pH-dependent mechanism of nitric oxide release in nitrophorins 2 and 4

A Fast and Interpretable Deep Learning Approach for Accurate Electrostatics-Driven pKa Predictions in Proteins

Existing computational methods for estimating pK_a values in proteins rely on theoretical approximations and lengthy computations. In this work, we use a data set of 6 million theoretically determined pK_a shifts to train deep learning models, which are shown to rival the physics-based predictors. These neural networks managed to infer the electrostatic contributions of different chemical groups and learned the importance of solvent exposure and close interactions, including hydrogen bonds. Although trained only using theoretical data, our pKAI+ model displayed the best accuracy in a test set of ∼750 experimental values. Inference times allow speedups of more than 1000× compared to physics-based methods. By combining speed, accuracy, and a reasonable understanding of the underlying physics, our models provide a game-changing solution for fast estimations of macroscopic pK_a values from ensembles of microscopic values as well as for many downstream applications such as molecular docking and constant-pH molecular dynamics simulations.

The Biologically Relevant Coordination Chemistry of Iron and Nitric Oxide: Electronic Structure and Reactivity

Nitric oxide (NO) is an important signaling molecule that is involved in a wide range of physiological and pathological events in biology. Metal coordination chemistry, especially with iron, is at the heart of many biological transformations involving NO. A series of heme proteins, nitric oxide synthases (NOS), soluble guanylate cyclase (sGC), and nitrophorins, are responsible for the biosynthesis, sensing, and transport of NO. Alternatively, NO can be generated from nitrite by heme- and copper-containing nitrite reductases (NIRs). The NO-bearing small molecules such as nitrosothiols and dinitrosyl iron complexes (DNICs) can serve as an alternative vehicle for NO storage and transport. Once NO is formed, the rich reaction chemistry of NO leads to a wide variety of biological activities including reduction of NO by heme or non-heme iron-containing NO reductases and protein post-translational modifications by DNICs. Much of our understanding of the reactivity of metal sites in biology with NO and the mechanisms of these transformations has come from the elucidation of the geometric and electronic structures and chemical reactivity of synthetic model systems, in synergy with biochemical and biophysical studies on the relevant proteins themselves. This review focuses on recent advancements from studies on proteins and model complexes that not only have improved our understanding of the biological roles of NO but also have provided foundations for biomedical research and for bio-inspired catalyst design in energy science.

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.

Use of gasotransmitters for the controlled release of polymer-based nitric oxide carriers in medical applications

Nitric Oxide (NO) is a small molecule gasotransmitter synthesized by nitric oxide synthase in almost all types of mammalian cells. NO is synthesized by NO synthase by conversion of l-arginine to l-citrulline in the human body. NO then stimulates soluble guanylate cyclase, from which various physiological functions are mediated in a concentration-dependent manner. High concentrations of NO induce apoptosis or antibacterial responses whereas low NO circulation leads to angiogenesis. The bidirectional effect of NO has attracted considerable attention, and efforts to deliver NO in a controlled manner, especially through polymeric carriers, has been the topic of much research. This naturally produced signaling molecule has stood out as a potentially more potent therapeutic agent compared to exogenously synthesized drugs. In this review, we will focus on past efforts of using the controlled release of NO via polymer-based materials to derive specific therapeutic results. We have also added studies and our future suggestions on co-delivery methods with other gasotransmitters as a step towards developing multifunctional carriers.

Nitrophorins and nitrobindins: structure and function

Classical all α-helical globins are present in all living organisms and are ordered in three lineages: (i) flavohemoglobins and single domain globins, (ii) protoglobins and globin coupled sensors and (iii) truncated hemoglobins, displaying the 3/3 or the 2/2 all α-helical fold. However, over the last two decades, all β-barrel and mixed α-helical-β-barrel heme-proteins displaying heme-based functional properties (e.g. ligand binding, transport and sensing) closely similar to those of all α-helical globins have been reported. Monomeric nitrophorins (NPs) and α1-microglobulin (α1-m), belonging to the lipocalin superfamily and nitrobindins (Nbs) represent prototypical heme-proteins displaying the all β-barrel and mixed α-helical-β-barrel folds. NPs are confined to the Reduviidae and Cimicidae families of Heteroptera, whereas α1-m and Nbs constitute heme-protein families spanning bacteria to Homo sapiens. The structural organization and the reactivity of the stable ferric solvent-exposed heme-Fe atom suggest that NPs and Nbs are devoted to NO transport, storage and sensing, whereas Hs-α1-m participates in heme metabolism. Here, the structural and functional properties of NPs and Nbs are reviewed in parallel with those of sperm whale myoglobin, which is generally taken as the prototype of monomeric globins.

