Insights on the mechanism of proton transfer reactions in amino acids

In the present work, the joint use of the potential energy, the reaction electronic flux profiles and NBO analysis along the intrinsic reaction coordinate within the framework of the reaction force analysis allows us to gain insights into the mechanism of the proton transfer process in amino acids. The reaction was studied in alanine and phenylalanine in the presence of a continuum and with addition of one water molecule acting as a bridge, the results were compared to those of tryptophan. The bridging water molecule stabilizes the zwitterionic form and increases the reaction barriers by a factor of two. This result is interpreted in terms of the energy required to bring the amino acid and the water molecule closer to each other and to promote the proton transfer through the reordering of the electron density. Furthermore, the bridging water molecule induces a concerted asynchronous double proton transfer, where the transfer of the carboxyl hydrogen atom is followed by the second proton transfer to the ammonium group. In addition, a second not intervening water molecule was added, which changes the proton acceptor and donor properties of the reactive water molecule modulating the reaction mechanism. The aforementioned methods allow us to identify the order of the transferred protons and the asynchronicity, thereby, evolving as promising tools to not only characterize but also manipulate reaction mechanisms.     

金催化炔基苯并二𫫇英环化合成8-羟基异香豆素的理论研究

Density functional theory（DFT） calculations were carried out on the gold-catalyzed cyclization of alkynyl benzodioxin to 8-hydroxy-isocoumarin reaction to show the molecular mechanism of the reaction. The conclusions obtained from this work are different from those in the previous experimental study. The results show that water molecule acts as both the reactant and the proton shuttle, and promotes the reaction with gold complexes under mild conditions. The nucleophilic addition site of water on the substrate is the C（sp³） atom on the side of the substrate far away from the oxabenzene ring, resulting in C（sp³）—O bond breaking in the substrate. The formation of new C—O bond and the cleavage of C—O bond in the substrate follow a step-by-step mechanism. The oxygen in the side-product acetone comes from the contribution of water in the reaction system. The regioselectivity of the reaction originates from the polarization of alkynyl π-electrons induced by substituents.     

C-H⋯O, O-H⋯C, and C-H⋯C interactions in complexes of carbocations and carboanions

Molecular complexes formed by different forms of carbocations (carbenium ions) and carboanions with water, acetylene, and methane molecules have been calculated by the MP2/6-311++G(2df,2pd) method. In complexes with water where the carbon atom of the carbocation (carboanion) acts as the proton donor (acceptor), the energies of the C-H?O and O-H?C hydrogen bonds turn out to be approximately the same being 13–20 kcal/mol for carbocation (carboanion) species differing in the valence state of the carbon atom. Two types of C-H?C interactions have been revealed depending on the charge at the bridging hydrogen atom, which is determined by the hybridization of the donor carbon atom. The C-H?C interaction energy in molecular complexes with the positively charged hydrogen atom (carboanion complexes with acetylene) is an order of magnitude higher than in the complexes where the bridging hydrogen atom has an excess of electron density (carbocation complexes with methane). In all the complexes under consideration, the covalent C-H bond involved in interaction is elongated, and the negative charge is transferred from the acceptor to the donor.     

Charge transport in poly-imidazole membranes: a fresh appraisal of the Grotthuss mechanism

A detailed theoretical investigation of the charge transport mechanism in poly(4-vinyl-imidazole) (P4VI), the parent polymer of a series of N-heterocyclic-based membranes used as an electrolyte in proton exchange membrane fuel cells, is presented. In particular, Density Functional Theory (DFT) results obtained for small model systems (protonated imidazole dimers and trimers) suggest that the commonly accepted conduction mechanism, based on a sequential proton transfer between imidazole moieties, could be impeded by the geometrical constraints imposed by the polymeric backbone. Indeed only one kind of proton transfer reaction is energetically allowed between adjacent imidazoles, so that a rotation of the protonated imidazole is required for a second proton transfer. Molecular dynamics simulations on a larger model (15 oligomers with an excess proton) show that the rotation of the imidazole carrying the excess proton is a soft large amplitude motion. These results allow us to propose a new proton conduction mechanism in P4VI, where a frustrated rotation of the protonated imidazole before each proton transfer reaction represents the rate-limiting step. Furthermore, in contrast with the Grotthuss proton transport mechanism in water, our results indicate that here it is the same proton which could be successively transferred. From a chemical point of view, these new insights into the mechanism are relevant for a rational design of modified azole-based systems for Proton Exchange Membrane Fuel Cells.     

