Ligand diffusion in globins: simulations versus experiment

Vibrational analyses of the reaction of oxymyoglobin with NO using a photolabile caged NO donor at cryogenic temperatures

The NO dioxygenation reaction catalyzed by heme-containing globin proteins is a crucial aerobic detoxification pathway. Accordingly, the second order reaction of NO with oxymyoglobin and oxyhemoglobin has been the focus of a large number of kinetic and spectroscopic studies. Stopped-flow and rapid-freeze-quench (RFQ) measurements have provided evidence for the formation of a Fe(III)-nitrato complex with millisecond lifetime prior to release of the nitrate product, but the temporal resolution of these techniques is insufficient for the characterization of precursor species. Most mechanistic models assume the formation of an initial Fe(III)-peroxynitrite species prior to homolytic cleavage of the OO bond and recombination of the resulting NO2 and Fe(IV)=O species. Here we report vibrational spectroscopy measurements for the reaction of oxymyoglobin with a photolabile caged NO donor at cryogenic temperatures. We show that this approach offers efficient formation and trapping of the Fe(III)-nitrato, enzyme-product, complex at 180 K. Resonance Raman spectra of the Fe(III)-nitrato complex trapped via RFQ in the liquid phase and photolabile NO release at cryogenic temperatures are indistinguishable, demonstrating the complementarity of these approaches. Caged NO is released by irradiation <180 K but diffusion into the heme pocket is fully inhibited. Our data provide no evidence for Fe(III)-peroxynitrite of Fe(IV)=O species, supporting low activation energies for the NO to nitrate conversion at the oxymyoglobin reaction site. Photorelease of NO at cryogenic temperatures allows monitoring of the reaction by transmittance FTIR which provides valuable quantitative information and promising prospects for the detection of protein sidechain reorganization events in NO-reacting metalloenzymes.

Ligand-Based Regulation of Dynamics and Reactivity of Hemoproteins

Hemoproteins include several heme-binding proteins with distinct structure and function. The presence of the heme group confers specific reactivity and spectroscopic properties to hemoproteins. In this review, we provide an overview of five families of hemoproteins in terms of dynamics and reactivity. First, we describe how ligands modulate cooperativity and reactivity in globins, such as myoglobin and hemoglobin. Second, we move on to another family of hemoproteins devoted to electron transport, such as cytochromes. Later, we consider heme-based reactivity in hemopexin, the main heme-scavenging protein. Then, we focus on heme-albumin, a chronosteric hemoprotein with peculiar spectroscopic and enzymatic properties. Eventually, we analyze the reactivity and dynamics of the most recently discovered family of hemoproteins, i.e., nitrobindins.

Computational Investigation of a Series of Small Molecules as Potential Compounds for Lysyl Hydroxylase-2 (LH2) Inhibition

The catalytic function of lysyl hydroxylase-2 (LH2), a member of the Fe(II)/αKG-dependent oxygenase superfamily, is to catalyze the hydroxylation of lysine to hydroxylysine in collagen, resulting in stable hydroxylysine aldehyde-derived collagen cross-links (HLCCs). Reports show that high amounts of LH2 lead to the accumulation of HLCCs, causing fibrosis and specific types of cancer metastasis. Some members of the Fe(II)/αKG-dependent family have also been reported to have intramolecular O2 tunnels, which aid in transporting one of the required cosubstrates into the active site. While LH2 can be a promising target to combat these diseases, efficacious inhibitors are still lacking. We have used computational simulations to investigate a series of 44 small molecules as lead compounds for LH2 inhibition. Tunneling analyses indicate the existence of several intramolecular tunnels. The lengths of the calculated O2-transporting tunnels in holoenzymes are relatively longer than those in the apoenzyme, suggesting that the ligands may affect the enzyme's structure and possibly block (at least partially) the tunnels. The sequence alignment analysis between LH enzymes from different organisms shows that all of the amino acid residues with the highest occurrence rate in the oxygen tunnels are conserved. Our results suggest that the enolate form of diketone compounds establishes stronger interactions with the Fe(II) in the active site. Branching the enolate compounds with functional groups such as phenyl and pyridinyl enhances the interaction with various residues around the active site. Our results provide information about possible leads for further LH2 inhibition design and development.

Structural and (Pseudo-)Enzymatic Properties of Neuroglobin: Its Possible Role in Neuroprotection

Neuroglobin (Ngb), the third member of the globin family, was discovered in human and murine brains in 2000. This monomeric globin is structurally similar to myoglobin (Mb) and hemoglobin (Hb) α and β subunits, but it hosts a bis-histidyl six-coordinated heme-Fe atom. Therefore, the heme-based reactivity of Ngb is modulated by the dissociation of the distal HisE7-heme-Fe bond, which reflects in turn the redox state of the cell. The high Ngb levels (~100–200 μM) present in the retinal ganglion cell layer and in the optic nerve facilitate the O₂ buffer and delivery. In contrast, the very low levels of Ngb (~1 μM) in most tissues and organs support (pseudo-)enzymatic properties including NO/O₂ metabolism, peroxynitrite and free radical scavenging, nitrite, hydroxylamine, hydrogen sulfide reduction, and the nitration of aromatic compounds. Here, structural and (pseudo-)enzymatic properties of Ngb, which are at the root of tissue and organ protection, are reviewed, envisaging a possible role in the protection from neuronal degeneration of the retina and the optic nerve.

