Exploring the early stages of the pH‐induced conformational change of influenza hemagglutinin

Hemagglutinin (HA) mediates the membrane fusion process of influenza virus through its pH‐induced conformational change. However, it remains challenging to study its structure reorganization pathways in atomic details. Here, we first applied continuous constant pH molecular dynamics approach to predict the pK_a values of titratable residues in H2 subtype HA. The calculated net‐charges in HA1 globular heads increase from 0e (pH 7.5) to +14e (pH 4.5), indicating that the charge repulsion drives the detrimerization of HA globular domains. In HA2 stem regions, critical pH sensors, such as Glu103₂, His18₁, and Glu89₁, are identified to facilitate the essential structural reorganizations in the fusing pathways, including fusion peptide release and interhelical loop transition. To probe the contribution of identified pH sensors and unveil the early steps of pH‐induced conformational change, we carried out conventional molecular dynamics simulations in explicit water with determined protonation state for each titratable residue in different environmental pH conditions. Particularly, energy barriers involving previously uncharacterized hydrogen bonds and hydrophobic interactions are identified in the fusion peptide release pathway. Nevertheless, comprehensive comparisons across HA family members indicate that different HA subtypes might employ diverse pH sensor groups along with different fusion pathways. Finally, we explored the fusion inhibition mechanism of antibody CR6261 and small molecular inhibitor TBHQ, and discovered a novel druggable pocket in H2 and H5 subtypes. Our results provide the underlying mechanism for the pH‐driven conformational changes and also novel insight for anti‐flu drug development. Proteins 2014; 82:2412–2428. © 2014 Wiley Periodicals, Inc.     

Redox-coupled proton pumping in cytochrome c oxidase: Further insights from computer simulation

The membrane-bound enzyme cytochrome c oxidase, the terminal member in the respiratory chain, converts oxygen into water and generates an electrochemical gradient by coupling the electron transfer to proton pumping across the membrane. Here we have investigated the dynamics of an excess proton and the surrounding protein environment near the active sites. The multi-state empirical valence bond (MS-EVB) molecular dynamics method was used to simulate the explicit dynamics of proton transfer through the critically important hydrophobic channel between Glu242 (bovine notation) and the D-propionate of heme a₃ (PRDa₃) for the first time. The results from these molecular dynamics simulations indicate that the PRDa₃ can indeed re-orientate and dissociate from Arg438, despite the high stability of such an ion pair, and has the ability to accept protons via bound water molecules. Any large conformational change of the adjacent heme a D-propionate group is, however, sterically blocked directly by the protein. Free energy calculations of the PRDa₃ side chain isomerization and the proton translocation between Glu242 and the PRDa₃ site have also been performed. The results exhibit a redox state-dependent dynamical behavior and indicate that reduction of the low-spin heme a may initiate internal transfer of the pumped proton from Glu242 to the PRDa₃ site.     

Dynamical study,hydrogen bond analysis,and constant pH simulations of the beta carbonic anhydrase of Methanobacterium thermoautotrophicum

Within the five classes (α, β, γ, δ, and ζ) of carbonic anhydrases (CAs) the first two, containing mammal and plant representatives, are the most studied among all CAs. In this study, we have focused our investigation on the beta-class carbonic anhydrase of Methanobacterium thermoautotrophicum. We investigated both the importance of the Asp-Arg dyad near the catalytic zinc-bound water and the possible roles that water molecules within the active site and residues near the entrance of the catalytic cleft have on the first step of the enzyme’s reaction mechanism. Hydrogen-bonding analysis of selected residues within the active site and constant pH replica exchange molecular dynamics constant pH replica exchange simulations were performed. The latter was done in order to evaluate the pK_a values of possible proton acceptors. We found an intricate hydrogen-bonding network involving two acidic residues within the active site, Asp16 and Asp34, and the catalytic water molecule. We also observed a very strong interaction between the zinc-bound water and residues Asp34 and Arg36. This interaction was not significantly affected by the change in the protonation state of both the catalytic water and aspartate residue 34. The pK_a analysis show that the effect of the R36A mutation affects not only the possible proton acceptors, but also the catalytic water itself.     

Resolving the Negative Potential Side (n-side) Water-accessible Proton Pathway of F-type ATP Synthase by Molecular Dynamics Simulations

