Identification of ligand-binding pathways in truncated hemoglobins using locally enhanced sampling molecular dynamics

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.

Substrate transport pathway inside outward open conformation of EmrD: a molecular dynamics simulation study

The EmrD transporter, which is a classical major facilitator superfamily (MFS) protein, can extrude a range of drug molecules out of E. coil. The drug molecules transport through the channel of MFS in an outward open state, an important issue in research about bacterial drug resistance, which however, is still unknown. In this paper, we construct a starting outward-open model of the EmrD transporter using a state transition method. The starting model is refined by a conventional molecular dynamics simulation. Locally enhanced sampling simulation (LES) is used to validate the outward-open model of EmrD. In the locally enhanced sampling simulation, ten substrates are placed along the channel of the outward-open EmrD, and these substrates are sampled in the outward-open center cavity. It is found that the translocation pathway of these substrates from the inside to the outside of the cell through the EmrD transporter is composed of two sub-pathways, one sub-pathway, including H2, H4, and H5, and another sub-pathway, including H8, H10, and H11. The results give us have a further insight to the ways of substrate translocation of an MFS protein. The model method is based on common features of an MFS protein, so this modeling method can be used to construct various MFS protein models which have a desired state with other conformations not known in the alternating-access mechanism.

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.

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.

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.

Volumetric properties underlying ligand binding in a monomeric hemoglobin: a high-pressure NMR study

The 2/2 hemoglobin of the cyanobacterium Synechococcus sp. PCC 7002, GlbN, coordinates the heme iron with two histidines and exists either with a b heme or with a covalently attached heme. The binding of exogenous ligands displaces the distal histidine and induces a conformational rearrangement involving the reorganization of internal void volumes. The formation of passageways within the resulting conformation is thought to facilitate ligand exchange and play a functional role. Here we monitored the perturbation induced by pressure on the ferric bis-histidine and cyanide-bound states of GlbN using (1)H-(15)N HSQC NMR spectroscopy. We inspected the outcome with a statistical analysis of 170 homologous 2/2 hemoglobin sequences. We found that the compression landscape of GlbN, as represented by the variation of an average chemical shift parameter, was highly sensitive to ligand swapping and heme covalent attachment. Stabilization of rare conformers was observed at high pressures and consistent with cavity redistribution upon ligand binding. In all states, the EF loop was found to be exceptionally labile to pressure, suggesting a functional role as a semi-flexible hinge between the adjacent helices. Finally, coevolved clusters presented a common pattern of compensating pressure responses. The high-pressure dissection combined with protein sequence analysis established locations with volumetric signatures relevant to residual communication of 2/2 hemoglobins. 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.

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.

Structure and dynamics of Mycobacterium tuberculosis truncated hemoglobin N: insights from NMR spectroscopy and molecular dynamics simulations

The potent nitric oxide dioxygenase (NOD) activity (trHbN-Fe²⁺-O₂ + (•)NO → trHbN-Fe³⁺-OH₂ + NO₃⁻) of Mycobacterium tuberculosis truncated hemoglobin N (trHbN) protects aerobic respiration from inhibition by (•)NO. The high activity of trHbN has been attributed in part to the presence of numerous short-lived hydrophobic cavities that allow partition and diffusion of the gaseous substrates (•)NO and O₂ to the active site. We investigated the relation between these cavities and the dynamics of the protein using solution NMR spectroscopy and molecular dynamics (MD). Results from both approaches indicate that the protein is mainly rigid with very limited motions of the backbone N-H bond vectors on the picoseconds-nanoseconds time scale, indicating that substrate diffusion and partition within trHbN may be controlled by side-chains movements. Model-free analysis also revealed the presence of slow motions (microseconds-milliseconds), not observed in MD simulations, for many residues located in helices B and G including the distal heme pocket Tyr33(B10). All currently known crystal structures and molecular dynamics data of truncated hemoglobins with the so-called pre-A N-terminal extension suggest a stable α-helical conformation that extends in solution. Moreover, a recent study attributed a crucial role to the pre-A helix for NOD activity. However, solution NMR data clearly show that in near-physiological conditions these residues do not adopt an α-helical conformation and are significantly disordered and that the helical conformation seen in crystal structures is likely induced by crystal contacts. Although this lack of order for the pre-A does not disagree with an important functional role for these residues, our data show that one should not assume an helical conformation for these residues in any functional interpretation. Moreover, future molecular dynamics simulations should not use an initial α-helical conformation for these residues in order to avoid a bias based on an erroneous initial structure for the N-termini residues. This work constitutes the first study of a truncated hemoglobin dynamics performed by solution heteronuclear relaxation NMR spectroscopy.

