Proton transfer reactions and hydrogen-bond networks in protein environments

Enhanced QM/MM sampling for free energy calculation of chemical reactions: A case study of double proton transfer

Free energy calculations for chemical reactions with a steep energy barrier require well defined reaction coordinates (RCs). However, when multiple parallel channels exist along selected RC, the application of conventional enhanced samplings is difficult to generate correct sampling within limited simulation time and thus cannot give correct prediction about the favorable pathways, the relative stability of multiple products or intermediates. Here, we implement the selective integrated tempering sampling (SITS) method with quantum mechanical and molecular mechanical (QM/MM) potential to investigate the chemical reactions in solution. The combined SITS-QM/MM scheme is used to identify possible reaction paths, intermediate and product states, and the free energy profiles for the different reaction paths. Two double proton transfer reactions were studied to validate the implemented method and simulation protocol, from which the independent and correlated proton transfer processes are identified in two representative systems, respectively. This protocol can be generalized to various kinds of chemical reactions for both academic studies and industry applications, such as in exploration and optimization of potential reactions in DNA encoded compound library and halogen or deuterium substitution of the hit discovery and lead optimization stages of drug design via providing a better understanding of the reaction mechanism along the designed chemical reaction pathways.

Evidence of proton-coupled mixed-valency by electrochemical behavior on transition metal complex dimers bridged by two Ag+ ions

H-Bonded metal complex dimers with reversible redox behaviour, which are connected by a low-barrier hydrogen bond (LBHB) with a very low energy barrier for proton transfer, can provide a unique mixed-valency state stabilized by the proton-coupled electron transfer (PCET) phenomenon. Using cyclic voltammetry measurements, newly prepared [ReIIICl2(PnPr3)2(Hbim)]2 (2) and [OsIIICl2(PnPr3)2(Hbim)]2 (3) existing as H-bonded dimers in a CH2Cl2 solution showed a four-step and four-electron transfer containing two mixed-valency states of ReIIReIII and ReIIIReIV, and OsIIOsIII and OsIIIOsVI, respectively. Furthermore, [ReIIICl2(PnPr3)2(Agbim)]2 (4) and [OsIIICl2(PnPr3)2(Agbim)]2 (5), bridged by two Ag+ ions instead of two H-bonding protons, were prepared, and their electrochemical behaviours changed to a two-step and four-electron transfer. It is clear that the H-bonded complex dimers 2 and 3, connected by an LBHB, can be electrochemically stabilized into unique pairs of mixed-valency states by PCET, and the H-bonding proton transfer also controls the electrochemical redox behaviour.

Reversible on/off switching of photochromic properties in N-salicylideneaniline co-crystals by heating and humidification

The photochromic properties of N-salicylideneaniline (SA) derivative co-crystals were switched on/off through the hydration–dehydration phase transition induced by humidification and heating.

Structural Stability of Peptidic His-Containing Proton Wire in Solution and in the Adsorbed State

Molecular wires are functional molecules applicable in the field of transfer processes in technological and biochemical applications. Besides molecular wires with the ability to transfer electrons, research is currently focused on molecular wires with high proton affinity and proton transfer ability. Recently, proposed peptidic proton wires (H wires) are one example. Their ability to mediate the transport of protons from aqueous solutions onto the surface of a Hg electrode in a catalytic hydrogen evolution reaction was investigated by constant-current chronopotentiometric stripping. However, elucidating the structure of H wires and rationalizing their stability are key requirements for their further research and application. In this article, we focus on the His (H) and Ala (A)-containing peptidic H wire A₃-(H-A₂)₆ in solution and after its immobilization onto the electrode surface in the presence of the secondary structure stabilizer 2,2,2-trifluoroethanol (TFE). We found that the solvent containing more than 25% of TFE stabilizes the helical structure of A₃-(H-A₂)₆ not only in solution but also in the adsorbed state. The TFE efficacy to stabilize α-helical structure was confirmed using high-resolution nuclear magnetic resonance, circular dichroism, and molecular dynamics simulation. Experimental and theoretical results indicated A₃-(H-A₂)₆ to be a high proton-affinity peptidic H wire with an α-helical structure stabilized by TFE, which was confirmed in a comparative study with hexahistidine as an example of a peptide with a definitely disordered and random coil structure. The results presented here could be used for further investigation of the peptidic H wires and for the application of electrochemical methods in the research of proton transfer phenomena in general.

A recurrent point mutation in PRKCA is a hallmark of chordoid gliomas

Chordoid glioma (ChG) is a characteristic, slow growing, and well-circumscribed diencephalic tumor, whose mutational landscape is unknown. Here we report the analysis of 16 ChG by whole-exome and RNA-sequencing. We found that 15 ChG harbor the same PRKCA^D463H mutation. PRKCA encodes the Protein kinase C (PKC) isozyme alpha (PKCα) and is mutated in a wide range of human cancers. However the hot spot PRKCA^D463H mutation was not described in other tumors. PRKCA^D463H is strongly associated with the activation of protein translation initiation (EIF2) pathway. PKCα^D463H mRNA levels are more abundant than wild-type PKCα transcripts, while PKCα^D463H is less stable than the PCKα^WT protein. Compared to PCKα^WT, the PKCα^D463H protein is depleted from the cell membrane. The PKCα^D463H mutant enhances proliferation of astrocytes and tanycytes, the cells of origin of ChG. In conclusion, our study identifies the hallmark mutation for chordoid gliomas and provides mechanistic insights on ChG oncogenesis.

