Periodicity in proton conduction along a H‐bonded chain. Application to biomolecules

Molecular complexes are constructed to simulate proton transfer channels of the influenza A virus and of the active site of carbonic anhydrase. These complexes consist of proton donor and acceptor groups connected by a chain of water molecules. Quantum chemical calculations on the methylimidazole(H⁺)? H₂O? CH₃COO^? model of the M2 virus channel indicate free translational motion of the water molecule between donor and acceptor, as well as concerted transfer of both H‐bond protons. The proton transfer barrier does not depend on the position of the bridged water molecule and varies linearly with the difference of electrostatic potentials between the donor and acceptor. When the water chain is elongated, and with various donor and acceptor models, periodicity appears in the H‐bond lengths and the progression of proton transfer in each link. This “wave” is shown to propagate along the chain, as it is driven by the displacement of a single proton. One can thereby estimate the velocity of the proton wave and proton conduction time. Computations are performed to examine the influence of immersing the system within a polarizable medium. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Exploration on regulating factors for proton transfer along hydrogen-bonded water chains.

Proton transfer along a single-file hydrogen-bonded water chain is elucidated with a special emphasis on the investigation of chain length, side water, and solvent effects, as well as the temperature and pressure dependences. The number of water molecules in the chain varies from one to nine. The proton can be transported to the acceptor fragment through the single-file hydrogen-bonded water wire which contains at most five water molecules. If the number of water molecule is more than five, the proton is trapped by the chain in the hydroxyl-centered H(7)O(3) (+) state. The farthest water molecule involved in the formation of H(7)O(3) (+) is the fifth one away from the donor fragment. These phenomena reappear in the molecular dynamics simulations. The energy of the system is reduced along with the proton conduction. The proton transfer mechanism can be altered by excess proton. The augmentation of the solvent dielectric constant weakens the stability of the system, but favors the proton transfer. NMR spin-spin coupling constants can be used as a criterion in judging whether the proton is transferred or not. The enhancement of temperature increases the thermal motion of the molecule, augments the internal energy of the system, and favors the proton transfer. The lengthening of the water wire increases the entropy of the system, concomitantly, the temperature dependence of the Gibbs free energy increases. The most favorable condition for the proton transfer along the H-bonded water wire is the four-water contained chain with side water attached near to the acceptor fragment in polar solvent under higher temperature.     

Interstitial water and the formation of low barrier hydrogen bonds: A computational model study

The B3LYP/D95+(d,p) analysis of the uncharged low barrier hydrogen bond (LBHB) between 4‐methyl‐1H‐imidazole (Mim) and acetic acid (HAc) shows that uncharged LBHBs can be formed either by adding three water molecules around the cluster or by placing the Mim–HAc pair in a dielectric environment created by a polarizable continuum model with a permittivity larger than 20.7. The permittivity of environment around uncharged LBHB can be lowered significantly by including water molecules into the system. A Mim–HAc LBHB stabilized with one water molecule observed in diethyl ether (ε = 4.34), with two water molecules in toluene (ε = 2.38), and with three water molecules in vacuo (ε = 1). Solvation models with different numbers of water molecules predict average differences in the proton affinities of the hydrogen bonded bases (ΔPA) for stable uncharged LBHB systems in vacuo to be 91.5 kcal/mol being different from the ΔPA values close to zero in charge‐assisted LBHB systems. The results clearly indicate that small amounts of interstitial water molecules at the active site of enzymes do not preclude the existence of LBHBs in biological catalysis. Our results also show that interstitial water molecules provide a useful clue in the search for uncharged LBHBs in an enzymatic environment and the number of water molecules can be used as a relative measure for the polarity around the direct environment of LBHBs. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

Potential Proton‐Release Channels in Bacteriorhodopsin

