Myoglobin: the hydrogen atom of biology and a paradigm of complexity

Computational infrared and two-dimensional infrared photon echo spectroscopy of both wild-type and double mutant myoglobin-CO proteins

The CO stretching mode of both wild-type and double mutant ( T67R / S92D ) MbCO (carbonmonoxymyoglobin) proteins is an ideal infrared (IR) probe for studying the local electrostatic environment inside the myoglobin heme pocket. Recently, to elucidate the conformational switching dynamics between two distinguishable states, extensive IR absorption, IR pump-probe, and two-dimensional (2D) IR spectroscopic studies for various mutant MbCO's have been performed by the Fayer group. They showed that the 2D IR spectroscopy of the double mutant, which has a peroxidase enzyme activity, reveals a rapid chemical exchange between two distinct states, whereas that of the wild-type does not. Despite the fact that a few simulation studies on these systems were already performed and reported, such complicated experimental results have not been fully reproduced nor described in terms of conformational state-to-state transition processes. Here, we first develop a distributed vibrational solvatochromic charge model for describing the CO stretch frequency shift reflecting local electric potential changes. Then, by carrying out molecular dynamic simulations of the two MbCO's and examining their CO frequency trajectories, it becomes possible to identify a proper reaction coordinate consisting of His64 imidazole ring rotation and its distance to the CO ligand. From the 2D surfaces of the resulting potential of mean forces, the spectroscopically distinguished A1 and A3 states of the wild-type as well as two more substates of the double mutant are identified and their vibrational frequencies and distributions are separately examined. Our simulated IR absorption and 2D IR spectra of the two MbCO's are directly compared with the previous experimental results reported by the Fayer group. The chemical exchange rate constants extracted from the two-state kinetic analyses of the simulated 2D IR spectra are in excellent agreement with the experimental values. On the basis of the quantitative agreement between the simulated spectra and experimental ones, we further examine the conformational differences in the heme pockets of the two proteins and show that the double mutation, T67R / S92D , suppresses the A1 population, restricts the imidazole ring rotation, and increases hydrogen-bond strength between the imidazole Nε-H and the oxygen atom of the CO ligand. It is believed that such delicate change of distal His64 imidazole ring dynamics induced by the double mutation may be responsible for its enhanced peroxidase catalytic activity as compared to the wild-type myoglobin.

Detecting repetitions and periodicities in proteins by tiling the structural space

The notion of energy landscapes provides conceptual tools for understanding the complexities of protein folding and function. Energy landscape theory indicates that it is much easier to find sequences that satisfy the “Principle of Minimal Frustration” when the folded structure is symmetric (Wolynes, P. G. Symmetry and the Energy Landscapes of Biomolecules. Proc. Natl. Acad. Sci. U.S.A.1996, 93, 14249–14255). Similarly, repeats and structural mosaics may be fundamentally related to landscapes with multiple embedded funnels. Here we present analytical tools to detect and compare structural repetitions in protein molecules. By an exhaustive analysis of the distribution of structural repeats using a robust metric, we define those portions of a protein molecule that best describe the overall structure as a tessellation of basic units. The patterns produced by such tessellations provide intuitive representations of the repeating regions and their association toward higher order arrangements. We find that some protein architectures can be described as nearly periodic, while in others clear separations between repetitions exist. Since the method is independent of amino acid sequence information, we can identify structural units that can be encoded by a variety of distinct amino acid sequences.

Ligand migration through hemeprotein cavities: insights from laser flash photolysis and molecular dynamics simulations

The presence of cavities and tunnels in the interior of proteins, in conjunction with the structural plasticity arising from the coupling to the thermal fluctuations of the protein scaffold, has profound consequences on the pathways followed by ligands moving through the protein matrix. In this perspective we discuss how quantitative analysis of experimental rebinding kinetics from laser flash photolysis, trapping of unstable conformational states by embedding proteins within the nanopores of silica gels, and molecular simulations can synergistically converge to gain insight into the migration mechanism of ligands. We show how the evaluation of the free energy landscape for ligand diffusion based on the outcome of computational techniques can assist the definition of sound reaction schemes, leading to a comprehensive understanding of the broad range of chemical events and time scales that encompass the transport of small ligands in hemeproteins.

On the role of thermal backbone fluctuations in myoglobin ligand gate dynamics