The rotation of F₁F_o-ATP synthase is powered by the proton motive force across the energy-transducing membrane. The protein complex functions like a turbine; the proton flow drives the rotation of the c-ring of the transmembrane F_o domain, which is coupled to the ATP-producing F₁ domain. The hairpin-structured c-protomers transport the protons by reversible protonation/deprotonation of a conserved Asp/Glu at the outer transmembrane helix (TMH). An open question is the proton transfer pathway through the membrane at atomic resolution. The protons are thought to be transferred via two half-channels to and from the conserved cAsp/Glu in the middle of the membrane. By molecular dynamics simulations of c-ring structures in a lipid bilayer, we mapped a water channel as one of the half-channels. We also analyzed the suppressor mutant cP24D/E61G in which the functional carboxylate is shifted to the inner TMH of the c-protomers. Current models concentrating on the “locked” and “open” conformations of the conserved carboxylate side chain are unable to explain the molecular function of this mutant. Our molecular dynamics simulations revealed an extended water channel with additional water molecules bridging the distance of the outer to the inner TMH. We suggest that the geometry of the water channel is an important feature for the molecular function of the membrane part of F₁F_o-ATP synthase. The inclination of the proton pathway isolates the two half-channels and may contribute to a favorable clockwise rotation in ATP synthesis mode.     

Prediction of Secondary Ionization of the Phosphate Group in Phosphotyrosine Peptides

A computational approach, based on a continuum molecular electrostatics model, for the calculation of the pK_a values of secondary ionization of the phosphate group in phenyl phosphate derivatives is described. The method uses the ESP atomic charges of the mono-anionic and di-anionic forms of the ionizable phosphate group, computed with the use of the density functional method, and applies a new concept of the model group, being the reference state for the pK_a calculations. Both conformational flexibility and tautomeric degrees of freedom are taken into account in the calculations. The method was parameterized using experimentally available pK_a values of four derivatives of phenyl phosphates, and phosphotyrosine. Subsequently this parameterization was used to predict pK_a of the phosphate group in a short peptide Gly-Gly-Tyr(P)-Ala, and in a longer peptide consisting of 12 residues, the latter in water, and in a complex with a protein—phospholipase. The agreement between the computed and the experimental pK_a values is better than ±0.3 pH units for the optimized solute dielectric constant of 11–13. This approach is promising and its extension to other phospho-amino acids is in progress.     

Dynamic water networks in cytochrome cbb3 oxidase

Heme-copper oxidases (HCOs) are terminal electron acceptors in aerobic respiration. They catalyze the reduction of molecular oxygen to water with concurrent pumping of protons across the mitochondrial and bacterial membranes. Protons required for oxygen reduction chemistry and pumping are transferred through proton uptake channels. Recently, the crystal structure of the first C-type member of the HCO superfamily was resolved [Buschmann et al. Science 329 (2010) 327–330], but crystallographic water molecules could not be identified. Here we have used molecular dynamics (MD) simulations, continuum electrostatic approaches, and quantum chemical cluster calculations to identify proton transfer pathways in cytochrome cbb₃. In MD simulations we observe formation of stable water chains that connect the highly conserved Glu323 residue on the proximal side of heme b₃ both with the N- and the P-sides of the membrane. We propose that such pathways could be utilized for redox-coupled proton pumping in the C-type oxidases. Electrostatics and quantum chemical calculations suggest an increased proton affinity of Glu323 upon reduction of high-spin heme b₃. Protonation of Glu323 provides a mechanism to tune the redox potential of heme b₃ with possible implications for proton pumping.     

Molecular dynamics study of the proposed proton transport pathways in [FeFe]-hydrogenase

Possible proton transport pathways in Clostridium pasteurianum (CpI) [FeFe]-hydrogenase were investigated with molecular dynamics simulations. This study was undertaken to evaluate the functional pathway and provide insight into the hydrogen bonding features defining an active proton transport pathway. Three pathways were evaluated, two of which consist of water wires and one of predominantly amino acid residues. Our simulations suggest that protons are not transported through water wires. Instead, the five-residue motif (Glu282, Ser319, Glu279, H₂O, Cys299) was found to be the likely pathway, consistent with previously made experimental observations. The pathway was found to have a persistent hydrogen bonded core (residues Cys299 to Ser319), with less persistent hydrogen bonds at the ends of the pathway for both H₂ release and H₂ uptake. Single site mutations of the four residues have been shown experimentally to deactivate the enzyme. The theoretical evaluation of these mutations demonstrates redistribution of the hydrogen bonds in the pathway, resulting in enzyme deactivation. Finally, coupling between the protein dynamics near the proton transport pathway and the redox partner binding regions was also found as a function of H₂ uptake and H₂ release states, which may be indicative of a correlation between proton and electron movement within the enzyme.     