The protein bacteriorhodopsin pumps protons across a bacterial membrane; its pumping cycle is triggered by the photoisomerization of a retinal cofactor and involves multiple proton‐transfer reactions between intermittent protonation sites. These transfers are either direct or mediated by hydrogen‐bonded networks, which may include internal water molecules. The terminal step of the proton‐transfer sequence is the proton release from a pocket near Glu194 and Glu204 to the extracellular bulk during the transition from the L to the M photointermediate states. The polar and charged side chains connecting these two regions in the crystal structures show no structural changes between the initial bR state and the L/M states, and no intermittent protonation changes have been detected so far in this region. Based on biomolecular simulations, we propose two potential proton‐release channels, which connect the release pocket to the extracellular medium. In simulations of the L photointermediate we observe bulk water entering these channels and forming transient hydrogen‐bonded networks, which could serve as fast deprotonation pathways from the release pocket to the bulk via a Grotthuss mechanism. For the first channel, we find that the triple Arg7, Glu9, and Tyr79 acts as a valve, thereby gating water uptake and release. The second channel has two release paths, which split at the position Asn76/Pro77 underneath the release group. Here, water molecules either exchange directly with the bulk or diffuse within the protein towards Arg 134/Lys129, where the exchange with the bulk occurs.     

A computational study of ultrafast acid dissociation and acid-base neutralization reactions. II. The relationship between the coordination state of solvent molecules and concerted versus sequential acid dissociation

We investigate the role played by the coordination state of pre-existing water wires during the dissociation of moderately strong acids by means of first-principles molecular dynamics calculations. By preparing 2,4,6-tricyanophenol (calc. pKa～0.5) in two different initial states, we are able to observe sequential as well as concerted trajectories of dissociation: On one hand, equilibrium dissociation takes place on a ～50 ps timescale; proton conduction occurs through three-coordinated water wires in this case, by means of sequential Grotthus hopping. On the other hand, by preparing 2,4,6-tricyanophenol in a hydration state inherited from that of equilibrated phenol (calc. pKa=7.6), the moderately strong acid finds itself in a presolvated state from which dissociation can take place on a ～1 ps timescale. In this case, concerted dissociation trajectories are observed, which consist of proton translocation through two intervening, four-coordinated, water molecules in 0.1-1.0 ps. The present results suggest that, in general, the mechanism of proton translocation depends on how the excess proton is injected into a hydrogen bond network. In particular, if the initial conditions favour proton release to a fourfold H-bonded water molecule, proton translocation by as much as 6-8 A? can take place on a sub-picosecond timescale.     

Cooperative interactions of hydrogen bonds in proton-transfer processes involving water molecules. Simulation of biochemical systems

Quantum-chemical calculations of molecular complexes simulating the proton channel of influenza A virus and the proton-transfer system of the active site of carboanhydrase enzyme were performed. These complexes comprise a proton-donor and a proton-acceptor groups bridged by a chain of water molecules. Calculations of the methylimidazole (H⁺)-H₂O-CH₃COO^? complex as a model of influenza M₂ virus revealed free translation motion of the water molecule between the donor and acceptor, as well as concerted proton transfer in both H bonds. The barrier for proton transfer is independent of the position of the bridging water molecule and varies linearly with the difference in the electrostatic potentials between the donor and acceptor. With elongation of the H-bond bridge between the donor and acceptor groups, the H-bond lengths and proton shifts in the chain links vary periodically. This process can be defined as an H-bond deformation wave (proton wave). It was shown that motion of one proton along the H bond is associated with vibrational motion of protons in other links, which results in wave propagation along the chain. The calculation results allowed the rate of the proton wave and the time of proton transfer from the donor to acceptor to be estimated.     

How Does a Membrane Protein Achieve a Vectorial Proton Transfer Via Water Molecules?

We present a detailed mechanism for the proton transfer from a protein‐bound protonated water cluster to the bulk water directed by protein side chains in the membrane protein bacteriorhodopsin. We use a combined approach of time‐resolved Fourier transform infrared spectroscopy, molecular dynamics simulations, and X‐ray structure analysis to elucidate the functional role of a hydrogen bond between Ser193 and Glu204. These two residues seal the internal protonated water cluster from the bulk water and the protein surface. During the photocycle of bacteriorhodopsin, a transient protonation of Glu204 leads to a breaking of this hydrogen bond. This breaking opens the gate to the extracellular bulk water, leading to a subsequent proton release from the protonated water cluster. We show in detail how the protein achieves vectorial proton transfer via protonated water clusters in contrast to random proton transfer in liquid water.     

