Ligand migration in the apolar tunnel of Cerebratulus lacteus mini-hemoglobin

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.

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.

Structure and function of haemoglobins

Haemoglobin (Hb) is widely known as the iron-containing protein in blood that is essential for O₂ transport in mammals. Less widely recognised is that erythrocyte Hb belongs to a large family of Hb proteins with members distributed across all three domains of life-bacteria, archaea and eukaryotes. This review, aimed chiefly at researchers new to the field, attempts a broad overview of the diversity, and common features, in Hb structure and function. Topics include structural and functional classification of Hbs; principles of O₂ binding affinity and selectivity between O₂/NO/CO and other small ligands; hexacoordinate (containing bis-imidazole coordinated haem) Hbs; bacterial truncated Hbs; flavohaemoglobins; enzymatic reactions of Hbs with bioactive gases, particularly NO, and protection from nitrosative stress; and, sensor Hbs. A final section sketches the evolution of work on the structural basis for allosteric O₂ binding by mammalian RBC Hb, including the development of newer kinetic models. Where possible, reference to historical works is included, in order to provide context for current advances in Hb research.

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.

Structural Plasticity in Globins: Role of Protein Dynamics in Defining Ligand Migration Pathways

Globins are a family of proteins characterized by the presence of the heme prosthetic group and involved in variety of biological functions in the cell. Due to their biological relevance and widespread distribution in all kingdoms of life, intense research efforts have been devoted to disclosing the relationships between structural features, protein dynamics, and function. Particular attention has been paid to the impact of differences in amino acid sequence on the topological features of docking sites and cavities and to the influence of conformational flexibility in facilitating the migration of small ligands through these cavities. Often, tunnels are carved in the interior of globins, and ligand exchange is regulated by gating residues. Understanding the subtle intricacies that relate the differences in sequence with the structural and dynamical features of globins with the ultimate aim of rationalizing the thermodynamics and kinetics of ligand binding continues to be a major challenge in the field. Due to the evolution of computational techniques, significant advances into our understanding of these questions have been made. In this review we focus our attention on the analysis of the ligand migration pathways as well as the function of the structural cavities and tunnels in a series of representative globins, emphasizing the synergy between experimental and theoretical approaches to gain a comprehensive knowledge into the molecular mechanisms of this diverse family of proteins.

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.

Understanding the kinetics of ligand binding to globins with molecular dynamics simulations: the necessity of multiple state models

Residue conformational changes and internal cavity migration processes play a key role in regulating the kinetics of ligand migration and binding events in globins. Molecular dynamics simulations have demonstrated their value in the study of these processes in different haemoglobins, but derivation of kinetic data demands the use of more complex techniques like enhanced sampling molecular dynamics methods. This review discusses the different methodologies that are currently applied to study the ligand migration process in globins and highlight those specially developed to derive kinetic data.

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.

Kinetic properties and heme pocket structure of two domains of the polymeric hemoglobin of Artemia in comparison with the native molecule

In this project, we studied some physicochemical properties of two different globin domains of the polymeric hemoglobin of the brine shrimp Artemia salina and compared them with those of the native molecule. Two domains (AsHbC1D1 and AsHbC1D5) were cloned and expressed in BL21(DE3)pLysS strain of Escherichia coli. The recombinant proteins as well as the native hemoglobin (AfHb) were purified from bacteria and frozen Artemia, respectively by standard chromatographic methods and assessed by SDS-PAGE. The heme environment of these proteins was studied by optical spectroscopy and ligand-binding kinetics (e.g. CO association and O2 binding affinity) were measured for the two recombinant proteins and the native hemoglobin. This indicates that the CO association rate for AsHbC1D1 is higher than that of AsHbC1D5 and AfHb, while the calculated P50 value for AsHbC1D1 is lower than that of AsHbC1D5 and AfHb. The geminate and bimolecular rebinding parameters indicate a significant difference between both domains. Moreover, EPR results showed that the heme pocket in AfHb is in a more closed conformation than the heme pocket in myoglobin. Finally, the reduction potential of -0.13V versus the standard hydrogen electrode was determined for AfHb by direct electrochemical measurements. It is about 0.06V higher than the potential of the single domain AsHbC1D5. This work shows that each domain in the hemoglobin of Artemia has different characteristics of ligand binding.

