Kinetic and dynamical properties of truncated hemoglobins of the Antarctic bacterium Pseudoalteromonas haloplanktis TAC125.

Due to the low temperature, the Antarctic marine environment is challenging for protein functioning. Cold-adapted organisms have evolved proteins endowed with higher flexibility and lower stability in comparison to their thermophilic homologs, resulting in enhanced reaction rates at low temperatures. The Antarctic bacterium Pseudoalteromonas haloplanktis TAC125 (PhTAC125) genome is one of the few examples of coexistence of multiple hemoglobin genes encoding, among others, two constitutively transcribed 2/2 hemoglobins (2/2Hbs), also named truncated Hbs (TrHbs), belonging to the Group II (or O), annotated as PSHAa0030 and PSHAa2217. In this work, we describe the ligand binding kinetics and their interrelationship with the dynamical properties of globin Ph-2/2HbO-2217 by combining experimental and computational approaches and implementing a new computational method to retrieve information from molecular dynamic trajectories. We show that our approach allows us to identify docking sites within the protein matrix that are potentially able to transiently accommodate ligands and migration pathways connecting them. Consistently with ligand rebinding studies, our modeling suggests that the distal heme pocket is connected to the solvent through a low energy barrier, while inner cavities play only a minor role in modulating rebinding kinetics.

Distal lysine (de)coordination in the algal hemoglobin THB1: A combined computer simulation and experimental study.

THB1 is a monomeric truncated hemoglobin from the green alga Chlamydomonas reinhardtii. In the absence of exogenous ligands and at neutral pH, the heme group of THB1 is coordinated by two protein residues, Lys53 and His77. THB1 is thought to function as a nitric oxide dioxygenase, and the distal binding of O2 requires the cleavage of the Fe-Lys53 bond accompanied by protonation and expulsion of the lysine from the heme cavity into the solvent. Nuclear magnetic resonance spectroscopy and crystallographic data have provided dynamic and structural insights of the process, but the details of the mechanism have not been fully elucidated. We applied a combination of computer simulations and site-directed mutagenesis experiments to shed light on this issue. Molecular dynamics simulations and hybrid quantum mechanics/molecular mechanics restrained optimizations were performed to explore the nature of the transition between the decoordinated and lysine-bound states of the ferrous heme in THB1. Lys49 and Arg52, which form ionic interactions with the heme propionates in the X-ray structure of lysine-bound THB1, were observed to assist in maintaining Lys53 inside the protein cavity and play a key role in the transition. Lys49Ala, Arg52Ala and Lys49Ala/Arg52Ala THB1 variants were prepared, and the consequences of the replacements on the Lys (de)coordination equilibrium were characterized experimentally for comparison with computational prediction. The results reinforced the dynamic role of protein-propionate interactions and strongly suggested that cleavage of the Fe-Lys53 bond and ensuing conformational rearrangement is facilitated by protonation of the amino group inside the distal cavity.

Physiological Response of Escherichia coli W3110 and BL21 to the Aerobic Expression of Vitreoscilla Hemoglobin.

The aerobic growth and metabolic performance of Escherichia coli strains BL21 and W3110 were studied when the Vitreoscilla hemoglobin (VHb) was constitutively expressed in the chromosome. When VHb was expressed, acetate production decreased in both strains and was nearly eliminated in BL21. Transcriptional levels of the glyoxylate shunt genes decreased in both strains when VHb was expressed. However, higher transcription of the α-ketoglutarate dehydrogenase genes were observed for W3110, while for BL21 transcription levels decreased. VHb expression reduced the transcription of the cytochrome bo3 genes only in BL21. These results are useful for better selecting a production host.

Engineering E. coli for improved microaerobic pDNA production.

Escherichia coli strains W3110 and BL21 were engineered for the production of plasmid DNA (pDNA) under aerobic and transitions to microaerobic conditions. The gene coding for recombinase A (recA) was deleted in both strains. In addition, the Vitreoscilla hemoglobin (VHb) gene (vgb) was chromosomally inserted and constitutively expressed in each E. coli recA mutant and wild type. The recA inactivation increased the supercoiled pDNA fraction (SCF) in both strains, while VHb expression improved the pDNA production in W3110, but not in BL21. Therefore, a codon-optimized version of vgb was inserted in strain BL21recA-, which, together with W3110recA-vgb+, was tested in cultures with shifts from aerobic to oxygen-limited regimes. VHb expression lowered the accumulation of fermentative by-products in both strains. VHb-expressing cells displayed higher oxidative activity as indicated by the Redox Sensor Green fluorescence, which was more intense in BL21 than in W3110. Furthermore, VHb expression did not change pDNA production in W3110, but decreased it in BL21. These results are useful for understanding the physiological effects of VHb expression in two industrially relevant E. coli strains, and for the selection of a host for pDNA production.

Design of a synthetic miniR1 plasmid and its production by engineered Escherichia coli.