Theoretical analysis of pyridine protonation in water clusters of increasing size.

The protonation of pyridine in water clusters as a function of the number of water molecules was theoretically analyzed as a prototypical case for the protonation of organic bases. We determined the variation of structural, bonding, and energetic properties on protonation, as well as the stabilization of the ionic species formed. Thus, we used supermolecular models in which pyridine interacts with clusters of up to five water molecules. For each complex, we determined the most stable unprotonated and protonated structures from a simulated annealing at the semi ab initio level. The structures were optimized at the B3LYP/cc-pVDZ level. We found that the hydroxyl group formed on protonation of pyridine abstracts a proton from the ortho-carbon atom of the pyridine ring. The "atoms in molecules" theory showed that this C-H group loses its covalent character. However, starting with clusters of four water molecules, the C-H bond recovers its covalent nature. This effect is associated with the presence of more than one ring between the water molecules and pyridine. These rings stabilize, by delocalization, the negative charge on the hydroxyl oxygen atom. Considering the protonation energy, we find that the protonated forms are increasingly stabilized with increasing size of the water cluster. When zero-point energy is included, the variation follows closely an exponential decrease with increasing number of water molecules. Analysis of the vibrational modes for the strongest bands in the IR spectra of the complexes suggests that the protonation of pyridine occurs by concerted proton transfers among the different water rings in the structure. Symmetric water stretching was found to be responsible for hydrogen transfer from the water molecule to the pyridine nitrogen atom.     

Different catalysis role of in‐loop and out‐of‐loop waters in assisting HNS/HSN proton transfer isomerizations: Bridging vs. surrounding effect

In this work, a density function theory (DFT) study is presented for the HNS/HSN isomerization assisted by 1–4 water molecules on the singlet state potential energy surface (PES). Two modes are considered to model the catalytic effect of these water molecules: (i) water molecule(s) participate directly in forming a proton transfer loop with HNS/HSN species, and (ii) water molecules are out of loop (referred to as out‐of‐loop waters) to assist the proton transfer. In the first mode, for the monohydration mechanism, the heat of reaction is 21.55 kcal · mol^?1 at the B3LYP/6‐311++G** level. The corresponding forward/backward barrier lowerings are obtained as 24.41/24.32 kcal · mol^?1 compared with the no‐water‐assisting isomerization barrier T (65.52/43.87 kcal · mol^?1). But when adding one water molecule on the HNS, there is another special proton‐transfer isomerization pathway with a transition state 10T′ in which the water is out of the proton transfer loop. The corresponding forward/backward barriers are 65.89/65.89 kcal · mol^?1. Clearly, this process is more difficult to follow than the R–T–P process. For the two‐water‐assisting mechanism, the heat of reaction is 19.61 kcal · mol^?1, and the forward/backward barriers are 32.27/12.66 kcal · mol^?1, decreased by 33.25/31.21 kcal · mol^?1 compared with T. For trihydration and tetrahydration, the forward/backward barriers decrease as 32.00/12.60 (30T) and 37.38/17.26 (40T) kcal · mol^?1, and the heat of reaction decreases by 19.39 and 19.23 kcal · mol^?1, compared with T, respectively. But, when four water molecules are involved in the reactant loop, the corresponding energy aspects increase compared with those of the trihydration. The forward/backward barriers are increased by 5.38 and 4.66 kcal · mol^?1 than the trihydration situation. In the second mode, the outer‐sphere water effect from the other water molecules directly H‐bonded to the loop is considered. When one to three water molecules attach to the looped water in one‐water in‐loop‐assisting proton transfer isomerization, their effects on the three energies are small, and the deviations are not more than 3 kcal · mol^?1 compared with the original monohydration‐assisting case. When adding one or two water molecules on the dihydration‐assisting mechanism, and increasing one water molecule on the trihydration, the corresponding energies also are not obviously changed. The results indicate that the forward/backward barriers for the three in‐loop water‐assisting case are the lowest, and the surrounding water molecules (out‐of‐loop) yield only a small effect. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Interaction between amphiphiles and water molecules in concentrated bilayer aqueous colloids

