Exploring the molecular basis of heme coordination in human neuroglobin

Neuroglobin (Ngb), a recently discovered ancient heme protein, presents the typical globin fold and is around 20% identical to myoglobin (Mb). In contrast with Mb, however, its heme is hexacoordinated (6c). It is expressed in the nervous system and has been the subject of numerous investigations in the last years, but its function is still unclear. The proposed roles include oxygen transport, reactive oxygen species (ROS) detoxification, hypoxia protection, and redox state sensing. All proposed functions require distal histidine dissociation from the heme to yield a reactive iron. With the aim of understanding the 6c to 5c transition, we have performed molecular dynamics simulations for ferrous Ngb in the 6c, 5c, and oxy states. We also computed free energy profiles associated with the transition employing an advanced sampling technique. Finally, we studied the effect of the redox state of CysCD7 and CysD5, which are known to form a disulfide bridge. Our results show that protein oxidation promotes a stabilization of the pentacoordinated species, thus favoring the protein to adopt the more reactive state and supporting the existence of a molecular mechanism whereby O2 would be released under hypoxic conditions, thereby suggesting an O(2) storage function for Ngb. Taken together, our results provide structural information not available experimentally which may shed light on the protein proposed functions, particularly as a redox sensor.     

Structural and kinetic characterization of myoglobins from eurythermal and stenothermal fish species

Teleost myoglobin (Mb) proteins from four fish species inhabiting different temperature environments were used to investigate the relationship between protein function and thermal stability. Mb was isolated from yellowfin tuna (homeothermal warm), mackerel (eurythermal warm), and the Antarctic teleost Notothenia coriiceps (stenothermal cold). Zebrafish (stenothermal tropical) myoglobin was expressed from cloned cDNA. These proteins differed in oxygen affinity, as measured by O2 dissociation rates and P50 values, and thermal stability as measured by autooxidation rates. Mackerel Mb had the highest P50 value at 25 degrees C (3.7 mmHg), corresponding to the lowest O2 affinity, followed by zebrafish (1.0 mmHg), yellowfin tuna (1.0 mmHg), and N. coriiceps (0.6 mmHg). Oxygen dissociation rates and Arrhenius plots were similar between all teleost species in this study, with the exception of mackerel myoglobin, which was two-fold faster at all temperatures tested. Myoglobin from the Antarctic teleost had the highest autooxidation rate (0.44 h(-1)), followed by mackerel (0.26 h(-1)), zebrafish (0.22 h(-1)), and yellowfin tuna (0.088 h(-1)). Primary structural analysis revealed residue differences distributed throughout the polypeptide sequences, making it difficult to identify, which, if any, residues contribute to structural flexibility. However, analysis of molecular dynamics trajectories indicates that Mb from the eurythermal mackerel is predicted to be the most flexible protein within the D loop and FG turn. At the same time, it has the lowest O2 affinity and the highest O2 dissociation rates when compared to myoglobins from teleosts that appear to be less flexible in our dynamics simulations.     

Rational design of a nitrite reductase based on myoglobin: a molecular modeling and dynamics simulation study

Myoglobin (Mb) is an ideal scaffold protein for rational protein design mimicking native enzymes. We recently designed a nitrite reductase (NiR) based on sperm whale Mb by introducing an additional distal histidine (Leu29 to His29 mutation) and generating a distal tyrosine (Phe43 to Tyr43 mutation) in the heme pocket, namely L29H/F43Y Mb, to mimic the active site of cytochrome cd (1) NiR from Ps. aeruginosa that contains two distal histidines and one distal tyrosine. The molecular modeling and dynamics simulation study herein revealed that L29H/F43Y Mb has the necessary structural features of native cytochrome cd (1) NiR and can provide comparable interactions with nitrite as in native NiRs, which provides rationality for the protein design and guides the protein engineering. Additionally, the present study provides an insight into the relatively low NiR activity of Mb in biological systems.     

Relaxations and fluctuations in myoglobin

