Automated computational screening of the thiol reactivity of substituted alkenes

AbstractElectrophilic olefins can react with the S–H moiety of cysteine side chains. The formation of a covalent adduct through this mechanism can result in the inhibition of an enzyme. The reactivity of an olefin towards cysteine depends on its functional groups. In this study, 325 reactions of thiol-Michael-type additions to olefins were modeled using density functional theory. All combinations of ethenes with hydrogen, methyl ester, amide, and cyano substituents were included. An automated workflow was developed to perform the construction, conformation search, minimization, and calculation of molecular properties for the reactant, carbanion intermediate, and thioether products for a model reaction of the addition of methanethiol to the electrophile. Known cysteine-reactive electrophiles present in the database were predicted to react exergonically with methanethiol through a carbanion with a stability in the 30–40 kcal mol−1 range. 13 other compounds in our database that are also present in the PubChem database have similar properties. Natural bond orbital parameters were computed and regression analysis was used to determine the relationship between properties of the olefin electronic structure and the product and intermediate stability. The stability of the intermediates is very sensitive to electronic effects on the carbon where the anionic charge is centered. The stability of the products is more sensitive to steric factors.

[1]  M T D Cronin,et al.  A review of the electrophilic reaction chemistry involved in covalent protein binding relevant to toxicity , 2011, Critical reviews in toxicology.

[2]  M. Potashman,et al.  Covalent modifiers: an orthogonal approach to drug design. , 2009, Journal of medicinal chemistry.

[3]  J. Hubbell,et al.  Systematic modulation of Michael-type reactivity of thiols through the use of charged amino acids. , 2001, Bioconjugate chemistry.

[4]  Giulio Rastelli,et al.  A computational workflow for the design of irreversible inhibitors of protein kinases , 2010, J. Comput. Aided Mol. Des..

[5]  Yanli Wang,et al.  PubChem: Integrated Platform of Small Molecules and Biological Activities , 2008 .

[6]  Scott B Ficarro,et al.  A structure-guided approach to creating covalent FGFR inhibitors. , 2010, Chemistry & biology.

[7]  Chris Morley,et al.  Open Babel: An open chemical toolbox , 2011, J. Cheminformatics.

[8]  P. Kollman,et al.  Automatic atom type and bond type perception in molecular mechanical calculations. , 2006, Journal of molecular graphics & modelling.

[9]  Shuo Zhou,et al.  CovalentDock: Automated covalent docking with parameterized covalent linkage energy estimation and molecular geometry constraints , 2013, J. Comput. Chem..

[10]  K. Houk,et al.  Transition states and energetics of nucleophilic additions of thiols to substituted α,β-unsaturated ketones: substituent effects involve enone stabilization, product branching, and solvation. , 2011, The Journal of organic chemistry.

[11]  Christopher N Bowman,et al.  Thiol-ene click chemistry. , 2010, Angewandte Chemie.

[12]  D. Nair,et al.  The Thiol-Michael Addition Click Reaction: A Powerful and Widely Used Tool in Materials Chemistry , 2014 .

[13]  M. Leippe,et al.  Novel peptidomimetics containing a vinyl ester moiety as highly potent and selective falcipain-2 inhibitors. , 2009, Journal of medicinal chemistry.

[14]  Michael W Deem,et al.  Parallel tempering: theory, applications, and new perspectives. , 2005, Physical chemistry chemical physics : PCCP.

[15]  A. Kamimura,et al.  Theoretical and experimental study on the stereoselectivity of Michael addition of alkoxide anion to nitro olefin , 1990 .

[16]  Nir London,et al.  Covalent Docking of Large Libraries for the Discovery of Chemical Probes , 2014, Nature chemical biology.

[17]  P. Kuzmič,et al.  Covalent EGFR inhibitor analysis reveals importance of reversible interactions to potency and mechanisms of drug resistance , 2013, Proceedings of the National Academy of Sciences.