Interaction between amphiphiles and water molecules in micelle or bilayer structure has been investigated using aqueous colloids of various amphiphiles through the rheological data and the spin-lattice relaxation timeT ₁ of the proton of water molecule.T ₁ of the water proton has been measured by the inversion recovery method and determined as a single exponential relaxation process.The chemical shift of the water proton is almost independent of the amphiphilic concentration; however, it shifts toward a higher magnetic field with increasing temperature in a way similar to that in pure water and in the amphiphilic aqueous systems. These facts mean that there is no significant difference in the magnetic field environment of the water protons in these systems.The water molecule is not necessarily bound in the fully developed micelle or bilayer (rod-like or lamella) structure which induces the high viscosity or high rigidity of the colloidal system. On the other hand, the water molecule is bound in the micelle colloids of amphoteric amphiphiles or amphiphiles whose molecular assembly creates a relatively strong electrostatic field. The activation entropy of the bound water is negative and this suggests that water molecules assume some ordered structure in the bound state.     

A first principles molecular dynamics study of the solvation structure and migration kinetics of an excess proton and a hydroxide ion in binary water-ammonia mixtures

We have investigated the solvation structure and migration kinetics of an excess proton and a hydroxide ion in water-ammonia mixed liquids of varying composition by means of ab initio molecular dynamics simulations. The excess proton is always found to be attached to an ammonia molecule to form the ammonium ion. Migration of the excess proton is found to occur very occasionally from one ammonia to the other but no proton transfer to a water molecule is observed during the entire simulations. Also, when the ammonium ion is solvated in water only, its hydrogen bond dynamics and rotation are found to occur at a faster rate than those in water-ammonia mixtures. For water-ammonia mixtures containing a proton less, the defect is found to stay like the hydroxide ion. For these systems, occasional proton transfer is found to occur only through the hydrogen bonded chains of water molecules in these water-ammonia mixtures. No proton transfer is found to take place from an ammonia molecule. The presence of ammonia molecules makes the realization of proper presolvated state of the hydroxide ion to accept a proton a more difficult process and, as a result, the rate of proton transfer and migration kinetics of the hydroxide ion in water-ammonia mixtures are found to be slower than that in liquid water and these rates are found to slow down further with increase of ammonia concentration.     

Computational Study of Proton Transfer in Tautomers of 3‐ and 5‐Hydroxypyrazole Assisted by Water

The tautomerism of 3‐ and 5‐hydroxypyrazole is studied at the B3LYP, CCSD and G3B3 computational levels, including the gas phase, PCM–water effects, and proton transfer assisted by water molecules. To understand the propensity of tautomerization, hydrogen‐bond acidity and basicity of neutral species is approached by means of correlations between donor/acceptor ability and H‐bond interaction energies. Tautomerism processes are highly dependent on the solvent environment, and a significant reduction of the transition barriers upon solvation is seen. In addition, the inclusion of a single water molecule to assist proton transfer decreases the barriers between tautomers. Although the second water molecule further reduces those barriers, its effect is less appreciable than the first one. Neutral species present more stable minima than anionic and cationic species, but relatively similar transition barriers to anionic tautomers.     

Perfluorobutane sulfonic acid hydration and interactions with O2 adsorbed on Pt3

