Structural determinants of ligand migration in Mycobacterium tuberculosis truncated hemoglobin O

Nitric oxide binding to ferrous nitrobindins: A computer simulation investigation

Nitrobindins (Nbs) represent an evolutionary conserved all-β-barrel heme-proteins displaying a highly solvent-exposed heme-Fe(III) atom, coordinated by a proximal His residue. Interestingly, even if the distal side is exposed to the solvent, the value of the second order rate constants for ligand binding to the ferrous derivative is almost one order of magnitude lower than those reported for myoglobins (Mbs). Noteworthy, nitric oxide binding to the sixth coordination position of the heme-Fe(II)-atom causes the cleavage or the severe weakening of the proximal His-Fe(II) bond. Here, we provide a computer simulation investigation to shed light on the molecular basis of ligand binding kinetics, by dissecting the ligand binding process into the ligand migration and the bond formation steps. Classical molecular dynamics simulations were performed employing a steered molecular dynamics approach and the Jarzinski equality to obtain ligand migration free energy profiles. The formation of the heme-Fe(II)-NO bond took into consideration the iron atom displacement from the heme plane. The ligand migration is almost unhindered, and the low rate constant for NO binding is due to the large displacement of the Fe(II) atom with respect to the heme plane responsible for the barrier for the Fe(II)-NO bond formation. In addition, we investigated the weakening and breaking of the proximal His-Fe(II) bond, observed experimentally upon NO binding, by means of a combination of classical molecular dynamics simulations and quantum-classical (QM-MM) optimizations. In both human and M. tuberculosis Nbs, a stable alternative conformation of the proximal His residue interacting with a network of water molecules was observed.

Truncated (2/2) hemoglobin: Unconventional structures and functional roles in vivo and in human pathogenesis

Truncated hemoglobins (trHbs) build a sub-class of the globin family, found in eubacteria, cyanobacteria, unicellular eukaryotes, and in higher plants; among these, selected human pathogens are found. The trHb fold is based on a 2/2 α-helical sandwich, consisting of a simplified and reduced-size version of the classical 3/3 α-helical sandwich of vertebrate and invertebrate globins. Phylogenetic analysis indicates that trHbs further branch into three groups: group I (or trHbN), group II (or trHbO), and group III (or trHbP), each group being characterized by specific structural features. Among these, a protein matrix tunnel, or a cavity system implicated in diatomic ligand diffusion through the protein matrix, is typical of group I and group II, respectively. In general, a highly intertwined network of hydrogen bonds stabilizes the heme bound ligand, despite variability of the heme distal residues in the different trHb groups. Notably, some organisms display genes from more than one trHb group, suggesting that trHbN, trHbO, and trHbP may support different functions in vivo, such as detoxification of reactive nitrogen and oxygen species, respiration, oxygen storage/sensoring, thus aiding survival of an invading microorganism. Here, structural features and proposed functions of trHbs from human pathogens are reviewed.

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.

Lessons learned about steered molecular dynamics simulations and free energy calculations

The calculation of free energy profiles is central in understanding differential enzymatic activity, for instance, involving chemical reactions that require QM-MM tools, ligand migration, and conformational rearrangements that can be modeled using classical potentials. The use of steered molecular dynamics (sMD) together with the Jarzynski equality is a popular approach in calculating free energy profiles. Here, we first briefly review the application of the Jarzynski equality to sMD simulations, then revisit the so-called stiff-spring approximation and the consequent expectation of Gaussian work distributions and, finally, reiterate the practical utility of the second-order cumulant expansion, as it coincides with the parametric maximum-likelihood estimator in this scenario. We illustrate this procedure using simulations of CO, both in aqueous solution and in a carbon nanotube as a model system for biologically relevant nanoheterogeneous environments. We conclude the use of the second-order cumulant expansion permits the use of faster pulling velocities in sMD simulations, without introducing bias due to large dispersion in the non-equilibrium work distribution.

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.

Tertiary and quaternary structural basis of oxygen affinity in human hemoglobin as revealed by multiscale simulations