We construct an energy function that describes the crystallographic structure of sperm whale myoglobin backbone. As a model in our construction, we use the Protein Data Bank entry 1ABS that has been measured at liquid helium temperature. Consequently, the thermal B-factor fluctuations are very small, which is an advantage in our construction. The energy function that we utilize resembles that of the discrete nonlinear Schrödinger equation. Likewise, ours supports topological solitons as local minimum energy configurations. We describe the 1ABS backbone in terms of topological solitons with a precision that deviates from 1ABS by an average root-mean-square distance, which is less than the experimentally observed Debye-Waller B-factor fluctuation distance. We then subject the topological multi-soliton solution to extensive numerical heating and cooling experiments, over a very wide range of temperatures. We concentrate in particular to temperatures above 300 K and below the Θ-point unfolding temperature, which is around 348 K. We confirm that the behavior of the topological multi-soliton is fully consistent with Anfinsen's thermodynamic principle, up to very high temperatures. We observe that the structure responds to an increase of temperature consistently in a very similar manner. This enables us to characterize the onset of thermally induced conformational changes in terms of three distinct backbone ligand gates. One of the gates is made of the helix F and the helix E. The two other gates are chosen similarly, when open they provide a direct access route for a ligand to reach the heme. We find that out of the three gates we investigate, the one which is formed by helices B and G is the most sensitive to thermally induced conformational changes. Our approach provides a novel perspective to the important problem of ligand entry and exit.

Non-covalent and covalent modifications modulate the reactivity of monomeric mammalian globins

Multimeric globins (e.g., hemoglobin) are considered to be the prototypes of allosteric enzymes, whereas monomeric globins (e.g., myoglobin; Mb) usually are assumed to be non-allosteric. However, the modulation of the functional properties of monomeric globins by non-covalent (or allosteric) and covalent modifications casts doubts on this general assumption. Here, we report examples referable to these two extreme mechanisms modulating the reactivity of three mammalian monomeric globins. Sperm whale Mb, which acts as a reserve supply of O2 and facilitates the O2 flux within a myocyte, displays the allosteric modulation of the O2 affinity on lactate, an obligatory product of glycolysis under anaerobic conditions, thus facilitating O2 diffusion to the mitochondria in supporting oxidative phosphorylation. Human neuroglobin (NGB), which appears to protect neurons from hypoxia in vitro and in vivo, undergoes hypoxia-dependent phosphorylation (i.e., covalent modulation) affecting the coordination equilibrium of the heme-Fe atom and, in turn, the heme-protein reactivity. This facilitates heme-Fe-ligand binding and enhances the rate of anaerobic nitrite reduction to form NO, thus contributing to cellular adaptation to hypoxia. The reactivity of human cytoglobin (CYGB), which has been postulated to protect cells against oxidative stress, depends on both non-covalent and covalent mechanisms. In fact, the heme reactivity of CYGB depends on the lipid, such as oleate, binding which stabilizes the penta-coordination geometry of the heme-Fe atom. Lastly, the reactivity of NGB and CYGB is modulated by the redox state of the intramolecular CysCD7/CysD5 and CysB2/CysE9 residue pairs, respectively, affecting the heme-Fe atom coordination state. In conclusion, the modulation of monomeric globins reactivity by non-covalent and covalent modifications appears a very widespread phenomenon, opening new perspectives in cell survival and protection. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.

Hydrophobic effect drives oxygen uptake in myoglobin via histidine E7

Since the elucidation of the myoglobin (Mb) structure, a histidine residue on the E helix (His-E7) has been proposed to act as a gate with an open or closed conformation controlling access to the active site. Although it is believed that at low pH, the His-E7 gate is in its open conformation, the full relationship between the His-E7 protonation state, its conformation, and ligand migration in Mb is hotly debated. We used molecular dynamics simulations to first address the effect of His-E7 protonation on its conformation. We observed the expected shift from the closed to the open conformation upon protonation, but more importantly, noted a significant difference between the conformations of the two neutral histidine tautomers. We further computed free energy profiles for oxygen migration in each of the possible His-E7 states as well as in two instructive Mb mutants: Ala-E7 and Trp-E7. Our results show that even in the closed conformation, the His-E7 gate does not create a large barrier to oxygen migration and permits oxygen entry with only a small rotation of the imidazole side chain and movement of the E helix. We identify, instead, a hydrophobic site in the E7 channel that can accommodate an apolar diatomic ligand and enhances ligand uptake particularly in the open His-E7 conformation. This rate enhancement is diminished in the closed conformation. Taken together, our results provide a new conceptual framework for the histidine gate hypothesis.

Carbon monoxide recombination dynamics in truncated hemoglobins studied with visible-pump midIR-probe spectroscopy

Carbon monoxide recombination dynamics upon photodissociation with visible light has been characterized by means of ultrafast visible-pump/MidIR probe spectroscopy for the truncated hemoglobins from Thermobifida fusca and Bacillus subtilis. Photodissociation has been induced by exciting the sample at two different wavelengths: 400 nm, corresponding to the heme absorption in the B-band, and 550 nm, in the Q-bands. The bleached iron-CO coordination band located at 1850-1950 cm(-1) and the free CO absorption band in the region 2050-2200 cm(-1) have been observed by probe pulses tuned in the appropriate infrared region. The kinetic traces measured at 1850-1950 cm(-1) reveal multiexponential subnanosecond dynamics that have been interpreted as arising from fast geminate recombination of the photolyzed CO. A compared analysis of the crystal structure of the two proteins reveals a similar structure of their distal heme pocket, which contains conserved polar and aromatic amino acid residues closely interacting with the iron ligand. Although fast geminate recombination is observed in both proteins, several kinetic differences can be evidenced, which can be interpreted in terms of a different structural flexibility of the corresponding heme distal pockets. The analysis of the free CO band-shape and of its dynamic evolution brings out novel features about the nature of the docking site inside the protein cavity.