The side chain of NAFION, a proton conductive membrane used as electrolyte in low-temperature fuel cells, is modeled with perfluorobutane sulfonic acid. Density functional theory is used to characterize structures and energetics of hydration of the model system interacting with a proton solvated with up to 24 water molecules and analyze interactions of some of these hydrated complexes with O(2) adsorbed on Pt(3). It is found that at least three water molecules are needed to ionize the sulfonic acid, and higher degrees of hydration induce the formation of cages where the water molecules are held together via complex hydrogen-bond networks. The interaction between the complex formed by the ionized acid and the hydrated proton, in contact with a bridge-adsorbed O(2)-Pt(3), promotes the protonation of the adsorbed O(2). Upon protonation, the O(2)-Pt(3) system evolves from hydrophobic to hydrophilic behavior, which may facilitate further interfacial contact.     

Long‐Range Electrostatics‐Induced Two‐Proton Transfer Captured by Neutron Crystallography in an Enzyme Catalytic Site

Neutron crystallography was used to directly locate two protons before and after a pH‐induced two‐proton transfer between catalytic aspartic acid residues and the hydroxy group of the bound clinical drug darunavir, located in the catalytic site of enzyme HIV‐1 protease. The two‐proton transfer is triggered by electrostatic effects arising from protonation state changes of surface residues far from the active site. The mechanism and pH effect are supported by quantum mechanics/molecular mechanics (QM/MM) calculations. The low‐pH proton configuration in the catalytic site is deemed critical for the catalytic action of this enzyme and may apply more generally to other aspartic proteases. Neutrons therefore represent a superb probe to obtain structural details for proton transfer reactions in biological systems at a truly atomic level.     

Theoretical study on the mechanism of the thermal decarboxylation of acrylic and benzoic acids. Models for aqueous solution

The mechanism of acrylic and benzoic acid decarboxylation in aqueous solution has been investigated by ab initio methods using the STO-3G and 3-21G basis sets. In those reactions, the solvent is represented successively by one and two water molecules. Their active participation as a proton relay in the chemical process is demonstrated by the large decrease in the activation energy with respect to the reaction studied in the absence of water. In the absence of any intermediate found along the reaction pathway, the proposed mechanism is the concerted process; the free acid being the species that undergoes decarboxylation via a pseudounimolecular mechanism by interaction with a chain of water molecules. At the transition state, the carboxylic hydrogen transfer to one water molecule, the reorganization of the chain of water molecules through which the proton is transferred and the cleavage of the C? C bond are much more advanced than the proton transfer from the last water molecule to the α-carbon atom of the carboxyl group.     

On the quantum nature of an excess proton in liquid hydrogen fluoride.

Liquid hydrogen fluoride consists of chains of hydrogen-bonded molecules. The nature of an excess proton in liquid HF, which is of interest not only for its own sake, but also in relation to super-acid chemistry and to its behavior in water, has been studied using computer simulations. The methodology employed is the density-functional-theory-based path-integral Car-Parrinello ab initio molecular dynamics. The excess proton, which formally exists as a H2F+ or a H2F2+ defect in an HF chain, is found to strongly perturb the chain to which it is attached. Moreover, due to large zero-point energy, the charge defect is largely delocalized over several HF molecules.     

Crystal and molecular structure of N-(benzylimidazolyl-2)-O-methylcarbamate salts

Single crystal X-ray diffraction has been applied to determine the structure of salts — formate and hydrochloride of N-(benzylimidazolyl-2)-O-methylcarbamate (BMC). The crystal structure of BMC formate is built by a molecule of a base and two formic acid molecules, one of them protonating a BMC molecule. Hydrogen bonds in the crystal form a weakly bound one-dimensional ribbon. BMC hydrochloride crystallizes as dihydrate. Two molecules of crystallization water and Cl ion make a robust H-bonded two-dimensional layer. BMC salts are formed through the protonation of N9 atom.     

Thermally Induced Intra‐Carboxyl Proton Shuttle in a Molecular Rack‐and‐Pinion Cascade Achieving Macroscopic Crystal Deformation

