Ligand migration through hemeprotein cavities: insights from laser flash photolysis and molecular dynamics simulations

Unusually Fast bis-Histidyl Coordination in a Plant Hemoglobin

The recently identified nonsymbiotic hemoglobin gene MtGlb1-2 of the legume Medicago truncatula possesses unique properties as it generates four alternative splice forms encoding proteins with one or two heme domains. Here we investigate the ligand binding kinetics of MtGlb1-2.1 and MtGlb1-2.4, bearing two hemes and one heme, respectively. Unexpectedly, the overall time-course of ligand rebinding was unusually fast. Thus, we complemented nanosecond laser flash photolysis kinetics with data collected with a hybrid femtosecond–nanosecond pump–probe setup. Most photodissociated ligands are rebound geminately within a few nanoseconds, which leads to rates of the bimolecular rebinding to pentacoordinate species in the 10⁸ M⁻¹s⁻¹ range. Binding of the distal histidine to the heme competes with CO rebinding with extremely high rates (k_h ~ 10⁵ s⁻¹). Histidine dissociation from the heme occurs with comparable rates, thus resulting in moderate equilibrium binding constants (K_H ~ 1). The rate constants for ligation and deligation of distal histidine to the heme are the highest reported for any plant or vertebrate globin. The combination of microscopic rates results in unusually high overall ligand binding rate constants, a fact that contributes to explaining at the mechanistic level the extremely high reactivity of these proteins toward the physiological ligands oxygen, nitric oxide and nitrite.

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.

Synthetic Fe/Cu Complexes: Toward Understanding Heme-Copper Oxidase Structure and Function

Heme-copper oxidases (HCOs) are terminal enzymes on the mitochondrial or bacterial respiratory electron transport chain, which utilize a unique heterobinuclear active site to catalyze the 4H+/4e- reduction of dioxygen to water. This process involves a proton-coupled electron transfer (PCET) from a tyrosine (phenolic) residue and additional redox events coupled to transmembrane proton pumping and ATP synthesis. Given that HCOs are large, complex, membrane-bound enzymes, bioinspired synthetic model chemistry is a promising approach to better understand heme-Cu-mediated dioxygen reduction, including the details of proton and electron movements. This review encompasses important aspects of heme-O2 and copper-O2 (bio)chemistries as they relate to the design and interpretation of small molecule model systems and provides perspectives from fundamental coordination chemistry, which can be applied to the understanding of HCO activity. We focus on recent advancements from studies of heme-Cu models, evaluating experimental and computational results, which highlight important fundamental structure-function relationships. Finally, we provide an outlook for future potential contributions from synthetic inorganic chemistry and discuss their implications with relevance to biological O2-reduction.

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.

Electrostatic Tuning of the Ligand Binding Mechanism by Glu27 in Nitrophorin 7

Nitrophorins (NP) 1-7 are NO-carrying heme proteins found in the saliva of the blood-sucking insect Rhodnius prolixus. The isoform NP7 displays peculiar properties, such as an abnormally high isoelectric point, the ability to bind negatively charged membranes, and a strong pH sensitivity of NO affinity. A unique trait of NP7 is the presence of Glu in position 27, which is occupied by Val in other NPs. Glu27 appears to be important for tuning the heme properties, but its influence on the pH-dependent NO release mechanism, which is assisted by a conformational change in the AB loop, remains unexplored. Here, in order to gain insight into the functional role of Glu27, we examine the effect of Glu27 → Val and Glu27 → Gln mutations on the ligand binding kinetics using CO as a model. The results reveal that annihilation of the negative charge of Glu27 upon mutation reduces the pH sensitivity of the ligand binding rate, a process that in turn depends on the ionization of Asp32. We propose that Glu27 exerts a through-space electrostatic action on Asp32, which shifts the pKa of the latter amino acid towards more acidic values thus reducing the pH sensitivity of the transition between open and closed states.

Characterization of the Heme Pocket Structure and Ligand Binding Kinetics of Non-symbiotic Hemoglobins from the Model Legume Lotus japonicus