Quantifying the importance of protein conformation on ligand migration in myoglobin

Myoglobin (Mb) is a model system for ligand binding and migration. The energy barriers (ΔG) for ligand migration in Mb have been studied in the past by experiment and theory and significant differences between different approaches were found. From experiment, it is known that Mb can assume a large number of conformational substates. In this work, these substates are investigated as a possible source of the differences in migration barriers. We show that the initial structure significantly affects the calculated ΔG for a particular transition and that fluctuations in barrier heights δΔG are of similar magnitude as the free energy barriers themselves. The sensitivity of ΔG to the initial structure is compared to other sources of errors. Different protein structures can affect the calculated ΔG by up to 4 kcal/mol, whereas differences between simple point charge models and more elaborate multipolar charge models for the CO-ligand are smaller by a factor of two. Analysis of the structural changes underlying the large effect of the conformational substate reveals the importance of coupling between protein and ligand motion for migration.

Neutron Scattering Reveals Enhanced Protein Dynamics in Concanavalin A Amyloid Fibrils

Protein aggregation is one of the most challenging topics in life sciences, and it is implicated in several human pathologies. The nature and the role of toxic species is highly debated, with amyloid fibrils being among the most relevant species for their peculiar structural and functional properties. Protein dynamics and in particular the ability to fluctuate through a large number of conformational substates are closely related to protein function. This Letter focuses on amyloid fibril dynamics, and, to our knowledge, it is the first neutron scattering study on a protein (Concanavalin A) isolated in its fibril state. Our results reveal enhanced atomic fluctuations in amyloid fibrils and indicate that the protein is "softer" in the fibril state with respect to the native and amorphous aggregate states. We discuss this finding in terms of a structural interpretation and suggest that the paradigm ordered structure ↔ lower flexibility can be questioned when considering the local fast side-chain protein dynamics.

Human serum albumin: from bench to bedside

Human serum albumin (HSA), the most abundant protein in plasma, is a monomeric multi-domain macromolecule, representing the main determinant of plasma oncotic pressure and the main modulator of fluid distribution between body compartments. HSA displays an extraordinary ligand binding capacity, providing a depot and carrier for many endogenous and exogenous compounds. Indeed, HSA represents the main carrier for fatty acids, affects pharmacokinetics of many drugs, provides the metabolic modification of some ligands, renders potential toxins harmless, accounts for most of the anti-oxidant capacity of human plasma, and displays (pseudo-)enzymatic properties. HSA is a valuable biomarker of many diseases, including cancer, rheumatoid arthritis, ischemia, post-menopausal obesity, severe acute graft-versus-host disease, and diseases that need monitoring of the glycemic control. Moreover, HSA is widely used clinically to treat several diseases, including hypovolemia, shock, burns, surgical blood loss, trauma, hemorrhage, cardiopulmonary bypass, acute respiratory distress syndrome, hemodialysis, acute liver failure, chronic liver disease, nutrition support, resuscitation, and hypoalbuminemia. Recently, biotechnological applications of HSA, including implantable biomaterials, surgical adhesives and sealants, biochromatography, ligand trapping, and fusion proteins, have been reported. Here, genetic, biochemical, biomedical, and biotechnological aspects of HSA are reviewed.

Conformational heterogeneity within the Michaelis complex of lactate dehydrogenase

A series of isotope edited IR measurements, both static as well as temperature jump relaxation spectroscopy, are performed on lactate dehydrogenase (LDH) to determine the ensemble of structures available to its Michaelis complex. There clearly has been a substantial reduction in the number of states available to the pyruvate substrate (as modeled by the substrate mimic, oxamate) and NADH when bound to protein compared to dissolved in solution, as determined by the bandwidths and positions of the critical C(2)═O band of the bound substrate mimic and the C(4)-H stretch of the NADH reduced nicotinamide group. Moreover, it is found that a strong ionic bond (characterized by a signature IR band discovered in this study) is formed between the carboxyl group of bound pyruvate with (presumably) Arg171, forming a strong "anchor" within the protein matrix. However, conformational heterogeneity within the Michaelis complex is found that has an impact on both catalytic efficiency and thermodynamics of the enzyme.

Dynamics of a myoglobin mutant enzyme: 2D IR vibrational echo experiments and simulations