[18]  Alessio Lodola,et al.  Quantum Mechanics/Molecular Mechanics Modeling of Covalent Addition between EGFR-Cysteine 797 and N-(4-Anilinoquinazolin-6-yl) Acrylamide , 2015, J. Chem. Inf. Model..

[19]  F. Weinhold,et al.  Natural population analysis , 1985 .

[20]  Guengerich Fp,et al.  In Vitro Metabolism of Acrylonitrile to 2-Cyanoethylene Oxide, Reaction with Glutathione, and Irreversible Binding to Proteins and Nucleic Acids , 1981 .

[21]  Jennifer M. Smith,et al.  Range-Separated DFT Functionals are Necessary to Model Thio-Michael Additions. , 2013, Journal of chemical theory and computation.

[22]  Michael M. Mysinger,et al.  Automated Docking Screens: A Feasibility Study , 2009, Journal of medicinal chemistry.

[23]  P. M. Edwards.  Neurotoxicity of acrylamide and its analogues and effects of these analogues and other agents on acrylamide neuropathy. , 1975, British journal of industrial medicine.

[24]  V. Barone,et al.  Toward reliable density functional methods without adjustable parameters: The PBE0 model , 1999 .

[25]  G. Schüürmann,et al.  Predicting Michael-acceptor reactivity and toxicity through quantum chemical transition-state calculations. , 2011, Organic & biomolecular chemistry.

[26]  A. Schäfer,et al.  Fully optimized contracted Gaussian basis sets of triple zeta valence quality for atoms Li to Kr , 1994 .

[27]  David Weininger,et al.  SMILES, a chemical language and information system. 1. Introduction to methodology and encoding rules , 1988, J. Chem. Inf. Comput. Sci..

[28]  B. Roux,et al.  Absolute hydration free energy scale for alkali and halide ions established from simulations with a polarizable force field. , 2006, The journal of physical chemistry. B.

[29]  Jianpeng Ma,et al.  CHARMM: The biomolecular simulation program , 2009, J. Comput. Chem..

[30]  T. H. Dunning.  Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen , 1989 .

[31]  Finn Verner Jensen,et al.  Bayesian networks , 1998, Data Mining and Knowledge Discovery Handbook.

[32]  J. Schwöbel,et al.  Prediction of michael-type acceptor reactivity toward glutathione. , 2010, Chemical research in toxicology.

[33]  U. Mrowietz,et al.  Reactivity of dimethyl fumarate and methylhydrogen fumarate towards glutathione and N-acetyl-L-cysteine--preparation of S-substituted thiosuccinic acid esters. , 2007, Bioorganic & medicinal chemistry.

[34]  L. Honigberg,et al.  Discovery of Selective Irreversible Inhibitors for Bruton’s Tyrosine Kinase , 2007, ChemMedChem.

[35]  S. Sieber,et al.  Electrophilic natural products and their biological targets. , 2012, Natural product reports.

[36]  J. Tomasi,et al.  Quantum mechanical continuum solvation models. , 2005, Chemical reviews.

[37]  J. Taunton,et al.  Targeting protein kinases with selective and semipromiscuous covalent inhibitors. , 2014, Methods in enzymology.

[38]  Junmei Wang,et al.  Development and testing of a general amber force field , 2004, J. Comput. Chem..

[39]  Benedetta Mennucci,et al.  Polarizable continuum model , 2012 .

[40]  Ramaswamy Nilakantan,et al.  Antitumor Activity of HKI-272, an Orally Active, Irreversible Inhibitor of the HER-2 Tyrosine Kinase , 2004, Cancer Research.

[41]  M. Dewar,et al.  Ground States of Molecules. 38. The MNDO Method. Approximations and Parameters , 1977 .

[43]  Jack Taunton,et al.  Reversible targeting of noncatalytic cysteines with chemically tuned electrophiles , 2012, Nature chemical biology.

[44]  M. Head‐Gordon,et al.  Systematic optimization of long-range corrected hybrid density functionals. , 2008, The Journal of chemical physics.