Oxygen migration pathways in NO-bound truncated hemoglobin

Multipolar Force Fields for Amide-I Spectroscopy from Conformational Dynamics of the Alanine Trimer

The dynamics and spectroscopy of N-methyl-acetamide (NMA) and trialanine in solution are characterized from molecular dynamics simulations using different energy functions, including a conventional point charge (PC)-based force field, one based on a multipolar (MTP) representation of the electrostatics, and a semiempirical DFT method. For the 1D infrared spectra, the frequency splitting between the two amide-I groups is 10 cm^-1 from the PC, 13 cm^-1 from the MTP, and 47 cm^-1 from self-consistent charge density functional tight-binding (SCC-DFTB) simulations, compared with 25 cm^-1 from experiment. The frequency trajectory required for the frequency fluctuation correlation function (FFCF) is determined from individual normal mode (INM) and full normal mode (FNM) analyses of the amide-I vibrations. The spectroscopy, time-zero magnitude of the FFCF C(t = 0), and the static component Δ₀² from simulations using MTP and analysis based on FNM are all consistent with experiments for (Ala)₃. Contrary to this, for the analysis excluding mode-mode coupling (INM), the FFCF decays to zero too rapidly and for simulations with a PC-based force field, the Δ₀² is too small by a factor of two compared with experiments. Simulations with SCC-DFTB agree better with experiment for these observables than those from PC-based simulations. The conformational ensemble sampled from simulations using PCs is consistent with the literature (including P_II, β, α_R, and α_L), whereas that covered by the MTP-based simulations is dominated by P_II with some contributions from β and α_R. This agrees with and confirms recently reported Bayesian-refined populations based on 1D infrared experiments. FNM analysis together with a MTP representation provides a meaningful model to correctly describe the dynamics of hydrated trialanine.

Allostery in the nitric oxide dioxygenase mechanism of flavohemoglobin

The substrates O₂ and NO cooperatively activate the NO dioxygenase function of Escherichia coli flavohemoglobin. Steady-state and transient kinetic measurements support a structure-based mechanistic model in which O₂ and NO movements and conserved amino acids at the E11, G8, E2, E7, B10, and F7 positions within the globin domain control activation. In the cooperative and allosteric mechanism, O₂ migrates to the catalytic heme site via a long hydrophobic tunnel and displaces LeuE11 away from the ferric iron, which forces open a short tunnel to the catalytic site gated by the ValG8/IleE15 pair and LeuE11. NO permeates this tunnel and leverages upon the gating side chains triggering the CD loop to furl, which moves the E and F-helices and switches an electron transfer gate formed by LysF7, GlnE7, and water. This allows FADH₂ to reduce the ferric iron, which forms the stable ferric–superoxide–TyrB10/GlnE7 complex. This complex reacts with internalized NO with a bimolecular rate constant of 10¹⁰ M⁻¹ s⁻¹ forming nitrate, which migrates to the CD loop and unfurls the spring-like structure. To restart the cycle, LeuE11 toggles back to the ferric iron. Actuating electron transfer with O₂ and NO movements averts irreversible NO poisoning and reductive inactivation of the enzyme. Together, structure snapshots and kinetic constants provide glimpses of intermediate conformational states, time scales for motion, and associated energies.

Structural and functional properties of Antarctic fish cytoglobins-1: Cold-reactivity in multi-ligand reactions

While the functions of the recently discovered cytoglobin, ubiquitously expressed in vertebrate tissues, remain uncertain, Antarctic fish provide unparalleled models to study novel protein traits that may arise from cold adaptation. We report here the spectral, ligand-binding and enzymatic properties (peroxynitrite isomerization, nitrite-reductase activity) of cytoglobin-1 from two Antarctic fish, Chaenocephalus aceratus and Dissostichus mawsoni, and present the crystal structure of D. mawsoni cytoglobin-1. The Antarctic cytoglobins-1 display high O₂ affinity, scarcely compatible with an O₂-supply role, a slow rate constant for nitrite-reductase activity, and do not catalyze peroxynitrite isomerization. Compared with mesophilic orthologues, the cold-adapted cytoglobins favor binding of exogenous ligands to the hexa-coordinated bis-histidyl species, a trait related to their higher rate constant for distal-His/heme-Fe dissociation relative to human cytoglobin. At the light of a remarkable 3D-structure conservation, the observed differences in ligand-binding kinetics may reflect Antarctic fish cytoglobin-1 specific features in the dynamics of the heme distal region and of protein matrix cavities, suggesting adaptation to functional requirements posed by the cold environment. Taken together, the biochemical and biophysical data presented suggest that in Antarctic fish, as in humans, cytoglobin-1 unlikely plays a role in O₂ transport, rather it may be involved in processes such as NO detoxification.