Role of the distal hydrogen-bonding network in regulating oxygen affinity in the truncated hemoglobin III from Campylobacter jejuni

Oxygen affinity in heme-containing proteins is determined by a number of factors, such as the nature and conformation of the distal residues that stabilize the heme bound-oxygen via hydrogen-bonding interactions. The truncated hemoglobin III from Campylobacter jejuni (Ctb) contains three potential hydrogen-bond donors in the distal site: TyrB10, TrpG8, and HisE7. Previous studies suggested that Ctb exhibits an extremely slow oxygen dissociation rate due to an interlaced hydrogen-bonding network involving the three distal residues. Here we have studied the structural and kinetic properties of the G8(WF) mutant of Ctb and employed state-of-the-art computer simulation methods to investigate the properties of the O(2) adduct of the G8(WF) mutant, with respect to those of the wild-type protein and the previously studied E7(HL) and/or B10(YF) mutants. Our data indicate that the unique oxygen binding properties of Ctb are determined by the interplay of hydrogen-bonding interactions between the heme-bound ligand and the surrounding TyrB10, TrpG8, and HisE7 residues.

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.

Structural characterization of the tunnels of Mycobacterium tuberculosis truncated hemoglobin N from molecular dynamics simulations

The structure of oxygenated trHbN from Mycobacterium tuberculosis shows an extended heme distal hydrogen-bond network that includes Tyr33(B10), Gln58(E11), and the bound O(2). In addition, trHbN structure shows a network of hydrophobic cavities organized in two orthogonal branches. In the present work, the structure and the dynamics of oxygenated and deoxygenated trHbN in explicit water was investigated from 100 ns molecular dynamics (MD) simulations. Results show that, depending on the presence or the absence of a coordinated O(2), the Tyr33(B10) and Gln58(E11) side chains adopt two different configurations in concert with hydrogen bond network rearrangement. In addition, our data indicate that Tyr33(B10) and Gln58(E11) control the dynamics of Phe62(E15). In deoxy-trHbN, Phe62(E15) is restricted to one conformation. Upon O(2) binding, the conformation of Gln58(E11) changes and residue Phe62(E15) fluctuates between two conformations. We also conducted a systematic study of trHbN tunnels by analyzing thousands of MD snapshots with CAVER. The results show that tunnel formation is the result of the dynamic reshaping of short-lived hydrophobic cavities. The analyses indicate that the presence of these cavities is likely linked to the rigid structure of trHbN and also reveal two tunnels, EH and BE, that link the protein surface to the buried distal heme pocket and not present in the crystallographic structure. The cavities are sufficiently large to accomodate and store ligands. Tunnel dynamics in trHbN was found to be controlled by the side-chain conformation of the Tyr33(B10), Gln58(E11), and Phe62(E15) residues. Importantly, in contrast to recently published works, our extensive systematic studies show that the presence or absence of a coordinated dioxygen does not control the opening of the long tunnel but rather the opening of the EH tunnel. In addition, the data lead to new and distinctly different conclusion on the impact of the Phe62(E15) residue on trHbN tunnels. We propose that the EH and the long tunnels are used for apolar ligands storage. The trajectories bring important new structural insights related to trHbN function and to ligand diffusion in proteins.