One major goal of biological physics is the discovery and understanding of the concepts and laws that govern biomolecules, in particular proteins. Since there exist at least 10(5) different proteins, the choice of a suitable prototype is necessary. Myoglobin (Mb) has for many years played the role of such a prototype. It appears to be simple enough so that many of its properties can be understood, yet it is complex enough to display many of the fascinating characteristics of biomolecules. One major achievement in the study of any protein would be the establishment of convincing connections among structure, kinetics, energy landscape, dynamics, and function. We believe that this goal has not yet been reached in any protein, but the present knowledge of Mb gives some hope that the end is near in this case. Here, we sketch some of the results that have been obtained in the past 50 or more years in the research on Mb, obtained by an army of investigators.     

Molecular dynamics simulations indicate that tyrosineB10 limits motions of distal histidine to regulate CO binding in soybean leghemoglobin

Myoglobin (Mb) uses strong electrostatic interaction in its distal heme pocket to regulate ligand binding. The mechanism of regulation of ligand binding in soybean leghemoglobin a (Lba) has been enigmatic and more so due to the absence of gaseous ligand bound atomic resolution three‐dimensional structure of the plant globin. While the 20‐fold higher oxygen affinity of Lba compared with Mb is required for its dual physiological function, the mechanism by which this high affinity is achieved is only emerging. Extensive mutational analysis combined with kinetic and CO‐FT‐IR spectroscopic investigation led to the hypothesis that Lba depended on weakened electrostatic interaction between distal HisE7 and bound ligand achieved by invoking B10Tyr, which itself hydrogen bonds with HisE7 thus restricting it in a single conformation detrimental to Mb‐like strong electrostatic interaction. Such theory has been re‐assessed here using CO‐Lba in silico model and molecular dynamics simulation. The investigation supports the presence of at least two major conformations of HisE7 in Lba brought about by imidazole ring flip, one of which makes hydrogen bonds effectively with B10Tyr affecting the former's ability to stabilize bound ligand, while the other does not. However, HisE7 in Lba has limited conformational freedom unlike high frequency of imidazole ring flips observed in Mb and in TyrB10Leu mutant of Lba. Thus, it appears that TyrB10 limits the conformational freedom of distal His in Lba, tuning down ligand dissociation rate constant by reducing the strength of hydrogen bonding to bound ligand, which the freedom of distal His of Mb allows. Proteins 2015; 83:1836–1848. © 2015 Wiley Periodicals, Inc.     

Electrochemical nitration of myoglobin at tyrosine 103: Structure and stability

Nitration in proteins is a physiologically relevant process and the formation of 3-nitrotyrosine was first proposed as an in vivo marker of the production of reactive nitrogen species in oxidative stress. No studies have been published on structural changes associated with nitration of myoglobin. To address this deficiency the electrochemical nitration of equine skeletal muscle (Mb) at amino acid tyrosine 103 has been investigated for the evaluation and characterization of structural and thermal stability changes. Y103 in Mb is one of the most exposed tyrosine residues and it is also close to the heme group. Effects of Y103 nitration on the secondary and tertiary structure of Y103 have been studied by UV–Vis, circular dichroism, fluorescence and NMR spectroscopy and by electrochemical studies. At physiological pH, subtle changes were observed involving slight loosening of the tertiary structure and conformational exchange processes. Thermal stability of the nitrated protein was found to be reduced by 5 °C for the nitrated Mb compared with the native Mb at physiological pH. Altogether, NMR data indicates that nitrated Mb has a very similar tertiary structure to that of native Mb, although with a slightly open conformation.     

Non-enzymatic glycation induces structural modifications of myoglobin

Increased glucose concentration in diabetes mellitus causes glycation of several proteins, leading to changes in their properties. Although glycation-induced functional modification of myoglobin is known, structural modification of the protein has not yet been reported. Here, we have studied glucose-modified structural changes of the heme protein. After in vitro glycation of metmyoglobin (Mb) by glucose at 25°C for 6 days, glycated myoglobin (GMb) and unchanged Mb have been separated by ion exchange (BioRex 70) chromatography, and their properties have been compared. Compared to Mb, GMb exhibits increased absorbance around 280 nm and enhanced fluorescence emission with excitation at 285 nm. Fluorescence quenching experiments of the proteins by acrylamide and KI indicate that more surface accessible tryptophan residues are exposed in GMb. CD spectroscopic study reveals a change in the secondary structure of GMb with decreased α-helix content. 1-anilino-naphthaline-8-sulfonate (ANS) binding with Mb and GMb indicates that glycation increases hydrophobicity of the heme protein. GMb appears to be less stable with respect to thermal denaturation and differential calorimetry experiments. Heme-globin linkage becomes weaker in GMb, as shown by spectroscopic and gel electrophoresis experiments. A correlation between glycation-induced structural and functional modifications of the heme protein has been suggested.     

