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

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.

Binding mechanism of disulfide species to ferric hemeproteins: The case of metmyoglobin

The interactions of the heme iron of hemeproteins with sulfide and disulfide compounds are of potential interest as physiological signaling processes. While the interaction with hydrogen sulfide has been described computationally and experimentally, the reaction with disulfide, and specifically the molecular mechanism for ligand binding has not been studied in detail. In this work, we study the association process for disulfane and its conjugate base disulfanide at different pH conditions. Additionally, by means of advanced sampling techniques based on multiple steered molecular dynamics, we provide free energy profiles for ligand migration for both acid/base species, showing a similar behavior to the previously reported for the related H2S/HS¯ pair. Finally, we studied the ligand interchange reaction (H2O/H2S, HS¯ and H2O/HSSH, HSS¯) by means of hybrid quantum mechanics-molecular mechanics calculations. We show that the anionic species are able to displace more efficiently the H2O bound to the iron, and that the H-bond network in the distal cavity can help the neutral species to perform the reaction. Altogether, we provide a molecular explanation for the experimental information and show that the global association process depends on a fine balance between the migration towards the active site and the ligand interchange reaction.

The role of the distal cavity in carbon monoxide stabilization in the coproheme decarboxylase enzyme from C. diphtheriae

This work focuses on the carbon monoxide adducts of the wild-type and selected variants of the coproheme decarboxylase from actinobacterial Corynebacterium diphtheriae complexed with coproheme, monovinyl monopropionyl deuteroheme (MMD), and heme b. The UV - vis and resonance Raman spectroscopies together with the molecular dynamics simulations clearly show that the wild-type coproheme-CO adduct is characterized by two CO conformers, one hydrogen-bonded to the distal H118 residue and the other showing a weak polar interaction with the distal cavity. Instead, upon conversion to heme b, i.e. after decarboxylation of propionates 2 and 4 and rotation by 90o of the porphyrin ring inside the cavity, CO probes a less polar environment. In the absence of the H118 residue, both coproheme and heme b complexes form only the non-H-bonded CO species. The unrotated MMD-CO adduct as observed in the H118F variant, confirms that decarboxylation of propionate 2 only, does not affect the heme cavity. The rupture of both the H-bonds involving propionates 2 and 4 destabilizes the porphyrin inside the cavity with the subsequent formation of a CO adduct in an open conformation. In addition, in this work we present data on CO binding to reversed heme b, obtained by hemin reconstitution of the H118A variant, and to heme d, obtained by addition of an excess of hydrogen peroxide. The results will be discussed and compared with those reported for the representatives of the firmicute clade.

Best Practices on QM/MM Simulations of Biological Systems

During the second half of the 20th century, following structural biology hallmark works on DNA and proteins, biochemists shifted their questions from "what does this molecule look like?" to "how does this process work?". Prompted by the theoretical and practical developments in computational chemistry, this led to the emergence of biomolecular simulations and, along with the 2013 Nobel Prize in Chemistry, to the development of hybrid QM/MM methods. QM/MM methods are necessary whenever the problem we want to address involves chemical reactivity and/or a change in the system's electronic structure, with archetypal examples being the studies of an enzyme's reaction mechanism and a metalloprotein's active site. In the last decades QM/MM methods have seen an increasing adoption driven by their incorporation in widely used biomolecular simulation software. However, properly setting up a QM/MM simulation is not an easy task, and several issues need to be properly addressed to obtain meaningful results. In the present work, we describe both the theoretical concepts and practical issues that need to be considered when performing QM/MM simulations. We start with a brief historical perspective on the development of these methods and describe when and why QM/MM methods are mandatory. Then we show how to properly select and analyze the performance of the QM level of theory, the QM system size, and the position and type of the boundaries. We show the relevance of performing prior QM model system (or QM cluster) calculations in a vacuum and how to use the corresponding results to adequately calibrate those derived from QM/MM. We also discuss how to prepare the starting structure and how to select an adequate simulation strategy, including those based on geometry optimizations as well as free energy methods. In particular, we focus on the determination of free energy profiles using multiple steered molecular dynamics (MSMD) combined with Jarzynski's equation. Finally, we describe the results for two illustrative and complementary examples: the reaction performed by chorismate mutase and the study of ligand binding to hemoglobins. Overall, we provide many practical recommendations (or shortcuts) together with important conceptualizations that we hope will encourage more and more researchers to incorporate QM/MM studies into their research projects.