Water-gated mechanism of proton translocation by cytochrome c oxidase

Cytochrome c oxidase is essential for aerobic life as a membrane-bound energy transducer. O₂ reduction at the haem a₃-Cu_B centre consumes electrons transferred via haem a from cytochrome c outside the membrane. Protons are taken up from the inside, both to form water and to be pumped across the membrane (M.K.F. Wikström, Nature 266 (1977) 271 [1]; M. Wikström, K. Krab, M. Saraste, Cytochrome Oxidase, A Synthesis, Academic Press, London, 1981 [2]). The resulting electrochemical proton gradient drives ATP synthesis (P. Mitchell, Chemiosmotic Coupling in Oxidative and Photosynthetic Phosphorylation, Glynn Research, Bodmin, UK, 1966 [3]). Here we present a molecular mechanism for proton pumping coupled to oxygen reduction that is based on the unique properties of water in hydrophobic cavities. An array of water molecules conducts protons from a conserved glutamic acid, either to the Δ-propionate of haem a₃ (pumping), or to haem a₃-Cu_B (water formation). Switching between these pathways is controlled by the redox-state-dependent electric field between haem a and haem a₃-Cu_B, which determines the water-dipole orientation, and therefore the proton transfer direction. Proton transfer via the propionate provides a gate to O₂ reduction. This pumping mechanism explains the unique arrangement of the metal cofactors in the structure. It is consistent with the large body of biochemical data, and is shown to be plausible by molecular dynamics simulations.     

Breaking the Carboxyl Rule: LYSINE 96 FACILITATES REPROTONATION OF THE SCHIFF BASE IN THE PHOTOCYCLE OF A RETINAL PROTEIN FROM EXIGUOBACTERIUM SIBIRICUM*?

A lysine instead of the usual carboxyl group is in place of the internal proton donor to the retinal Schiff base in the light-driven proton pump of Exiguobacterium sibiricum (ESR). The involvement of this lysine in proton transfer is indicated by the finding that its substitution with alanine or other residues slows reprotonation of the Schiff base (decay of the M intermediate) by more than 2 orders of magnitude. In these mutants, the rate constant of the M decay linearly decreases with a decrease in proton concentration, as expected if reprotonation is limited by the uptake of a proton from the bulk. In wild type ESR, M decay is biphasic, and the rate constants are nearly pH-independent between pH 6 and 9. Proton uptake occurs after M formation but before M decay, which is especially evident in D₂O and at high pH. Proton uptake is biphasic; the amplitude of the fast phase decreases with a pK_a of 8.5 ± 0.3, which reflects the pK_a of the donor during proton uptake. Similarly, the fraction of the faster component of M decay decreases and the slower one increases, with a pK_a of 8.1 ± 0.2. The data therefore suggest that the reprotonation of the Schiff base in ESR is preceded by transient protonation of an initially unprotonated donor, which is probably the ϵ-amino group of Lys-96 or a water molecule in its vicinity, and it facilitates proton delivery from the bulk to the reaction center of the protein.     

Network analysis of a proposed exit pathway for protons to the P-side of cytochrome c oxidase

