Hydrogen bonding in nucleosides and nucleotides

An analysis of the hydrogen bonding in 76 nucleoside and 11 nucleotide crystal structures shows that the hydrogen bond lengths fall into well-defined categories according to the nature of the donor or acceptor groups. The shortest bonds are those involving P---OH or O=P groups. For donor groups, the sequence in bond lengths is

P—OH<C—OH< N−H<O_w(H)—H<N(H)—H<C−H

There are ten examples of two centre

H—H…O

bonds, which are comparable in length with P---OH …O bonds. The acceptor seqeunce is

O=P<OH₂<OH₂<O=CO(H)C<N N(H₂)C<Cl^…<O<S=C

The number of three-centre bonds, about 24%, is comparable to that observed in the carbohydrates and the amino acids. Most hydrogen bonds are involved in short finite chains. Only in the nucleotides are cyclic hydrogen bonding schemes observed.     

Proton transfer along the external proton channel of bacteriorhodopsin

The process of proton transfer along a proton channel is considered using bacteriorhodopsin as a model system, for which a large body of experimental data is available. The possible amino acid composition of the external proton half-channel of bacteriorhodopsin and the stepwise scheme of proton transfer consistent with experimental data are proposed. The rate of proton transfer between fixed centers is assessed for certain regions of this channel for which spectroscopic data are available.     

The mechanism of proton transfer between adjacent sites on the molecular surface

The surface of a protein, or a membrane, is spotted with a multitude of proton binding sites, some of which are only few Å apart. When a proton is released from one site, it propagates through the water by a random walk under the bias of the local electrostatic potential determined by the distribution of the charges on the protein. Eventually, the released protons are dispersed in the bulk, but during the first few nanoseconds after the dissociation, the protons can be trapped by encounter with nearby acceptor sites. While the study of this reaction on the surface of a protein suffers from experimental and theoretical difficulties, it can be investigated with simple model compounds like derivatives of fluorescein. In the present study, we evaluate the mechanism of proton transfer reactions that proceed, preferentially, inside the Coulomb cage of the dye molecules. Kinetic analysis of the measured dynamics reveals the role of the dimension of the Coulomb cage on the efficiency of the reaction and how the ordering of the water molecules by the dye affects the kinetic isotope effect.     

Unraveling the mechanism of proton translocation in the extracellular half‐channel of bacteriorhodopsin

Bacteriorhodopsin, a light activated protein that creates a proton gradient in halobacteria, has long served as a simple model of proton pumps. Within bacteriorhodopsin, several key sites undergo protonation changes during the photocycle, moving protons from the higher pH cytoplasm to the lower pH extracellular side. The mechanism underlying the long‐range proton translocation between the central (the retinal Schiff base SB216, D85, and D212) and exit clusters (E194 and E204) remains elusive. To obtain a dynamic view of the key factors controlling proton translocation, a systematic study using molecular dynamics simulation was performed for eight bacteriorhodopsin models varying in retinal isomer and protonation states of the SB216, D85, D212, and E204. The side‐chain orientation of R82 is determined primarily by the protonation states of the residues in the EC. The side‐chain reorientation of R82 modulates the hydrogen‐bond network and consequently possible pathways of proton transfer. Quantum mechanical intrinsic reaction coordinate calculations of proton‐transfer in the methyl guanidinium‐hydronium‐hydroxide model system show that proton transfer via a guanidinium group requires an initial geometry permitting proton donation and acceptance by the same amine. In all the bacteriorhodopsin models, R82 can form proton wires with both the CC and the EC connected by the same amine. Alternatively, rare proton wires for proton transfer from the CC to the EC without involving R82 were found in an O′ state where the proton on D85 is transferred to D212. Proteins 2016; 84:639–654. © 2016 Wiley Periodicals, Inc.     

Identification of proton transfer pathways in the X-ray crystal structure of the bacterial reaction center from Rhodobacter sphaeroides