Chordoid glioma is a slow growing diencephalic tumor whose mutational landscape is poorly characterized. Here, the authors perform whole-exome and RNA-sequencing and find that 15 of 16 chordoid glioma cases studied harbor the same PRKCA mutation which results in enhanced proliferation.

Negative influence of pKa on activation energy barrier: A case study for double proton transfer reaction in inorganic acid dimers

Strength of acid can be determined by means of pK_a value. Attempts have been made to find a relationship between pK_a and activation energy barrier for a double proton transfer (DPT) reaction in inorganic acid dimers. Negative influence of pK_a is observed on activation energy (E_a ) which is contrary to the general convention of pK_a . Four different levels of theories with two different basis sets have been used to calculate the activation energy barrier of the DPT reaction in inorganic acid dimers. A model based on first and second order polynomial has been created to find the relationship between activation energy for DPT reaction. © 2018 Wiley Periodicals, Inc.

QM/MM free energy Simulations of an efficient Gluten Hydrolase (Kuma030) Implicate for a Reactant-State Based Protein-Design Strategy for General Acid/Base Catalysis

It is a grand attraction for contemporary biochemists to computationally design enzymes for novel chemical transformation or improved catalytic efficiency. Rosetta by Baker et al. is no doubt the leading software in the protein design society. Generally, optimization of the transition state (TS) is part of the Rosetta’s protocol to enhance the catalytic efficiency of target enzymes, since TS stabilization is the determining factor for catalytic efficiency based on the TS theory (TST). However, it is confusing that optimization of the reactant state (RS) also results in significant improvement of catalytic efficiency in some cases, such as design of gluten hydrolase (Kuma030). Therefore, it is interesting to uncover underlying reason why a better binding in the RS leading to an increased k_cat. In this study, the combined quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) and free energy (PMF) simulations, pK_a calculation, and the statistical analysis such as the ANOVA test were carried out to shed light on the interesting but elusive question. By integration of our computational results and general acid/base theory, we answered the question why optimization of RS stabilization leads to a better TS stabilization in the general acid/base catalysis. In addition, a new and simplified protein-design strategy is proposed for the general acid/base catalysis. The idea, that application of traditional well-defined enzyme mechanism to protein design strategy, would be a great help for methodology development of protein design.

Local energy decomposition analysis of hydrogen-bonded dimers within a domain-based pair natural orbital coupled cluster study

The local energy decomposition (LED) analysis allows for a decomposition of the accurate domain-based local pair natural orbital CCSD(T) [DLPNO-CCSD(T)] energy into physically meaningful contributions including geometric and electronic preparation, electrostatic interaction, interfragment exchange, dynamic charge polarization, and London dispersion terms. Herein, this technique is employed in the study of hydrogen-bonding interactions in a series of conformers of water and hydrogen fluoride dimers. Initially, DLPNO-CCSD(T) dissociation energies for the most stable conformers are computed and compared with available experimental data. Afterwards, the decay of the LED terms with the intermolecular distance (r) is discussed and results are compared with the ones obtained from the popular symmetry adapted perturbation theory (SAPT). It is found that, as expected, electrostatic contributions slowly decay for increasing r and dominate the interaction energies in the long range. London dispersion contributions decay as expected, as r-6. They significantly affect the depths of the potential wells. The interfragment exchange provides a further stabilizing contribution that decays exponentially with the intermolecular distance. This information is used to rationalize the trend of stability of various conformers of the water and hydrogen fluoride dimers.

Voltage and pH sensing by the voltage-gated proton channel, HV1

Voltage-gated proton channels are unique ion channels, membrane proteins that allow protons but no other ions to cross cell membranes. They are found in diverse species, from unicellular marine life to humans. In all cells, their function requires that they open and conduct current only under certain conditions, typically when the electrochemical gradient for protons is outwards. Consequently, these proteins behave like rectifiers, conducting protons out of cells. Their activity has electrical consequences and also changes the pH on both sides of the membrane. Here we summarize what is known about the way these proteins sense the membrane potential and the pH inside and outside the cell. Currently, it is hypothesized that membrane potential is sensed by permanently charged arginines (with very high pKa) within the protein, which results in parts of the protein moving to produce a conduction pathway. The mechanism of pH sensing appears to involve titratable side chains of particular amino acids. For this purpose their pKa needs to be within the operational pH range. We propose a 'counter-charge' model for pH sensing in which electrostatic interactions within the protein are selectively disrupted by protonation of internally or externally accessible groups.

O2 evolution and recovery of the water-oxidizing enzyme

