A Multilevel Strategy for the Exploration of the Conformational Flexibility of Small Molecules.

Predicting the conformational preferences of flexible compounds is still a challenging problem with important implications in areas such as molecular recognition and drug design. In this work, we describe a multilevel strategy to explore the conformational preferences of molecules. The method relies on the predominant-state approximation, which partitions the conformational space into distinct conformational wells. Moreover, it combines low-level (LL) methods for sampling the conformational minima and high-level (HL) techniques for calibrating their relative stability. In the implementation used in this study, the LL sampling is performed with the semiempirical RM1 Hamiltonian, and solvent effects are included using the RM1 version of the MST continuum solvation model. The HL refinement of the conformational wells is performed by combining geometry optimizations of the minima at the B3LYP (gas phase) or MST-B3LYP (solution) level, followed by single point MP2 computations using Dunning's augmented basis sets. Then, the effective free energy of a conformational well is estimated by combining the MP2 energy, supplemented with the MST-B3LYP solvation free energy for a conformational search in solution, with the local curvature of the well sampled at the semiempirical level. Applications of this strategy involve the exploration of the conformational preferences of 1,2-dichloroethane and neutral histamine in both the gas phase and water solution. Finally, the multilevel strategy is used to estimate the reorganization cost required for selecting the bioactive conformation of HIV reverse transcriptase inhibitors, which is estimated to be at most 1.3 kcal/mol.

[1]  Kevin E. Riley,et al.  Critical Assessment of the Performance of Density Functional Methods for Several Atomic and Molecular Properties. , 2007, Journal of chemical theory and computation.

[2]  W. Richards,et al.  Conformation of histamine derivatives. 1. Application of molecular orbital calculations and nuclear magnetic resonance spectroscopy. , 1973, Journal of medicinal chemistry.

[3]  C. Ganellin.  The tautomer ratio of histamine , 1973, The Journal of pharmacy and pharmacology.

[4]  R. R. Ison,et al.  Solution conformations of histamine and some related derivatives. , 1973, Journal of medicinal chemistry.

[5]  Stefan Grimme,et al.  Semiempirical GGA‐type density functional constructed with a long‐range dispersion correction , 2006, J. Comput. Chem..

[6]  Gisbert Schneider,et al.  Protein-ligand interactions from molecular recognition to drug design , 2003 .

[7]  J A McCammon,et al.  Accommodating protein flexibility in computational drug design. , 2000, Molecular pharmacology.

[8]  Xicheng Wang,et al.  Cyndi: a multi-objective evolution algorithm based method for bioactive molecular conformational generation , 2009, BMC Bioinformatics.

[9]  Stephen R. Johnson,et al.  Molecular properties that influence the oral bioavailability of drug candidates. , 2002, Journal of medicinal chemistry.

[10]  F. J. Luque,et al.  Protein flexibility and ligand recognition: challenges for molecular modeling. , 2011, Current topics in medicinal chemistry.

[11]  Michael K. Gilson,et al.  ''Mining minima'': Direct computation of conformational free energy , 1997 .

[12]  P. Charifson,et al.  Conformational analysis of drug-like molecules bound to proteins: an extensive study of ligand reorganization upon binding. , 2004, Journal of medicinal chemistry.

[13]  P. Ugliengo,et al.  B3LYP augmented with an empirical dispersion term (B3LYP-D*) as applied to molecular crystals , 2008 .

[14]  Minoru Katō,et al.  Raman study of the trans–gauche conformational equilibrium of 1,2-dichloroethane in water: Experimental evidence for the hydrophobic effect , 1999 .

[15]  I. Tuñón,et al.  Structural and vibrational study of the tautomerism of histamine free-base in solution. , 2003, Journal of the American Chemical Society.

[16]  D. Truhlar,et al.  The M06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, and transition elements: two new functionals and systematic testing of four M06-class functionals and 12 other functionals , 2008 .

[17]  Martin Stahl,et al.  Small Molecule Conformational Preferences Derived from Crystal Structure Data. A Medicinal Chemistry Focused Analysis , 2008, J. Chem. Inf. Model..

[18]  C. Dobson,et al.  Acid-base and tautomeric equilibriums in the solid state: nitrogen-15 NMR spectroscopy of histidine and imidazole , 1982 .

[19]  Keith T. Butler,et al.  Toward accurate relative energy predictions of the bioactive conformation of drugs , 2009, J. Comput. Chem..

[20]  F. Javier Luque,et al.  Extension of the MST continuum solvation model to the RM1 semiempirical hamiltonian , 2008, J. Comput. Chem..

[21]  P. Godfrey,et al.  Proportions of Species Observed in Jet Spectroscopy−Vibrational Energy Effects: Histamine Tautomers and Conformers , 1998 .

[22]  J. Tomasi,et al.  Solvent Effects on trans/gauche Conformational Equilibria of Substituted Chloroethanes: a Polarizable Continuum Model Study , 2001 .

[23]  K. Mikkelsen,et al.  A coupled-cluster solvent reaction field method , 1999 .

[24]  E. Vilaseca.  Solvent effect on conformational equilibrium: A Monte Carlo study of 1,2‐dichloroethane in carbon tetrachloride , 1996 .

[25]  William L. Jorgensen,et al.  Importance of Polarization for Dipolar Solutes in Low-Dielectric Media: 1,2-Dichloroethane and Water in Cyclohexane , 1995 .

[26]  Stefan Grimme,et al.  Accurate description of van der Waals complexes by density functional theory including empirical corrections , 2004, J. Comput. Chem..

[27]  M. Orozco,et al.  Solvent effects on tautomerism equilibria in heterocycles , 1997 .