Structural features that have important implications for the fundamental process of transmembrane proton transfer are examined in the recently published high resolution atomic structures of the reaction center (RC) from Rhodobacter sphaeroides in the dark adapted state (DQ_AQ_B) and the charged separated state (D⁺Q_AQ_B ^–); the latter is the active state for proton transfer to the semiquinone. The structures have been determined at 2.2 Å and 2.6 Å resolution, respectively, as reported by Stowell et al. (1997) [Science 276: 812–816]. Three possible proton transfer pathways (P1, P2, P3) consisting of water molecules and/or protonatable residues were identified which connect the Q_B binding region with the cytoplasmic exposed surface at Asp H224 & Asp M240 (P1), Tyr M3 (P2) and Asp M17 (P3). All three represent possible pathways for proton transfer into the RC. P1 contains an uninterrupted chain of water molecules. This path could, in addition, facilitate the exchange of quinone for quinol during the photocycle by allowing water to move into and out of the binding pocket. Located near these pathways is a cluster of electrostatically interacting acid residues (Asp-L213, Glu-H173, Asp-M17, Asp H124, Asp-L210 and Asp H170) each being within 4.5 Å of a neighboring carboxylic acid or a bridging water molecule. This cluster could serve as an internal proton reservoir facilitating fast protonation of Q_B ^– that could occur at a rate greater than that attainable by proton uptake from solution.     

Theory of passive proton conductance in lipid bilayers

The large permeability of lipid bilayers to protons compared to other small ions calls for a special proton transport mechanism. At the present time, only mechanisms involving transient hydrogen-bonded chains of water can account for the experimental result that the conductance is nearly independent of pH. Three models involving transient hydrogen-bonded chains are discussed, including an outline of the kinetic calculations that lead to predictions of current versus voltage drop and current versus pH differences. These calculations can be compared to experiment to determine which, if any, of these models pertains to lipid bilayers.     

Resolution of electron and proton transfer events in the electrochromism associated with quinone reduction in bacterial reaction centers

We have measured the electrochromic response of the bacteriopheophytin, BPh, and bacteriochlorophyll, BChl, cofactors during the Q_A ^–Q_B Q_AQ_B ^– electron transfer in chromatophores of Rhodobacter (Rb.) capsulatus and Rb. sphaeroides. The electrochromic response rises faster in chromatophores and is more clearly biexponential than it is in isolated reaction centers. The chromatophore spectra can be interpreted in terms of a clear kinetic separation between fast electron transfer and slower non-electron transfer events such as proton transfer or protein relaxation. The electrochromic response to electron transfer exhibits rise times of about 4 µs (70%) and 40 µs (30%) in Rb. capsulatus and 4 µs (60%) and 80 µs (40%) in Rb. sphaeroides. The BPh absorption band is shifted to nearly equivalent positions in the Q_A ^– and nascent Q_B ^– states, indicating that the electrochromic perturbation of BPh absorption from the newly formed Q_B ^– state is comparable to that of Q_A ^– . Subsequently, partial attenuation of the Q_B ^– electrochromism occurs with a time constant on the order of 200 µs. This can be attributed to partial charge compensation by H⁺ (or other counter ion) movement into the Q_B pocket. Electron transfer events were found to be slower in detergent isolated RCs than in chromatophores, more nearly monoexponential, and overlap H⁺ transfer, suggesting that a change in rate-limiting step has occurred upon detergent solubilization.     

Kinetics and thermodynamics of proton transfer to CpRu(dppe)H: Via dihydrogen bonding and (η-H2)-complex to the dihydride

The interaction between Cp^∗RuH(dppe) and a series of proton donors (HA) of increasing strength: CFH₂CH₂OH (MFE), CF₃CH₂OH (TFE), (CF₃)₂CHOH (HFIP), p-nitrophenol, CF₃COOH and HBF₄ has been investigated spectroscopically by variable-temperature IR, UV-Vis, and NMR spectroscopy in solvents of differing polarity (n-hexane, dichloromethane and their mixture). The low-temperature IR study shows the establishment of a hydrogen-bond which involves the hydride ligand as the proton accepting site. The basicity factor E_j for the hydride was found to be 1.39. All techniques indicate that an equilibrium exists between the dihydrogen-bonded complex and the cationic dihydrogen complex, [Cp^∗Ru(η²-H₂)(dppe)]⁺, the formation of which is shown here for the first time. The proton transfer from HFIP is characterized by ΔH^° = −8.1 ± 0.6 kcal mol⁻¹ and ΔS^° = −17 ± 3 eu. The activation parameters for the subsequent irreversible isomerization leading to the classical dihydride complex, [Cp^∗Ru(H)₂(dppe)]⁺, are ΔH^‡ = 20.9 ± 0.8 kcal mol⁻¹ and ΔS^‡ = 9 ± 3 eu as determined from ¹H NMR spectroscopy for protonation by HBF₄. Computational results at the DFT/B3PW91 level confirm the experimentally observed hydride basicity increase on descending the Group from iron to ruthenium and also the formation of the non-classical complex as an intermediate, prior to the thermodynamically favored dihydride.     

Contribution of proton linkage to the thermodynamic stability of the major cold-shock protein of Escherichia coli CspA

