Potential energy function for cation–peptide interactions: An ab initio study

Abstract

A potential energy function is developed to represent the interaction of small monovalent cations, Li', Na', and K' , with the backbone of polypeptides. The results are based on ab initio calculations up to the 6-31G* level of the interactions of the ions with acetamide and N-methylacetamide. Basis set superposition errors are corrected with the counterpoise method. A systematic overestimate of the bond polarities is taken into account by an empirical scaling procedure that uses the ratio of the experimental to ab initio dipole moment. The calculated binding energies obtained with this procedure show consistent convergence with different basis sets and are in good agreement with experimental data on cation-water and cation-dimethylformamide systems. Investigations of the calculated ab initio potential energy surface indicate that the cation-peptide interaction is dominated by electrostatics and includes a nonnegligible contribution from polarization of the peptide group by the ion. The induced polarization results in a steeper-than-Coulombic interaction and cannot be described by fixed ion-peptide partial charges electrostatics. Atomic polarizabilities located on the atoms of the ligand molecule are introduced to account for the induced polarization in the empirical energy function. A -l/r4 attractive interaction appears in the potential function. The resulting radial and angular dependence of the potential energy surface is well reproduced. 0 1995

GPU-Accelerated Molecular Dynamics and Free Energy Methods in Amber18: Performance Enhancements and New Features

“…30 The 12-6-4 potentials are derived from regular 12–6 potentials by adding r − 4 terms, which has been proposed by Roux and Karplus. 57,58 …”

Dynamical behavior of one-dimensional water molecule chains in zeolites: Nanosecond time-scale molecular dynamics simulations of bikitaite

“…The special features of Li ϩ -water interactions were recently discussed by Lyubartsev et al 29 in an ab initio MD simulation study of lithium salt aqueous solutions. Therefore, a new model fitting both ab initio results [30][31][32][33][34] for the Li ϩ -water system and experimental data for bikitaite was performed. The final form of the Li ϩ -water potential function reads…”

“…The values of the parameters are reported in Table I, and the optimized quantities of the isolated Li ϩ -H 2 O system are compared with the ab initio results [30][31][32][33][34] in Table II. As in our previous papers, [12][13][14] water was assumed to interact with Si and Al atoms via a Coulomb potential only, and the potential functions between the oxygen atoms of the framework and the oxygen or hydrogen atoms of water were derived from a simplified form of the corresponding O-O and O-H ones for water-water interactions.…”

A QM–MM interface between CHARMM and TURBOMOLE: Implementation and application to systems in bulk phase and biologically active systems

“…2) was partitioned in two ways: in one, the QM water is hydrogen bond donating ("Donor" in Table 1); in the other, it is hydrogen bond accepting ("Acceptor" in Table 1). For the MM water the TIP3P parameters 54,55 were used that reflect the average amount of polarization experienced in water. Additionally, the whole water dimer was calculated quantum mechanically.…”