Design, 3D QSAR modeling and docking of TGF-β type I inhibitors to target cancer

Transforming growth factor-β (TGF-β) family members plays a vital role in regulating hormonal function, bone formation, tissue remodeling, and erythropoiesis, cell growth and apoptosis. TGF-β super-family members mediate signal transduction via serine/threonine kinase receptors located on the cell membrane. Variation in expression of the TGF-β type I and II receptors in the cancer cells compromise its tumor suppressor activities which direct to oncogenic functions. The present study was aimed to screen the potent TGF-β type I inhibitors through atom based 3D-QSAR and pharmacophore modelling. For this purpose, we have chosen known TGF-β type I inhibitors from the binding database. The PHASE module of Schrodinger identified the best Pharmacophore model which includes three hydrogen bond acceptors (A), one hydrophobic region (H), and one ring (R) as the structural features. The top pharmacophore model AAAHR.27 was chosen with the R2 value of 0.94 and validated externally using molecules of the test set. Moreover the AAAHR.27 model underwent virtual screening using the molecules from the NCI, ZINC and Maybridge database. The screened molecules were further filtered using molecular docking and ADME properties prediction. Additionally, the 7 lead molecules were compared with a newly identified compound "SB431542" (well known TGF-β type I receptor inhibitor) and Galunisertib, a small molecule inhibitor of TGF-β type I, under clinical development (phase II trials) using the docking score and other binding properties. Also a top scored screened molecule from our study has been compared and confirmed using molecular dynamic simulation studies. By this way, we have obtained 7 distinct drug-like TGF-β type I inhibitors which can be beneficial in suppressing cancers reported with up-regulation of TGF-β type I. This result highlights the guidelines for designing molecules with TGF-β Type I inhibitory properties.

[1]  David E. Shaw,et al.  PHASE: a new engine for pharmacophore perception, 3D QSAR model development, and 3D database screening: 1. Methodology and preliminary results , 2006, J. Comput. Aided Mol. Des..

[2]  Brian K. Shoichet,et al.  ZINC - A Free Database of Commercially Available Compounds for Virtual Screening , 2005, J. Chem. Inf. Model..

[3]  Donald B. Boyd,et al.  Successes of Computer‐Assisted Molecular Design , 2007 .

[4]  Scott R Granter,et al.  EWI-2 negatively regulates TGF-β signaling leading to altered melanoma growth and metastasis , 2015, Cell Research.

[5]  Dae-Kee Kim,et al.  Synthesis and biological evaluation of 4(5)-(6-alkylpyridin-2-yl)imidazoles as transforming growth factor-beta type 1 receptor kinase inhibitors. , 2007, Journal of medicinal chemistry.

[6]  David Calkins,et al.  Towards the comprehensive, rapid, and accurate prediction of the favorable tautomeric states of drug-like molecules in aqueous solution , 2010, J. Comput. Aided Mol. Des..

[7]  Anita B. Roberts,et al.  Tumor suppressor activity of the TGF-β pathway in human cancers , 1996 .

[8]  Ruth Lehr,et al.  Identification of novel inhibitors of the transforming growth factor beta1 (TGF-beta1) type 1 receptor (ALK5). , 2002, Journal of medicinal chemistry.

[9]  Baidurya Bhattacharya,et al.  Nonergodicity of the Nose-Hoover chain thermostat in computationally achievable time. , 2014, Physical review. E, Statistical, nonlinear, and soft matter physics.

[10]  Zhili Zuo,et al.  Docking study and three-dimensional quantitative structure-activity relationship (3D-QSAR) analyses of transforming growth factor-ß type I receptor kinase inhibitors , 2009 .

[11]  Hege S. Beard,et al.  Glide: a new approach for rapid, accurate docking and scoring. 2. Enrichment factors in database screening. , 2004, Journal of medicinal chemistry.

[12]  Shanthi Nagarajan,et al.  Homology modeling, molecular dynamics, e-pharmacophore mapping and docking study of Chikungunya virus nsP2 protease , 2011, Journal of Molecular Modeling.

[13]  Oliver Eickelberg,et al.  The impact of TGF-β on lung fibrosis: from targeting to biomarkers. , 2012, Proceedings of the American Thoracic Society.

[14]  A. Hata,et al.  Targeting the TGFβ signalling pathway in disease , 2012, Nature Reviews Drug Discovery.

[15]  Ying E. Zhang,et al.  Smad-dependent and Smad-independent pathways in TGF-β family signalling , 2003, Nature.

[16]  J. Massagué,et al.  How cells read TGF-β signals , 2000, Nature Reviews Molecular Cell Biology.