Diffracted X-ray tracking method for recording single-molecule protein motions

BackgroundProteins change their conformation depending on function. Although a vast number of static pictures of proteins have been accumulated, information regarding their dynamics in function is limited. Diffracted X-ray tracking (DXT) is a good candidate to obtain the missing data.Scope of reviewA gold nanocrystal was attached to the target protein as a probe and the motion of the X-ray diffraction spots from the crystal corresponded to the motion of the target. Although it has advantages of high temporal (sub-millisecond) and spatial (approximately 0.1°) resolutions, it is not extensively utilized. This review focused on its effective application from a user's perspective. We also present an example with the KcsA channel and the status of recent developments to show the future possibilities of the method.Major conclusionsDXT is a powerful method to investigate intramolecular structural changes. For instance, in the KcsA channel, the method revealed a wave of conformational changes transmitted from the gate region to the end of the molecule. The method is continuously being developed, and users can choose an appropriate measurement system depending on the condition of their sample.General significanceRevealing the protein structural changes with respect to function is an important frontier. The most distinctive feature of the DXT method is that both high temporal and spatial resolutions are achievable, and it is possible to track the motions of multiple molecules at the same time. This feature is an advantage for screening molecules associated with the target proteins (e.g., ligands and medicines).     

Identification of Selection Footprints on the X Chromosome in Pig

Identifying footprints of selection can provide a straightforward insight into the mechanism of artificial selection and further dig out the causal genes related to important traits. In this study, three between-population and two within-population approaches, the Cross Population Extend Haplotype Homozygosity Test (XPEHH), the Cross Population Composite Likelihood Ratio (XPCLR), the F-statistics (Fst), the Integrated Haplotype Score (iHS) and the Tajima''s D, were implemented to detect the selection footprints on the X chromosome in three pig breeds using Illumina Porcine60K SNP chip. In the detection of selection footprints using between-population methods, 11, 11 and 7 potential selection regions with length of 15.62 Mb, 12.32 Mb and 9.38 Mb were identified in Landrace, Chinese Songliao and Yorkshire by XPEHH, respectively, and 16, 13 and 17 potential selection regions with length of 15.20 Mb, 13.00 Mb and 19.21 Mb by XPCLR, 4, 2 and 4 potential selection regions with length of 3.20 Mb, 1.60 Mb and 3.20 Mb by Fst. For within-population methods, 7, 10 and 9 potential selection regions with length of 8.12 Mb, 8.40 Mb and 9.99 Mb were identified in Landrace, Chinese Songliao and Yorkshire by iHS, and 4, 3 and 2 potential selection regions with length of 3.20 Mb, 2.40 Mb and 1.60 Mb by Tajima''s D. Moreover, the selection regions from different methods were partly overlapped, especially the regions around 22 ∼25 Mb were detected under selection in Landrace and Yorkshire while no selection in Chinese Songliao by all three between-population methods. Only quite few overlap of selection regions identified by between-population and within-population methods were found. Bioinformatics analysis showed that the genes relevant with meat quality, reproduction and immune were found in potential selection regions. In addition, three out of five significant SNPs associated with hematological traits reported in our genome-wide association study were harbored in potential selection regions.     

Neutron Spin-Echo Studies of Hemoglobin and Myoglobin: Multiscale Internal Dynamics

