ForceGen 3D structure and conformer generation: from small lead-like molecules to macrocyclic drugs

We introduce the ForceGen method for 3D structure generation and conformer elaboration of drug-like small molecules. ForceGen is novel, avoiding use of distance geometry, molecular templates, or simulation-oriented stochastic sampling. The method is primarily driven by the molecular force field, implemented using an extension of MMFF94s and a partial charge estimator based on electronegativity-equalization. The force field is coupled to algorithms for direct sampling of realistic physical movements made by small molecules. Results are presented on a standard benchmark from the Cambridge Crystallographic Database of 480 drug-like small molecules, including full structure generation from SMILES strings. Reproduction of protein-bound crystallographic ligand poses is demonstrated on four carefully curated data sets: the ConfGen Set (667 ligands), the PINC cross-docking benchmark (1062 ligands), a large set of macrocyclic ligands (182 total with typical ring sizes of 12–23 atoms), and a commonly used benchmark for evaluating macrocycle conformer generation (30 ligands total). Results compare favorably to alternative methods, and performance on macrocyclic compounds approaches that observed on non-macrocycles while yielding a roughly 100-fold speed improvement over alternative MD-based methods with comparable performance.

[1]  J. Gasteiger,et al.  ITERATIVE PARTIAL EQUALIZATION OF ORBITAL ELECTRONEGATIVITY – A RAPID ACCESS TO ATOMIC CHARGES , 1980 .

[2]  Nicolas Foloppe,et al.  Tackling the conformational sampling of larger flexible compounds and macrocycles in pharmacology and drug discovery. , 2013, Bioorganic & medicinal chemistry.

[3]  Thomas A. Halgren,et al.  Merck molecular force field. II. MMFF94 van der Waals and electrostatic parameters for intermolecular. interactions , 1996, J. Comput. Chem..

[4]  W. Goddard,et al.  Charge equilibration for molecular dynamics simulations , 1991 .

[5]  Thomas A. Halgren,et al.  Merck molecular force field. IV. conformational energies and geometries for MMFF94 , 1996, J. Comput. Chem..

[6]  Roger Fletcher,et al.  A Rapidly Convergent Descent Method for Minimization , 1963, Comput. J..

[7]  Wilfried J. Mortier,et al.  Electronegativity-equalization method for the calculation of atomic charges in molecules , 1986 .

[8]  Zhong-Zhi Yang,et al.  General atom-bond electronegativity equalization method and its application in prediction of charge distributions in polypeptide , 2000 .

[9]  John Mullay,et al.  A simple method for calculating atomic charges in charged molecular systems of biochemical interest , 1991 .

[10]  Satoshi Omura,et al.  Creation of Customized Bioactivity within a 14-Membered Macrolide Scaffold: Design, Synthesis, and Biological Evaluation Using a Family-18 Chitinase. , 2015, Journal of medicinal chemistry.

[11]  Huichun Zhu,et al.  Discovery of (10R)-7-Amino-12-Fluoro-2,10,16-Trimethyl-15-Oxo-10,15,16,17-Tetrahydro-2H-8,4-(Metheno)Pyrazolo[4,3-H][2,5,11]Benzoxadiazacyclotetradecine-3-Carbonitrile (Pf-06463922), a Macrocyclic Inhibitor of Alk/Ros1 with Pre-Clinical Brain Exposure and Broad Spectrum Potency Against Alk-Resistant , 2014 .

[12]  Karl-Heinz Altmann,et al.  Pyridomycin bridges the NADH- and substrate-binding pockets of the enoyl reductase InhA. , 2014, Nature chemical biology.

[13]  Ajay N. Jain Surflex: fully automatic flexible molecular docking using a molecular similarity-based search engine. , 2003, Journal of medicinal chemistry.

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

[15]  Jürgen Bajorath,et al.  Predicting bioactive conformations and binding modes of macrocycles , 2016, Journal of Computer-Aided Molecular Design.

[16]  Hong Cao,et al.  The Molecular Basis of Drug Resistance against Hepatitis C Virus NS3/4A Protease Inhibitors , 2012, PLoS pathogens.

[17]  Wei Liu,et al.  Discovery of (10R)-7-amino-12-fluoro-2,10,16-trimethyl-15-oxo-10,15,16,17-tetrahydro-2H-8,4-(metheno)pyrazolo[4,3-h][2,5,11]-benzoxadiazacyclotetradecine-3-carbonitrile (PF-06463922), a macrocyclic inhibitor of anaplastic lymphoma kinase (ALK) and c-ros oncogene 1 (ROS1) with preclinical brain expos , 2014, Journal of medicinal chemistry.

