1
|
Gu X, Aranganathan A, Tiwary P. Empowering AlphaFold2 for protein conformation selective drug discovery with AlphaFold2-RAVE. eLife 2024; 13:RP99702. [PMID: 39240197 PMCID: PMC11379456 DOI: 10.7554/elife.99702] [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] [Indexed: 09/07/2024] Open
Abstract
Small-molecule drug design hinges on obtaining co-crystallized ligand-protein structures. Despite AlphaFold2's strides in protein native structure prediction, its focus on apo structures overlooks ligands and associated holo structures. Moreover, designing selective drugs often benefits from the targeting of diverse metastable conformations. Therefore, direct application of AlphaFold2 models in virtual screening and drug discovery remains tentative. Here, we demonstrate an AlphaFold2-based framework combined with all-atom enhanced sampling molecular dynamics and Induced Fit docking, named AF2RAVE-Glide, to conduct computational model-based small-molecule binding of metastable protein kinase conformations, initiated from protein sequences. We demonstrate the AF2RAVE-Glide workflow on three different mammalian protein kinases and their type I and II inhibitors, with special emphasis on binding of known type II kinase inhibitors which target the metastable classical DFG-out state. These states are not easy to sample from AlphaFold2. Here, we demonstrate how with AF2RAVE these metastable conformations can be sampled for different kinases with high enough accuracy to enable subsequent docking of known type II kinase inhibitors with more than 50% success rates across docking calculations. We believe the protocol should be deployable for other kinases and more proteins generally.
Collapse
Affiliation(s)
- Xinyu Gu
- Institute for Physical Science and Technology, University of Maryland, College Park, United States
- University of Maryland Institute for Health Computing, Bethesda, United States
| | - Akashnathan Aranganathan
- Institute for Physical Science and Technology, University of Maryland, College Park, United States
- Biophysics Program, University of Maryland, College Park, United States
| | - Pratyush Tiwary
- Institute for Physical Science and Technology, University of Maryland, College Park, United States
- University of Maryland Institute for Health Computing, Bethesda, United States
- Department of Chemistry and Biochemistry, University of Maryland, College Park, United States
| |
Collapse
|
2
|
Gu X, Aranganathan A, Tiwary P. Empowering AlphaFold2 for protein conformation selective drug discovery with AlphaFold2-RAVE. ARXIV 2024:arXiv:2404.07102v3. [PMID: 38659642 PMCID: PMC11042445] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Subscribe] [Scholar Register] [Indexed: 04/26/2024]
Abstract
Small molecule drug design hinges on obtaining co-crystallized ligand-protein structures. Despite AlphaFold2's strides in protein native structure prediction, its focus on apo structures overlooks ligands and associated holo structures. Moreover, designing selective drugs often benefits from the targeting of diverse metastable conformations. Therefore, direct application of AlphaFold2 models in virtual screening and drug discovery remains tentative. Here, we demonstrate an AlphaFold2 based framework combined with all-atom enhanced sampling molecular dynamics and induced fit docking, named AF2RAVE-Glide, to conduct computational model based small molecule binding of metastable protein kinase conformations, initiated from protein sequences. We demonstrate the AF2RAVE-Glide workflow on three different protein kinases and their type I and II inhibitors, with special emphasis on binding of known type II kinase inhibitors which target the metastable classical DFG-out state. These states are not easy to sample from AlphaFold2. Here we demonstrate how with AF2RAVE these metastable conformations can be sampled for different kinases with high enough accuracy to enable subsequent docking of known type II kinase inhibitors with more than 50% success rates across docking calculations. We believe the protocol should be deployable for other kinases and more proteins generally.
Collapse
Affiliation(s)
- Xinyu Gu
- Institute for Physical Science and Technology, University of Maryland, College Park, Maryland 20742, USA
- University of Maryland Institute for Health Computing, Bethesda, United States
| | - Akashnathan Aranganathan
- Institute for Physical Science and Technology, University of Maryland, College Park, Maryland 20742, USA
- Biophysics Program, University of Maryland, College Park 20742, USA
| | - Pratyush Tiwary
- Institute for Physical Science and Technology, University of Maryland, College Park, Maryland 20742, USA
- Department of Chemistry and Biochemistry, University of Maryland, College Park 20742, USA
- University of Maryland Institute for Health Computing, Bethesda, United States
| |
Collapse
|
3
|
Gu S, Yang Y, Zhao Y, Qiu J, Wang X, Tong HHY, Liu L, Wan X, Liu H, Hou T, Kang Y. Evaluation of AlphaFold2 Structures for Hit Identification across Multiple Scenarios. J Chem Inf Model 2024; 64:3630-3639. [PMID: 38630855 DOI: 10.1021/acs.jcim.3c01976] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/19/2024]
Abstract
The introduction of AlphaFold2 (AF2) has sparked significant enthusiasm and generated extensive discussion within the scientific community, particularly among drug discovery researchers. Although previous studies have addressed the performance of AF2 structures in virtual screening (VS), a more comprehensive investigation is still necessary considering the paramount importance of structural accuracy in drug design. In this study, we evaluate the performance of AF2 structures in VS across three common drug discovery scenarios: targets with holo, apo, and AF2 structures; targets with only apo and AF2 structures; and targets exclusively with AF2 structures. We utilized both the traditional physics-based Glide and the deep-learning-based scoring function RTMscore to rank the compounds in the DUD-E, DEKOIS 2.0, and DECOY data sets. The results demonstrate that, overall, the performance of VS on AF2 structures is comparable to that on apo structures but notably inferior to that on holo structures across diverse scenarios. Moreover, when a target has solely AF2 structure, selecting the holo structure of the target from different subtypes within the same protein family produces comparable results with the AF2 structure for VS on the data set of the AF2 structures, and significantly better results than the AF2 structures on its own data set. This indicates that utilizing AF2 structures for docking-based VS may not yield most satisfactory outcomes, even when solely AF2 structures are available. Moreover, we rule out the possibility that the variations in VS performance between the binding pockets of AF2 and holo structures arise from the differences in their biological assembly composition.
Collapse
Affiliation(s)
- Shukai Gu
- Faculty of Applied Science, Macao Polytechnic University, Macao 999078, SAR, China
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, Zhejiang, China
| | - Yuwei Yang
- Faculty of Applied Science, Macao Polytechnic University, Macao 999078, SAR, China
| | - Yihao Zhao
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, Zhejiang, China
| | - Jiayue Qiu
- Faculty of Applied Science, Macao Polytechnic University, Macao 999078, SAR, China
| | - Xiaorui Wang
- State Key Laboratory of Quality Re-search in Chinese Medicine, Macau Institute for Applied Research in Medicine and Health, Macau University of Science and Technology, Macao 999078, China
| | - Henry Hoi Yee Tong
- Faculty of Applied Science, Macao Polytechnic University, Macao 999078, SAR, China
| | - Liwei Liu
- Advanced Computing and Storage Laboratory, Central Research Institute, 2012 Laboratories, Huawei Technologies Co., Ltd., Nanjing 210000, Jiangsu, China
| | - Xiaozhe Wan
- Advanced Computing and Storage Laboratory, Central Research Institute, 2012 Laboratories, Huawei Technologies Co., Ltd., Nanjing 210000, Jiangsu, China
| | - Huanxiang Liu
- Faculty of Applied Science, Macao Polytechnic University, Macao 999078, SAR, China
| | - Tingjun Hou
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, Zhejiang, China
| | - Yu Kang
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, Zhejiang, China
| |
Collapse
|
4
|
Zhang H, Lu C, Yao Q, Jiao Q. In silico study to identify novel NEK7 inhibitors from natural sources by a combination strategy. Mol Divers 2024:10.1007/s11030-024-10838-4. [PMID: 38598164 DOI: 10.1007/s11030-024-10838-4] [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: 09/02/2023] [Accepted: 03/06/2024] [Indexed: 04/11/2024]
Abstract
Cancer poses a significant global health challenge and significantly contributes to mortality. NEK7, related to the NIMA protein kinase family, plays a crucial role in spindle assembly and cell division. The dysregulation of NEK7 is closely linked to the onset and progression of various cancers, especially colon and breast cancer, making it a promising target for cancer therapy. Nevertheless, the shortage of high-quality NEK7 inhibitors highlights the need for new therapeutic strategies. In this study, we utilized a multidisciplinary approach, including virtual screening, molecular docking, pharmacokinetics, molecular dynamics simulations (MDs), and MM/PBSA calculations, to evaluate natural compounds as NEK7 inhibitors comprehensively. Through various docking strategies, we identified three natural compounds: (-)-balanol, digallic acid, and scutellarin. Molecular docking revealed significant interactions at residues such as GLU112 and ALA114, with docking scores of -15.054, -13.059, and -11.547 kcal/mol, respectively, highlighting their potential as NEK7 inhibitors. MDs confirmed the stability of these compounds at the NEK7-binding site. Hydrogen bond analysis during simulations revealed consistent interactions, supporting their strong binding capacity. MM/PBSA analysis identified other crucial amino acids contributing to binding affinity, including ILE20, VAL28, ILE75, LEU93, ALA94, LYS143, PHE148, LEU160, and THR161, crucial for stabilizing the complex. This research demonstrated that these compounds exceeded dabrafenib in binding energy, according to MM/PBSA calculations, underscoring their effectiveness as NEK7 inhibitors. ADME/T predictions showed lower oral toxicity for these compounds, suggesting their potential for further development. This study highlights the promise of these natural compounds as bases for creating more potent derivatives with significant biological activities, paving the way for future experimental validation.
Collapse
Affiliation(s)
- Heng Zhang
- State Key Laboratory of Pharmaceutical Biotechnology, School of Life Sciences, Nanjing University, Nanjing, 210023, China
| | - Chenhong Lu
- State Key Laboratory of Pharmaceutical Biotechnology, School of Life Sciences, Nanjing University, Nanjing, 210023, China
| | - Qilong Yao
- State Key Laboratory of Pharmaceutical Biotechnology, School of Life Sciences, Nanjing University, Nanjing, 210023, China
| | - Qingcai Jiao
- State Key Laboratory of Pharmaceutical Biotechnology, School of Life Sciences, Nanjing University, Nanjing, 210023, China.
| |
Collapse
|
5
|
Yang H, Kang M, Jang S, Baek SY, Kim J, Kim GU, Kim D, Ha J, Kim JS, Jung C, Kim NJ, Cho SY, Shin WH, Lee J, Ko J, Lee A, Keum G, Lee S, Kang T. Discovery of thiophen-2-ylmethylene bis-dimedone derivatives as novel WRN inhibitors for treating cancers with microsatellite instability. Bioorg Med Chem 2024; 100:117588. [PMID: 38295487 DOI: 10.1016/j.bmc.2024.117588] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2023] [Revised: 11/23/2023] [Accepted: 01/08/2024] [Indexed: 02/02/2024]
Abstract
Microsatellite instability (MSI) is a hypermutable condition caused by DNA mismatch repair system defects, contributing to the development of various cancer types. Recent research has identified Werner syndrome ATP-dependent helicase (WRN) as a promising synthetic lethal target for MSI cancers. Herein, we report the first discovery of thiophen-2-ylmethylene bis-dimedone derivatives as novel WRN inhibitors for MSI cancer therapy. Initial computational analysis and biological evaluation identified a new scaffold for a WRN inhibitor. Subsequent SAR study led to the discovery of a highly potent WRN inhibitor. Furthermore, we demonstrated that the optimal compound induced DNA damage and apoptotic cell death in MSI cancer cells by inhibiting WRN. This study provides a new pharmacophore for WRN inhibitors, emphasizing their therapeutic potential for MSI cancers.
Collapse
Affiliation(s)
- Hwasun Yang
- Brain Science Institute, Korea Institute of Science and Technology, Seoul 02792, Republic of Korea; Department of Chemistry, Korea University, Seoul 02841, Republic of Korea
| | - Miso Kang
- Brain Science Institute, Korea Institute of Science and Technology, Seoul 02792, Republic of Korea; Department of Fundamental Pharmaceutical Sciences, Graduate School, Kyung Hee University, Seoul 02447, Republic of Korea
| | - Seonyeong Jang
- Brain Science Institute, Korea Institute of Science and Technology, Seoul 02792, Republic of Korea; Department of Biotechnology, College of Life Sciences and Biotechnology, Korea University, Seoul 02841, Republic of Korea
| | - Soo Yeon Baek
- Brain Science Institute, Korea Institute of Science and Technology, Seoul 02792, Republic of Korea
| | - Jiwon Kim
- Department of Biomedical Sciences, Seoul National University College of Medicine, Seoul 03080, Republic of Korea
| | - Gyeong Un Kim
- Brain Science Institute, Korea Institute of Science and Technology, Seoul 02792, Republic of Korea; Division of Bio-Medical Science & Technology, KIST School, University of Science and Technology, Seoul 02792, Republic of Korea
| | - Dongwoo Kim
- College of Pharmacy, Seoul National University, Seoul 08826, Republic of Korea
| | - Junsu Ha
- Arontier Co., Ltd., Seoul 06735, Republic of Korea
| | - Jong Seung Kim
- Department of Chemistry, Korea University, Seoul 02841, Republic of Korea
| | - Cheulhee Jung
- Department of Biotechnology, College of Life Sciences and Biotechnology, Korea University, Seoul 02841, Republic of Korea
| | - Nam-Jung Kim
- Department of Fundamental Pharmaceutical Sciences, Graduate School, Kyung Hee University, Seoul 02447, Republic of Korea
| | - Sung-Yup Cho
- Department of Biomedical Sciences, Seoul National University College of Medicine, Seoul 03080, Republic of Korea; Medical Research Center, Genomic Medicine Institute, Seoul National University College of Medicine, Seoul 03080, Republic of Korea; Cancer Research Institute, Seoul National University, Seoul 03080, Republic of Korea
| | - Woong-Hee Shin
- Arontier Co., Ltd., Seoul 06735, Republic of Korea; Department of Medicine, Korea University College of Medicine, Seoul 02708, Republic of Korea
| | - Juyong Lee
- Arontier Co., Ltd., Seoul 06735, Republic of Korea; Research Institute of Pharmaceutical Science, Seoul National University, Seoul 08826, Republic of Korea; Molecular Medicine and Biopharmaceutical Sciences, Graduate School of Convergence Science and Technology, Seoul National University, Seoul 08826, Republic of Korea
| | - Junsu Ko
- Arontier Co., Ltd., Seoul 06735, Republic of Korea
| | - Ansoo Lee
- Brain Science Institute, Korea Institute of Science and Technology, Seoul 02792, Republic of Korea; Division of Bio-Medical Science & Technology, KIST School, University of Science and Technology, Seoul 02792, Republic of Korea
| | - Gyochang Keum
- Brain Science Institute, Korea Institute of Science and Technology, Seoul 02792, Republic of Korea; Division of Bio-Medical Science & Technology, KIST School, University of Science and Technology, Seoul 02792, Republic of Korea
| | - Sanghee Lee
- Brain Science Institute, Korea Institute of Science and Technology, Seoul 02792, Republic of Korea; Department for HY-KIST Bio-convergence, Hanyang University, Seoul 04763, Republic of Korea.
| | - Taek Kang
- Brain Science Institute, Korea Institute of Science and Technology, Seoul 02792, Republic of Korea.
| |
Collapse
|
6
|
Bouguerra OM, Wahab RA, Huyop F, Al-Fakih AM, Mahmood WMAW, Mahat NA, Sabullah MK. An Overview of Crosslinked Enzyme Aggregates: Concept of Development and Trends of Applications. Appl Biochem Biotechnol 2024:10.1007/s12010-023-04809-y. [PMID: 38180645 DOI: 10.1007/s12010-023-04809-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 12/09/2023] [Indexed: 01/06/2024]
Abstract
Enzymes are commonly used as biocatalysts for various biological and chemical processes in industrial applications. However, their limited operational stability, catalytic efficiency, poor reusability, and high-cost hamper further industrial usage. Thus, crosslinked enzyme aggregates (CLEAs) are developed as a better enzyme immobilization tool to extend the enzymes' operational stability. This immobilization method is appealing because it is simpler due to the absence of ballast and permits the collective use of crude enzyme cocktails. CLEAs, so far, have been successfully developed using a variety of enzymes, viz., hydrolases, proteases, amidases, lipases, esterases, and oxidoreductase. Recent years have seen the emergence of novel strategies for preparing better CLEAs, which include the combi- and multi-CLEAs, magnetics CLEAs, and porous CLEAs for various industrial applications, viz., laundry detergents, organic synthesis, food industries, pharmaceutical applications, oils, and biodiesel production. To better understand the different strategies for CLEAs' development, this review explores these strategies and highlights the relevant concerns in designing innovative CLEAs. This article also details the challenges faced during CLEAs preparation and solutions for overcoming them. Finally, the trending strategies to improve the preparation of CLEAs alongside their industrial application trends are also discussed.
Collapse
Affiliation(s)
- Oumaima Maroua Bouguerra
- Department of Bioscience, Faculty of Science, Universiti Teknologi Malaysia, UTM, 81310, Johor Bahru, Johor, Malaysia
| | - Roswanira Abdul Wahab
- Department of Chemistry, Faculty of Science, Universiti Teknologi Malaysia, UTM, 81310, Johor Bahru, Johor, Malaysia.
- Advanced Membrane Technology Research Centre (AMTEC), Universiti Teknologi Malaysia, UTM, 81310, Johor Bahru, Malaysia.
| | - Fahrul Huyop
- Department of Bioscience, Faculty of Science, Universiti Teknologi Malaysia, UTM, 81310, Johor Bahru, Johor, Malaysia
| | - Abdo Mohammed Al-Fakih
- Department of Chemistry, Faculty of Science, Universiti Teknologi Malaysia, UTM, 81310, Johor Bahru, Johor, Malaysia
| | - Wan Muhd Asyraf Wan Mahmood
- Centre of Foundation Studies, Dengkil Campus, Universiti Teknologi MARA (UiTM) Selangor Branch, 43800, Dengkil, Selangor, Malaysia
| | - Naji Arafat Mahat
- Department of Chemistry, Faculty of Science, Universiti Teknologi Malaysia, UTM, 81310, Johor Bahru, Johor, Malaysia
| | - Mohd Khalizan Sabullah
- Faculty of Science and Natural Resources, Universiti Malaysia Sabah, Jalan UMS, 88400, Kota Kinabalu, Sabah, Malaysia.
| |
Collapse
|
7
|
Nicoli A, Weber V, Bon C, Steuer A, Gustincich S, Gainetdinov RR, Lang R, Espinoza S, Di Pizio A. Structure-Based Discovery of Mouse Trace Amine-Associated Receptor 5 Antagonists. J Chem Inf Model 2023; 63:6667-6680. [PMID: 37847527 PMCID: PMC10647090 DOI: 10.1021/acs.jcim.3c00755] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/18/2023] [Indexed: 10/18/2023]
Abstract
Trace amine-associated receptors (TAARs) were discovered in 2001 as new members of class A G protein-coupled receptors (GPCRs). With the only exception of TAAR1, TAAR members (TAAR2-9, also known as noncanonical olfactory receptors) were originally described exclusively in the olfactory epithelium and believed to mediate the innate perception of volatile amines. However, most noncanonical olfactory receptors are still orphan receptors. Given its recently discovered nonolfactory expression and therapeutic potential, TAAR5 has been the focus of deorphanization campaigns that led to the discovery of a few druglike antagonists. Here, we report four novel TAAR5 antagonists identified through high-throughput screening, which, along with the four ligands published in the literature, constituted our starting point to design a computational strategy for the identification of TAAR5 ligands. We developed a structure-based virtual screening protocol that allowed us to identify three new TAAR5 antagonists with a hit rate of 10%. Despite lacking an experimental structure, we accurately modeled the TAAR5 binding site by integrating comparative sequence- and structure-based analyses of serotonin receptors with homology modeling and side-chain optimization. In summary, we have identified seven new TAAR5 antagonists that could serve as lead candidates for the development of new treatments for depression, anxiety, and neurodegenerative diseases.
Collapse
Affiliation(s)
- Alessandro Nicoli
- Leibniz
Institute for Food Systems Biology at the Technical University of
Munich, 85354 Freising, Germany
- Chemoinformatics
and Protein Modelling, Department of Molecular Life Sciences, School
of Life Sciences, Technical University of
Munich, 85354 Freising, Germany
| | - Verena Weber
- Leibniz
Institute for Food Systems Biology at the Technical University of
Munich, 85354 Freising, Germany
- Institute
for Advanced Simulations (IAS)-5/Institute for Neuroscience and Medicine
(INM)-9, Forschungszentrum Jülich, 52428 Jülich, Germany
- Faculty
of Mathematics, Computer Science and Natural Sciences, RWTH Aachen, Aachen, 52062 Germany
| | - Carlotta Bon
- Istituto
Italiano di Tecnologia, 16163 Genova, Italy
| | - Alexandra Steuer
- Leibniz
Institute for Food Systems Biology at the Technical University of
Munich, 85354 Freising, Germany
- Chemoinformatics
and Protein Modelling, Department of Molecular Life Sciences, School
of Life Sciences, Technical University of
Munich, 85354 Freising, Germany
| | | | - Raul R. Gainetdinov
- Institute
of Translational Biomedicine and Saint Petersburg University Hospital,
Saint Petersburg State University, Saint Petersburg 199034, Russia
| | - Roman Lang
- Leibniz
Institute for Food Systems Biology at the Technical University of
Munich, 85354 Freising, Germany
| | - Stefano Espinoza
- Istituto
Italiano di Tecnologia, 16163 Genova, Italy
- Dipartimento
di Scienze della Salute, Università
del Piemonte Orientale, 28100 Novara, Italy
| | - Antonella Di Pizio
- Leibniz
Institute for Food Systems Biology at the Technical University of
Munich, 85354 Freising, Germany
- Chemoinformatics
and Protein Modelling, Department of Molecular Life Sciences, School
of Life Sciences, Technical University of
Munich, 85354 Freising, Germany
| |
Collapse
|
8
|
Nguyen LP, Khan RA, Kang S, Lee H, Hwang JI, Kim HR. Discovery of Chemical Scaffolds as Lysophosphatidic Acid Receptor 1 Antagonists: Virtual Screening, In Vitro Validation, and Molecular Dynamics Analysis. ACS OMEGA 2023; 8:40375-40386. [PMID: 37929144 PMCID: PMC10620911 DOI: 10.1021/acsomega.3c04798] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/04/2023] [Accepted: 10/06/2023] [Indexed: 11/07/2023]
Abstract
Lysophosphatidic acid receptor 1 (LPAR1) is an emerging therapeutic target for numerous human diseases including fibrosis. However, the limited number of available core structures of LPAR1 antagonists has prompted the need for novel chemical templates. In this study, we conducted a high-throughput virtual screening to discover potential new scaffolds. We tested three existing crystal structures alongside an AlphaFold model to evaluate their suitability in structure-based virtual screening, finding that the crystal structures show superior performance compared with the predictive model. Furthermore, we also found that enhancing the precision in the screening process did not necessarily improve the enrichment of hits. From the screening campaign, we identified five structures that were validated using an LPAR1-dependent calcium flux assay. To gain a deeper insight into the protein-ligand interaction, we extensively analyzed the binding modes of these compounds using in silico techniques, laying the groundwork for the discovery of novel LPAR1 antagonists.
Collapse
Affiliation(s)
- Lan Phuong Nguyen
- Department of Biomedical Sciences,
College of Medicine, Korea University, Seoul 02841, Republic of Korea
| | - Rasel Ahmed Khan
- Department of Biomedical Sciences,
College of Medicine, Korea University, Seoul 02841, Republic of Korea
| | - Soomin Kang
- Department of Biomedical Sciences,
College of Medicine, Korea University, Seoul 02841, Republic of Korea
| | - Hobin Lee
- Department of Biomedical Sciences,
College of Medicine, Korea University, Seoul 02841, Republic of Korea
| | - Jong-Ik Hwang
- Department of Biomedical Sciences,
College of Medicine, Korea University, Seoul 02841, Republic of Korea
| | - Hong-Rae Kim
- Department of Biomedical Sciences,
College of Medicine, Korea University, Seoul 02841, Republic of Korea
| |
Collapse
|
9
|
Guterres H, Im W. CHARMM-GUI-Based Induced Fit Docking Workflow to Generate Reliable Protein-Ligand Binding Modes. J Chem Inf Model 2023; 63:4772-4779. [PMID: 37462607 PMCID: PMC10428204 DOI: 10.1021/acs.jcim.3c00416] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2023] [Indexed: 08/15/2023]
Abstract
Molecular docking is a preferred method to predict ligand binding modes and their binding energy to target protein receptors, which is critical in early phase structure-based drug discovery. However, there is a persistent challenge in docking that can be attributed to the induced fit effect, as receptor binding sites undergo induced fit conformational changes upon ligand binding to achieve better binding modes. In this work, based on CHARMM-GUI LBS Finder& Refiner and High-Throughput Simulator, we present a straightforward CHARMM-GUI induced fit docking (CGUI-IFD) workflow to generate reliable protein-ligand binding modes. The CGUI-IFD workflow generates an ensemble of receptor binding site conformations through ligand-binding site (LBS) refinement, runs rigid receptor docking, and performs high-throughput molecular dynamics (MD) simulations of protein-ligand complex structures in explicit solvents. The results are evaluated based on the ligand root-mean-square deviation (RMSD)-based binding stability and the molecular mechanics generalized Born surface area binding energy. For a benchmark test, we used 258 cross-docking protein-ligand pairs across 41 target proteins from the Schrodinger IFD-MD data set. The application of CGUI-IFD on this data set shows 80% success rate (within 2.5 Å RMSD from the experimental structures). We expect that the CGUI-IFD workflow can be useful to generate reliable ligand binding modes for cross-docking cases.
Collapse
Affiliation(s)
- Hugo Guterres
- Departments of Biological
Sciences, Chemistry, Bioengineering, and Computer Science and Engineering, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| | - Wonpil Im
- Departments of Biological
Sciences, Chemistry, Bioengineering, and Computer Science and Engineering, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| |
Collapse
|
10
|
Liessmann F, Künze G, Meiler J. Improving the Modeling of Extracellular Ligand Binding Pockets in RosettaGPCR for Conformational Selection. Int J Mol Sci 2023; 24:7788. [PMID: 37175495 PMCID: PMC10178219 DOI: 10.3390/ijms24097788] [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: 04/03/2023] [Revised: 04/19/2023] [Accepted: 04/22/2023] [Indexed: 05/15/2023] Open
Abstract
G protein-coupled receptors (GPCRs) are the largest class of drug targets and undergo substantial conformational changes in response to ligand binding. Despite recent progress in GPCR structure determination, static snapshots fail to reflect the conformational space of putative binding pocket geometries to which small molecule ligands can bind. In comparative modeling of GPCRs in the absence of a ligand, often a shrinking of the orthosteric binding pocket is observed. However, the exact prediction of the flexible orthosteric binding site is crucial for adequate structure-based drug discovery. In order to improve ligand docking and guide virtual screening experiments in computer-aided drug discovery, we developed RosettaGPCRPocketSize. The algorithm creates a conformational ensemble of biophysically realistic conformations of the GPCR binding pocket between the TM bundle, which is consistent with a knowledge base of expected pocket geometries. Specifically, tetrahedral volume restraints are defined based on information about critical residues in the orthosteric binding site and their experimentally observed range of Cα-Cα-distances. The output of RosettaGPCRPocketSize is an ensemble of binding pocket geometries that are filtered by energy to ensure biophysically probable arrangements, which can be used for docking simulations. In a benchmark set, pocket shrinkage observed in the default RosettaGPCR was reduced by up to 80% and the binding pocket volume range and geometric diversity were increased. Compared to models from four different GPCR homology model databases (RosettaGPCR, GPCR-Tasser, GPCR-SSFE, and GPCRdb), the here-created models showed more accurate volumes of the orthosteric pocket when evaluated with respect to the crystallographic reference structure. Furthermore, RosettaGPCRPocketSize was able to generate an improved realistic pocket distribution. However, while being superior to other homology models, the accuracy of generated model pockets was comparable to AlphaFold2 models. Furthermore, in a docking benchmark using small-molecule ligands with a higher molecular weight between 400 and 700 Da, a higher success rate in creating native-like binding poses was observed. In summary, RosettaGPCRPocketSize can generate GPCR models with realistic orthosteric pocket volumes, which are useful for structure-based drug discovery applications.
Collapse
Affiliation(s)
- Fabian Liessmann
- Institute for Drug Discovery, Medical Faculty, Leipzig University, 04103 Leipzig, Germany; (F.L.)
| | - Georg Künze
- Institute for Drug Discovery, Medical Faculty, Leipzig University, 04103 Leipzig, Germany; (F.L.)
| | - Jens Meiler
- Institute for Drug Discovery, Medical Faculty, Leipzig University, 04103 Leipzig, Germany; (F.L.)
- Department of Chemistry, Vanderbilt University, Nashville, TN 37235, USA
- Center for Structural Biology, Vanderbilt University, Nashville, TN 37235, USA
- Center for Scalable Data Analytics and Artificial Intelligence, Leipzig University, 04105 Leipzig, Germany
| |
Collapse
|
11
|
Scardino V, Di Filippo JI, Cavasotto CN. How good are AlphaFold models for docking-based virtual screening? iScience 2023; 26:105920. [PMID: 36686396 PMCID: PMC9852548 DOI: 10.1016/j.isci.2022.105920] [Citation(s) in RCA: 43] [Impact Index Per Article: 43.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2022] [Revised: 11/12/2022] [Accepted: 12/28/2022] [Indexed: 12/31/2022] Open
Abstract
A crucial component in structure-based drug discovery is the availability of high-quality three-dimensional structures of the protein target. Whenever experimental structures were not available, homology modeling has been, so far, the method of choice. Recently, AlphaFold (AF), an artificial-intelligence-based protein structure prediction method, has shown impressive results in terms of model accuracy. This outstanding success prompted us to evaluate how accurate AF models are from the perspective of docking-based drug discovery. We compared the high-throughput docking (HTD) performance of AF models with their corresponding experimental PDB structures using a benchmark set of 22 targets. The AF models showed consistently worse performance using four docking programs and two consensus techniques. Although AlphaFold shows a remarkable ability to predict protein architecture, this might not be enough to guarantee that AF models can be reliably used for HTD, and post-modeling refinement strategies might be key to increase the chances of success.
Collapse
Affiliation(s)
- Valeria Scardino
- Meton AI, Inc, Wilmington, DE 19801, USA
- Austral Institute for Applied Artificial Intelligence, Universidad Austral, Pilar, Buenos Aires, Argentina
| | - Juan I. Di Filippo
- Austral Institute for Applied Artificial Intelligence, Universidad Austral, Pilar, Buenos Aires, Argentina
- Computational Drug Design and Biomedical Informatics Laboratory, Instituto de Investigaciones en Medicina Traslacional (IIMT), Universidad Austral-CONICET, Pilar, Buenos Aires, Argentina
| | - Claudio N. Cavasotto
- Austral Institute for Applied Artificial Intelligence, Universidad Austral, Pilar, Buenos Aires, Argentina
- Computational Drug Design and Biomedical Informatics Laboratory, Instituto de Investigaciones en Medicina Traslacional (IIMT), Universidad Austral-CONICET, Pilar, Buenos Aires, Argentina
- Facultad de Ciencias Biomédicas, and Facultad de Ingeniería, Universidad Austral, Pilar, Buenos Aires, Argentina
| |
Collapse
|
12
|
Zhang J, Li H, Zhao X, Wu Q, Huang SY. Holo Protein Conformation Generation from Apo Structures by Ligand Binding Site Refinement. J Chem Inf Model 2022; 62:5806-5820. [PMID: 36342197 DOI: 10.1021/acs.jcim.2c00895] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
An important part in structure-based drug design is the selection of an appropriate protein structure. It has been revealed that a holo protein structure that contains a well-defined binding site is a much better choice than an apo structure in structure-based drug discovery. Therefore, it is valuable to obtain a holo-like protein conformation from apo structures in the case where no holo structure is available. Meeting the need, we present a robust approach to generate reliable holo-like structures from apo structures by ligand binding site refinement with restraints derived from holo templates with low homology. Our method was tested on a test set of 32 proteins from the DUD-E data set and compared with other approaches. It was shown that our method successfully refined the apo structures toward the corresponding holo conformations for 23 of 32 proteins, reducing the average all-heavy-atom RMSD of binding site residues by 0.48 Å. In addition, when evaluated against all the holo structures in the protein data bank, our method can improve the binding site RMSD for 14 of 19 cases that experience significant conformational changes. Furthermore, our refined structures also demonstrate their advantages over the apo structures in ligand binding mode predictions by both rigid docking and flexible docking and in virtual screening on the database of active and decoy ligands from the DUD-E. These results indicate that our method is effective in recovering holo-like conformations and will be valuable in structure-based drug discovery.
Collapse
Affiliation(s)
- Jinze Zhang
- School of Physics, Huazhong University of Science and Technology, Wuhan430074, Hubei, P. R. China
| | - Hao Li
- School of Physics, Huazhong University of Science and Technology, Wuhan430074, Hubei, P. R. China
| | - Xuejun Zhao
- School of Physics, Huazhong University of Science and Technology, Wuhan430074, Hubei, P. R. China
| | - Qilong Wu
- School of Physics, Huazhong University of Science and Technology, Wuhan430074, Hubei, P. R. China
| | - Sheng-You Huang
- School of Physics, Huazhong University of Science and Technology, Wuhan430074, Hubei, P. R. China
| |
Collapse
|
13
|
Zemla AT, Allen JE, Kirshner D, Lightstone FC. PDBspheres: a method for finding 3D similarities in local regions in proteins. NAR Genom Bioinform 2022; 4:lqac078. [PMID: 36225529 PMCID: PMC9549786 DOI: 10.1093/nargab/lqac078] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2022] [Revised: 08/06/2022] [Accepted: 09/29/2022] [Indexed: 11/05/2022] Open
Abstract
We present a structure-based method for finding and evaluating structural similarities in protein regions relevant to ligand binding. PDBspheres comprises an exhaustive library of protein structure regions ('spheres') adjacent to complexed ligands derived from the Protein Data Bank (PDB), along with methods to find and evaluate structural matches between a protein of interest and spheres in the library. PDBspheres uses the LGA (Local-Global Alignment) structure alignment algorithm as the main engine for detecting structural similarities between the protein of interest and template spheres from the library, which currently contains >2 million spheres. To assess confidence in structural matches, an all-atom-based similarity metric takes side chain placement into account. Here, we describe the PDBspheres method, demonstrate its ability to detect and characterize binding sites in protein structures, show how PDBspheres-a strictly structure-based method-performs on a curated dataset of 2528 ligand-bound and ligand-free crystal structures, and use PDBspheres to cluster pockets and assess structural similarities among protein binding sites of 4876 structures in the 'refined set' of the PDBbind 2019 dataset.
Collapse
Affiliation(s)
- Adam T Zemla
- To whom correspondence should be addressed. Tel: +1 925 423 5571; Fax: +1 925 423 6437;
| | - Jonathan E Allen
- Global Security Computing Applications, Lawrence Livermore National Laboratory, Livermore, CA 94550, USA
| | - Dan Kirshner
- Biosciences and Biotechnology Division, Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, CA 94550, USA
| | - Felice C Lightstone
- Biosciences and Biotechnology Division, Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, CA 94550, USA
| |
Collapse
|
14
|
Jaafar NR, Jailani N, Rahman RA, Öner ET, Murad AMA, Illias RM. Protein surface engineering and interaction studies of maltogenic amylase towards improved enzyme immobilisation. Int J Biol Macromol 2022; 213:70-82. [DOI: 10.1016/j.ijbiomac.2022.05.169] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2021] [Revised: 05/14/2022] [Accepted: 05/24/2022] [Indexed: 11/30/2022]
|
15
|
Heat Shock Protein 90 (HSP90) Inhibitors as Anticancer Medicines: A Review on the Computer-Aided Drug Discovery Approaches over the Past Five Years. COMPUTATIONAL AND MATHEMATICAL METHODS IN MEDICINE 2022; 2022:2147763. [PMID: 35685897 PMCID: PMC9173959 DOI: 10.1155/2022/2147763] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 10/12/2021] [Revised: 05/08/2022] [Accepted: 05/19/2022] [Indexed: 12/24/2022]
Abstract
Cancer is a disease caused by the uncontrolled, abnormal growth of cells in different anatomic sites. In 2018, it was predicted that the worldwide cancer burden would rise to 18.1 million new cases and 9.6 million deaths. Anticancer compounds, often known as chemotherapeutic medicines, have gained much interest in recent cancer research. These medicines work through various biological processes in targeting cells at various stages of the cell's life cycle. One of the most significant roadblocks to developing anticancer drugs is that traditional chemotherapy affects normal cells and cancer cells, resulting in substantial side effects. Recently, advancements in new drug development methodologies and the prediction of the targeted interatomic and intermolecular ligand interaction sites have been beneficial. This has prompted further research into developing and discovering novel chemical species as preferred therapeutic compounds against specific cancer types. Identifying new drug molecules with high selectivity and specificity for cancer is a prerequisite in the treatment and management of the disease. The overexpression of HSP90 occurs in patients with cancer, and the HSP90 triggers unstable harmful kinase functions, which enhance carcinogenesis. Therefore, the development of potent HSP90 inhibitors with high selectivity and specificity becomes very imperative. The activities of HSP90 as chaperones and cochaperones are complex due to the conformational dynamism, and this could be one of the reasons why no HSP90 drugs have made it beyond the clinical trials. Nevertheless, HSP90 modulations appear to be preferred due to the competitive inhibition of the targeted N-terminal adenosine triphosphate pocket. This study, therefore, presents an overview of the various computational models implored in the development of HSP90 inhibitors as anticancer medicines. We hereby suggest an extensive investigation of advanced computational modelling of the three different domains of HSP90 for potent, effective inhibitor design with minimal off-target effects.
Collapse
|
16
|
Guterres H, Park SJ, Cao Y, Im W. CHARMM-GUI Ligand Designer for Template-Based Virtual Ligand Design in a Binding Site. J Chem Inf Model 2021; 61:5336-5342. [PMID: 34757752 DOI: 10.1021/acs.jcim.1c01156] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Rational drug design involves a task of finding ligands that would bind to a specific target protein. This work presents CHARMM-GUI Ligand Designer that is an intuitive and interactive web-based tool to design virtual ligands that match the shape and chemical features of a given protein binding site. Ligand Designer provides ligand modification capabilities with 3D visualization that allow researchers to modify and redesign virtual ligands while viewing how the protein-ligand interactions are affected. Virtual ligands can also be parameterized for further molecular dynamics (MD) simulations and free energy calculations. Using 8 targets from 8 different protein classes in the directory of useful decoys, enhanced (DUD-E) data set, we show that Ligand Designer can produce similar ligands to the known active ligands in the crystal structures. Ligand Designer also produces stable protein-ligand complex structures when tested using short MD simulations. We expect that Ligand Designer can be a useful and user-friendly tool to design small molecules in any given potential ligand binding site on a protein of interest.
Collapse
Affiliation(s)
- Hugo Guterres
- Departments of Biological Sciences, Chemistry, Bioengineering, and Computer Science and Engineering, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| | - Sang-Jun Park
- Departments of Biological Sciences, Chemistry, Bioengineering, and Computer Science and Engineering, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| | - Yiwei Cao
- Departments of Biological Sciences, Chemistry, Bioengineering, and Computer Science and Engineering, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| | - Wonpil Im
- Departments of Biological Sciences, Chemistry, Bioengineering, and Computer Science and Engineering, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| |
Collapse
|
17
|
Guterres H, Park SJ, Zhang H, Im W. CHARMM-GUI LBS Finder & Refiner for Ligand Binding Site Prediction and Refinement. J Chem Inf Model 2021; 61:3744-3751. [PMID: 34296608 DOI: 10.1021/acs.jcim.1c00561] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
A protein performs its task by binding a variety of ligands in its local region that is also known as the ligand-binding-site (LBS). Therefore, accurate prediction, characterization, and refinement of LBS can facilitate protein functional annotations and structure-based drug design. In this work, we present CHARMM-GUI LBS Finder & Refiner (https://www.charmm-gui.org/input/lbsfinder) that predicts potential LBS, offers interactive features for local LBS structure analysis, and prepares various molecular dynamics (MD) systems and inputs by setting up distance restraint potentials for LBS structure refinement. LBS Finder & Refiner supports 5 different commonly used simulation programs, such as NAMD, AMBER, GROMACS, GENESIS, and OpenMM, for LBS structure refinement together with hydrogen mass repartitioning. The capability of LBS Finder & Refiner is illustrated through LBS structure predictions and refinements of 48 modeled and 20 apo benchmark target proteins. Overall, successful LBS structure predictions and refinements are seen in our benchmark tests. We hope that LBS Finder & Refiner is useful to predict, characterize, and refine potential LBS on any given protein of interest.
Collapse
Affiliation(s)
- Hugo Guterres
- Departments of Biological Sciences, Chemistry, Bioengineering, and Computer Science and Engineering, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| | - Sang-Jun Park
- Departments of Biological Sciences, Chemistry, Bioengineering, and Computer Science and Engineering, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| | - Han Zhang
- Departments of Biological Sciences, Chemistry, Bioengineering, and Computer Science and Engineering, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| | - Wonpil Im
- Departments of Biological Sciences, Chemistry, Bioengineering, and Computer Science and Engineering, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| |
Collapse
|
18
|
Singh N, Villoutreix BO. Resources and computational strategies to advance small molecule SARS-CoV-2 discovery: Lessons from the pandemic and preparing for future health crises. Comput Struct Biotechnol J 2021; 19:2537-2548. [PMID: 33936562 PMCID: PMC8074526 DOI: 10.1016/j.csbj.2021.04.059] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/04/2021] [Revised: 04/22/2021] [Accepted: 04/24/2021] [Indexed: 12/11/2022] Open
Abstract
There is an urgent need to identify new therapies that prevent SARS-CoV-2 infection and improve the outcome of COVID-19 patients. This pandemic has thus spurred intensive research in most scientific areas and in a short period of time, several vaccines have been developed. But, while the race to find vaccines for COVID-19 has dominated the headlines, other types of therapeutic agents are being developed. In this mini-review, we report several databases and online tools that could assist the discovery of anti-SARS-CoV-2 small chemical compounds and peptides. We then give examples of studies that combined in silico and in vitro screening, either for drug repositioning purposes or to search for novel bioactive compounds. Finally, we question the overall lack of discussion and plan observed in academic research in many countries during this crisis and suggest that there is room for improvement.
Collapse
Affiliation(s)
- Natesh Singh
- Université de Paris, Inserm UMR 1141 NeuroDiderot, Robert-Debré Hospital, 75019 Paris, France
| | - Bruno O. Villoutreix
- Université de Paris, Inserm UMR 1141 NeuroDiderot, Robert-Debré Hospital, 75019 Paris, France
| |
Collapse
|