1
|
Wong CF. 15 Years of molecular simulation of drug-binding kinetics. Expert Opin Drug Discov 2023; 18:1333-1348. [PMID: 37789731 PMCID: PMC10926948 DOI: 10.1080/17460441.2023.2264770] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2023] [Accepted: 09/26/2023] [Indexed: 10/05/2023]
Abstract
INTRODUCTION Drug-binding kinetics has been increasingly recognized as an important factor to be considered in drug discovery. Long residence time could prolong the action of some drugs while produce toxicity on others. Early evaluation of the binding kinetics of drug candidates could reduce attrition rate late in the drug discovery process. Computational prediction of drug-binding kinetics is useful as compounds can be evaluated even before they are made. However, simulation of drug-binding kinetics is a challenging problem because of the long-time scale involved. Nevertheless, significant progress has been made. AREAS COVERED This review illustrates the rapid evolution of qualitative to quantitative molecular dynamics-based methods that have been developed over the last 15 years. EXPERT OPINION The development of new methods based on molecular dynamics simulations now enables computation of absolute association/dissociation rate constants. Cheaper methods capable of identifying candidates with fast or slow binding kinetics, or rank-ordering rate constants are also available. Together, these methods have generated useful insights into the molecular mechanisms of drug-binding kinetics, and the design of drug candidates with therapeutically favorable kinetics. Although predicting absolute rate constants is still expensive and challenging, rapid improvement is expected in the coming years with the continuing refinement of current technologies, development of new methodologies, and the utilization of machine learning.
Collapse
Affiliation(s)
- Chung F Wong
- Department of Chemistry and Biochemistry, University of Missouri-St. Louis, St. Louis, MO, USA
| |
Collapse
|
2
|
Zhou F, Yin S, Xiao Y, Lin Z, Fu W, Zhang YJ. Structure-Kinetic Relationship for Drug Design Revealed by a PLS Model with Retrosynthesis-Based Pre-Trained Molecular Representation and Molecular Dynamics Simulation. ACS OMEGA 2023; 8:18312-18322. [PMID: 37251166 PMCID: PMC10210189 DOI: 10.1021/acsomega.3c02294] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 04/05/2023] [Accepted: 04/26/2023] [Indexed: 05/31/2023]
Abstract
Drug design based on kinetic properties is growing in application. Here, we applied retrosynthesis-based pre-trained molecular representation (RPM) in machine learning (ML) to train 501 inhibitors of 55 proteins and successfully predicted the dissociation rate constant (koff) values of 38 inhibitors from an independent dataset for the N-terminal domain of heat shock protein 90α (N-HSP90). Our RPM molecular representation outperforms other pre-trained molecular representations such as GEM, MPG, and general molecular descriptors from RDKit. Furthermore, we optimized the accelerated molecular dynamics to calculate the relative retention time (RT) for the 128 inhibitors of N-HSP90 and obtained the protein-ligand interaction fingerprints (IFPs) on their dissociation pathways and their influencing weights on the koff value. We observed a high correlation among the simulated, predicted, and experimental -log(koff) values. Combining ML, molecular dynamics (MD) simulation, and IFPs derived from accelerated MD helps design a drug for specific kinetic properties and selectivity profiles to the target of interest. To further validate our koff predictive ML model, we tested our model on two new N-HSP90 inhibitors, which have experimental koff values and are not in our ML training dataset. The predicted koff values are consistent with experimental data, and the mechanism of their kinetic properties can be explained by IFPs, which shed light on the nature of their selectivity against N-HSP90 protein. We believe that the ML model described here is transferable to predict koff of other proteins and will enhance the kinetics-based drug design endeavor.
Collapse
|
3
|
Wang J, Do HN, Koirala K, Miao Y. Predicting Biomolecular Binding Kinetics: A Review. J Chem Theory Comput 2023; 19:2135-2148. [PMID: 36989090 DOI: 10.1021/acs.jctc.2c01085] [Citation(s) in RCA: 8] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/30/2023]
Abstract
Biomolecular binding kinetics including the association (kon) and dissociation (koff) rates are critical parameters for therapeutic design of small-molecule drugs, peptides, and antibodies. Notably, the drug molecule residence time or dissociation rate has been shown to correlate with their efficacies better than binding affinities. A wide range of modeling approaches including quantitative structure-kinetic relationship models, Molecular Dynamics simulations, enhanced sampling, and Machine Learning has been developed to explore biomolecular binding and dissociation mechanisms and predict binding kinetic rates. Here, we review recent advances in computational modeling of biomolecular binding kinetics, with an outlook for future improvements.
Collapse
Affiliation(s)
- Jinan Wang
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, United States
| | - Hung N Do
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, United States
| | - Kushal Koirala
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, United States
| | - Yinglong Miao
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, United States
| |
Collapse
|
4
|
Stegemann S, Moreton C, Svanbäck S, Box K, Motte G, Paudel A. Trends in oral small-molecule drug discovery and product development based on product launches before and after the Rule of Five. Drug Discov Today 2023; 28:103344. [PMID: 36442594 DOI: 10.1016/j.drudis.2022.103344] [Citation(s) in RCA: 19] [Impact Index Per Article: 19.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2022] [Revised: 07/28/2022] [Accepted: 09/01/2022] [Indexed: 11/26/2022]
Abstract
In 1997, the 'Rule of Five' (Ro5) suggested physicochemical limitations for orally administered drugs, based on the analysis of chemical libraries from the early 1990s. In this review, we report on the trends in oral drug product development by analyzing products launched between 1994 and 1997 and between 2013 and 2019. Our analysis confirmed that most new oral drugs are within the Ro5 descriptors; however, the number of new drug products of drugs with molecular weight (MW) and calculated partition coefficient (clogP) beyond the Ro5 has slightly increased. Analysis revealed that there is no single scientific or technological reason for this trend, but that it likely results from incremental advances are being made in molecular biology, target diversity, drug design, medicinal chemistry, predictive modeling, drug metabolism and pharmacokinetics, and drug delivery.
Collapse
Affiliation(s)
- Sven Stegemann
- Institute for Process and Particle Engineering, Graz University of Technology, Inffeldgasse 13, 8010 Graz, Austria.
| | | | - Sami Svanbäck
- The Solubility Company Ltd, Viikinkaari 4, 00790 Helsinki, Finland
| | - Karl Box
- Pion Inc. (UK) Ltd, Forest Row, UK
| | - Geneviève Motte
- JEN Pharma Consulting, 182 Rue Henri Latour, 1450 Chastre, Belgium
| | - Amrit Paudel
- Institute for Process and Particle Engineering, Graz University of Technology, Inffeldgasse 13, 8010 Graz, Austria; Research Center Pharmaceutical Engineering GmbH, Inffeldgasse 13, 8010 Graz, Austria
| |
Collapse
|
5
|
Zhou Y, Wu J, Xue G, Li J, Jiang L, Huang M. Structural study of the uPA-nafamostat complex reveals a covalent inhibitory mechanism of nafamostat. Biophys J 2022; 121:3940-3949. [PMID: 36039386 PMCID: PMC9674978 DOI: 10.1016/j.bpj.2022.08.034] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/04/2022] [Revised: 08/02/2022] [Accepted: 08/24/2022] [Indexed: 11/23/2022] Open
Abstract
Nafamostat mesylate (NM) is a synthetic compound that inhibits various serine proteases produced during the coagulation cascade and inflammation. Previous studies showed that NM was a highly safe drug for the treatment of different cancers, but the precise functions and mechanisms of NM are not clear. In this study, we determined a series of crystal structures of NM and its hydrolysates in complex with a serine protease (urokinase-type plasminogen activator [uPA]). These structures reveal that NM was cleaved by uPA and that a hydrolyzed product (4-guanidinobenzoic acid [GBA]) remained covalently linked to Ser195 of uPA, and the other hydrolyzed product (6-amidino-2-naphthol [6A2N]) released from uPA. Strikingly, in the inactive uPA (uPA-S195A):NM structure, the 6A2N side of intact NM binds to the specific pocket of uPA. Molecular dynamics simulations and end-point binding free-energy calculations show that the conf1 of NM (6A2N as P1 group) in the uPA-S195A:NM complex may be more stable than conf2 of NM (GBA as P1 group). Moreover, in the structure of uPA:NM complex, the imidazole group of His57 flips further away from Ser195 and disrupts the stable canonical catalytic triad conformation. These results not only reveal the inhibitory mechanism of NM as an efficient serine protease inhibitor but also might provide the structural basis for the further development of serine protease inhibitors.
Collapse
Affiliation(s)
- Yang Zhou
- College of Chemistry, Fuzhou University, Fuzhou, Fujian, P.R. China
| | - Juhong Wu
- College of Chemistry, Fuzhou University, Fuzhou, Fujian, P.R. China
| | - Guangpu Xue
- College of Chemistry, Fuzhou University, Fuzhou, Fujian, P.R. China
| | - Jinyu Li
- College of Chemistry, Fuzhou University, Fuzhou, Fujian, P.R. China
| | - Longguang Jiang
- College of Chemistry, Fuzhou University, Fuzhou, Fujian, P.R. China; Fujian Key Laboratory of Marine Enzyme Engineering, Fuzhou University, Fuzhou, Fujian, P.R. China.
| | - Mingdong Huang
- College of Chemistry, Fuzhou University, Fuzhou, Fujian, P.R. China.
| |
Collapse
|
6
|
Chang YF, Chen SY, Lee CC, Chen J, Lai CS. Easy and Rapid Approach to Obtaining the Binding Affinity of Biomolecular Interactions Based on the Deep Learning Boost. Anal Chem 2022; 94:10427-10434. [PMID: 35837692 DOI: 10.1021/acs.analchem.2c01620] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Recently, the deep learning (DL) dimension of artificial intelligence has received much attention from biochemical researchers and thus has gradually become the key approach adopted in the area of biosensing applications. Studies have shown that the use of DL techniques for sensing can not only shorten the time of data analysis but also significantly increase the accuracy of data analysis and prediction, resulting in the performance improvement of biosensing systems in comparison to conventional methods. However, obtaining reliable equilibrium and rate constants of biomolecular interactions during the detection process remains difficult and time-consuming to date. In this study, we propose a transformed model based on the deep transfer learning and sequence-to-sequence autoencoder that can successfully transfer the SPR sensorgram to the protein-binding constants, that is, the association rate constant (ka) and dissociation rate constant (kd), which provide crucial information to understand the mechanisms of drug action and the functional structures of biomolecules. Experimentally, we first trained and tested the pre-trained model using the Langmuir model which generated ideal SPR sensorgrams and then we fine-tuned the pre-trained model through the augmented SPR sensorgrams which were synthesized by using the synthesized minority oversampling technique (SMOTE) through the moderate-scale experiment. Next, the fine-tuned model was inputted with a short experimental SPR sensorgram that only needs 110 s, and the sensorgram was directly transformed into a reconstructed ideal sensorgram. Finally, the binding kinetic constants, that is, ka and kd, as outputs, were obtained through fitting the reconstructed ideal sensorgram. The results showed that the prediction errors of ka and kd obtained by our model were less than 12 and 24%, respectively. Based on the convenience, accuracy, and reliability of the proposed DL approach, we believe our strategy significantly boosts the feasibility to monitor the binding affinity of antibodies online during production.
Collapse
Affiliation(s)
- Ying-Feng Chang
- Artificial Intelligence Research Center, Chang Gung University, Kweishan District, Taoyuan City 33302, Taiwan
| | - Sin-You Chen
- Artificial Intelligence Research Center, Chang Gung University, Kweishan District, Taoyuan City 33302, Taiwan
| | - Chi-Ching Lee
- Artificial Intelligence Research Center, Chang Gung University, Kweishan District, Taoyuan City 33302, Taiwan.,Department of Computer Science and Information Engineering, Chang Gung University, Kweishan District, Taoyuan City 33302, Taiwan.,Genomic Medicine Core Laboratory, Chang Gung Memorial Hospital, Kweishan District, Taoyuan City 33305, Taiwan
| | - Jenhui Chen
- Artificial Intelligence Research Center, Chang Gung University, Kweishan District, Taoyuan City 33302, Taiwan.,Department of Computer Science and Information Engineering, Chang Gung University, Kweishan District, Taoyuan City 33302, Taiwan.,Division of Breast Surgery and General Surgery, Department of Surgery, Chang Gung Memorial Hospital, Linkou Branch, Kweishan District, Taoyuan City 33305, Taiwan.,Department of Electronic Engineering, Ming Chi University of Technology, Taishan District, New Taipei City 24301, Taiwan
| | - Chao-Sung Lai
- Artificial Intelligence Research Center, Chang Gung University, Kweishan District, Taoyuan City 33302, Taiwan.,Department of Electronic Engineering, Chang Gung University, Kweishan District, Taoyuan City 33302, Taiwan.,Center for Biomedical Engineering, Chang Gung University, Kweishan District, Taoyuan City 33302, Taiwan.,Department of Nephrology, Chang Gung Memorial Hospital, Kweishan District, Taoyuan City 33305, Taiwan.,Department of Materials Engineering, Ming Chi University of Technology, Taishan District, New Taipei City 24301, Taiwan
| |
Collapse
|
7
|
Penna E, Niso M, Podlewska S, Volpicelli F, Crispino M, Perrone-Capano C, Bojarski AJ, Lacivita E, Leopoldo M. In Vitro and In Silico Analysis of the Residence Time of Serotonin 5-HT 7 Receptor Ligands with Arylpiperazine Structure: A Structure-Kinetics Relationship Study. ACS Chem Neurosci 2022; 13:497-509. [PMID: 35099177 DOI: 10.1021/acschemneuro.1c00710] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
During the last decade, the kinetics of drug-target interaction has received increasing attention as an important pharmacological parameter in the drug development process. Several studies have suggested that the lipophilicity of a molecule can play an important role. To date, this aspect has been studied for several G protein-coupled receptors (GPCRs) ligands but not for the 5-HT7 receptor (5-HT7R), a GPCR proposed as a valid therapeutic target in neurodevelopmental and neuropsychiatric disorders associated with abnormal neuronal connectivity. In this study, we report on structure-kinetics relationships of a set of arylpiperazine-based 5-HT7R ligands. We found that it is not the overall lipophilicity of the molecule that influences drug-target interaction kinetics but rather the position of polar groups within the molecule. Next, we performed a combination of molecular docking studies and molecular dynamics simulations to gain insights into structure-kinetics relationships. These studies did not suggest specific contact patterns between the ligands and the receptor-binding site as determinants for compounds kinetics. Finally, we compared the abilities of two 5-HT7R agonists with similar receptor-binding affinities and different residence times to stimulate the 5-HT7R-mediated neurite outgrowth in mouse neuronal primary cultures and found that the compounds induced the effect with different timing. This study provides the first insights into the binding kinetics of arylpiperazine-based 5-HT7R ligands that can be helpful to design new 5-HT7R ligands with fine-tuning of the kinetic profile.
Collapse
Affiliation(s)
- Eduardo Penna
- Department of Biology, University of Naples Federico II, via Cintia 26, 80126 Naples, Italy
- Biofordrug srl, via Dante 99, 70019 Triggiano (Bari), Italy
| | - Mauro Niso
- Dipartimento di Farmacia-Scienze del Farmaco, Università degli Studi di Bari Aldo Moro, via Orabona 4, 70125 Bari, Italy
| | - Sabina Podlewska
- Maj Institute of Pharmacology, Polish Academy of Sciences, 12 Smętna Street, 31-343 Kraków, Poland
| | - Floriana Volpicelli
- Department of Pharmacy, School of Medicine and Surgery, University of Naples Federico II, via Domenico Montesano 49, 80131 Naples, Italy
| | - Marianna Crispino
- Department of Biology, University of Naples Federico II, via Cintia 26, 80126 Naples, Italy
| | - Carla Perrone-Capano
- Department of Pharmacy, School of Medicine and Surgery, University of Naples Federico II, via Domenico Montesano 49, 80131 Naples, Italy
- Institute of Genetics and Biophysics “Adriano Buzzati Traverso”, National Research Council (CNR), via Pietro Castellino 111, 80131 Naples, Italy
| | - Andrzej J. Bojarski
- Maj Institute of Pharmacology, Polish Academy of Sciences, 12 Smętna Street, 31-343 Kraków, Poland
| | - Enza Lacivita
- Dipartimento di Farmacia-Scienze del Farmaco, Università degli Studi di Bari Aldo Moro, via Orabona 4, 70125 Bari, Italy
| | - Marcello Leopoldo
- Dipartimento di Farmacia-Scienze del Farmaco, Università degli Studi di Bari Aldo Moro, via Orabona 4, 70125 Bari, Italy
| |
Collapse
|
8
|
Zhang Q, Zhao N, Meng X, Yu F, Yao X, Liu H. The prediction of protein-ligand unbinding for modern drug discovery. Expert Opin Drug Discov 2021; 17:191-205. [PMID: 34731059 DOI: 10.1080/17460441.2022.2002298] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Abstract
INTRODUCTION Drug-target thermodynamic and kinetic information have perennially important roles in drug design. The prediction of protein-ligand unbinding, which can provide important kinetic information, in experiments continues to face great challenges. Uncovering protein-ligand unbinding through molecular dynamics simulations has become efficient and inexpensive with the progress and enhancement of computing power and sampling methods. AREAS COVERED In this review, various sampling methods for protein-ligand unbinding and their basic principles are firstly briefly introduced. Then, their applications in predicting aspects of protein-ligand unbinding, including unbinding pathways, dissociation rate constants, residence time and binding affinity, are discussed. EXPERT OPINION Although various sampling methods have been successfully applied in numerous systems, they still have shortcomings and deficiencies. Most enhanced sampling methods require researchers to possess a wealth of prior knowledge of collective variables or reaction coordinates. In addition, most systems studied at present are relatively simple, and the study of complex systems in real drug research remains greatly challenging. Through the combination of machine learning and enhanced sampling methods, prediction accuracy can be further improved, and some problems encountered in complex systems also may be solved.
Collapse
Affiliation(s)
| | - Nannan Zhao
- School of Pharmacy, Lanzhou University, Lanzhou, China
| | - Xiaoxiao Meng
- School of Pharmacy, Lanzhou University, Lanzhou, China
| | - Fansen Yu
- School of Pharmacy, Lanzhou University, Lanzhou, China
| | - Xiaojun Yao
- College of Chemistry and Chemical Engineering, Lanzhou University, Lanzhou, China.,Dr. Neher's Biophysics Laboratory for Innovative Drug Discovery, State Key Laboratory of Quality Research in Chinese Medicine, Macau University of Science and Technology, Macau, China
| | - Huanxiang Liu
- School of Pharmacy, Lanzhou University, Lanzhou, China
| |
Collapse
|
9
|
Decherchi S, Cavalli A. Thermodynamics and Kinetics of Drug-Target Binding by Molecular Simulation. Chem Rev 2020; 120:12788-12833. [PMID: 33006893 PMCID: PMC8011912 DOI: 10.1021/acs.chemrev.0c00534] [Citation(s) in RCA: 130] [Impact Index Per Article: 32.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2020] [Indexed: 12/19/2022]
Abstract
Computational studies play an increasingly important role in chemistry and biophysics, mainly thanks to improvements in hardware and algorithms. In drug discovery and development, computational studies can reduce the costs and risks of bringing a new medicine to market. Computational simulations are mainly used to optimize promising new compounds by estimating their binding affinity to proteins. This is challenging due to the complexity of the simulated system. To assess the present and future value of simulation for drug discovery, we review key applications of advanced methods for sampling complex free-energy landscapes at near nonergodicity conditions and for estimating the rate coefficients of very slow processes of pharmacological interest. We outline the statistical mechanics and computational background behind this research, including methods such as steered molecular dynamics and metadynamics. We review recent applications to pharmacology and drug discovery and discuss possible guidelines for the practitioner. Recent trends in machine learning are also briefly discussed. Thanks to the rapid development of methods for characterizing and quantifying rare events, simulation's role in drug discovery is likely to expand, making it a valuable complement to experimental and clinical approaches.
Collapse
Affiliation(s)
- Sergio Decherchi
- Computational
and Chemical Biology, Fondazione Istituto
Italiano di Tecnologia, 16163 Genoa, Italy
| | - Andrea Cavalli
- Computational
and Chemical Biology, Fondazione Istituto
Italiano di Tecnologia, 16163 Genoa, Italy
- Department
of Pharmacy and Biotechnology, University
of Bologna, 40126 Bologna, Italy
| |
Collapse
|
10
|
Amangeldiuly N, Karlov D, Fedorov MV. Baseline Model for Predicting Protein–Ligand Unbinding Kinetics through Machine Learning. J Chem Inf Model 2020; 60:5946-5956. [DOI: 10.1021/acs.jcim.0c00450] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Affiliation(s)
- Nurlybek Amangeldiuly
- Center for Data-Intensive Science and Engineering, Skolkovo Institute of Science and Technology, Moscow 121205, Russia
| | - Dmitry Karlov
- Center for Data-Intensive Science and Engineering, Skolkovo Institute of Science and Technology, Moscow 121205, Russia
| | - Maxim V. Fedorov
- Center for Data-Intensive Science and Engineering, Skolkovo Institute of Science and Technology, Moscow 121205, Russia
- Department of Physics, Scottish Universities Physics Alliance (SUPA), University of Strathclyde, Glasgow G4 0NG, U.K
| |
Collapse
|
11
|
Huang S, Chen L, Mei H, Zhang D, Shi T, Kuang Z, Heng Y, Xu L, Pan X. In Silico Prediction of the Dissociation Rate Constants of Small Chemical Ligands by 3D-Grid-Based VolSurf Method. Int J Mol Sci 2020; 21:ijms21072456. [PMID: 32252223 PMCID: PMC7177943 DOI: 10.3390/ijms21072456] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/17/2020] [Accepted: 03/30/2020] [Indexed: 12/26/2022] Open
Abstract
Accumulated evidence suggests that binding kinetic properties—especially dissociation rate constant or drug-target residence time—are crucial factors affecting drug potency. However, quantitative prediction of kinetic properties has always been a challenging task in drug discovery. In this study, the VolSurf method was successfully applied to quantitatively predict the koff values of the small ligands of heat shock protein 90α (HSP90α), adenosine receptor (AR) and p38 mitogen-activated protein kinase (p38 MAPK). The results showed that few VolSurf descriptors can efficiently capture the key ligand surface properties related to dissociation rate; the resulting models demonstrated to be extremely simple, robust and predictive in comparison with available prediction methods. Therefore, it can be concluded that the VolSurf-based prediction method can be widely applied in the ligand-receptor binding kinetics and de novo drug design researches.
Collapse
Affiliation(s)
- Shuheng Huang
- Key Laboratory of Biorheological Science and Technology (Ministry of Education), Chongqing University, Chongqing 400044, China; (S.H.); (L.C.)
- College of Bioengineering, Chongqing University, Chongqing 400044, China; (D.Z.); (T.S.); (Z.K.); (Y.H.); (L.X.)
| | - Linxin Chen
- Key Laboratory of Biorheological Science and Technology (Ministry of Education), Chongqing University, Chongqing 400044, China; (S.H.); (L.C.)
- College of Bioengineering, Chongqing University, Chongqing 400044, China; (D.Z.); (T.S.); (Z.K.); (Y.H.); (L.X.)
| | - Hu Mei
- Key Laboratory of Biorheological Science and Technology (Ministry of Education), Chongqing University, Chongqing 400044, China; (S.H.); (L.C.)
- College of Bioengineering, Chongqing University, Chongqing 400044, China; (D.Z.); (T.S.); (Z.K.); (Y.H.); (L.X.)
- Correspondence: (H.M.); (X.P.); Tel.: +86-23-65112677 (H.M.)
| | - Duo Zhang
- College of Bioengineering, Chongqing University, Chongqing 400044, China; (D.Z.); (T.S.); (Z.K.); (Y.H.); (L.X.)
| | - Tingting Shi
- College of Bioengineering, Chongqing University, Chongqing 400044, China; (D.Z.); (T.S.); (Z.K.); (Y.H.); (L.X.)
| | - Zuyin Kuang
- College of Bioengineering, Chongqing University, Chongqing 400044, China; (D.Z.); (T.S.); (Z.K.); (Y.H.); (L.X.)
| | - Yu Heng
- College of Bioengineering, Chongqing University, Chongqing 400044, China; (D.Z.); (T.S.); (Z.K.); (Y.H.); (L.X.)
| | - Lei Xu
- College of Bioengineering, Chongqing University, Chongqing 400044, China; (D.Z.); (T.S.); (Z.K.); (Y.H.); (L.X.)
| | - Xianchao Pan
- College of Bioengineering, Chongqing University, Chongqing 400044, China; (D.Z.); (T.S.); (Z.K.); (Y.H.); (L.X.)
- Department of Medicinal Chemistry, College of Pharmacy, Southwest Medical University, Luzhou 646000, China
- Correspondence: (H.M.); (X.P.); Tel.: +86-23-65112677 (H.M.)
| |
Collapse
|
12
|
Norval LW, Krämer SD, Gao M, Herz T, Li J, Rath C, Wöhrle J, Günther S, Roth G. KOFFI and Anabel 2.0-a new binding kinetics database and its integration in an open-source binding analysis software. DATABASE-THE JOURNAL OF BIOLOGICAL DATABASES AND CURATION 2020; 2019:5585575. [PMID: 31608948 PMCID: PMC6790968 DOI: 10.1093/database/baz101] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 04/25/2019] [Revised: 07/12/2019] [Accepted: 07/23/2019] [Indexed: 12/31/2022]
Abstract
The kinetics of featured interactions (KOFFI) database is a novel tool and resource for binding kinetics data from biomolecular interactions. While binding kinetics data are abundant in literature, finding valuable information is a laborious task. We used text extraction methods to store binding rates (association, dissociation) as well as corresponding meta-information (e.g. methods, devices) in a novel database. To date, over 270 articles were manually curated and binding data on over 1705 interactions was collected and stored in the (KOFFI) database. Moreover, the KOFFI database application programming interface was implemented in Anabel (open-source software for the analysis of binding interactions), enabling users to directly compare their own binding data analyses with related experiments described in the database.
Collapse
Affiliation(s)
- Leo William Norval
- ZBSA Center for Biological Systems Analysis, Albert-Ludwigs-University Freiburg, Habsburgerstrasse 49, D-79104 Freiburg, Germany.,Institute of Pharmaceutical Sciences, Pharmaceutical Bioinformatics, Albert-Ludwigs-University Freiburg, Hermann-Herder-Straße 9, D-79104 Freiburg, Germany
| | - Stefan Daniel Krämer
- ZBSA Center for Biological Systems Analysis, Albert-Ludwigs-University Freiburg, Habsburgerstrasse 49, D-79104 Freiburg, Germany.,Faculty for Biology, Albert-Ludwigs-University Freiburg, Schaenzlestrasse 1, D-79104 Freiburg, Germany
| | - Mingjie Gao
- Institute of Pharmaceutical Sciences, Pharmaceutical Bioinformatics, Albert-Ludwigs-University Freiburg, Hermann-Herder-Straße 9, D-79104 Freiburg, Germany
| | - Tobias Herz
- ZBSA Center for Biological Systems Analysis, Albert-Ludwigs-University Freiburg, Habsburgerstrasse 49, D-79104 Freiburg, Germany.,Faculty for Biology, Albert-Ludwigs-University Freiburg, Schaenzlestrasse 1, D-79104 Freiburg, Germany
| | - Jianyu Li
- Institute of Pharmaceutical Sciences, Pharmaceutical Bioinformatics, Albert-Ludwigs-University Freiburg, Hermann-Herder-Straße 9, D-79104 Freiburg, Germany
| | - Christin Rath
- ZBSA Center for Biological Systems Analysis, Albert-Ludwigs-University Freiburg, Habsburgerstrasse 49, D-79104 Freiburg, Germany.,Faculty for Biology, Albert-Ludwigs-University Freiburg, Schaenzlestrasse 1, D-79104 Freiburg, Germany.,BioCopy GmbH, Spechtweg 25, D-79110 Freiburg, Germany.,BIOSS Center for Biological Signalling Studies, Albert-Ludwigs-University Freiburg, Schänzlestrasse 18, D-79104 Freiburg, Germany
| | - Johannes Wöhrle
- ZBSA Center for Biological Systems Analysis, Albert-Ludwigs-University Freiburg, Habsburgerstrasse 49, D-79104 Freiburg, Germany.,IMTEK Department of Microsystems Engineering, Albert-Ludwigs-University Freiburg, Georges-Köhler-Allee 103, D-79110 Freiburg, Germany
| | - Stefan Günther
- Institute of Pharmaceutical Sciences, Pharmaceutical Bioinformatics, Albert-Ludwigs-University Freiburg, Hermann-Herder-Straße 9, D-79104 Freiburg, Germany
| | - Günter Roth
- ZBSA Center for Biological Systems Analysis, Albert-Ludwigs-University Freiburg, Habsburgerstrasse 49, D-79104 Freiburg, Germany.,Faculty for Biology, Albert-Ludwigs-University Freiburg, Schaenzlestrasse 1, D-79104 Freiburg, Germany.,BioCopy GmbH, Spechtweg 25, D-79110 Freiburg, Germany.,BIOSS Center for Biological Signalling Studies, Albert-Ludwigs-University Freiburg, Schänzlestrasse 18, D-79104 Freiburg, Germany
| |
Collapse
|
13
|
Kokh DB, Kaufmann T, Kister B, Wade RC. Machine Learning Analysis of τRAMD Trajectories to Decipher Molecular Determinants of Drug-Target Residence Times. Front Mol Biosci 2019; 6:36. [PMID: 31179286 PMCID: PMC6543870 DOI: 10.3389/fmolb.2019.00036] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/27/2019] [Accepted: 05/02/2019] [Indexed: 12/25/2022] Open
Abstract
Drug-target residence times can impact drug efficacy and safety, and are therefore increasingly being considered during lead optimization. For this purpose, computational methods to predict residence times, τ, for drug-like compounds and to derive structure-kinetic relationships are desirable. A challenge for approaches based on molecular dynamics (MD) simulation is the fact that drug residence times are typically orders of magnitude longer than computationally feasible simulation times. Therefore, enhanced sampling methods are required. We recently reported one such approach: the τRAMD procedure for estimating relative residence times by performing a large number of random acceleration MD (RAMD) simulations in which ligand dissociation occurs in times of about a nanosecond due to the application of an additional randomly oriented force to the ligand. The length of the RAMD simulations is used to deduce τ. The RAMD simulations also provide information on ligand egress pathways and dissociation mechanisms. Here, we describe a machine learning approach to systematically analyze protein-ligand binding contacts in the RAMD trajectories in order to derive regression models for estimating τ and to decipher the molecular features leading to longer τ values. We demonstrate that the regression models built on the protein-ligand interaction fingerprints of the dissociation trajectories result in robust estimates of τ for a set of 94 drug-like inhibitors of heat shock protein 90 (HSP90), even for the compounds for which the length of the RAMD trajectories does not provide a good estimation of τ. Thus, we find that machine learning helps to overcome inaccuracies in the modeling of protein-ligand complexes due to incomplete sampling or force field deficiencies. Moreover, the approach facilitates the identification of features important for residence time. In particular, we observed that interactions of the ligand with the sidechain of F138, which is located on the border between the ATP binding pocket and a hydrophobic transient sub-pocket, play a key role in slowing compound dissociation. We expect that the combination of the τRAMD simulation procedure with machine learning analysis will be generally applicable as an aid to target-based lead optimization.
Collapse
Affiliation(s)
- Daria B Kokh
- Molecular and Cellular Modeling Group, Heidelberg Institute for Theoretical Studies (HITS), Heidelberg, Germany
| | - Tom Kaufmann
- Molecular and Cellular Modeling Group, Heidelberg Institute for Theoretical Studies (HITS), Heidelberg, Germany.,Department of Biosciences, Heidelberg University, Heidelberg, Germany
| | - Bastian Kister
- Molecular and Cellular Modeling Group, Heidelberg Institute for Theoretical Studies (HITS), Heidelberg, Germany.,Department of Biosciences, Heidelberg University, Heidelberg, Germany
| | - Rebecca C Wade
- Molecular and Cellular Modeling Group, Heidelberg Institute for Theoretical Studies (HITS), Heidelberg, Germany.,Zentrum für Molekulare Biologie der Universität Heidelberg, DKFZ-ZMBH Alliance, Heidelberg University, Heidelberg, Germany.,Interdisciplinary Center for Scientific Computing (IWR), Heidelberg University, Heidelberg, Germany.,Department of Physics, Heidelberg University, Heidelberg, Germany
| |
Collapse
|
14
|
Zhang D, Huang S, Mei H, Kevin M, Shi T, Chen L. Protein-ligand interaction fingerprints for accurate prediction of dissociation rates of p38 MAPK Type II inhibitors. Integr Biol (Camb) 2019; 11:53-60. [PMID: 30855664 DOI: 10.1093/intbio/zyz004] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2018] [Revised: 11/20/2018] [Accepted: 02/01/2019] [Indexed: 12/22/2022]
Abstract
Binding/unbinding kinetics are key determinants of drug potencies. However, there are still a lot of challenges in predicting kinetic properties during early-stage drug development. In this work, position-restrained molecular dynamics simulations combined with energy decomposition were applied to extract protein-ligand interaction (PLI) fingerprints along the unbinding pathway of 20 p38 mitogen-activated protein kinase (p38 MAPK) Type II inhibitors. The results showed that the electrostatic and/or van der Waals interaction fingerprints at three key positions can be used for accurate prediction of the dissociation rate constants (koff) of p38 MAPK Type II inhibitors. The strategy proposed in this paper can provide not only an efficient method of predicting the dissociation rates of the p38 MAPK Type II inhibitors, but also the atom-level mechanism of enthalpy-driven unbinding process.
Collapse
Affiliation(s)
- Duo Zhang
- Key Laboratory of Biorheological Science and Technology (Ministry of Education), Chongqing University, Chongqing, China
- College of Bioengineering, Chongqing University, Chongqing, China
| | - Shuheng Huang
- Key Laboratory of Biorheological Science and Technology (Ministry of Education), Chongqing University, Chongqing, China
- College of Bioengineering, Chongqing University, Chongqing, China
| | - Hu Mei
- Key Laboratory of Biorheological Science and Technology (Ministry of Education), Chongqing University, Chongqing, China
- College of Bioengineering, Chongqing University, Chongqing, China
| | | | - Tingting Shi
- College of Bioengineering, Chongqing University, Chongqing, China
| | - Linxin Chen
- College of Bioengineering, Chongqing University, Chongqing, China
| |
Collapse
|
15
|
|
16
|
Marques SM, Bednar D, Damborsky J. Computational Study of Protein-Ligand Unbinding for Enzyme Engineering. Front Chem 2019; 6:650. [PMID: 30671430 PMCID: PMC6331733 DOI: 10.3389/fchem.2018.00650] [Citation(s) in RCA: 30] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2018] [Accepted: 12/13/2018] [Indexed: 12/28/2022] Open
Abstract
The computational prediction of unbinding rate constants is presently an emerging topic in drug design. However, the importance of predicting kinetic rates is not restricted to pharmaceutical applications. Many biotechnologically relevant enzymes have their efficiency limited by the binding of the substrates or the release of products. While aiming at improving the ability of our model enzyme haloalkane dehalogenase DhaA to degrade the persistent anthropogenic pollutant 1,2,3-trichloropropane (TCP), the DhaA31 mutant was discovered. This variant had a 32-fold improvement of the catalytic rate toward TCP, but the catalysis became rate-limited by the release of the 2,3-dichloropropan-1-ol (DCP) product from its buried active site. Here we present a computational study to estimate the unbinding rates of the products from DhaA and DhaA31. The metadynamics and adaptive sampling methods were used to predict the relative order of kinetic rates in the different systems, while the absolute values depended significantly on the conditions used (method, force field, and water model). Free energy calculations provided the energetic landscape of the unbinding process. A detailed analysis of the structural and energetic bottlenecks allowed the identification of the residues playing a key role during the release of DCP from DhaA31 via the main access tunnel. Some of these hot-spots could also be identified by the fast CaverDock tool for predicting the transport of ligands through tunnels. Targeting those hot-spots by mutagenesis should improve the unbinding rates of the DCP product and the overall catalytic efficiency with TCP.
Collapse
Affiliation(s)
- Sérgio M. Marques
- Loschmidt Laboratories, Department of Experimental Biology and Research Centre for Toxic Compounds in the Environment RECETOX, Faculty of Science, Masaryk University, Brno, Czechia
- International Clinical Research Center, St. Anne's University Hospital Brno, Brno, Czechia
| | - David Bednar
- Loschmidt Laboratories, Department of Experimental Biology and Research Centre for Toxic Compounds in the Environment RECETOX, Faculty of Science, Masaryk University, Brno, Czechia
- International Clinical Research Center, St. Anne's University Hospital Brno, Brno, Czechia
| | - Jiri Damborsky
- Loschmidt Laboratories, Department of Experimental Biology and Research Centre for Toxic Compounds in the Environment RECETOX, Faculty of Science, Masaryk University, Brno, Czechia
- International Clinical Research Center, St. Anne's University Hospital Brno, Brno, Czechia
| |
Collapse
|
17
|
Jagger BR, Lee CT, Amaro RE. Quantitative Ranking of Ligand Binding Kinetics with a Multiscale Milestoning Simulation Approach. J Phys Chem Lett 2018; 9:4941-4948. [PMID: 30070844 PMCID: PMC6443090 DOI: 10.1021/acs.jpclett.8b02047] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/07/2023]
Abstract
Efficient prediction and ranking of small molecule binders by their kinetic ( kon and koff) and thermodynamic ( Δ G) properties can be a valuable metric for drug lead optimization, as these quantities are often indicators of in vivo efficacy. We have previously described a hybrid molecular dynamics, Brownian dynamics, and milestoning model, Simulation Enabled Estimation of Kinetic Rates (SEEKR), that can predict kon's, koff's, and Δ G's. Here we demonstrate the effectiveness of this approach for ranking a series of seven small molecule compounds for the model system, β-cyclodextrin, based on predicted kon's and koff's. We compare our results using SEEKR to experimentally determined rates as well as rates calculated using long time scale molecular dynamics simulations and show that SEEKR can effectively rank the compounds by koff and Δ G with reduced computational cost. We also provide a discussion of convergence properties and sensitivities of calculations with SEEKR to establish "best practices" for its future use.
Collapse
Affiliation(s)
- Benjamin R Jagger
- Department of Chemistry and Biochemistry , University of California, San Diego , 9500 Gilman Drive , La Jolla , California 92093-0340 , United States
| | - Christopher T Lee
- Department of Chemistry and Biochemistry , University of California, San Diego , 9500 Gilman Drive , La Jolla , California 92093-0340 , United States
| | - Rommie E Amaro
- Department of Chemistry and Biochemistry , University of California, San Diego , 9500 Gilman Drive , La Jolla , California 92093-0340 , United States
| |
Collapse
|
18
|
You W, Chang CEA. Role of Molecular Interactions and Protein Rearrangement in the Dissociation Kinetics of p38α MAP Kinase Type-I/II/III Inhibitors. J Chem Inf Model 2018; 58:968-981. [PMID: 29620886 PMCID: PMC5975198 DOI: 10.1021/acs.jcim.7b00640] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Understanding the governing factors of fast or slow inhibitor binding/unbinding assists in developing drugs with preferred kinetic properties. For inhibitors with the same binding affinity targeting different binding sites of the same protein, the kinetic behavior can profoundly differ. In this study, we investigated unbinding kinetics and mechanisms of fast (type-I) and slow (type-II/III) binders of p38α mitogen-activated protein kinase, where the crystal structures showed that type-I and type-II/III inhibitors bind to pockets with different conformations of the Asp-Phe-Gly (DFG) motif. The work used methods that combine conventional molecular dynamics (MD), accelerated molecular dynamics (AMD) simulations, and the newly developed pathway search guided by internal motions (PSIM) method to find dissociation pathways. The study focuses on revealing key interactions and molecular rearrangements that hinder ligand dissociation by using umbrella sampling and post-MD processing to examine changes in free energy during ligand unbinding. As anticipated, the initial dissociation steps all require breaking interactions that appeared in crystal structures of the bound complexes. Interestingly, for type-I inhibitors such as SB2, p38α keeps barrier-free conformational fluctuation in the ligand-bound complex and during ligand dissociation. In contrast, with a type-II/III inhibitor such as BIRB796, with the rearrangements of p38α in its bound state, ligand unbinding features energetically unfavorable protein-ligand concerted movement. Our results also show that the type-II/III inhibitors preferred dissociation pathways through the allosteric channel, which is consistent with an existing publication. The study suggests that the level of required protein rearrangement is one major determining factor of drug binding kinetics in p38α systems, providing useful information for development of inhibitors.
Collapse
Affiliation(s)
- Wanli You
- Department of Chemistry, University of California at Riverside, Riverside, California 92521, United States
| | - Chia-en A. Chang
- Department of Chemistry, University of California at Riverside, Riverside, California 92521, United States
| |
Collapse
|
19
|
Bruce NJ, Ganotra GK, Kokh DB, Sadiq SK, Wade RC. New approaches for computing ligand-receptor binding kinetics. Curr Opin Struct Biol 2017; 49:1-10. [PMID: 29132080 DOI: 10.1016/j.sbi.2017.10.001] [Citation(s) in RCA: 102] [Impact Index Per Article: 14.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2017] [Revised: 09/22/2017] [Accepted: 10/01/2017] [Indexed: 02/08/2023]
Abstract
The recent and growing evidence that the efficacy of a drug can be correlated to target binding kinetics has seeded the development of a multitude of novel methods aimed at computing rate constants for receptor-ligand binding processes, as well as gaining an understanding of the binding and unbinding pathways and the determinants of structure-kinetic relationships. These new approaches include various types of enhanced sampling molecular dynamics simulations and the combination of energy-based models with chemometric analysis. We assess these approaches in the light of the varying levels of complexity of protein-ligand binding processes.
Collapse
Affiliation(s)
- Neil J Bruce
- Molecular and Cellular Modeling Group, Heidelberg Institute for Theoretical Studies, Schloss-Wolfsbrunnenweg 35, 69118 Heidelberg, Germany
| | - Gaurav K Ganotra
- Molecular and Cellular Modeling Group, Heidelberg Institute for Theoretical Studies, Schloss-Wolfsbrunnenweg 35, 69118 Heidelberg, Germany; Interdisciplinary Center for Scientific Computing (IWR), Heidelberg University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Daria B Kokh
- Molecular and Cellular Modeling Group, Heidelberg Institute for Theoretical Studies, Schloss-Wolfsbrunnenweg 35, 69118 Heidelberg, Germany
| | - S Kashif Sadiq
- Molecular and Cellular Modeling Group, Heidelberg Institute for Theoretical Studies, Schloss-Wolfsbrunnenweg 35, 69118 Heidelberg, Germany
| | - Rebecca C Wade
- Molecular and Cellular Modeling Group, Heidelberg Institute for Theoretical Studies, Schloss-Wolfsbrunnenweg 35, 69118 Heidelberg, Germany; Center for Molecular Biology (ZMBH), DKFZ-ZMBH Alliance, Heidelberg University, Im Neuenheimer Feld 282, 69120 Heidelberg, Germany; Interdisciplinary Center for Scientific Computing (IWR), Heidelberg University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany.
| |
Collapse
|
20
|
Nickkholgh B, Sittadjody S, Rothberg MB, Fang X, Li K, Chou JW, Hawkins GA, Balaji K. Beta-catenin represses protein kinase D1 gene expression by non-canonical pathway through MYC/MAX transcription complex in prostate cancer. Oncotarget 2017; 8:78811-78824. [PMID: 29108267 PMCID: PMC5668000 DOI: 10.18632/oncotarget.20229] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2017] [Accepted: 07/09/2017] [Indexed: 12/13/2022] Open
Abstract
Down regulation of Protein Kinase D1 (PrKD1), a novel serine threonine kinase, in prostate, gastric, breast and colon cancers in humans leads to disease progression. While the down regulation of PrKD1 by DNA methylation in gastric cancer and by nuclear beta-catenin in colon cancer has been shown, the regulatory mechanisms in other cancers are unknown. Because we had demonstrated that PrKD1 is the only known kinase to phosphorylate threonine 120 (T120) of beta-catenin in prostate cancer resulting in increased nuclear beta-catenin, we explored the role of beta-catenin in gene regulation of PrKD1. An initial CHIP assay identified potential binding sites for beta-catenin in and downstream of PrKD1 promoter and sequencing confirmed recruitment of beta-catenin to a 166 base pairs sequence upstream of exon 2. Co-transfection studies with PrKD1-promoter-reporter suggested that beta-catenin represses PrKD1 promoter. Efforts to identify transcription factors that mediate the co-repressor effects of beta-catenin identified recruitment of both MYC and its obligate heterodimer MAX to the same binding site as beta-catenin on the PrKD1 promoter site. Moreover, treatment with MYC inhibitor rescued the co-repressor effect of beta-catenin on PrKD1 gene expression. Prostate specific knock out of PrKD1 in transgenic mice demonstrated increased nuclear expression of beta-catenin validating the in vitro studies. Functional studies showed that nuclear translocation of beta-catenin as a consequence of PrKD1 down regulation, increases AR transcriptional activity with attendant downstream effects on androgen responsive genes. In silico human gene expression analysis confirmed the down regulation of PrKD1 in metastatic prostate cancer correlated inversely with the expression of MAX, but not MYC, and positively with MXD1, a competing heterodimer of MAX, suggesting that the dimerization of MAX with either MYC or MXD1 regulates PrKD1 gene expression. The study has identified a novel auto-repressive loop that perpetuates PrKD1 down regulation through beta-catenin/MYC/MAX protein complex.
Collapse
Affiliation(s)
- Bita Nickkholgh
- Department of Physiology-Pharmacology, Wake Forest University School of Medicine, Winston-Salem, NC, USA
- Wake Forest Institute for Regenerative Medicine (WFRM), Wake Forest University School of Medicine, Winston-Salem, NC, USA
| | - Sivanandane Sittadjody
- Wake Forest Institute for Regenerative Medicine (WFRM), Wake Forest University School of Medicine, Winston-Salem, NC, USA
| | | | - Xiaolan Fang
- Wake Forest Institute for Regenerative Medicine (WFRM), Wake Forest University School of Medicine, Winston-Salem, NC, USA
- Department of Cancer Biology, School of Medicine, Wake Forest University, Winston-Salem, NC, USA
- Current/Present address: Clinical Bioinformatics, New York Genome Center, New York, NY, USA
| | - Kunzhao Li
- Biology Department, Wake Forest University, Winston-Salem, NC, USA
| | - Jeff W. Chou
- Department of Biostatistical Sciences, Comprehensive Cancer Center, Wake Forest Baptist Health, Winston-Salem, NC, USA
| | - Gregory A. Hawkins
- Center for Genomics and Personalized Medicine and WFB Comprehensive Cancer Center, Winston-Salem, NC, USA
| | - K.C. Balaji
- Wake Forest Institute for Regenerative Medicine (WFRM), Wake Forest University School of Medicine, Winston-Salem, NC, USA
- Department of Urology, Wake Forest Baptist Health, Winston Salem, NC, USA
- W.G.(Bill) Hefner Veterans Administration Medical Center, Salisbury, NC, USA
- Department of Cancer Biology, School of Medicine, Wake Forest University, Winston-Salem, NC, USA
| |
Collapse
|
21
|
Abstract
Systems pharmacology aims to holistically understand mechanisms of drug actions to support drug discovery and clinical practice. Systems pharmacology modeling (SPM) is data driven. It integrates an exponentially growing amount of data at multiple scales (genetic, molecular, cellular, organismal, and environmental). The goal of SPM is to develop mechanistic or predictive multiscale models that are interpretable and actionable. The current explosions in genomics and other omics data, as well as the tremendous advances in big data technologies, have already enabled biologists to generate novel hypotheses and gain new knowledge through computational models of genome-wide, heterogeneous, and dynamic data sets. More work is needed to interpret and predict a drug response phenotype, which is dependent on many known and unknown factors. To gain a comprehensive understanding of drug actions, SPM requires close collaborations between domain experts from diverse fields and integration of heterogeneous models from biophysics, mathematics, statistics, machine learning, and semantic webs. This creates challenges in model management, model integration, model translation, and knowledge integration. In this review, we discuss several emergent issues in SPM and potential solutions using big data technology and analytics. The concurrent development of high-throughput techniques, cloud computing, data science, and the semantic web will likely allow SPM to be findable, accessible, interoperable, reusable, reliable, interpretable, and actionable.
Collapse
Affiliation(s)
- Lei Xie
- Department of Computer Science, Hunter College, The City University of New York, New York, NY 10065; .,The Graduate Center, The City University of New York, New York, NY 10016
| | - Eli J Draizen
- National Center for Biotechnology Information, National Library of Medicine, National Institutes of Health, Bethesda, Maryland 20894; .,Program in Bioinformatics, Boston University, Boston, Massachusetts 02215
| | - Philip E Bourne
- National Center for Biotechnology Information, National Library of Medicine, National Institutes of Health, Bethesda, Maryland 20894; .,Office of the Director, National Institutes of Health, Bethesda, Maryland 20894
| |
Collapse
|
22
|
Qu S, Huang S, Pan X, Yang L, Mei H. Constructing Interconsistent, Reasonable, and Predictive Models for Both the Kinetic and Thermodynamic Properties of HIV-1 Protease Inhibitors. J Chem Inf Model 2016; 56:2061-2068. [PMID: 27624663 DOI: 10.1021/acs.jcim.6b00326] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
Accumulated evidence suggests that the in vivo biological potency of a ligand is more strongly correlated with the binding/unbinding kinetics than the equilibrium thermodynamics of the protein-ligand interaction (PLI). However, the existing experimental and computational techniques are largely insufficient and limited in large-scale measurements or accurate predictions of the kinetic properties of PLI. In this work, elaborate efforts have been made to develop interconsistent, reasonable, and predictive models of the association rate constant (kon), dissociation rate constant (koff), and equilibrium dissociation constant (KD) of a series of HIV protease inhibitors with different structural skeletons. The results showed that nine Volsurf descriptors derived from water (OH2) and hydrophobic (DRY) probes are key molecular determinants for the kinetic and thermodynamic properties of HIV-1 protease inhibitors. To the best of our knowledge, this is the first time that interconsistent and reasonable models with strong prediction power have been established for both the kinetic and thermodynamic properties of HIV protease inhibitors.
Collapse
Affiliation(s)
- Sujun Qu
- Key Laboratory of Biorheological Science and Technology, Ministry of Education, Chongqing University , Chongqing 400044, China
| | - Shuheng Huang
- College of Bioengineering, Chongqing University , Chongqing 400044, China
| | - Xianchao Pan
- College of Bioengineering, Chongqing University , Chongqing 400044, China
| | - Li Yang
- Key Laboratory of Biorheological Science and Technology, Ministry of Education, Chongqing University , Chongqing 400044, China.,College of Bioengineering, Chongqing University , Chongqing 400044, China
| | - Hu Mei
- Key Laboratory of Biorheological Science and Technology, Ministry of Education, Chongqing University , Chongqing 400044, China.,College of Bioengineering, Chongqing University , Chongqing 400044, China
| |
Collapse
|
23
|
Hart T, Xie L. Providing data science support for systems pharmacology and its implications to drug discovery. Expert Opin Drug Discov 2016; 11:241-56. [PMID: 26689499 DOI: 10.1517/17460441.2016.1135126] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
INTRODUCTION The conventional one-drug-one-target-one-disease drug discovery process has been less successful in tracking multi-genic, multi-faceted complex diseases. Systems pharmacology has emerged as a new discipline to tackle the current challenges in drug discovery. The goal of systems pharmacology is to transform huge, heterogeneous, and dynamic biological and clinical data into interpretable and actionable mechanistic models for decision making in drug discovery and patient treatment. Thus, big data technology and data science will play an essential role in systems pharmacology. AREAS COVERED This paper critically reviews the impact of three fundamental concepts of data science on systems pharmacology: similarity inference, overfitting avoidance, and disentangling causality from correlation. The authors then discuss recent advances and future directions in applying the three concepts of data science to drug discovery, with a focus on proteome-wide context-specific quantitative drug target deconvolution and personalized adverse drug reaction prediction. EXPERT OPINION Data science will facilitate reducing the complexity of systems pharmacology modeling, detecting hidden correlations between complex data sets, and distinguishing causation from correlation. The power of data science can only be fully realized when integrated with mechanism-based multi-scale modeling that explicitly takes into account the hierarchical organization of biological systems from nucleic acid to proteins, to molecular interaction networks, to cells, to tissues, to patients, and to populations.
Collapse
Affiliation(s)
- Thomas Hart
- a The Rockefeller University , New York , NY , USA.,b Department of Biological Sciences, Hunter College , The City University of New York , New York , NY , USA
| | - Lei Xie
- c Department of Computer Science, Hunter College , The City University of New York , New York , NY , USA.,d The Graduate Center , The City University of New York , New York , NY , USA
| |
Collapse
|