Neutron spin-echo spectroscopy was used to study structural fluctuations that occur in hemoglobin (Hb) and myoglobin (Mb) in solution. Using neutron spin-echo data up to a very high momentum transfer q (∼ 0.62 Å^− 1), we characterized the internal dynamics of these proteins at the levels of dynamic pair correlation function and self-correlation function in the time range of several picoseconds to a few nanoseconds. In the same protein solution, data transition from pair correlation motion to self-correlation motion as the momentum transfer q increases. At low q, coherent scattering dominates; at high q, observations are largely due to incoherent scattering. The low q data were interpreted in terms of an effective diffusion coefficient; on the other hand, the high q data were interpreted in terms of mean square displacements. Comparison of data from the two homologous proteins collected at different temperatures and protein concentrations was used to assess the contributions made by translational and rotational diffusion and internal modes of motion to the data. The temperature dependence of decay times can be attributed to changes in the viscosity and temperature of the solvent, as predicted by the Stokes-Einstein relationship. This is true for contributions from both diffusive and internal modes of motion, indicating an intimate relationship between the internal dynamics of the proteins and the viscosity of the solvent. Viscosity change associated with protein concentration can account for changes in diffusion observed at different concentrations, but is apparently not the only factor involved in the changes in internal dynamics observed with change in protein concentration. Data collected at high q indicate that internal modes in Mb are generally faster than those in Hb, perhaps due to the greater surface-to-volume ratio of Mb and the fact that surface groups tend to exhibit faster motion than buried groups. Comparison of data from Hb and data from Mb at low q indicates an unexpectedly rapid motion of Hb αβ dimers relative to one another. Dynamic motion of subunits is increasingly perceived as important to the allosteric behavior of Hb. Our data demonstrate that this motion is highly sensitive to protein concentration, temperature, and solvent viscosity, indicating that great care needs to be exercised in interpreting its effect on protein function.     

Coupling of protein and environment fluctuations

We review the concepts of protein dynamics developed over the last 35years and extend applications of the unified model of protein dynamics to heat flow and spatial fluctuations in hydrated myoglobin (Mb) powders. Differential scanning calorimetry (DSC) and incoherent neutron scattering (INS) data on hydration Mb powders are explained by the temperature-dependence of the hydration-shell β(h) process measured by dielectric relaxation spectroscopy (DRS). The unified model explains the temperature dependence of DSC and INS data as a kinetic effect due to a fixed experimental time window and a broad distribution of hydration-shell β(h) fluctuation rates. We review the slaving of large scale protein motions to the bulk solvent α process, and the metastability of Mb molecules in glass forming bulk solvent at low temperatures. This article is part of a Special Issue entitled: "Protein Dynamics: Experimental and Computational Approaches".     

Dynamics comparison of two myoglobins with a distinct heme active site

Sperm whale myoglobin (swMb) is a well-studied heme protein, both experimentally and theoretically. Comparatively, little attention has been paid to another member of Mb family, Aplysia limacina myoglobin (apMb). swMb and apMb have the same overall structure and perform the same biological function, i.e., O₂ carrier, while using a distinct heme active site. To provide insights into the structure-function relationship for these two Mbs, we herein made a comparison in terms of their dynamics properties using molecular dynamics simulation. We analyzed the overall structure and protein motions, as well as the intramolecular contacts, namely salt-bridges and hydrogen bonds, especially the interactions between the heme propionate groups and the polypeptide chain. The internal cavities in apMb were also compared to the well-known xenon and other cavities in swMb. Based on current simulations, we propose a unique “arginine gate” for apMb, which has a similar function to the histidine gate observed for swMb in previous studies. This study provides valuable information for understanding the homology of heme proteins, and also aids in rational design of structural and functional heme proteins by alternating the heme active site.     

Evolutionarily conserved linkage between enzyme fold, flexibility, and catalysis