Kinetic mechanisms for O2 binding to myoglobins and hemoglobins

Antonini and Brunori's 1971 book "Hemoglobin and Myoglobin in Their Reactions with Ligands" was a truly remarkable publication that summarized almost 100 years of research on O₂ binding to these globins. Over the ensuing 50 years, ultra-fast laser photolysis techniques, high-resolution and time resolved X-ray crystallography, molecular dynamics simulations, and libraries of recombinant myoglobin (Mb) and hemoglobin (Hb) variants have provided structural interpretations of O₂ binding to these proteins. The resultant mechanisms provide quantitative descriptions of the stereochemical factors that govern overall affinity, including proximal and distal steric restrictions that affect iron reactivity and favorable positive electrostatic interactions that preferentially stabilize bound O₂. The pathway for O₂ uptake and release by Mb and subunits of Hb has been mapped by screening libraries of site-directed mutants in laser photolysis experiments. O₂ enters mammalian Mb and the α and β subunits of human HbA through a channel created by upward and outward rotation of the distal His at the E7 helical position, is non-covalently captured in the interior of the distal cavity, and then internally forms a bond with the heme Fe(II) atom. O₂ dissociation is governed by disruption of hydrogen bonding interactions with His (E7), breakage of the Fe(II)-O₂ bond, and then competition between rebinding and escape through the E7-gate. The structural features that govern the rates of both the individual steps and overall reactions have been determined and provide the framework for: (1) defining the physiological functions of specific globins and their evolution; (2) understanding the clinical features of hemoglobinopathies; and (3) designing safer and more efficient acellular hemoglobin-based oxygen carriers (HBOCs) for transfusion therapy, organ preservation, and other commercially relevant O₂ transport and storage processes.

Lessons Learned from 50 Years of Hemoglobin Research: Unstirred and Cell-Free Layers, Electrostatics, Baseball Gloves, and Molten Globules

Significance: Over the past 50 years, the mechanisms for O₂ storage and transport have been determined quantitatively on distance scales from millimeters to tenths of nanometers and timescales from seconds to picoseconds. Recent Advances: In this review, I have described four key conclusions from work done by my group and our close colleagues. (i) O₂ uptake by mammalian red cells is limited by diffusion through unstirred water layers adjacent to the cell surface and across cell-free layers adjacent to vessel walls. (ii) In most vertebrates, hemoglobins (Hbs) and myoglobins (Mbs), the distal histidine at the E7 helical position donates a strong hydrogen bond to bound O₂, which selectively enhances O₂ affinity, prevents carbon monoxide poisoning, and markedly slows autoxidation. (iii) O₂ binding to mammalian Hbs and Mbs occurs by migration of the ligand through a channel created by upward rotation of the His(E7) side chain, capture in the empty space of the distal pocket, and then coordination with the ferroprotoporphyrin IX (heme) iron atom. (iv) The assembly of Mbs and Hbs occurs by formation of molten globule intermediates, in which the N- and C-terminal helices have almost fully formed secondary structures, but the heme pockets are disordered and followed by high-affinity binding of heme. Critical Issues: These conclusions indicate that there are often compromises between O₂ transport function, holoprotein stability, and the efficiency of assembly. Future Directions: However, the biochemical mechanisms underlying these conclusions provide the framework for understanding globin evolution in greater detail and for engineering more efficient and stable globins.

Hemoglobin: Structure, Function and Allostery

This chapter reviews how allosteric (heterotrophic) effectors and natural mutations impact hemoglobin (Hb) primary physiological function of oxygen binding and transport. First, an introduction about the structure of Hb is provided, including the ensemble of tense and relaxed Hb states and the dynamic equilibrium of Hb multistate. This is followed by a brief review of Hb variants with altered Hb structure and oxygen binding properties. Finally, a review of different endogenous and exogenous allosteric effectors of Hb is presented with particular emphasis on the atomic interactions of synthetic ligands with altered allosteric function of Hb that could potentially be harnessed for the treatment of diseases.

Dynamical comparison between myoglobin and hemoglobin

Myoglobin and hemoglobin are globular hemeproteins, when the former is a monomer and the latter a heterotetramer. Despite the structural similarity of myoglobin to α and β subunits of hemoglobin, there is a functional difference between the two proteins, owing to the quaternary structure of hemoglobin. The effect of the quaternary structure of hemoglobin on the intrinsic dynamics of its subunits is explored by dynamical comparison of the two proteins. Anisotropic Network Model modes of motion were calculated for hemoglobin and myoglobin. Dynamical comparison between the proteins was performed using global and local Anisotropic Network Model mode alignment algorithms based on the algorithms of Smith-Waterman and Needleman-Wunsch for sequence comparison. The results indicate that the quaternary structure of hemoglobin substantially alters the intrinsic dynamics of its subunits, an effect that may contribute to the functional difference between the two proteins. Local dynamics similarity between the proteins is still observed at the major exit route of the ligand.

Atomic resolution mechanism of ligand binding to a solvent inaccessible cavity in T4 lysozyme

