Conformational dependence of charges in protein simulations

We have studied the conformational dependence of molecular mechanics atomic charges for proteins by calculating the charges fitted to the quantum mechanical (QM) electrostatic potential (ESP) for all atoms in complexes between avidin and seven biotin analogues for 20 snapshots from molecular dynamics simulations. We have studied how various other charge sets reproduce those charges. The QM charges, even if averaged over all snapshots or all residues, in general have a larger magnitude than standard Amber charges, indicating that the restraint toward zero in the restrained ESP method is too strong. This has a significant influence on the electrostatic conformational energies and the interaction energy between the biotin ligand and the protein, giving a difference between the QM and Amber charges of 43 and 8 kJ/mol for the negatively charged and neutral biotin analogues, respectively (3-4%). However, this energy difference is strongly reduced if the solvation energy (calculated by the Poisson-Boltzmann or Generalized Born methods) is added, viz., to 7 kJ/mol for charged and 3 kJ/mol for uncharged ligand. In fact, charges need to be recalculated with a QM method only for residues within 7 or 4 A of the ligand, if the error should be less than 4 kJ/mol. Unfortunately, the QM charges do not give significantly better MM/PBSA estimates of ligand-binding affinities than standard Amber charges.     

A practical procedure for the determination of electrostatic charges of large molecules

Summary A practical procedure for the precise determination of electrostatic charges, which are evaluated by fitting the rigorous quantum mechanical molecular electrostatic potential to a monopole-monopole expression, is presented. The proposal of this procedure arises from the study of the minimum requirements necessary to obtain reliable electrostatic charges. Such a study is focused on: (i) the dependence of the electrostatic charges on the set of points where the quantum mechanical and the monopole-monopole molecular electrostatic potentials are fitted; thus, both the influence of the number of points and their distribution in layers located out of the van der Waals radii of the atoms are examined, and (ii) the reliability of the use of fractional models for the evaluation of electrostatic charges of large molecules. Results point out that the optimum number of points is defined by a density of points ranging from 0.45 to 0.60 points/Ų when four layers (separated by 0.2 Å) are considered. Nevertheless, the use of only two layers (separated by 0.4 Å) for large molecules is recommended, thus enabling one to obtain reliable charges at a reduced computational cost. Moreover, results justify the use of fractional models for the determination of electrostatic charges of extremely large molecules, even when aromatic structures exist.     

Generation of OPLS-like charges from molecular electrostatic potential using restraints

A new method for generating atom-centered charges for use in condensed phase computer simulations is presented, which is based on a restrained electrostatic potential (RESP) procedure. Charges are calculated from a least-squares fit to the quantum mechanical electrostatic potential with a restraint applied to reduce their magnitude. The restraint developed here offers advantages over that used in RESP. The magnitude of the restraint is optimized to yield charges as close to the equivalent OPLS values as possible while still reproducing the molecule's electrostatic potential. A cross-validation analysis is used to show that the restraint is insensitive to the selection of OPLS molecules from which it is derived. Thus, with this method, OPLS-like charges may be produced from the electrostatic potential for atom types not in the OPLS force field. In addition, the restraint is shown to reduce the conformational dependence of the charges. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 483–498, 1999     

Reaction path potential for complex systems derived from combined ab initio quantum mechanical and molecular mechanical calculations

