skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Exploring the folding energy landscapes of heme proteins using a hybrid AWSEM-heme model
Abstract Heme is an active center in many proteins. Here we explore computationally the role of heme in protein folding and protein structure. We model heme proteins using a hybrid model employing the AWSEM Hamiltonian, a coarse-grained forcefield for the protein chain along with AMBER, an all-atom forcefield for the heme. We carefully designed transferable force fields that model the interactions between the protein and the heme. The types of protein–ligand interactions in the hybrid model include thioester covalent bonds, coordinated covalent bonds, hydrogen bonds, and electrostatics. We explore the influence of different types of hemes (heme b and heme c) on folding and structure prediction. Including both types of heme improves the quality of protein structure predictions. The free energy landscape shows that both types of heme can act as nucleation sites for protein folding and stabilize the protein folded state. In binding the heme, coordinated covalent bonds and thioester covalent bonds for heme c drive the heme toward the native pocket. The electrostatics also facilitates the search for the binding site.  more » « less
Award ID(s):
2019745
PAR ID:
10363122
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
Springer Science + Business Media
Date Published:
Journal Name:
Journal of Biological Physics
Volume:
48
Issue:
1
ISSN:
0092-0606
Page Range / eLocation ID:
p. 37-53
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Schneidman-Duhovny, Dina (Ed.)
    We present OpenAWSEM and Open3SPN2, new cross-compatible implementations of coarse-grained models for protein (AWSEM) and DNA (3SPN2) molecular dynamics simulations within the OpenMM framework. These new implementations retain the chemical accuracy and intrinsic efficiency of the original models while adding GPU acceleration and the ease of forcefield modification provided by OpenMM’s Custom Forces software framework. By utilizing GPUs, we achieve around a 30-fold speedup in protein and protein-DNA simulations over the existing LAMMPS-based implementations running on a single CPU core. We showcase the benefits of OpenMM’s Custom Forces framework by devising and implementing two new potentials that allow us to address important aspects of protein folding and structure prediction and by testing the ability of the combined OpenAWSEM and Open3SPN2 to model protein-DNA binding. The first potential is used to describe the changes in effective interactions that occur as a protein becomes partially buried in a membrane. We also introduced an interaction to describe proteins with multiple disulfide bonds. Using simple pairwise disulfide bonding terms results in unphysical clustering of cysteine residues, posing a problem when simulating the folding of proteins with many cysteines. We now can computationally reproduce Anfinsen’s early Nobel prize winning experiments by using OpenMM’s Custom Forces framework to introduce a multi-body disulfide bonding term that prevents unphysical clustering. Our protein-DNA simulations show that the binding landscape is funneled towards structures that are quite similar to those found using experiments. In summary, this paper provides a simulation tool for the molecular biophysics community that is both easy to use and sufficiently efficient to simulate large proteins and large protein-DNA systems that are central to many cellular processes. These codes should facilitate the interplay between molecular simulations and cellular studies, which have been hampered by the large mismatch between the time and length scales accessible to molecular simulations and those relevant to cell biology. 
    more » « less
  2. Abstract Charged residues on the surface of proteins are critical for both protein stability and interactions. However, many proteins contain binding regions with a high net charge that may destabilize the protein but are useful for binding to oppositely charged targets. We hypothesized that these domains would be marginally stable, as electrostatic repulsion would compete with favorable hydrophobic collapse during folding. Furthermore, by increasing the salt concentration, we predict that these protein folds would be stabilized by mimicking some of the favorable electrostatic interactions that take place during target binding. We varied the salt and urea concentrations to probe the contributions of electrostatic and hydrophobic interactions for the folding of the yeast SH3 domain found in Abp1p. The SH3 domain was significantly stabilized with increased salt concentrations due to Debye–Huckel screening and a nonspecific territorial ion‐binding effect. Molecular dynamics and NMR show that sodium ions interact with all 15 acidic residues but do little to change backbone dynamics or overall structure. Folding kinetics experiments show that the addition of urea or salt primarily affects the folding rate, indicating that almost all the hydrophobic collapse and electrostatic repulsion occur in the transition state. After the transition state formation, modest yet favorable short‐range salt bridges are formed along with hydrogen bonds, as the native state fully folds. Thus, hydrophobic collapse offsets electrostatic repulsion to ensure this highly charged binding domain can still fold and be ready to bind to its charged peptide targets, a property that is likely evolutionarily conserved over 1 billion years. 
    more » « less
  3. In our continued investigations of microbial globins, we solved the structure of a truncated hemoglobin from Shewanella benthica, an obligate psychropiezophilic bacterium. The distal side of the heme active site is lined mostly with hydrophobic residues, with the exception of a tyrosine, Tyr34 (CD1) and a histidine, His24 (B13). We found that purified SbHbN, when crystallized in the ferric form with polyethylene glycol as precipitant, turned into a green color over weeks. The electron density obtained from the green crystals accommodated a trans heme d, a chlorin-type derivative featuring a γ-spirolactone and a vicinal hydroxyl group on a pyrroline ring. In solution, exposure of the protein to one equivalent of hydrogen peroxide resulted in a similar green color change, but caused by the formation of multiple products. These were oxidation species released on protein denaturation, likely including heme d, and a species with heme covalently attached to the polypeptide. The Tyr34Phe replacement prevented the formation of both heme d and the covalent linkage. The ready modification of heme b by SbHbN expands the range of chemistries supported by the globin fold and offers a route to a novel heme cofactor. 
    more » « less
  4. Champion, Patricia A (Ed.)
    ABSTRACT Cellular life relies on enzymes that require metals, which must be acquired from extracellular sources. Bacteria utilize surface and secreted proteins to acquire such valuable nutrients from their environment. These include the cargo proteins of the type eleven secretion system (T11SS), which have been connected to host specificity, metal homeostasis, and nutritional immunity evasion. This Sec-dependent, Gram-negative secretion system is encoded by organisms throughout the phylum Proteobacteria, including human pathogensNeisseria meningitidis, Proteus mirabilis, Acinetobacter baumannii,andHaemophilus influenzae. Experimentally verified T11SS-dependent cargo includetransferrin-bindingprotein B (TbpB), the hemophilin homologshemereceptorprotein C (HrpC),hemophilinA(HphA), the immune evasion proteinfactor-H bindingprotein (fHbp), and the host symbiosis factornematodeintestinallocalization protein C (NilC). Here, we examined the specificity of T11SS systems for their cognate cargo proteins using taxonomically distributed homolog pairs of T11SS and hemophilin cargo and explored the ligand binding ability of those hemophilin cargo homologs.In vivoexpression inEscherichia coliof hemophilin homologs revealed that each is secreted in a specific manner by its cognate T11SS protein. Sequence analysis and structural modeling suggest that all hemophilin homologs share an N-terminal ligand-binding domain with the same topology as the ligand-binding domains of theHaemophilus haemolyticusheme binding protein (Hpl) and HphA. We term this signature feature of this group of proteins the hemophilin ligand-binding domain. Network analysis of hemophilin homologs revealed five subclusters and representatives from four of these showed variable heme-binding activities, which, combined with sequence-structure variation, suggests that hemophilins are diversifying in function.IMPORTANCEThe secreted protein hemophilin and its homologs contribute to the survival of several bacterial symbionts within their respective host environments. Here, we compared taxonomically diverse hemophilin homologs and their paired Type 11 secretion systems (T11SS) to determine if heme binding and T11SS secretion are conserved characteristics of this family. We establish the existence of divergent hemophilin sub-families and describe structural features that contribute to distinct ligand-binding behaviors. Furthermore, we demonstrate that T11SS are specific for their cognate hemophilin family cargo proteins. Our work establishes that hemophilin homolog-T11SS pairs are diverging from each other, potentially evolving into novel ligand acquisition systems that provide competitive benefits in host niches. 
    more » « less
  5. Abstract Engineered heme proteins were developed to possess numerous excellent biocatalytic nitrenoid C−H functionalizations. Computational approaches such as density functional theory (DFT), hybrid quantum mechanics/molecular mechanics (QM/MM), and molecular dynamics (MD) calculations were employed to help understand some important mechanistic aspects of these heme nitrene transfer reactions. This review summarizes advances of computational reaction pathway results of these biocatalytic intramolecular and intermolecular C−H aminations/amidations, focusing on mechanistic origins of reactivity, regioselectivity, enantioselectivity, diastereoselectivity as well as effects of substrate substituent, axial ligand, metal center, and protein environment. Some important common and distinctive mechanistic features of these reactions were also described with brief outlook of future development. 
    more » « less