In photosystem II, light-induced water oxidation occurs at the Mn₄CaO₅ cluster. Here we demonstrate proton releases, dioxygen formation, and substrate water incorporation in response to Mn₄CaO₅ oxidation in the protein environment, using a quantum mechanical/molecular mechanical approach and molecular dynamics simulations. In S₂, H₂O at the W1 site forms a low-barrier H-bond with D1-Asp61. In the S₂-to-S₃ transition, oxidation of O_W1H^– to O_W1^•–, concerted proton transfer from O_W1H^– to D1-Asp61, and binding of a water molecule W_n-W1 at O_W1^•– are observed. In S₄, W_n-W1 facilitates oxo-oxyl radical coupling between O_W1^•– and corner μ-oxo O4. Deprotonation via D1-Asp61 leads to formation of O_W1=O4. As O_W1=O4 moves away from Mn, H₂O at W539 is incorporated into the vacant O4 site of the O₂-evolved Mn₄CaO₄ cluster, forming a μ-oxo bridge (Mn3–O_W539–Mn4) in an exergonic process. Simultaneously, W_n-W1 is incorporated as W1, recovering the Mn₄CaO₅ cluster.

Water splitting during photosynthesis results in the combination of two oxygen atoms to form O₂. Here, based on computational simulations, the authors develop a possible mechanism for this reaction, which is different from the mechanisms previous studies have suggested.

The roles of buckled geometry and water environment in the excitonic properties of graphitic C3N4

The exciton plays a crucial role in two-dimensional materials involved in photocatalytic water splitting, where its properties are determined not only by the material itself, but also by the surrounding water environment. By the framework of many-body perturbation theory, we investigated the excitonic effects in pure and water-adsorbed g-C₃N₄. It is shown that the excitonic properties are very sensitive to the geometry of g-C₃N₄ and the adsorption of water molecules. Firstly, the optical band gap, i.e. the first bright excitonic energy of pure g-C₃N₄ decreases remarkably from a high symmetry planar structure (3.8 eV) to a P1 buckled configuration (2.7 eV). Secondly, the hydrogen bonds between water and g-C₃N₄ induce the generation of interface excitons. With a reduced binding energy (at least 0.2 eV), interface excitons can contribute to a more efficient separation of electrons and holes. Our work provides an insight into the excitation mechanism of 2D photocatalysts in a real environment.

Structural Articulation of Biochemical Reactions Using Restrained Geometries and Topology Switching

A strategy named "restrained geometries and topology switching" (RGATS) is presented to obtain detailed trajectories for complex biochemical reactions using molecular mechanics (MM) methods. It enables prediction of realistic dynamical pathways for chemical reactions, especially for accurately characterizing the structural adjustments of highly complex environments to any proximal biochemical reaction. It can be used to generate reactive conformations, model stepwise or concerted reactions in complex environments, and probe the influence of changes in the environment. Its ability to take reactively nonoptimal conformations and generate favorable starting conformations for a biochemical reaction is illustrated for a proton transfer between two model compounds. Its ability to study concerted reactions in explicit solvent is illustrated using proton transfers between an ammonium ion and two conserved histidines in an ammonia transporter channel embedded in a lipid membrane. Its ability to characterize the changes induced by subtle differences in the active site environment is illustrated using nucleotide addition by a DNA polymerase in the presence of two versus three Mg²⁺ ions. RGATS can be employed within any MM program and requires no additional software implementation. This allows the full assortment of computational methods implemented in all available MM programs to be used to tackle virtually any question about biochemical reactions that is answerable without using a quantum mechanical (QM) model. It can also be applied to generate reasonable starting structures for more detailed and expensive QM or QM/MM methods. In particular, this strategy enables rapid prediction of reactant, intermediary, or product state structures in any macromolecular context, with the only requirement being that the structure in any one of these states is either known or can be accurately modeled.

Understanding proton capture and cation-induced dimerization of [Ag29(BDT)12]3−clusters by ion mobility mass spectrometry

We present a unique reaction of [Ag₂₉(BDT)₁₂]³⁻cluster with protons and dimerization of the cluster induced by alkali metal ions.

MolProbity: More and better reference data for improved all-atom structure validation

This paper describes the current update on macromolecular model validation services that are provided at the MolProbity website, emphasizing changes and additions since the previous review in 2010. There have been many infrastructure improvements, including rewrite of previous Java utilities to now use existing or newly written Python utilities in the open-source CCTBX portion of the Phenix software system. This improves long-term maintainability and enhances the thorough integration of MolProbity-style validation within Phenix. There is now a complete MolProbity mirror site at http://molprobity.manchester.ac.uk. GitHub serves our open-source code, reference datasets, and the resulting multi-dimensional distributions that define most validation criteria. Coordinate output after Asn/Gln/His "flip" correction is now more idealized, since the post-refinement step has apparently often been skipped in the past. Two distinct sets of heavy-atom-to-hydrogen distances and accompanying van der Waals radii have been researched and improved in accuracy, one for the electron-cloud-center positions suitable for X-ray crystallography and one for nuclear positions. New validations include messages at input about problem-causing format irregularities, updates of Ramachandran and rotamer criteria from the million quality-filtered residues in a new reference dataset, the CaBLAM Cα-CO virtual-angle analysis of backbone and secondary structure for cryoEM or low-resolution X-ray, and flagging of the very rare cis-nonProline and twisted peptides which have recently been greatly overused. Due to wide application of MolProbity validation and corrections by the research community, in Phenix, and at the worldwide Protein Data Bank, newly deposited structures have continued to improve greatly as measured by MolProbity's unique all-atom clashscore.

Neutron crystallography of photoactive yellow protein reveals unusual protonation state of Arg52 in the crystal

