Showing papers by "Vincent Zoete published in 2004"

Double proton transfer in the isolated and DNA-embedded guanine-cytosine base pair

Abstract: The energetics and dynamics of double proton transfer (DPT) is investigated theoretically for the Watson-Crick conformation of the guanine-cytosine (GC) base pair. Using semiempirical density functional theory the isolated and DNA-embedded GC pair is considered. Differences in the energetics and dynamics of DPT thus addresses the question of how relevant studies of isolated base pairs are for the understanding of processes occurring in DNA. Two-dimensional potential energy surfaces involving the transferring hydrogen atoms and the proton donors and acceptors are presented for both systems. The DPT reaction is accompanied by a contraction of the distance between the two bases with virtually identical energetic barriers being 18.8 and 18.7 kcal/mol for the isolated and DNA-embedded system, respectively. However, the transition state for DPT in the DNA-embedded GC pair is offset by 0.1 A to larger N-H separation compared to the isolated GC pair. Using activated ab initio molecular dynamics, DPT is readily observed for the isolated base pair with a minimal amount of 21.4 kcal/mol of initial average kinetic energy along the DPT normal mode vector. On a time scale of approximately 100 fs DPT has occurred and the excess energy is redistributed. For the DNA-embedded GC pair considerably more kinetic energy is required (30.0 kcal/mol) for DPT and the process is completed within one hydrogen vibration. The relevance of studies of isolated base pairs and base pair analogs in regard of reactions or properties involving DNA is discussed.

Investigation of glucose binding sites on insulin

Abstract: Possible insulin binding sites for D-glucose have been investigated theoretically by docking and molecular dynamics (MD) simulations. Two different docking programs for small molecules were used; Multiple Copy Simultaneous Search (MCSS) and Solvation Energy for Exhaustive Docking (SEED) programs. The configurations resulting from the MCSS search were evaluated with a scoring function developed to estimate the binding free energy. SEED calculations were performed using various values for the dielectric constant of the solute. It is found that scores emphasizing non-polar interactions gave a preferential binding site in agreement with that inferred from recent fluorescence and NMR NOESY experiments. The calculated binding affinity of −1.4 to −3.5 kcal/mol is within the measured range of −2.0 ± 0.5 kcal/mol. The validity of the binding site is suggested by the dynamical stability of the bound glucose when examined with MD simulations with explicit solvent. Alternative binding sites were found in the simulations and their relative stabilities were estimated. The motions of the bound glucose during molecular dynamics simulations are correlated with the motions of the insulin side chains that are in contact with it and with larger scale insulin motions. These results raise the question of whether glucose binding to insulin could play a role in its activity. The results establish the complementarity of molecular dynamics simulations and normal mode analyses with the search for binding sites proposed with small molecule docking programs. Proteins 2004. © 2004 Wiley-Liss, Inc.

On the influence of semirigid environments on proton transfer along molecular chains.

Abstract: The dynamics of proton transfer along ammonia chains (chemical composition N(x)H(+)(3x+1), x=2, 4, and 6) in a constraining environment is investigated by ab initio molecular dynamics simulations. A carbon nanotube of defined length and diameter is used as an idealized constraining environment such that the ammonia chain is forced to maintain its quasilinear geometry. It is found that, although the energetics of proton transport shows considerable energetic barriers, proton translocation along the wire is possible at finite temperature for all chain lengths studied. The proton transport involves rotational reorientation of the proton-carrying ammonia molecule. High level ab initio calculations (MP2/aug-cc-pVTZ) yield barriers for internal rotation of 9.1 kcal/mol for NH(4) (+)-NH(3) and 11.7 kcal/mol for OH(3) (+)-OH(2), respectively. The infrared spectrum calculated from the dipole-dipole autocorrelation function shows distinct spectral features in the regions (2000-3000 cm(-1)) where the NHN proton transfer mode is expected to absorb. Assigning moderate opposite total charges between 0.002 and 0.2e to the carbon atoms at the end caps of the nanotube leads to a considerable speedup of the proton transfer.