Ab Initio Stochastic Optimization of Conformational and Many-Body Degrees of Freedom

Moderate to large size molecules in solution have complex energy surfaces due to intramolecular (conformational) and intermolecular (many-body) interactions. The first principles Monte Carlo (FPMC) method, previously shown to effectively locate minimum-energy structures for systems with only many-body complexity, has been extended to address conformational flexibility by adding three new Monte Carlo move types. The primary advantage of the FPMC method is the ability to efficiently locate minimum energy structures of molecules with conformational flexibility in the presence of explicit solvent molecules using highly accurate quantum chemical calculations. The additions to FPMC were validated by studying conformers of glycerol, glyceraldehyde, and a large humic acid monomer unit. The structure of glyceraldehyde in the presence of one and two water molecules was also explored to demonstrate the power of FPMC to study systems with both conformational and many-body degrees of freedom.

[1]  C. D. Gelatt,et al.  Optimization by Simulated Annealing , 1983, Science.

[2]  C. Steelink,et al.  Humic Acids: Marvelous Products of Soil Chemistry , 2001 .

[3]  A. Leach.  A Survey of Methods for Searching the Conformational Space of Small and Medium‐Sized Molecules , 1992 .

[4]  T. Cieplak,et al.  A New Effective Algorithm for the Unambiguous Identification of the Stereochemical Characteristics of Compounds During Their Registration in Databases , 2001, Molecules : A Journal of Synthetic Chemistry and Natural Product Chemistry.

[5]  Jay W. Ponder,et al.  Accurate modeling of the intramolecular electrostatic energy of proteins , 1995, J. Comput. Chem..

[6]  M. Ľożyński,et al.  MP2 and Density Functional Studies of Hydrogen Bonding in Model Trioses: d-(+)-Glyceraldehyde and Dihydroxyacetone , 1997 .

[7]  Jay W. Ponder,et al.  Analysis and Application of Potential Energy Smoothing and Search Methods for Global Optimization , 1998 .

[8]  J. Stewart.  Optimization of parameters for semiempirical methods I. Method , 1989 .

[9]  J. Ponder,et al.  Calculation of the reaction field due to off-center point multipoles , 1997 .

[10]  R. Cramer,et al.  Validation of the general purpose tripos 5.2 force field , 1989 .

[11]  Jay W. Ponder,et al.  Algorithms for calculating excluded volume and its derivatives as a function of molecular conformation and their use in energy minimization , 1991 .

[12]  Masaaki Kojima,et al.  Calculation of stacking geometries by a modified MM2 method , 2001, J. Comput. Chem..

[13]  M. Karplus,et al.  CHARMM: A program for macromolecular energy, minimization, and dynamics calculations , 1983 .

[14]  J. Ponder,et al.  An efficient newton‐like method for molecular mechanics energy minimization of large molecules , 1987 .

[15]  H. Scheraga,et al.  Monte Carlo-minimization approach to the multiple-minima problem in protein folding. , 1987, Proceedings of the National Academy of Sciences of the United States of America.

[16]  S. J. Singer,et al.  Computational analysis of the potential energy surfaces of glycerol in the gas and aqueous phases: effects of level of theory, basis set, and solvation on strongly intramolecularly hydrogen-bonded systems. , 2001, Journal of the American Chemical Society.

[17]  Linda J. Broadbelt,et al.  Many-Body Optimization Using an Ab Initio Monte Carlo Method , 2003, J. Chem. Inf. Comput. Sci..

[18]  S. Jansen,et al.  Conformational Modeling of a New Building Block of Humic Acid: Approaches to the Lowest Energy Conformer , 1999 .

[19]  J. Ponder,et al.  Protein structure prediction using a combination of sequence homology and global energy minimization: II. Energy functions , 1998 .

[20]  S. D. Cunningham,et al.  SEQUESTRATION OF HYDROPHOBIC ORGANIC CONTAMINANTS BY GEOSORBENTS , 1997 .

[21]  Ab initio GB study of conformational space of d-glyceraldehyde in aqueous solution , 2000 .

[22]  Peter A. Kollman,et al.  AMBER, a package of computer programs for applying molecular mechanics, normal mode analysis, molecular dynamics and free energy calculations to simulate the structural and energetic properties of molecules , 1995 .

[23]  D. Ferguson,et al.  A new approach to probing conformational space with molecular mechanics: Random incremental pulse search , 1989 .

[24]  Edward P. Jaeger,et al.  Conformational searching methods for small molecules. III. Study of stochastic methods available in SYBYL and MACROMODEL , 1996, J. Comput. Chem..

[25]  W. Buermann,et al.  Tight metal binding by humic acids and its role in biomineralization , 1997 .

[26]  Norman L. Allinger,et al.  Molecular mechanics. The MM3 force field for hydrocarbons. 1 , 1989 .

[27]  Donald G. Truhlar,et al.  MODEL FOR AQUEOUS SOLVATION BASED ON CLASS IV ATOMIC CHARGES AND FIRST SOLVATION SHELL EFFECTS , 1996 .