Ionization basis for activation of enzymes soluble in ionic liquids

BACKGROUND

The complex interactions between electrolytes and proteins have been studied for more than a century. However, understanding is not yet complete and does not provide a basis for predicting the activity of enzymes in ionic media. The use of ionic liquids (ILs) as reaction medium has opened up new opportunities for better understanding of the mechanism of enzymatic catalysis. Although a number of properties of ILs have been correlated with enzyme function, these relationships are not completely understood at a molecular level.

METHODS

We propose that ILs must be able to promote ionization of protein ionizable groups in order to dissolve active enzymes. The biocompatible IL need to possess a functional group with large donor number and acceptor number in both cationic and anionic units, each of which is based on a high dielectric constant lead structure. We designed and synthesized two series of ILs and determined their ionizing-dissociating abilities and activities of lipases soluble in these new ILs.

RESULTS

The results showed that the ionizing-dissociating abilities of ILs paralleled the catalytic activity trend of lipases dissolved in the ILs. The activities of lipases soluble in the newly designed ILs were comparable to those in water.

CONCLUSIONS

We can conclude that ionizing-dissociating abilities of an IL can be used as a basis for predicting the activity of enzymes soluble in the IL. General significance Ionization basis for activation of enzymes gives a deeper understanding of the behavior of enzymes in non-aqueous media at a molecular level.

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.

Underlying thermodynamics of pH-dependent allostery

Understanding the effects of coupling protein protonation and conformational states is critical to the development of drugs targeting pH sensors and to the rational engineering of pH switches. In this work, we address this issue by performing a comprehensive study of the pH-regulated switch from the closed to the open conformation in nitrophorin 4 (NP4) that determines its pH-dependent activity. Our calculations show that D30 is the only amino acid that has two significantly different pKas in the open and closed conformations, confirming its critical role in regulating pH-dependent behavior. In addition, we describe the free-energy landscape of the conformational change as a function of pH, obtaining accurate estimations of free-energy barriers and equilibrium constants using different methods. The underlying thermodynamic model of the switch workings suggests the possibility of tuning the observed pKa only through the conformational equilibria, keeping the same conformation-specific pKas, as evidenced by the proposed K125L mutant. Moreover, coupling between the protonation and conformational equilibria results in efficient regulation and pH-sensing around physiological pH values only for some combinations of protonation and conformational equilibrium constants, placing constraints on their possible values and leaving a narrow space for protein molecular evolution. The calculations and analysis presented here are of general applicability and provide a guide as to how more complex systems can be studied, offering insight into how pH-regulated allostery works of great value for designing drugs that target pH sensors and for rational engineering of pH switches beyond the common histidine trigger.

Recent development and application of constant pH molecular dynamics

Solution pH is a critical environmental factor for chemical and biological processes. Over the last decade, significant efforts have been made in the development of constant pH molecular dynamics (pHMD) techniques for gaining detailed insights into pH-coupled dynamical phenomena. In this article we review the advancement of this field in the past five years, placing a special emphasis on the development of the all-atom continuous pHMD technique. We discuss various applications, including the prediction of pK_a shifts for proteins, nucleic acids and surfactant assemblies, elucidation of pH-dependent population shifts, protein-protein and protein-RNA binding, as well as the mechanisms of pH-dependent self-assembly and phase transitions of surfactants and peptides. We also discuss future directions for the further improvement of the pHMD techniques.

NMR studies of the dynamics of nitrophorin 2 bound to nitric oxide