Myoglobin (Mb) double mutant T67R/S92D displays peroxidase enzymatic activity in contrast to the wild type protein. The CO adduct of T67R/S92D shows two CO absorption bands corresponding to the A(1) and A(3) substates. The equilibrium protein dynamics for the two distinct substates of the Mb double mutant are investigated by using two-dimensional infrared (2D IR) vibrational echo spectroscopy and molecular dynamics (MD) simulations. The time-dependent changes in the 2D IR vibrational echo line shapes for both of the substates are analyzed using the center line slope (CLS) method to obtain the frequency-frequency correlation function (FFCF). The results for the double mutant are compared to those from the wild type Mb. The experimentally determined FFCF is compared to the FFCF obtained from molecular dynamics simulations, thereby testing the capacity of a force field to determine the amplitudes and time scales of protein structural fluctuations on fast time scales. The results provide insights into the nature of the energy landscape around the free energy minimum of the folded protein structure.

Mechanisms of peroxynitrite interactions with heme proteins

Oxygenated heme proteins are known to react rapidly with nitric oxide (NO) to produce peroxynitrite (PN) at the heme site. This process could lead either to attenuation of the effects of NO or to nitrosative protein damage. PN is a powerful nitrating and oxidizing agent that has been implicated in a variety of cell injuries. Accordingly, it is important to delineate the nature and variety of reaction mechanisms of PN interactions with heme proteins. In this Forum, we survey the range of reactions of PN with heme proteins, with particular attention to myoglobin and cytochrome c. While these two proteins are textbook paradigms for oxygen binding and electron transfer, respectively, both have recently been shown to have other important functions that involve NO and PN. We have recently described direct evidence that ferrylmyolgobin (ferrylMb) and nitrogen dioxide (NO(2)) are both produced during the reaction of PN and metmyolgobin (metMb) (Su, J.; Groves, J. T. J. Am. Chem. Soc. 2009, 131, 12979-12988). Kinetic evidence indicates that these products evolve from the initial formation of a caged radical intermediate [Fe(IV) horizontal lineO.NO(2)]. This caged pair reacts mainly via internal return with a rate constant k(r) to form metMb and nitrate in an oxygen-rebound scenario. Detectable amounts of ferrylMb are observed by stopped-flow spectrophotometry, appearing at a rate consistent with the rate, k(obs), of heme-mediated PN decomposition. Freely diffusing NO(2), which is liberated concomitantly from the radical pair (k(e)), preferentially nitrates myoglobin Tyr103 and added fluorescein. For cytochrome c, Raman spectroscopy has revealed that a substantial fraction of cytochrome c converts to a beta-sheet structure, at the expense of turns and helices at low pH (Balakrishnan, G.; Hu, Y.; Oyerinde, O. F.; Su, J.; Groves, J. T.; Spiro, T. G. J. Am. Chem. Soc., 2007, 129, 504-505). It is proposed that a short beta-sheet segment, comprising residues 37-39 and 58-61, extends itself into the large 37-61 loop when the latter is destabilized by protonation of H26, which forms an anchoring hydrogen bond to loop residue P44. This conformation change ruptures the Met80-Fe bond, as revealed by changes in ligation-sensitive Raman bands. It also induces peroxidase activity with the same temperature profile. This process is suggested to model the apoptotic peroxidation of cardiolipin by cytochrome c.

Protein dynamics and ligand migration interplay as studied by computer simulation

Since proteins are dynamic systems in living organisms, the employment of methodologies contemplating this crucial characteristic results fundamental to allow revealing several aspects of their function. In this work, we present results obtained using classical mechanical atomistic simulation tools applied to understand the connection between protein dynamics and ligand migration. Firstly, we will present a review of the different sampling schemes used in the last years to obtain both ligand migration pathways and the thermodynamic information associated with the process. Secondly, we will focus on representative examples in which the schemes previously presented are employed, concerning the following: i) ligand migration, tunnels, and cavities in myoglobin and neuroglobin; ii) ligand migration in truncated hemoglobin members; iii) NO escape and conformational changes in nitrophorins; iv) ligand selectivity in catalase and hydrogenase; and v) larger ligand migration: the P450 and haloalkane dehalogenase cases. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

Computational study of the resistance shown by the subtype B/HIV-1 protease to currently known inhibitors