Ligand Binding Rate Constants in Heme Proteins Using Markov State Models and Molecular Dynamics Simulations

Computer simulation studies of the molecular basis for ligand migration in proteins allow the description of key events such as the transition between docking sites, displacement of existing ligands and solvent molecules, and open/closure of specific "gates", among others. In heme proteins, ligand migration from the solvent to the active site preludes the binding to the heme iron and triggers different functions. In this work, molecular dynamics simulations, a Markov State Model of migration and empirical kinetic equations are combined to study the migration of O₂ and NO in two truncated hemoglobins of Mycobacterium tuberculosis (Mt-TrHbN and Mt-TrHbO). For Mt-TrHbN, we show that the difference in the association constant in the oxy and deoxy states relies mainly in the displacement of water molecules anchored in the distal cavity in the deoxy form. The results here provide a valuable approach to study ligand migration in globins.

Kinetic Analysis and Structural Interpretation of Competitive Ligand Binding for NO Dioxygenation in Truncated Hemoglobin N

The conversion of nitric oxide (NO) into nitrate (NO₃^- ) by dioxygenation protects cells from lethal NO. Starting from NO-bound heme, the first step in converting NO into benign NO₃^- is the ligand exchange reaction FeNO+O₂ →FeO₂ +NO, which is still poorly understood at a molecular level. For wild-type (WT) truncated hemoglobin N (trHbN) and its Y33A mutant, the calculated barriers for the exchange reaction differ by 1.5 kcal mol^-1 , compared with 1.7 kcal mol^-1 from experiment. It is directly confirmed that the ligand exchange reaction is rate-limiting in trHbN and that entropic contributions account for 75 % of the difference between the WT and the mutant. Residues Tyr 33, Phe 46, Val 80, His 81, and Gln 82 surrounding the active site are expected to control the reaction path. By comparison with electronic structure calculations, the transition state separating the two ligand-bound states was assigned to a ² A state.

Implications of short time scale dynamics on long time processes

This review provides a comprehensive overview of the structural dynamics in topical gas- and condensed-phase systems on multiple length and time scales. Starting from vibrationally induced dissociation of small molecules in the gas phase, the question of vibrational and internal energy redistribution through conformational dynamics is further developed by considering coupled electron/proton transfer in a model peptide over many orders of magnitude. The influence of the surrounding solvent is probed for electron transfer to the solvent in hydrated I-. Next, the dynamics of a modified PDZ domain over many time scales is analyzed following activation of a photoswitch. The hydration dynamics around halogenated amino acid side chains and their structural dynamics in proteins are relevant for iodinated TyrB26 insulin. Binding of nitric oxide to myoglobin is a process for which experimental and computational analyses have converged to a common view which connects rebinding time scales and the underlying dynamics. Finally, rhodopsin is a paradigmatic system for multiple length- and time-scale processes for which experimental and computational methods provide valuable insights into the functional dynamics. The systems discussed here highlight that for a comprehensive understanding of how structure, flexibility, energetics, and dynamics contribute to functional dynamics, experimental studies in multiple wavelength regions and computational studies including quantum, classical, and more coarse grained levels are required.

Migration of small ligands in globins: Xe diffusion in truncated hemoglobin N

In heme proteins, the efficient transport of ligands such as NO or O₂ to the binding site is achieved via ligand migration networks. A quantitative assessment of ligand diffusion in these networks is thus essential for a better understanding of the function of these proteins. For this, Xe migration in truncated hemoglobin N (trHbN) of Mycobacterium Tuberculosis was studied using molecular dynamics simulations. Transitions between pockets of the migration network and intra-pocket relaxation occur on similar time scales (10 ps and 20 ps), consistent with low free energy barriers (1-2 kcal/mol). Depending on the pocket from where Xe enters a particular transition, the conformation of the side chains lining the transition region differs which highlights the coupling between ligand and protein degrees of freedom. Furthermore, comparison of transition probabilities shows that Xe migration in trHbN is a non-Markovian process. Memory effects arise due to protein rearrangements and coupled dynamics as Xe moves through it.

Binding and transport of ligands in proteins is essential, in particular in globular proteins which often exhibit internal cavities. In truncated Hemoglobin N (trHbN) these cavities are arranged as a network with particular connectivities. The present work supports the notion that ligand diffusion in trHbN is an active process and coupled to the protein dynamics. Furthermore, transition probabilities between neighboring pockets depend on the location from where the ligand entered the transition, which is typical for non-Markovian processes. Hence, ligand migration in trHbN exhibits memory effects due to dynamical coupling between the protein and ligand motion.

Drug targeting of heme proteins in Mycobacterium tuberculosis