The Rhodnius nitrophorins are β-barrel proteins of the lipocalin fold with a heme protruding from the open end of the barrel. They are found in the saliva of the blood-sucking insect Rhodnius prolixus, which synthesizes and stores nitric oxide (NO) in the salivary glands, where NO is bound to iron. NO is released by dilution and an increase in pH when the insect spits its saliva into the tissues of a victim, to aid in obtaining a blood meal. In the adult insect, there are four nitrophorins, NP1-NP4. At pH 7.3, NP4 releases NO 17 times faster than NP2 does, as measured by stopped-flow kinetics. A number of crystal structures of the least abundant protein, NP4, are available. These structures have been used to propose that two loops between adjacent β-strands at the front opening of the protein, the A-B and G-H loops, determine the rate of NO release. To learn how the protein loops contribute to the release of NO for each of the nitrophorins, the dynamics of these proteins are being studied in our laboratory. In this work, the NP2-NO complex has been investigated by nuclear magnetic resonance relaxation measurements to probe the picosecond-to-nanosecond and microsecond-to-millisecond time scale motions at three pH values, 5.0, 6.5, and 7.3. It is found that at pH 5.0 and 6.5, the NP2-NO complex is rigid and only a few residues in the loop regions show dynamics, while at pH 7.3, somewhat more dynamics, particularly of the A-B loop, are observed. Comparison to other lipocalins shows that all are relatively rigid, and that the dynamics of lipocalins in general are much more subtle than those of mainly α-helical proteins.

pH-dependent picosecond structural dynamics in the distal pocket of nitrophorin 4 investigated by 2D IR spectroscopy

Nitrophorin 4 (NP4) belongs to a family of pH-sensitive, nitric oxide (NO) transporter proteins that undergo a large structural change from a closed to an open conformation at high pH to allow for NO delivery. Measuring the pH-dependent structural dynamics in NP4-NO around the ligand binding site is crucial for developing a mechanistic understanding of NO binding and release. In this study, we use coherent two-dimensional infrared (2D IR) spectroscopy to measure picosecond structural dynamics sampled by the nitrosyl stretch in NP4-NO as a function of pH at room temperature. Our results show that both the closed and open conformers of the protein are present at low (pD 5.1) and high (pD 7.9) pH conditions. The closed and open conformers are characterized by two frequencies of the nitrosyl stretching vibration labeled A0 and A1, respectively. Analysis of the 2D IR line shapes reveals that at pD 5.1, the closed conformer experiences structural fluctuations arising from solvation dynamics on a ∼3 ps time scale. At pD 7.9, both the open and closed conformers exhibit fluctuations on a ∼1 ps time scale. At both pD conditions, the closed conformers maintain a static distribution of structures within the experimental time window of 100 ps. This is in contrast to the open conformer, which is able to interconvert among its substates on a ∼100 ps time scale. Our results directly measure the time scales of solvation dynamics in the distal pocket, the flexibility of the open conformation at high pH, and the rigidity of the closed conformers at both pH conditions. We discuss how the pH-dependent equilibrium structural fluctuations of the nitrosyl ligand measured in this study are related to the uptake and delivery of nitric oxide in NP4.

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.

Protonation and pK changes in protein-ligand binding

Formation of protein-ligand complexes causes various changes in both the receptor and the ligand. This review focuses on changes in pK and protonation states of ionizable groups that accompany protein-ligand binding. Physical origins of these effects are outlined, followed by a brief overview of the computational methods to predict them and the associated corrections to receptor-ligand binding affinities. Statistical prevalence, magnitude and spatial distribution of the pK and protonation state changes in protein-ligand binding are discussed in detail, based on both experimental and theoretical studies. While there is no doubt that these changes occur, they do not occur all the time; the estimated prevalence varies, both between individual complexes and by method. The changes occur not only in the immediate vicinity of the interface but also sometimes far away. When receptor-ligand binding is associated with protonation state change at particular pH, the binding becomes pH dependent: we review the interplay between sub-cellular characteristic pH and optimum pH of receptor-ligand binding. It is pointed out that there is a tendency for protonation state changes upon binding to be minimal at physiologically relevant pH for each complex (no net proton uptake/release), suggesting that native receptor-ligand interactions have evolved to reduce the energy cost associated with ionization changes. As a result, previously reported statistical prevalence of these changes - typically computed at the same pH for all complexes - may be higher than what may be expected at optimum pH specific to each complex. We also discuss whether proper account of protonation state changes appears to improve practical docking and scoring outcomes relevant to structure-based drug design. An overview of some of the existing challenges in the field is provided in conclusion.