Because of its high pK_a, arginine (Arg) is believed to be protonated even in the hydrophobic environment of the protein interior. However, our neutron crystallographic structure of photoactive yellow protein, a light sensor, demonstrated that Arg52 adopts an electrically neutral form. We also showed that the hydrogen bond between the chromophore and Glu46 is a so-called low barrier hydrogen bond (LBHB). Because both the neutral Arg and LBHB are unusual in proteins, these observations remain controversial. To validate our findings, we carried out neutron crystallographic analysis of the E46Q mutant of PYP. The resultant structure revealed that the proportion of the cationic form is higher in E46Q than in WT, although the cationic and neutral forms of Arg52 coexist in E46Q. These observations were confirmed by the occupancy of the deuterium atom bound to the N _η1 atom combined with an alternative conformation of the N_(η2)D₂ group comprising sp² hybridisation. Based on these results, we propose that the formation of the LBHB decreases the proton affinity of Arg52, stabilizing the neutral form in the crystal.

High Pressure Spectroscopic Investigation on Proton Transfer in Squaric Acid and 4,4'-Bipyridine Co-crystal

In attempt to the obtain detailed geometric information of proton transfer compound (subsequently denote as SQBP) formed between squaric acid (SQ)and 4,4′-bipyridine(BP), and to investigate the mechanisms of pressure-induced double proton transfer and related structural phase transition, we carried out in-situ high pressure Raman spectroscopy of SQBP up to 20 GPa. A solid-solid phase transition together with double proton transfer phenomenon was confirmed by Raman spectroscopy at about 1.5 GPa, and the activation of C = O stretching mode in Raman spectra indicates a square-ring structure of SQ with four symmetric C = O bond formation. These results are further supported by first-principals calculations and in-situ high pressure infrared absorption spectroscopy. Additionally, Raman intensity analysis suggests that a higher-order phase transition with planar BP molecular structure occurred in the pressure range of 3~6 GPa. As a result, the π electron delocalization effect in BP dominated the intensity enhancement of C = O stretching mode in SQ. To the best of our knowledge, this is the first time observation of the intensity enhancement of proton donor’s normal modes induced by proton acceptor’s π electron delocalization.

The Existence of an Isolated Hydronium Ion in the Interior of Proteins

Neutron diffraction analysis studies reported an isolated hydronium ion (H₃O⁺) in the interior of d‐xylose isomerase (XI) and phycocyanobilin‐ferredoxin oxidoreductase (PcyA). H₃O⁺ forms hydrogen bonds (H‐bonds) with two histidine side‐chains and a backbone carbonyl group in PcyA, whereas H₃O⁺ forms H‐bonds with three acidic residues in XI. Using a quantum mechanical/molecular mechanical (QM/MM) approach, we analyzed stabilization of H₃O⁺ by the protein environment. QM/MM calculations indicated that H₃O⁺ was unstable in the PcyA crystal structure, releasing a proton to an H‐bond partner His88, producing H₂O and protonated His88. On the other hand, H₃O⁺ was stable in the XI crystal structure. H‐bond partners of isolated H₃O⁺ would be practically limited to acidic residues such as aspartic and glutamic acids in the protein environment.

NMR Spectroscopic Characterization of Charge Assisted Strong Hydrogen Bonds in Brønsted Acid Catalysis

