Analysis of Fluctuation in the Heme-Binding Pocket and Heme Distortion in Hemoglobin and Myoglobin

Heme is located in the active site of proteins and has diverse and important biological functions, such as electron transfer and oxygen transport and/or storage. The distortion of heme porphyrin is considered an important factor for the diverse functions of heme because it correlates with the physical properties of heme, such as oxygen affinity and redox potential. Therefore, clarification of the relationship between heme distortion and the protein environment is crucial in protein science. Here, we analyzed the fluctuation in heme distortion in the protein environment for hemoglobin and myoglobin using molecular dynamics (MD) simulations and quantum mechanical (QM) calculations. We also investigated the protein structures of hemoglobin and myoglobin stored in Protein Data Bank and found that heme is distorted along the doming mode, which correlates with its oxygen affinity, more prominently in the protein environment than in the isolated state, and the magnitude of distortion is different between hemoglobin and myoglobin. This tendency was also observed in the results of MD simulations and QM calculations. These results suggest that heme distortion is affected by its protein environment and fluctuates around its fitted conformation, leading to physical properties that are appropriate for protein functions.

Ultrafast dynamics of heme distortion in the O2-sensor of a thermophilic anaerobe bacterium

Heme-Nitric oxide and Oxygen binding protein domains (H-NOX) are found in signaling pathways of both prokaryotes and eukaryotes and share sequence homology with soluble guanylate cyclase, the mammalian NO receptor. In bacteria, H-NOX is associated with kinase or methyl accepting chemotaxis domains. In the O2-sensor of the strict anaerobe Caldanaerobacter tengcongensis (Ct H-NOX) the heme appears highly distorted after O2 binding, but the role of heme distortion in allosteric transitions was not yet evidenced. Here, we measure the dynamics of the heme distortion triggered by the dissociation of diatomics from Ct H-NOX using transient electronic absorption spectroscopy in the picosecond to millisecond time range. We obtained a spectroscopic signature of the heme flattening upon O2 dissociation. The heme distortion is immediately (<1 ps) released after O2 dissociation to produce a relaxed state. This heme conformational change occurs with different proportions depending on diatomics as follows: CO < NO < O2. Our time-resolved data demonstrate that the primary structural event of allostery is the heme distortion in the Ct H-NOX sensor, contrastingly with hemoglobin and the human NO receptor, in which the primary structural events are respectively the motion of the proximal histidine and the rupture of the iron-histidine bond.

Interactive regulation between aliphatic hydroxylation and aromatic hydroxylation of thaxtomin D in TxtC: a theoretical investigation

ABSTRACTTxtC is an unusual bifunctional cytochrome P450 that is able to perform sequential aliphatic and aromatic hydroxylation of the diketopiperazine substrate thaxtomin D in two remote sites to produce thaxtomin A. Though the X-ray structure of TxtC complexed with thaxtomin D revealed a binding mode for its aromatic hydroxylation, the preferential hydroxylation site is aliphatic C14. It is thus intriguing to unravel how TxtC accomplishes such two-step catalytic hydroxylation on distinct aliphatic and aromatic carbons and why the aliphatic site is preferred in the hydroxylation step. In this work, by employing molecular docking and molecular dynamics (MD) simulation, we revealed that thaxtomin D could adopt two different conformations in the TxtC active site, which were equal in energy with either the aromatic C-H or aliphatic C14-H laying towards the active Cpd I oxyferryl moiety. Further ONIOM calculations indicated that the energy barrier for the rate-limiting hydroxylation step on the aliphatic C14 site was 8.9 kcal/mol more favorable than that on the aromatic C20 site. The hydroxyl group on the monohydroxylated intermediate thaxtomin B C14 site formed hydrogen bonds with Ser280 and Thr385, which induced the L-Phe moiety to rotate around the Cβ−Cγ bond of the 4-nitrotryptophan moiety. Thus, it adopted an energy favorable conformation with aromatic C20 adjacent to the oxyferryl moiety. In addition, the hydroxyl group induced solvent water molecules to enter the active site, which propelled thaxtomin B towards the heme plane and resulted in heme distortion. Based on this geometrical layout, the rate-limiting aromatic hydroxylation energy barrier decreased to 15.4 kcal/mol, which was comparable to that of the thaxtomin D aliphatic hydroxylation process. Our calculations indicated that heme distortion lowered the energy level of the lowest Cpd I α-vacant orbital, which promoted electron transfer in the rate-limiting thaxtomin B aromatic hydroxylation step in TxtC.

Iron transitions during activation of allosteric heme proteins in cell signaling

Allosteric heme proteins can fulfill a very large number of different functions, thanks to the remarkable chemical versatility of the heme through the entire living kingdom. The involved heme properties include: changes of iron coordination and redox state, heme macrocycle electrostatic interactions and heme distortion.

Fractional transfer of a free unpaired electron to overcome energy barriers in the formation of Fe4+ from Fe3+ during the core contraction of macrocycles: implication for heme distortion

The free unpaired electron in Fe3+ ions cannot be directly removed, and needs a transfer pathway with at least four steps to overcome the high energy barriers to form Fe4+ ions.