DNA charge transport leading to disulfide bond formation

Here, we show that DNA-mediated charge transport (CT) can lead to the oxidation of thiols to form disulfide bonds in DNA. DNA assemblies were prepared possessing anthraquinone (AQ) as a photooxidant spatially separated on the duplex from two SH groups incorporated into the DNA backbone. Upon AQ irradiation, HPLC analysis reveals DNA ligated through a disulfide. The reaction efficiency is seen to vary in assemblies containing intervening DNA mismatches, confirming that the reaction is DNA-mediated. Interestingly, one intervening mismatch near the thiols promotes an increase in efficiency, which we attribute to increased base dynamics. Hence, here, where the reaction is on the backbone rather than within the base stack, stacking perturbations do not necessarily lead to an inhibitory effect on DNA CT.     

Properties of phosphorothioate DNA analogs. An ab initio study of prototype model linkages derived from dimethyl-phosphate anion

Ab initio HF and MP2 calculations on prototype model linkages of phosphorothioate DNA backbones illuminate the effects of phosphorothioation on electronic and structural properties of DNA backbone. The replacement of a bridging oxygen atom by sulfur in the phosphodiester linkage is energetically favored over that of replacement of a non-bridging oxygen atom. In phosphorothioate derivatives containing the P(OS)_nb moiety, the non-bridging oxygen atom always bears a higher negative charge than the non-bridging sulfur. Additional calculations on protonated (neutral) adducts suggest that phosphorothioation of the phosphodiester linkage lowers its proton affinity. Moreover, protonation of the non-bridging oxygen atom at phosphorous is favored over the protonation of the non-bridging sulfur atom for linkages containing the P(OS)_nb moiety. The ab initio calculated structural parameters are compared to the available crystallographic data of small phosphorothioate molecules and phosphorothioate oligodeoxynucleotides. These results have implications upon the biological activity of phosphorothioate DNA analogs.     

A theoretical study on the water-mediated asynchronous addition between urea and formaldehyde

The reaction between urea and formaldehyde in water solution was theoretically investigated by using B3LYP and MP2 methods.It was found that the addition of the nitrogen atom in urea to the carbonyl group in formaldehyde precedes the proton transfer and the proton migration from water to the carbonyl group occurs before the proton abstraction from the nitrogen.With one or two water molecules involved in the TS.the activation energy barrier is lowered compared to the TS of the mechanism with no water participation.The energy change along the reaction coordinate clearly shows that a zwitterionic-like intermediate does not exist on the PES.The reaction between urea and formaldehyde occurs in a concerted mechanism but with asynchronous characters.This is different from the stepwise mechanism recently found for the amination reactions of formaldehyde.     

Computer simulation of the chemical catalysis of DNA polymerases: discriminating between alternative nucleotide insertion mechanisms for T7 DNA polymerase

Understanding the chemical step in the catalytic reaction of DNA polymerases is essential for elucidating the molecular basis of the fidelity of DNA replication. The present work evaluates the free energy surface for the nucleotide transfer reaction of T7 polymerase by free energy perturbation/empirical valence bond (FEP/EVB) calculations. A key aspect of the enzyme simulation is a comparison of enzymatic free energy profiles with the corresponding reference reactions in water using the same computational methodology, thereby enabling a quantitative estimate for the free energy of the nucleotide insertion reaction. The reaction is driven by the FEP/EVB methodology between valence bond structures representing the reactant, pentacovalent intermediate, and the product states. This pathway corresponds to three microscopic chemical steps, deprotonation of the attacking group, a nucleophilic attack on the P(alpha) atom of the dNTP substrate, and departure of the leaving group. Three different mechanisms for the first microscopic step, the generation of the RO(-) nucleophile from the 3'-OH hydroxyl of the primer, are examined: (i) proton transfer to the bulk solvent, (ii) proton transfer to one of the ionic oxygens of the P(alpha) phosphate group, and (iii) proton transfer to the ionized Asp654 residue. The most favorable reaction mechanism in T7 pol is predicted to involve the proton transfer to Asp654. This finding sheds light on the long standing issue of the actual role of conserved aspartates. The structural preorganization that helps to catalyze the reaction is also considered and analyzed. The overall calculated mechanism consists of three subsequent steps with a similar activation free energy of about 12 kcal/mol. The similarity of the activation barriers of the three microscopic chemical steps indicates that the T7 polymerase may select against the incorrect dNTP substrate by raising any of these barriers. The relative height of these barriers comparing right and wrong dNTP substrates should therefore be a primary focus of future computational studies of the fidelity of DNA polymerases.     

