Structure-function relationship in a variant hemoglobin: a combined computational-experimental approach

Dynamics of camel and human hemoglobin revealed by molecular simulations

Hemoglobin is one of the most widely studied proteins genetically, biochemically, and structurally. It is an oxygen carrying tetrameric protein that imparts the characteristic red color to blood. Each chain of hemoglobin harbors a heme group embedded in a hydrophobic pocket. Several studies have investigated structural variations present in mammalian hemoglobin and their functional implications. However, camel hemoglobin has not been thoroughly explored, especially from a structural perspective. Importantly, very little is known about how the heme group interacts with hemoglobin under varying conditions of osmolarity and temperature. Several experimental studies have indicated that the tense (T) state is more stable than the relaxed (R) state of hemoglobin under normal physiological conditions. Despite the fact that R state is less stable than the T state, no extensive structural dynamics studies have been performed to investigate global quaternary transitions of R state hemoglobin under normal physiological conditions. To evaluate this, several 500 ns all-atom molecular dynamics simulations were performed to get a deeper understanding of how camel hemoglobin behaves under stress, which it is normally exposed to, when compared to human hemoglobin. Notably, camel hemoglobin was more stable under physiological stress when compared to human hemoglobin. Additionally, when compared to camel hemoglobin, cofactor-binding regions of hemoglobin also exhibited more fluctuations in human hemoglobin under the conditions studied. Several differences were observed between the residues of camel and human hemoglobin that interacted with heme. Importantly, distal residues His58 of α hemoglobin and His63 of β hemoglobin formed more sustained interactions, especially at higher temperatures, in camel hemoglobin. These residues are important for oxygen binding to hemoglobin. Thus, this work provides insights into how camel and human hemoglobin differ in their interactions under stress.

Epistasis constrains mutational pathways of hemoglobin adaptation in high-altitude pikas

A fundamental question in evolutionary genetics concerns the roles of mutational pleiotropy and epistasis in shaping trajectories of protein evolution. This question can be addressed most directly by using site-directed mutagenesis to explore the mutational landscape of protein function in experimentally defined regions of sequence space. Here, we evaluate how pleiotropic trade-offs and epistatic interactions influence the accessibility of alternative mutational pathways during the adaptive evolution of hemoglobin (Hb) function in high-altitude pikas (Mammalia: Lagomorpha). By combining ancestral protein resurrection with a combinatorial protein-engineering approach, we examined the functional effects of sequential mutational steps in all possible pathways that produced an increased Hb–O₂ affinity. These experiments revealed that the effects of mutations on Hb–O₂ affinity are highly dependent on the temporal order in which they occur: Each of three β-chain substitutions produced a significant increase in Hb–O₂ affinity on the ancestral genetic background, but two of these substitutions produced opposite effects when they occurred as later steps in the pathway. The experiments revealed pervasive epistasis for Hb–O₂ affinity, but affinity-altering mutations produced no significant pleiotropic trade-offs. These results provide insights into the properties of adaptive substitutions in naturally evolved proteins and suggest that the accessibility of alternative mutational pathways may be more strongly constrained by sign epistasis for positively selected biochemical phenotypes than by antagonistic pleiotropy.

Structural and Functional Characterization of a New Double Variant Haemoglobin (HbG-Philadelphia/Duarte α(2)β(2))

WE REPORT THE FIRST CASE OF COSEGREGATION OF TWO HAEMOGLOBINS (HBS): HbG-Philadelphia [α68(E17)Asn → Lys] and HbDuarte [β62(E6)Ala → Pro]. The proband is a young patient heterozygous also for β°-thalassaemia. We detected exclusively two haemoglobin variants: HbDuarte and HbG-Philadelphia/Duarte. Functional study of the new double variant HbG-Philadelphia/Duarte exhibited an increase in oxygen affinity, with a slight decrease of cooperativity and Bohr effect. This functional behaviour is attributed to β62Ala → Pro instead of α68Asn → Lys substitution. Indeed, HbG-Philadelphia isolated in our laboratory from blood cells donor carrier for this variant is not affected by any functional modification, whereas purified Hb Duarte showed functional properties very similar to the double variant. NMR and MD simulation studies confirmed that the presence of Pro instead of Ala at the β62 position produces displacement of the E helix and modifications of the tertiary structure. The substitution α68(E17)Asn → Lys does not cause significant structural and dynamical modifications of the protein. A possible structure-based rational of substitution effects is suggested.