Kinetics and computational studies of ligand migration in nitrophorin 7 and its Δ1-3 mutant

Nitrophorins (NPs) are nitric oxide (NO)-carrying heme proteins found in the saliva of the blood-sucking insect Rhodnius prolixus. Though NP7 exhibits a large sequence resemblance with other NPs, two major differential features are the ability to interact with negatively charged cell surfaces and the presence of a specific N-terminus composed of three extra residues (Leu1-Pro2-Gly3). The aim of this study is to examine the influence of the N-terminus on the ligand binding, and the topological features of inner cavities in closed and open states of NP7, which can be associated to the protein structure at low and high pH, respectively. Laser flash photolysis measurements of the CO rebinding kinetics to NP7 and its variant NP7(Δ1-3), which lacks the three extra residues at the N-terminus, exhibit a similar pattern and support the existence of a common kinetic mechanism for ligand migration and binding. This is supported by the existence of a common topology of inner cavities, which consists of two docking sites in the heme pocket and a secondary site at the back of the protein. The ligand exchange between these cavities is facilitated by an additional site, which can be transiently occupied by the ligand in NP7, although it is absent in NP4. These features provide a basis to explain the enhanced internal gas hosting capacity found experimentally in NP7 and the absence of ligand rebinding from secondary sites in NP4. The current data allow us to speculate that the processes of docking to cell surfaces and NO release may be interconnected in NP7, thereby efficiently releasing NO into a target cell. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.

Comparative study of enzyme activity and heme reactivity in Drosophila melanogaster and Homo sapiens cystathionine β-synthases

Cystathionine β-synthase (CBS) is the first and rate-limiting enzyme in the transsulfuration pathway, which is critical for the synthesis of cysteine from methionine in eukaryotes. CBS uses coenzyme pyridoxal 5'-phosphate (PLP) for catalysis, and S-adenosylmethionine regulates the activity of human CBS, but not yeast CBS. Human and fruit fly CBS contain heme; however, the role for heme is not clear. This paper reports biochemical and spectroscopic characterization of CBS from fruit fly Drosophila melanogaster (DmCBS) and the CO/NO gas binding reactions of DmCBS and human CBS. Like CBS enzymes from lower organisms (e.g., yeast), DmCBS is intrinsically highly active and is not regulated by AdoMet. The DmCBS heme coordination environment, the reactivity, and the accompanying effects on enzyme activity are similar to those of human CBS. The DmCBS heme bears histidine and cysteine axial ligands, and the enzyme becomes inactive when the cysteine ligand is replaced. The Fe(II) heme in DmCBS is less stable than that in human CBS, undergoing more facile reoxidation and ligand exchange. In both CBS proteins, the overall stability of the protein is correlated with the heme oxidation state. Human and DmCBS Fe(II) hemes react relatively slowly with CO and NO, and the rate of the CO binding reaction is faster at low pH than at high pH. Together, the results suggest that heme incorporation and AdoMet regulation in CBS are not correlated, possibly providing two independent means for regulating the enzyme.

pH-Dependent conformational changes in proteins and their effect on experimental pK(a)s: the case of Nitrophorin 4