Intrinsic proton-donating power of zinc-bound water in a carbonic anhydrase active site model estimated by NMR

Using liquid-state NMR spectroscopy we have estimated the proton-donating ability of Zn-bound water in organometallic complexes designed as models for the active site of the metalloenzyme carbonic anhydrase (CA). This ability is important for the understanding of the enzyme reaction mechanism. The desired information was obtained by (1)H and (15)N NMR at 180 K of solutions of [Tp(Ph,Me)ZnOH] [1, Tp(Ph,Me) = tris(2-methyl-4-phenylpyrazolyl)hydroborate] in CD(2)Cl(2), in the absence and presence of the proton donors (C(6)F(5))(3)BOH(2) [aquatris(pentafluorophenyl)boron] and Col-H(+) (2,4,6-trimethylpyridine-H(+)). Col-H(+) forms a strong OHN hydrogen bond with 1, where the proton is located closer to nitrogen than to oxygen. (C(6)F(5))(3)BOH(2), which exhibits a pK(a) value of 1 in water, also forms a strong hydrogen bond with 1, where the proton is shifted slightly across the hydrogen-bond center toward the Zn-bound oxygen. Finally, a complex between Col and (C(6)F(5))(3)BOH(2) was identified, exhibiting a zwitterionic OHN hydrogen bond, where H is entirely shifted to nitrogen. The comparison with complexes of Col with carboxylic acids studied previously suggests that, surprisingly, the Zn-bound water exhibits in an aprotic environment a similar proton-donating ability as a carboxylic acid characterized in water by a pK(a) of 2.2 ± 0.6. This value is much smaller than the value of 9 found for [Zn(OH(2))(6)](2+) in water and those between 5 and 8 reported for different forms of CA. Implications for the biological function of CA are discussed.     

Synthesis and Characterization of Sulfonated Polyphosphazene-graft-Polystyrene Copolymersfor Proton Exchange Membranes简

A novel series of polyphosphazene-grafl-polystyrene （PP-g-PS） copolymers were successfully prepared by atom transfer radical polymerization （ATRP） of styrene monomers and brominated poly（bis（4-methylphenoxy）phosphazene） macroinitiator. The graft density and the graft length could be regulated by changing the bromination degree of the macroinitiator and the ATRP reaction time, respectively. The PP-g-PS copolymers readily underwent a regioselective sulfonation reaction, which occurred preferentially at the polystyrene sites, producing the sulfonated PP-g-PS copolymers with a range of ion exchange capacities. The resulting sulfonated PP-g-PS membranes prepared by solution casting showed high water uptake, low water swelling and considerable proton conductivity. They also exhibited good oxidative stability and high resistance to methanol crossover. Morphological studies of the membranes by transmission electron microscopy showed clear nanophase-separated structures resulted from hydrophobic polyphosphazene backbone and hydrophilic polystyrene sulfonic acid segments, indicating the formation of proton transferring tunnels. Therefore, these sulfonated copolymers may be candidate materials for proton exchange membranes in direct methanol fuel cell （DMFC） applications.     

The N6-methyl group of adenine further increases the BI stability of DNA compared to C5-methyl groups

Methylated DNA bases are natural modifications which play an important role in protein-DNA interactions. Recent experimental and theoretical results have shown an influence of the base modification on the conformational behavior of the DNA backbone. MD simulations of four different B-DNA dodecamers (d(GC)(6), d(AT)(6), d(G(5mCG)(5)C), and d(A(T6mA)(5)T)) have been performed with the aim to examine the influence of methyl groups on the B-DNA backbone behavior. An additional control simulation of d(AU)(6) has also been performed to examine the further influence of the C5-methyl group in thymine. Methyl groups in the major groove (as in C5-methylcytosine, thymine, or N6-methyladenine) decrease the BII substate population of RpY steps. Due to methylation a clearer distinction of the BI substate stability between YpR and RpY (CpG/GpC or TpA/ApT) steps arises. A positive correlation between the BII substate population and base stacking distances is seen only for poly(GC). A methyl group added into the major groove increases mean water residence times around the purine N7 atom, which may stabilize the BI substate by improving the hydration network between the DNA backbone and the major groove. The N6-methyl group also forms a water molecule bridge between the N6 and O4 atoms, and thus further stabilizes the BI substate.     