Plant hemoglobins (Hbs) are found in nodules of legumes and actinorhizal plants but also in non-symbiotic organs of monocots and dicots. Non-symbiotic Hbs (nsHbs) have been classified into two phylogenetic groups. Class 1 nsHbs show an extremely high O2 affinity and are induced by hypoxia and nitric oxide (NO), whereas class 2 nsHbs have moderate O2 affinity and are induced by cold and cytokinins. The functions of nsHbs are still unclear, but some of them rely on the capacity of hemes to bind diatomic ligands and catalyze the NO dioxygenase (NOD) reaction (oxyferrous Hb + NO → ferric Hb + nitrate). Moreover, NO may nitrosylate Cys residues of proteins. It is therefore important to determine the ligand binding properties of the hemes and the role of Cys residues. Here, we have addressed these issues with the two class 1 nsHbs (LjGlb1-1 and LjGlb1-2) and the single class 2 nsHb (LjGlb2) of Lotus japonicus, which is a model legume used to facilitate the transfer of genetic and biochemical information into crops. We have employed carbon monoxide (CO) as a model ligand and resonance Raman, laser flash photolysis, and stopped-flow spectroscopies to unveil major differences in the heme environments and ligand binding kinetics of the three proteins, which suggest non-redundant functions. In the deoxyferrous state, LjGlb1-1 is partially hexacoordinate, whereas LjGlb1-2 shows complete hexacoordination (behaving like class 2 nsHbs) and LjGlb2 is mostly pentacoordinate (unlike other class 2 nsHbs). LjGlb1-1 binds CO very strongly by stabilizing it through hydrogen bonding, but LjGlb1-2 and LjGlb2 show lower CO stabilization. The changes in CO stabilization would explain the different affinities of the three proteins for gaseous ligands. These affinities are determined by the dissociation rates and follow the order LjGlb1-1 > LjGlb1-2 > LjGlb2. Mutations LjGlb1-1 C78S and LjGlb1-2 C79S caused important alterations in protein dynamics and stability, indicating a structural role of those Cys residues, whereas mutation LjGlb1-1 C8S had a smaller effect. The three proteins and their mutant derivatives exhibited similarly high rates of NO consumption, which were due to NOD activity of the hemes and not to nitrosylation of Cys residues.

Photoacoustic calorimetry studies of CO photo-dissociation from chloramine-T modified horse heart cytochrome-c

Treatment of horse heart Cytochrome-c (Cc) with N-chloro-4-toluosulfonamide (Chloramine-t, CT) results in the oxidation of methionine (Met) residues to the corresponding sulfoxide including the distal heme ligand, Met80. The resulting Fe-sulfoxide coordination is sufficiently labile in the ferrous form to be displaced by gaseous ligands, including CO. Photolysis of the CO-CT-Cc complex provides an opportunity to examine ligand binding dynamics that are associated with a relatively rigid distal heme pocket. In this work, photoacoustic calorimetry (PAC) was utilized to obtain the kinetics as well as enthalpy and molar volume changes subsequent to CO photo-dissociation from CO-CT-Cc. Previous photolysis studies of CO-CT-Cc have led to a proposed model for ligand recombination in which the Met80-sulfoxide and CO recombine with the heme on relatively slow timescales (50 μs and ∼500 μs, respectively). The PAC data presented here reveals two additional kinetic phases with lifetimes of <20 ns and 534 ± 75 ns. The fast phase (<20 ns) is associated with an ΔH of 44 ± 5 kcal mol^-1 and ΔV of -0.5 ± 0.5 mL mol^-1, whereas the slower phase (534 ns) is associated with a small ΔH of 2 ± 3 kcal mol^-1 and ΔV of 1 ± 0.5 mL mol^-1.

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.

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.

Docking and migration of carbon monoxide in nitrogenase: the case for gated pockets from infrared spectroscopy and molecular dynamics

Evidence of a CO docking site near the FeMo cofactor in nitrogenase has been obtained by Fourier transform infrared spectroscopy-monitored low-temperature photolysis. We investigated the possible migration paths for CO from this docking site using molecular dynamics calculations. The simulations support the notion of a gas channel with multiple internal pockets from the active site to the protein exterior. Travel between pockets is gated by the motion of protein residues. Implications for the mechanism of nitrogenase reactions with CO and N2 are discussed.

Structure and dynamics of the membrane attaching nitric oxide transporter nitrophorin 7

Nitrophorins represent a unique class of heme proteins that are able to perform the delicate transportation and release of the free-radical gaseous messenger nitric oxide (NO) in a pH-triggered manner. Besides its ability to bind to phospholipid membranes, the N-terminus contains an additional Leu-Pro-Gly stretch, which is a unique sequence trait, and the heme cavity is significantly altered with respect to other nitrophorins. These distinctive features encouraged us to solve the X-ray crystallographic structures of NP7 at low and high pH and bound with different heme ligands (nitric oxide, histamine, imidazole). The overall fold of the lipocalin motif is well preserved in the different X-ray structures and resembles the fold of other nitrophorins. However, a chain-like arrangement in the crystal lattice due to a number of head-to-tail electrostatic stabilizing interactions is found in NP7. Furthermore, the X-ray structures also reveal ligand-dependent changes in the orientation of the heme, as well as in specific interactions between the A-B and G-H loops, which are considered to be relevant for the biological function of nitrophorins. Fast and ultrafast laser triggered ligand rebinding experiments demonstrate the pH-dependent ligand migration within the cavities and the exit route. Finally, the topological distribution of pockets located around the heme as well as from inner cavities present at the rear of the protein provides a distinctive feature in NP7, so that while a loop gated exit mechanism to the solvent has been proposed for most nitrophorins, a more complex mechanism that involves several interconnected gas hosting cavities is proposed for NP7.

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.