Hydrogen bonding plays a crucial role in Brønsted acid catalysis. However, the hydrogen bond properties responsible for the activation of the substrate are still under debate. Here, we report an in depth study of the properties and geometries of the hydrogen bonds in (R)-TRIP imine complexes (TRIP: 3,3'-Bis(2,4,6-triisopropylphenyl)-1,1'-binaphthyl-2,2'-diylhydrogen phosphate). From NMR spectroscopic investigations 1H and 15N chemical shifts, a Steiner-Limbach correlation, a deuterium isotope effect as well as quantitative values of 1JNH,2hJPH and 3hJPN were used to determine atomic distances (rOH, rNH, rNO) and geometry information. Calculations at SCS-MP2/CBS//TPSS-D3/def2-SVP-level of theory provided potential surfaces, atomic distances and angles. In addition, scalar coupling constants were computed at TPSS-D3/IGLO-III. The combined experimental and theoretical data reveal mainly ion pair complexes providing strong hydrogen bonds with an asymmetric single well potential. The geometries of the hydrogen bonds are not affected by varying the steric or electronic properties of the aromatic imines. Hence, the strong hydrogen bond reduces the degree of freedom of the substrate and acts as a structural anchor in the (R)-TRIP imine complex.

The Low Barrier Hydrogen Bond in the Photoactive Yellow Protein: A Vacuum Artifact Absent in the Crystal and Solution

There has been considerable debate on the existence of a low-barrier hydrogen bond (LBHB) in the photoactive yellow protein (PYP). The debate was initially triggered by the neutron diffraction study of Yamaguchi et al. ( Proc. Natl. Acad. Sci., U. S. A. , 2009 , 106 , 440 - 444 ) who suggested a model in which a neutral Arg52 residue triggers the formation of the LBHB in PYP. Here, we present an alternative model that is consistent within the error margins of the Yamaguchi structure factors. The model explains an increased hydrogen bond length without nuclear quantum effects and for a protonated Arg52. We tested both models by calculations under crystal, solution, and vacuum conditions. Contrary to the common assumption in the field, we found that a single PYP in vacuum does not provide an accurate description of the crystal conditions but instead introduces strong artifacts, which favor a LBHB and a large ¹H NMR chemical shift. Our model of the crystal environment was found to stabilize the two Arg52 hydrogen bonds and crystal water positions for the protonated Arg52 residue in free MD simulations and predicted an Arg52 pK_a upshift with respect to PYP in solution. The crystal and solution environments resulted in almost identical ¹H chemical shifts that agree with NMR solution data. We also calculated the effect of the Arg52 protonation state on the LBHB in 3D nuclear equilibrium density calculations. Only the charged crystal structure in vacuum supports a LBHB if Arg52 is neutral in PYP at the previously reported level of theory ( J. Am. Chem. Soc. , 2014 , 136 , 3542 - 3552 ). We attribute the anomalies in the interpretation of the neutron data to a shift of the potential minimum, which does not involve nuclear quantum effects and is transferable beyond the Yamaguchi structure.

Mechanisms for Two-Step Proton Transfer Reactions in the Outward-Facing Form of MATE Transporter

Bacterial pathogens or cancer cells can acquire multidrug resistance, which causes serious clinical problems. In cells with multidrug resistance, various drugs or antibiotics are extruded across the cell membrane by multidrug transporters. The multidrug and toxic compound extrusion (MATE) transporter is one of the five families of multidrug transporters. MATE from Pyrococcus furiosus uses H(+) to transport a substrate from the cytoplasm to the outside of a cell. Crystal structures of MATE from P. furiosus provide essential information on the relevant H(+)-binding sites (D41 and D184). Hybrid quantum mechanical/molecular mechanical simulations and continuum electrostatic calculations on the crystal structures predict that D41 is protonated in one structure (Straight) and, both D41 and D184 protonated in another (Bent). All-atom molecular dynamics simulations suggest a dynamic equilibrium between the protonation states of the two aspartic acids and that the protonation state affects hydration in the substrate binding cavity and lipid intrusion in the cleft between the N- and C-lobes. This hypothesis is examined in more detail by quantum mechanical/molecular mechanical calculations on snapshots taken from the molecular dynamics trajectories. We find the possibility of two proton transfer (PT) reactions in Straight: the 1st PT takes place between side-chains D41 and D184 through a transient formation of low-barrier hydrogen bonds and the 2nd through another H(+) from the headgroup of a lipid that intrudes into the cleft resulting in a doubly protonated (both D41 and D184) state. The 1st PT affects the local hydrogen bond network and hydration in the N-lobe cavity, which would impinge on the substrate-binding affinity. The 2nd PT would drive the conformational change from Straight to Bent. This model may be applicable to several prokaryotic H(+)-coupled MATE multidrug transporters with the relevant aspartic acids.

Atomistic determinants of co-enzyme Q reduction at the Qi-site of the cytochrome bc1 complex

The cytochrome (cyt) bc₁ complex is an integral component of the respiratory electron transfer chain sustaining the energy needs of organisms ranging from humans to bacteria. Due to its ubiquitous role in the energy metabolism, both the oxidation and reduction of the enzyme’s substrate co-enzyme Q has been studied vigorously. Here, this vast amount of data is reassessed after probing the substrate reduction steps at the Q_i-site of the cyt bc₁ complex of Rhodobacter capsulatus using atomistic molecular dynamics simulations. The simulations suggest that the Lys251 side chain could rotate into the Q_i-site to facilitate binding of half-protonated semiquinone – a reaction intermediate that is potentially formed during substrate reduction. At this bent pose, the Lys251 forms a salt bridge with the Asp252, thus making direct proton transfer possible. In the neutral state, the lysine side chain stays close to the conserved binding location of cardiolipin (CL). This back-and-forth motion between the CL and Asp252 indicates that Lys251 functions as a proton shuttle controlled by pH-dependent negative feedback. The CL/K/D switching, which represents a refinement to the previously described CL/K pathway, fine-tunes the proton transfer process. Lastly, the simulation data was used to formulate a mechanism for reducing the substrate at the Q_i-site.

Proton transfer aiding phase transitions in oxalic acid dihydrate under pressure

Oxalic acid dihydrate, an important molecular solid in crystal chemistry, ecology and physiology, has been studied for nearly 100 years now. The most debated issues regarding its proton dynamics have arisen due to an unusually short hydrogen bond between the acid and water molecules. Using combined in situ spectroscopic studies and first-principles simulations at high pressures, we show that the structural modification associated with this hydrogen bond is much more significant than ever assumed. Initially, under pressure, proton migration takes place along this strong hydrogen bond at a very low pressure of 2 GPa. This results in the protonation of water with systematic formation of dianionic oxalate and hydronium ion motifs, thus reversing the hydrogen bond hierarchy in the high pressure phase II. The resulting hydrogen bond between a hydronium ion and a carboxylic group shows remarkable strengthening under pressure, even in the pure ionic phase III. The loss of cooperativity of hydrogen bonds leads to another phase transition at ∼ 9 GPa through reorientation of other hydrogen bonds. The high pressure phase IV is stabilized by a strong hydrogen bond between the dominant CO2 and H2O groups of oxalate and hydronium ions, respectively. These findings suggest that oxalate systems may provide useful insights into proton transfer reactions and assembly of simple molecules under extreme conditions.

H-Bond: Τhe Chemistry-Biology H-Bridge

H-bonding, as a non covalent stabilizing interaction of diverse nature, has a central role in the structure, function and dynamics of chemical and biological processes, pivotal to molecular recognition and eventually to drug design. Types of conventional and non conventional (H-H, dihydrogen, H- π, CH- π, anti- , proton coordination and H-S) H-bonding interactions are discussed as well as features emerging from their interplay, such as cooperativity (σ- and π-) effects and allostery. Its utility in many applications is described. Catalysis, proton and electron transfer processes in various materials or supramolecular architectures of preorganized hosts for guest binding, are front-line technology. The H-bond-related concept of proton transfer (PT) addresses energy issues or deciphering the mechanism of many natural and synthetic processes. PT is also of paramount importance in the functions of cells and is assisted by large complex proteins embedded in membranes. Both intermolecular and intramolecular PT in H-bonded systems has received attention, theoretically and experimentally, using prototype molecules. It is found in rearrangement reactions, protein functions, and enzyme reactions or across proton channels and pumps. Investigations on the competition between intra- and intermolecular H bonding are discussed. Of particular interest is the H-bond furcation, a common phenomenon in protein-ligand binding. Multiple H-bonding (H-bond furcation) is observed in supramolecular structures.

Intramolecular, Exciplex-Mediated, Proton-Coupled, Charge-Transfer Processes in N,N-Dimethyl-3-(1-pyrenyl)propan-1-ammonium Cations: Influence of Anion, Solvent Polarity, and Temperature

An intramolecular exciplex-mediated, proton-coupled, charge-transfer (PCCT) process has been investigated for a series of N,N-dimethyl-3-(1-pyrenyl)propan-1-ammonium cations with different anions (PyS) in solvents of low to intermediate polarity over a wide temperature range. Solvent mediates both the equilibrium between conformations of the cation that place the pyrenyl and ammonium groups in proximity (conformation C) or far from each other (conformation O) and the ability of the ammonium group to transfer a proton adiabatically in the PyS excited singlet state. Thus, exciplex emission, concurrent with the PCCT process, was observed only in hydrogen-bond accepting solvents of relatively low polarity (tetrahydrofuran, ethyl acetate, and 1,4-dioxane) and not in dichloromethane. From the exciplex emission and other spectroscopic and thermodynamic data, the acidity of the ammonium group in conformation C of the excited singlet state of PyS (pKa*) has been estimated to be ca. -3.4 in tetrahydrofuran. The ratios between the intensities of emission from the exciplex and the locally excited state (IEx/ILE) appear to be much more dependent on the nature of the anion than are the rates of exciplex formation and decay, although the excited state data do not provide a quantitative measure of the anion effect on the C-O equilibrium. The activation energies associated with exciplex formation in THF are calculated to be 0.08 to 0.15 eV lower than for the neutral amine, N,N-dimethyl-3-(1-pyrenyl)propan-1-amine. Decay of the exciplexes formed from the deprotonation of PyS is hypothesized to occur through charge-recombination processes. To our knowledge, this is the first example in which photoacidity and intramolecular exciplex formation (i.e., a PCCT reaction) are coupled.

pK(a) of a Proton-Conducting Water Chain in Photosystem II

Recent high-resolution crystal structures of the water-oxidizing enzyme photosystem II (PSII) show that O4 of the catalytic Mn4CaO5 cluster forms an H-bond with a water molecule W539, which belongs to a chain of water molecules (O4-water chain). Oxidation of Mn4CaO5 to S1 resulted in elongation of the O-H bonds and decrease in pKa(O-H/O(-)) in the [O4-H···OW539-H···OW538-H···OW393] region along the O4-water chain. In S1, removal of all water molecules from the O4-water chain, except W539, resulted in a significant pKa upshift at O4; this suggests that the proton-conducting water chain serves as a conducting media for protons and significantly decreases the donor pKa, leading to a downhill proton transfer. The absence of a corresponding proton-conducting channel is disadvantageous for release of protons from the proton-releasing site, as in the case of O5 that has no H-bond partner.

Experimental Binding Energies in Supramolecular Complexes

On the basis of many literature measurements, a critical overview is given on essential noncovalent interactions in synthetic supramolecular complexes, accompanied by analyses with selected proteins. The methods, which can be applied to derive binding increments for single noncovalent interactions, start with the evaluation of consistency and additivity with a sufficiently large number of different host-guest complexes by applying linear free energy relations. Other strategies involve the use of double mutant cycles, of molecular balances, of dynamic combinatorial libraries, and of crystal structures. Promises and limitations of these strategies are discussed. Most of the analyses stem from solution studies, but a few also from gas phase. The empirically derived interactions are then presented on the basis of selected complexes with respect to ion pairing, hydrogen bonding, electrostatic contributions, halogen bonding, π-π-stacking, dispersive forces, cation-π and anion-π interactions, and contributions from the hydrophobic effect. Cooperativity in host-guest complexes as well as in self-assembly, and entropy factors are briefly highlighted. Tables with typical values for single noncovalent free energies and polarity parameters are in the Supporting Information.

Biexponential photon antibunching: recombination kinetics within the Förster-cycle in DMSO

Time-resolved experiments with pulsed-laser excitation are the standard approach to map the dynamic evolution of excited states, but ground-state kinetics remain hidden or require pump-dump-probe schemes. Here, we exploit the so-called photon antibunching, a purely quantum-optical effect related to single molecule detection to assess the rate constants for a chemical reaction in the electronic ground state. The measurement of the second-order correlation function g((2)), i.e. the evaluation of inter-photon arrival times, is applied to the reprotonation in a Förster-cycle. We find that the antibunching of three different photoacids in the aprotic solvent DMSO significantly differs from the behavior in water. The longer decay constant of the biexponential antibunching tl is linked to the bimolecular reprotonation kinetics of the fully separated ion-pair, independent of the acidic additives. The value of the corresponding bimolecular rate constant, kp = 4 × 10(9) M(-1) s(-1), indicates diffusion-controlled reprotonation. The analysis of tl also allows for the extraction of the separation yield of proton and the conjugated base after excitation and amounts to approximately 15%. The shorter time component ts is connected to the decay of the solvent-separated ion pair. The associated time constant for geminate reprotonation is approximately 3 ± 1 ns in agreement with independent tcspc experiments. These experiments verify that the transfer of quantum-optical experiments to problems in chemistry enables mechanistic conclusions which are hardly accessible by other methods.

Effect of the Hydrogen Bond on Photochemical Synthesis of Silver Nanoparticles

The effect of a hydrogen bond on the photochemical synthesis of silver nanoparticles has been investigated via experimental and theoretical methods. In a benzophenone system, the photochemical synthesis process includes two steps, which are that hydrogen abstraction reaction and the following reduction reaction. We found that for the first step, an intermolecular hydrogen bond enhances the proton transfer. The efficiency of hydrogen abstraction increases with the hydrogen bond strength. For the second step, the hydrogen-bonded ketyl radical complex shows higher reducibility than the ketyl radical. The inductively coupled plasma-optical emission spectroscopy (ICP-OES) measurement exhibits a 2.49 times higher yield of silver nanoparticles in the hydrogen bond ketyl radical complex system than that for the ketyl radical system. Theoretical calculations show that the hydrogen bond accelerates electron transfer from the ketyl radical to the silver ion by raising the SOMO energy of the ketyl radical; thus, the SOMO-LUMO interaction is more favorable.

Energetics of proton release on the first oxidation step in the water-oxidizing enzyme

In photosystem II (PSII), the Mn₄CaO₅ cluster catalyses the water splitting reaction. The crystal structure of PSII shows the presence of a hydrogen-bonded water molecule directly linked to O4. Here we show the detailed properties of the H-bonds associated with the Mn₄CaO₅ cluster using a quantum mechanical/molecular mechanical approach. When O4 is taken as a μ-hydroxo bridge acting as a hydrogen-bond donor to water539 (W539), the S₀ redox state best describes the unusually short O4–O_W539 distance (2.5 Å) seen in the crystal structure. We find that in S₁, O4 easily releases the proton into a chain of eight strongly hydrogen-bonded water molecules. The corresponding hydrogen-bond network is absent for O5 in S₁. The present study suggests that the O4-water chain could facilitate the initial deprotonation event in PSII. This unexpected insight is likely to be of real relevance to mechanistic models for water oxidation.

The availability of crystal structures of photosystem II opens up the possibility of gaining insights into its mechanism. Here, the authors use a computational approach and propose a deprotonation event at O4 followed by long-range proton-transfer along a chain of strongly bonded water molecules.

Accurate ab initio calculations of O-HO and O-H(-)O proton chemical shifts: towards elucidation of the nature of the hydrogen bond and prediction of hydrogen bond distances

The inability to determine precisely the location of labile protons in X-ray molecular structures has been a key barrier to progress in many areas of molecular sciences. We report an approach for predicting hydrogen bond distances beyond the limits of X-ray crystallography based on accurate ab initio calculations of O-HO proton chemical shifts, using a combination of DFT and contactor-like polarizable continuum model (PCM). Very good linear correlation between experimental and computed (at the GIAO/B3LYP/6-311++G(2d,p) level of theory) chemical shifts were obtained with a large set of 43 compounds in CHCl3 exhibiting intramolecular O-HO and intermolecular and intramolecular ionic O-H(-)O hydrogen bonds. The calculated OH chemical shifts exhibit a strong linear dependence on the computed (O)HO hydrogen bond length, in the region of 1.24 to 1.85 Å, of -19.8 ppm Å(-1) and -20.49 ppm Å(-1) with optimization of the structures at the M06-2X/6-31+G(d) and B3LYP/6-31+G(d) level of theory, respectively. A Natural Bond Orbitals (NBO) analysis demonstrates a very good linear correlation between the calculated (1)H chemical shifts and (i) the second-order perturbation stabilization energies, corresponding to charge transfer between the oxygen lone pairs and σ antibonding orbital and (ii) Wiberg bond order of the O-HO and O-H(-)O hydrogen bond. Accurate ab initio calculations of O-HO and O-H(-)O (1)H chemical shifts can provide improved structural and electronic description of hydrogen bonding and a highly accurate measure of distances of short and strong hydrogen bonds.

A Practical Quantum Mechanics Molecular Mechanics Method for the Dynamical Study of Reactions in Biomolecules

Quantum mechanics/molecular mechanics (QM/MM) methods are excellent tools for the modeling of biomolecular reactions. Recently, we have implemented a new QM/MM method (Fireball/Amber), which combines an efficient density functional theory method (Fireball) and a well-recognized molecular dynamics package (Amber), offering an excellent balance between accuracy and sampling capabilities. Here, we present a detailed explanation of the Fireball method and Fireball/Amber implementation. We also discuss how this tool can be used to analyze reactions in biomolecules using steered molecular dynamics simulations. The potential of this approach is shown by the analysis of a reaction catalyzed by the enzyme triose-phosphate isomerase (TIM). The conformational space and energetic landscape for this reaction are analyzed without a priori assumptions about the protonation states of the different residues during the reaction. The results offer a detailed description of the reaction and reveal some new features of the catalytic mechanism. In particular, we find a new reaction mechanism that is characterized by the intramolecular proton transfer from O1 to O2 and the simultaneous proton transfer from Glu 165 to C2.

Unraveling ultrafast photoinduced proton transfer dynamics in a fluorescent protein biosensor for Ca(2+) imaging

Imaging Ca(2+) dynamics in living systems holds great potential to advance neuroscience and cellular biology. G-GECO1.1 is an intensiometric fluorescent protein Ca(2+) biosensor with a Thr-Tyr-Gly chromophore. The protonated chromophore emits green upon photoexcitation via excited-state proton transfer (ESPT). Upon Ca(2+) binding, a significant population of the chromophores becomes deprotonated. It remains elusive how the chromophore structurally evolves prior to and during ESPT, and how it is affected by Ca(2+) . We use femtosecond stimulated Raman spectroscopy to dissect ESPT in both the Ca(2+) -free and bound states. The protein chromophores exhibit a sub-200 fs vibrational frequency shift due to coherent small-scale proton motions. After wavepackets move out of the Franck-Condon region, ESPT gets faster in the Ca(2+) -bound protein, indicative of the formation of a more hydrophilic environment. These results reveal the governing structure-function relationship of Ca(2+) -sensing protein biosensors.

Excited state structural events of a dual-emission fluorescent protein biosensor for Ca²⁺ imaging studied by femtosecond stimulated Raman spectroscopy

Fluorescent proteins (FPs) are luminescent biomolecules that emit characteristic hues upon irradiation. A group of calmodulin (CaM)-green FP (GFP) chimeras have been previously engineered to enable the optical detection of calcium ions (Ca(2+)). We investigate one of these genetically encoded Ca(2+) biosensors for optical imaging (GECOs), GEM-GECO1, which fluoresces green without Ca(2+) but blue with Ca(2+), using femtosecond stimulated Raman spectroscopy (FSRS). The time-resolved FSRS data (<800 cm(-1)) reveal that initial structural evolution following 400 nm photoexcitation involves small-scale coherent proton motions on both ends of the chromophore two-ring system with a <250 fs time constant. Upon Ca(2+) binding, the chromophore adopts a more twisted conformation in the protein pocket with increased hydrophobicity, which inhibits excited-state proton transfer (ESPT) by effectively trapping the protonated chromophore in S1. Both the chromophore photoacidity and local environment form the ultrafast structural dynamics basis for the dual-emission properties of GEM-GECO1. Its photochemical transformations along multidimensional reaction coordinates are evinced by distinct stages of FSRS spectral evolution, particularly related to the ∼460 and 504 cm(-1) modes. The direct observation of lower frequency modes provides crucial information about the nuclear motions preceding ESPT, which enriches our understanding of photochemistry and enables the rational design of new biosensors.

Revealing the functional states in the active site of BLUF photoreceptors from electrochromic shift calculations

Photoexcitation with blue light of the flavin chromophore in BLUF photoreceptors induces a switch into a metastable signaling state that is characterized by a red-shifted absorption maximum. The red shift is due to a rearrangement in the hydrogen bond pattern around Gln63 located in the immediate proximity of the isoalloxazine ring system of the chromophore. There is a long-lasting controversy between two structural models, named Q63A and Q63J in the literature, on the local conformation of the residues Gln63 and Tyr21 in the dark state of the photoreceptor. As regards the mechanistic details of the light-activation mechanism, rotation of Gln63 is opposed by tautomerism in the Q63A and Q63J models, respectively. We provide a structure-based simulation of electrochromic shifts of the flavin chromophore in the wild type and in various site-directed mutants. The excellent overall agreement between experimental and computed data allows us to evaluate the two structural models. Compelling evidence is obtained that the Q63A model is incorrect, whereas the Q63J is fully consistent with the present computations. Finally, we confirm independently that a keto-enol tautomerization of the glutamine at position 63, which was proposed as molecular mechanism for the transition between the dark and the light-adapted state, explains the measured 10 to 15 nm red shift in flavin absorption between these two states of the protein. We believe that the accurateness of our results provides evidence that the BLUF photoreceptors absorption is fine-tuned through electrostatic interactions between the chromophore and the protein matrix, and finally that the simplicity of our theoretical model is advantageous as regards easy reproducibility and further extensions.