Human hemoglobin (Hb) is a benchmark protein of structural biology that shaped our view of allosterism over 60 years ago, with the introduction of the MWC model based on Perutz structures of the oxy(R) and deoxy(T) states and the more recent Tertiary Two-State model that proposed the existence of individual subunit states -“r” and “t”-, whose structure is yet unknown. Cooperative oxygen binding is essential for Hb function, and despite decades of research there are still open questions related to how tertiary and quaternary changes regulate oxygen affinity. In the present work, we have determined the free energy profiles of oxygen migration and for HisE7 gate opening, with QM/MM calculations of the oxygen binding energy in order to address the influence of tertiary differences in the control of oxygen affinity. Our results show that in the α subunit the low to high affinity transition is achieved by a proximal effect that mostly affects oxygen dissociation and is the driving force of the allosteric transition, while in the β subunit the affinity change results from a complex interplay of proximal and distal effects, including an increase in the HE7 gate opening, that as shown by free energy profiles promotes oxygen uptake.

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.

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.

Structural and Functional Significance of the N- and C-Terminal Appendages in Arabidopsis Truncated Hemoglobin

Plant hemoglobins constitute three distinct groups: symbiotic, nonsymbiotic, and truncated hemoglobins. Structural investigation of symbiotic and nonsymbiotic (class I) hemoglobins revealed the presence of a vertebrate-like 3/3 globin fold in these proteins. In contrast, plant truncated hemoglobins are similar to bacterial truncated hemoglobins with a putative 2/2 α-helical globin fold. While multiple structures have been reported for plant hemoglobins of the first two categories, for plant truncated globins only one structure has been reported of late. Here, we report yet another crystal structure of the truncated hemoglobin from Arabidopsis thaliana (AHb3) with two water molecules in the heme pocket, of which one is distinctly coordinated to the heme iron, unlike the only available crystal structure of AHb3 with a hydroxyl ligand. AHb3 was monomeric in its crystallographic asymmetric unit; however, dimer was evident in the crystallographic symmetry, and the globin indeed existed as a stable dimer in solution. The tertiary structure of the protein exhibited a bacterial-like 2/2 α-helical globin fold with an additional N-terminal α-helical extension and disordered C-termini. To address the role of these extended termini in AHb3, which is yet unknown, N- and C-terminal deletion mutants were created and characterized and molecular dynamics simulations performed. The C-terminal deletion had an insignificant effect on most properties but perturbed the dimeric equilibrium of AHb3 and significantly influenced azide binding kinetics in the ferric state. These results along with the disordered nature of the C-terminus indicated its putative role in intramolecular or intermolecular interactions probably regulating protein-ligand and protein-protein interactions. While the N-terminal deletion did not change the overall globin fold, stability, or ligand binding kinetics, it seemed to have influenced coordination at the heme iron, the hydration status of the active site, and the quaternary structure of AHb3. Evidence indicated that the N-terminus is the predominant factor regulating the quaternary interaction appropriate to physiological requirements, dynamics of the side chains in the heme pocket, and tunnel organization in the protein matrix.

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.

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.

Structural flexibility of the heme cavity in the cold-adapted truncated hemoglobin from the Antarctic marine bacterium Pseudoalteromonas haloplanktis TAC125

Truncated hemoglobins build one of the three branches of the globin protein superfamily. They display a characteristic two-on-two α-helical sandwich fold and are clustered into three groups (I, II and III) based on distinct structural features. Truncated hemoglobins are present in eubacteria, cyanobacteria, protozoa and plants. Here we present a structural, spectroscopic and molecular dynamics characterization of a group-II truncated hemoglobin, encoded by the PSHAa0030 gene from Pseudoalteromonas haloplanktis TAC125 (Ph-2/2HbO), a cold-adapted Antarctic marine bacterium hosting one flavohemoglobin and three distinct truncated hemoglobins. The Ph-2/2HbO aquo-met crystal structure (at 2.21 Å resolution) shows typical features of group-II truncated hemoglobins, namely the two-on-two α-helical sandwich fold, a helix Φ preceding the proximal helix F, and a heme distal-site hydrogen-bonded network that includes water molecules and several distal-site residues, including His(58)CD1. Analysis of Ph-2/2HbO by electron paramagnetic resonance, resonance Raman and electronic absorption spectra, under varied solution conditions, shows that Ph-2/2HbO can access diverse heme ligation states. Among these, detection of a low-spin heme hexa-coordinated species suggests that residue Tyr(42)B10 can undergo large conformational changes in order to act as the sixth heme-Fe ligand. Altogether, the results show that Ph-2/2HbO maintains the general structural features of group-II truncated hemoglobins but displays enhanced conformational flexibility in the proximity of the heme cavity, a property probably related to the functional challenges, such as low temperature, high O2 concentration and low kinetic energy of molecules, experienced by organisms living in the Antarctic environment.

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.