Human immunodeficiency virus type 1 protease (HIV-1 PR) is an essential enzyme in the HIV-1 life cycle. As such, this protein represents a major drug target in AIDS therapy, but emerging resistance to antiretroviral inhibitor cocktails, caused by high viral mutation rates, represents a significant challenge in AIDS treatment. Many mutations are not located within the active site or binding pocket, nor they do significantly modify the three-dimensional structural organization of the enzyme; hence, the mechanism(s) by which they alter inhibitor affinity for the protease remains uncertain. In this article, we present an all-atom computational analysis of the dynamic residue-residue coordination between the active site residues and the rest of the protein and of the energetic properties of different HIV-1 PR complexes. We analyze both the wild-type form and mutated forms that induce drug resistance. In particular, the results show differences between the wild type and the mutants in their mechanism of dynamic coordination, in the signal propagation between the active site residues and the rest of the protein, and in the energy networks responsible for the stabilization of the bound inhibitor conformation. Finally, we propose a dynamic and energetic explanation for HIV-1 protease drug resistance, and, through this model, we identify a possible new site that could be helpful in the design of a new family of HIV-1 PR allosteric inhibitors.

Reversibly bound chloride in the atrial natriuretic peptide receptor hormone-binding domain: possible allosteric regulation and a conserved structural motif for the chloride-binding site

The binding of atrial natriuretic peptide (ANP) to its receptor requires chloride, and it is chloride concentration dependent. The extracellular domain (ECD) of the ANP receptor (ANPR) contains a chloride near the ANP-binding site, suggesting a possible regulatory role. The bound chloride, however, is completely buried in the polypeptide fold, and its functional role has remained unclear. Here, we have confirmed that chloride is necessary for ANP binding to the recombinant ECD or the full-length ANPR expressed in CHO cells. ECD without chloride (ECD(-)) did not bind ANP. Its binding activity was fully restored by bromide or chloride addition. A new X-ray structure of the bromide-bound ECD is essentially identical to that of the chloride-bound ECD. Furthermore, bromide atoms are localized at the same positions as chloride atoms both in the apo and in the ANP-bound structures, indicating exchangeable and reversible halide binding. Far-UV CD and thermal unfolding data show that ECD(-) largely retains the native structure. Sedimentation equilibrium in the absence of chloride shows that ECD(-) forms a strongly associated dimer, possibly preventing the structural rearrangement of the two monomers that is necessary for ANP binding. The primary and tertiary structures of the chloride-binding site in ANPR are highly conserved among receptor-guanylate cyclases and metabotropic glutamate receptors. The chloride-dependent ANP binding, reversible chloride binding, and the highly conserved chloride-binding site motif suggest a regulatory role for the receptor bound chloride. Chloride-dependent regulation of ANPR may operate in the kidney, modulating ANP-induced natriuresis.

Myoglobin strikes back

Over the last half century, myoglobin (Mb) has been an excellent model system to test a number of concepts, theories, and new experimental methods that proved valuable to investigate protein structure, function, evolution, and dynamics. Mb's function, most often considered just an oxygen repository, has considerably diversified over the last 15 years, especially because it was shown to have a role in the biochemistry of quenching and synthesizing nitric oxide in the red muscle, thereby protecting the cell. To tackle protein's structural dynamics by innovative biophysical methods, Mb has been the best prototype; laser flash technology made it possible to obtain molecular movies by time-resolved Laue crystallography (with ps resolution). This approach unveiled the complexity of the energy landscape and the structural basis of the stretched interconversion between conformational substates of a protein.

Direct detection of the oxygen rebound intermediates, ferryl Mb and NO2, in the reaction of metmyoglobin with peroxynitrite

Oxygenated hemoproteins are known to react rapidly with nitric oxide (NO) to produce peroxynitrite (PN) at the heme site. This process could lead either to attenuation of the effects of NO or to nitrosative protein damage. Peroxynitrite is a powerful nitrating and oxidizing agent that has been implicated in a variety of cell injuries. Accordingly, it is important to delineate the nature and variety of reaction mechanisms of PN reactions with heme proteins. Here, we present direct evidence that ferrylMb and NO(2) are both produced during the reaction of PN and metmyoglobin (metMb). Kinetic evidence indicates that these products evolve from initial formation of a caged radical intermediate [Fe(IV)=O *NO(2)]. This caged pair reacts mainly via internal return with a rate constant k(r) to form metMb and nitrate in an oxygen rebound scenario. Detectable amounts of ferrylMb are observed by stopped-flow spectrophotometry, appearing at a rate consistent with the rate, k(obs), of heme-mediated PN decomposition. Freely diffusing NO(2), which is liberated concomitantly from the radical pair (k(e)), preferentially nitrates Tyr103 in horse heart myoglobin. The ratio of the rates of in-cage rebound and cage escape, k(r)/k(e), was found to be approximately 10 by examining the nitration yields of fluorescein, an external NO(2) trap. This rebound/escape model for the metMb/PN interaction is analogous to the behavior of alkyl hyponitrites and the well-studied geminate recombination processes of deoxymyoglobin with O(2), CO, and NO. The scenario is also similar to the stepwise events of substrate hydroxylation by cytochrome P450 and other oxygenases. It is likely, therefore, that the reaction of metMb with ONOO(-) and that of oxyMb with NO proceed through the same [Fe(IV)=O *NO(2)] caged radical intermediate and lead to similar outcomes. The results indicate that while oxyMb may reduce the concentration of intracellular NO, it would not eliminate the formation of NO(2) as a decomposition product of peroxynitrite.