Breathing motions of a respiratory protein revealed by molecular dynamics simulations

Internal cavities, which are central to the biological functions of myoglobin, are exploited by gaseous ligands (e.g., O(2), NO, CO, etc.) to migrate inside the protein matrix. At present, it is not clear whether the ligand makes its own way inside the protein or instead the internal cavities are an intrinsic feature of myoglobin. To address this issue, standard molecular dynamics simulations were performed on horse-heart met-myoglobin with no ligand migrating inside the protein matrix. To reveal intrinsic internal pathways, the use of a statistical approach was applied to the cavity calculation, with special emphasis on the major pathway from the distal pocket to Xe1. Our study points out the remarkable dynamical behavior of Xe4, whose "breathing motions" may facilitate migration of ligands through the distal region. Additionally, our results highlight a two-way path for a ligand to diffuse through the proximal region, possibly allowing an alternative route in case Xe1 is occupied. Finally, our approach has led us to the identification of key residues, such as leucines, that may work as switches between cavities.

Evidences of xenon-induced structural changes in the active site of cyano-metmyoglobins: a 1H NMR study

Using xenon atoms as a biomolecular probe raises the concern of whether they may influence in some way the molecular and electronic structure of the system under study. In this paper, the relevance of guest-host interactions in xenon complexes with paramagnetic myoglobins (Mbs) is thoroughly analyzed, and the issue about the use of xenon to detect and characterize voids within flexible biomolecules is critically discussed. A detailed 1H NMR study useful for describing the hydrophobic cavities close to the active site of low-spin ferric myoglobins with respect to their interaction with the xenon atom is presented. The method is subsequently validated by the analysis of Xe-Mb with two different myoglobins, extracted from horse and pig. These myoglobins differ by 14 amino acids. One of these, Ile142 in horse Mb, is located in the proximal cavity, which is the main xenon binding site in horse Mb, and is replaced by Met142 in pig Mb. We demonstrated specific behaviors associated with the capacity of each of the two myoglobins to bind xenon and provided site-specific information on the host-guest interaction. Moreover, 1H NMR measurements produce a picture of xenon-related local distortions of the protein, associated with a functionally relevant residue located right at the active site, the proximal hystidine E7(His93). According to the 1H NMR data, xenon induces the tilt of the residue His93 relative to the heme plane and consequently causes an alteration of the magnetic axes. Similar conclusions are obtained both for pig cyano-myoglobin and for horse cyano-myoglobin, the structural deformation being in the former of minor entity.

Haemoglobin polymorphisms affect the oxygen-binding properties in Atlantic cod populations

A major challenge in evolutionary biology is to identify the genes underlying adaptation. The oxygen-transporting haemoglobins directly link external conditions with metabolic needs and therefore represent a unique system for studying environmental effects on molecular evolution. We have discovered two haemoglobin polymorphisms in Atlantic cod populations inhabiting varying temperature and oxygen regimes in the North Atlantic. Three-dimensional modelling of the tetrameric haemoglobin structure demonstrated that the two amino acid replacements Met55beta1Val and Lys62beta1Ala are located at crucial positions of the alpha1beta1 subunit interface and haem pocket, respectively. The replacements are proposed to affect the oxygen-binding properties by modifying the haemoglobin quaternary structure and electrostatic feature. Intriguingly, the same molecular mechanism for facilitating oxygen binding is found in avian species adapted to high altitudes, illustrating convergent evolution in water- and air-breathing vertebrates to reduction in environmental oxygen availability. Cod populations inhabiting the cold Arctic waters and the low-oxygen Baltic Sea seem well adapted to these conditions by possessing the high oxygen affinity Val55-Ala62 haplotype, while the temperature-insensitive Met55-Lys62 haplotype predominates in the southern populations. The distinct distributions of the functionally different haemoglobin variants indicate that the present biogeography of this ecologically and economically important species might be seriously affected by global warming.