The diversity of 2/2 (truncated) globins

Small size globins that have been defined as 'truncated haemoglobins' or as '2/2 haemoglobins' have increasingly been discovered in microorganisms since the early 1990s. Analysis of amino acid sequences allowed to distinguish three groups that collect proteins with specific and common structural properties. All three groups display 3D structures that are based on four main α-helices, which are a subset of the conventional eight-helices globin fold. Specific features, such as the presence of protein matrix tunnels that are held to promote diffusion of functional ligands to/from the haem, distinguish members of the three groups. Haem distal sites vary for their accessibility, local structures, polarity, and ligand stabilization mechanisms, suggesting functional roles that are related to O2/NO chemistry. In a few cases, such activities have been proven in vitro and in vivo through deletion mutants. The issue of 2/2 haemoglobin varied biological functions throughout the three groups remains however fully open.

Haemoglobins of Mycobacteria: structural features and biological functions

The genus Mycobacterium is comprised of Gram-positive bacteria occupying a wide range of natural habitats and includes species that range from severe intracellular pathogens to economically useful and harmless microbes. The recent upsurge in the availability of microbial genome data has shown that genes encoding haemoglobin-like proteins are ubiquitous among Mycobacteria and that multiple haemoglobins (Hbs) of different classes may be present in pathogenic and non-pathogenic species. The occurrence of truncated haemoglobins (trHbs) and flavohaemoglobins (flavoHbs) showing distinct haem active site structures and ligand-binding properties suggests that these Hbs may be playing diverse functions in the cellular metabolism of Mycobacteria. TrHbs and flavoHbs from some of the severe human pathogens such as Mycobacterium tuberculosis and Mycobacterium leprae have been studied recently and their roles in effective detoxification of reactive nitrogen and oxygen species, electron cycling, modulation of redox state of the cell and facilitation of aerobic respiration have been proposed. This multiplicity in the function of Hbs may aid these pathogens to cope with various environmental stresses and survive during their intracellular regime. This chapter provides recent updates on genomic, structural and functional aspects of Mycobacterial Hbs to address their role in Mycobacteria.

Ligand uptake modulation by internal water molecules and hydrophobic cavities in hemoglobins

Internal water molecules play an active role in ligand uptake regulation, since displacement of retained water molecules from protein surfaces or cavities by incoming ligands can promote favorable or disfavorable effects over the global binding process. Detection of these water molecules by X-ray crystallography is difficult given their positional disorder and low occupancy. In this work, we employ a combination of molecular dynamics simulations and ligand rebinding over a broad time range to shed light into the role of water molecules in ligand migration and binding. Computational studies on the unliganded structure of the thermostable truncated hemoglobin from Thermobifida fusca (Tf-trHbO) show that a water molecule is in the vicinity of the iron heme, stabilized by WG8 with the assistance of YCD1, exerting a steric hindrance for binding of an exogenous ligand. Mutation of WG8 to F results in a significantly lower stabilization of this water molecule and in subtle dynamical structural changes that favor ligand binding, as observed experimentally. Water is absent from the fully hydrophobic distal cavity of the triple mutant YB10F-YCD1F-WG8F (3F), due to the lack of residues capable of stabilizing it nearby the heme. In agreement with these effects on the barriers for ligand rebinding, over 97% of the photodissociated ligands are rebound within a few nanoseconds in the 3F mutant case. Our results demonstrate the specific involvement of water molecules in shaping the energetic barriers for ligand migration and binding.