Proteins are intrinsically flexible molecules. The role of internal motions in a protein''s designated function is widely debated. The role of protein structure in enzyme catalysis is well established, and conservation of structural features provides vital clues to their role in function. Recently, it has been proposed that the protein function may involve multiple conformations: the observed deviations are not random thermodynamic fluctuations; rather, flexibility may be closely linked to protein function, including enzyme catalysis. We hypothesize that the argument of conservation of important structural features can also be extended to identification of protein flexibility in interconnection with enzyme function. Three classes of enzymes (prolyl-peptidyl isomerase, oxidoreductase, and nuclease) that catalyze diverse chemical reactions have been examined using detailed computational modeling. For each class, the identification and characterization of the internal protein motions coupled to the chemical step in enzyme mechanisms in multiple species show identical enzyme conformational fluctuations. In addition to the active-site residues, motions of protein surface loop regions (>10 Å away) are observed to be identical across species, and networks of conserved interactions/residues connect these highly flexible surface regions to the active-site residues that make direct contact with substrates. More interestingly, examination of reaction-coupled motions in non-homologous enzyme systems (with no structural or sequence similarity) that catalyze the same biochemical reaction shows motions that induce remarkably similar changes in the enzyme–substrate interactions during catalysis. The results indicate that the reaction-coupled flexibility is a conserved aspect of the enzyme molecular architecture. Protein motions in distal areas of homologous and non-homologous enzyme systems mediate similar changes in the active-site enzyme–substrate interactions, thereby impacting the mechanism of catalyzed chemistry. These results have implications for understanding the mechanism of allostery, and for protein engineering and drug design.

Author''s Summary

Enzymes are nature''s molecular machines that catalyze biochemical reactions with remarkable efficiency. Recent evidence suggests that enzyme function may involve not only direct structural interactions between the enzyme and its substrate, but also internal motions of the enzyme itself. Here, we describe a computational investigation of three classes of enzymes that catalyze completely different biochemical reactions. Remarkably, the mobile enzyme regions and the nature of these motions are the same across species ranging from single-celled organisms to complex life-forms. Also surprisingly, non-homologous enzymes that catalyze the same chemical reaction but do not share sequence or structural similarity reveal a similar impact of enzyme motions on their reaction mechanisms. Flexible enzyme regions are found to be connected by conserved networks of coupled interactions that connect surface regions to active-site residues. These networks may provide a mechanism for the solvent on an enzyme''s surface to couple to the reaction catalyzed by the enzyme. These results have implications for understanding the mechanism of allostery (long-range effects), and for protein engineering and drug design.     

Theoretical studies of viral capsid proteins

Recent results in structural biology and increases in computer power have prompted initial theoretical studies on capsids of nonenveloped icosahedral viruses. The macromolecular assembly of 60 to 180 protein copies into a protein shell results in a structure of considerable size for molecular dynamics simulations. Nonetheless, progress has been made in examining these capsid assemblies from molecular dynamics calculations and kinetic models. The goals of these studies are to understand capsid function and structural properties, including quarternary structural stability, effects of antiviral compounds that bind the capsid and the self-assembly process. The insight that can be gained from the detailed information provided by simulations is demonstrated in studies of human rhinovirus; an entropic basis for the antiviral activity of hydrophobic compounds, predicted from calculated compressibility values, has been corroborated by experimental measurements on poliovirus.     

Robust classification of protein variation using structural modelling and large-scale data integration

Existing methods for interpreting protein variation focus on annotating mutation pathogenicity rather than detailed interpretation of variant deleteriousness and frequently use only sequence-based or structure-based information. We present VIPUR, a computational framework that seamlessly integrates sequence analysis and structural modelling (using the Rosetta protein modelling suite) to identify and interpret deleterious protein variants. To train VIPUR, we collected 9477 protein variants with known effects on protein function from multiple organisms and curated structural models for each variant from crystal structures and homology models. VIPUR can be applied to mutations in any organism''s proteome with improved generalized accuracy (AUROC .83) and interpretability (AUPR .87) compared to other methods. We demonstrate that VIPUR''s predictions of deleteriousness match the biological phenotypes in ClinVar and provide a clear ranking of prediction confidence. We use VIPUR to interpret known mutations associated with inflammation and diabetes, demonstrating the structural diversity of disrupted functional sites and improved interpretation of mutations associated with human diseases. Lastly, we demonstrate VIPUR''s ability to highlight candidate variants associated with human diseases by applying VIPUR to de novo variants associated with autism spectrum disorders.     

Time-resolved resonance Raman study on ultrafast structural relaxation and vibrational cooling of photodissociated carbonmonoxy myoglobin