Ligand binding sites in proteins are often localized to deeply buried cavities, inaccessible to bulk solvent. Yet, in many cases binding of cognate ligands occurs rapidly. An intriguing system is presented by the L99A cavity mutant of T4 Lysozyme (T4L L99A) that rapidly binds benzene (~106 M-1s-1). Although the protein has long served as a model system for protein thermodynamics and crystal structures of both free and benzene-bound T4L L99A are available, the kinetic pathways by which benzene reaches its solvent-inaccessible binding cavity remain elusive. The current work, using extensive molecular dynamics simulation, achieves this by capturing the complete process of spontaneous recognition of benzene by T4L L99A at atomistic resolution. A series of multi-microsecond unbiased molecular dynamics simulation trajectories unequivocally reveal how benzene, starting in bulk solvent, diffuses to the protein and spontaneously reaches the solvent inaccessible cavity of T4L L99A. The simulated and high-resolution X-ray derived bound structures are in excellent agreement. A robust four-state Markov model, developed using cumulative 60 μs trajectories, identifies and quantifies multiple ligand binding pathways with low activation barriers. Interestingly, none of these identified binding pathways required large conformational changes for ligand access to the buried cavity. Rather, these involve transient but crucial opening of a channel to the cavity via subtle displacements in the positions of key helices (helix4/helix6, helix7/helix9) leading to rapid binding. Free energy simulations further elucidate that these channel-opening events would have been unfavorable in wild type T4L. Taken together and via integrating with results from experiments, these simulations provide unprecedented mechanistic insights into the complete ligand recognition process in a buried cavity. By illustrating the power of subtle helix movements in opening up multiple pathways for ligand access, this work offers an alternate view of ligand recognition in a solvent-inaccessible cavity, contrary to the common perception of a single dominant pathway for ligand binding.

Identification of the binding site between bovine serum albumin and ultrasmall SiC fluorescent biomarkers

Combined experimental and theoretical studies propose a delayed diffusion model for describing the interaction between ultrasmall NPs and proteins.

Computational investigation of O2 diffusion through an intra-molecular tunnel in AlkB; influence of polarization on O2 transport

E. Coli AlkB catalyzes the direct dealkylation of various alkylated bases in damaged DNA. The diffusion of molecular oxygen to the active site in AlkB is an essential step for the oxidative dealkylation activity. Despite detailed studies on the stepwise oxidation mechanism of AlkB, there is no conclusive picture of how O2 molecules reach the active site of the protein. Yu et al. (Nature, 439, 879) proposed the existence of an intra-molecular tunnel based on their initial crystal structures of AlkB. We have employed computational simulations to investigate possible migration pathways inside AlkB for O2 molecules. Extensive molecular dynamics (MD) simulations, including explicit ligand sampling and potential of mean force (PMF) calculations, have been performed to provide a microscopic description of the O2 delivery pathway in AlkB. Analysis of intra-molecular tunnels using the CAVER software indicates two possible pathways for O2 to diffuse into the AlkB active site. Explicit ligand sampling simulations suggests that only one of these tunnels provides a viable route. The free energy path for an oxygen molecule to travel along each of these tunnels has been determined with AMBER and AMOEBA. Both PMFs indicate passive transport of O2 from the surface of the protein. However, the inclusion of explicit polarization shows a very large barrier for diffusion of the co-substrate out of the active site, compared with the non-polarizable potential. In addition, our results suggest that the mutation of a conserved residue along the tunnel, Y178, has dramatic effects on the dynamics of AlkB and on the transport of O2 along the tunnel.

Thermodynamics of camphor migration in cytochrome P450cam by atomistic simulations

Understanding the mechanisms of ligand binding to enzymes is of paramount importance for the design of new drugs. Here, we report on the use of a novel biased molecular dynamics (MD) methodology to study the mechanism of camphor binding to cytochrome P450cam. Microsecond-long MD simulations allowed us to observe reaction coordinates characterizing ligand diffusion from the active site of cytochrome P450cam to solvent via three egress routes. These atomistic simulations were used to estimate thermodynamic quantities along the reaction coordinates and indicate diverse binding configurations. The results suggest that the diffusion of camphor along the pathway near the substrate recognition site (SRS) is thermodynamically preferred. In addition, we show that the diffusion near the SRS is triggered by a transition from a heterogeneous collection of closed ligand-bound conformers to the basin comprising the open conformations of cytochrome P450cam. The conformational change accompanying this switch is characterized by the retraction of the F and G helices and the disorder of the B' helix. These results are corroborated by experimental studies and provide detailed insight into ligand binding and conformational behavior of the cytochrome family. The presented methodology is general and can be applied to other ligand-protein systems.

Privileged Electrophile Sensors: A Resource for Covalent Drug Development

This Perspective delineates how redox signaling affects the activity of specific enzyme isoforms and how this property may be harnessed for rational drug design. Covalent drugs have resurged in recent years and several reports have extolled the general virtues of developing irreversible inhibitors. Indeed, many modern pharmaceuticals contain electrophilic appendages. Several invoke a warhead that hijacks active-site nucleophiles whereas others take advantage of spectator nucleophilic side chains that do not participate in enzymatic chemistry, but are poised to bind/react with electrophiles. The latest data suggest that innate electrophile sensing-which enables rapid reaction with an endogenous signaling electrophile-is a quintessential resource for the development of covalent drugs. For instance, based on recent work documenting isoform-specific electrophile sensing, isozyme non-specific drugs may be converted to isozyme-specific analogs by hijacking privileged first-responder electrophile-sensing cysteines. Because this approach targets functionally relevant cysteines, we can simultaneously harness previously untapped moonlighting roles of enzymes linked to redox sensing.

