1
|
Keretsu S, Ghosh S, Cho SJ. Computer aided designing of novel pyrrolopyridine derivatives as JAK1 inhibitors. Sci Rep 2021; 11:23051. [PMID: 34845259 PMCID: PMC8630053 DOI: 10.1038/s41598-021-02364-2] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2021] [Accepted: 11/10/2021] [Indexed: 02/07/2023] Open
Abstract
Janus kinases (JAKs) are a family of non-receptor kinases that play a key role in cytokine signaling and their aberrant activities are associated with the pathogenesis of various immune diseases. The JAK1 isoform plays an essential role in the types 1 and II interferon signaling and elicits signals from the interleukin-2, interleukin-4, gp130, and class 2 receptor families. It is ubiquitously expressed in humans and its overexpression has been linked with autoimmune diseases such as myeloproliferative neoplasm. Although JAK1 inhibitors such as Tofacitinib have been approved for medical use, the low potency and off-target effects of these inhibitors have limited their use and calls for the development of novel JAK1 inhibitors. In this study, we used computational methods on a series of pyrrolopyridine derivatives to design new JAK1 inhibitors. Molecular docking and molecular dynamics simulation methods were used to study the protein-inhibitor interactions. 3D-quantitative structure–activity relationship models were developed and were used to predict the activity of newly designed compounds. Free energy calculation methods were used to study the binding affinity of the inhibitors with JAK1. Of the designed compounds, seventeen of the compounds showed a higher binding energy value than the most active compound in the dataset and at least six of the compounds showed higher binding energy value than the pan JAK inhibitor Tofacitinib. The findings made in this study could be utilized for the further development of JAK1 inhibitors.
Collapse
Affiliation(s)
- Seketoulie Keretsu
- Department of Biomedical Sciences, College of Medicine, Chosun University, 375 Seosuk-dong, Dong-gu, Gwangju, 501-759, Republic of Korea
| | - Suparna Ghosh
- Department of Biomedical Sciences, College of Medicine, Chosun University, 375 Seosuk-dong, Dong-gu, Gwangju, 501-759, Republic of Korea
| | - Seung Joo Cho
- Department of Biomedical Sciences, College of Medicine, Chosun University, 375 Seosuk-dong, Dong-gu, Gwangju, 501-759, Republic of Korea. .,Department of Cellular Molecular Medicine, College of Medicine, Chosun University, Gwangju, 501-759, Republic of Korea.
| |
Collapse
|
2
|
Lin FY, MacKerell AD. Polarizable Empirical Force Field for Halogen-Containing Compounds Based on the Classical Drude Oscillator. J Chem Theory Comput 2018; 14:1083-1098. [PMID: 29357257 PMCID: PMC5811359 DOI: 10.1021/acs.jctc.7b01086] [Citation(s) in RCA: 36] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
The quality of the force field is crucial to ensure the accuracy of simulations used in molecular modeling, including computer-aided drug design (CADD). To perform more accurate modeling and simulations of halogenated molecules, in this study the polarizable force field based on the classical Drude oscillator model was extended to both aliphatic and aromatic systems using halogenated ethane and benzene model compounds for the halogens F, Cl, Br, and I. The force field parameters were optimized targeting quantum mechanical dipole moments, water interactions, and molecular polarizabilities as well as experimental observables, including enthalpies of vaporization, molecular volumes, hydration free energies, and dielectric constants. The developed halogenated polarizable force field is capable of reproducing QM relative energies and geometries of both halogen bonds and halogen-hydrogen bond donor interactions at an unprecedented level due to the inclusion of a virtual particle and anisotropic atomic polarizability on the halogen and, notably, the inclusion of Lennard-Jones parameters on the halogen Drude particle. The model was validated on the basis of its ability to accurately reproduce pure solvent properties for halogenated naphthalenes and alkanes, including species analogous to those used as refrigerants. Accordingly, it is anticipated that the model will be applicable for the study of halogenated derivatives in CADD as well as in other chemical and biophysical studies.
Collapse
Affiliation(s)
- Fang-Yu Lin
- Computer-Aided Drug Design Center, Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, MD 21201, USA
| | - Alexander D. MacKerell
- Computer-Aided Drug Design Center, Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, MD 21201, USA
| |
Collapse
|
3
|
Hernández González JE, Hernández Alvarez L, Pascutti PG, Valiente PA. Predicting binding modes of reversible peptide-based inhibitors of falcipain-2 consistent with structure-activity relationships. Proteins 2017; 85:1666-1683. [DOI: 10.1002/prot.25322] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2017] [Revised: 05/05/2017] [Accepted: 05/17/2017] [Indexed: 12/17/2022]
Affiliation(s)
- Jorge Enrique Hernández González
- Departamento de Física; Instituto de Biociências, Letras e Ciências Exatas, Universidade Estadual Paulista Júlio de Mesquita Filho; Rua Cristóvão Colombo, 2265, Jardim Nazareth, São José do Rio Preto São Paulo CEP 15054-000 Brazil
- Centro de Estudios de Proteínas, Facultad de Biología, Universidad de La Habana; Calle 25 No. 455, entre J e I, Vedado, Plaza de la Revolución La Habana CP 10400 Cuba
| | - Lilian Hernández Alvarez
- Departamento de Física; Instituto de Biociências, Letras e Ciências Exatas, Universidade Estadual Paulista Júlio de Mesquita Filho; Rua Cristóvão Colombo, 2265, Jardim Nazareth, São José do Rio Preto São Paulo CEP 15054-000 Brazil
| | - Pedro Geraldo Pascutti
- Laboratório de Dinâmica e Modelagem Molecular, Instituto de Biofísica Carlos Chagas Filho, Universidade Federal do Rio de Janeiro, Ave. Carlos Chagas Filho, 373, CCS-Bloco D sala 30, Cidade Universitária Ilha de Fundão; Rio de Janeiro CEP 21941-902 Brazil
| | - Pedro A. Valiente
- Centro de Estudios de Proteínas, Facultad de Biología, Universidad de La Habana; Calle 25 No. 455, entre J e I, Vedado, Plaza de la Revolución La Habana CP 10400 Cuba
| |
Collapse
|
4
|
Zhang L, Li Y, Yuan Y, Jiang Y, Guo Y, Li M, Pu X. Molecular mechanism of carbon nanotube to activate Subtilisin Carlsberg in polar and non-polar organic media. Sci Rep 2016; 6:36838. [PMID: 27874101 PMCID: PMC5118797 DOI: 10.1038/srep36838] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/23/2016] [Accepted: 10/17/2016] [Indexed: 02/04/2023] Open
Abstract
In the work, we mainly used molecular dynamics (MD) simulation and protein structure network (PSN) to study subtilisin Carlsberg (SC) immobilized onto carbon nanotube (CNT) in water, acetonitrile and heptane solvents, in order to explore activation mechanism of enzymes in non-aqueous media. The result indicates that the affinity of SC with CNT follows the decreasing order of water > acetonitrile > heptane. The overall structure of SC and the catalytic triad display strong robustness to the change of environments, responsible for the activity retaining. However, the distances between two β-strands of substrate-binding pocket are significantly expanded by the immobilization in the increasing order of water < acetonitrile < heptane, contributing to the highest substrate-binding energy in heptane media. PSN analysis further reveals that the immobilization enhances structural communication paths to the substrate-binding pocket, leading to its larger change than the free-enzymes. Interestingly, the increase in the number of the pathways upon immobilization is not dependent on the absorbed extent but the desorbed one, indicating significant role of shifting process of experimental operations in influencing the functional region. In addition, some conserved and important hot-residues in the paths are identified, providing molecular information for functional modification.
Collapse
Affiliation(s)
- Liyun Zhang
- Faculty of Chemistry, Sichuan University, Chengdu 610064, People's Republic of China
| | - Yuzhi Li
- Faculty of Chemistry, Sichuan University, Chengdu 610064, People's Republic of China
| | - Yuan Yuan
- College of Management, Southwest University for Nationalities, Chengdu 610041, People's Republic of China
| | - Yuanyuan Jiang
- Faculty of Chemistry, Sichuan University, Chengdu 610064, People's Republic of China
| | - Yanzhi Guo
- Faculty of Chemistry, Sichuan University, Chengdu 610064, People's Republic of China
| | - Menglong Li
- Faculty of Chemistry, Sichuan University, Chengdu 610064, People's Republic of China
| | - Xuemei Pu
- Faculty of Chemistry, Sichuan University, Chengdu 610064, People's Republic of China
| |
Collapse
|
5
|
Probing immobilization mechanism of alpha-chymotrypsin onto carbon nanotube in organic media by molecular dynamics simulation. Sci Rep 2015; 5:9297. [PMID: 25787884 PMCID: PMC4365409 DOI: 10.1038/srep09297] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/12/2014] [Accepted: 02/26/2015] [Indexed: 11/16/2022] Open
Abstract
The enzyme immobilization has been adopted to enhance the activity and stability of enzymes in non-aqueous enzymatic catalysis. However, the activation and stabilization mechanism has been poorly understood on experiments. Thus, we used molecular dynamics simulation to study the adsorption of α-chymotrypsin (α-ChT) on carbon nanotube (CNT) in aqueous solution and heptane media. The results indicate that α-ChT has stronger affinity with CNT in aqueous solution than in heptane media, as confirmed by more adsorption atoms, larger contact area and higher binding free energies. Although the immobilization causes significant structure deviations from the crystal one, no significant changes in secondary structure of the enzyme upon adsorption are observed in the two media. Different from aqueous solution, the stabilization effects on some local regions far from the surface of CNT were observed in heptane media, in particular for S1 pocket, which should contribute to the preservation of specificity reported by experiments. Also, CNT displays to some extent stabilization role in retaining the catalytic H-bond network of the active site in heptane media, which should be associated with the enhanced activity of enzymes. The observations from the work can provide valuable information for improving the catalytic properties of enzymes in non-aqueous media.
Collapse
|
6
|
Homeyer N, Stoll F, Hillisch A, Gohlke H. Binding Free Energy Calculations for Lead Optimization: Assessment of Their Accuracy in an Industrial Drug Design Context. J Chem Theory Comput 2014; 10:3331-44. [PMID: 26588302 DOI: 10.1021/ct5000296] [Citation(s) in RCA: 119] [Impact Index Per Article: 11.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
Correctly ranking compounds according to their computed relative binding affinities will be of great value for decision making in the lead optimization phase of industrial drug discovery. However, the performance of existing computationally demanding binding free energy calculation methods in this context is largely unknown. We analyzed the performance of the molecular mechanics continuum solvent, the linear interaction energy (LIE), and the thermodynamic integration (TI) approach for three sets of compounds from industrial lead optimization projects. The data sets pose challenges typical for this early stage of drug discovery. None of the methods was sufficiently predictive when applied out of the box without considering these challenges. Detailed investigations of failures revealed critical points that are essential for good binding free energy predictions. When data set-specific features were considered accordingly, predictions valuable for lead optimization could be obtained for all approaches but LIE. Our findings lead to clear recommendations for when to use which of the above approaches. Our findings also stress the important role of expert knowledge in this process, not least for estimating the accuracy of prediction results by TI, using indicators such as the size and chemical structure of exchanged groups and the statistical error in the predictions. Such knowledge will be invaluable when it comes to the question which of the TI results can be trusted for decision making.
Collapse
Affiliation(s)
- Nadine Homeyer
- Institute for Pharmaceutical and Medicinal Chemistry, Department of Mathematics and Natural Sciences, Heinrich-Heine-University Düsseldorf , Universitätsstr.1, 40225 Düsseldorf, Germany
| | - Friederike Stoll
- Global Drug Discovery, Medicinal Chemistry, Bayer Pharma AG , Aprather Weg 18A, 42113 Wuppertal, Germany
| | - Alexander Hillisch
- Global Drug Discovery, Medicinal Chemistry, Bayer Pharma AG , Aprather Weg 18A, 42113 Wuppertal, Germany
| | - Holger Gohlke
- Institute for Pharmaceutical and Medicinal Chemistry, Department of Mathematics and Natural Sciences, Heinrich-Heine-University Düsseldorf , Universitätsstr.1, 40225 Düsseldorf, Germany
| |
Collapse
|
7
|
Slynko I, Scharfe M, Rumpf T, Eib J, Metzger E, Schüle R, Jung M, Sippl W. Virtual screening of PRK1 inhibitors: ensemble docking, rescoring using binding free energy calculation and QSAR model development. J Chem Inf Model 2014; 54:138-50. [PMID: 24377786 DOI: 10.1021/ci400628q] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Abstract
Protein kinase C Related Kinase 1 (PRK1) has been shown to be involved in the regulation of androgen receptor signaling and has been identified as a novel potential drug target for prostate cancer therapy. Since there is no PRK1 crystal structure available to date, multiple PRK1 homology models were generated in order to address the protein flexibility. An in-house library of compounds tested on PRK1 was docked into the ATP binding site of the generated models. In most cases a correct pose of the inhibitors could be identified by ensemble docking, while there is still a challenge of finding a reasonable scoring function that is able to rank compounds according to their biological activity. We estimated the binding free energy for our data set of structurally diverse PRK1 inhibitors using the MM-PB(GB)SA and QM/MM-GBSA methods. The obtained results demonstrate that a correlation between calculated binding free energies and experimental IC50 values was found to be usually higher than using docking scores. Furthermore, the developed approach was tested on a set of diverse PRK1 inhibitors taken from literature, which resulted in a significant correlation. The developed method is computationally inexpensive and can be applied as a postdocking filter in virtual screening as well as for optimization of PRK1 inhibitors in order to prioritize compounds for further biological characterization.
Collapse
Affiliation(s)
- Inna Slynko
- Institute of Pharmacy, MLU Halle-Wittenberg , 06120 Halle (Saale), Germany
| | | | | | | | | | | | | | | |
Collapse
|
8
|
Meng Y, Yuan Y, Zhu Y, Guo Y, Li M, Wang Z, Pu X, Jiang L. Effects of organic solvents and substrate binding on trypsin in acetonitrile and hexane media. J Mol Model 2013; 19:3749-66. [DOI: 10.1007/s00894-013-1900-2] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2013] [Accepted: 05/27/2013] [Indexed: 11/29/2022]
|
9
|
Computation of relative binding free energy for an inhibitor and its analogs binding with Erk kinase using thermodynamic integration MD simulation. J Comput Aided Mol Des 2012; 26:1159-69. [DOI: 10.1007/s10822-012-9606-6] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2012] [Accepted: 09/10/2012] [Indexed: 01/17/2023]
|
10
|
Muzzioli E, Del Rio A, Rastelli G. Assessing protein kinase selectivity with molecular dynamics and mm-pbsa binding free energy calculations. Chem Biol Drug Des 2011; 78:252-9. [PMID: 21585710 DOI: 10.1111/j.1747-0285.2011.01140.x] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
An application of molecular dynamics and molecular mechanics Poisson-Boltzmann surface area techniques to the prediction of protein kinase inhibitor selectivity is presented. A highly active and selective ERK2 inhibitor was placed in equivalent orientations in five different protein kinases (SRC, LCK, GSK3, JNK3 and Aurora-A). Binding free energies were then computed with the molecular mechanics Poisson-Boltzmann surface area approach using 15 nanosecond fully solvated molecular dynamics trajectories of the corresponding protein-ligand complexes. The results show correlation with experimentally determined selectivities and provide useful insights into the underlying structural determinants for selectivity.
Collapse
Affiliation(s)
- Elena Muzzioli
- Dipartimento di Scienze Farmaceutiche, Università di Modena e Reggio Emilia, Via Campi 183, 41125 Modena, Italy
| | | | | |
Collapse
|
11
|
A computational workflow for the design of irreversible inhibitors of protein kinases. J Comput Aided Mol Des 2010; 24:183-94. [DOI: 10.1007/s10822-010-9324-x] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2009] [Accepted: 03/05/2010] [Indexed: 01/09/2023]
|
12
|
Del Rio A, Barbosa AJM, Caporuscio F, Mangiatordi GF. CoCoCo: a free suite of multiconformational chemical databases for high-throughput virtual screening purposes. MOLECULAR BIOSYSTEMS 2010; 6:2122-8. [DOI: 10.1039/c0mb00039f] [Citation(s) in RCA: 45] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|