Combined ab initio quantum mechanical and molecular mechanical calculations have been widely used for modeling chemical reactions in complex systems such as enzymes, with most applications being based on the determination of a minimum energy path connecting the reactant through the transition state to the product in the enzyme environment. However, statistical mechanics sampling and reaction dynamics calculations with a combined ab initio quantum mechanical (QM) and molecular mechanical (MM) potential are still not feasible because of the computational costs associated mainly with the ab initio quantum mechanical calculations for the QM subsystem. To address this issue, a reaction path potential energy surface is developed here for statistical mechanics and dynamics simulation of chemical reactions in enzymes and other complex systems. The reaction path potential follows the ideas from the reaction path Hamiltonian of Miller, Handy and Adams for gas phase chemical reactions but is designed specifically for large systems that are described with combined ab initio quantum mechanical and molecular mechanical methods. The reaction path potential is an analytical energy expression of the combined quantum mechanical and molecular mechanical potential energy along the minimum energy path. An expansion around the minimum energy path is made in both the nuclear and the electronic degrees of freedom for the QM subsystem internal energy, while the energy of the subsystem described with MM remains unchanged from that in the combined quantum mechanical and molecular mechanical expression and the electrostatic interaction between the QM and MM subsystems is described as the interaction of the MM charges with the QM charges. The QM charges are polarizable in response to the changes in both the MM and the QM degrees of freedom through a new response kernel developed in the present work. The input data for constructing the reaction path potential are energies, vibrational frequencies, and electron density response properties of the QM subsystem along the minimum energy path, all of which can be obtained from the combined quantum mechanical and molecular mechanical calculations. Once constructed, it costs much less for its evaluation. Thus, the reaction path potential provides a potential energy surface for rigorous statistical mechanics and reaction dynamics calculations of complex systems. As an example, the method is applied to the statistical mechanical calculations for the potential of mean force of the chemical reaction in triosephosphate isomerase.     

Application of error-ranked singular value decomposition for the determination of potential-derived atomic-centered point charges

The present work provides a detailed investigation on the use of singular value decomposition (SVD) to solve the linear least-squares problem (LLS) for the purposes of obtaining potential-derived atom-centered point charges (PD charges) from the ab initio molecular electrostatic potential (V(QM)). Given the SVD of any PD charge calculation LLS problem, it was concluded that (1) all singular vectors are not necessary to obtain the optimal set of PD charges and (2) the most effective set of singular vectors do not necessarily correspond to those with the largest singular values. It is shown that the efficient use of singular vectors can provide statistically well-defined PD charges when compared with conventional PD charge calculation methods without sacrificing the agreement with V(QM). As can be expected, the methodology outlined here is independent of the algorithm for sampling V(QM) as well as the basis set used to calculate V(QM). An algorithm is provided to select the best set of singular vectors used for optimal PD charge calculations. To minimize the subjective comparisons of different PD charge sets, we also provide an objective criterion for determining if two sets of PD charges are significantly different from one another.     

Atomic charges derived from semiempirical methods

It is demonstrated that semiempirical methods give electrostatic potential (ESP) derived atomic point charges that are in reasonable agreement with ab initio ESP charges. Furthermore, we find that MNDO ESP charges are superior to AM1 ESP charges in correlating with ESP charges derived from the 6-31G* basis set. Thus, it is possible to obtain 6-31G* quality point charges by simply scaling MNDO ESP charges. The charges are scaled in a linear (y = Mx) manner to conserve charge. In this way researchers desiring to carry out force field simulations or minimizations can obtain charges by using MNDO, which requires much less computer time than the corresponding 6-31G* calculation.     

Effective atomic charges in alanine dipeptide

A new set of effective atomic charges of different conformers of alanine dipeptide is presented. These charges are obtained by fitting the electrostatic potential resulting from the ab initio SCF wave function of the system obtained in a 6-31G basis set. A specific fit procedure is used providing charges weakly dependent on the fit points as well as on the geometry of the molecule. It is shown that these charges retain a reasonable chemical meaning. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 473–482, 1999     

On the use of AM1 and MNDO wave functions to compute accurate electrostatic charges