[18]  D. Hazuda,et al.  MK-7009, a Potent and Selective Inhibitor of Hepatitis C Virus NS3/4A Protease , 2009, Antimicrobial Agents and Chemotherapy.

[19]  K Osterlund,et al.  Unexpected binding mode of a cyclic sulfamide HIV-1 protease inhibitor. , 1997, Journal of medicinal chemistry.

[20]  Thomas A. Halgren,et al.  Merck molecular force field. V. Extension of MMFF94 using experimental data, additional computational data, and empirical rules , 1996, J. Comput. Chem..

[21]  Ajay N. Jain,et al.  Knowledge-guided docking: accurate prospective prediction of bound configurations of novel ligands using Surflex-Dock , 2015, Journal of Computer-Aided Molecular Design.

[22]  F. Allen The Cambridge Structural Database: a quarter of a million crystal structures and rising. , 2002, Acta crystallographica. Section B, Structural science.

[23]  Charlotte M. Deane,et al.  Freely Available Conformer Generation Methods: How Good Are They? , 2012, J. Chem. Inf. Model..

[24]  Li Xing,et al.  Binding site elucidation and structure guided design of macrocyclic IL-17A antagonists , 2016, Scientific Reports.

[25]  J. Gasteiger,et al.  Automatic generation of 3D-atomic coordinates for organic molecules , 1990 .

[26]  Atsushi Suda,et al.  Design and synthesis of novel macrocyclic 2-amino-6-arylpyrimidine Hsp90 inhibitors. , 2012, Bioorganic & medicinal chemistry letters.

[27]  Paul Labute,et al.  LowModeMD - Implicit Low-Mode Velocity Filtering Applied to Conformational Search of Macrocycles and Protein Loops , 2010, J. Chem. Inf. Model..

[28]  Benjamin A. Ellingson,et al.  Conformer Generation with OMEGA: Algorithm and Validation Using High Quality Structures from the Protein Databank and Cambridge Structural Database , 2010, J. Chem. Inf. Model..

[29]  Michael K. Gilson,et al.  Fast Assignment of Accurate Partial Atomic Charges: An Electronegativity Equalization Method that Accounts for Alternate Resonance Forms , 2003, J. Chem. Inf. Comput. Sci..

[30]  Thomas A. Halgren Merck molecular force field. III. Molecular geometries and vibrational frequencies for MMFF94 , 1996, J. Comput. Chem..

[31]  Matthew P Jacobson,et al.  Exhaustive Conformational Sampling of Complex Fused Ring Macrocycles Using Inverse Kinematics. , 2016, Journal of chemical theory and computation.

[32]  L. Van,et al.  Structural outline of the detailed mechanism for elongation factor Ts-mediated guanine nucleotide exchange on elongation factor Tu. , 2015, Journal of structural biology.

[33]  Ajay N. Jain Effects of protein conformation in docking: improved pose prediction through protein pocket adaptation , 2009, J. Comput. Aided Mol. Des..

[34]  Thomas A. Halgren MMFF VI. MMFF94s option for energy minimization studies , 1999, J. Comput. Chem..

[35]  Thomas A. Halgren Merck molecular force field. I. Basis, form, scope, parameterization, and performance of MMFF94 , 1996, J. Comput. Chem..

[36]  Ajay N. Jain,et al.  Robust ligand-based modeling of the biological targets of known drugs. , 2006, Journal of medicinal chemistry.

[37]  Ajay N. Jain Morphological similarity: A 3D molecular similarity method correlated with protein-ligand recognition , 2000, J. Comput. Aided Mol. Des..

[38]  Ajay N. Jain,et al.  Prediction of Off-Target Drug Effects Through Data Fusion , 2014, Pacific Symposium on Biocomputing.

[39]  Anders Dahlgren,et al.  Highly potent macrocyclic BACE-1 inhibitors incorporating a hydroxyethylamine core: design, synthesis and X-ray crystal structures of enzyme inhibitor complexes. , 2012, Bioorganic & medicinal chemistry.

[40]  Harold A. Scheraga,et al.  Determination of net atomic charges using a modified partial equalization of orbital electronegativity method. 1. Application to neutral molecules as models for polypeptides , 1990 .

[41]  Wilfried J. Mortier,et al.  Electronegativity-equalization method for the calculation of atomic charges in molecules , 1986 .