TB, caused by the human pathogen Mycobacterium tuberculosis (Mtb), causes more deaths than any other infectious disease. Iron is crucial for Mtb to infect the host and to sustain infection, with Mtb encoding large numbers of iron-binding proteins. Many of these are hemoproteins with key roles, including defense against oxidative stress, cellular signaling and regulation, host cholesterol metabolism, and respiratory processes. Various heme enzymes in Mtb are validated drug targets and/or products of genes essential for bacterial viability or survival in the host. Here, we review the structure, function, and druggability of key Mtb heme enzymes and strategies used for their inhibition.

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.

A comparative analysis of clustering algorithms: O2 migration in truncated hemoglobin I from transition networks

The ligand migration network for O2-diffusion in truncated Hemoglobin N is analyzed based on three different clustering schemes. For coordinate-based clustering, the conventional k-means and the kinetics-based Markov Clustering (MCL) methods are employed, whereas the locally scaled diffusion map (LSDMap) method is a collective-variable-based approach. It is found that all three methods agree well in their geometrical definition of the most important docking site, and all experimentally known docking sites are recovered by all three methods. Also, for most of the states, their population coincides quite favourably, whereas the kinetics of and between the states differs. One of the major differences between k-means and MCL clustering on the one hand and LSDMap on the other is that the latter finds one large primary cluster containing the Xe1a, IS1, and ENT states. This is related to the fact that the motion within the state occurs on similar time scales, whereas structurally the state is found to be quite diverse. In agreement with previous explicit atomistic simulations, the Xe3 pocket is found to be a highly dynamical site which points to its potential role as a hub in the network. This is also highlighted in the fact that LSDMap cannot identify this state. First passage time distributions from MCL clusterings using a one- (ligand-position) and two-dimensional (ligand-position and protein-structure) descriptor suggest that ligand- and protein-motions are coupled. The benefits and drawbacks of the three methods are discussed in a comparative fashion and highlight that depending on the questions at hand the best-performing method for a particular data set may differ.

Vibrational relaxation and energy migration of N-methylacetamide in water: the role of nonbonded interactions

Nonequilibrium molecular dynamics (MD) simulations together with physics-based force fields are used to follow energy flow between vibrationally excited N-methylacetamide (NMA) and water. The simulations are carried out with a previously validated force field for NMA, based on a multipolar representation of the electrostatics, and with a new fluctuating point charge model. For the water solvent, a flexible and a rigid model was employed to distinguish between the role of inter- and intramolecular degrees of freedom. On a 10 ps time scale about 90% of the available energy goes into the solvent. The remaining energy resides within internal NMA-degrees of freedom from where energy flow takes place on longer time scales. The total amount of energy transferred to the solvent on the 10 ps time scale does not depend on whether the water molecules are rigid or flexible during the simulations. Vibrational energy relaxation time scales include two regimes: one on the several 100 fs time scale and a longer one, ranging from 6 to 10 ps. This longer time scale agrees with previous simulations but overestimates the experimentally determined relaxation time by a factor of 2, which can be explained by the classical treatment of the vibrations. Including a previously determined quantum correction factor brings the long time scale into quite favorable agreement with experiment. Coupling to the bending vibration of the water molecules in H-bonding contact with the excited C═O chromophore is substantial. The equilibrium and nonequilibrium distribution of the bending angles of the water molecules in contact with the local oscillator are non-Gaussian, and one approaches the other on the subpicosecond time scale. Analysis of the water velocity distribution suggests that the C═O vibrational energy relaxes into the solvent water shells in an impulsive fashion on a picosecond time scale.

Coupled protein-ligand dynamics in truncated hemoglobin N from atomistic simulations and transition networks

BACKGROUND

The nature of ligand motion in proteins is difficult to characterize directly using experiment. Specifically, it is unclear to what degree these motions are coupled.

METHODS

All-atom simulations are used to sample ligand motion in truncated Hemoglobin N. A transition network analysis including ligand- and protein-degrees of freedom is used to analyze the microscopic dynamics.

RESULTS

Clustering of two different subsets of MD trajectories highlights the importance of a diverse and exhaustive description to define the macrostates for a ligand-migration network. Monte Carlo simulations on the transition matrices from one particular clustering are able to faithfully capture the atomistic simulations. Contrary to clustering by ligand positions only, including a protein degree of freedom yields considerably improved coarse grained dynamics. Analysis with and without imposing detailed balance agree closely which suggests that the underlying atomistic simulations are converged with respect to sampling transitions between neighboring sites.

CONCLUSIONS

Protein and ligand dynamics are not independent from each other and ligand migration through globular proteins is not passive diffusion.

GENERAL SIGNIFICANCE

Transition network analysis is a powerful tool to analyze and characterize the microscopic dynamics in complex systems. This article is part of a Special Issue entitled Recent developments of molecular dynamics.

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.

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.