Cytochrome c Oxidase (CcO) reduces O₂, the terminal electron acceptor, to water in the aerobic, respiratory electron transport chain. The energy released by O₂ reductions is stored by removing eight protons from the high pH, N-side, of the membrane with four used for chemistry in the active site and four pumped to the low pH, P-side. The proton transfers must occur along controllable proton pathways that prevent energy dissipating movement towards the N-side. The CcO N-side has well established D- and K-channels to deliver protons to the protein interior. The P-side has a buried core of hydrogen-bonded protonatable residues designated the Proton Loading Site cluster (PLS cluster) and many protonatable residues on the P-side surface, providing no obvious unique exit. Hydrogen bond pathways were identified in Molecular Dynamics (MD) trajectories of Rb. sphaeroides CcO prepared in the P_R state with the heme a₃ propionate and Glu286 in different protonation states. Grand Canonical Monte Carlo sampling of water locations, polar proton positions and residue protonation states in trajectory snapshots identify a limited number of water mediated, proton paths from PLS cluster to the surface via a (P-exit) cluster of residues. Key P-exit residues include His93, Ser168, Thr100 and Asn96. The hydrogen bonds between PLS cluster and P-exit clusters are mediated by a water wire in a cavity centered near Thr100, whose hydration can be interrupted by a hydrophobic pair, Leu255B (near Cu_A) and Ile99. Connections between the D channel and PLS via Glu286 are controlled by a second, variably hydrated cavity.

Constant pH molecular dynamics (CpHMD) and molecular docking studies of CquiOBP1 pH-induced ligand releasing mechanism

The odorant binding protein of Culex quinquefasciatus (CquiOBP1), expressed on the insect antenna, is crucial for the investigation of trapping baited with oviposition semi-chemicals and controlling mosquito populations. The acidic titratable residues pKa prediction and the ligand binding poses investigation in two systems (pH 7 and pH 5) are studied by constant pH molecular dynamics (CpHMD) and molecular docking methods. Research results reveal that the change of the protonation states would disrupt some important H-bonds, such as Asp 66-Asp 70, Glu 105-Asn 102, etc. The cleavage of these H-bonds leads to the movement of the relative position of hydrophobic tunnel, N- and C- termini loops and pH-sensing triad (His23-Tyr54-Val125) in acid solution. Ligand MOP has lower affinity and shows different binding poses to protein CquiOBP1 at pH 5. This ligand may be released from another tunnel between helices α3 and α4 in acidic environment. However, it would bind to the protein with high affinity in neutral environment. This work could provide more penetrating understanding of the pH-induced ligand-releasing mechanism.     

Dynamics of the glutamic acid 242 side chain in cytochrome c oxidase

In many cytochrome c oxidases glutamic acid 242 is required for proton transfer to the binuclear heme a₃/Cu_B site, and for proton pumping. When present, the side chain of Glu-242 is orientated “down” towards the proton-transferring D-pathway in all available crystal structures. A nonpolar cavity “above” Glu-242 is empty in these structures. Yet, proton transfer from Glu-242 to the binuclear site, and for proton-pumping, is well established, and the cavity has been proposed to at least transiently contain water molecules that would mediate proton transfer. Such proton transfer has been proposed to require isomerisation of the Glu-242 side chain into an “up” position pointing towards the cavity. Here, we have explored the molecular dynamics of the protonated Glu-242 side chain. We find that the “up” position is preferred energetically when the cavity contains four water molecules, but the “down” position is favoured with less water. We conclude that the cavity might be deficient in water in the crystal structures, possibly reflecting the “resting” state of the enzyme, and that the “up/down” equilibrium of Glu-242 may be coupled to the presence of active-site water molecules produced by O₂ reduction.     

Complexation of aluminium with (−) -epigallocathechin gallate studied by spectrophotometry

Complexation equilibria between Al(III) and (−)-epigallocathechin gallate (EGCG) in the presence of acetate buffer have been studied by spectrophotometry. The method is based on the competition between EGCG and buffer ligands for Al(III) ions. The apparent formation constant of the EGCG complex for Al(III), which could be determined by measuring the absorbance of the free EGCG, decreased with increasing acetate ion concentration at a fixed pH. This phenomenon has been quantitatively investigated and both types of complexes (EGCG and acetate) could be analyzed. The apparent formation constant of Al(III) complex with EGCG also decreased with decreasing pH at a fixed acetate ion concentration. The pH dependence of the apparent formation constant indicates the 1:1 competition between metal ions and hydrogen ions for the binding site of EGCG. The intrinsic formation constant of Al-EGCG complex, the proton association constant of EGCG and the formation constant of Al-acetate complex are found to be log K_Al-EGCG = 7.6 ± 0.1, log K_H-EGCG = 7.65 ± 0.03 and log K_Al-acetate = 2.07 ± 0.05 by a graphical analysis.     

Effect of dielectric constant of medium on protonation equilibria of ethylenediamine

Abstract