The stability of protein is defined not only by the hydrogen bonding, hydrophobic effect, van der Waals interactions, and salt bridges. Additional, much more subtle contributions to protein stability can arise from surface residues that change their properties upon unfolding. The recombinant major cold shock protein of Escherichia coli CspA an all-beta protein unfolds reversible in a two-state manner, and behaves in all other respects as typical globular protein. However, the enthalpy of CspA unfolding strongly depends on the pH and buffer composition. Detailed analysis of the unfolding enthalpies as a function of pH and buffers with different heats of ionization shows that CspA unfolding in the pH range 5.5-9.0 is linked to protonation of an amino group. This amino group appears to be the N-terminal alpha-amino group of the CspA molecule. It undergoes a 1.6 U shift in pKa values between native and unfolded states. Although this shift in pKa is expected to contribute approximately 5 kJ/mol to CspA stabilization energy the experimentally observed stabilization is only approximately 1 kJ/mol. This discrepancy is related to a strong enthalpy-entropy compensation due, most likely, to the differences in hydration of the protonated and deprotonated forms of the alpha-amino group.     

Molecular dynamics simulation to investigate anhydrous phosphoric acid-doped polybenzimidazole

The study conducted a molecular dynamics simulation based on a condensed-phase optimised molecular potentials for atomistic simulation studies force field model to investigate an anhydrous system of phosphoric acid-doped polybenzimidazole (poly[2,2′-(m-phenylene)-5,5′-bibenzimidazole], PBI). Intermolecular pair correlation functions and corresponding coordination numbers were calculated to research the strengths for various types of hydrogen bonding. The results display that the strengths of the hydrogen bonding interactions are in the order of o1–h pair > o2–h pair > n2a–h pair > n3a–h pair, and most protons are located around the neighbourhood of H₂PO₄^– rather than that of PBI. The proton conductivities are 3.86 × 10^?3 S cm^?1 at 298 K and 8.50 × 10^?3 S cm^?1 at 413 K. Moreover, the value obtained from our simulation system at 413 K is within the same order of magnitude as the experimentally measured value 0.012 S cm^?1 at 420% doping level. The distribution of proton displacement exhibits that the displacement of most protons is about 1.25–2.5 Å. The displacement is over 3.0 Å only for a fraction of protons. In addition, the greatest displacement can approach 4.595 Å. The trajectory analyses of protons show that the most possible mechanisms of proton transfer come from three ways: (a) between two H₂PO₄^– anions, (b) between H₂PO₄^– anions and benzimidazole moieties and (c) between two benzimidazole moieties. The dynamics of polymer motion was studied by the trajectory analyses of ring flips. The large amplitude flips of rings in the polymer chains were found in the system. The flips between benzene and benzimidazole are more frequent than that between benzimidazole moieties.     

Flash-induced proton transfer in photosynthetic bacteria

A proton electrochemical potential across the membranes of photosynthetic purple bacteria is established by a light-driven proton pump mechanism: the absorbed light in the reaction center initiates electron transfer which is coupled to the vectorial displacement of protons from the cytoplasm to the periplasm. The stoichiometry and kinetics of proton binding and release can be tracked directly by electric (glass electrodes), spectrophotometric (pH indicator dyes) and conductimetric techniques. The primary step in the formation of the transmembrane chemiosmotic potential is the uptake of two protons by the doubly reduced secondary quinone in the reaction center and the subsequent exchange of hydroquinol for quinone from the membrane quinone-pool. However, the proton binding associated with singly reduced promary and/or secondary quinones of the reaction center is substoichiometric, pH-dependent and its rate is electrostatically enhanced but not diffusion limited. Molecular details of protonation are discussed based on the crystallographic structure of the reaction center of purple bacteriaRb. sphaeroides andRps. viridis, structure-based molecular (electrostatic) calculations and mutagenesis directed at protonatable amino acids supposed to be involved in proton conduction pathways.     

Synthesis and characterization of a proton transfer salt between 2,6-pyridinedicarboxylic acid and 2-aminobenzothiazole,and its complexes and their inhibition studies on carbonic anhydrase isoenzymes

Abstract

A novel proton transfer compound (HABT)⁺(Hdipic)^? (1) obtained from ABT and H₂dipic and its metal complexes (2–5) have been prepared and characterized by spectroscopic techniques. Single crystal X-ray diffraction method has also been applied to 2 and 5. While complex 2 has a distorted octahedral conformation, 5 exhibits a distorted square pyramidal structure. The structures of 3 and 4 might be proposed as octahedral according to experimental data. All compounds were also evaluated for their in vitro inhibition effects on hCA I and II for their hydratase and esterase activities. Although there is no inhibition for hydratase activities, all compounds have inhibited the esterase activities of hCA I and II. The comparison of the inhibition studies of 1–5 to parent compounds indicates that 1–5 have superior inhibitory effects. The inhibition effects of 2–5 are also compared to inhibitory properties of the metal complexes of ABT and H₂dipic, revealing an improved transfection profile.     