Osmolyte-induced perturbations of hydrogen bonding between hydration layer waters: correlation with protein conformational changes

Gadolinium vibronic sideband luminescence spectroscopy (GVSBLS) is used to probe osmolyte-induced changes in the hydrogen bond strength between first and second shell waters on the surface of free Gd(3+) and Gd(3+) coordinated to EDTA and to structured calcium binding peptides in solution. In parallel, Raman is used to probe the corresponding impact of the same set of osmolytes on hydrogen bonding among waters in the bulk phase. Increasing concentration of added urea is observed to progressively weaken the hydrogen bonding within the hydration layer but has minimal observed impact on bulk water. In contrast, polyols are observed to enhance hydrogen bonding in both the hydration layer and the bulk with the amplitude being polyol dependent with trehalose being more effective than sucrose, glucose, or glycerol. The observed patterns indicate that the size and properties of the osmolyte as well as the local architecture of the specific surface site of hydration impact preferential exclusion effects and local hydrogen bond strength. Correlation of the vibronic spectra with CD measurements on the peptides as a function of added osmolytes shows an increase in secondary structure with added polyols and that the progressive weakening of the hydrogen bonding upon addition of urea first increases water occupancy within the peptide and only subsequently does the peptide unfold. The results support models in which the initial steps in the unfolding process involve osmolyte-induced enhancement of water occupancy within the interior of the protein.

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.

Charge density-dependent modifications of hydration shell waters by Hofmeister ions

Gadolinium (Gd(3+)) vibronic sideband luminescence spectroscopy (GVSBLS) is used to probe, as a function of added Hofmeister series salts, changes in the OH stretching frequency derived from first-shell waters of aqueous Gd(3+) and of Gd(3+) coordinated to three different types of molecules: (i) a chelate (EDTA), (ii) structured peptides (mSE3/SE2) of the lanthanide-binding tags (LBTs) family with a single high-affinity binding site, and (iii) a calcium-binding protein (calmodulin) with four binding sites. The vibronic sideband (VSB) corresponding to the OH stretching mode of waters coordinated to Gd(3+), whose frequency is inversely correlated with the strength of the hydrogen bonding to neighboring waters, exhibits an increase in frequency when Gd(3+) becomes coordinated to either EDTA, calmodulin, or mSE3 peptide. In all of these cases, the addition of cation chloride or acetate salts to the solution increases the frequency of the vibronic band originating from the OH stretching mode of the coordinated waters in a cation- and concentration-dependent fashion. The cation dependence of the frequency increase scales with charge density of the cations, giving rise to an ordering consistent with the Hofmeister ordering. On the other hand, water Raman spectroscopy shows no significant change upon addition of these salts. Additionally, it is shown that the cation effect is modulated by the specific anion used. The results indicate a mechanism of action for Hofmeister series ions in which hydrogen bonding among hydration shell waters is modulated by several factors. High charge density cations sequester waters in a configuration that precludes strong hydrogen bonding to neighboring waters. Under such conditions, anion effects emerge as anions compete for hydrogen-bonding sites with the remaining free waters on the surface of the hydration shell. The magnitude of the anion effect is both cation and Gd(3+)-binding site specific.

Time-resolved flow-flash FT-IR difference spectroscopy: the kinetics of CO photodissociation from myoglobin revisited

Fourier-transform infrared (FT-IR) difference spectroscopy has been proven to be a significant tool in biospectroscopy. In particular, the step-scan technique monitors structural and electronic changes at time resolutions down to a few nanoseconds retaining the multiplex advantage of FT-IR. For the elucidation of the functional mechanisms of proteins, this technique is currently limited to repetitive systems undergoing a rapid photocycle. To overcome this obstacle, we developed a flow-flash experiment in a miniaturised flow channel which was integrated into a step-scan FT-IR spectroscopic setup. As a proof of principle, we studied the rebinding reaction of CO to myoglobin after photodissociation. The use of microfluidics reduced the sample consumption drastically such that a typical step-scan experiment takes only a few 10 ml of a millimolar sample solution, making this method particularly interesting for the investigation of biological samples that are only available in small quantities. Moreover, the flow cell provides the unique opportunity to assess the reaction mechanism of proteins that cycle slowly or react irreversibly. We infer that this novel approach will help in the elucidation of molecular reactions as complex as those of vectorial ion transfer in membrane proteins. The potential application to the oxygen splitting reaction of cytochrome c oxidase is discussed.

Design and engineering of an O(2) transport protein