[28]  M. Frisch,et al.  Ab Initio Calculation of Vibrational Absorption and Circular Dichroism Spectra Using Density Functional Force Fields , 1994 .

[29]  C. Hunter,et al.  Quantifying intermolecular interactions: guidelines for the molecular recognition toolbox. , 2004, Angewandte Chemie.

[30]  M. Gilson,et al.  Ligand configurational entropy and protein binding , 2007, Proceedings of the National Academy of Sciences.

[31]  Ronald M. Levy,et al.  Conformational populations of ligand‐sized molecules by replica exchange molecular dynamics and temperature reweighting , 2009, J. Comput. Chem..

[32]  E. Vilaseca,et al.  Solvent Effect on the Conformational Equilibrium of 1,2-Dichloroethane in Water. The Role of Solute Polarization , 2004 .

[33]  B. Noszál,et al.  Concentration and basicity of histamine rotamers , 2002 .

[34]  Julian Tirado-Rives,et al.  Contribution of conformer focusing to the uncertainty in predicting free energies for protein-ligand binding. , 2006, Journal of medicinal chemistry.

[35]  Modesto Orozco,et al.  Prediction of Conformational Free Energy Differences of Solutes in Solution∶ An MC-MST Study , 2002 .

[36]  Graham J. Durant,et al.  Theoretical Analyses of the Tautomeric and Conformational Equilibria of Histamine and (.alpha.R,.beta.S)-.alpha.,.beta.-Dimethylhistamine in the Gas Phase and Aqueous Solution , 1994 .

[37]  Yan Zhao,et al.  Exchange-correlation functional with broad accuracy for metallic and nonmetallic compounds, kinetics, and noncovalent interactions. , 2005, The Journal of chemical physics.

[38]  Thom Vreven,et al.  Chapter 3 Hybrid Methods: ONIOM(QM:MM) and QM/MM , 2006 .

[39]  Claudio N. Cavasotto,et al.  Quantum mechanical binding free energy calculation for phosphopeptide inhibitors of the Lck SH2 domain , 2011, J. Comput. Chem..

[40]  András Aszódi,et al.  Analyzing the performance of conformational search programs on compound databases. , 2007, Journal of molecular graphics & modelling.

[41]  Woody Sherman,et al.  ConfGen: A Conformational Search Method for Efficient Generation of Bioactive Conformers , 2010, J. Chem. Inf. Model..

[42]  F. Javier Luque,et al.  Performance of the IEF-MST solvation continuum model in the SAMPL2 blind test prediction of hydration and tautomerization free energies , 2010, J. Comput. Aided Mol. Des..

[43]  F. Javier Luque,et al.  Extension of the MST model to the IEF formalism: HF and B3LYP parametrizations , 2005 .

[44]  Thierry Langer,et al.  Comparative Performance Assessment of the Conformational Model Generators Omega and Catalyst: A Large-Scale Survey on the Retrieval of Protein-Bound Ligand Conformations , 2006, J. Chem. Inf. Model..

[45]  Eric J. Martin,et al.  Conformational Sampling of Bioactive Molecules: A Comparative Study , 2007, J. Chem. Inf. Model..

[46]  Donald G Truhlar,et al.  Free-energy surfaces for liquid-phase reactions and their use to study the border between concerted and nonconcerted alpha,beta-elimination reactions of esters and thioesters. , 2010, Journal of the American Chemical Society.

[47]  Alfredo Mayall Simas,et al.  RM1: A reparameterization of AM1 for H, C, N, O, P, S, F, Cl, Br, and I , 2006, J. Comput. Chem..

[48]  Claudio N. Cavasotto and Narender Singh.  Docking and High Throughput Docking: Successes and the Challenge of Protein Flexibility , 2008 .

[49]  L. Mario Amzel,et al.  Impact of linker strain and flexibility in the design of a fragment-based inhibitor , 2009, Nature chemical biology.

[50]  Jonas Boström,et al.  Reproducing the conformations of protein-bound ligands: A critical evaluation of several popular conformational searching tools , 2001, J. Comput. Aided Mol. Des..

[51]  M. Murcko,et al.  Solvent Effects on 1,2-Dihaloethane Gauche/Trans Ratios , 1995 .

[52]  Jonas Boström,et al.  Assessing the performance of OMEGA with respect to retrieving bioactive conformations. , 2003, Journal of molecular graphics & modelling.

[53]  M. Gilson,et al.  Calculation of cyclodextrin binding affinities: energy, entropy, and implications for drug design. , 2004, Biophysical journal.

[54]  B. Vogelsanger,et al.  Rotational spectra of biomolecules : histamine , 1991 .

[55]  G. Klebe,et al.  Approaches to the description and prediction of the binding affinity of small-molecule ligands to macromolecular receptors. , 2002, Angewandte Chemie.

[56]  Jiabo Li,et al.  CAESAR: A New Conformer Generation Algorithm Based on Recursive Buildup and Local Rotational Symmetry Consideration , 2007, J. Chem. Inf. Model..

[57]  Michiko Amano,et al.  Novel Method for the Evaluation of 3D Conformation Generators , 2009, J. Chem. Inf. Model..

[58]  John D. Roberts,et al.  A nitrogen-15 nuclear magnetic resonance study of the acid-base and tautomeric equilibriums of 4-substituted imidazoles and its relevance to the catalytic mechanism of .alpha.-lytic protease , 1982 .

[59]  W. Bachovchin,et al.  Direct observation of the tautomeric forms of histidine in nitrogen-15 NMR spectra at low temperatures. Comments on intramolecular hydrogen bonding and on tautomeric equilibrium constants , 1993 .