The effect of dielectric constant of medium on protonation equilibria has been studied by determining protonation constants of ethylenediamine pH metrically in various concentrations (0–60%v/v) of acetoni-trile– and ethylene glycol–water mixtures, at an ionic strength of 0.16mol L^?1 and at 303.0 K. MINIQUAD75 computer program has been used for the calculation of protonation constants. Linear and non-linear variations of step-wise protonation constants with reciprocal of dielectric constant of the solvent mixtures have been attributed to the dominance of the electrostatic and non-electrostatic forces, respectively. The trend is explained on the basis of solute–solute and solute–solvent interactions, solvation, proton transfer processes and dielectric constants of the media.     

Morphological and transport characteristics of swollen chitosan‐based proton exchange membranes studied by molecular modeling

Chitosan biopolymer has been extensively applied in direct methanol fuel cells (DMFCs) as a potential replacement to conventional Nafion membrane for its considerably reduced methanol crossover. Here, we computationally explored the influences of methanol concentration, temperature, and pH parameters upon the nanostructure and dynamics, particularly the methanol crossover, in chitosan proton‐exchange membrane (PEM) through molecular dynamics simulations. Theoretical results demonstrated the increased swelling and radius of gyration of chitosan chains at higher concentrations. Structural examinations further revealed that an increase in methanol loading weakened the water interactions with chitosan functionalities (amine? NH₂, hydroxyl? OH, and methoxy? CH₂OH) whereas improved the methanol affinities toward chitosan, reflecting higher methanol sorption capability of chitosan at enhanced concentrations. Additionally, it was found that interactions between solvents and chitosan strengthened under acidic pH conditions on account of amine protonation. The water diffusivity inside the swollen chitosan diminished by increasing CH₃OH uptake, and in contrast diffusivity of methanol was noted to enhance. Furthermore, it was observed that an enhancement in temperature or a decrease in pH intensified solvent mobility. These insights imply that supplying methanol‐concentrated and/or acidic feed solutions into DMFCs based on chitosan PEMs could lower membrane performance due to the significant methanol transport dynamics.     

The role of Arginine 38 in horseradish peroxidase enzyme revisited: A computational investigation

Molecular dynamics simulations on hydrogen peroxide complex with wild-type (WT) and Arg38Leu mutated (R38L) Horseradish Peroxidase (HRP) were carried out over nanoseconds timescale in water solution at 300 K. Comparison of the results provides interesting insights about the role of highly conserved Arg38 and His42 residues in the chemical features of HRP, underlying its biological activity which initiates with Compound0 (Cpd0). In the WT-HRP enzyme current molecular dynamics simulations show, for the first time, that Arg38 residue: i) prevents the entrance of water inside the reaction cavity, hence providing a hydrophobic reactive scenario, ii) it maintains the distance between His42 and heme–H₂O₂ complex suitable for the occurrence of proton transfer reaction leading, thereafter, to heme–H₂O₂ disruption according to Poulos-Kraut mechanism. On the other hand, R38L mutant can be considered as a “wet enzyme” where the presence of water solvent molecules in the heme reaction pocket, unfavoring the initial heme–H₂O₂ complex formation, decreases the catalytic efficiency in agreement with experimental kinetics measurements. Furthermore, we note that Arg38Leu mutation pushes the His42 residue far from the heme–H₂O₂ complex, making unlikely a direct proton transfer and suggesting that, in the mutant, a solvent water molecule could be involved in the first step of the Poulos-Kraut mechanism.     

pH replica-exchange method based on discrete protonation states

We propose a new algorithm for obtaining proton titration curves of ionizable residues. The algorithm is a pH replica-exchange method (PHREM), which is based on the constant pH algorithm of Mongan et al. (J Comput Chem 2004;25:2038-2048). In the original replica-exchange method, simulations of different replicas are performed at different temperatures, and the temperatures are exchanged between the replicas. In our PHREM, simulations of different replicas are performed at different pH values, and the pHs are exchanged between the replicas. The PHREM was applied to a blocked amino acid and to two protein systems (snake cardiotoxin and turkey ovomucoid third domain), in conjunction with a generalized Born implicit solvent. The performance and accuracy of this algorithm and the original constant pH method (PHMD) were compared. For a single set of simulations at different pHs, the use of PHREM yields more accurate Hill coefficients of titratable residues. By performing multiple sets of constant pH simulations started with different initial states, the accuracy of predicted pK(a) values and Hill coefficients obtained with PHREM and PHMD methods becomes comparable. However, the PHREM algorithm exhibits better samplings of the protonation states of titratable residues and less scatter of the titration points and thus better precision of measured pK(a) values and Hill coefficients. In addition, PHREM exhibits faster convergence of individual simulations than the original constant pH algorithm.     