Mycobacterium tuberculosis DosS binds H2S through its Fe3+ heme iron to regulate the DosR dormancy regulon

Mycobacterium tuberculosis (Mtb) senses and responds to host-derived gasotransmitters NO and CO via heme-containing sensor kinases DosS and DosT and the response regulator DosR. Hydrogen sulfide (H2S) is an important signaling molecule in mammals, but its role in Mtb physiology is unclear. We have previously shown that exogenous H2S can modulate expression of genes in the Dos dormancy regulon via an unknown mechanism(s). Here, we test the hypothesis that Mtb senses and responds to H2S via the DosS/T/R system. Using UV-Vis and EPR spectroscopy, we show that H2S binds directly to the ferric (Fe3+) heme of DosS (KDapp = 5.30 μM) but not the ferrous (Fe2+) form. No interaction with DosT(Fe2+-O2) was detected. We found that the binding of sulfide can slowly reduce the DosS heme iron to the ferrous form. Steered Molecular Dynamics simulations show that H2S, and not the charged HS- species, can enter the DosS heme pocket. We also show that H2S increases DosS autokinase activity and subsequent phosphorylation of DosR, and H2S-mediated increases in Dos regulon gene expression is lost in Mtb lacking DosS. Finally, we demonstrate that physiological levels of H2S in macrophages can induce DosR regulon genes via DosS. Overall, these data reveal a novel mechanism whereby Mtb senses and responds to a third host gasotransmitter, H2S, via DosS(Fe3+). These findings highlight the remarkable plasticity of DosS and establish a new paradigm for how bacteria can sense multiple gasotransmitters through a single heme sensor kinase.

Data for molecular dynamics simulations of Escherichia coli cytochrome bd oxidase with the Amber force field

Cytochrome bd-type quinol oxidase is an important metalloenzyme that allows many bacteria to survive in low oxygen conditions. Since bd oxidase is found in many prokaryotes but not in eukaryotes, it has emerged as a promising bacterial drug target. Examples of organisms containing bd oxidases include the Mycobacterium tuberculosis (Mtb) bacterium that causes tuberculosis (TB) in humans, the Vibrio cholerae bacterium that causes cholera, the Pseudomonas aeruginosa bacterium that contributes to antibiotic resistance and sepsis, and the Campylobacter jejuni bacterium that causes food poisoning. Escherichia coli (E. coli) is another organism exhibiting the cytochrome bd oxidase. Since it has the highest sequence identity to Mtb (36%) and we are ultimately interested in finding drug targets for TB, we have built parameters for the E. coli bd oxidase (Protein Data Bank ID number: 6RKO) that are compatible with the all-atom Amber ff14SB force field for molecular dynamics (MD) simulations. Specifically, we built parameters for the three heme cofactors present in all species of bacterial cytochrome bd-type oxidases (heme b558 , heme b595 , and heme d) along with their axial ligands. This data report includes the parameter and library files that can be used with Amber's LEaP program to generate input files for MD simulations using the Amber software package. We also provide the PDB data files of the initial model both by itself and solvated with TIP3P water molecules and counterions.

Heme-Thiolate Perturbation in Cystathionine β-Synthase by Mercury Compounds