Ligand diffusion in proteins via enhanced sampling in molecular dynamics

Computational simulations in biophysics describe the dynamics and functions of biological macromolecules at the atomic level. Among motions particularly important for life are the transport processes in heterogeneous media. The process of ligand diffusion inside proteins is an example of a complex rare event that can be modeled using molecular dynamics simulations. The study of physical interactions between a ligand and its biological target is of paramount importance for the design of novel drugs and enzymes. Unfortunately, the process of ligand diffusion is difficult to study experimentally. The need for identifying the ligand egress pathways and understanding how ligands migrate through protein tunnels has spurred the development of several methodological approaches to this problem. The complex topology of protein channels and the transient nature of the ligand passage pose difficulties in the modeling of the ligand entry/escape pathways by canonical molecular dynamics simulations. In this review, we report a methodology involving a reconstruction of the ligand diffusion reaction coordinates and the free-energy profiles along these reaction coordinates using enhanced sampling of conformational space. We illustrate the above methods on several ligand-protein systems, including cytochromes and G-protein-coupled receptors. The methods are general and may be adopted to other transport processes in living matter.

Temperature Sensitive Singlet Oxygen Photosensitization by LOV-Derived Fluorescent Flavoproteins

Optogenetic sensitizers that selectively produce a given reactive oxygen species (ROS) constitute a promising tool for studying cell signaling processes with high levels of spatiotemporal control. However, to harness the full potential of this tool for live cell studies, the photophysics of currently available systems need to be explored further and optimized. Of particular interest in this regard, are the flavoproteins miniSOG and SOPP, both of which (1) contain the chromophore flavin mononucleotide, FMN, in a LOV-derived protein enclosure, and (2) photosensitize the production of singlet oxygen, O₂(a¹Δ_g). Here we present an extensive experimental study of the singlet and triplet state photophysics of FMN in SOPP and miniSOG over a physiologically relevant temperature range. Although changes in temperature only affect the singlet excited state photophysics slightly, the processes that influence the deactivation of the triplet excited state are more sensitive to temperature. Most notably, for both proteins, the rate constant for quenching of ³FMN by ground state oxygen, O₂(X³Σ_g^-), increases ∼10-fold upon increasing the temperature from 10 to 43 °C, while the oxygen-independent channels of triplet state deactivation are less affected. As a consequence, this increase in temperature results in higher yields of O₂(a¹Δ_g) formation for both SOPP and miniSOG. We also show that the quantum yields of O₂(a¹Δ_g) production by both miniSOG and SOPP are mainly limited by the fraction of FMN triplet states quenched by O₂(X³Σ_g^-). The results presented herein provide a much-needed quantitative framework that will facilitate the future development of optogenetic ROS sensitizers.

Exploration of gated ligand binding recognizes an allosteric site for blocking FABP4-protein interaction

Fatty acid binding protein 4 (FABP4), reversibly binding to fatty acids and other lipids with high affinities, is a potential target for treatment of cancers. The binding site of FABP4 is buried in an interior cavity and thereby ligand binding/unbinding is coupled with opening/closing of FABP4. It is a difficult task both experimentally and computationally to illuminate the entry or exit pathway, especially with the conformational gating. In this report we combine extensive computer simulations, clustering analysis, and the Markov state model to investigate the binding mechanism of FABP4 and troglitazone. Our simulations capture spontaneous binding and unbinding events as well as the conformational transition of FABP4 between the open and closed states. An allosteric binding site on the protein surface is recognized for the development of novel FABP4 inhibitors. The binding affinity is calculated and compared with the experimental value. The kinetic analysis suggests that ligand residence on the protein surface may delay the binding process. Overall, our results provide a comprehensive picture of ligand diffusion on the protein surface, ligand migration into the buried cavity, and the conformational change of FABP4 at an atomic level.

A molecule for all seasons: The heme

If life without heme-Fe were at all possible, it would definitely be different. Indeed this complex and versatile iron-porphyrin macrocycle upon binding to different “globins” yields hemeproteins crucial to sustain a variety of vital functions, generally classified, for convenience, in a limited number of functional families. Over-and-above the array of functions briefly outlined below, the spectacular progress in molecular genetics seen over the last 30 years led to the discovery of many hitherto unknown novel hemeproteins in prokaryotes and eukaryotes. Here, we highlight a few basic aspects of the chemistry of the hemeprotein universe, in particular those that are relevant to the control of heme-Fe reactivity and specialization, as sculpted by a variety of interactions with the protein moiety.

Benchmarking Rapid TLES Simulations of Gas Diffusion in Proteins: Mapping O2 Migration and Escape in Myoglobin as a Case Study