The principles of natural protein engineering are obscured by overlapping functions and complexity accumulated through natural selection and evolution. Completely artificial proteins offer a clean slate on which to define and test these protein engineering principles, while recreating and extending natural functions. Here we introduce this method with the design of an oxygen transport protein, akin to human neuroglobin. Beginning with a simple and unnatural helix-forming sequence with just three different amino acids, we assembled a four-helix bundle, positioned histidines to bis-histidine ligate haems, and exploited helical rotation and glutamate burial on haem binding to introduce distal histidine strain and facilitate O(2) binding. For stable oxygen binding without haem oxidation, water is excluded by simple packing of the protein interior and loops that reduce helical-interface mobility. O(2) affinities and exchange timescales match natural globins with distal histidines, with the remarkable exception that O(2) binds tighter than CO.

An integrated model for enzyme catalysis emerges from studies of hydrogen tunneling

The origins of the enormous rate accelerations brought about by enzymes are discussed. The focus is on enzymatic C-H activation, which has been shown to take place via tunneling. Four enzyme systems illustrate the impact of site-specific mutagenesis, changes in temperature or changes in protein solvation on the tunneling properties. A model emerges in which conformational sampling is required to access a subset of protein conformers where the H-donor and acceptor undergo a close approach. The evidence for an inverse relationship between protein flexibility and active site compression is likely to extend to all classes of enzyme catalysts.

Modeling signal propagation mechanisms and ligand-based conformational dynamics of the Hsp90 molecular chaperone full-length dimer

Hsp90 is a molecular chaperone essential for protein folding and activation in normal homeostasis and stress response. ATP binding and hydrolysis facilitate Hsp90 conformational changes required for client activation. Hsp90 plays an important role in disease states, particularly in cancer, where chaperoning of the mutated and overexpressed oncoproteins is important for function. Recent studies have illuminated mechanisms related to the chaperone function. However, an atomic resolution view of Hsp90 conformational dynamics, determined by the presence of different binding partners, is critical to define communication pathways between remote residues in different domains intimately affecting the chaperone cycle. Here, we present a computational analysis of signal propagation and long-range communication pathways in Hsp90. We carried out molecular dynamics simulations of the full-length Hsp90 dimer, combined with essential dynamics, correlation analysis, and a signal propagation model. All-atom MD simulations with timescales of 70 ns have been performed for complexes with the natural substrates ATP and ADP and for the unliganded dimer. We elucidate the mechanisms of signal propagation and determine “hot spots” involved in interdomain communication pathways from the nucleotide-binding site to the C-terminal domain interface. A comprehensive computational analysis of the Hsp90 communication pathways and dynamics at atomic resolution has revealed the role of the nucleotide in effecting conformational changes, elucidating the mechanisms of signal propagation. Functionally important residues and secondary structure elements emerge as effective mediators of communication between the nucleotide-binding site and the C-terminal interface. Furthermore, we show that specific interdomain signal propagation pathways may be activated as a function of the ligand. Our results support a “conformational selection model” of the Hsp90 mechanism, whereby the protein may exist in a dynamic equilibrium between different conformational states available on the energy landscape and binding of a specific partner can bias the equilibrium toward functionally relevant complexes.

Dynamic processes underlie the functions of all proteins. Hence, to understand, control, and design protein functions in the cell, we need to unravel the basic principles of protein dynamics. This is fundamental in studying the mechanisms of a specific class of proteins known as molecular chaperones, which oversee the correct conformational maturation of other proteins. In particular, molecular chaperones of the stress response machinery have become the focus of intense research, because their upregulation is responsible for the ability of tumor cells to cope with unfavorable environments. This is largely centered on the expression and function of the molecular chaperone Hsp90, which has provided an attractive target for therapeutic intervention in cancer. Experiments have shown that the chaperone functions through a nucleotide-directed conformational cycle. Here, we show that it is possible to identify the effects of nucleotide-related chemical differences on functionally relevant motions at the atomic level of resolution. The protein may fluctuate at equilibrium among different available dynamic states, and binding of a specific partner may shift the equilibrium toward the thermodynamically most stable complexes. These results provide us with important mechanistic insight for the identification of new regulatory sites and the design of possible new drugs.

Free-energy landscape, principal component analysis, and structural clustering to identify representative conformations from molecular dynamics simulations: the myoglobin case

Several molecular dynamics (MD) simulations were used to sample conformations in the neighborhood of the native structure of holo-myoglobin (holo-Mb), collecting trajectories spanning 0.22 micros at 300 K. Principal component (PCA) and free-energy landscape (FEL) analyses, integrated by cluster analysis, which was performed considering the position and structures of the individual helices of the globin fold, were carried out. The coherence between the different structural clusters and the basins of the FEL, together with the convergence of parameters derived by PCA indicates that an accurate description of the Mb conformational space around the native state was achieved by multiple MD trajectories spanning at least 0.14 micros. The integration of FEL, PCA, and structural clustering was shown to be a very useful approach to gain an overall view of the conformational landscape accessible to a protein and to identify representative protein substates. This method could be also used to investigate the conformational and dynamical properties of Mb apo-, mutant, or delete versions, in which greater conformational variability is expected and, therefore identification of representative substates from the simulations is relevant to disclose structure-function relationship.