The acid-base behavior of amino acids is an important subject of study due to their prominent role in enzyme catalysis, substrate binding and protein structure. Due to interactions with the protein environment, their pK_as can be shifted from their solution values and, if a protein has two stable conformations, it is possible for a residue to have different “microscopic”, conformation-dependent pK_a values. In those cases, interpretation of experimental measurements of the pK_a is complicated by the coupling between pH, protonation state and protein conformation. We explored these issues using Nitrophorin 4 (NP4), a protein that releases NO in a pH sensitive manner. At pH 5.5 NP4 is in a closed conformation where NO is tightly bound, while at pH 7.5 Asp30 becomes deprotonated, causing the conformation to change to an open state from which NO can easily escape. Using constant pH molecular dynamics we found two distinct microscopic Asp30 pK_as: 8.5 in the closed structure and 4.3 in the open structure. Using a four-state model, we then related the obtained microscopic values to the experimentally observed “apparent” pK_a, obtaining a value of 6.5, in excellent agreement with experimental data. This value must be interpreted as the pH at which the closed to open population transition takes place. More generally, our results show that it is possible to relate microscopic structure dependent pKa values to experimentally observed ensemble dependent apparent pK_as and that the insight gained in the relatively simple case of NP4 can be useful in several more complex cases involving a pH dependent transition, of great biochemical interest.

The interaction of an amino acid with its protein environment can result in an acid-base behavior that is very different from what would be observed in solution. This environment can be greatly altered when the protein changes conformation. As a result, the amino acid will have two different “microscopic” pK_a values. Nitrophorin 4 is a good case study to explore this behavior, because it undergoes a pH-dependent conformational change that is well characterized experimentally. Using computer simulation tools, we found that the key titratable Aspartic acid 30, has two very different microscopic pK_as: 4.3 and 8.5, which are significantly different to the observed transition pK_a in solution. However, using a simple model, we were able to understand how this causes the conformational change to take place at pH∼6.5, as measured experimentally. The insight gained in this relatively simple case can be useful in other more complex cases where the apparent pK_a is also a result of the interplay of different conformations where some amino acids experience very different environments.

Assignment of the 1H NMR resonances of protein residues in close proximity to the heme of the nitrophorins: similarities and differences among the four proteins from the saliva of the adult blood-sucking insect Rhodnius prolixus

The nuclear Overhauser effects (NOEs) observed between heme substituent protons and a small number of nearby protein side chain protons in the water-elimination Fourier transform NOE spectroscopy (WEFT-NOESY) spectra of high- and low-spin wild-type nitrophorin (NP) 2 and its ligand complexes have been analyzed and compared with those observed for the same complexes of wild-type NP3. These assignments were made on naturally abundant isotope samples, with the most useful protein side chains being those of Ile120, Leu122, and Leu132 for NP2 and NP3, and Thr121, Leu123, and Leu133 for NP1 and NP4. It is found that the NOEs observed are identical, with extremely similar protein side chain proton chemical shifts. This is strong evidence that the structure of NP3, for which no X-ray crystal structures are available, is essentially identical to that of NP2, at least near the heme binding pocket. Similarly, the NOEs observed between heme substituents and protein side chains for NP1 and NP4 also indicate that the structures of the protein having both A and B heme orientations are very similar to each other, as well as to the proteins with major B heme orientation of NP2 and NP3. These A and B connectivities can be seen, even though the two heme orientations have similar populations in NP1 and NP4, which complicates the analysis of the NOESY spectra. The histamine complex of wild-type NP2 shows significant shifts of the Leu132 side chain protons relative to all other ligand complexes of NP1-NP4 because of the perturbation of the structure near Leu132 caused by the histamine's side chain ammonium hydrogen bond to the Asp29 side chain carboxylate.

Heterogeneous kinetics of the carbon monoxide association and dissociation reaction to nitrophorin 4 and 7 coincide with structural heterogeneity of the gate-loop

