Matthew Clark

Abstract A molecular mechanics force field implemented in the Sybyl program is described along with a statistical evaluation of its efficiency on a variety of compounds by analysis of internal coordinates and thermodynamic barriers. The goal of the force field is to provide good quality geometries and relative energies for a large variety of organic molecules by energy minimization. Performance in protein modeling was tested by minimizations starting from crystallographic coordinates for three cyclic hexapeptides in the crystal lattice with rms movements of 0.019 angstroms, 2.06 degrees, and 6.82 degrees for bond lengths, angles, and torsions, respectively, and an rms movement of 0.16 angstroms for heavy atoms. Isolated crambin was also analyzed with rms movements of 0.025 angstroms, 2.97 degrees, and 13.0 degrees for bond lengths, angles, and torsions respectively, and an rms movement of 0.42 angstroms for heavy atoms. Accuracy in calculating thermodynamic barriers was tested for 17 energy differences between conformers, 12 stereoisomers, and 15 torsional barriers. The rms errors were 0.8, 1.7, and 1.13 kcal/mol, respectively, for the three tests. Performance in general purpose applications was assessed by minimizing 76 diverse complex organic crystal structures, with and without randomization by coordinate truncation, with rms movements of 0.025 angstroms, 2.50 degrees, and 9.54 degrees for bond lengths, angles and torsions respectively, and an average rms movement of 0.192 angstroms for heavy atoms.

Abstract The frequency of chance correlation using partial least squares (PLS) has been measured experimentally for variously dimensioned data, comprising either completely random numbers, random numbers containing a perfect correlation within, and CoMFA field descriptors. This frequency, much lower than that for stepwise multiple regression, is maximal for datasets in which the number of descriptors equals the number of compounds, and surprisingly decreases indefinitely as the number of descriptors becomes much greater than the number of compounds. However, perfect correlations involving descriptor subsets are not detected by PLS if the number of irrelevant descriptors is excessive. In CoMFA applications, the probability of chance correlation is usually negligible. For example with 21 compounds a crossvalidated r 2 value greater than 0.25 will occur by chance in less than 5% of trials.

We describe a new interface of the GPU parallelized TeraChem electronic structure package and the Amber molecular dynamics package for quantum mechanical (QM) and mixed QM and molecular mechanical (MM) molecular dynamics simulations. This QM/MM interface is used for computation of the absorption spectra of the photoactive yellow protein (PYP) chromophore in vacuum, aqueous solution, and protein environments. The computed excitation energies of PYP require a very large QM region (hundreds of atoms) covalently bonded to the chromophore in order to achieve agreement with calculations that treat the entire protein quantum mechanically. We also show that 40 or more surrounding water molecules must be included in the QM region in order to obtain converged excitation energies of the solvated PYP chromophore. These results indicate that large QM regions (with hundreds of atoms) are a necessity in QM/MM calculations.