Standard molecular dynamics (MD) simulations of gas diffusion consume considerable computational time and resources even for small proteins. To combat this, temperature-controlled locally enhanced sampling (TLES) examines multiple diffusion trajectories per simulation by accommodating multiple noninteracting copies of a gas molecule that diffuse independently, while the protein and water molecules experience an average interaction from all copies. Furthermore, gas migration within a protein matrix can be accelerated without altering protein dynamics by increasing the effective temperature of the TLES copies. These features of TLES enable rapid simulations of gas diffusion within a protein matrix at significantly reduced (∼98%) computational cost. However, the results of TLES and standard MD simulations have not been systematically compared, which limits the adoption of the TLES approach. We address this drawback here by benchmarking TLES against standard MD in the simulation of O2 diffusion in myoglobin (Mb) as a case study since this model system has been extensively characterized. We find that 2 ns TLES and 108 ns standard simulations map the same network of diffusion tunnels in Mb and uncover the same docking sites, barriers, and escape portals. We further discuss the influence of simulation time as well as the number of independent simulations on the O2 population density within the diffusion tunnels and on the sampling of Mb's conformational space as revealed by principal component analysis. Overall, our comprehensive benchmarking reveals that TLES is an appropriate and robust tool for the rapid mapping of gas diffusion in proteins when the kinetic data provided by standard MD are not required. Furthermore, TLES provides explicit ligand diffusion pathways, unlike most rapid methods.

A quantitative model for oxygen uptake and release in a family of hemeproteins

MOTIVATION

Hemeproteins have many diverse functions that largely depend on the rate at which they uptake or release small ligands, like oxygen. These proteins have been extensively studied using either simulations or experiments, albeit only qualitatively and one or two proteins at a time.

RESULTS

We present a physical-chemical model, which uses data obtained exclusively from computer simulations, to describe the uptake and release of oxygen in a family of hemeproteins, called truncated hemoglobins (trHbs). Through a rigorous statistical analysis we demonstrate that our model successfully recaptures all the reported experimental oxygen association and dissociation kinetic rate constants, thus allowing us to establish the key factors that determine the rates at which these hemeproteins uptake and release oxygen. We found that internal tunnels as well as the distal site water molecules control ligand uptake, whereas oxygen stabilization by distal site residues controls ligand release. Because these rates largely determine the functions of these hemeproteins, these approaches will also be important tools in characterizing the trHbs members with unknown functions.

CONTACT

lboechi@ic.fcen.uba.ar

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Oxygen diffusion pathways in a cofactor-independent dioxygenase

A combination of experimental and computational studies reveals the main O₂ diffusion pathways, providing insight into how cofactor-independent oxygenases control stereospecificity and prevent oxidative inactivation.

Molecular oxygen plays an important role in a wide variety of enzymatic reactions. Through recent research efforts combining computational and experimental methods a new view of O₂ diffusion is emerging, where specific channels guide O₂ to the active site. The focus of this work is DpgC, a cofactor-independent oxygenase. Molecular dynamics simulations, together with mutagenesis experiments and xenon-binding data, reveal that O₂ reaches the active site of this enzyme using three main pathways and four different access points. These pathways connect a series of dynamic hydrophobic pockets, concentrating O₂ at a specific face of the enzyme substrate. Extensive molecular dynamics simulations provide information about which pathways are more frequently used. This data is consistent with the results of kinetic measurements on mutants and is difficult to obtain using computational cavity-location methods. Taken together, our results reveal that although DpgC is rare in its ability of activating O₂ in the absence of cofactors or metals, the way O₂ reaches the active site is similar to that reported for other O₂-using proteins: multiple access channels are available, and the architecture of the pathway network can provide regio- and stereoselectivity. Our results point to the existence of common themes in O₂ access that are conserved among very different types of proteins.

Docking and migration of carbon monoxide in nitrogenase: the case for gated pockets from infrared spectroscopy and molecular dynamics

Evidence of a CO docking site near the FeMo cofactor in nitrogenase has been obtained by Fourier transform infrared spectroscopy-monitored low-temperature photolysis. We investigated the possible migration paths for CO from this docking site using molecular dynamics calculations. The simulations support the notion of a gas channel with multiple internal pockets from the active site to the protein exterior. Travel between pockets is gated by the motion of protein residues. Implications for the mechanism of nitrogenase reactions with CO and N2 are discussed.

Identification of Mutational Hot Spots for Substrate Diffusion: Application to Myoglobin

The pathways by which small molecules (substrates or inhibitors) access active sites are a key aspect of the function of enzymes and other proteins. A key problem in designing or altering such proteins is to identify sites for mutation that will have the desired effect on the substrate transport properties. While specific access channels have been invoked in the past, molecular simulations suggest that multiple routes are possible, complicating the analysis. This complexity, however, can be captured by a Markov State Model (MSM) of the ligand diffusion process. We have developed a sensitivity analysis of the resulting rate matrix, which identifies the locations where mutations should have the largest effect on the diffusive on rate. We apply this method to myoglobin, which is the best characterized example both from experiment and simulation. We validate the approach by translating the sensitivity parameter obtained from this method into the CO binding rates in myoglobin upon mutation, resulting in a semi-quantitative correlation with experiments. The model is further validated against an explicit simulation for one of the experimental mutants.

Effects of mutations on the molecular dynamics of oxygen escape from the dimeric hemoglobin of Scapharca inaequivalvis