NO is an important signaling molecule in human tissue. However, the mechanisms by which this molecule is controlled and directed are currently little understood. Nitrophorins (NPs) comprise a group of ferriheme proteins originating from blood-sucking insects that are tailored to protect and deliver NO via coordination to and release from the heme iron. Therefore, the kinetics of the association and dissociation reactions were studied in this work using the ferroheme-CO complexes of NP4, NP4(D30N), and NP7 as isoelectronic models for the ferriheme-NO complexes. The kinetic measurements performed by nanosecond laser-flash-photolysis and stopped-flow are accompanied by resonance Raman and FT-IR spectroscopy to characterize the carbonyl species. Careful analysis of the CO rebinding kinetics reveals that in NP4 and, to a larger extent, NP7 internal gas binding cavities are located, which temporarily trap photodissociated ligands. Moreover, changes in the free energy barriers throughout the rebinding and release pathway upon increase of the pH are surprisingly small in case of NP4. Also in case of NP4, a heterogeneous kinetic trace is obtained at pH 7.5, which corresponds to the presence of two carbonyl species in the heme cavity that are seen in vibrational spectroscopy and that are due to the change of the distal heme pocket polarity. Quantification of the two species from FT-IR spectra allowed the fitting of the kinetic traces as two processes, corresponding to the previously reported open and closed conformation of the A-B and G-H loops. With the use of the A-B loop mutant NP4(D30N), it was confirmed that the kinetic heterogeneity is controlled by pH through the disruption of the H-bond between the Asp30 side chain and the Leu130 backbone carbonyl. Overall, this first study on the slow phase of the dynamics of diatomic gas molecule interaction with NPs comprises an important experimental contribution for the understanding of the dynamics involved in the binding/release processes of NO/CO in NPs.

Carbon Dioxide Migration Pathways in Proteins

Some of the most important biological processes, such as carbon fixation, are dependent on protein-gas interactions. The motion of CO2 through the enzyme phosphoenolpyruvate carboxykinase was investigated using extensive all-atom molecular dynamics simulations. Three discrete migration pathways were located, suggesting the protein directs the movement of CO2. The chemical nature of these pathways is discussed, as are their biotechnological ramifications.

Nuclear inelastic scattering and Mössbauer spectroscopy as local probes for ligand binding modes and electronic properties in proteins: vibrational behavior of a ferriheme center inside a β-barrel protein

In this work, we present a study of the influence of the protein matrix on its ability to tune the binding of small ligands such as NO, cyanide (CN(-)), and histamine to the ferric heme iron center in the NO-storage and -transport protein Nitrophorin 2 (NP2) from the salivary glands of the blood-sucking insect Rhodnius prolixus. Conventional Mössbauer spectroscopy shows a diamagnetic ground state of the NP2-NO complex and Type I and II electronic ground states of the NP2-CN(-) and NP2-histamine complex, respectively. The change in the vibrational signature of the protein upon ligand binding has been monitored by Nuclear Inelastic Scattering (NIS), also called Nuclear Resonant Vibrational Spectroscopy (NRVS). The NIS data thus obtained have also been calculated by quantum mechanical (QM) density functional theory (DFT) coupled with molecular mechanics (MM) methods. The calculations presented here show that the heme ruffling in NP2 is a consequence of the interaction with the protein matrix. Structure optimizations of the heme and its ligands with DFT retain the characteristic saddling and ruffling only if the protein matrix is taken into account. Furthermore, simulations of the NIS data by QM/MM calculations suggest that the pH dependence of the binding of NO, but not of CN(-) and histamine, might be a consequence of the protonation state of the heme carboxyls.

NMR studies of nitrophorin distal pocket side chain effects on the heme orientation and seating of NP2 as compared to NP1

The nitrophorins (NP) of the adult blood-sucking insect Rhodnius prolixus fall into two pairs based on sequence identity (NP1,4 (90%) and NP2,3 (79%)), which differ significantly in the size of side chains of residues which contact the heme. These residues include those in the distal pocket of NP2 (I120) and NP1 (T121) and the "belt" that surrounds the heme of NP2 (S40, F42), and NP1(A42, L44). To determine the importance of these residues and others conserved or very similar for the two pairs, including L122(123), L132(133), appropriate mutants of NP2 and NP1 have been prepared and studied by (1)H NMR spectroscopy. Wild-type NP2 has heme orientation ratio (A:B) of 1:8 at equilibrium, while wild-type NP1 has A:B ~1:1 at equilibrium. Another difference between NP2 and NP1 is in the heme seating with regard to His57(59). It is found that among the distal pocket residues investigated, the residue most responsible for heme orientation and seating is I120(T121). F42(L44) and L106(F107) may also be important, but must be investigated in greater detail.

Histidine-mediated pH-sensitive regulation of M-ficolin:GlcNAc binding activity in innate immunity examined by molecular dynamics simulations