A new strategy to evaluate accurate electrostatic charges from semiempirical wave functions is reported. The rigorous quantum mechanical molecular electrostatic potentials computed from both MNDO and AM1 wave functions are fitted to the point-charge molecular electrostatic potential to obtain the electrostatic charges. The reliability of this strategy is tested by comparing the semiempirical electrostatic charges for 21 molecules with the semiempirical Mulliken charges and with the ab initio STO-3G and 6-31G* electrostatic charges. The ability of the dipoles derived from the semiempirical electrostatic and Mulliken charges as well as from the SCF charge distributions to reproduce the ab initio 6-31G* electrostatic dipoles and the gas phase experimental values is determined. The statistical analysis clearly point out the goodness of the semiempirical electrostatic charges, specially when the MNDO method is used. The excellent relationships found between the MNDO and 6-31G* electrostatic charges permit to define a scaling factor which allows to accurately reproduce the 6-31G* electrostatic charge distribution as well as the experimental dipoles from the semiempirical electrostatic charges.     

Communication: Quantum polarized fluctuating charge model: a practical method to include ligand polarizability in biomolecular simulations

We present a simple and practical method to include ligand electronic polarization in molecular dynamics (MD) simulation of biomolecular systems. The method involves periodically spawning quantum mechanical (QM) electrostatic potential (ESP) calculations on an extra set of computer processors using molecular coordinate snapshots from a running parallel MD simulation. The QM ESPs are evaluated for the small-molecule ligand in the presence of the electric field induced by the protein, solvent, and ion charges within the MD snapshot. Partial charges on ligand atom centers are fit through the multi-conformer restrained electrostatic potential (RESP) fit method on several successive ESPs. The RESP method was selected since it produces charges consistent with the AMBER/GAFF force-field used in the simulations. The updated charges are introduced back into the running simulation when the next snapshot is saved. The result is a simulation whose ligand partial charges continuously respond in real-time to the short-term mean electrostatic field of the evolving environment without incurring additional wall-clock time. We show that (1) by incorporating the cost of polarization back into the potential energy of the MD simulation, the algorithm conserves energy when run in the microcanonical ensemble and (2) the mean solvation free energies for 15 neutral amino acid side chains calculated with the quantum polarized fluctuating charge method and thermodynamic integration agree better with experiment relative to the Amber fixed charge force-field.     

Molecular electrostatic potentials and partial atomic charges from correlated wave functions: Applications to the electronic ground and excited states of 3-methylindole

Procedures have been developed to generate molecular electrostatic potentials based on correlated wave function from ab initio or semiempirical electronic structure programs. A new algorithm for point-wise sampling of the potential is described and used to obtain partial atomic charges via a linear, least squares fit between classical and quantum mechanical electrostatic potentials. The proposed sampling algorithm is efficient and promises to introduce less rotational variance in the potential derived partial charges than algorithms applied previously. Electrostatic potentials and fitted atomic charges from ab initio (HF/6–31G* and MP2/6-31G*) and semiempirical (INDO/S; HF, SECI, and SDCI) wave functions are presented for the electronic ground (S₀) and excited (¹L_b, ¹L_a) states of 3-methylindole. © 1992 by John Wiley & Sons, Inc.     

Force-field parametrization of retro-inverso modified residues: Development of torsional and electrostatic parameters

Summary Torsional and the electrostatic parameters for molecular mechanics studies of retro-inverso modified peptides have been developed using quantum mechanical calculations. The resulting parameters have been compared with those calculated for conventional peptides. Rotational profiles, which were obtained spanning the corresponding dihedral angle, were corrected by removing the energy contributions associated to changes in interactions different from torsion under study. For this purpose, the torsional energy associated to each point of the profiles was estimated as the corresponding quantum mechanical energy minus the bonding and nonbonding energy contributions produced by the perturbations that the variation of the spanned dihedral angle causes in the bond distances, bond angles and the other dihedral angles. These energies were calculated using force-field expressions. The corrected profiles were fitted to a three-term Fourier expansion to derive the torsional parameters. Atomic charges for retro-inverso modified residues were derived from the rigorously calculated quantum mechanical electrostatic potential. Furthermore, the reliability of electrostatic models based on geometry-dependent charges and fixed charges has been examined.     

Importance of accurate charges in molecular docking: quantum mechanical/molecular mechanical (QM/MM) approach