Like many hemoglobins, the structure of the dimeric hemoglobin from the clam Scapharca inaequivalvis is a “closed bottle” since there is no direct tunnel from the oxygen binding site on the heme to the solvent.  The proximal histidine faces the dimer interface, which consists of the E and F helicies.  This is significantly different from tetrameric vertebrate hemoglobins and brings the heme groups near the subunit interface. The subunit interface is also characterized by an immobile, hydrogen-bonded network of water molecules.  Although there is data which is consistent with the histidine gate pathway for ligand escape, these aspects of the structure would seem to make that pathway less likely. Locally enhanced sampling molecular dynamics are used here to suggest alternative pathways in the wild-type and six mutant proteins. In most cases the point mutations change the selection of exit routes observed in the simulations. Exit via the histidine gate is rarely seem although oxygen molecules do occasionally cross over the interface from one subunit to the other. The results suggest that changes in flexibility and, in some cases, creation of new cavities can explain the effects of the mutations on ligand exit paths.

Full kinetics of CO entry, internal diffusion, and exit in myoglobin from transition-path theory simulations

We use Markovian milestoning molecular dynamics (MD) simulations on a tessellation of the collective variable space for CO localization in myoglobin to estimate the kinetics of entry, exit, and internal site-hopping. The tessellation is determined by analysis of the free-energy surface in that space using transition-path theory (TPT), which provides criteria for defining optimal milestones, allowing short, independent, cell-constrained MD simulations to provide properly weighted kinetic data. We coarse grain the resulting kinetic model at two levels: first, using crystallographically relevant internal cavities and their predicted interconnections and solvent portals; and second, as a three-state side-path scheme inspired by similar models developed from geminate recombination experiments. We show semiquantitative agreement with experiment on entry and exit rates and in the identification of the so-called "histidine gate" at position 64 through which ≈90% of flux between solvent and the distal pocket passes. We also show with six-dimensional calculations that the minimum free-energy pathway of escape through the histidine gate is a "knock-on" mechanism in which motion of the ligand and the gate are sequential and interdependent. In total, these results suggest that such TPT simulations are indeed a promising approach to overcome the practical time-scale limitations of MD to allow reliable estimation of transition mechanisms and rates among metastable states.

Ligand uptake in Mycobacterium tuberculosis truncated hemoglobins is controlled by both internal tunnels and active site water molecules

Mycobacterium tuberculosis, the causative agent of human tuberculosis, has two proteins belonging to the truncated hemoglobin (trHb) family. Mt-trHbN presents well-defined internal hydrophobic tunnels that allow O 2 and •NO to migrate easily from the solvent to the active site, whereas Mt-trHbO possesses tunnels interrupted by a few bulky residues, particularly a tryptophan at position G8. Differential ligand migration rates allow Mt-trHbN to detoxify •NO, a crucial step for pathogen survival once under attack by the immune system, much more efficiently than Mt-trHbO. In order to investigate the differences between these proteins, we performed experimental kinetic measurements, •NO decomposition, as well as molecular dynamics simulations of the wild type Mt-trHbN and two mutants, VG8F and VG8W. These mutations affect both the tunnels accessibility as well as the affinity of distal site water molecules, thus modifying the ligand access to the iron. We found that a single mutation allows Mt-trHbN to acquire ligand migration rates comparable to those observed for Mt-trHbO, confirming that ligand migration is regulated by the internal tunnel architecture as well as by water molecules stabilized in the active site.

Ligand uptake in Mycobacterium tuberculosis truncated hemoglobins is controlled by both internal tunnels and active site water molecules

Mycobacterium tuberculosis, the causative agent of human tuberculosis, has two proteins belonging to the truncated hemoglobin (trHb) family. Mt-trHbN presents well-defined internal hydrophobic tunnels that allow O 2 and (•)NO to migrate easily from the solvent to the active site, whereas Mt-trHbO possesses tunnels that are partially blocked by a few bulky residues, particularly a tryptophan at position G8. Differential ligand migration rates allow Mt-trHbN to detoxify (•)NO, a crucial step for pathogen survival once under attack by the immune system, much more efficiently than Mt-trHbO. In order to investigate the differences between these proteins, we performed experimental kinetic measurements, (•)NO decomposition, as well as molecular dynamics simulations of the wild type Mt-trHbN and two mutants, VG8F and VG8W. These mutations introduce modifications in both tunnel topologies and affect the incoming ligand capacity to displace retained water molecules at the active site. We found that a single mutation allows Mt-trHbN to acquire ligand migration rates comparable to those observed for Mt-trHbO, confirming that ligand migration is regulated by the internal tunnel architecture as well as by water molecules stabilized in the active site.

Bridging Theory and Experiment to Address Structural Properties of Truncated Haemoglobins: Insights from Thermobifida fusca HbO