Proton transport via dynamic molecules is ubiquitous in chemistry and biology. However, its use as a switching mechanism for properties in functional molecular assemblies is far less common. In this study, we demonstrate how an intra‐carboxyl proton shuttle can be generated in a molecular assembly akin to a rack‐and‐pinion cascade via a thermally induced single‐crystal‐to‐single‐crystal phase transition. In a triply interpenetrated supramolecular organic framework (SOF), a 4,4′‐azopyridine (azpy) molecule connects to two biphenyl‐3,3′,5,5′‐tetracarboxylic acid (H₄BPTC) molecules to form a functional molecular system with switchable mechanical properties. A temperature change reversibly triggers a molecular movement akin to a rack‐and‐pinion cascade, which mainly involves 1) an intra‐carboxyl proton shuttle coupled with tilting of the azo molecules and azo pedal motion and 2) H₄BPTC translation. Moreover, both the molecular motions are collective, and being propagated across the entire framework, leading to a macroscopic crystal expansion and contraction.     

Proton dissociation and transfer in hydrated phosphoric acid clusters

The dynamics and mechanisms of proton dissociation and transfer in hydrated phosphoric acid (H₃PO₄) clusters under excess proton conditions were studied based on the concept of presolvation using the H₃PO₄–H₃O⁺–nH₂O complexes (n = 1–3) as the model systems and ab initio calculations and Born–Oppenheimer molecular dynamics (BOMD) simulations at the RIMP2/TZVP level as model calculations. The static results showed that the smallest, most stable intermediate complex for proton dissociation (n = 1) is formed in a low local‐dielectric constant environment (e.g., ε = 1), whereas proton transfer from the first to the second hydration shell is driven by fluctuations in the number of water molecules in a high local‐dielectric constant environment (e.g., ε = 78) through the Zundel complex in a linear H‐bond chain (n = 3). The two‐dimensional potential energy surfaces (2D‐PES) of the intermediate complex (n = 1) suggested three characteristic vibrational and ¹H NMR frequencies associated with a proton moving on the oscillatory shuttling and structural diffusion paths, which can be used to monitor the dynamics of proton dissociation in the H‐bond clusters. The BOMD simulations over the temperature range of 298–430 K validated the proposed proton dissociation and transfer mechanisms by showing that good agreement between the theoretical and experimental data can be achieved with the proposed rate‐determining processes. The theoretical results suggest the roles played by the polar solvent and iterate that insights into the dynamics and mechanisms of proton transfer in the protonated H‐bond clusters can be obtained from intermediate complexes provided that an appropriate presolvation model is selected and that all of the important rate‐determining processes are included in the model calculations. © 2015 Wiley Periodicals, Inc.     

Energetics during proton transfer process in hydrated zündel ion complex

Zündel ion (H₅O) is one of the two important structures formed during the proton transfer process in aqueous system. This work reports microsolvation of Zündel ion using density functional theory based B3LYP method with aug‐cc‐pVTZ basis set. Interaction of Zündel ion with four water molecules in its first solvation shell is studied using many‐body analysis approach. A change in many‐body energies and their contribution to the binding energy of a complex during the proton transfer process from donor to acceptor water molecule in Zündel ion‐4H₂O complex is obtained. For the hydrated Zündel ion complex, the contribution from total two‐body, three‐body, four‐body, five‐body, and relaxation energy to the binding energy is 84.7, 14, 6.87, 1.6, and 4%, respectively, at B3LYP/aug‐cc‐pVTZ level. Relaxation energy and total five‐body energy have repulsive contribution to the binding energy of a hydrated Zündel ion complex. It is found that the relaxation energy and binding energy of a Zündel‐4H₂O complex is the maximum and minimum, respectively, when a shared proton is at equal distance from oxygen atom of donor and acceptor water molecules. A significant change in two‐body, three‐body, and four‐body energies for which Zündel ion is one of the many‐body terms is observed during the proton transfer process. A change in total two‐body, total three‐body, total four‐body, and relaxation energy is about 2.6, 1.8, 0.4, and 1.1%, respectively, during the proton transfer process. A change in two‐body, three‐body, and four‐body interaction energies between water molecules is very small during the proton transfer process. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012