The extent to which accuracy of electric charges plays a role in protein-ligand docking is investigated through development of a docking algorithm, which incorporates quantum mechanical/molecular mechanical (QM/MM) calculations. In this algorithm, fixed charges of ligands obtained from force field parameterization are replaced by QM/MM calculations in the protein environment, treating only the ligands as the quantum region. The algorithm is tested on a set of 40 cocrystallized structures taken from the Protein Data Bank (PDB) and provides strong evidence that use of nonfixed charges is important. An algorithm, dubbed "Survival of the Fittest" (SOF) algorithm, is implemented to incorporate QM/MM charge calculations without any prior knowledge of native structures of the complexes. Using an iterative protocol, this algorithm is able in many cases to converge to a nativelike structure in systems where redocking of the ligand using a standard fixed charge force field exhibits nontrivial errors. The results demonstrate that polarization effects can play a significant role in determining the structures of protein-ligand complexes, and provide a promising start towards the development of more accurate docking methods for lead optimization applications.     

Ligand atom partial charges assignment for complementary electrostatic potentials

Summary The design of molecules to fit into the active site of receptors is a rapidly developing area of pharmacology and medicinal chemistry. A good ligand needs a suitable geometry and also appropriate electrostatic properties. The electrostatic properties of the ligand should complement those of the receptor. We present a method for the assignment of atom-centred point charges for a ligand, based on the electrostatic potential of the receptor. These point charges are chosen to give the best possible complementarity to the receptor electrostatic potential over the van der Waals surface of the ligand. We demonstrate that point charges can be chosen to give good electrostatic complementarity, and suggest that a molecule with similar electrostatic properties should bind well to the receptor.     

Comparison of atomic charges derived via different procedures

Atomic charges were obtained from ab initio molecular orbital calculations using a variety of procedures to compare them and assess their utility. Two procedures based on the molecular orbitals were examined, the Mulliken population analysis and the Weinhold–Reed Natural Population Analysis. Two procedures using the charge density distribution were included; the Hirshfeld procedure and Bader's Atoms in Molecules method. Charges also were derived by fitting the electrostatic potential (CHELPG) and making use of the atomic polar tensors (GAPT). The procedures were first examined for basis set independence, and then applied to a group of hydrocarbons. The dipole moments for these molecules were computed from the various atomic charges and compared to the total SCF dipole moments. This was followed by an examination of a series of substituted methanes, simple hydrides, and a group of typical organic compounds such as carbonyl derivatives, nitriles, and nitro compounds. In some cases, the ability of the charges to reproduce electrostatic potentials was examined. © John Wiley & Sons, Inc.     

Fast,efficient generation of high‐quality atomic charges. AM1‐BCC model: I. Method

The AM1‐BCC method quickly and efficiently generates high‐quality atomic charges for use in condensed‐phase simulations. The underlying features of the electron distribution including formal charge and delocalization are first captured by AM1 atomic charges for the individual molecule. Bond charge corrections (BCCs), which have been parameterized against the HF/6‐31G* electrostatic potential (ESP) of a training set of compounds containing relevant functional groups, are then added using a formalism identical to the consensus BCI (bond charge increment) approach. As a proof of the concept, we fit BCCs simultaneously to 45 compounds including O‐, N‐, and S‐containing functionalities, aromatics, and heteroaromatics, using only 41 BCC parameters. AM1‐BCC yields charge sets of comparable quality to HF/6‐31G* ESP‐derived charges in a fraction of the time while reducing instabilities in the atomic charges compared to direct ESP‐fit methods. We then apply the BCC parameters to a small “test set” consisting of aspirin, d ‐glucose, and eryodictyol; the AM1‐BCC model again provides atomic charges of quality comparable with HF/6‐31G* RESP charges, as judged by an increase of only 0.01 to 0.02 atomic units in the root‐mean‐square (RMS) error in ESP. Based on these encouraging results, we intend to parameterize the AM1‐BCC model to provide a consistent charge model for any organic or biological molecule. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 132–146, 2000