Cystathionine β-synthase (CBS) is an enzyme involved in sulfur metabolism that catalyzes the pyridoxal phosphate-dependent condensation of homocysteine with serine or cysteine to form cystathionine and water or hydrogen sulfide (H₂S), respectively. CBS possesses a b-type heme coordinated by histidine and cysteine. Fe(III)-CBS is inert toward exogenous ligands, while Fe(II)-CBS is reactive. Both Fe(III)- and Fe(II)-CBS are sensitive to mercury compounds. In this study, we describe the kinetics of the reactions with mercuric chloride (HgCl₂) and p-chloromercuribenzoic acid. These reactions were multiphasic and resulted in five-coordinate CBS lacking thiolate ligation, with six-coordinate species as intermediates. Computational QM/MM studies supported the feasibility of formation of species in which the thiolate is proximal to both the iron ion and the mercury compound. The reactions of Fe(II)-CBS were faster than those of Fe(III)-CBS. The observed rate constants of the first phase increased hyperbolically with concentration of the mercury compounds, with limiting values of 0.3-0.4 s^-1 for Fe(III)-CBS and 40 ± 4 s^-1 for Fe(II)-CBS. The data were interpreted in terms of alternative models of conformational selection or induced fit. Exposure of Fe(III)-CBS to HgCl₂ led to heme release and activity loss. Our study reveals the complexity of the interactions between mercury compounds and CBS.

The Role of Host-Generated H2S in Microbial Pathogenesis: New Perspectives on Tuberculosis

For centuries, hydrogen sulfide (H₂S) was considered primarily as a poisonous gas and environmental hazard. However, with the discovery of prokaryotic and eukaryotic enzymes for H₂S production, breakdown, and utilization, H₂S has emerged as an important signaling molecule in a wide range of physiological and pathological processes. Hence, H₂S is considered a gasotransmitter along with nitric oxide (•NO) and carbon monoxide (CO). Surprisingly, despite having overlapping functions with •NO and CO, the role of host H₂S in microbial pathogenesis is understudied and represents a gap in our knowledge. Given the numerous reports that followed the discovery of •NO and CO and their respective roles in microbial pathogenesis, we anticipate a rapid increase in studies that further define the importance of H₂S in microbial pathogenesis, which may lead to new virulence paradigms. Therefore, this review provides an overview of sulfide chemistry, enzymatic production of H₂S, and the importance of H₂S in metabolism and immunity in response to microbial pathogens. We then describe our current understanding of the role of host-derived H₂S in tuberculosis (TB) disease, including its influences on host immunity and bioenergetics, and on Mycobacterium tuberculosis (Mtb) growth and survival. Finally, this review discusses the utility of H₂S-donor compounds, inhibitors of H₂S-producing enzymes, and their potential clinical significance.

Conformational Flexibility Drives Cold Adaptation in Pseudoalteromonas haloplanktis TAC125 Globins

Significance: Temperature is one of the most important drivers in shaping protein adaptations. Many biochemical and physiological processes are influenced by temperature. Proteins and enzymes from organisms living at low temperature are less stable in comparison to high-temperature adapted proteins. The lower stability is generally due to greater conformational flexibility. Recent Advances: Adaptive changes in the structure of cold-adapted proteins may occur at subunit interfaces, distant from the active site, thus producing energy changes associated with conformational transitions transmitted to the active site by allosteric modulation, valid also for monomeric proteins in which tertiary structural changes may play an essential role. Critical Issues: Despite efforts, the current experimental and computational methods still fail to produce general principles on protein evolution, since many changes are protein and species dependent. Environmental constraints or other biological cellular signals may override the ancestral information included in the structure of the protein, thus introducing inaccuracy in estimates and predictions on the evolutionary adaptations of proteins in response to cold adaptation. Future Directions: In this review, we describe the studies and approaches used to investigate stability and flexibility in the cold-adapted globins of the Antarctic marine bacterium Pseudoalteromonas haloplanktis TAC125. In fact, future research directions will be prescient on more detailed investigation of cold-adapted proteins and the role of fluctuations between different conformational states.

Distinctive structural properties of THB11, a pentacoordinate Chlamydomonas reinhardtii truncated hemoglobin with N- and C-terminal extensions