A theoretical study of the mechanism of the nucleotidyl transfer reaction catalyzed by yeast RNA polymerase II

The mechanism of the nucleotidyl transfer reaction catalyzed by yeast RNA polymerase II has been investigated using molecular mechanics and quantum mechanics methods.Molecular dynamics(MD) simulations were carried out using the TIP3 water model and generalized solvent boundary potential(GSBP) by CHARMM based on the X-ray crystal structure.Two models of the ternary elongation complex were constructed based on CHARMM MD calculations.All the species including reactants,transition states,intermediates,and products were optimized using the DFT-PBE method coupled with the basis set DZVP and the auxiliary basis set GEN-A2.Three pathways were explored using the DFT method.The most favorable reaction pathway involves indirect proton migration from the RNA primer 3’-OH to the oxygen atom of-phosphate via a solvent water molecule,proton rotation from the oxygen atom of-phosphate to the-phosphate side,the RNA primer 3’-O nucleophilic attack on the-phosphorus atom,and P-O bond breakage.The corresponding reaction potential profile was obtained.The rate limiting step,with a barrier height of 21.5 kcal/mol,is the RNA primer 3’-O nucleophilic attack,rather than the commonly considered proton transfer process.A high-resolution crystal structure including crystallographic water molecules is required for further studies.     

Protonation of the proximal histidine ligand in heme peroxidases

The heme peroxidases have a histidine group as the axial ligand of iron. This ligand forms a hydrogen bond to an aspartate carboxylate group by the other nitrogen atom in the side chain. The aspartate is not present in the globins and it has been suggested that it gives an imidazolate character to the histidine ligand. Quantum chemical calculations have indicated that the properties of the heme site strongly depend on the position of the proton in this hydrogen bond. Therefore, we have studied the location of this proton in all intermediates in the reaction mechanism, using a set of different quantum mechanical and combined experimental and computational methods. Quantum refinements of a crystal structure of the resting FeIII state in yeast cytochrome c peroxidase show that the geometric differences of the two states are so small that it cannot be unambiguously decided where the proton is in the crystal structure. Vacuum calculations indicate that the position of the proton is sensitive to the surroundings and to the side chains of the porphyrin ring. Combined quantum and molecular mechanics (QM/MM) calculations indicate that the proton prefers to reside on the His ligand in all states in the reaction mechanism of the peroxidases. QM/MM free energy perturbations confirm these results, but reduce the energy difference between the two states to 12-44 kJ/mol.     

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.     

Theoretical study of the acid-promoted hydrolysis of oxazolin-5-one: a microhydration model

The acid-promoted hydrolysis of 2,4,4-trimethyloxazolin-5-one (TMO) is studied employing the density functional theory (B3LYP) method in conjunction with the 6-31++G(d,p) basis set. Two types of reaction mechanism, N-protonated and O-protonated, are considered, involving protonation at the nitrogen and carbonyl oxygen of TMO to activate the C2 and C5 atoms, respectively, in favor of attack by water molecules. In the N-protonated pathway, the nucleophilic water molecule attacks the activated C2 atom, with a proton transfer from the water molecule to the oxygen atom attached to C2 and the fission of the C2-O bond, leading to a cis ring-opening product (N-acyl-alpha-amino isobutyric acid). While, in the O-protonated pathway, the nucleophilic water molecule attacks the activated carbonyl C5 atom, accompanied by a proton transfer from the water molecule toward the nitrogen atom of oxazole ring and the cleavage of C5-O bond; as a result, a corresponding trans product is generated. The water-assisted hydrolysis reactions are also examined together. A local microhydration model, in which an extra water molecule was added to obtain a continuous H-bond network around the reaction centers, was adopted to mimic the system for the two types of reaction processes. In addition, bulk solvent effect is introduced by use of the conductor-like polarizable continuum model (CPCM). Our computational results in kinetics and thermodynamics clearly manifest that the O-protonated pathway with the nucleophilic attack at the carbonyl C5 atom is more favorable than the N-protonated one, in nice agreement with the available experimental conclusion.     