BACKGROUND

M-ficolin, a pathogen recognition molecule in the innate immune system, binds sugar residues including N-acetyl-D-glucosamine (GlcNAc), which is displayed on invading microbes and on apoptotic cells. The cis and trans Asp282-Cys283 peptide bond in the M-ficolin, which was found to occur at neutral and acidic pH in crystal structures, has been suggested to represent binding and non-binding activity, respectively. A detailed understanding of the pH-dependent conformational changes in M-ficolin and pH-mediated discrimination mechanism of GlcNAc-binding activity are crucial to both immune-surveillance and clearance of apoptotic cells.

METHODOLOGY/PRINCIPAL FINDINGS

By immunodetection analysis, we found that the pH-sensitive binding of GlcNAc is regulated by a conformational equilibrium between the active and inactive states of M-ficolin. We performed constant pH molecular dynamics (MD) simulation at a series of pH values to explore the pH effect on the cis-trans isomerization of the Asp282-Cys283 peptide bond in the M-ficolin fibrinogen-like domain (FBG). Analysis of the hydrogen bond occupancy of wild type FBG compared with three His mutants (H251A, H284A and H297A) corroborates that His284 is indispensible for pH-dependent binding. H251A formed new but weaker hydrogen bonds with GlcNAc. His297, unlike the other two His mutants, is more dependent on the solution pH and also contributes to cis-trans isomerization of the Asp282-Cys283 peptide bond in weak basic solution.

CONCLUSIONS/SIGNIFICANCE

Constant pH MD simulation indicated that the cis active isomer of Asp282-Cys283 peptide bond was predominant around neutral pH while the trans bond gradually prevailed towards acidic environment. The protonation of His284 was found to be associated with the trans-to-cis isomerization of Asp282-Cys283 peptide bond which dominantly regulates the GlcNAc binding. Our MD simulation approach provides an insight into the pH-sensitive proteins and hence, ligand binding activity.

Nitrophorins: nitrite disproportionation reaction and other novel functionalities of insect heme-based nitric oxide transport proteins

Nitrophorins (NPs) comprise a unique class of heme proteins used by the blood-sucking insect Rhodnius prolixus to deliver the signaling gas molecule NO into the blood vessel of a host during feeding. Upon NO release, histamine can be scavenged by coordination to the heme iron. Although the protein is of similar size as the mammalian globin monomers and shares the same cofactor and proximal histidine coordination, nitrophorin structure, in contrast, is almost entirely composed of a β-barrel. Comparison of the NO and histamine association constants with the concentrations of both compounds invivo raises concerns about the very simple ligand release model in case of at least some of the NPs. Therefore, novel functionalities of the NPs were sought. As a result, catalysis of the nitrite disproportionation reaction was found, which leads to the formation of NO with nitrite as the sole substrate. This is the first example of a ferriheme protein that can perform this reaction. Furthermore, although NPs stabilize the ferriheme state, a peroxidase reactivity of the cofactor involving the higher oxidation state iron (Compound I/II) was studied with the potential to catalyze the oxidation of histamine and norepinephrine. In contrast to many other heme proteins including the globins, the ferroheme state was found to be extremely sensitive to O(2) , which is a consequence of the much lower reduction potential of the NPs, so that the 1-electron reduction of O(2) to O (•-)(2) becomes a thermodynamically favored process. Altogether, the detailed study of the NPs gives insight into the structure-function relationships required for the targeted delivery of diatomic gas molecules in biology. Moreover, the comparison of the structure-function relationships of the NPs (NO transporters) with those of the globins (O(2) transporters) will help to elucidate the architectural requirement for the respective tasks.

Accurate proton affinity and gas-phase basicity values for molecules important in biocatalysis

Benchmark quantum calculations of proton affinities and gas-phase basicities of molecules relevant to biochemical processes, particularly acid/base catalysis, are presented and compared for a variety of multilevel and density functional quantum models. Included are nucleic acid bases in both keto and enol tautomeric forms, ribose in B-form and A-form sugar pucker conformations, amino acid side chains and backbone molecules, and various phosphates and phosphoranes, including thio substitutions. This work presents a high-level thermodynamic characterization of biologically relevant protonation states and provides a benchmark database for development of next-generation semiempirical and approximate density functional quantum models and parametrization of methods to predict pK(a) values and relative solvation energies.

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.