A synthetic plasmid consisting of the minimal elements for replication control of the R1 replicon and kanamycin resistance marker, which was named pminiR1, was developed. pminiR1 production was tested at 30 °C under aerobic and microaerobic conditions in Escherichia coli W3110 recA- (W1). The plasmid DNA yields from biomass (YpDNA/X) were only 0.06 ± 0.02 and 0.22 ± 0.11 mg/g under aerobic and microaerobic conditions, respectively. As an option to increase YpDNA/X values, pminiR1 was introduced in an engineered E. coli strain expressing the Vitreoscilla hemoglobin inserted in chromosome (W12). The YpDNA/X values using strain W12 increased to 0.85 ± 0.05 and 1.53 ± 0.14 mg/g under aerobic and microaerobic conditions, respectively. pminiR1 production in both strains was compared with that of pUC57Kan at 37 °C under aerobic and microaerobic conditions. The YpDNA/X values for pminiR1 using strain W12 were 6.25 ± 0.16 and 9.27 ± 0.95 mg/g under aerobic and microaerobic conditions, respectively. Such yields were similar to those obtained for plasmid pUC57Kan using strain W12 (6.9 ± 0.64 and 10.85 ± 1.06 mg/g for aerobic and microaerobic cultures, respectively). Therefore, the synthetic minimal plasmid based on the R1 replicon is a valuable alternative to pUC plasmids for biotechnological applications.

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.

Molecular basis of thermal stability in truncated (2/2) hemoglobins.

BACKGROUND: Understanding the molecular mechanism through which proteins are functional at extreme high and low temperatures is one of the key issues in structural biology. To investigate this phenomenon, we have focused on two instructive truncated hemoglobins from Thermobifida fusca (Tf-trHbO) and Mycobacterium tuberculosis (Mt-trHbO); although the two proteins are structurally nearly identical, only the former is stable at high temperatures. METHODS: We used molecular dynamics simulations at different temperatures as well as thermal melting profile measurements of both wild type proteins and two mutants designed to interchange the amino acid residue, either Pro or Gly, at E3 position. RESULTS: The results show that the presence of a Pro at the E3 position is able to increase (by 8°) or decrease (by 4°) the melting temperature of Mt-trHbO and Tf-trHbO, respectively. We observed that the ProE3 alters the structure of the CD loop, making it more flexible. CONCLUSIONS: This gain in flexibility allows the protein to concentrate its fluctuations in this single loop and avoid unfolding. The alternate conformations of the CD loop also favor the formation of more salt-bridge interactions, together augmenting the protein's thermostability. GENERAL SIGNIFICANCE: These results indicate a clear structural and dynamical role of a key residue for thermal stability in truncated hemoglobins.

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.

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

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

[Construction of conjugal transfer system of Streptomyces venezuelae var. Qinlingensis and the expression of vitreoscilla hemoglobin].

OBJECTIVE: To Construct of conjugal transfer system of Streptomyces venezuelae var. Qinlingensis and express of vitreoscilla hemoglobin. METHODS: Intergeneric genetic transfer system was based upon plasmid pSET152 and pHZ1358 from Donor E. coli ET12567 (pUZ8002) to Streptomyces venezuelae var. Qinlingensis. RESULTS: Our data showed that Intergeneric genetic transfer system was demonstrated and optimized. Integrating plasmid pJD100 containing Perm E and construction gene of vhb, constructed by SOE-PCR, was transformed to ET12567 (pUZ8002), and then transfered to Streptomyces venezuelae var. Qinlingensis by conjugation. The integrating expression of vhb gene in Streptomyces venezuelae var. Qinlingensis was verified by PCR and CO binding diference spectrum. CONCLUSION: In summary, intergeneric genetic transfer system was demonstrated and optimized, and this is the first report to express vhb gene in Streptomyces venezuelae var. Qinlingensis.

Unraveling the molecular basis for ligand binding in truncated hemoglobins: the trHbO Bacillus subtilis case.

Truncated hemoglobins (trHbs) are heme proteins present in bacteria, unicellular eukaryotes, and higher plants. Their tertiary structure consists in a 2-over-2 helical sandwich, which display typically an inner tunnel/cavity system for ligand migration and/or storage. The microorganism Bacillus subtilis contains a peculiar trHb, which does not show an evident tunnel/cavity system connecting the protein active site with the solvent, and exhibits anyway a very high oxygen association rate. Moreover, resonant Raman results of CO bound protein, showed that a complex hydrogen bond network exists in the distal cavity, making it difficult to assign unambiguously the residues involved in the stabilization of the bound ligand. To understand these experimental results with atomistic detail, we performed classical molecular dynamics simulations of the oxy, carboxy, and deoxy proteins. The free energy profiles for ligand migration suggest that there is a key residue, GlnE11, that presents an alternate conformation, in which a wide ligand migration tunnel is formed, consistently with the kinetic data. This tunnel is topologically related to the one found in group I trHbs. On the other hand, the results for the CO and O(2) bound protein show that GlnE11 is directly involved in the stabilization of the cordinated ligand, playing a similar role as TyrB10 and TrpG8 in other trHbs. Our results not only reconcile the structural data with the kinetic information, but also provide additional insight into the general behaviour of trHbs. Proteins 2010. (c) 2009 Wiley-Liss, Inc.

Spectroscopic characterization of a truncated hemoglobin from the nitrogen-fixing bacterium Herbaspirillum seropedicae.

