1
|
Amezcua M, Setiadi J, Mobley DL. The SAMPL9 host-guest blind challenge: an overview of binding free energy predictive accuracy. Phys Chem Chem Phys 2024; 26:9207-9225. [PMID: 38444308 PMCID: PMC10954238 DOI: 10.1039/d3cp05111k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2023] [Accepted: 02/03/2024] [Indexed: 03/07/2024]
Abstract
We report the results of the SAMPL9 host-guest blind challenge for predicting binding free energies. The challenge focused on macrocycles from pillar[n]-arene and cyclodextrin host families, including WP6, and bCD and HbCD. A variety of methods were used by participants to submit binding free energy predictions. A machine learning approach based on molecular descriptors achieved the highest accuracy (RMSE of 2.04 kcal mol-1) among the ranked methods in the WP6 dataset. Interestingly, predictions for WP6 obtained via docking tended to outperform all methods (RMSE of 1.70 kcal mol-1), most of which are MD based and computationally more expensive. In general, methods applying force fields achieved better correlation with experiments for WP6 opposed to the machine learning and docking models. In the cyclodextrin-phenothiazine challenge, the ATM approach emerged as the top performing method with RMSE less than 1.86 kcal mol-1. Correlation metrics of ranked methods in this dataset were relatively poor compared to WP6. We also highlight several lessons learned to guide future work and help improve studies on the systems discussed. For example, WP6 may be present in other microstates other than its -12 state in the presence of certain guests. Machine learning approaches can be used to fine tune or help train force fields for certain chemistry (i.e. WP6-G4). Certain phenothiazines occupy distinct primary and secondary orientations, some of which were considered individually for accurate binding free energies. The accuracy of predictions from certain methods while starting from a single binding pose/orientation demonstrates the sensitivity of calculated binding free energies to the orientation, and in some cases the likely dominant orientation for the system. Computational and experimental results suggest that guest phenothiazine core traverses both the secondary and primary faces of the cyclodextrin hosts, a bulky cationic side chain will primarily occupy the primary face, and the phenothiazine core substituent resides at the larger secondary face.
Collapse
Affiliation(s)
- Martin Amezcua
- Department of Pharmaceutical Sciences, University of California, Irvine, Irvine, California 92697, USA.
| | - Jeffry Setiadi
- Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California, San Diego, La Jolla, California 92093, USA
| | - David L Mobley
- Department of Pharmaceutical Sciences, University of California, Irvine, Irvine, California 92697, USA.
- Department of Chemistry, University of California, Irvine, Irvine, California 92697, USA
| |
Collapse
|
2
|
Shrivastava A, Magani SKJ, Lokhande KB, Chintakhindi M, Singh A. Exploring the role of TLK2 mutation in tropical calcific pancreatitis: an in silico and molecular dynamics simulation study. J Biomol Struct Dyn 2024:1-20. [PMID: 38500246 DOI: 10.1080/07391102.2024.2329797] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/08/2023] [Accepted: 03/06/2024] [Indexed: 03/20/2024]
Abstract
Tropical calcific pancreatitis (TCP) is a juvenile form of non-alcoholic chronic pancreatitis seen exclusively in tropical countries. The disease poses a high risk of complications, including pancreatic diabetes and cancer, leading to significant mortality due to poor diagnosis and ineffective treatments. This study employed whole exome sequencing (WES) of 5 TCP patient samples to identify genetic variants associated with TCP. Advanced computational techniques were used to gain atomic-level insights into disease progression, including microsecond-scale long MD simulations and essential dynamics. In silico virtual screening was performed to identify potential therapeutic compounds targeting the mutant protein using the Asinex and DrugBank compound library. WES analysis predicted several single nucleotide variants (SNVs) associated with TCP, including a novel missense variant (c.T1802A or p.V601E) in the TLK2 gene. Computational analysis revealed that the p.V601E mutation significantly affected the structure of the TLK2 kinase domain and its conformational dynamics, altering the interaction profile between ATP and the binding pocket. These changes could impact TLK2's kinase activity and functions, potentially correlating with TCP progression. Promising lead compounds that selectively bind to the TLK2 mutant protein were identified, offering potential for therapeutic interventions in TCP. These findings hold great potential for future research.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Ashish Shrivastava
- Translational Bioinformatics and Computational Genomics Research Lab, Department of Life Sciences, Shiv Nadar Institution of Eminence, Gautam Buddha Nagar, UP, India
| | - Sri Krishna Jayadev Magani
- Cancer Biology Lab, Department of Life Sciences, Shiv Nadar Institution of Eminence, Gautam Buddha Nagar, UP, India
| | - Kiran Bharat Lokhande
- Translational Bioinformatics and Computational Genomics Research Lab, Department of Life Sciences, Shiv Nadar Institution of Eminence, Gautam Buddha Nagar, UP, India
| | | | - Ashutosh Singh
- Translational Bioinformatics and Computational Genomics Research Lab, Department of Life Sciences, Shiv Nadar Institution of Eminence, Gautam Buddha Nagar, UP, India
| |
Collapse
|
3
|
Zhang X, Guseinov AA, Jenkins L, Li K, Tikhonova IG, Milligan G, Zhang C. Structural basis for the ligand recognition and signaling of free fatty acid receptors. SCIENCE ADVANCES 2024; 10:eadj2384. [PMID: 38198545 PMCID: PMC10780892 DOI: 10.1126/sciadv.adj2384] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/14/2023] [Accepted: 12/05/2023] [Indexed: 01/12/2024]
Abstract
Free fatty acid receptors 1 to 4 (FFA1 to FFA4) are class A G protein-coupled receptors (GPCRs). FFA1 to FFA3 share substantial sequence similarity, whereas FFA4 is unrelated. However, FFA1 and FFA4 are activated by long-chain fatty acids, while FFA2 and FFA3 respond to short-chain fatty acids generated by intestinal microbiota. FFA1, FFA2, and FFA4 are potential drug targets for metabolic and inflammatory conditions. Here, we determined the active structures of FFA1 and FFA4 bound to docosahexaenoic acid, FFA4 bound to the synthetic agonist TUG-891, and butyrate-bound FFA2, each complexed with an engineered heterotrimeric Gq protein (miniGq), by cryo-electron microscopy. Together with computational simulations and mutagenesis studies, we elucidated the similarities and differences in the binding modes of fatty acid ligands to their respective GPCRs. Our findings unveiled distinct mechanisms of receptor activation and G protein coupling. We anticipate that these outcomes will facilitate structure-based drug development and underpin future research on this group of GPCRs.
Collapse
Affiliation(s)
- Xuan Zhang
- Department of Pharmacology and Chemical Biology, School of Medicine, University of Pittsburgh, Pittsburgh, PA 15261, USA
| | - Abdul-Akim Guseinov
- School of Pharmacy, Medical Biology Centre, Queen’s University Belfast, Belfast BT9 7BL, Northern Ireland, UK
| | - Laura Jenkins
- Centre for Translational Pharmacology, School of Molecular Biosciences, College of Medical, Veterinary and Life Sciences, University of Glasgow, Glasgow G12 8QQ, Scotland, UK
| | - Kunpeng Li
- Cryo-EM Core Facility, Case Western Reserve University, Cleveland, OH 44106, USA
| | - Irina G. Tikhonova
- School of Pharmacy, Medical Biology Centre, Queen’s University Belfast, Belfast BT9 7BL, Northern Ireland, UK
| | - Graeme Milligan
- Centre for Translational Pharmacology, School of Molecular Biosciences, College of Medical, Veterinary and Life Sciences, University of Glasgow, Glasgow G12 8QQ, Scotland, UK
| | - Cheng Zhang
- Department of Pharmacology and Chemical Biology, School of Medicine, University of Pittsburgh, Pittsburgh, PA 15261, USA
| |
Collapse
|
4
|
Zhang X, Guseinov AA, Jenkins L, Li K, Tikhonova IG, Milligan G, Zhang C. Structural basis for the ligand recognition and signaling of free fatty acid receptors. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.08.20.553924. [PMID: 37662198 PMCID: PMC10473637 DOI: 10.1101/2023.08.20.553924] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 09/05/2023]
Abstract
Free fatty acid receptors 1-4 (FFA1-4) are class A G protein-coupled receptors (GPCRs). FFA1-3 share substantial sequence similarity whereas FFA4 is unrelated. Despite this FFA1 and FFA4 are activated by the same range of long chain fatty acids (LCFAs) whilst FFA2 and FFA3 are instead activated by short chain fatty acids (SCFAs) generated by the intestinal microbiota. Each of FFA1, 2 and 4 are promising targets for novel drug development in metabolic and inflammatory conditions. To gain insights into the basis of ligand interactions with, and molecular mechanisms underlying activation of, FFAs by LCFAs and SCFAs, we determined the active structures of FFA1 and FFA4 bound to the polyunsaturated LCFA docosahexaenoic acid (DHA), FFA4 bound to the synthetic agonist TUG-891, as well as SCFA butyrate-bound FFA2, each complexed with an engineered heterotrimeric Gq protein (miniGq), by cryo-electron microscopy. Together with computational simulations and mutagenesis studies, we elucidated the similarities and differences in the binding modes of fatty acid ligands with varying chain lengths to their respective GPCRs. Our findings unveil distinct mechanisms of receptor activation and G protein coupling. We anticipate that these outcomes will facilitate structure-based drug development and underpin future research to understand allosteric modulation and biased signaling of this group of GPCRs.
Collapse
Affiliation(s)
- Xuan Zhang
- Department of Pharmacology and Chemical Biology, School of Medicine, University of Pittsburgh, Pittsburgh, PA15261, USA
| | - Abdul-Akim Guseinov
- School of Pharmacy, Medical Biology Centre, Queen’s University Belfast, Belfast BT9 7BL, Northern Ireland, United Kingdom
| | - Laura Jenkins
- Centre for Translational Pharmacology, School of Molecular Biosciences, College of Medical, Veterinary and Life Sciences, University of Glasgow, Glasgow G12 8QQ, Scotland, United Kingdom
| | - Kunpeng Li
- Cryo-EM core facility, Case Western Reserve University, OH44106, USA
| | - Irina G. Tikhonova
- School of Pharmacy, Medical Biology Centre, Queen’s University Belfast, Belfast BT9 7BL, Northern Ireland, United Kingdom
| | - Graeme Milligan
- Centre for Translational Pharmacology, School of Molecular Biosciences, College of Medical, Veterinary and Life Sciences, University of Glasgow, Glasgow G12 8QQ, Scotland, United Kingdom
| | - Cheng Zhang
- Department of Pharmacology and Chemical Biology, School of Medicine, University of Pittsburgh, Pittsburgh, PA15261, USA
| |
Collapse
|
5
|
An J, Guan J, Nie Y. Semi-Rational Design of L-Isoleucine Dioxygenase Generated Its Activity for Aromatic Amino Acid Hydroxylation. Molecules 2023; 28:molecules28093750. [PMID: 37175159 PMCID: PMC10180240 DOI: 10.3390/molecules28093750] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/11/2023] [Revised: 04/15/2023] [Accepted: 04/24/2023] [Indexed: 05/15/2023] Open
Abstract
Fe (II)-and 2-ketoglutarate-dependent dioxygenases (Fe (II)/α-KG DOs) have been applied to catalyze hydroxylation of amino acids. However, the Fe (II)/α-KG DOs that have been developed and characterized are not sufficient. L-isoleucine dioxygenase (IDO) is an Fe (II)/α-KG DO that specifically catalyzes the formation of 4-hydroxyisoleucine (4-HIL) from L-isoleucine (L-Ile) and exhibits a substrate specificity toward L-aliphatic amino acids. To expand the substrate spectrum of IDO toward aromatic amino acids, in this study, we analyzed the regularity of the substrate spectrum of IDO using molecular dynamics (MD) simulation and found that the distance between Fe2+, C2 of α-KG and amino acid chain's C4 may be critical for regulating the substrate specificity of the enzyme. The mutation sites (Y143, S153 and R227) were also subjected to single point saturation mutations based on polarity pockets and residue free energy contributions. It was found that Y143D, Y143I and S153A mutants exhibited catalytic L-phenylalanine activity, while Y143I, S153A, S153Q and S153Y exhibited catalytic L-homophenylalanine activity. Consequently, this study extended the substrate spectrum of IDO with aromatic amino acids and enhanced its application property.
Collapse
Affiliation(s)
- Jianhong An
- School of Biotechnology and Key Laboratory of Industrial Biotechnology, Ministry of Education, Jiangnan University, 1800 Lihu Road, Wuxi 214122, China
- International Joint Research Laboratory for Brewing Microbiology and Applied Enzymology, Jiangnan University, 1800 Lihu Road, Wuxi 214122, China
- Eye Hospital, Wenzhou Medical University, 270 Xueyuan Road, Wenzhou 325000, China
| | - Jiaojiao Guan
- School of Biotechnology and Key Laboratory of Industrial Biotechnology, Ministry of Education, Jiangnan University, 1800 Lihu Road, Wuxi 214122, China
| | - Yao Nie
- School of Biotechnology and Key Laboratory of Industrial Biotechnology, Ministry of Education, Jiangnan University, 1800 Lihu Road, Wuxi 214122, China
- International Joint Research Laboratory for Brewing Microbiology and Applied Enzymology, Jiangnan University, 1800 Lihu Road, Wuxi 214122, China
| |
Collapse
|
6
|
Liu X, Zheng L, Cong Y, Gong Z, Yin Z, Zhang JZH, Liu Z, Sun Z. Comprehensive evaluation of end-point free energy techniques in carboxylated-pillar[6]arene host-guest binding: II. regression and dielectric constant. J Comput Aided Mol Des 2022; 36:879-894. [PMID: 36394776 DOI: 10.1007/s10822-022-00487-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/26/2022] [Accepted: 10/29/2022] [Indexed: 11/18/2022]
Abstract
End-point free energy calculations as a powerful tool have been widely applied in protein-ligand and protein-protein interactions. It is often recognized that these end-point techniques serve as an option of intermediate accuracy and computational cost compared with more rigorous statistical mechanic models (e.g., alchemical transformation) and coarser molecular docking. However, it is observed that this intermediate level of accuracy does not hold in relatively simple and prototypical host-guest systems. Specifically, in our previous work investigating a set of carboxylated-pillar[6]arene host-guest complexes, end-point methods provide free energy estimates deviating significantly from the experimental reference, and the rank of binding affinities is also incorrectly computed. These observations suggest the unsuitability and inapplicability of standard end-point free energy techniques in host-guest systems, and alteration and development are required to make them practically usable. In this work, we consider two ways to improve the performance of end-point techniques. The first one is the PBSA_E regression that varies the weights of different free energy terms in the end-point calculation procedure, while the second one is considering the interior dielectric constant as an additional variable in the end-point equation. By detailed investigation of the calculation procedure and the simulation outcome, we prove that these two treatments (i.e., regression and dielectric constant) are manipulating the end-point equation in a somehow similar way, i.e., weakening the electrostatic contribution and strengthening the non-polar terms, although there are still many detailed differences between these two methods. With the trained end-point scheme, the RMSE of the computed affinities is improved from the standard ~ 12 kcal/mol to ~ 2.4 kcal/mol, which is comparable to another altered end-point method (ELIE) trained with system-specific data. By tuning PBSA_E weighting factors with the host-specific data, it is possible to further decrease the prediction error to ~ 2.1 kcal/mol. These observations along with the extremely efficient optimized-structure computation procedure suggest the regression (i.e., PBSA_E as well as its GBSA_E extension) as a practically applicable solution that brings end-point methods back into the library of usable tools for host-guest binding. However, the dielectric-constant-variable scheme cannot effectively minimize the experiment-calculation discrepancy for absolute binding affinities, but is able to improve the calculation of affinity ranks. This phenomenon is somehow different from the protein-ligand case and suggests the difference between host-guest and biomacromolecular (protein-ligand and protein-protein) systems. Therefore, the spectrum of tools usable for protein-ligand complexes could be unsuitable for host-guest binding, and numerical validations are necessary to screen out really workable solutions in these 'prototypical' situations.
Collapse
Affiliation(s)
- Xiao Liu
- School of Mathematics, Physics and Statistics, Shanghai University of Engineering Science, Shanghai, 201620, China.
| | - Lei Zheng
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai, 200062, China
| | - Yalong Cong
- School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, 200062, China
| | - Zhihao Gong
- School of Micro-Nano Electronics, Zhejiang University, Hangzhou, 310027, China.,Hangzhou Global Scientific and Technological Innovation Center, Zhejiang University, Hangzhou, 310027, China
| | - Zhixiang Yin
- School of Mathematics, Physics and Statistics, Shanghai University of Engineering Science, Shanghai, 201620, China
| | - John Z H Zhang
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai, 200062, China. .,School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, 200062, China. .,Shenzhen Institute of Advanced Technology, Chinese Academy of Sciences, Shenzhen, Guangdong, China. .,Department of Chemistry, New York University, NY, NY, 10003, USA.
| | - Zhirong Liu
- College of Chemistry and Molecular Engineering, Peking University, Beijing, 100871, China
| | - Zhaoxi Sun
- College of Chemistry and Molecular Engineering, Peking University, Beijing, 100871, China.
| |
Collapse
|
7
|
Liu X, Zheng L, Qin C, Zhang JZH, Sun Z. Comprehensive evaluation of end-point free energy techniques in carboxylated-pillar[6]arene host-guest binding: I. Standard procedure. J Comput Aided Mol Des 2022; 36:735-752. [PMID: 36136209 DOI: 10.1007/s10822-022-00475-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/06/2022] [Accepted: 09/06/2022] [Indexed: 10/14/2022]
Abstract
Despite the massive application of end-point free energy methods in protein-ligand and protein-protein interactions, computational understandings about their performance in relatively simple and prototypical host-guest systems are limited. In this work, we present a comprehensive benchmark calculation with standard end-point free energy techniques in a recent host-guest dataset containing 13 host-guest pairs involving the carboxylated-pillar[6]arene host. We first assess the charge schemes for solutes by comparing the charge-produced electrostatics with many ab initio references, in order to obtain a preliminary albeit detailed view of the charge quality. Then, we focus on four modelling details of end-point free energy calculations, including the docking procedure for the generation of initial condition, the charge scheme for host and guest molecules, the water model used in explicit-solvent sampling, and the end-point methods for free energy estimation. The binding thermodynamics obtained with different modelling schemes are compared with experimental references, and some practical guidelines on maximizing the performance of end-point methods in practical host-guest systems are summarized. Further, we compare our simulation outcome with predictions in the grand challenge and discuss further developments to improve the prediction quality of end-point free energy methods. Overall, unlike the widely acknowledged applicability in protein-ligand binding, the standard end-point calculations cannot produce useful outcomes in host-guest binding and thus are not recommended unless alterations are performed.
Collapse
Affiliation(s)
- Xiao Liu
- School of Mathematics, Physics and Statistics, Shanghai University of Engineering Science, Shanghai, 201620, China.
| | - Lei Zheng
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai, 200062, China
| | - Chu Qin
- School of Mathematics, Physics and Statistics, Shanghai University of Engineering Science, Shanghai, 201620, China
| | - John Z H Zhang
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai, 200062, China.,School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, 200062, China.,Shenzhen Institute of Advanced Technology, Chinese Academy of Sciences, Shenzhen, Guangdong, China.,Department of Chemistry, New York University, New York, NY, 10003, USA
| | - Zhaoxi Sun
- College of Chemistry and Molecular Engineering, Peking University, Beijing, 100871, China.
| |
Collapse
|
8
|
Affinity prediction using deep learning based on SMILES input for D3R grand challenge 4. J Comput Aided Mol Des 2022; 36:225-235. [PMID: 35314897 DOI: 10.1007/s10822-022-00448-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2021] [Accepted: 03/08/2022] [Indexed: 10/18/2022]
Abstract
Modern molecular docking comprises the prediction of pose and affinity. Prediction of docking poses is required for affinity prediction when three-dimensional coordinates of the ligand have not been provided. However, a large number of feature engineering is required for existing methods. In addition, there is a need for a robust model for the sequential combination of pose and affinity prediction due to the probabilistic deviation of the ligand position issue. We propose a pipeline using a bipartite graph neural network and transfer learning trained on a re-docking dataset. We evaluated our model on the released data from drug design data resource grand challenge 4 (D3R GC4). The two target protein data provided by the challenge have different patterns. The model outperformed the best participant by 9% on the BACE target protein from stage 2. Further, our model showed competitive performance on the CatS target protein.
Collapse
|
9
|
Zhu YX, Sheng YJ, Ma YQ, Ding HM. Assessing the Performance of Screening MM/PBSA in Protein-Ligand Interactions. J Phys Chem B 2022; 126:1700-1708. [PMID: 35188781 DOI: 10.1021/acs.jpcb.1c09424] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
Accurate calculation of the binding free energies between a protein and a ligand is the primary objective of structure-based drug design, but it still remains a challenging problem. In this work, we apply the screening molecular mechanics/Poisson Boltzmann surface area (MM/PBSA) method to calculate the binding affinity of protein-ligand interactions. Our results show that the performance of the screening MM/PBSA is better than that of the standard MM/PBSA, especially in a charged-ligand system. In addition, we also investigate the effect of the solute dielectric constant on the results, and find that the optimal solute dielectric constants are different between the neutral-ligand system and the charged-ligand system. Moreover, we also evaluate the effect of the atomic-charge methods on the performance of the screening MM/PBSA. The present study demonstrates that the screening MM/PBSA should be a reliable method for calculating binding energy of biosystems.
Collapse
Affiliation(s)
- Yu-Xin Zhu
- Center for Soft Condensed Matter Physics and Interdisciplinary Research, School of Physical Science and Technology, Soochow University, Suzhou 215006, China
| | - Yan-Jing Sheng
- Center for Soft Condensed Matter Physics and Interdisciplinary Research, School of Physical Science and Technology, Soochow University, Suzhou 215006, China
| | - Yu-Qiang Ma
- National Laboratory of Solid State Microstructures and Department of Physics, Collaborative Innovation Center of Advanced Microstructures, Nanjing University, Nanjing 210093, China
| | - Hong-Ming Ding
- Center for Soft Condensed Matter Physics and Interdisciplinary Research, School of Physical Science and Technology, Soochow University, Suzhou 215006, China
| |
Collapse
|
10
|
King E, Aitchison E, Li H, Luo R. Recent Developments in Free Energy Calculations for Drug Discovery. Front Mol Biosci 2021; 8:712085. [PMID: 34458321 PMCID: PMC8387144 DOI: 10.3389/fmolb.2021.712085] [Citation(s) in RCA: 49] [Impact Index Per Article: 16.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2021] [Accepted: 07/27/2021] [Indexed: 01/11/2023] Open
Abstract
The grand challenge in structure-based drug design is achieving accurate prediction of binding free energies. Molecular dynamics (MD) simulations enable modeling of conformational changes critical to the binding process, leading to calculation of thermodynamic quantities involved in estimation of binding affinities. With recent advancements in computing capability and predictive accuracy, MD based virtual screening has progressed from the domain of theoretical attempts to real application in drug development. Approaches including the Molecular Mechanics Poisson Boltzmann Surface Area (MM-PBSA), Linear Interaction Energy (LIE), and alchemical methods have been broadly applied to model molecular recognition for drug discovery and lead optimization. Here we review the varied methodology of these approaches, developments enhancing simulation efficiency and reliability, remaining challenges hindering predictive performance, and applications to problems in the fields of medicine and biochemistry.
Collapse
Affiliation(s)
- Edward King
- Department of Molecular Biology and Biochemistry, University of California, Irvine, CA, United States
| | - Erick Aitchison
- Department of Molecular Biology and Biochemistry, University of California, Irvine, CA, United States
| | - Han Li
- Department of Chemical and Biomolecular Engineering, University of California, Irvine, CA, United States
| | - Ray Luo
- Department of Molecular Biology and Biochemistry, University of California, Irvine, CA, United States
- Department of Chemical and Biomolecular Engineering, University of California, Irvine, CA, United States
- Department of Materials Science and Engineering, University of California, Irvine, CA, United States
- Department of Biomedical Engineering, University of California, Irvine, CA, United States
| |
Collapse
|
11
|
Ji B, He X, Zhai J, Zhang Y, Man VH, Wang J. Machine learning on ligand-residue interaction profiles to significantly improve binding affinity prediction. Brief Bioinform 2021; 22:6184410. [PMID: 33758923 DOI: 10.1093/bib/bbab054] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/08/2020] [Revised: 01/06/2021] [Accepted: 02/02/2021] [Indexed: 01/01/2023] Open
Abstract
Structure-based virtual screenings (SBVSs) play an important role in drug discovery projects. However, it is still a challenge to accurately predict the binding affinity of an arbitrary molecule binds to a drug target and prioritize top ligands from an SBVS. In this study, we developed a novel method, using ligand-residue interaction profiles (IPs) to construct machine learning (ML)-based prediction models, to significantly improve the screening performance in SBVSs. Such a kind of the prediction model is called an IP scoring function (IP-SF). We systematically investigated how to improve the performance of IP-SFs from many perspectives, including the sampling methods before interaction energy calculation and different ML algorithms. Using six drug targets with each having hundreds of known ligands, we conducted a critical evaluation on the developed IP-SFs. The IP-SFs employing a gradient boosting decision tree (GBDT) algorithm in conjunction with the MIN + GB simulation protocol achieved the best overall performance. Its scoring power, ranking power and screening power significantly outperformed the Glide SF. First, compared with Glide, the average values of mean absolute error and root mean square error of GBDT/MIN + GB decreased about 38 and 36%, respectively. Second, the mean values of squared correlation coefficient and predictive index increased about 225 and 73%, respectively. Third, more encouragingly, the average value of the areas under the curve of receiver operating characteristic for six targets by GBDT, 0.87, is significantly better than that by Glide, which is only 0.71. Thus, we expected IP-SFs to have broad and promising applications in SBVSs.
Collapse
Affiliation(s)
- Beihong Ji
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, PA 15261, USA
| | - Xibing He
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, PA 15261, USA
| | - Jingchen Zhai
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, PA 15261, USA
| | - Yuzhao Zhang
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, PA 15261, USA
| | - Viet Hoang Man
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, PA 15261, USA
| | - Junmei Wang
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, PA 15261, USA
| |
Collapse
|
12
|
Ji B, He X, Zhang Y, Zhai J, Man VH, Liu S, Wang J. Incorporating structural similarity into a scoring function to enhance the prediction of binding affinities. J Cheminform 2021; 13:11. [PMID: 33588902 PMCID: PMC7884591 DOI: 10.1186/s13321-021-00493-4] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2020] [Accepted: 02/03/2021] [Indexed: 11/23/2022] Open
Abstract
In this study, we developed a novel algorithm to improve the screening performance of an arbitrary docking scoring function by recalibrating the docking score of a query compound based on its structure similarity with a set of training compounds, while the extra computational cost is neglectable. Two popular docking methods, Glide and AutoDock Vina were adopted as the original scoring functions to be processed with our new algorithm and similar improvement performance was achieved. Predicted binding affinities were compared against experimental data from ChEMBL and DUD-E databases. 11 representative drug receptors from diverse drug target categories were applied to evaluate the hybrid scoring function. The effects of four different fingerprints (FP2, FP3, FP4, and MACCS) and the four different compound similarity effect (CSE) functions were explored. Encouragingly, the screening performance was significantly improved for all 11 drug targets especially when CSE = S4 (S is the Tanimoto structural similarity) and FP2 fingerprint were applied. The average predictive index (PI) values increased from 0.34 to 0.66 and 0.39 to 0.71 for the Glide and AutoDock vina scoring functions, respectively. To evaluate the performance of the calibration algorithm in drug lead identification, we also imposed an upper limit on the structural similarity to mimic the real scenario of screening diverse libraries for which query ligands are general-purpose screening compounds and they are not necessarily structurally similar to reference ligands. Encouragingly, we found our hybrid scoring function still outperformed the original docking scoring function. The hybrid scoring function was further evaluated using external datasets for two systems and we found the PI values increased from 0.24 to 0.46 and 0.14 to 0.42 for A2AR and CFX systems, respectively. In a conclusion, our calibration algorithm can significantly improve the virtual screening performance in both drug lead optimization and identification phases with neglectable computational cost.![]()
Collapse
Affiliation(s)
- Beihong Ji
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, PA, 15261, USA
| | - Xibing He
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, PA, 15261, USA
| | - Yuzhao Zhang
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, PA, 15261, USA
| | - Jingchen Zhai
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, PA, 15261, USA
| | - Viet Hoang Man
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, PA, 15261, USA
| | - Shuhan Liu
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, PA, 15261, USA
| | - Junmei Wang
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, PA, 15261, USA.
| |
Collapse
|
13
|
Hao D, He X, Ji B, Zhang S, Wang J. How Well Does the Extended Linear Interaction Energy Method Perform in Accurate Binding Free Energy Calculations? J Chem Inf Model 2020; 60:6624-6633. [PMID: 33213150 DOI: 10.1021/acs.jcim.0c00934] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
With continually increased computer power, molecular mechanics force field-based approaches, such as the endpoint methods of molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) and molecular mechanics generalized Born surface area (MM-GBSA), have been routinely applied in both drug lead identification and optimization. However, the MM-PB/GBSA method is not as accurate as the pathway-based alchemical free energy methods, such as thermodynamic integration (TI) or free energy perturbation (FEP). Although the pathway-based methods are more rigorous in theory, they suffer from slow convergence and computational cost. Moreover, choosing adequate perturbation routes is also crucial for the pathway-based methods. Recently, we proposed a new method, coined extended linear interaction energy (ELIE) method, to overcome some disadvantages of the MM-PB/GBSA method to improve the accuracy of binding free energy calculation. In this work, we have systematically assessed this approach using in total 229 protein-ligand complexes for eight protein targets. Our results showed that ELIE performed much better than the molecular docking and MM-PBSA method in terms of root-mean-square error (RMSE), correlation coefficient (R), predictive index (PI), and Kendall's τ. The mean values of PI, R, and τ are 0.62, 0.58, and 0.44 for ELIE calculations. We also explored the impact of the length of simulation, ranging from 1 to 100 ns, on the performance of binding free energy calculation. In general, extending simulation length up to 25 ns could significantly improve the performance of ELIE, while longer molecular dynamics (MD) simulation does not always perform better than short MD simulation. Considering both the computational efficiency and achieved accuracy, ELIE is adequate in filling the gap between the efficient docking methods and computationally demanding alchemical free energy methods. Therefore, ELIE provides a practical solution for the routine ranking of compounds in lead optimization.
Collapse
Affiliation(s)
- Dongxiao Hao
- School of Physics, Xi'an Jiaotong University, Xi'an 710049, China.,Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Xibing He
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Beihong Ji
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Shengli Zhang
- School of Physics, Xi'an Jiaotong University, Xi'an 710049, China
| | - Junmei Wang
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| |
Collapse
|
14
|
He X, Man VH, Yang W, Lee TS, Wang J. A fast and high-quality charge model for the next generation general AMBER force field. J Chem Phys 2020; 153:114502. [PMID: 32962378 DOI: 10.1063/5.0019056] [Citation(s) in RCA: 202] [Impact Index Per Article: 50.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Abstract
The General AMBER Force Field (GAFF) has been broadly used by researchers all over the world to perform in silico simulations and modelings on diverse scientific topics, especially in the field of computer-aided drug design whose primary task is to accurately predict the affinity and selectivity of receptor-ligand binding. The atomic partial charges in GAFF and the second generation of GAFF (GAFF2) were originally developed with the quantum mechanics derived restrained electrostatic potential charge, but in practice, users usually adopt an efficient charge method, Austin Model 1-bond charge corrections (AM1-BCC), based on which, without expensive ab initio calculations, the atomic charges could be efficiently and conveniently obtained with the ANTECHAMBER module implemented in the AMBER software package. In this work, we developed a new set of BCC parameters specifically for GAFF2 using 442 neutral organic solutes covering diverse functional groups in aqueous solution. Compared to the original BCC parameter set, the new parameter set significantly reduced the mean unsigned error (MUE) of hydration free energies from 1.03 kcal/mol to 0.37 kcal/mol. More excitingly, this new AM1-BCC model also showed excellent performance in the solvation free energy (SFE) calculation on diverse solutes in various organic solvents across a range of different dielectric constants. In this large-scale test with totally 895 neutral organic solvent-solute systems, the new parameter set led to accurate SFE predictions with the MUE and the root-mean-square-error of 0.51 kcal/mol and 0.65 kcal/mol, respectively. This newly developed charge model, ABCG2, paved a promising path for the next generation GAFF development.
Collapse
Affiliation(s)
- Xibing He
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, USA
| | - Viet H Man
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, USA
| | - Wei Yang
- Department of Chemistry and Biochemistry and Institute of Molecular Biophysics, Florida State University, Tallahassee, Florida 32306, USA
| | - Tai-Sung Lee
- Laboratory for Biomolecular Simulation Research, Center for Integrative Proteomics Research, and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854, USA
| | - Junmei Wang
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, USA
| |
Collapse
|
15
|
Abstract
![]()
Accurate determination
of the binding affinity of the ligand to
the receptor remains a difficult problem in computer-aided drug design.
Here, we study and compare the efficiency of Jarzynski’s equality
(JE) combined with steered molecular dynamics and the linear interaction
energy (LIE) method by assessing the binding affinity of 23 small
compounds to six receptors, including β-lactamase, thrombin,
factor Xa, HIV-1 protease (HIV), myeloid cell leukemia-1, and cyclin-dependent
kinase 2 proteins. It was shown that Jarzynski’s nonequilibrium
binding free energy ΔGneqJar correlates with the available
experimental data with the correlation levels R =
0.89, 0.86, 0.83, 0.80, 0.83, and 0.81 for six data sets, while for
the binding free energy ΔGLIE obtained
by the LIE method, we have R = 0.73, 0.80, 0.42,
0.23, 0.85, and 0.01. Therefore, JE is recommended to be used for
ranking binding affinities as it provides accurate and robust results.
In contrast, LIE is not as reliable as JE, and it should be used with
caution, especially when it comes to new systems.
Collapse
Affiliation(s)
- Kiet Ho
- Institute for Computational Sciences and Technology, Quang Trung Software City, SBI Building, Tan Chanh Hiep Ward, District 12, Ho Chi Minh City, Vietnam
| | - Duc Toan Truong
- Institute for Computational Sciences and Technology, Quang Trung Software City, SBI Building, Tan Chanh Hiep Ward, District 12, Ho Chi Minh City, Vietnam.,Department of Theoretical Physics, Faculty of Physics and Engineering Physics, Ho Chi Minh University of Science, Ho Chi Minh City, Vietnam
| | - Mai Suan Li
- Institute of Physics, Polish Academy of Sciences, Al. Lotnikow 32/46, 02-668 Warsaw, Poland
| |
Collapse
|
16
|
Wang E, Liu H, Wang J, Weng G, Sun H, Wang Z, Kang Y, Hou T. Development and Evaluation of MM/GBSA Based on a Variable Dielectric GB Model for Predicting Protein–Ligand Binding Affinities. J Chem Inf Model 2020; 60:5353-5365. [DOI: 10.1021/acs.jcim.0c00024] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Affiliation(s)
- Ercheng Wang
- Hangzhou Institute of Innovative Medicine, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou Zhejiang 310058, China
| | - Hui Liu
- Hangzhou Institute of Innovative Medicine, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou Zhejiang 310058, China
| | - Junmei Wang
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Gaoqi Weng
- Hangzhou Institute of Innovative Medicine, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou Zhejiang 310058, China
| | - Huiyong Sun
- Hangzhou Institute of Innovative Medicine, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou Zhejiang 310058, China
| | - Zhe Wang
- Hangzhou Institute of Innovative Medicine, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou Zhejiang 310058, China
| | - Yu Kang
- Hangzhou Institute of Innovative Medicine, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou Zhejiang 310058, China
| | - Tingjun Hou
- Hangzhou Institute of Innovative Medicine, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou Zhejiang 310058, China
| |
Collapse
|
17
|
He X, Liu S, Lee TS, Ji B, Man VH, York DM, Wang J. Fast, Accurate, and Reliable Protocols for Routine Calculations of Protein-Ligand Binding Affinities in Drug Design Projects Using AMBER GPU-TI with ff14SB/GAFF. ACS OMEGA 2020; 5:4611-4619. [PMID: 32175507 PMCID: PMC7066661 DOI: 10.1021/acsomega.9b04233] [Citation(s) in RCA: 67] [Impact Index Per Article: 16.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/10/2019] [Accepted: 02/13/2020] [Indexed: 05/12/2023]
Abstract
Accurate prediction of the absolute or relative protein-ligand binding affinity is one of the major tasks in computer-aided drug design projects, especially in the stage of lead optimization. In principle, the alchemical free energy (AFE) methods such as thermodynamic integration (TI) or free-energy perturbation (FEP) can fulfill this task, but in practice, a lot of hurdles prevent them from being routinely applied in daily drug design projects, such as the demanding computing resources, slow computing processes, unavailable or inaccurate force field parameters, and difficult and unfriendly setting up and post-analysis procedures. In this study, we have exploited practical protocols of applying the CPU (central processing unit)-TI and newly developed GPU (graphic processing unit)-TI modules and other tools in the AMBER software package, combined with ff14SB/GAFF1.8 force fields, to conduct efficient and accurate AFE calculations on protein-ligand binding free energies. We have tested 134 protein-ligand complexes in total for four target proteins (BACE, CDK2, MCL1, and PTP1B) and obtained overall comparable performance with the commercial Schrodinger FEP+ program (WangJ. Am. Chem. Soc.2015, 137, 2695-2703). The achieved accuracy fits within the requirements for computations to generate effective guidance for experimental work in drug lead optimization, and the needed wall time is short enough for practical application. Our verified protocol provides a practical solution for routine AFE calculations in real drug design projects.
Collapse
Affiliation(s)
- Xibing He
- Department
of Pharmaceutical Sciences and Computational Chemical Genomics Screening
Center, School of Pharmacy, University of
Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Shuhan Liu
- Department
of Pharmaceutical Sciences and Computational Chemical Genomics Screening
Center, School of Pharmacy, University of
Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Tai-Sung Lee
- Laboratory
for Biomolecular Simulation Research, Center for Integrative Proteomics
Research, and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854, United States
| | - Beihong Ji
- Department
of Pharmaceutical Sciences and Computational Chemical Genomics Screening
Center, School of Pharmacy, University of
Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Viet H. Man
- Department
of Pharmaceutical Sciences and Computational Chemical Genomics Screening
Center, School of Pharmacy, University of
Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Darrin M. York
- Laboratory
for Biomolecular Simulation Research, Center for Integrative Proteomics
Research, and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854, United States
| | - Junmei Wang
- Department
of Pharmaceutical Sciences and Computational Chemical Genomics Screening
Center, School of Pharmacy, University of
Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
- E-mail: . Phone: (412) 383-3268. Fax: (412) 383-7436
| |
Collapse
|
18
|
Parks CD, Gaieb Z, Chiu M, Yang H, Shao C, Walters WP, Jansen JM, McGaughey G, Lewis RA, Bembenek SD, Ameriks MK, Mirzadegan T, Burley SK, Amaro RE, Gilson MK. D3R grand challenge 4: blind prediction of protein-ligand poses, affinity rankings, and relative binding free energies. J Comput Aided Mol Des 2020; 34:99-119. [PMID: 31974851 PMCID: PMC7261493 DOI: 10.1007/s10822-020-00289-y] [Citation(s) in RCA: 60] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/05/2019] [Accepted: 01/13/2020] [Indexed: 12/11/2022]
Abstract
The Drug Design Data Resource (D3R) aims to identify best practice methods for computer aided drug design through blinded ligand pose prediction and affinity challenges. Herein, we report on the results of Grand Challenge 4 (GC4). GC4 focused on proteins beta secretase 1 and Cathepsin S, and was run in an analogous manner to prior challenges. In Stage 1, participant ability to predict the pose and affinity of BACE1 ligands were assessed. Following the completion of Stage 1, all BACE1 co-crystal structures were released, and Stage 2 tested affinity rankings with co-crystal structures. We provide an analysis of the results and discuss insights into determined best practice methods.
Collapse
Affiliation(s)
- Conor D Parks
- Drug Design Data Resource, University of California, San Diego, La Jolla, CA, 92093, USA
| | - Zied Gaieb
- Drug Design Data Resource, University of California, San Diego, La Jolla, CA, 92093, USA
| | - Michael Chiu
- Drug Design Data Resource, University of California, San Diego, La Jolla, CA, 92093, USA
| | - Huanwang Yang
- RCSB Protein Data Bank, Institute for Quantitative Biomedicine, Rutgers, The State University of New Jersey, New Brunswick, NJ, 08903, USA
- San Diego Supercomputer Center, University of California, San Diego, La Jolla, CA, 92093, USA
| | - Chenghua Shao
- RCSB Protein Data Bank, Institute for Quantitative Biomedicine, Rutgers, The State University of New Jersey, New Brunswick, NJ, 08903, USA
- San Diego Supercomputer Center, University of California, San Diego, La Jolla, CA, 92093, USA
| | | | - Johanna M Jansen
- Novartis Institutes for BioMedical Research, Emeryville, CA, 94608, USA
| | | | - Richard A Lewis
- Novartis Institutes for BioMedical Research, Novartis Pharma AG, 4002, Basel, Switzerland
| | | | | | | | - Stephen K Burley
- RCSB Protein Data Bank, Institute for Quantitative Biomedicine, Rutgers, The State University of New Jersey, New Brunswick, NJ, 08903, USA
- San Diego Supercomputer Center, University of California, San Diego, La Jolla, CA, 92093, USA
| | - Rommie E Amaro
- Drug Design Data Resource, University of California, San Diego, La Jolla, CA, 92093, USA.
- Department of Chemistry and Biochemistry, UC San Diego, La Jolla, CA, 92093-0340, USA.
| | - Michael K Gilson
- Drug Design Data Resource, University of California, San Diego, La Jolla, CA, 92093, USA.
- Skaggs School of Pharmacy and Pharmaceutical Sciences, UC San Diego, 9500 Gilman Drive, MC0751, La Jolla, CA, 92093, USA.
| |
Collapse
|
19
|
Kadukova M, Chupin V, Grudinin S. Docking rigid macrocycles using Convex-PL, AutoDock Vina, and RDKit in the D3R Grand Challenge 4. J Comput Aided Mol Des 2019; 34:191-200. [PMID: 31784861 DOI: 10.1007/s10822-019-00263-3] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2019] [Accepted: 11/22/2019] [Indexed: 12/15/2022]
Abstract
The D3R Grand Challenge 4 provided a brilliant opportunity to test macrocyclic docking protocols on a diverse high-quality experimental data. We participated in both pose and affinity prediction exercises. Overall, we aimed to use an automated structure-based docking pipeline built around a set of tools developed in our team. This exercise again demonstrated a crucial importance of the correct local ligand geometry for the overall success of docking. Starting from the second part of the pose prediction stage, we developed a stable pipeline for sampling macrocycle conformers. This resulted in the subangstrom average precision of our pose predictions. In the affinity prediction exercise we obtained average results. However, we could improve these when using docking poses submitted by the best predictors. Our docking tools including the Convex-PL scoring function are available at https://team.inria.fr/nano-d/software/.
Collapse
Affiliation(s)
- Maria Kadukova
- Univ. Grenoble Alpes, CNRS, Inria, Grenoble INP, LJK, 38000, Grenoble, France
- Moscow Institute of Physics and Technology, Dolgoprudny, Russia, 141700
| | - Vladimir Chupin
- Moscow Institute of Physics and Technology, Dolgoprudny, Russia, 141700
| | - Sergei Grudinin
- Univ. Grenoble Alpes, CNRS, Inria, Grenoble INP, LJK, 38000, Grenoble, France.
| |
Collapse
|
20
|
Huang S, Cao Y. Correlation of cathepsin S with coronary stenosis degree, carotid thickness, blood pressure, glucose and lipid metabolism and vascular endothelial function in atherosclerosis. Exp Ther Med 2019; 19:61-66. [PMID: 31853273 PMCID: PMC6909596 DOI: 10.3892/etm.2019.8222] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2019] [Accepted: 09/12/2019] [Indexed: 02/06/2023] Open
Abstract
Correlation of cathepsin S with coronary stenosis degree, carotid thickness, blood pressure, glucose and lipid metabolism, and vascular endothelial function in patients with atherosclerosis was investigated. Data from 120 patients with increased cathepsin S levels (increased group) and 120 subjects with normal cathepsin S levels (normal group) were retrospectively analyzed. The serum cathepsin S level and Gensini score were compared between the healthy subjects and patients with coronary atherosclerotic heart disease, and the correlation between serum cathepsin S level and Gensini score was analyzed. The carotid thickness, mean arterial pressure, and indexes related to glucose and lipid metabolism, as well as vascular endothelial function were compared between the two groups. The correlation of the serum cathepsin S level with carotid intima-media thickness (IMT), mean arterial pressure, fasting blood glucose, total cholesterol (TC) and nitric oxide (NO) was also investigated. Patients with multi-vessel coronary artery disease (CAD) had higher serum cathepsin S level than those with double-vessel and single-vessel disease, and higher level than that of healthy subjects. The Gensini score of patients with multi-vessel CAD was higher than that of patients with double-vessel and single-vessel disease, as well as that of healthy subjects. The serum cathepsin S level was positively correlated with the Gensini score. Patients with increased cathepsin S level had greater IMT, higher mean arterial pressure, fasting blood glucose, fasting insulin (FINS), triglyceride (TG), TC, and endothelin-1 (ET-1), however, lower NO level than those of healthy subjects. The serum cathepsin S level was positively correlated with IMT, mean arterial pressure, fasting blood glucose, and TC, however, it was negatively correlated with the NO level. In conclusion, as the serum cathepsin S level is elevated, the coronary stenosis is aggravated, the carotid thickness and blood pressure are increased, and the glucose and lipid metabolism, as well as vascular endothelial function are significantly abnormal.
Collapse
Affiliation(s)
- Shengyang Huang
- Cardiovascular Department, The Third Xiangya Hospital of Central South University, Changsha, Hunan 410000, P.R. China
| | - Yu Cao
- Cardiovascular Department, The Third Xiangya Hospital of Central South University, Changsha, Hunan 410000, P.R. China
| |
Collapse
|
21
|
Liu S, He X, Man VH, Ji B, Liu J, Wang J. New application of in silico methods in identifying mechanisms of action and key components of anti-cancer herbal formulation YIV-906 (PHY906). Phys Chem Chem Phys 2019; 21:23501-23513. [PMID: 31617551 DOI: 10.1039/c9cp03803e] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
YIV-906 (formally PHY906, KD018) is a four-herb formulation that is currently being developed to improve the therapeutic index and ameliorate the side effects of many chemotherapeutic drugs including sorafenib, irinotecan, and capecitabine. However, as a promising anti-cancer adjuvant, the molecular mechanism of action of YIV-906 remains unrevealed due to its multi-component and multi-target features. Since YIV-906 has been shown to induce apoptosis and autophagy in cancer cells through modulating the negative regulators of ERK1/2, namely DUSPs, it is of great interest to elucidate the key components that cause the therapeutic effect of YIV-906. In this work, we investigated the mechanism of YIV-906 inhibiting DUSPs, using a broad spectrum of molecular modelling techniques, including molecular docking, molecular dynamics (MD) simulations, and binding free energy calculations. In total, MD simulations and binding free energy calculations were performed for 99 DUSP-ligand complexes. We found that some herbal components or their metabolites could inhibit DUSPs. Based on the docking scores and binding free energies, the sulfation and glucuronidation metabolites of the S ingredient in YIV-906 play a leading role in inhibiting DUSPs, although several original herbal chemicals with carboxyl groups from the P and Z ingredients also make contributions to this inhibitory effect. It is not a surprise that the electrostatic interaction plays the dominant role in the ligand binding process, given the fact that several charged residues reside in the binding pockets of DUSPs. Our MD simulation results demonstrate that the sulfate moieties and carboxyl moieties of the advantageous ligands from YIV-906 can occupy the enzymes' catalytic sites, mimicking the endogenous phosphate substrates of DUSPs. As such, the ligand binding can inhibit the association of DUSPs and ERK1/2, which in turn reduces the dephosphorylation of ERK1/2 and causes cell cycle arrest in the tumor. Our modelling study provides useful insights into the rational design of highly potent anti-cancer drugs targeting DUSPs. Finally, we have demonstrated that multi-scale molecular modelling techniques are able to elucidate molecular mechanisms involving complex molecular systems.
Collapse
Affiliation(s)
- Shuhan Liu
- School of Pharmacy, Computational Chemical Genomics Screening Center, University of Pittsburgh, 3501 Terrace St, Pittsburgh, Pennsylvania 15261, USA.
| | | | | | | | | | | |
Collapse
|
22
|
Rifai EA, van Dijk M, Vermeulen NPE, Yanuar A, Geerke DP. A Comparative Linear Interaction Energy and MM/PBSA Study on SIRT1-Ligand Binding Free Energy Calculation. J Chem Inf Model 2019; 59:4018-4033. [PMID: 31461271 PMCID: PMC6759767 DOI: 10.1021/acs.jcim.9b00609] [Citation(s) in RCA: 50] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2019] [Indexed: 12/25/2022]
Abstract
Binding free energy (ΔGbind) computation can play an important role in prioritizing compounds to be evaluated experimentally on their affinity for target proteins, yet fast and accurate ΔGbind calculation remains an elusive task. In this study, we compare the performance of two popular end-point methods, i.e., linear interaction energy (LIE) and molecular mechanics/Poisson-Boltzmann surface area (MM/PBSA), with respect to their ability to correlate calculated binding affinities of 27 thieno[3,2-d]pyrimidine-6-carboxamide-derived sirtuin 1 (SIRT1) inhibitors with experimental data. Compared with the standard single-trajectory setup of MM/PBSA, our study elucidates that LIE allows to obtain direct ("absolute") values for SIRT1 binding free energies with lower compute requirements, while the accuracy in calculating relative values for ΔGbind is comparable (Pearson's r = 0.72 and 0.64 for LIE and MM/PBSA, respectively). We also investigate the potential of combining multiple docking poses in iterative LIE models and find that Boltzmann-like weighting of outcomes of simulations starting from different poses can retrieve appropriate binding orientations. In addition, we find that in this particular case study the LIE and MM/PBSA models can be optimized by neglecting the contributions from electrostatic and polar interactions to the ΔGbind calculations.
Collapse
Affiliation(s)
- Eko Aditya Rifai
- AIMMS
Division of Molecular and Computational Toxicology, Department of
Chemistry and Pharmaceutical Sciences, Vrije
Universiteit Amsterdam, De Boelelaan 1108, 1081 HZ Amsterdam, The Netherlands
| | - Marc van Dijk
- AIMMS
Division of Molecular and Computational Toxicology, Department of
Chemistry and Pharmaceutical Sciences, Vrije
Universiteit Amsterdam, De Boelelaan 1108, 1081 HZ Amsterdam, The Netherlands
| | - Nico P. E. Vermeulen
- AIMMS
Division of Molecular and Computational Toxicology, Department of
Chemistry and Pharmaceutical Sciences, Vrije
Universiteit Amsterdam, De Boelelaan 1108, 1081 HZ Amsterdam, The Netherlands
| | - Arry Yanuar
- Faculty
of Pharmacy, Universitas Indonesia, Depok 16424, Indonesia
| | - Daan P. Geerke
- AIMMS
Division of Molecular and Computational Toxicology, Department of
Chemistry and Pharmaceutical Sciences, Vrije
Universiteit Amsterdam, De Boelelaan 1108, 1081 HZ Amsterdam, The Netherlands
| |
Collapse
|
23
|
Wang E, Sun H, Wang J, Wang Z, Liu H, Zhang JZH, Hou T. End-Point Binding Free Energy Calculation with MM/PBSA and MM/GBSA: Strategies and Applications in Drug Design. Chem Rev 2019; 119:9478-9508. [DOI: 10.1021/acs.chemrev.9b00055] [Citation(s) in RCA: 578] [Impact Index Per Article: 115.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Ercheng Wang
- Hangzhou Institute of Innovative Medicine, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China
| | - Huiyong Sun
- Hangzhou Institute of Innovative Medicine, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China
| | - Junmei Wang
- Department of Pharmaceutical Sciences, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Zhe Wang
- Hangzhou Institute of Innovative Medicine, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China
| | - Hui Liu
- Hangzhou Institute of Innovative Medicine, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China
| | - John Z. H. Zhang
- Shanghai Engineering Research Center of Molecular Therapeutics & New Drug Development, Shanghai Key Laboratory of Green Chemistry & Chemical Process, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
- NYU−ECNU Center for Computational Chemistry, NYU Shanghai, Shanghai 200122, China
- Department of Chemistry, New York University, New York, New York 10003, United States
- Collaborative Innovation Center of Extreme Optics, Shanxi University, Taiyuan, Shanxi 030006, China
| | - Tingjun Hou
- Hangzhou Institute of Innovative Medicine, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China
| |
Collapse
|
24
|
Bitencourt-Ferreira G, Veit-Acosta M, de Azevedo WF. Van der Waals Potential in Protein Complexes. Methods Mol Biol 2019; 2053:79-91. [PMID: 31452100 DOI: 10.1007/978-1-4939-9752-7_6] [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] [Indexed: 12/13/2022]
Abstract
Van der Waals forces are determinants of the formation of protein-ligand complexes. Physical models based on the Lennard-Jones potential can estimate van der Waals interactions with considerable accuracy and with a computational complexity that allows its application to molecular docking simulations and virtual screening of large databases of small organic molecules. Several empirical scoring functions used to evaluate protein-ligand interactions approximate van der Waals interactions with the Lennard-Jones potential. In this chapter, we present the main concepts necessary to understand van der Waals interactions relevant to molecular recognition of a ligand by the binding pocket of a protein target. We describe the Lennard-Jones potential and its application to calculate potential energy for an ensemble of structures to highlight the main features related to the importance of this interaction for binding affinity.
Collapse
Affiliation(s)
- Gabriela Bitencourt-Ferreira
- Escola de Ciências da Saúde, Pontifícia Universidade Católica do Rio Grande do Sul-PUCRS, Porto Alegre, RS, Brazil
| | - Martina Veit-Acosta
- Escola de Ciências da Saúde, Pontifícia Universidade Católica do Rio Grande do Sul-PUCRS, Porto Alegre, RS, Brazil
| | - Walter Filgueira de Azevedo
- Escola de Ciências da Saúde, Pontifícia Universidade Católica do Rio Grande do Sul-PUCRS, Porto Alegre, RS, Brazil.
| |
Collapse
|