Variations in clique and community patterns in protein structures during allosteric communication: investigation of dynamically equilibrated structures of methionyl tRNA synthetase complexes

The allosteric concept has played a key role in understanding the biological functions of proteins. The rigidity or plasticity and the conformational population are the two important ideas invoked in explaining the allosteric effect. Although molecular insights have been gained from a large number of structures, a precise assessment of the ligand-induced conformational changes in proteins at different levels, ranging from gross topology to intricate details, remains a challenge. In this study, we have explored the conformational changes in the complexes of methionyl tRNA synthetase (MetRS) through novel network parameters such as cliques and communities, which identify the rigid regions in the protein structure networks (PSNs) constructed from the noncovalent interactions of amino acid side chains. MetRS belongs to the aminoacyl tRNA synthetase (aaRS) family that plays a crucial role in the translation of genetic code. These enzymes are modular with distinct domains from which extensive genetic, kinetic, and structural data are available, highlighting the role of interdomain communication. The network parameters evaluated here on the conformational ensembles of MetRS complexes, generated from molecular dynamics simulations, have enabled us to understand the interdomain communication in detail. Additionally, the characterization of conformational changes in terms of cliques and communities has also become possible, which had eluded conventional analyses. Furthermore, we find that most of the residues participating in cliques and communities are strikingly different from those that take part in long-range communication. The cliques and communities evaluated here for the first time on PSNs have beautifully captured the local geometries in detail within the framework of global topology. Here the allosteric effect is revealed at the residue level via identification of the important residues specific for structural rigidity and functional flexibility in MetRS. This ought to enhance our understanding of the functioning of aaRS in general.

Atomic level computational identification of ligand migration pathways between solvent and binding site in myoglobin

Myoglobin is a globular protein involved in oxygen storage and transport. No consensus yet exists on the atomic level mechanism by which oxygen and other small nonpolar ligands move between the myoglobin's buried heme, which is the ligand binding site, and surrounding solvent. This study uses room temperature molecular dynamics simulations to provide a complete atomic level picture of ligand migration in myoglobin. Multiple trajectories--providing a cumulative total of 7 micros of simulation--are analyzed. Our simulation results are consistent with and tie together previous experimental findings. Specifically, we characterize: (i) Explicit full trajectories in which the CO ligand shuttles between the internal binding site and the solvent and (ii) pattern and structural origins of transient voids available for ligand migration. The computations are performed both in sperm whale myoglobin wild-type and in sperm whale V68F myoglobin mutant, which is experimentally known to slow ligand-binding kinetics. On the basis of these independent, but mutually consistent ligand migration and transient void computations, we find that there are two discrete dynamical pathways for ligand migration in myoglobin. Trajectory hops between these pathways are limited to two bottleneck regions. Ligand enters and exits the protein matrix in common identifiable portals on the protein surface. The pathways are located in the "softer" regions of the protein matrix and go between its helices and in its loop regions. Localized structural fluctuations are the primary physical origin of the simulated CO migration pathways inside the protein.

Crystallographic studies of Xe and Kr binding within the large internal cavity of cytochrome ba3 from Thermus thermophilus: structural analysis and role of oxygen transport channels in the heme-Cu oxidases

Cytochrome ba3 is a cytochrome c oxidase from the plasma membrane of Thermus thermophilus and is the preferred terminal enzyme of cellular respiration at low dioxygen tensions. Using cytochrome ba 3 crystals pressurized at varying conditions under Xe or Kr gas, and X-ray data for six crystals, we identify the relative affinities of Xe and Kr atoms for as many as seven distinct binding sites. These sites track a continuous, Y-shaped channel, 18-20 A in length, lined by hydrophobic residues, which leads from the surface of the protein where two entrance holes, representing the top of the Y, connect the bilayer to the a3-CuB center at the base of the Y. Considering the increased affinity of O2 for hydrophobic environments, the hydrophobic nature of the channel, its orientation within the bilayer, its connection to the active site, its uniform diameter, its virtually complete occupation by Xe, and its isomorphous presence in the native enzyme, we infer that the channel is a diffusion pathway for O2 into the dinuclear center of cytochrome ba3. These observations provide a basis for analyzing similar channels in other oxidases of known structure, and these structures are discussed in terms of mechanisms of O2 transport in biological systems, details of CO binding to and egress from the dinuclear center, the bifurcation of the oxygen-in and water-out pathways, and the possible role of the oxygen channel in aerobic thermophily.