Hemoglobins (Hbs) utilize heme b as a cofactor and are found in all kingdoms of life. The current knowledge reveals an enormous variability of Hb primary sequences, resulting in topological, biochemical and physiological individuality. As Hbs appear to modulate their reactivities through specific combinations of structural features, predicting the characteristics of a given Hb is still hardly possible. The unicellular green alga Chlamydomonas reinhardtii contains 12 genes encoding diverse Hbs of the truncated lineage, several of which possess extended N- or C-termini of unknown function. Studies on some of the Chlamydomonas Hbs revealed yet unpredictable structural and biochemical variations, which, along with a different expression of their genes, suggest diverse physiological roles. Chlamydomonas thus represents a promising system to analyze the diversification of Hb structure, biochemistry and physiology. Here, we report the crystal structure, resolved to 1.75 Å, of the heme-binding domain of cyanomet THB11 (Cre16.g662750), one of the pentacoordinate algal Hbs, which offer a free Fe-coordination site in the reduced state. The overall fold of THB11 is conserved, but individual features such as a kink in helix E, a tilted heme plane and a clustering of methionine residues at a putative tunnel exit appear to be unique. Both N- and C-termini promote the formation of oligomer mixtures, and the absence of the C terminus results in reduced nitrite reduction rates. This work widens the structural and biochemical knowledge on the 2/2Hb family and suggests that the N- and C-terminal extensions of the Chlamydomonas 2/2Hbs modulate their reactivity by intermolecular interactions.

Hemeproteins as Targets for Sulfide Species

Significance: Sulfides are endogenous and ubiquitous signaling species that share the hemeproteins as biochemical targets with O₂, nitric oxide, and carbon monoxide. The description of the binding mechanisms is mandatory to anticipate the biochemical relevance of the interaction. Recent Advances: The binding of sulfide to ferric hemeproteins has been described in more than 40 systems, including native proteins, mutants, and model systems. Mechanisms of sulfide binding to ferric hemeproteins have been examined by a combination of kinetic and computational experiments. The distal control of the association process, dissected into the migration of the ligand to the active site and the binding event, reveals that neutral hydrogen sulfide (H₂S) reaches the active site and is the predominant binding ligand, while the HS^- is excluded by the protein matrix. Experiments with model compounds, devoid of a protein scaffold, reveal that both H₂S and HS^- can bind the ferric heme if accessing the site. A critical role of the proximal ligand in the prevention of the metal-centered reduction has been experimentally assessed. For metmyoglobin and methemoglobin, the coordination of sulfide leads to noncanonical functions: sulfide storage and its oxidative detoxification have been evidenced under physiological and excess sulfide concentrations, respectively. Critical Issues: The bound species is suggested to predominate in the monoprotonated form, although spectroscopic evidence is pending. Future Directions: A description of the role of hemeproteins as biochemical targets for inorganic sulfide requires understanding the reactivity of bound sulfide, for example: the metal-centered reduction, the reaction with excess sulfide, oxidants, or other gasotransmitters, among other biomolecules.

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.

Cold-Adaptation Signatures in the Ligand Rebinding Kinetics to the Truncated Hemoglobin of the Antarctic Bacterium Pseudoalteromonas haloplanktis TAC125

Cold-adapted organisms have evolved proteins endowed with higher flexibility and lower stability in comparison to their thermophilic homologues, resulting in enhanced reaction rates at low temperatures. In this context, protein-bound water molecules were suggested to play a major role, and their weaker interactions at protein active sites have been associated with cold adaptation. In this work, we tested this hypothesis on truncated hemoglobins (a family of microbial heme-proteins of yet-unclear function) applying molecular dynamics simulations and ligand-rebinding kinetics on a protein from the Antarctic bacterium Pseudoalteromonas haloplanktis TAC125 in comparison with its thermophilic Thermobifida fusca homologue. The CO rebinding kinetics of the former highlight several geminate phases, with an unusually long-lived geminate intermediate. An articulated tunnel with at least two distinct docking sites was identified by analysis of molecular dynamics simulations and was suggested to be at the origin of the unusual geminate rebinding phase. Water molecules are present in the distal pocket, but their stabilization by TrpG8, TyrB10, and HisCD1 is much weaker than in thermophilic Thermobifida fusca truncated hemoglobin, resulting in a faster geminate rebinding. Our results support the hypothesis that weaker water-molecule interactions at the reaction site are associated with cold adaptation.

Computational approaches for deciphering the equilibrium and kinetic properties of iron transport proteins