Heme-based globin-coupled oxygen sensors: linking oxygen binding to functional regulation of diguanylate cyclase, histidine kinase, and methyl-accepting chemotaxis

An emerging class of novel heme-based oxygen sensors containing a globin fold binds and senses environmental O2 via a heme iron complex. Structure-function relationships of oxygen sensors containing a heme-bound globin fold are different from those containing heme-bound PAS and GAF folds. It is thus worth reconsidering from an evolutionary perspective how heme-bound proteins with a globin fold similar to that of hemoglobin and myoglobin could act as O2 sensors. Here, we summarize the molecular mechanisms of heme-based oxygen sensors containing a globin fold in an effort to shed light on the O2-sensing properties and O2-stimulated catalytic enhancement observed for these proteins.

Single vs. multiple ligand pathways in globins: a computational view

Diatomic ligand migration in globins has been the subject of numerous studies. Still, a consensus picture for the ligand entrance is not clear, with a growing concern among experimental researchers that computational simulations always show multiple pathways for any globin. Modeling non-biased ligand entrance from conventional molecular dynamics techniques, however, has shown to be difficult (and expensive). Here we use our Monte Carlo methodology, capable of freely mapping ligand diffusion and the description of rare events, to two well-studied systems: myoglobin and the mini-hemoglobin from the sea worm Cerebratulus lacteus. Our results clearly show that the simulations are specific to the system providing a different trend in the entrance pathway, as expected from experiments. While Mb presents multiple entrance pathways, populating the well-known xenon cavities, in CerHb the ligand enters the protein only by one apolar channel. Most of the trajectories (64%) visiting myoglobin's active site though, are gated by the distal histidine. Such detailed information, accessible through the state of the art algorithms in PELE, is computationally inexpensive and available to all non-profit researchers. 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.

Ligand binding to heme proteins: a comparison of cytochrome c variants with globins

We have studied the binding of carbon monoxide (CO) in mutants of Cyt c having its methionine at position 80 replaced by alanine, aspartate, and arginine, so that the sixth coordination is available for ligand binding. We have employed Fourier transform infrared (FTIR) photolysis difference spectroscopy to examine interactions of the heme-bound and photolyzed CO (and also nitric oxide, NO) in the small heme pocket created by the mutations. By using FTIR temperature derivative spectroscopy (TDS) and nanosecond flash photolysis, the enthalpy barrier distributions for CO rebinding were determined. In flash photolysis experiments, the majority of ligands rebind to the heme iron on picosecond time scales so that only the high-barrier tail of the distributions is visible on the nanosecond scale. By continuous wave excitation prior to TDS characterization of the barriers, however, each Cyt c molecule is photoexcited multiple times and complete photodissociation can be achieved, which likely arises from a rotation of the CO within the heme pocket so that the oxygen faces the heme iron. Apparently, reorientation prior to rebinding constitutes an additional and significant contribution to the rebinding barrier. Our experiments reveal that the compact, rigid structure of Cyt c offers no alternative binding sites for photodissociated ligands in the protein matrix. A comparison of ligand binding in these Cyt c mutants and hemoglobins underscores the importance of internal ligand docking sites and ligand migration routes for conveying a ligand binding function to heme proteins.

Determination of ligand pathways in globins: apolar tunnels versus polar gates

