Ab initio path-integral molecular dynamics

A path-integral molecular dynamics technique for strongly interacting atoms using ab initio potentials derived from density functional theory is implemented. This allows the efficient inclusion of nuclear quantum dispersion in ab initio simulations at finite temperatures. We present an application to the quantum cluster H ₅ ⁺ .     

Ab initio molecular dynamics-based assignment of the protonation state of pepstatin A/HIV-1 protease cleavage site

A recent 13C NMR experiment (Smith et al. Nature Struct. Biol. 1996, 3, 946-950) on the Asp 25-Asp25' dyad in pepstatin A/HIV-1 protease measured two separate resonance lines, which were interpreted as being a singly protonated dyad. We address this issue by performing ab initio molecular dynamics calculations on models for this site accompanied by calculations of 13C NMR chemical shifts and isotopic shifts. We find that already on the picosecond time-scale the model proposed by Smith et al. is not stable and evolves toward a different monoprotonated form whose NMR pattern differs from the experimental one. We suggest, instead, a different protonation state in which both aspartic groups are protonated. Despite the symmetric protonation state, the calculated 13C NMR properties are in good agreement with the experiment. We rationalize this result using a simple valence bond model, which explains the chemical inequality of the two C sites. The model calculations, together with our calculations on the complex, allow also the rationalization of 13C NMR properties on other HIV-1 PR/inhibitor complexes. Both putative binding of the substrate to the free enzyme, which has the dyad singly protonated (Piana, S.; Carloni, P. Proteins: Struct., Funct., Genet. 2000, 39, 26-36), and pepstatin A binding to the diprotonated form are consistent with the inverse solvent isotope effect on the onset of inhibition of pepsin by pepstatin and the kinetic iso-mechanism proposed for aspartic proteases (Cho, T.-K.; Rebholz, K.; Northrop, D.B. Biochemistry 1994, 33, 9637-9642).     

Nonperiodic boundary conditions for solvated systems

The simulation of charged and/or strongly polar solutes represents a challenge for standard molecular-dynamics techniques. The use of periodic boundary conditions (PBCs) leads to artifacts due to the interaction between two replicas in the presence of the long-range Coulomb forces. A way to avoid these problems is the use of nonperiodic boundary conditions. A possible realization is to consider a finite system, a sphere, embedded in a reaction field described by the method of the images. In the present work the modified image approximation has been implemented in a molecular-dynamics code and optimized for the use of two standard solvents, water and acetonitrile. The methodology has then been applied to investigate the conformational changes in water-solvated alanine dipeptide. The free-energy surface calculated with this method is comparable to that obtained with PBC.     

Grid-free DFT implementation of local and gradient-corrected XC functionals

Following an approach to density functional theory calculations based on the matrix representation of operators, we implemented a scheme as an alternative to traditional grid-based methods. These techniques allow integrals over exchange-correlation operators to be evaluated through matrix manipulations. Both local and gradient-corrected functionals can be treated in a similar way. After deriving all the required expressions, selected examples with various functionals are given. Received: 7 March 1998 / Accepted: 21 May 1998 / Published on line: 6 August 1998