With the advances in three-dimensional structure determination techniques, high quality structures of the iron transport proteins transferrin and the bacterial ferric binding protein (FbpA) have been deposited in the past decade. These are proteins of relatively large size, and developments in hardware and software have only recently made it possible to study their dynamics using standard computational resources. We review computational techniques towards understanding the equilibrium and kinetic properties of iron transport proteins under different environmental conditions. At the level of detail that requires quantum chemical treatments, the octahedral geometry around iron has been scrutinized and it has been established that the iron coordinating tyrosines are in an unusual deprotonated state. At the atomistic level, both the N-lobe and the full bilobal structure of transferrin have been studied under varying conditions of pH, ionic strength and binding of other metal ions by molecular dynamics (MD) simulations. These studies have allowed questions to be answered, among others, on the function of second shell residues in iron release, the role of synergistic anions in preparing the active site for iron binding, and the differences between the kinetics of the N- and the C-lobe. MD simulations on FbpA have led to the detailed observation of the binding kinetics of phosphate to the apo form, and to the conformational preferences of the holo form under conditions mimicking the environmental niches provided by the periplasmic space. To study the dynamics of these proteins with their receptors, one must resort to coarse-grained methodologies, since these systems are prohibitively large for atomistic simulations. A study of the complex of human transferrin (hTf) with its pathogenic receptor by such methods has revealed a potential mechanistic explanation for the defense mechanism that arises in evolutionary warfare. Meanwhile, the motions in the transferrin receptor bound hTf have been shown to disfavor apo hTf dissociation, explaining why the two proteins remain in complex during the recycling process from the endosome to the cell surface. Open problems and possible technological applications related to metal ion binding-release in iron transport proteins that may be handled by hybrid use of quantum mechanical, MD and coarse-grained approaches are discussed.

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.

Water exit pathways and proton pumping mechanism in B-type cytochrome c oxidase from molecular dynamics simulations

Cytochrome c oxidase (CcO) is a vital enzyme that catalyzes the reduction of molecular oxygen to water and pumps protons across mitochondrial and bacterial membranes. While proton uptake channels as well as water exit channels have been identified for A-type CcOs, the means by which water and protons exit B-type CcOs remain unclear. In this work, we investigate potential mechanisms for proton transport above the dinuclear center (DNC) in ba3-type CcO of Thermus thermophilus. Using long-time scale, all-atom molecular dynamics (MD) simulations for several relevant protonation states, we identify a potential mechanism for proton transport that involves propionate A of the active site heme a3 and residues Asp372, His376 and Glu126(II), with residue His376 acting as the proton-loading site. The proposed proton transport process involves a rotation of residue His376 and is in line with experimental findings. We also demonstrate how the strength of the salt bridge between residues Arg225 and Asp287 depends on the protonation state and that this salt bridge is unlikely to act as a simple electrostatic gate that prevents proton backflow. We identify two water exit pathways that connect the water pool above the DNC to the outer P-side of the membrane, which can potentially also act as proton exit transport pathways. Importantly, these water exit pathways can be blocked by narrowing the entrance channel between residues Gln151(II) and Arg449/Arg450 or by obstructing the entrance through a conformational change of residue Tyr136, respectively, both of which seem to be affected by protonation of residue His376.

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.

Data for molecular dynamics simulations of B-type cytochrome c oxidase with the Amber force field

Cytochrome c oxidase (CcO) is a vital enzyme that catalyzes the reduction of molecular oxygen to water and pumps protons across mitochondrial and bacterial membranes. This article presents parameters for the cofactors of ba₃-type CcO that are compatible with the all-atom Amber ff12SB and ff14SB force fields. Specifically, parameters were developed for the Cu_A pair, heme b, and the dinuclear center that consists of heme a₃ and Cu_B bridged by a hydroperoxo group. The data includes geometries in XYZ coordinate format for cluster models that were employed to compute proton transfer energies and derive bond parameters and point charges for the force field using density functional theory. Also included are the final parameter files that can be employed with the Amber leap program to generate input files for molecular dynamics simulations with the Amber software package. Based on the high resolution (1.8 Å) X-ray crystal structure of the ba₃-type CcO from Thermus thermophilus (Protein Data Bank ID number PDB: 3S8F), we built a model that is embedded in a POPC lipid bilayer membrane and solvated with TIP3P water molecules and counterions. We provide PDB data files of the initial model and the equilibrated model that can be used for further studies.

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.