Although molecular dynamics simulations suggest multiple interior pathways for O(2) entry into and exit from globins, most experiments indicate well defined single pathways. In 2001, we highlighted the effects of large-to-small amino acid replacements on rates for ligand entry and exit onto the three-dimensional structure of sperm whale myoglobin. The resultant map argued strongly for ligand movement through a short channel from the heme iron to solvent that is gated by the distal histidine (His-64(E7)) near the solvent edge of the porphyrin ring. In this work, we have applied the same mutagenesis mapping strategy to the neuronal mini-hemoglobin from Cerebratulus lacteus (CerHb), which has a large internal tunnel from the heme iron to the C-terminal ends of the E and H helices, a direction that is 180° opposite to the E7 channel. Detailed comparisons of the new CerHb map with expanded results for Mb show unambiguously that the dominant (>90%) ligand pathway in CerHb is through the internal tunnel, and the major (>75%) ligand pathway in Mb is through the E7 gate. These results demonstrate that: 1) mutagenesis mapping can identify internal pathways when they exist; 2) molecular dynamics simulations need to be refined to address discrepancies with experimental observations; and 3) alternative pathways have evolved in globins to meet specific physiological demands.

High resolution crystal structures of the Cerebratulus lacteus mini-Hb in the unligated and carbomonoxy states

The nerve tissue mini-hemoglobin from Cerebratulus lacteus (CerHb) displays an essential globin fold hosting a protein matrix tunnel held to allow traffic of small ligands to and from the heme. CerHb heme pocket hosts the distal TyrB10/GlnE7 pair, normally linked to low rates of O₂ dissociation and ultra-high O₂ affinity. However, CerHb affinity for O₂ is similar to that of mammalian myoglobins, due to a dynamic equilibrium between high and low affinity states driven by the ability of ThrE11 to orient the TyrB10 OH group relative to the heme ligand. We present here the high resolution crystal structures of CerHb in the unligated and carbomonoxy states. Although CO binds to the heme with an orientation different from the O₂ ligand, the overall binding schemes for CO and O₂ are essentially the same, both ligands being stabilized through a network of hydrogen bonds based on TyrB10, GlnE7, and ThrE11. No dramatic protein structural changes are needed to support binding of the ligands, which can freely reach the heme distal site through the apolar tunnel. A lack of main conformational changes between the heme-unligated and -ligated states grants stability to the folded mini-Hb and is a prerequisite for fast ligand diffusion to/from the heme.

Modulating distal cavities in the α and β subunits of human HbA reveals the primary ligand migration pathway

The free volume in the active site of human HbA plays a crucial role in governing the bimolecular rates of O(2), CO, and NO binding, the fraction of geminate ligand recombination, and the rate of NO dioxygenation by the oxygenated complex. We have decreased the size of the distal pocket by mutating Leu(B10), Val(E11), and Leu(G8) to Phe and Trp and that of other more internal cavities by filling them with Xe at high gas pressures. Increasing the size of the B10 side chain reduces bimolecular rates of ligand binding nearly 5000-fold and inhibits CO geminate recombination due to both reduction of the capture volume in the distal pocket and direct steric hindrance of Fe-ligand bond formation. Phe and Trp(E11) mutations also cause a decrease in distal pocket volume but, at the same time, increase access to the Fe atom because of the loss of the γ2 CH(3) group of the native Val(E11) side chain. The net result of these E11 substitutions is a dramatic increase in the rate of geminate recombination because dissociated CO is sequestered close to the Fe atom and can rapidly rebind without steric resistance. However, the bimolecular rate constants for binding of ligand to the Phe and Trp(E11) mutants are decreased 5-30-fold, because of a smaller capture volume. Geminate and bimolecular kinetic parameters for Phe and Trp(G8) mutants are similar to those for the native HbA subunits because the aromatic rings at this position cause little change in distal pocket volume and because ligands do not move past this position into the globin interior of wild-type HbA subunits. The latter conclusion is verified by the observation that Xe binding to the α and β Hb subunits has little effect on either geminate or bimolecular ligand rebinding. All of these experimental results argue strongly against alternative ligand migration pathways that involve movements through the protein interior in HbA. Instead, ligands appear to enter through the His(E7) gate and are captured directly in the distal cavity.