QM/MM study of the C-C coupling reaction mechanism of CYP121, an essential cytochrome p450 of Mycobacterium tuberculosis

Among 20 p450s of Mycobacterium tuberculosis (Mt), CYP121 has received an outstanding interest, not only due to its essentiality for bacterial viability but also because it catalyzes an unusual carbon-carbon coupling reaction. Based on the structure of the substrate bound enzyme, several reaction mechanisms were proposed involving first Tyr radical formation, second Tyr radical formation, and C-C coupling. Key and unknown features, being the nature of the species that generate the first and second radicals, and the role played by the protein scaffold each step. In the present work we have used classical and quantum based computer simulation methods to study in detail its reaction mechanism. Our results show that substrate binding promotes formation of the initial oxy complex, Compound I is the responsible for first Tyr radical formation, and that the second Tyr radical is formed subsequently, through a PCET reaction, promoted by the presence of key residue Arg386. The final C-C coupling reaction possibly occurs in bulk solution, thus yielding the product in one oxygen reduction cycle. Our results thus contribute to a better comprehension of MtCYP121 reaction mechanism, with direct implications for inhibitor design, and also contribute to our general understanding of these type of enzymes.

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.

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.

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.

The key role of water in the dioxygenase function of Escherichia coli flavohemoglobin

Flavohemoglobins (FHbs) are members of the globin superfamily, widely distributed among prokaryotes and eukaryotes that have been shown to carry out nitric oxide dioxygenase (NOD) activity. In prokaryotes, such as Escherichia coli, NOD activity is a defence mechanism against the NO release by the macrophages of the hosts' immune system during infection. Because of that, FHbs have been studied thoroughly and several drugs have been developed in an effort to fight infectious processes. Nevertheless, the protein's structural determinants involved in the NOD activity are still poorly understood. In this context, the aim of the present work is to unravel the molecular basis of FHbs structural dynamics-to-function relationship using state of the art computer simulation tools. In an effort to fulfill this goal, we studied three key processes that determine NOD activity, namely i) ligand migration into the active site ii) stabilization of the coordinated oxygen and iii) intra-protein electron transfer (ET). Our results allowed us to determine key factors related to all three processes like the presence of a long hydrophobic tunnel for ligand migration, the presence of a water mediated hydrogen bond to stabilize the coordinated oxygen and therefore achieve a high affinity, and the best possible ET paths between the FAD and the heme, where water molecules play an important role. Taken together the presented results close an important gap in our understanding of the wide and diverse globin structural-functional relationships.

Following ligand migration pathways from picoseconds to milliseconds in type II truncated hemoglobin from Thermobifida fusca

CO recombination kinetics has been investigated in the type II truncated hemoglobin from Thermobifida fusca (Tf-trHb) over more than 10 time decades (from 1 ps to ∼100 ms) by combining femtosecond transient absorption, nanosecond laser flash photolysis and optoacoustic spectroscopy. Photolysis is followed by a rapid geminate recombination with a time constant of ∼2 ns representing almost 60% of the overall reaction. An additional, small amplitude geminate recombination was identified at ∼100 ns. Finally, CO pressure dependent measurements brought out the presence of two transient species in the second order rebinding phase, with time constants ranging from ∼3 to ∼100 ms. The available experimental evidence suggests that the two transients are due to the presence of two conformations which do not interconvert within the time frame of the experiment. Computational studies revealed that the plasticity of protein structure is able to define a branched pathway connecting the ligand binding site and the solvent. This allowed to build a kinetic model capable of describing the complete time course of the CO rebinding kinetics to Tf-trHb.

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.

Blocking the gate to ligand entry in human hemoglobin