Ultrafast dynamics of diatomic ligand binding to nitrophorin 4

Nitrophorin 4 (NP4) is a heme protein that stores and delivers nitric oxide (NO) through pH-sensitive conformational change. This protein uses the ferric state of a highly ruffled heme to bind NO tightly at low pH and release it at high pH. In this work, the rebinding kinetics of NO and CO to NP4 are investigated as a function of iron oxidation state and the acidity of the environment. The geminate recombination process of NO to ferrous NP4 at both pH 5 and pH 7 is dominated by a single approximately 7 ps kinetic phase that we attribute to the rebinding of NO directly from the distal pocket. The lack of pH dependence explains in part why NP4 cannot use the ferrous state to fulfill its function. The kinetic response of ferric NP4NO shows two distinct phases. The relative geminate amplitude of the slower phase increases dramatically as the pH is raised from 5 to 8. We assign the fast phase of NO rebinding to a conformation of the ferric protein with a closed hydrophobic pocket. The slow phase is assigned to the protein in an open conformation with a more hydrophilic heme pocket environment. Analysis of the ultrafast kinetics finds the equilibrium off-rate of NO to be proportional to the open state population as well as the pH-dependent amplitude of escape from the open pocket. When both factors are considered, the off-rate increases by more than an order of magnitude as the pH changes from 5 to 8. The recombination of CO to ferrous NP4 is observed to have a large nonexponential geminate amplitude with rebinding time scales of approximately 10(-11)-10(-9) s at pH 5 and approximately 10(-10)-10(-8) s at pH 7. The nonexponential CO rebinding kinetics at both pH 5 and pH 7 are accounted for using a simple model that has proven effective for understanding CO binding in a variety of other heme systems (Ye, X.; et al. Proc. Natl. Acad. Sci. U.S.A. 2007, 104, 14682).

1H and 13C NMR spectroscopic studies of the ferriheme resonances of three low-spin complexes of wild-type nitrophorin 2 and nitrophorin 2(V24E) as a function of pH

The ferriheme resonances of the low-spin (S = 1/2) complexes of wild-type (wt) nitrophorin 2 (NP2) and its heme pocket mutant NP2(V24E) with imidazole (ImH), histamine (Hm), and cyanide (CN⁻) as the sixth ligand have been investigated by NMR spectroscopy as a function of pH (4.0–7.5). For the three wt NP2 complexes, the ratio of the two possible heme orientational isomers, A and B, remains almost unchanged (ratio of A:B approximately 1:6 to 1:5) over this wide pH range. However, strong chemical exchange cross peaks appear in the nuclear Overhauser effect spectroscopy/exchange spectroscopy (NOESY/EXSY) spectra for the heme methyl resonances at low pH (pH* 4.0–5.5), which indicate chemical exchange between two species. We have shown these to be two different exogenous ImH or Hm orientations that are denoted B and B′, with the ImH plane nearly parallel and perpendicular to the ImH plane of the protein-provided His57, respectively. The wt NP2–CN complex also shows EXSY cross peaks due to chemical exchange, which is shown to be a result of interchange between two ruffling distortions of the heme. The same ruffling distortion interchange is also responsible for the ImH and Hm chemical exchange. For the three NP2(V24E) ligand complexes, no EXSY cross peaks are observed, but the A:B ratios change dramatically with pH. The fact that heme favors the A orientation highly for NP2(V24E) at low pH as compared with wt NP2 is believed to be due to the steric effect of the V24E mutation. The existence of the B′ species at lower pH for wt NP2 complexes and the increase in A heme orientation at lower pH for NP2(V24E) are believed to be a result of a change in structure near Glu53 when it is protonated at low pH. ¹H{¹³C} heteronuclear multiple quantum coherence (HMQC) spectra are very helpful for the assignment of heme and nearby protein side chain resonances.