The Herbaspirillum seropedicae genome sequence encodes a truncated hemoglobin typical of group II (Hs-trHb1) members of this family. We show that His-tagged recombinant Hs-trHb1 is monomeric in solution, and its optical spectrum resembles those of previously reported globins. NMR analysis allowed us to assign heme substituents. All data suggest that Hs-trHb1 undergoes a transition from an aquomet form in the ferric state to a hexacoordinate low-spin form in the ferrous state. The close positions of Ser-E7, Lys-E10, Tyr-B10, and His-CD1 in the distal pocket place them as candidates for heme coordination and ligand regulation. Peroxide degradation kinetics suggests an easy access to the heme pocket, as the protein offered no protection against peroxide degradation when compared with free heme. The high solvent exposure of the heme may be due to the presence of a flexible loop in the access pocket, as suggested by a structural model obtained by using homologous globins as templates. The truncated hemoglobin described here has unique features among truncated hemoglobins and may function in the facilitation of O(2) transfer and scavenging, playing an important role in the nitrogen-fixation mechanism.

Structural determinants of ligand migration in Mycobacterium tuberculosis truncated hemoglobin O.

Mycobacterium tuberculosis is the causative agent of human tuberculosis, one of the most prevalent infectious diseases in the world. Its genome hosts the glbN and glbO genes coding for two proteins, truncated hemoglobin N (trHbN) and truncated hemoglobin O (trHbO), that belong to different groups (I and II, respectively) of the recently discovered trHb family of hemeproteins. The different expression pattern and kinetics rates constants for ligand association and NO oxidation rate suggest different functions for these proteins. Previous experimental and theoretical studies showed that, in trHbs, ligand migration along the internal tunnel cavity system is a key issue in determining the ligand-binding characteristics. The X-ray structure of trHbO has been solved and shows several internal cavities and secondary-docking sites. In this work, we present an extensive investigation of the tunnel/cavity system ofM. tuberculosis trHbO by means of computer-simulation techniques. We have computed the free-energy profiles for ligand migration along three found tunnels in the oxy and deoxy w.t. and mutant trHbO proteins. Our results show that multiple-ligand migration paths are possible and that several conserved residues such as TrpG8 play a key role in the ligand-migration regulation.

Nitric oxide reactivity with globins as investigated through computer simulation.

This chapter reviews the application of classical and quantum-mechanical atomistic simulation tools used in the investigation of several relevant issues in nitric oxide reactivity with globins and presents different simulation strategies based on classical force fields: standard molecular dynamics, essential dynamics, umbrella sampling, multiple steering molecular dynamics, and a novel technique for exploring the protein energy landscape. It also presents hybrid quantum-classical schemes as a tool to obtain relevant information regarding binding energies and chemical reactivity of globins. As illustrative examples, investigations of the structural flexibility, ligand migration profiles, oxygen affinity, and reactivity toward nitric oxide of truncated hemoglobin N of Mycobacterium tuberculosis are presented.

Oxygen affinity controlled by dynamical distal conformations: the soybean leghemoglobin and the Paramecium caudatum hemoglobin cases.

The binding of diatomic ligands, such as O(2), NO, and CO, to heme proteins is a process intimately related with their function. In this work, we analyzed by means of a combination of classical Molecular Dynamics (MD) and Hybrid Quantum-Classical (QM/MM) techniques the existence of multiple conformations in the distal site of heme proteins and their influence on oxygen affinity regulation. We considered two representative examples: soybean leghemoglobin (Lba) and Paramecium caudatum truncated hemoglobin (PcHb). The results presented in this work provide a molecular interpretation for the kinetic, structural, and mutational data that cannot be obtained by assuming a single distal conformation.

Theoretical study of the truncated hemoglobin HbN: exploring the molecular basis of the NO detoxification mechanism.

Mycobacterium tuberculosis is the causative agent of human tuberculosis. The nitric oxide reaction with oxy-truncated hemoglobin N (trHbN) has been proposed to be responsible for the resistance mechanism by which this microorganism can evade the toxic effects of NO. In this work, we explore the molecular basis of the NO detoxification mechanism using a combination of classical and hybrid quantum-classical (QM-MM) simulation techniques. We have investigated the structural flexibility of the protein, the ligand affinity properties, and the nitric oxide reaction with coordinated O2. The analysis of the classical MD trajectory allowed us to identify Phe62 as the gate of the main channel for ligand diffusion to the active site. Moreover, the opening of the channel stems from the interplay between collective backbone motions and local rearrangements in the side chains of the residues that form the bottleneck of the tunnel. Even though the protein environment is not found to make a significant contribution to the heme moiety catalyzed reaction, the binding site influences the physiological function of the enzyme at three different levels. First, by isolating the intermediates formed in the reaction, it prevents nondesired reactions from proceeding. Second, it modulates the ligand (O2, NO) affinity of the protein, which can be ascribed to both distal and proximal effects. Finally, the stabilization of the Tyr33-Gln58 pair upon O2 binding might alter the essential dynamics of the protein, leading in turn to a mechanism for ligand-induced regulation.