A localized small structural change is converted to a higher order conformational change of protein and extends to a mesoscopic scale to induce a physiological function. To understand such features of protein, ultrafast dynamics of myoglobin (Mb) following photolysis of carbon monoxide were investigated. Recent results are summarized here with a stress on structural and vibrational energy relaxation. The core expansion of heme takes place within 2 ps but the out of plane displacement of the heme iron and the accompanying protein conformational change occur in 10 and 100 s of the picosecond regimes, respectively. Unexpectedly, it was found from UV resonance Raman spectra that Trp7 in the N-terminal region and Tyr151 in the C-terminal region undergo appreciable structural changes upon ligand binding-dissociation while Tyr104, Tyr146, and Trp14 do not. Because of the communication between the movements of these surface residues and the heme iron, the rate of spectral change of the iron-histidine (Fe- His) stretching band after CO photodissociation is influenced by the viscosity of solvent. Temporal changes of the anti-Stokes Raman intensity demonstrated immediate generation of vibrationally excited heme upon photodissociation and its decay with a time constant of 1-2 ps.     

Molecular Dynamics Simulations of a Membrane Protein/Amphipol Complex

Amphipathic polymers known as “amphipols” provide a highly stabilizing environment for handling membrane proteins in aqueous solutions. A8-35, an amphipol with a polyacrylate backbone and hydrophobic grafts, has been extensively characterized and widely employed for structural and functional studies of membrane proteins using biochemical and biophysical approaches. Given the sensitivity of membrane proteins to their environment, it is important to examine what effects amphipols may have on the structure and dynamics of the proteins they complex. Here we present the first molecular dynamics study of an amphipol-stabilized membrane protein, using Escherichia coli OmpX as a model. We begin by describing the structure of the complexes formed by supplementing OmpX with increasing amounts of A8-35, in order to determine how the amphipol interacts with the transmembrane and extramembrane surfaces of the protein. We then compare the dynamics of the protein in either A8-35, a detergent, or a lipid bilayer. We find that protein dynamics on all accessible length scales is restrained by A8-35, which provides a basis to understanding some of the stabilizing and functional effects of amphipols that have been experimentally observed.     

Hemoprotein time-resolved X-ray crystallography

In the last decade the role of structural dynamics in controlling protein function was actively investigated using new and advanced experimental approaches. In particular, time resolved crystallography, despite some practical difficulties, is being used extensively to complement the study of protein structure-function relationships with information on the dynamics, based on experimental evidence. Here we present a short overview of the results obtained on dynamical properties of myoglobin and homologous hemoproteins, where the photosensitive heme-Fe--ligand bond has allowed transient intermediates to be studied by different flash photolysis methods coupled to Laue X-ray diffraction, thus highlighting some of the dynamical events that characterize diffusion of a diatomic ligand to/from the heme in model hemoproteins.     

Serum Factor protecting against Experimental Arthritis

POLYARTHRITIS induced in rats by the administration of complete Freund's adjuvant is an immunologically mediated disease which can be passively transferred with cells¹. Using this model, we have recently described² the development of a non-complement-dependent antibody specific for lipoprotein antigenic components of Mycobacterium butyricum (Mb). The antibody response precedes both the development of the disease and purified protein derivative (PPD) skin reactivity and is not inhibited by effective doses of immunosuppressive agents that prevent disease and modify PPD reactivity. These observations suggested that the Mb antibody may have a regulatory function in the disease process.     

Picosecond time-resolved X-ray crystallography: probing protein function in real time

A detailed mechanistic understanding of how a protein functions requires knowledge not only of its static structure, but also how its conformation evolves as it executes its function. The recent development of picosecond time-resolved X-ray crystallography has allowed us to visualize in real time and with atomic detail the conformational evolution of a protein. Here, we report the photolysis-induced structural evolution of wild-type and L29F myoglobin over times ranging from 100 ps to 3 micros. The sub-ns structural rearrangements that accompany ligand dissociation in wild-type and the mutant form differ dramatically, and lead to vastly different ligand migration dynamics. The correlated protein displacements provide a structural explanation for the kinetic differences. Our observation of functionally important protein motion on the sub-ns time scale was made possible by the 150-ps time resolution of the measurement, and demonstrates that picosecond dynamics are relevant to protein function. To visualize subtle structural changes without modeling, we developed a novel method for rendering time-resolved electron density that depicts motion as a color gradient across the atom or group of atoms that move. A sequence of these time-resolved images have been stitched together into a movie, which allows one to literally "watch" the protein as it executes its function.