Hydration effects on the reaction with an open-shell transition state: QM/MM-ER study for the dehydration reaction of alcohol in hot water

The reaction mechanism for the dehydration of 1,4-butanediol in hot water has been investigated by means of the hybrid quantum mechanical/molecular mechanical approach combined with the theory of energy representation (QM/MM-ER). We have assumed that the proton transfers along the hydrogen bonds of the water molecules catalyze the reaction, where the transition state (TS) forms a singlet biradical electronic structure. It has been revealed by the simulation that the biradical electronic state at the TS changes to zwitterionic structure in solution due to the hydration of the polar solvent. Such the electronic structure change gives rise to the substantial stabilization of the TS in hot water. As a result, the water-catalytic path becomes more favorable in aqueous solution than another possible path that proceeds without proton transfers as opposed to the reaction mechanism in the gas phase. Furthermore, the activation free energy computed by the present method is in excellent agreement with the experimental result.     

Deprotonation of a histidine residue in aqueous solution using constrained ab initio molecular dynamics

The dissociation of a weak acid - a histidine residue - in water was investigated by means of constrained Car-Parrinello ab initio molecular dynamics. Both linear and coordination constraints were employed, and the structural, electronic, and dynamical transformations along the respective reaction coordinates were analyzed. The calculated potentials of mean force for the dissociation of a hydrogen atom from the Nepsilon and Ndelta positions of the imidazole ring reveal that protonated forms are approximately 9.0-9.5 kcal/mol more stable than the deprotonated. This result seems to agree well with the experimental estimate based on pKa. A possible transition state for the deprotonation has also been identified. Analysis of the electron localization function indicates that the proton transfer along the selected reaction path is not a fully concerted process.     

用密度泛函方法研究质子化苯基丙酮-水团簇

The structure and growth trend of the protonated acetophenone-water clusters have been investigated using the DFT-B3LYP method combined with the standard 6-31+G(d,p) basis set. In order to obtain more accurate single-point energy the B3LYP/6-311++G(3df,2p) method was adapted. The results show that the formation of H+C8H8O-H2O is a barrierless reaction process and the equilibrium distance between the proton and the O atom in C8H8O molecule is 1.015 A. For H+C8H8O-(H2O)n(n=1,2,3) clusters, the proton lies between the acetophenone molecule C8H8O and the water molecule H2O. The distance between the proton and the O atom of the C8H8O molecule increased from n=1 to n=3; C8H8O-H+-H2O can be regarded as an solvation shell. For H+C8H8O (H2O)n (n=4,5,6,7,8) clusters, the proton lies between the two H2O molecules forming a H5O2+ structure, C8H8O-H5O2+ is an important structure, which the other H2O molecules will attack from different sides.     

Catalytic mechanism of RNA backbone cleavage by ribonuclease H from quantum mechanics/molecular mechanics simulations

We use quantum mechanics/molecular mechanics simulations to study the cleavage of the ribonucleic acid (RNA) backbone catalyzed by ribonuclease H. This protein is a prototypical member of a large family of enzymes that use two-metal catalysis to process nucleic acids. By combining Hamiltonian replica exchange with a finite-temperature string method, we calculate the free energy surface underlying the RNA-cleavage reaction and characterize its mechanism. We find that the reaction proceeds in two steps. In a first step, catalyzed primarily by magnesium ion A and its ligands, a water molecule attacks the scissile phosphate. Consistent with thiol-substitution experiments, a water proton is transferred to the downstream phosphate group. The transient phosphorane formed as a result of this nucleophilic attack decays by breaking the bond between the phosphate and the ribose oxygen. In the resulting intermediate, the dissociated but unprotonated leaving group forms an alkoxide coordinated to magnesium ion B. In a second step, the reaction is completed by protonation of the leaving group, with a neutral Asp132 as a likely proton donor. The overall reaction barrier of ～15 kcal mol(-1), encountered in the first step, together with the cost of protonating Asp132, is consistent with the slow measured rate of ～1-100/min. The two-step mechanism is also consistent with the bell-shaped pH dependence of the reaction rate. The nonmonotonic relative motion of the magnesium ions along the reaction pathway agrees with X-ray crystal structures. Proton-transfer reactions and changes in the metal ion coordination emerge as central factors in the RNA-cleavage reaction.