His(E7) to Trp replacements in HbA lead to markedly biphasic bimolecular CO rebinding after laser photolysis. For isolated mutant subunits, the fraction of fast phase increases with increasing [CO], suggesting a competition between binding to an open conformation with an empty E7 channel and relaxation to blocked or closed, slowly reacting states. The rate of conformational relaxation of the open state is ∼18,000 s(-1) in α subunits and ∼10-fold faster in β subunits, ∼175,000 s(-1). Crystal structures were determined for tetrameric α(WT)β(Trp-63) HbCO, α(Trp-58)β(WT) deoxyHb, and Trp-64 deoxy- and CO-Mb as controls. In Trp-63(E7) βCO, the indole side chain is located in the solvent interface, blocking entry into the E7 channel. Similar blocked Trp-64(E7) conformations are observed in the mutant Mb crystal structures. In Trp-58(E7) deoxy-α subunits, the indole side chain fills both the channel and the distal pocket, forming a completely closed state. The bimolecular rate constant for CO binding, k'(CO), to the open conformations of both mutant Hb subunits is ∼80-90 μm(-1) s(-1), whereas k'(CO) for the completely closed states is 1000-fold slower, ∼0.08 μm(-1) s(-1). A transient intermediate with k'(CO) ≈ 0.7 μm(-1) s(-1) is observed after photolysis of Trp-63(E7) βCO subunits and indicates that the indole ring blocks the entrance to the E7 channel, as observed in the crystal structures of Trp(E7) deoxyMb and βCO subunits. Thus, either blocking or completely filling the E7 channel dramatically slows bimolecular binding, providing strong evidence that the E7 channel is the major pathway (≥90%) for ligand entry in human hemoglobin.

Exploring the molecular basis of human manganese superoxide dismutase inactivation mediated by tyrosine 34 nitration

Manganese Superoxide Dismutase (MnSOD) is an essential mitochondrial antioxidant enzyme that protects organisms against oxidative damage, dismutating superoxide radical (O₂(.)⁻) into H₂O₂ and O₂. The active site of the protein presents a Mn ion in a distorted trigonal-bipyramidal environment, coordinated by H26, H74, H163, D159 and one ⁻OH ion or H₂O molecule. The catalytic cycle of the enzyme is a "ping-pong" mechanism involving Mn³+/Mn²+. It is known that nitration of Y34 is responsible for enzyme inactivation, and that this protein oxidative modification is found in tissues undergoing inflammatory and degenerative processes. However, the molecular basis about MnSOD tyrosine nitration affects the protein catalytic function is mostly unknown. In this work we strongly suggest, using computer simulation tools, that Y34 nitration affects protein function by restricting ligand access to the active site. In particular, deprotonation of 3-nitrotyrosine increases drastically the energetic barrier for ligand entry due to the absence of the proton. Our results for the WT and selected mutant proteins confirm that the phenolic moiety of Y34 plays a key role in assisting superoxide migration.

Heme pocket structural properties of a bacterial truncated hemoglobin from Thermobifida fusca

An acidic surface variant (ASV) of the "truncated" hemoglobin from Thermobifida fusca was designed with the aim of creating a versatile globin scaffold endowed with thermostability and a high level of recombinant expression in its soluble form while keeping the active site unmodified. This engineered protein was obtained by mutating the surface-exposed residues Phe107 and Arg91 to Glu. Molecular dynamics simulations showed that the mutated residues remain solvent-exposed, not affecting the overall protein structure. Thus, the ASV was used in a combinatorial mutagenesis of the distal heme pocket residues in which one, two, or three of the conserved polar residues [TyrB10(54), TyrCD1(67), and TrpG8(119)] were substituted with Phe. Mutants were characterized by infrared and resonance Raman spectroscopy and compared with the wild-type protein. Similar Fe-proximal His stretching frequencies suggest that none of the mutations alters the proximal side of the heme cavity. Two conformers were observed in the spectra of the CO complexes of both wild-type and ASV protein: form 1 with ν(FeC) and ν(CO) at 509 and 1938 cm(-1) and form 2 with ν(FeC) and ν(CO) at 518 and 1920 cm(-1), respectively. Molecular dynamics simulations were performed for the wild-type and ASV forms, as well as for the TyrB10 mutant. The spectroscopic and computational results demonstrate that CO interacts with TrpG8 in form 1 and interacts with both TrpG8 and TyrCD1 in form 2. TyrB10 does not directly interact with the bound CO.

The peculiar heme pocket of the 2/2 hemoglobin of cold-adapted Pseudoalteromonas haloplanktis TAC125