In this chapter, we will discuss the paradigmatic case of Thermobifida fusca (Tf-trHb) HbO in its ferrous and ferric states and its behaviour towards a battery of possible ligands. This choice was dictated by the fact that it has been one of the most extensively studied truncated haemoglobins, both in terms of spectroscopic and molecular dynamics studies. Tf-trHb typifies the structural properties of group II trHbs, as the active site is characterized by a highly polar distal environment in which TrpG8, TyrCD1, and TyrB10 provide three potential H-bond donors in the distal cavity capable of stabilizing the incoming ligands. The role of these residues in key topological positions, and their interplay with the iron-bound ligands, has been addressed in studies carried out on the CO, F(-), OH(-), CN(-), and HS(-) adducts formed with the wild-type protein and a combinatorial set of mutants, in which the distal polar residues, TrpG8, TyrCD1, and TyrB10, have been singly, doubly, or triply replaced by a Phe residue. In this context, such a complete analysis provides an excellent benchmark for the investigation of the relationship between protein structure and function, allowing one to translate physicochemical properties of the active site into the observed functional behaviour. Tf-trHb will be compared with other members of the group II trHbs and, more generally, with members of the other trHb subgroups.

Principal Component Analysis reveals correlation of cavities evolution and functional motions in proteins

Protein conformation has been recognized as the key feature determining biological function, as it determines the position of the essential groups specifically interacting with substrates. Hence, the shape of the cavities or grooves at the protein surface appears to drive those functions. However, only a few studies describe the geometrical evolution of protein cavities during molecular dynamics simulations (MD), usually with a crude representation. To unveil the dynamics of cavity geometry evolution, we developed an approach combining cavity detection and Principal Component Analysis (PCA). This approach was applied to four systems subjected to MD (lysozyme, sperm whale myoglobin, Dengue envelope protein and EF-CaM complex). PCA on cavities allows us to perform efficient analysis and classification of the geometry diversity explored by a cavity. Additionally, it reveals correlations between the evolutions of the cavities and structures, and can even suggest how to modify the protein conformation to induce a given cavity geometry. It also helps to perform fast and consensual clustering of conformations according to cavity geometry. Finally, using this approach, we show that both carbon monoxide (CO) location and transfer among the different xenon sites of myoglobin are correlated with few cavity evolution modes of high amplitude. This correlation illustrates the link between ligand diffusion and the dynamic network of internal cavities.

Solvation and cavity occupation in biomolecules

BACKGROUND

Solvation density locations are important for protein dynamics and structure. Knowledge of the preferred hydration sites at biomolecular interfaces and those in the interior of cavities can enhance understanding of structure and function. While advanced X-ray diffraction methods can provide accurate atomic structures for proteins, that technique is challenged when it comes to providing accurate hydration structures, especially for interfacial and cavity bound solvent molecules.

METHODS

Advances in integral equation theories which include more accurate methods for calculating the long-ranged Coulomb interaction contributions to the three-dimensional distribution functions make it possible to calculate angle dependent average solvent structure, accurately, around and inside irregular molecular conformations. The proximal radial distribution method provides another approximate method to determine average solvent structures for biomolecular systems based on a proximal or near neighbor solvent distribution that can be constructed from previously collected solvent distributions. These two approximate methods, along with all-atom molecular dynamics simulations are used to determine the solvent density inside the myoglobin heme cavity.

DISCUSSION AND RESULTS

Myoglobin is a good test system for these methods because the cavities are many and one is large, tens of Å(3), but is shown to have only four hydration sites. These sites are not near neighbors which implies that the large cavity must have more than one way in and out.

CONCLUSIONS

Our results show that main solvation sites are well reproduced by all three methods. The techniques also produce a clearly identifiable solvent pathway into the interior of the protein.

GENERAL SIGNIFICANCE

The agreement between molecular dynamics and less computationally demanding approximate methods is encouraging. This article is part of a Special Issue entitled Recent developments of molecular dynamics.

Internal water and microsecond dynamics in myoglobin

Myoglobin (Mb) binds diatomic ligands, like O₂, CO, and NO, in a cavity that is only transiently accessible. Crystallography and molecular simulations show that the ligands can migrate through an extensive network of transiently connected cavities but disagree on the locations and occupancy of internal hydration sites. Here, we use water ²H and ¹⁷O magnetic relaxation dispersion (MRD) to characterize the internal water molecules in Mb under physiological conditions. We find that equine carbonmonoxy Mb contains 4.5 ± 1.0 ordered internal water molecules with a mean survival time of 5.6 ± 0.5 μs at 25 °C. The likely locations of these water molecules are the four polar hydration sites, including one of the xenon-binding cavities, that are fully occupied in all high-resolution crystal structures of equine Mb. The finding that water escapes from these sites, located 17–31 Å apart in the protein, on the same μs time scale suggests a global exchange mechanism. We propose that this mechanism involves transient penetration of the protein by H-bonded water chains. Such a mechanism could play a functional role by eliminating trapped ligands. In addition, the MRD results indicate that 2 or 3 of the 11 histidine residues of equine Mb undergo intramolecular hydrogen exchange on a μs time scale.

Functional consequences of the open distal pocket of dehaloperoxidase-hemoglobin observed by time-resolved X-ray crystallography