Evolutionary and Functional Relationships in the Truncated Hemoglobin Family

Predicting function from sequence is an important goal in current biological research, and although, broad functional assignment is possible when a protein is assigned to a family, predicting functional specificity with accuracy is not straightforward. If function is provided by key structural properties and the relevant properties can be computed using the sequence as the starting point, it should in principle be possible to predict function in detail. The truncated hemoglobin family presents an interesting benchmark study due to their ubiquity, sequence diversity in the context of a conserved fold and the number of characterized members. Their functions are tightly related to O2 affinity and reactivity, as determined by the association and dissociation rate constants, both of which can be predicted and analyzed using in-silico based tools. In the present work we have applied a strategy, which combines homology modeling with molecular based energy calculations, to predict and analyze function of all known truncated hemoglobins in an evolutionary context. Our results show that truncated hemoglobins present conserved family features, but that its structure is flexible enough to allow the switch from high to low affinity in a few evolutionary steps. Most proteins display moderate to high oxygen affinities and multiple ligand migration paths, which, besides some minor trends, show heterogeneous distributions throughout the phylogenetic tree, again suggesting fast functional adaptation. Our data not only deepens our comprehension of the structural basis governing ligand affinity, but they also highlight some interesting functional evolutionary trends.

The N-terminal pre-A region of Mycobacterium tuberculosis 2/2HbN promotes NO-dioxygenase activity

A unique defense mechanisms by which Mycobacterium tuberculosis protects itself from nitrosative stress is based on the O2 -dependent NO-dioxygenase (NOD) activity of truncated hemoglobin 2/2HbN (Mt2/2HbN). The NOD activity largely depends on the efficiency of ligand migration to the heme cavity through a two-tunnel (long and short) system; recently, it was also correlated with the presence at the Mt2/2HbN N-terminus of a short pre-A region, not conserved in most 2/2HbNs, whose deletion results in a drastic reduction of NO scavenging. In the present study, we report the crystal structure of Mt2/2HbN-ΔpreA, lacking the pre-A region, at a resolution of 1.53 Å. We show that removal of the pre-A region results in long range effects on the protein C-terminus, promoting the assembly of a stable dimer, both in the crystals and in solution. In the Mt2/2HbN-ΔpreA dimer, access of heme ligands to the short tunnel is hindered. Molecular dynamics simulations show that the long tunnel branch is the only accessible pathway for O2 -ligand migration to/from the heme, and that the gating residue Phe(62)E15 partly restricts the diameter of the tunnel. Accordingly, kinetic measurements indicate that the kon value for peroxynitrite isomerization by Mt2/2HbN-ΔpreA-Fe(III) is four-fold lower relative to the full-length protein, and that NO scavenging by Mt2/2HbN-ΔpreA-Fe(II)-O2 is reduced by 35-fold. Therefore, we speculate that Mt2/2HbN evolved to host the pre-A region as a mechanism for preventing dimerization, thus reinforcing the survival of the microorganism against the reactive nitrosative stress in macrophages.

DATABASE

Coordinates and structure factors have been deposited in the Protein Data Bank under accession number 5AB8.

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 stability of myoglobin and glycomyoglobin: a comparative molecular dynamics simulation study

Glycoproteins are formed as the result of enzymatic glycosylation or chemical glycation in the body, and produced in vitro in industrial processes. The covalently attached carbohydrate molecule(s) confer new properties to the protein, including modified stability. In the present study, the structural stability of a glycoprotein form of myoglobin, bearing a glucose unit in the N-terminus, has been compared with its native form by the use of molecular dynamics simulation. Both structures were subjected to temperatures of 300 and 500 K in an aqueous environment for 10 ns. Changes in secondary structures and RMSD were then assessed. An overall higher stability was detected for glycomyoglobin, for which the most stable segments/residues were highlighted and compared with the native form. The simple addition of a covalently bound glucose is suggested to exert its stabilizing effect via increased contacts with surrounding water molecules, as well as a different pattern of interactions with neighbor residues.

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.