[17]  Sahaya Asirvatham,et al.  Quantitative structure activity relationships studies of non-steroidal anti-inflammatory drugs: A review , 2016, Arabian Journal of Chemistry.

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

[19]  D. Kingsley,et al.  The TGF-beta superfamily: new members, new receptors, and new genetic tests of function in different organisms. , 1994, Genes & development.

[20]  W. L. Jorgensen,et al.  The OPLS [optimized potentials for liquid simulations] potential functions for proteins, energy minimizations for crystals of cyclic peptides and crambin. , 1988, Journal of the American Chemical Society.

[21]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[22]  A. De,et al.  In silico Analysis of 3D-QSAR and Molecular Docking for Bcl-2 Inhibitors to Potential Anticancer Drug Development , 2013 .

[23]  Richard A. Friesner,et al.  Flexible ligand docking with Glide. , 2007, Current protocols in bioinformatics.

[24]  Jian Luo,et al.  Inhibition of breast cancer metastases by a novel inhibitor of TGFβ receptor 1. , 2013, Journal of the National Cancer Institute.

[25]  W. L. Jorgensen,et al.  Temperature dependence of TIP3P, SPC, and TIP4P water from NPT Monte Carlo simulations: Seeking temperatures of maximum density , 1998 .

[26]  I. Kuntz,et al.  Automated docking with grid‐based energy evaluation , 1992 .

[27]  Alexander Tropsha,et al.  Best Practices for QSAR Model Development, Validation, and Exploitation , 2010, Molecular informatics.

[28]  T. Cheatham,et al.  Convergence and reproducibility in molecular dynamics simulations of the DNA duplex d(GCACGAACGAACGAACGC). , 2015, Biochimica et biophysica acta.

[29]  K. Muthusamy,et al.  Ligand-based Pharmacophore Modeling; Atom-based 3D-QSAR Analysis and Molecular Docking Studies of Phosphoinositide-Dependent Kinase-1 Inhibitors , 2012, Indian journal of pharmaceutical sciences.

[30]  Françoise Gellibert,et al.  Identification of 1,5-naphthyridine derivatives as a novel series of potent and selective TGF-beta type I receptor inhibitors. , 2004, Journal of medicinal chemistry.

[31]  Ya-wen Wang,et al.  Pharmacophore modeling and 3D-QSAR analysis of phosphoinositide 3-kinase p110α inhibitors , 2010, Journal of molecular modeling.

[32]  Werner J Geldenhuys,et al.  3D-QSAR and docking studies on transforming growth factor (TGF)-beta receptor 1 antagonists. , 2010, Bioorganic & medicinal chemistry letters.

[33]  R. McAnulty,et al.  TGF-beta: its role in asthma and therapeutic potential. , 2006, Current drug targets.

[34]  Richard A. Friesner,et al.  Comparative Performance of Several Flexible Docking Programs and Scoring Functions: Enrichment Studies for a Diverse Set of Pharmaceutically Relevant Targets , 2007, J. Chem. Inf. Model..

[35]  Steven L Dixon,et al.  PHASE: A Novel Approach to Pharmacophore Modeling and 3D Database Searching , 2006, Chemical biology & drug design.

[36]  J. Massagué TGF-beta signal transduction. , 1998, Annual review of biochemistry.

[37]  A. Tedgui,et al.  Inhibition of Transforming Growth Factor- Signaling Accelerates Atherosclerosis and Induces an Unstable Plaque Phenotype in Mice , 2001 .

[38]  George W. A. Milne,et al.  National Cancer Institute Drug Information System 3D Database , 1994, J. Chem. Inf. Comput. Sci..

[39]  David A. Evans,et al.  3D QSAR Methods: Phase and Catalyst Compared , 2007, J. Chem. Inf. Model..

[40]  Tao Wang,et al.  Quantitative structure–activity relationship: promising advances in drug discovery platforms , 2015, Expert opinion on drug discovery.

[41]  Woody Sherman,et al.  Protein and ligand preparation: parameters, protocols, and influence on virtual screening enrichments , 2013, Journal of Computer-Aided Molecular Design.

[42]  K. Saravanan,et al.  Comparison of Molecular Docking and Molecular Dynamics Simulations of 1,3-Thiazin-4-One with MDM2 Protein , 2015, International Letters of Chemistry, Physics and Astronomy.

[43]  Xin Wen,et al.  BindingDB: a web-accessible database of experimentally determined protein–ligand binding affinities , 2006, Nucleic Acids Res..

[44]  Ester Piek,et al.  Role of Transforming Growth Factor-β Signaling in Cancer , 2000 .