Using time-resolved X-ray crystallography, we contrast a bifunctional dehaloperoxidase-hemoglobin (DHP) with previously studied examples of myoglobin and hemoglobin to understand the functional role of the distal pocket of globins. One key functional difference between DHP and other globins is the requirement that H2O2 enter the distal pocket of oxyferrous DHP to displace O2 from the heme Fe atom and thereby activate the heme for the peroxidase function. The open architecture of DHP permits more than one molecule to simultaneously enter the distal pocket of the protein above the heme to facilitate the unique peroxidase cycle starting from the oxyferrous state. The time-resolved X-ray data show that the distal pocket of DHP lacks a protein valve found in the two other globins that have been studied previously. The photolyzed CO ligand trajectory in DHP does not have a docking site; rather, the CO moves immediately to the Xe-binding site. From there, CO can escape but can also recombine an order of magnitude more rapidly than in other globins. The contrast with DHP dynamics and function more precisely defines the functional role of the multiple conformational states of myoglobin. Taken together with the high reduction potential of DHP, the open distal site helps to explain how a globin can also function as a peroxidase.

Uncovering a dynamically formed substrate access tunnel in carbon monoxide dehydrogenase/acetyl-CoA synthase

The transport of small ligands to active sites of proteins is the basis of vital processes in biology such as enzymatic catalysis and cell signaling, but also of more destructive ones including enzyme inhibition and oxidative damage. Here, we show how a diffusion-reaction model solved by means of molecular dynamics and density functional theory calculations provides novel insight into the transport of small ligands in proteins. In particular, we unravel the existence of an elusive, dynamically formed gas channel, which CO2 takes to diffuse from the solvent to the active site (C-cluster) of the bifunctional multisubunit enzyme complex carbon monoxide dehydrogenase/acetyl-CoA synthase (CODH/ACS). Two cavities forming this channel are temporarily created by protein fluctuations and are not apparent in the X-ray structures. The ligand transport is controlled by two residues at the end of this tunnel, His113 and His116, and occurs on the same time scale on which chemical binding to the active site takes place (0.1-1 ms), resulting in an overall binding rate on the second time scale. We find that upon reduction of CO2 to CO, the newly formed Fe-hydroxy ligand greatly strengthens the hydrogen-bond network, preventing CO from exiting the protein through the same way that CO2 takes to enter the protein. This is the basis for directional transport of CO from the production site (C-cluster of CODH subunit) to the utilization site (A-cluster of ACS subunit). In view of these results, a general picture emerges of how large proteins guide small ligands toward their active sites.

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.

Comparative analysis of inner cavities and ligand migration in non-symbiotic AHb1 and AHb2

This study reports a comparative analysis of the topological properties of inner cavities and the intrinsic dynamics of non-symbiotic hemoglobins AHb1 and AHb2 from Arabidopsis thaliana. The two proteins belong to the 3/3 globin fold and have a sequence identity of about 60%. However, it is widely assumed that they have distinct physiological roles. In order to investigate the structure-function relationships in these proteins, we have examined the bis-histidyl and ligand-bound hexacoordinated states by atomistic simulations using in silico structural models. The results allow us to identify two main pathways to the distal cavity in the bis-histidyl hexacoordinated proteins. Nevertheless, a larger accessibility to small gaseous molecules is found in AHb2. This effect can be attributed to three factors: the mutation Leu35(AHb1)→Phe32(AHb2), the enhanced flexibility of helix B, and the more favorable energetic profile for ligand migration to the distal cavity. The net effect of these factors would be to facilitate the access of ligands, thus compensating the preference for the fully hexacoordination of AHb2, in contrast to the equilibrium between hexa- and pentacoordinated species in AHb1. On the other hand, binding of the exogenous ligand introduces distinct structural changes in the two proteins. A well-defined tunnel is formed in AHb1, which might be relevant to accomplish the proposed NO detoxification reaction. In contrast, no similar tunnel is found in AHb2, which can be ascribed to the reduced flexibility of helix E imposed by the larger number of salt bridges compared to AHb1. This feature would thus support the storage and transport functions proposed for AHb2. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.

Small ligand-globin interactions: reviewing lessons derived from computer simulation

In this work we review the application of classical and quantum-mechanical atomistic computer simulation tools to the investigation of small ligand interaction with globins. In the first part, studies of ligand migration, with its connection to kinetic association rate constants (kon), are presented. In the second part, we review studies for a variety of ligands such as O2, NO, CO, HS(-), F(-), and NO2(-) showing how the heme structure, proximal effects, and the interactions with the distal amino acids can modulate protein ligand binding. The review presents mainly results derived from our previous works on the subject, in the context of other theoretical and experimental studies performed by others. The variety and extent of the presented data yield a clear example of how computer simulation tools have, in the last decade, contributed to our deeper understanding of small ligand interactions with globins. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.

Molecular dynamics study of hell's gate globin I (HGbI) from a methanotrophic extremophile: oxygen migration through a large cavity

Hell's gate globin I (HGbI), a heme-containing protein from the extremophile Methylacidiphilum infernorum, has fast oxygen-binding/slow release characteristics due to its distal residues Gln and Tyr. The combination of Gln/Tyr distal iron coordination, adaptation to extreme environmental conditions, and lack of a D helix suggests that ligand migration in HGbI differs from other previously studied globins. Locally enhanced molecular dynamics trajectories of oxygen migration indicate a large internal cavity. This may increase the tendency of oxygen to exit from portals other than the most direct exit from the space near the heme. Oxygen may reside transiently in shallow surface depressions around the exits. Such surface trapping may enhance both oxygen uptake by increasing contact time between molecules, and decrease release by increasing the probability of oxygen reentry from the vicinity of the portal.