Ground‐ and excited‐state proton transfer and antioxidant activity of 3‐hydroxyflavone in egg yolk phosphatidylcholine liposomes: absorption and fluorescence spectroscopic studies

Excited‐state intramolecular proton transfer (ESIPT) and dual luminescence behaviour of 3‐hydroxyflavone (3‐HF) have been utilized to monitor its binding to liposomal membranes prepared from egg yolk phosphatydilcholine (EYPC). Additionally, absorption spectrophotometric assay has been performed to evaluate the antioxidant activity of 3‐HF against lipid peroxidation in this membrane system. When 3‐HF molecules are partitioned into EYPC liposomes, a weak long‐wavelength absorption band with λ ～410 nm appears in addition to the principal absorption at ～λ = 345 nm. Selective excitation of the 410 nm band produces the characteristic emission (λ～460 nm) of the ground‐state anionic species, whereas excitation at the higher energy absorption band leads to dual emission with predominatly ESIPT tautomer fluorescence (λ = 528 nm). Both ESIPT tautomer and the anionic species exhibit fairly high fluorescence anisotropy (r) values (r = 0.122 and 0.180, respectively). Biexponential fluorescence decay kinetics are observed for the ESIPT tautomer as well as the ground‐state anionic forms, indicating heterogeneity in the microenvironments of the corresponding emitting species. Furthermore, we demonstrate that lipid peroxidation of EYPC liposomes is significantly inhibited upon 3‐HF binding, suggesting that 3‐HF can be potentially useful as an inhibitor of peroxidative damage of cell membranes. Copyright © 2008 John Wiley & Sons, Ltd.     

Potentials of mean force in acidic proton transfer reactions in constrained geometries

Free energy barriers associated with the transfer of an excess proton in water and related to the potentials of mean force in proton transfer episodes have been computed in a wide range of thermodynamic states, from low-density amorphous ices to high-temperature liquids under the critical point for unconstrained and constrained systems. The latter were represented by set-ups placed inside hydrophobic graphene slabs at the nanometric scale allocating a few water layers, namely one or two in the narrowest case. Water–proton and carbon–proton forces were modelled with a Multi-State Empirical Valence Bond method. As a general trend, a competition between the effects of confinement and temperature is observed on the local hydrogen-bonded structures around the lone proton and, consequently, on the mean force exerted by its environment on the water molecule carrying the proton. Free energy barriers estimated from the computed potentials of mean force tend to rise with the combined effect of increasing temperatures and the packing effect due to a larger extent of hydrophobic confinement. The main reason observed for such enhancement of the free energy barriers was the breaking of the second coordination shell around the lone proton.     

The role of hydrogen bonding in the enzymatic reaction catalyzed by HIV-1 protease

The hydrogen-bond network in various stages of the enzymatic reaction catalyzed by HIV-1 protease was studied through quantum-classical molecular dynamics simulations. The approximate valence bond method was applied to the active site atoms participating directly in the rearrangement of chemical bonds. The rest of the protein with explicit solvent was treated with a classical molecular mechanics model. Two possible mechanisms were studied, general-acid/general-base (GA/GB) with Asp 25 protonated at the inner oxygen, and a direct nucleophilic attack by Asp 25. Strong hydrogen bonds leading to spontaneous proton transfers were observed in both reaction paths. A single-well hydrogen bond was formed between the peptide nitrogen and outer oxygen of Asp 125. The proton was diffusely distributed with an average central position and transferred back and forth on a picosecond scale. In both mechanisms, this interaction helped change the peptide-bond hybridization, increased the partial charge on peptidyl carbon, and in the GA/GB mechanism, helped deprotonate the water molecule. The inner oxygens of the aspartic dyad formed a low-barrier, but asymmetric hydrogen bond; the proton was not positioned midway and made a slightly elongated covalent bond, transferring from one to the other aspartate. In the GA/GB mechanism both aspartates may help deprotonate the water molecule. We observed the breakage of the peptide bond and found that the protonation of the peptidyl amine group was essential for the peptide-bond cleavage. In studies of the direct nucleophilic mechanism, the peptide carbon of the substrate and oxygen of Asp 25 approached as close as 2.3 A.     