Mechanistic insight into the enzymatic reduction of truncated hemoglobin N of Mycobacterium tuberculosis: role of the CD loop and pre-A motif in electron cycling

Many pathogenic microorganisms have evolved hemoglobin-mediated nitric oxide (NO) detoxification mechanisms, where a globin domain in conjunction with a partner reductase catalyzes the conversion of toxic NO to innocuous nitrate. The truncated hemoglobin HbN of Mycobacterium tuberculosis displays a potent NO dioxygenase activity despite lacking a reductase domain. The mechanism by which HbN recycles itself during NO dioxygenation and the reductase that participates in this process are currently unknown. This study demonstrates that the NADH-ferredoxin/flavodoxin system is a fairly efficient partner for electron transfer to HbN with an observed reduction rate of 6.2 μM/min(-1), which is nearly 3- and 5-fold faster than reported for Vitreoscilla hemoglobin and myoglobin, respectively. Structural docking of the HbN with Escherichia coli NADH-flavodoxin reductase (FdR) together with site-directed mutagenesis revealed that the CD loop of the HbN forms contacts with the reductase, and that Gly(48) may have a vital role. The donor to acceptor electron coupling parameters calculated using the semiempirical pathway method amounts to an average of about 6.4 10(-5) eV, which is lower than the value obtained for E. coli flavoHb (8.0 10(-4) eV), but still supports the feasibility of an efficient electron transfer. The deletion of Pre-A abrogated the heme iron reduction by FdR in the HbN, thus signifying its involvement during intermolecular interactions of the HbN and FdR. The present study, thus, unravels a novel role of the CD loop and Pre-A motif in assisting the interactions of the HbN with the reductase and the electron cycling, which may be vital for its NO-scavenging function.

A model for the flexibility of the distal histidine in dehaloperoxidase-hemoglobin A based on X-ray crystal structures of the carbon monoxide adduct

Dehaloperoxidase hemoglobin A (DHP A) is a multifunctional hemoglobin that appears to have evolved oxidative pathways for the degradation of xenobiotics as a protective function that complements the oxygen transport function. DHP A possesses at least two internal binding sites, one for substrates and one for inhibitors, which include various halogenated phenols and indoles. Herein, we report the X-ray crystallographic structure of the carbonmonoxy complex (DHPCO). Unlike other DHP structures with 6-coordinated heme, the conformation of the distal histidine (H55) in DHPCO is primarily external or solvent exposed, despite the fact that the heme Fe is 6-coordinated. As observed generally in globins, DHP exhibits two distal histidine conformations (one internal and one external). In previous structural studies, we have shown that the distribution of H55 conformations is weighted strongly toward the external position when the DHP heme Fe is 5-coordinated. The large population of the external conformation of the distal histidine observed in DHPCO crystals at pH 6.0 indicates that some structural factor in DHP must account for the difference from other globins, which exhibit a significant external conformation only when pH < 4.5. While the original hypothesis suggested that interaction with a heme-Fe-bound ligand was the determinant of H55 conformation, the current study forces a refinement of that hypothesis. The external or open conformation of H55 is observed to have interactions with two propionate groups in heme, at distances of 3.82 and 2.73 Å, respectively. A relatively weak hydrogen bonding interaction between H55 and CO, combined with strong interactions with heme propionate (position 6), is hypothesized to strengthen the external conformation of H55. Density function theory (DFT) calculations were conducted to test whether there is a weaker hydrogen bond interaction between H55 and heme bonded CO or O₂. Molecular dynamics simulations were conducted to examine how the tautomeric forms of H55 affect the dynamic motions of the distal histidine that govern the switching between open and closed conformations. The calculations support the modified hypothesis suggesting a competition between the strength of interactions with heme ligand and the heme propionates as the factors that determine the conformation of the distal histidine.

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.