The genome of the cold-adapted bacterium Pseudoalteromonas haloplanktis TAC125 contains multiple genes encoding three distinct monomeric hemoglobins exhibiting a 2/2 α-helical fold. In the present work, one of these hemoglobins is studied by resonance Raman, electronic absorption and electronic paramagnetic resonance spectroscopies, kinetic measurements, and different bioinformatic approaches. It is the first cold-adapted bacterial hemoglobin to be characterized. The results indicate that this protein belongs to the 2/2 hemoglobin family, Group II, characterized by the presence of a tryptophanyl residue on the bottom of the heme distal pocket in position G8 and two tyrosyl residues (TyrCD1 and TyrB10). However, unlike other bacterial hemoglobins, the ferric state, in addition to the aquo hexacoordinated high-spin form, shows multiple hexacoordinated low-spin forms, where either TyrCD1 or TyrB10 can likely coordinate the iron. This is the first example in which both TyrCD1 and TyrB10 are proposed to be the residues that are alternatively involved in heme hexacoordination by endogenous ligands.

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.

The first step of the dioxygenation reaction carried out by tryptophan dioxygenase and indoleamine 2,3-dioxygenase as revealed by quantum mechanical/molecular mechanical studies

Tryptophan dioxygenase (TDO) and indoleamine 2,3-dioxygenase (IDO) are two heme-containing enzymes which catalyze the conversion of L: -tryptophan to N-formylkynurenine (NFK). In mammals, TDO is mostly expressed in liver and is involved in controlling homeostatic serum tryptophan concentrations, whereas IDO is ubiquitous and is involved in modulating immune responses. Previous studies suggested that the first step of the dioxygenase reaction involves the deprotonation of the indoleamine group of the substrate by an evolutionarily conserved distal histidine residue in TDO and the heme-bound dioxygen in IDO. Here, we used classical molecular dynamics and hybrid quantum mechanical/molecular mechanical methods to evaluate the base-catalyzed mechanism. Our data suggest that the deprotonation of the indoleamine group of the substrate by either histidine in TDO or heme-bound dioxygen in IDO is not energetically favorable. Instead, the dioxygenase reaction can be initiated by a direct attack of heme-bound dioxygen on the C(2)=C(3) bond of the indole ring, leading to a protein-stabilized 2,3-alkylperoxide transition state and a ferryl epoxide intermediate, which subsequently recombine to generate NFK. The novel sequential two-step oxygen addition mechanism is fully supported by our recent resonance Raman data that allowed identification of the ferryl intermediate (Lewis-Ballester et al. in Proc Natl Acad Sci USA 106:17371-17376, 2009). The results reveal the subtle differences between the TDO and IDO reactions and highlight the importance of protein matrix in modulating stereoelectronic factors for oxygen activation and the stabilization of both transition and intermediate states.

Ligand migration in the truncated hemoglobin-II from Mycobacterium tuberculosis: the role of G8 tryptophan

Resonance Raman studies show that the heme-bound CO in trHbO, a truncated-II hemoglobin from Mycobacterium tuberculosis, is exposed to an environment with a positive electrostatic potential. The mutation of Trp(G8), an absolutely conserved residue in group II and III truncated hemoglobins, to Phe introduces two new Fe-CO conformers, both of which exhibit reduced electrostatic potentials. Computer simulations reveal that the structural perturbation is a result of the increased flexibility of the Tyr(CD1) and Leu(E11) side chains due to the reduction of the size of the G8 residue. Laser flash photolysis studies show that the G8 mutation induces 1) the presence of two new geminate recombination phases, one with a rate faster than the time resolution of our instrument and the other with a rate 13-fold slower than that of the wild type protein, and 2) the reduction of the total geminate recombination yield from 86 to 62% and the increase in the bimolecular recombination rate by a factor of 530. Computer simulations uncover that the photodissociated ligand migrates between three distal temporary docking sites before it subsequently rebinds to the heme iron or ultimately escapes into the solvent via a hydrophobic tunnel. The calculated energy profiles associated with the ligand migration processes are in good agreement with the experimental observations. The results highlight the importance of the Trp(G8) in regulating ligand migration in trHbO, underscoring its pivotal role in the structural and functional properties of the group II and III truncated hemoglobins.