Hydrogen bonding is the prime determinant of carboxyl pKa values at the N-termini of alpha-helices

Experimentally determined mean pK(a) values of carboxyl residues located at the N-termini of alpha-helices are lower than their overall mean values. Here, we perform three types of analyses to account for this phenomenon. We estimate the magnitude of the helix macrodipole to determine its potential role in lowering carboxyl pK(a) values at the N-termini. No correlation between the magnitude of the macrodipole and the pK(a) values is observed. Using the pK(a) program propKa we compare the molecular surroundings of 18 N-termini carboxyl residues versus 233 protein carboxyl groups from a previously studied database. Although pK(a) lowering interactions at the N-termini are similar in nature to those encountered in other protein regions, pK(a) lowering backbone and side-chain hydrogen bonds appear in greater number at the N-termini. For both Asp and Glu, there are about 0.5 more hydrogen bonds per residue at the N-termini than in other protein regions, which can be used to explain their lower than average pK(a) values. Using a QM-based pK(a) prediction model, we investigate the chemical environment of the two lowest Asp and the two lowest Glu pK(a) values at the N-termini so as to quantify the effect of various pK(a) determinants. We show that local interactions suffice to account for the acidity of carboxyl residues at the N-termini. The effect of the helix dipole on carboxyl pK(a) values, if any, is marginal. Backbone amide hydrogen bonds constitute the single biggest contributor to the lowest carboxyl pK(a) values at the N-termini. Their estimated pK(a) lowering effects range from about 1.0 to 1.9 pK(a) units.     

Conotoxins as sensors of local pH and electrostatic potential in the outer vestibule of the sodium channel

We examined the block of voltage-dependent rat skeletal muscle sodium channels by derivatives of mu-conotoxin GIIIA (muCTX) having either histidine, glutamate, or alanine residues substituted for arginine-13. Toxin binding and dissociation were observed as current fluctuations from single, batrachotoxin-treated sodium channels in planar lipid bilayers. R13X derivatives of muCTX only partially block the single-channel current, enabling us to directly monitor properties of both muCTX-bound and -unbound states under different conditions. The fractional residual current through the bound channel changes with pH according to a single-site titration curve for toxin derivatives R13E and R13H, reflecting the effect of changing the charge on residue 13, in the bound state. Experiments with R13A provided a control reflecting the effects of titration of all residues on toxin and channel other than toxin residue 13. The apparent pKs for the titration of residual conductance are shifted 2-3 pH units positive from the nominal pK values for histidine and glutamate, respectively, and from the values for these specific residues, determined in the toxin molecule in free solution by NMR measurements. Toxin affinity also changes dramatically as a function of pH, almost entirely due to changes in the association rate constant, kon. Interpreted electrostatically, our results suggest that, even in the presence of the bound cationic toxin, the channel vestibule strongly favors cation entry with an equivalent local electrostatic potential more negative than -100 mV at the level of the "outer charged ring" formed by channel residues E403, E758, D1241, and D1532. Association rates are apparently limited at a transition state where the pK of toxin residue 13 is closer to the solution value than in the bound state. The action of these unique peptides can thus be used to sense the local environment in the ligand--receptor complex during individual molecular transitions and defined conformational states.     

Protein electron transfer (mechanism and reproductive toxicity): iminium, hydrogen bonding, homoconjugation, amino acid side chains (redox and charged), and cell signaling

This contribution presents novel biochemical perspectives of protein electron transfer (ET) with focus on the iminium nature of the peptide link, along with relationships to reproductive toxicity. The favorable influence of hydrogen bonding on protein ET has been widely documented. Hydrogen bonding of the zwitterionic peptide enhances iminium character. A wide array of such bonding agents is available in vivo, with many reports on the peptide link itself. ET proceeds along the backbone, due in part, to homoconjugation. Redox amino acids (AAs), mainly tyrosine (Tyr), tryptophan (Typ), histidine (His), cysteine (Cys), disulfide, and methionine (Met), are involved in the competing processes for radical formation: direct hydrogen atom abstraction versus electron and proton loss. It appears that the radical or radical cation generated during the redox process is capable of interacting with n-electrons of the backbone. Beneficial effects of cationic AAs impact the conduction process. A relationship apparently exists involving cell signaling, protein conduction, and radicals or electrons. In addition, the link between protein ET and reproductive toxicity is examined. A key element is the role of reactive oxygen species (ROS) generated by protein ET. There is extensive evidence for involvement of ROS in generation of birth defects. The radical species arise in protein mainly by ET transformations by enzymes, as illustrated in the case of alcoholism.