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
Differential responses in the core, active site and peripheral regions of cytochrome c peroxidase to extreme pressure and temperature
In consideration of life in extreme environments, the effects of hydrostatic pressure on proteins at the atomic level have drawn substantial interest. Large deviations of temperature and pressure from ambient conditions can shift the free energy landscape of proteins to reveal otherwise lowly populated structural states and even promote unfolding. We report the crystal structure of the heme-containing peroxidase, cytochrome c peroxidase (CcP) at 1.5 and 3.0 kbar and make comparisons to structures determined at 1.0 bar and cryo-temperatures (100 K). Compressibility plateaus after 1.5 kbar and pressure produces anisotropic changes in CcP. CcP responds to pressure with volume declines at the periphery of the protein where B-factors are relatively high but maintains nearly intransient core structure and active site channels. Compression at the surface affects neither alternate side-chain conformers nor B-factors. Thus, packing in the core, which resembles a crystalline solid, limits motion and protects the active site, whereas looser packing at the surface preserves side-chain dynamics. Changes in active-site solvation and heme ligation reveal pressure sensitivity to protein-ligand interactions and reveal a potential docking site for the substrate peroxide.
more »
« less
- Award ID(s):
- 2129728
- PAR ID:
- 10526429
- Publisher / Repository:
- bioRxiv
- Date Published:
- Format(s):
- Medium: X
- Institution:
- bioRxiv
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Bernd Reif (Ed.)It has long been known that the alteration of protein side chains that occlude or expose the heme cofactor to water can greatly affect the stability of the oxyferrous heme state. Here, we demonstrate that the rate of dynamically driven water penetration into the core of an artificial oxygen transport protein also correlates with oxyferrous state lifetime by reducing global dynamics, without altering the structure of the active site, via the simple linking of the two monomers in a homodimeric artificial oxygen transport protein using a glycine-rich loop. The tethering of these two helices does not significantly affect the active site structure, pentacoordinate heme-binding affinity, reduction potential, or gaseous ligand affinity. It does, however, significantly reduce the hydration of the protein core, as demonstrated by resonance Raman spectroscopy, backbone amide hydrogen exchange, and pKa shifts in buried histidine side chains. This further destabilizes the charge-buried entatic state and nearly triples the oxyferrous state lifetime. These data are the first direct evidence that dynamically driven water penetration is a rate-limiting step in the oxidation of these complexes. It furthermore demonstrates that structural rigidity that limits water penetration is a critical design feature in metalloenzyme construction and provides an explanation for both the failures and successes of earlier attempts to create oxygen-binding proteins.more » « less
-
We present a comprehensive investigation on the different role of CO in carboxy- neuroglobin i) as ligand of the heme group in the active site forming a bond with the heme iron and ii) dissociated from the heme group but still trapped inside the active site, focusing on two specific orientations, one with CO perpendicular to the plane defined by the distal histidine of the enzyme (form A) and one with CO located parallel to that plane (form B). Our study includes wild type carboxy-neuroglobin and nine known protein mutations. Considering that the distal histidine interacting with the heme group can adapt two different tautomeric forms and the two possible orientations of the dissociated CO, a total of 36 protein systems were analyzed in this study. Fully optimized geometries and vibrational frequencies were calculated at the QM/MM level, followed by the local mode analysis, to decode CO bond properties. The intrinsic bond strengths derived from the local mode analysis, complemented with NBO and QTAIM data, reveal that the strength of the CO bond, in the hexacoordinate (where CO is a ligand of the heme group) and pentacoordinate (where CO is dissociated from the heme group) scenarios, is dominated by through bond and through space charge transfer between CO and Fe, fine-tuned by electrostatic and dispersion interactions with the side chain amino acids in the distal heme pocket. Suggestions are made as to advise on how protein modifications can influence the molecular properties of the coordinated or dissociated CO, which could serve the fine-tuning of existing and the design of new neuroglobin models with specific FeC and CO bond strengths.more » « less
-
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
-
The primary and secondary coordination spheres of metal binding sites in metalloproteins have been investigated extensively, leading to the creation of high-performing functional metalloproteins; however, the impact of the overall structure of the protein scaffold on the unique properties of metalloproteins has rarely been studied. A primary example is the binuclear CuA center, an electron transfer cupredoxin domain of photosynthetic and respiratory complexes and, recently, a protein co-regulated with particulate methane and ammonia monooxygenases. The redox potential, Cu–Cu spectroscopic features, and a valence delocalized state of CuA are difficult to reproduce in synthetic models, and every artificial protein CuA center to-date has used a modified cupredoxin. Here we present a fully functional CuA center designed in a structurally non-homologous protein, cytochrome c peroxidase (CcP), by only two mutations (CuACcP). We demonstrate with UV-visible absorption, resonance Raman, and MCD spectroscopy that CuACcP is valence delocalized. CW and pulsed (HYSCORE) X-band EPR show it has a highly compact gz area and small Az hyperfine principal value with g and A tensors that resemble axially perturbed CuA. Stopped-flow kinetics found that CuA formation proceeds through a single T2Cu intermediate. The reduction potential of CuACcP is comparable to native CuA and can transfer electrons to a physiological redox partner. We built a structural model of the designed Cu binding site from EXAFS and validated it by mutation of coordinating Cys and His residues, revealing that a triad of residues (R48C, W51C, and His52) rigidly arranged on one α-helix is responsible for chelating the first Cu atom and that His175 stabilizes the binuclear complex by rearrangement of the CcP heme-coordinating helix. This design is a demonstration that a highly conserved protein fold is not uniquely necessary to induce certain characteristic physical and chemical properties to a metal redox center.more » « less
An official website of the United States government