Transient Binding of CO to CuB in Cytochrome c Oxidase Is Dynamically Linked to Structural Changes around a Carboxyl Group: A Time-Resolved Step-Scan Fourier Transform Infrared Investigation

The redox-driven proton pump cytochrome c oxidase is that enzymatic machinery of the respiratory chain that transfers electrons from cytochrome c to molecular oxygen and thereby splits molecular oxygen to form water. To investigate the reaction mechanism of cytochrome c oxidase on the single vibrational level, we used time-resolved step-scan Fourier transform infrared spectroscopy and studied the dynamics of the reduced enzyme after photodissociation of bound carbon monoxide across the midinfrared range (2300-950 cm⁻¹). Difference spectra of the bovine complex were obtained at -20°C with 5 μs time resolution. The data demonstrate a dynamic link between the transient binding of CO to Cu_B and changes in hydrogen bonding at the functionally important residue E(I-286). Variation of the pH revealed that the pK_a of E(I-286) is >9.3 in the fully reduced CO-bound oxidase. Difference spectra of cytochrome c oxidase from beef heart are compared with those of the oxidase isolated from Rhodobacter sphaeroides. The bacterial enzyme does not show the environmental change in the vicinity of E(I-286) upon CO dissociation. The characteristic band shape appears, however, in redox-induced difference spectra of the bacterial enzyme but is absent in redox-induced difference spectra of mammalian enzyme. In conclusion, it is demonstrated that the dynamics of a large protein complex such as cytochrome c oxidase can be resolved on the single vibrational level with microsecond Fourier transform infrared spectroscopy. The applied methodology provides the basis for future investigations of the physiological reaction steps of this important enzyme.     

High resolution nuclear magnetic resonance study of the histidine--aspartate hydrogen bond in chymotrypsin and chymotrypsinogen

A high resolution proton nuclear magnetic resonance study of chymotrypsin A_δ and Chymotrypsinogen A in water has shown a single resonance at very low magnetic fields (− 18 to − 15 p.p.m. relative to dimethyl-silapentane-sulfonate). From its pH dependence (pK = 7·2) and response to chemical modification the resonance has been assigned to the hydrogen-bonded proton between His-57 and Asp-102.     

Constant pH molecular dynamics of proteins in explicit solvent with proton tautomerism

pH is a ubiquitous regulator of biological activity, including protein‐folding, protein‐protein interactions, and enzymatic activity. Existing constant pH molecular dynamics (CPHMD) models that were developed to address questions related to the pH‐dependent properties of proteins are largely based on implicit solvent models. However, implicit solvent models are known to underestimate the desolvation energy of buried charged residues, increasing the error associated with predictions that involve internal ionizable residue that are important in processes like hydrogen transport and electron transfer. Furthermore, discrete water and ions cannot be modeled in implicit solvent, which are important in systems like membrane proteins and ion channels. We report on an explicit solvent constant pH molecular dynamics framework based on multi‐site λ‐dynamics (CPHMD^MSλD). In the CPHMD^MSλD framework, we performed seamless alchemical transitions between protonation and tautomeric states using multi‐site λ‐dynamics, and designed novel biasing potentials to ensure that the physical end‐states are predominantly sampled. We show that explicit solvent CPHMD^MSλD simulations model realistic pH‐dependent properties of proteins such as the Hen‐Egg White Lysozyme (HEWL), binding domain of 2‐oxoglutarate dehydrogenase (BBL) and N‐terminal domain of ribosomal protein L9 (NTL9), and the pK_a predictions are in excellent agreement with experimental values, with a RMSE ranging from 0.72 to 0.84 pK_a units. With the recent development of the explicit solvent CPHMD^MSλD framework for nucleic acids, accurate modeling of pH‐dependent properties of both major class of biomolecules—proteins and nucleic acids is now possible. Proteins 2014; 82:1319–1331. © 2013 Wiley Periodicals, Inc.