1
|
Zhao M, Yu W, MacKerell AD. Enhancing SILCS-MC via GPU Acceleration and Ligand Conformational Optimization with Genetic and Parallel Tempering Algorithms. J Phys Chem B 2024. [PMID: 39031121 DOI: 10.1021/acs.jpcb.4c03045] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/22/2024]
Abstract
In the domain of computer-aided drug design, achieving precise and accurate estimates of ligand-protein binding is paramount in the context of screening extensive drug libraries and performing ligand optimization. A fundamental aspect of the SILCS (site identification by ligand competitive saturation) methodology lies in the generation of comprehensive 3D free-energy functional group affinity maps (FragMaps), encompassing the entirety of the target molecule structure. These FragMaps offer an intricate landscape of functional group affinities across the protein, bilayer, or RNA, acting as the basis for subsequent SILCS-Monte Carlo (MC) simulations wherein ligands are docked to the target molecule. To augment the efficiency and breadth of ligand sampling capabilities, we implemented an improved SILCS-MC methodology. By harnessing the parallel computing capability of GPUs, our approach facilitates concurrent calculations over multiple ligands and binding sites, markedly enhancing the computational efficiency. Moreover, the integration of a genetic algorithm (GA) with MC allows us to employ an evolutionary approach to perform ligand sampling, assuring enhanced convergence characteristics. In addition, the potential utility of parallel tempering (PT) to improve sampling was investigated. Implementation of SILCS-MC on GPU architecture is shown to accelerate the speed of SILCS-MC calculations by over 2-orders of magnitude. Use of GA and PT yield improvements over Markov-chain MC, increasing the precision of the resultant docked orientations and binding free energies, though the extent of improvements is relatively small. Accordingly, significant improvements in speed are obtained through the GPU implementation with minor improvements in the precision of the docking obtained via the tested GA and PT algorithms.
Collapse
Affiliation(s)
- Mingtian Zhao
- Computer Aided Drug Design Center, Department of Pharmaceutical Sciences, University of Maryland, School of Pharmacy, 20 Penn Street, Baltimore, Maryland 21201, United States
| | - Wenbo Yu
- Computer Aided Drug Design Center, Department of Pharmaceutical Sciences, University of Maryland, School of Pharmacy, 20 Penn Street, Baltimore, Maryland 21201, United States
| | - Alexander D MacKerell
- Computer Aided Drug Design Center, Department of Pharmaceutical Sciences, University of Maryland, School of Pharmacy, 20 Penn Street, Baltimore, Maryland 21201, United States
| |
Collapse
|
2
|
Ge Y, Pande V, Seierstad MJ, Damm-Ganamet KL. Exploring the Application of SiteMap and Site Finder for Focused Cryptic Pocket Identification. J Phys Chem B 2024; 128:6233-6245. [PMID: 38904218 DOI: 10.1021/acs.jpcb.4c00664] [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: 06/22/2024]
Abstract
The characterization of cryptic pockets has been elusive, despite substantial efforts. Computational modeling approaches, such as molecular dynamics (MD) simulations, can provide atomic-level details of binding site motions and binding pathways. However, the time scale that MD can achieve at a reasonable cost often limits its application for cryptic pocket identification. Enhanced sampling techniques can improve the efficiency of MD simulations by focused sampling of important regions of the protein, but prior knowledge of the simulated system is required to define the appropriate coordinates. In the case of a novel, unknown cryptic pocket, such information is not available, limiting the application of enhanced sampling techniques for cryptic pocket identification. In this work, we explore the ability of SiteMap and Site Finder, widely used commercial packages for pocket identification, to detect focus points on the protein and further apply other advanced computational methods. The information gained from this analysis enables the use of computational modeling, including enhanced MD sampling techniques, to explore potential cryptic binding pockets suggested by SiteMap and Site Finder. Here, we examined SiteMap and Site Finder results on 136 known cryptic pockets from a combination of the PocketMiner dataset (a recently curated set of cryptic pockets), the Cryptosite Set (a classic set of cryptic pockets), and Natural killer group 2D (NKG2D, a protein target where a cryptic pocket is confirmed). Our findings demonstrate the application of existing, well-studied tools in efficiently mapping potential regions harboring cryptic pockets.
Collapse
Affiliation(s)
- Yunhui Ge
- Computer-Aided Drug Design, Therapeutics Discovery, Janssen Research & Development, 3210 Merryfield Row, San Diego, California 92121, United States
| | - Vineet Pande
- Computer-Aided Drug Design, Therapeutics Discovery, Janssen Research & Development, Turnhoutseweg 30, 2340 Beerse, Belgium
| | - Mark J Seierstad
- Computer-Aided Drug Design, Therapeutics Discovery, Janssen Research & Development, 3210 Merryfield Row, San Diego, California 92121, United States
| | - Kelly L Damm-Ganamet
- Computer-Aided Drug Design, Therapeutics Discovery, Janssen Research & Development, 3210 Merryfield Row, San Diego, California 92121, United States
| |
Collapse
|
3
|
Chen JN, Dai B, Wu YD. Probability Density Reweighting of High-Temperature Molecular Dynamics. J Chem Theory Comput 2024; 20:4977-4985. [PMID: 38758038 DOI: 10.1021/acs.jctc.3c01423] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/18/2024]
Abstract
Molecular dynamics (MD) simulation is a popular method for elucidating the structures and functions of biomolecules. However, exploring the conformational space, especially for large systems with slow transitions, often requires enhanced sampling methods. Although conducting MD at high temperatures provides a straightforward approach, resulting conformational ensembles diverge significantly from those at low temperatures. To address this discrepancy, we propose a novel probability density-based reweighting (PDR) method. PDR exhibits robust performance across four distinct systems, including a miniprotein, a cyclic peptide, a protein loop, and a protein-peptide complex. It accurately restores the conformational distributions at high temperatures to those at low temperatures. Additionally, we apply PDR to reweight previously studied high-T MD simulations of 12 protein-peptide complexes, enabling a comprehensive investigation of the conformational space of protein-peptide complexes.
Collapse
Affiliation(s)
- Jia-Nan Chen
- Laboratory of Computational Chemistry and Drug Design, State Key Laboratory of Chemical Oncogenomics, Peking University Shenzhen Graduate School, Shenzhen 518055, China
| | - Botao Dai
- Laboratory of Computational Chemistry and Drug Design, State Key Laboratory of Chemical Oncogenomics, Peking University Shenzhen Graduate School, Shenzhen 518055, China
| | - Yun-Dong Wu
- Laboratory of Computational Chemistry and Drug Design, State Key Laboratory of Chemical Oncogenomics, Peking University Shenzhen Graduate School, Shenzhen 518055, China
- Shenzhen Bay Laboratory, Shenzhen 518132, China
- College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| |
Collapse
|
4
|
Sinha K, Basu I, Shah Z, Shah S, Chakrabarty S. Leveraging Bidirectional Nature of Allostery To Inhibit Protein-Protein Interactions (PPIs): A Case Study of PCSK9-LDLR Interaction. J Chem Inf Model 2024; 64:3923-3932. [PMID: 38615325 DOI: 10.1021/acs.jcim.4c00294] [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/16/2024]
Abstract
The protein PCSK9 (proprotein convertase subtilisin/Kexin type 9) negatively regulates the recycling of LDLR (low-density lipoprotein receptor), leading to an elevated plasma level of LDL. Inhibition of PCSK9-LDLR interaction has emerged as a promising therapeutic strategy to manage hypercholesterolemia. However, the large interaction surface area between PCSK9 and LDLR makes it challenging to identify a small molecule competitive inhibitor. An alternative strategy would be to identify distal cryptic sites as targets for allosteric inhibitors that can remotely modulate PCSK9-LDLR interaction. Using several microseconds long molecular dynamics (MD) simulations, we demonstrate that on binding with LDLR, there is a significant conformational change (population shift) in a distal loop (residues 211-222) region of PCSK9. Consistent with the bidirectional nature of allostery, we establish a clear correlation between the loop conformation and the binding affinity with LDLR. Using a thermodynamic argument, we establish that the loop conformations predominantly present in the apo state of PCSK9 would have lower LDLR binding affinity, and they would be potential targets for designing allosteric inhibitors. We elucidate the molecular origin of the allosteric coupling between this loop and the LDLR binding interface in terms of the population shift in a set of salt bridges and hydrogen bonds. Overall, our work provides a general strategy toward identifying allosteric hotspots: compare the conformational ensemble of the receptor between the apo and bound states of the protein and identify distal conformational changes, if any. The inhibitors should be designed to bind and stabilize the apo-specific conformations.
Collapse
Affiliation(s)
- Krishnendu Sinha
- Department of Chemical and Biological Sciences, S. N. Bose National Centre for Basic Sciences, Kolkata 700 106, India
| | - Ipsita Basu
- Department of Chemical and Biological Sciences, S. N. Bose National Centre for Basic Sciences, Kolkata 700 106, India
| | - Zacharia Shah
- Hingez Therapeutics Inc., 8000 Towers Crescent Drive, STE 1331, Vienna, Virginia 22182, United States
| | - Salim Shah
- Hingez Therapeutics Inc., 8000 Towers Crescent Drive, STE 1331, Vienna, Virginia 22182, United States
| | - Suman Chakrabarty
- Department of Chemical and Biological Sciences, S. N. Bose National Centre for Basic Sciences, Kolkata 700 106, India
| |
Collapse
|
5
|
Zhang M, Chen T, Lu X, Lan X, Chen Z, Lu S. G protein-coupled receptors (GPCRs): advances in structures, mechanisms, and drug discovery. Signal Transduct Target Ther 2024; 9:88. [PMID: 38594257 PMCID: PMC11004190 DOI: 10.1038/s41392-024-01803-6] [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: 08/15/2023] [Revised: 02/19/2024] [Accepted: 03/13/2024] [Indexed: 04/11/2024] Open
Abstract
G protein-coupled receptors (GPCRs), the largest family of human membrane proteins and an important class of drug targets, play a role in maintaining numerous physiological processes. Agonist or antagonist, orthosteric effects or allosteric effects, and biased signaling or balanced signaling, characterize the complexity of GPCR dynamic features. In this study, we first review the structural advancements, activation mechanisms, and functional diversity of GPCRs. We then focus on GPCR drug discovery by revealing the detailed drug-target interactions and the underlying mechanisms of orthosteric drugs approved by the US Food and Drug Administration in the past five years. Particularly, an up-to-date analysis is performed on available GPCR structures complexed with synthetic small-molecule allosteric modulators to elucidate key receptor-ligand interactions and allosteric mechanisms. Finally, we highlight how the widespread GPCR-druggable allosteric sites can guide structure- or mechanism-based drug design and propose prospects of designing bitopic ligands for the future therapeutic potential of targeting this receptor family.
Collapse
Affiliation(s)
- Mingyang Zhang
- Key Laboratory of Protection, Development and Utilization of Medicinal Resources in Liupanshan Area, Ministry of Education, Peptide & Protein Drug Research Center, School of Pharmacy, Ningxia Medical University, Yinchuan, 750004, China
- Medicinal Chemistry and Bioinformatics Center, Shanghai Jiao Tong University School of Medicine, Shanghai, 200025, China
| | - Ting Chen
- Department of Cardiology, Changzheng Hospital, Affiliated to Naval Medical University, Shanghai, 200003, China
| | - Xun Lu
- Medicinal Chemistry and Bioinformatics Center, Shanghai Jiao Tong University School of Medicine, Shanghai, 200025, China
| | - Xiaobing Lan
- Key Laboratory of Protection, Development and Utilization of Medicinal Resources in Liupanshan Area, Ministry of Education, Peptide & Protein Drug Research Center, School of Pharmacy, Ningxia Medical University, Yinchuan, 750004, China
| | - Ziqiang Chen
- Department of Orthopedics, Changhai Hospital, Affiliated to Naval Medical University, Shanghai, 200433, China.
| | - Shaoyong Lu
- Key Laboratory of Protection, Development and Utilization of Medicinal Resources in Liupanshan Area, Ministry of Education, Peptide & Protein Drug Research Center, School of Pharmacy, Ningxia Medical University, Yinchuan, 750004, China.
- Medicinal Chemistry and Bioinformatics Center, Shanghai Jiao Tong University School of Medicine, Shanghai, 200025, China.
| |
Collapse
|
6
|
Andreev G, Kovalenko M, Bozdaganyan ME, Orekhov PS. Colabind: A Cloud-Based Approach for Prediction of Binding Sites Using Coarse-Grained Simulations with Molecular Probes. J Phys Chem B 2024; 128:3211-3219. [PMID: 38514440 DOI: 10.1021/acs.jpcb.3c07853] [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: 03/23/2024]
Abstract
Binding site prediction is a crucial step in understanding protein-ligand and protein-protein interactions (PPIs) with broad implications in drug discovery and bioinformatics. This study introduces Colabind, a robust, versatile, and user-friendly cloud-based approach that employs coarse-grained molecular dynamics simulations in the presence of molecular probes, mimicking fragments of drug-like compounds. Our method has demonstrated high effectiveness when validated across a diverse range of biological targets spanning various protein classes, successfully identifying orthosteric binding sites, as well as known druggable allosteric or PPI sites, in both experimentally determined and AI-predicted protein structures, consistently placing them among the top-ranked sites. Furthermore, we suggest that careful inspection of the identified regions with a high affinity for specific probes can provide valuable insights for the development of pharmacophore hypotheses. The approach is available at https://github.com/porekhov/CG_probeMD.
Collapse
Affiliation(s)
- Georgy Andreev
- Insilico Medicine AI Ltd., Masdar City 145748, United Arab Emirates
| | - Max Kovalenko
- Division of Scientific Computing, Department of Information Technology, Uppsala University, Uppsala 752 37, Sweden
| | | | - Philipp S Orekhov
- Faculty of Biology, Shenzhen MSU-BIT University, Shenzhen 518172, China
| |
Collapse
|
7
|
Khan O, Jones G, Lazou M, Joseph-McCarthy D, Kozakov D, Beglov D, Vajda S. Expanding FTMap for Fragment-Based Identification of Pharmacophore Regions in Ligand Binding Sites. J Chem Inf Model 2024; 64:2084-2100. [PMID: 38456842 DOI: 10.1021/acs.jcim.3c01969] [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] [Indexed: 03/09/2024]
Abstract
The knowledge of ligand binding hot spots and of the important interactions within such hot spots is crucial for the design of lead compounds in the early stages of structure-based drug discovery. The computational solvent mapping server FTMap can reliably identify binding hot spots as consensus clusters, free energy minima that bind a variety of organic probe molecules. However, in its current implementation, FTMap provides limited information on regions within the hot spots that tend to interact with specific pharmacophoric features of potential ligands. E-FTMap is a new server that expands on the original FTMap protocol. E-FTMap uses 119 organic probes, rather than the 16 in the original FTMap, to exhaustively map binding sites, and identifies pharmacophore features as atomic consensus sites where similar chemical groups bind. We validate E-FTMap against a set of 109 experimentally derived structures of fragment-lead pairs, finding that highly ranked pharmacophore features overlap with the corresponding atoms in both fragments and lead compounds. Additionally, comparisons of mapping results to ensembles of bound ligands reveal that pharmacophores generated with E-FTMap tend to sample highly conserved protein-ligand interactions. E-FTMap is available as a web server at https://eftmap.bu.edu.
Collapse
Affiliation(s)
- Omeir Khan
- Department of Chemistry, Boston University, Boston, Massachusetts 02215, United States
| | - George Jones
- Department of Applied Mathematics and Statistics, Stony Brook University, Stony Brook, New York 11794, United States
| | - Maria Lazou
- Department of Biomedical Engineering, Boston University, Boston, Massachusetts 02215, United States
| | - Diane Joseph-McCarthy
- Department of Biomedical Engineering, Boston University, Boston, Massachusetts 02215, United States
| | - Dima Kozakov
- Department of Applied Mathematics and Statistics, Stony Brook University, Stony Brook, New York 11794, United States
- Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, New York 11794, United States
| | - Dmitri Beglov
- Department of Biomedical Engineering, Boston University, Boston, Massachusetts 02215, United States
- Acpharis Inc., Holliston, Massachusetts 01746, United States
| | - Sandor Vajda
- Department of Chemistry, Boston University, Boston, Massachusetts 02215, United States
- Department of Biomedical Engineering, Boston University, Boston, Massachusetts 02215, United States
| |
Collapse
|
8
|
DasGupta D, Mehrani R, Carlson HA, Sharma S. Identifying Potential Ligand Binding Sites on Glycogen Synthase Kinase 3 Using Atomistic Cosolvent Simulations. ACS APPLIED BIO MATERIALS 2024; 7:588-595. [PMID: 37141501 DOI: 10.1021/acsabm.2c01079] [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: 05/06/2023]
Abstract
Glycogen synthase kinase 3 β (GSK3β) is a serine/threonine kinase that phosphorylates several protein substrates in crucial cell signaling pathways. Owing to its therapeutic importance, there is a need to develop GSK3β inhibitors with high specificity and potency. One approach is to find small molecules that can allosterically bind to the GSK3β protein surface. We have employed fully atomistic mixed-solvent molecular dynamics (MixMD) simulations to identify three plausible allosteric sites on GSK3β that can facilitate the search for allosteric inhibitors. Our MixMD simulations narrow down the allosteric sites to precise regions on the GSK3β surface, thereby improving upon the previous predictions of the locations of these sites.
Collapse
Affiliation(s)
- Debarati DasGupta
- College of Pharmacy, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Ramin Mehrani
- Department of Mechanical Engineering, Ohio University, Athens, Ohio 45701, United States
| | - Heather A Carlson
- College of Pharmacy, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Sumit Sharma
- Department of Chemical and Biomolecular Engineering, Ohio University, Athens, Ohio 45701, United States
| |
Collapse
|
9
|
Porcelli F, Casavola AR, Grottesi A, Schiumarini D, Avaldi L. Probing the conformational dynamics of an Ago-RNA complex in water/methanol solution. Phys Chem Chem Phys 2024; 26:2497-2508. [PMID: 38170800 DOI: 10.1039/d3cp05530b] [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: 01/05/2024]
Abstract
Argonaute (Ago) proteins mediate target recognition guiding miRNA to bind complementary mRNA primarily in the seed region. However, additional pairing can occur beyond the seed, forming a supplementary duplex that can contribute to the guide-target affinity. In order to shed light on the connection, between protein-RNA interactions and miRNA-mRNA seed and supplementary duplex mobility, we carried out molecular dynamics simulations at the microsecond time-scale using a different approach compared to the ones normally used. Until now, theoretical investigations with classical MD on Ago-RNA complexes have been focused primarily on pure water solvent, which mimics the natural environment of biological molecules. Here, we explored the conformational space of a human Ago2 (hAgo2) bound to the seed + supplementary miRNA-mRNA duplex, using the solvent environment as a molecular probe. MD simulations have been performed in a mixture of water/MeOH at a molar ratio of 70 : 30 as well as in pure water for comparison. Our findings revealed that the mixed solvent promotes protein RNA association, principally enhancing salt-linkages between basic amino acid side-chains and acidic phosphates of the sugar-phosphate backbone. The primary effect registered was the restriction of supplementary duplex flexibility and the stabilization of the miRNA 3' terminus. Interestingly, we observed that the influence of the solvent appears to have almost no impact on the conformation of the seed duplex.
Collapse
Affiliation(s)
- Francesco Porcelli
- CNR-Istituto di Struttura della Materia, Area della Ricerca di Roma 1, CP 10 Monterotondo Scalo, Italy.
| | - Anna Rita Casavola
- CNR-Istituto di Struttura della Materia, Area della Ricerca di Roma 1, CP 10 Monterotondo Scalo, Italy.
| | | | - Donatella Schiumarini
- CNR-Istituto di Struttura della Materia, Area della Ricerca di Roma 1, CP 10 Monterotondo Scalo, Italy.
| | - Lorenzo Avaldi
- CNR-Istituto di Struttura della Materia, Area della Ricerca di Roma 1, CP 10 Monterotondo Scalo, Italy.
| |
Collapse
|
10
|
Kotev M, Diaz Gonzalez C. Molecular Dynamics and Other HPC Simulations for Drug Discovery. Methods Mol Biol 2024; 2716:265-291. [PMID: 37702944 DOI: 10.1007/978-1-0716-3449-3_12] [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] [Indexed: 09/14/2023]
Abstract
High performance computing (HPC) is taking an increasingly important place in drug discovery. It makes possible the simulation of complex biochemical systems with high precision in a short time, thanks to the use of sophisticated algorithms. It promotes the advancement of knowledge in fields that are inaccessible or difficult to access through experimentation and it contributes to accelerating the discovery of drugs for unmet medical needs while reducing costs. Herein, we report how computational performance has evolved over the past years, and then we detail three domains where HPC is essential. Molecular dynamics (MD) is commonly used to explore the flexibility of proteins, thus generating a better understanding of different possible approaches to modulate their activity. Modeling and simulation of biopolymer complexes enables the study of protein-protein interactions (PPI) in healthy and disease states, thus helping the identification of targets of pharmacological interest. Virtual screening (VS) also benefits from HPC to predict in a short time, among millions or billions of virtual chemical compounds, the best potential ligands that will be tested in relevant assays to start a rational drug design process.
Collapse
Affiliation(s)
- Martin Kotev
- Evotec SE, Integrated Drug Discovery, Molecular Architects, Campus Curie, Toulouse, France
| | | |
Collapse
|
11
|
Talevi A. Computer-Aided Drug Discovery and Design: Recent Advances and Future Prospects. Methods Mol Biol 2024; 2714:1-20. [PMID: 37676590 DOI: 10.1007/978-1-0716-3441-7_1] [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] [Indexed: 09/08/2023]
Abstract
Computer-aided drug discovery and design involve the use of information technologies to identify and develop, on a rational ground, chemical compounds that align a set of desired physicochemical and biological properties. In its most common form, it involves the identification and/or modification of an active scaffold (or the combination of known active scaffolds), although de novo drug design from scratch is also possible. Traditionally, the drug discovery and design processes have focused on the molecular determinants of the interactions between drug candidates and their known or intended pharmacological target(s). Nevertheless, in modern times, drug discovery and design are conceived as a particularly complex multiparameter optimization task, due to the complicated, often conflicting, property requirements.This chapter provides an updated overview of in silico approaches for identifying active scaffolds and guiding the subsequent optimization process. Recent groundbreaking advances in the field have also analyzed the integration of state-of-the-art machine learning approaches in every step of the drug discovery process (from prediction of target structure to customized molecular docking scoring functions), integration of multilevel omics data, and the use of a diversity of computational approaches to assist target validation and assess plausible binding pockets.
Collapse
Affiliation(s)
- Alan Talevi
- Laboratory of Bioactive Compound Research and Development (LIDeB), Faculty of Exact Sciences, National University of La Plata (UNLP), La Plata, Argentina.
- Argentinean National Council of Scientific and Technical Research (CONICET), La Plata, Argentina.
| |
Collapse
|
12
|
Kudo G, Yanagisawa K, Yoshino R, Hirokawa T. AAp-MSMD: Amino Acid Preference Mapping on Protein-Protein Interaction Surfaces Using Mixed-Solvent Molecular Dynamics. J Chem Inf Model 2023; 63:7768-7777. [PMID: 38085669 PMCID: PMC10751795 DOI: 10.1021/acs.jcim.3c01677] [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: 10/18/2023] [Revised: 11/28/2023] [Accepted: 11/29/2023] [Indexed: 12/26/2023]
Abstract
Peptides have attracted much attention recently owing to their well-balanced properties as drugs against protein-protein interaction (PPI) surfaces. Molecular simulation-based predictions of binding sites and amino acid residues with high affinity to PPI surfaces are expected to accelerate the design of peptide drugs. Mixed-solvent molecular dynamics (MSMD), which adds probe molecules or fragments of functional groups as solutes to the hydration model, detects the binding hotspots and cryptic sites induced by small molecules. The detection results vary depending on the type of probe molecule; thus, they provide important information for drug design. For rational peptide drug design using MSMD, we proposed MSMD with amino acid residue probes, named amino acid probe-based MSMD (AAp-MSMD), to detect hotspots and identify favorable amino acid types on protein surfaces to which peptide drugs bind. We assessed our method in terms of hotspot detection at the amino acid probe level and binding free energy prediction with amino acid probes at the PPI site for the complex structure that formed the PPI. In hotspot detection, the max-spatial probability distribution map (max-PMAP) obtained from AAp-MSMD detected the PPI site, to which each type of amino acid can bind favorably. In the binding free energy prediction using amino acid probes, ΔGFE obtained from AAp-MSMD roughly estimated the experimental binding affinities from the structure-activity relationship. AAp-MSMD, with amino acid probes, provides estimated binding sites and favorable amino acid types at the PPI site of a target protein.
Collapse
Affiliation(s)
- Genki Kudo
- Physics
Department, Graduate School of Pure and Applied Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba 305-8571, Ibaraki Japan
| | - Keisuke Yanagisawa
- Department
of Computer Science, School of Computing, Tokyo Institute of Technology, 2-12-1 Ookayama, Meguro 152-8550, Tokyo Japan
- Middle
Molecule IT-based Drug Discovery Laboratory, Tokyo Institute of Technology, 2-12-1 Ookayama, Meguro 152-8550, Tokyo Japan
| | - Ryunosuke Yoshino
- Faculty
of Medicine, University of Tsukuba, 1-1-1 Tennodai, Tsukuba 305-8575, Ibaraki Japan
- Transborder
Medical Research Center, University of Tsukuba, 1-1-1 Tennodai, Tsukuba 305-8577, Ibaraki Japan
| | - Takatsugu Hirokawa
- Faculty
of Medicine, University of Tsukuba, 1-1-1 Tennodai, Tsukuba 305-8575, Ibaraki Japan
- Transborder
Medical Research Center, University of Tsukuba, 1-1-1 Tennodai, Tsukuba 305-8577, Ibaraki Japan
| |
Collapse
|
13
|
Barrera-Téllez FJ, Prieto-Martínez FD, Hernández-Campos A, Martínez-Mayorga K, Castillo-Bocanegra R. In Silico Exploration of the Trypanothione Reductase (TryR) of L. mexicana. Int J Mol Sci 2023; 24:16046. [PMID: 38003236 PMCID: PMC10671491 DOI: 10.3390/ijms242216046] [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: 08/17/2023] [Revised: 10/23/2023] [Accepted: 10/31/2023] [Indexed: 11/26/2023] Open
Abstract
Human leishmaniasis is a neglected tropical disease which affects nearly 1.5 million people every year, with Mexico being an important endemic region. One of the major defense mechanisms of these parasites is based in the polyamine metabolic pathway, as it provides the necessary compounds for its survival. Among the enzymes in this route, trypanothione reductase (TryR), an oxidoreductase enzyme, is crucial for the Leishmania genus' survival against oxidative stress. Thus, it poses as an attractive drug target, yet due to the size and features of its catalytic pocket, modeling techniques such as molecular docking focusing on that region is not convenient. Herein, we present a computational study using several structure-based approaches to assess the druggability of TryR from L. mexicana, the predominant Leishmania species in Mexico, beyond its catalytic site. Using this consensus methodology, three relevant pockets were found, of which the one we call σ-site promises to be the most favorable one. These findings may help the design of new drugs of trypanothione-related diseases.
Collapse
Affiliation(s)
- Francisco J. Barrera-Téllez
- Departamento de Farmacia, Facultad de Química, Universidad Nacional Autónoma de México, Ciudad de México 04510, Mexico
| | - Fernando D. Prieto-Martínez
- Instituto de Química, Unidad Mérida, Universidad Nacional Autónoma de México, Carretera Mérida-Tetiz, Km. 4.5, Ucú 97357, Mexico
| | - Alicia Hernández-Campos
- Departamento de Farmacia, Facultad de Química, Universidad Nacional Autónoma de México, Ciudad de México 04510, Mexico
| | - Karina Martínez-Mayorga
- Instituto de Investigaciones en Matemáticas Aplicadas y en Sistemas, Unidad Mérida, Universidad Nacional Autónoma de México, Sierra Papacal, Mérida 97302, Mexico
| | - Rafael Castillo-Bocanegra
- Departamento de Farmacia, Facultad de Química, Universidad Nacional Autónoma de México, Ciudad de México 04510, Mexico
| |
Collapse
|
14
|
Yang WC, Gong DH, Hong Wu, Gao YY, Hao GF. Grasping cryptic binding sites to neutralize drug resistance in the field of anticancer. Drug Discov Today 2023; 28:103705. [PMID: 37453458 DOI: 10.1016/j.drudis.2023.103705] [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: 02/27/2023] [Revised: 06/09/2023] [Accepted: 07/10/2023] [Indexed: 07/18/2023]
Abstract
Drug resistance is a significant obstacle to successful cancer treatment. The utilization and development of cryptic binding sites (CBSs) in proteins involved in cancer-related drug-resistance (CRDR) could help to overcome that drug resistance. However, there is no comprehensive review of the successful use of CBSs in addressing CRDR. Here, we have systematically summarized and analyzed the opportunities and challenges of using CBSs in addressing CRDR and revealed the key role that CBSs have in targeting CRDR. First, we have identified the CRDR targets and the corresponding CBSs. Second, we discuss the mechanisms by which CBSs can overcome CRDR. Finally, we have provided examples of successful CBS applications in addressing CRDR. We hope that this approach will provide guidance to biologists and chemists in effectively utilizing CBSs for the development of new drugs to alleviate CRDR.
Collapse
Affiliation(s)
- Wei-Cheng Yang
- National Key Laboratory of Green Pesticide, Key Laboratory of Green Pesticide and Agricultural Bioengineering, Ministry of Education, Center for Research and Development of Fine Chemicals, Guizhou University, Guiyang 550025, China
| | - Dao-Hong Gong
- National Key Laboratory of Green Pesticide, Key Laboratory of Green Pesticide and Agricultural Bioengineering, Ministry of Education, Center for Research and Development of Fine Chemicals, Guizhou University, Guiyang 550025, China
| | - Hong Wu
- National Key Laboratory of Green Pesticide, Key Laboratory of Green Pesticide and Agricultural Bioengineering, Ministry of Education, Center for Research and Development of Fine Chemicals, Guizhou University, Guiyang 550025, China
| | - Yang-Yang Gao
- National Key Laboratory of Green Pesticide, Key Laboratory of Green Pesticide and Agricultural Bioengineering, Ministry of Education, Center for Research and Development of Fine Chemicals, Guizhou University, Guiyang 550025, China.
| | - Ge-Fei Hao
- National Key Laboratory of Green Pesticide, Key Laboratory of Green Pesticide and Agricultural Bioengineering, Ministry of Education, Center for Research and Development of Fine Chemicals, Guizhou University, Guiyang 550025, China; National Key Laboratory of Green Pesticide, Central China Normal University, Wuhan 430079, China.
| |
Collapse
|
15
|
La Sala G, Pfleger C, Käck H, Wissler L, Nevin P, Böhm K, Janet JP, Schimpl M, Stubbs CJ, De Vivo M, Tyrchan C, Hogner A, Gohlke H, Frolov AI. Combining structural and coevolution information to unveil allosteric sites. Chem Sci 2023; 14:7057-7067. [PMID: 37389247 PMCID: PMC10306073 DOI: 10.1039/d2sc06272k] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/14/2022] [Accepted: 06/02/2023] [Indexed: 07/01/2023] Open
Abstract
Understanding allosteric regulation in biomolecules is of great interest to pharmaceutical research and computational methods emerged during the last decades to characterize allosteric coupling. However, the prediction of allosteric sites in a protein structure remains a challenging task. Here, we integrate local binding site information, coevolutionary information, and information on dynamic allostery into a structure-based three-parameter model to identify potentially hidden allosteric sites in ensembles of protein structures with orthosteric ligands. When tested on five allosteric proteins (LFA-1, p38-α, GR, MAT2A, and BCKDK), the model successfully ranked all known allosteric pockets in the top three positions. Finally, we identified a novel druggable site in MAT2A confirmed by X-ray crystallography and SPR and a hitherto unknown druggable allosteric site in BCKDK validated by biochemical and X-ray crystallography analyses. Our model can be applied in drug discovery to identify allosteric pockets.
Collapse
Affiliation(s)
- Giuseppina La Sala
- Medicinal Chemistry, Research and Early Development, Cardiovascular, Renal and Metabolism (CVRM), BioPharmaceuticals R&D, AstraZeneca Gothenburg Sweden
| | - Christopher Pfleger
- Mathematisch-Naturwissenschaftliche Fakultät, Institut für Pharmazeutische und Medizinische Chemie, Heinrich-Heine-Universität Düsseldorf 40225 Düsseldorf Germany
| | - Helena Käck
- Mechanistic and Structural Biology, Discovery Sciences, BioPharmaceuticals R&D, AstraZeneca Gothenburg Sweden
| | - Lisa Wissler
- Mechanistic and Structural Biology, Discovery Sciences, BioPharmaceuticals R&D, AstraZeneca Gothenburg Sweden
| | - Philip Nevin
- Discovery Biology, Discovery Sciences, BioPharmaceuticals R&D, AstraZeneca Gothenburg Sweden
| | - Kerstin Böhm
- Discovery Biology, Discovery Sciences, BioPharmaceuticals R&D, AstraZeneca Gothenburg Sweden
| | - Jon Paul Janet
- Medicinal Chemistry, Research and Early Development, Cardiovascular, Renal and Metabolism (CVRM), BioPharmaceuticals R&D, AstraZeneca Gothenburg Sweden
| | - Marianne Schimpl
- Mechanistic and Structural Biology, Discovery Sciences, BioPharmaceuticals R&D, AstraZeneca Cambridge UK
| | - Christopher J Stubbs
- Mechanistic and Structural Biology, Discovery Sciences, BioPharmaceuticals R&D, AstraZeneca Cambridge UK
| | - Marco De Vivo
- Laboratory of Molecular Modeling and Drug Design, Istituto Italiano di Tecnologia Via Morego 30 16163 Genoa Italy
| | - Christian Tyrchan
- Medicinal Chemistry, Research and Early Development, Respiratory & Immunology (R&I), BioPharmaceuticals R&D, AstraZeneca Gothenburg Sweden
| | - Anders Hogner
- Medicinal Chemistry, Research and Early Development, Cardiovascular, Renal and Metabolism (CVRM), BioPharmaceuticals R&D, AstraZeneca Gothenburg Sweden
| | - Holger Gohlke
- Mathematisch-Naturwissenschaftliche Fakultät, Institut für Pharmazeutische und Medizinische Chemie, Heinrich-Heine-Universität Düsseldorf 40225 Düsseldorf Germany
- John von Neumann Institute for Computing (NIC), Jülich Supercomputing Centre (JSC), Institute of Biological Information Processing (IBI-7: Structural Biochemistry), Institute of Bio- and Geosciences (IBG-4: Bioinformatics) Forschungszentrum Jülich GmbH 52425 Jülich Germany
| | - Andrey I Frolov
- Medicinal Chemistry, Research and Early Development, Cardiovascular, Renal and Metabolism (CVRM), BioPharmaceuticals R&D, AstraZeneca Gothenburg Sweden
| |
Collapse
|
16
|
Chan WKB, Carlson HA, Traynor JR. Application of Mixed-Solvent Molecular Dynamics Simulations for Prediction of Allosteric Sites on G Protein-Coupled Receptors. Mol Pharmacol 2023; 103:274-285. [PMID: 36868791 PMCID: PMC10166447 DOI: 10.1124/molpharm.122.000612] [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: 08/16/2022] [Revised: 12/30/2022] [Accepted: 01/23/2023] [Indexed: 03/05/2023] Open
Abstract
The development of small molecule allosteric modulators acting at G protein-coupled receptors (GPCRs) is becoming increasingly attractive. Such compounds have advantages over traditional drugs acting at orthosteric sites on these receptors, in particular target specificity. However, the number and locations of druggable allosteric sites within most clinically relevant GPCRs are unknown. In the present study, we describe the development and application of a mixed-solvent molecular dynamics (MixMD)-based method for the identification of allosteric sites on GPCRs. The method employs small organic probes with druglike qualities to identify druggable hotspots in multiple replicate short-timescale simulations. As proof of principle, we first applied the method retrospectively to a test set of five GPCRs (cannabinoid receptor type 1, C-C chemokine receptor type 2, M2 muscarinic receptor, P2Y purinoceptor 1, and protease-activated receptor 2) with known allosteric sites in diverse locations. This resulted in the identification of the known allosteric sites on these receptors. We then applied the method to the μ-opioid receptor. Several allosteric modulators for this receptor are known, although the binding sites for these modulators are not known. The MixMD-based method revealed several potential allosteric sites on the mu-opioid receptor. Implementation of the MixMD-based method should aid future efforts in the structure-based drug design of drugs targeting allosteric sites on GPCRs. SIGNIFICANCE STATEMENT: Allosteric modulation of G protein-coupled receptors (GPCRs) has the potential to provide more selective drugs. However, there are limited structures of GPCRs bound to allosteric modulators, and obtaining such structures is problematic. Current computational methods utilize static structures and therefore may not identify hidden or cryptic sites. Here we describe the use of small organic probes and molecular dynamics to identify druggable allosteric hotspots on GPCRs. The results reinforce the importance of protein dynamics in allosteric site identification.
Collapse
Affiliation(s)
- Wallace K B Chan
- Department of Pharmacology and Edward F. Domino Research Center (W.K.B.C., J.R.T.) and Department of Medicinal Chemistry (H.A.C., J.R.T.), University of Michigan, Ann Arbor, Michigan
| | - Heather A Carlson
- Department of Pharmacology and Edward F. Domino Research Center (W.K.B.C., J.R.T.) and Department of Medicinal Chemistry (H.A.C., J.R.T.), University of Michigan, Ann Arbor, Michigan
| | - John R Traynor
- Department of Pharmacology and Edward F. Domino Research Center (W.K.B.C., J.R.T.) and Department of Medicinal Chemistry (H.A.C., J.R.T.), University of Michigan, Ann Arbor, Michigan
| |
Collapse
|
17
|
Masters MR, Mahmoud AH, Wei Y, Lill MA. Deep Learning Model for Efficient Protein-Ligand Docking with Implicit Side-Chain Flexibility. J Chem Inf Model 2023; 63:1695-1707. [PMID: 36916514 DOI: 10.1021/acs.jcim.2c01436] [Citation(s) in RCA: 15] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/16/2023]
Abstract
Protein-ligand docking is an essential tool in structure-based drug design with applications ranging from virtual high-throughput screening to pose prediction for lead optimization. Most docking programs for pose prediction are optimized for redocking to an existing cocrystallized protein structure, ignoring protein flexibility. In real-world drug design applications, however, protein flexibility is an essential feature of the ligand-binding process. Flexible protein-ligand docking still remains a significant challenge to computational drug design. To target this challenge, we present a deep learning (DL) model for flexible protein-ligand docking based on the prediction of an intermolecular Euclidean distance matrix (EDM), making the typical use of iterative search algorithms obsolete. The model was trained on a large-scale data set of protein-ligand complexes and evaluated on independent test sets. Our model generates high quality poses for a diverse set of protein and ligand structures and outperforms comparable docking methods.
Collapse
Affiliation(s)
- Matthew R Masters
- Department of Pharmaceutical Sciences, University of Basel, Klingelbergstrasse 50, 4056 Basel, Switzerland
| | - Amr H Mahmoud
- Department of Pharmaceutical Sciences, University of Basel, Klingelbergstrasse 50, 4056 Basel, Switzerland
| | - Yao Wei
- Department of Pharmaceutical Sciences, University of Basel, Klingelbergstrasse 50, 4056 Basel, Switzerland
| | - Markus A Lill
- Department of Pharmaceutical Sciences, University of Basel, Klingelbergstrasse 50, 4056 Basel, Switzerland
| |
Collapse
|
18
|
Is the Stalk of the SARS-CoV-2 Spike Protein Druggable? Viruses 2022; 14:v14122789. [PMID: 36560795 PMCID: PMC9786045 DOI: 10.3390/v14122789] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2022] [Revised: 11/28/2022] [Accepted: 12/12/2022] [Indexed: 12/23/2022] Open
Abstract
The spike protein is key to SARS-CoV-2 high infectivity because it facilitates the receptor binding domain (RBD) encounter with ACE2. As targeting subunit S1 has not yet delivered an ACE2-binding inhibitor, we have assessed the druggability of the conserved segment of the spike protein stalk within subunit S2 by means of an integrated computational approach that combines the molecular docking of an optimized library of fragments with high-throughput molecular dynamics simulations. The high propensity of the spike protein to mutate in key regions that are responsible for the recognition of the human angiotensin-converting enzyme 2 (hACE2) or for the recognition of antibodies, has made subunit S1 of the spike protein difficult to target. Despite the inherent flexibility of the stalk region, our results suggest two hidden interhelical binding sites, whose accessibility is only partially hampered by glycan residues.
Collapse
|
19
|
Discovery of the Cryptic Sites of SARS-CoV-2 Papain-like Protease and Analysis of Its Druggability. Int J Mol Sci 2022; 23:ijms231911265. [PMID: 36232570 PMCID: PMC9569941 DOI: 10.3390/ijms231911265] [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: 07/30/2022] [Revised: 08/30/2022] [Accepted: 09/02/2022] [Indexed: 11/17/2022] Open
Abstract
In late 2019, a new coronavirus (CoV) caused the outbreak of a deadly respiratory disease, resulting in the COVID-19 pandemic. In view of the ongoing pandemic, there is an immediate need to find drugs to treat patients. SARS-CoV-2 papain-like cysteine protease (PLpro) not only plays an important role in the pathogenesis of the virus but is also a target protein for the development of inhibitor drugs. Therefore, to develop targeted inhibitors, it is necessary to analyse and verify PLpro sites and explore whether there are other cryptic binding pockets with better activity. In this study, first, we detected the site of the whole PLpro protein by sitemap of Schrödinger (version 2018), the cavity of LigBuilder V3, and DeepSite, and roughly judged the possible activated binding site area. Then, we used the mixed solvent dynamics simulation (MixMD) of probe molecules to induce conformational changes in the protein to find the possible cryptic active sites. Finally, the TRAPP method was used to predict the druggability of cryptic pockets and analyse the changes in the physicochemical properties of residues around these sites. This work will help promote the research of SARS-CoV-2 PLpro inhibitors.
Collapse
|
20
|
A pocket-based 3D molecule generative model fueled by experimental electron density. Sci Rep 2022; 12:15100. [PMID: 36068257 PMCID: PMC9448726 DOI: 10.1038/s41598-022-19363-6] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2022] [Accepted: 08/29/2022] [Indexed: 11/08/2022] Open
Abstract
We report for the first time the use of experimental electron density (ED) as training data for the generation of drug-like three-dimensional molecules based on the structure of a target protein pocket. Similar to a structural biologist building molecules based on their ED, our model functions with two main components: a generative adversarial network (GAN) to generate the ligand ED in the input pocket and an ED interpretation module for molecule generation. The model was tested on three targets: a kinase (hematopoietic progenitor kinase 1), protease (SARS-CoV-2 main protease), and nuclear receptor (vitamin D receptor), and evaluated with a reference dataset composed of over 8000 compounds that have their activities reported in the literature. The evaluation considered the chemical validity, chemical space distribution-based diversity, and similarity with reference active compounds concerning the molecular structure and pocket-binding mode. Our model can generate molecules with similar structures to classical active compounds and novel compounds sharing similar binding modes with active compounds, making it a promising tool for library generation supporting high-throughput virtual screening. The ligand ED generated can also be used to support fragment-based drug design. Our model is available as an online service to academic users via https://edmg.stonewise.cn/#/create .
Collapse
|
21
|
Wakefield AE, Kozakov D, Vajda S. Mapping the binding sites of challenging drug targets. Curr Opin Struct Biol 2022; 75:102396. [PMID: 35636004 PMCID: PMC9790766 DOI: 10.1016/j.sbi.2022.102396] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/17/2022] [Revised: 04/20/2022] [Accepted: 04/25/2022] [Indexed: 02/03/2023]
Abstract
An increasing number of medically important proteins are challenging drug targets because their binding sites are too shallow or too polar, are cryptic and thus not detectable without a bound ligand or located in a protein-protein interface. While such proteins may not bind druglike small molecules with sufficiently high affinity, they are frequently druggable using novel therapeutic modalities. The need for such modalities can be determined by experimental or computational fragment based methods. Computational mapping by mixed solvent molecular dynamics simulations or the FTMap server can be used to determine binding hot spots. The strength and location of the hot spots provide very useful information for selecting potentially successful approaches to drug discovery.
Collapse
Affiliation(s)
- Amanda E. Wakefield
- Department of Biomedical Engineering, Boston University, Boston, Massachusetts 02215,Department of Chemistry, Boston University, Boston, Massachusetts 02215
| | - Dima Kozakov
- Department of Applied Mathematics and Statistics, Stony Brook University, Stony Brook, New York, USA,Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, New York, USA NY, USA
| | - Sandor Vajda
- Department of Biomedical Engineering, Boston University, Boston, Massachusetts 02215,Department of Chemistry, Boston University, Boston, Massachusetts 02215
| |
Collapse
|
22
|
Mayol GF, Defelipe LA, Arcon JP, Turjanski AG, Marti MA. Solvent Sites Improve Docking Performance of Protein–Protein Complexes and Protein–Protein Interface-Targeted Drugs. J Chem Inf Model 2022; 62:3577-3588. [DOI: 10.1021/acs.jcim.2c00264] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Gonzalo F. Mayol
- Departamento de Química Biológica, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires (FCEyN-UBA) e Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales (IQUIBICEN) CONICET, Pabellòn 2 de Ciudad Universitaria, Ciudad de Buenos Aires C1428EHA, Argentina
| | - Lucas A. Defelipe
- Departamento de Química Biológica, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires (FCEyN-UBA) e Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales (IQUIBICEN) CONICET, Pabellòn 2 de Ciudad Universitaria, Ciudad de Buenos Aires C1428EHA, Argentina
- European Molecular Biology Laboratory - Hamburg Unit, Notkestrasse 85, Hamburg 22607, Germany
| | - Juan Pablo Arcon
- Departamento de Química Biológica, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires (FCEyN-UBA) e Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales (IQUIBICEN) CONICET, Pabellòn 2 de Ciudad Universitaria, Ciudad de Buenos Aires C1428EHA, Argentina
- Institute for Research in Biomedicine (IRB), 08028 Barcelona, Spain
- The Barcelona Institute of Science and Technology, 08036 Barcelona, Spain
| | - Adrian G. Turjanski
- Departamento de Química Biológica, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires (FCEyN-UBA) e Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales (IQUIBICEN) CONICET, Pabellòn 2 de Ciudad Universitaria, Ciudad de Buenos Aires C1428EHA, Argentina
| | - Marcelo A. Marti
- Departamento de Química Biológica, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires (FCEyN-UBA) e Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales (IQUIBICEN) CONICET, Pabellòn 2 de Ciudad Universitaria, Ciudad de Buenos Aires C1428EHA, Argentina
| |
Collapse
|
23
|
Lal Gupta P, Carlson HA. Cosolvent Simulations with Fragment-Bound Proteins Identify Hot Spots to Direct Lead Growth. J Chem Theory Comput 2022; 18:3829-3844. [PMID: 35533286 DOI: 10.1021/acs.jctc.1c01054] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
In drug design, chemical groups are sequentially added to improve a weak-binding fragment into a tight-binding lead molecule. Often, the direction to make these additions is unclear, and there are numerous chemical modifications to choose. Lead development can be guided by crystal structures of the fragment-bound protein, but this alone is unable to capture structural changes like closing or opening of the binding site and any side-chain movements. Accounting for adaptation of the site requires a dynamic approach. Here, we use molecular dynamics calculations of small organic solvents with protein-fragment pairs to reveal the nearest "hot spots". These close hot spots show the direction to make appropriate additions and suggest types of chemical modifications that could improve binding affinity. Mixed-solvent molecular dynamics (MixMD) is a cosolvent simulation technique that is well established for finding binding "hot spots" in active sites and allosteric sites of proteins. We simulated 20 fragment-bound and apo forms of key pharmaceutical targets to map out hot spots for potential lead space. Furthermore, we analyzed whether the presence of a fragment facilitates the probes' binding in the lead space, a type of binding cooperativity. To the best of our knowledge, this is the first use of cosolvent MD conducted with bound inhibitors in the simulation. Our work provides a general framework to extract molecular features of binding sites to choose chemical groups for growing lead molecules. Of the 20 systems, 17 systems were well mapped by MixMD. For the three not-mapped systems, two had lead growth out into solution away from the protein, and the third had very small modifications which indicated no nearby hot spots. Therefore, our lack of mapping in three systems was appropriate given the experimental data (true-negative cases). The simulations are run for very short time scales, making this method tractable for use in the pharmaceutical industry.
Collapse
Affiliation(s)
- Pancham Lal Gupta
- Department of Medicinal Chemistry, College of Pharmacy, 428 Church Street, Ann Arbor, Michigan 48109-1065, United States
| | - Heather A Carlson
- Department of Medicinal Chemistry, College of Pharmacy, 428 Church Street, Ann Arbor, Michigan 48109-1065, United States
| |
Collapse
|
24
|
Zang Y, Wang H, Kang Y, Zhang J, Li X, Zhang L, Yang Z, Zhang S. TAB1 binding induced p38α conformation change: an accelerated molecular dynamics simulation study. Phys Chem Chem Phys 2022; 24:10506-10513. [PMID: 35441632 DOI: 10.1039/d2cp00144f] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
p38α mitogen-activated protein kinase (MAPK) undergoes autophosphorylation induced by the binding of TGFβ-activated kinase 1 binding protein 1 (TAB1) in myocardial ischemia. Investigation of the conformational transformations in p38α triggered by TAB1 binding is motivated by the need to find selective p38α activation inhibitors to treat myocardial ischemia. Herein, the conformational transformations of p38α were studied via all-atom accelerated molecular dynamics simulations and principal component analysis. With the binding of TAB1, the conformational changes of p38α auto-activation were characterized by the movement of the activation loop (A-loop) away from the αG helix toward the αF, αE helixes and L16-loop. In addition, a diverse intermediate state with an extensional and phosphorylated A-loop different from the transition intermediate state was explored. The conformational changes, including the A-loop alpha-structure breaking and the stronger hydrogen bond network formation, are accompanied by the extension of the A-loop and more intramolecular interactions in p38α. TAB1 correlates with other regions of p38α that are distal from the TAB1-binding site, including the A-loop, αC helix, and L16-loop, which regulates the intramolecular correlation of p38α. And, the phosphorylation further enhances the correlations between the A-loop and the other regions of p38α. The correlation results imply the regulation process of p38α conformational transformations. These findings will improve our understanding of the autophosphorylation of kinase and facilitate the development of selective inhibitors for the treatment of ischemic injury.
Collapse
Affiliation(s)
- Yongjian Zang
- MOE Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, School of Physics, Xi'an Jiaotong University, Xi'an 710049, China.
| | - He Wang
- MOE Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, School of Physics, Xi'an Jiaotong University, Xi'an 710049, China.
| | - Ying Kang
- MOE Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, School of Physics, Xi'an Jiaotong University, Xi'an 710049, China.
| | - Jianwen Zhang
- MOE Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, School of Physics, Xi'an Jiaotong University, Xi'an 710049, China.
| | - Xuhua Li
- MOE Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, School of Physics, Xi'an Jiaotong University, Xi'an 710049, China.
| | - Lei Zhang
- MOE Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, School of Physics, Xi'an Jiaotong University, Xi'an 710049, China.
| | - Zhiwei Yang
- MOE Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, School of Physics, Xi'an Jiaotong University, Xi'an 710049, China.
| | - Shengli Zhang
- MOE Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, School of Physics, Xi'an Jiaotong University, Xi'an 710049, China.
| |
Collapse
|
25
|
Alvarez-Garcia D, Schmidtke P, Cubero E, Barril X. Extracting Atomic Contributions to Binding Free Energy Using Molecular Dynamics Simulations with Mixed Solvents (MDmix). Curr Drug Discov Technol 2022; 19:62-68. [PMID: 34951392 PMCID: PMC9906626 DOI: 10.2174/1570163819666211223162829] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2021] [Revised: 09/28/2021] [Accepted: 10/05/2021] [Indexed: 11/22/2022]
Abstract
BACKGROUND Mixed solvents MD (MDmix) simulations have proved to be a useful and increasingly accepted technique with several applications in structure-based drug discovery. One of the assumptions behind the methodology is the transferability of free energy values from the simulated cosolvent molecules to larger drug-like molecules. However, the binding free energy maps (ΔGbind) calculated for the different moieties of the cosolvent molecules (e.g. a hydroxyl map for the ethanol) are largely influenced by the rest of the solvent molecule and do not reflect the intrinsic affinity of the moiety in question. As such, they are hardly transferable to different molecules. METHOD To achieve transferable energies, we present here a method for decomposing the molecular binding free energy into accurate atomic contributions. RESULT We demonstrate with two qualitative visual examples how the corrected energy maps better match known binding hotspots and how they can reveal hidden hotspots with actual drug design potential. CONCLUSION Atomic decomposition of binding free energies derived from MDmix simulations provides transferable and quantitative binding free energy maps.
Collapse
Affiliation(s)
- Daniel Alvarez-Garcia
- Gain Therapeutics, Parc Cientific de Barcelona, Baldiri Reixac 10, 08029 Barcelona, Spain
| | - Peter Schmidtke
- Facultat de Farmacia, Universitat de Barcelona, Av. Joan XXIII 27-31, 08028 Barcelona, Spain;,Current address: Discngine, 79 Avenue Ledru Rollin, 75012 Paris, France;
| | - Elena Cubero
- Gain Therapeutics, Parc Cientific de Barcelona, Baldiri Reixac 10, 08029 Barcelona, Spain
| | - Xavier Barril
- Gain Therapeutics, Parc Cientific de Barcelona, Baldiri Reixac 10, 08029 Barcelona, Spain;,Facultat de Farmacia, Universitat de Barcelona, Av. Joan XXIII 27-31, 08028 Barcelona, Spain;,Catalan Institution for Research and Advanced Studies (ICREA), Passeig Lluis Companys 23, 08010 Barcelona, Spain,Address correspondence to this author at the Gain Therapeutics, Parc Cientific de Barcelona, Baldiri Reixac 10, 08029 Barcelona, Spain; E-mail:
| |
Collapse
|
26
|
Tze-Yang Ng J, Tan YS. Accelerated Ligand-Mapping Molecular Dynamics Simulations for the Detection of Recalcitrant Cryptic Pockets and Occluded Binding Sites. J Chem Theory Comput 2022; 18:1969-1981. [PMID: 35175753 DOI: 10.1021/acs.jctc.1c01177] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The identification and characterization of binding sites is a critical component of structure-based drug design (SBDD). Probe-based/cosolvent molecular dynamics (MD) methods that allow for protein flexibility have been developed to predict ligand binding sites. However, cryptic pockets that appear only upon ligand binding and occluded binding sites with no access to the solvent pose significant challenges to these methods. Here, we report the development of accelerated ligand-mapping MD (aLMMD), which combines accelerated MD with LMMD, for the detection of these challenging binding sites. The method was validated on five proteins with what we term "recalcitrant" cryptic pockets, which are deeply buried pockets that require extensive movement of the protein backbone to expose, and three proteins with occluded binding sites. In all the cases, aLMMD was able to detect and sample the binding sites. Our results suggest that aLMMD could be used as a general approach for the detection of such elusive binding sites in protein targets, thus providing valuable information for SBDD.
Collapse
Affiliation(s)
- Justin Tze-Yang Ng
- Bioinformatics Institute, Agency for Science, Technology and Research (A*STAR), 30 Biopolis Street, #07-01 Matrix, Singapore 138671, Singapore
| | - Yaw Sing Tan
- Bioinformatics Institute, Agency for Science, Technology and Research (A*STAR), 30 Biopolis Street, #07-01 Matrix, Singapore 138671, Singapore
| |
Collapse
|
27
|
DasGupta D, Chan WKB, Carlson HA. Computational Identification of Possible Allosteric Sites and Modulators of the SARS-CoV-2 Main Protease. J Chem Inf Model 2022; 62:618-626. [PMID: 35107014 DOI: 10.1021/acs.jcim.1c01223] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
In this study, we target the main protease (Mpro) of the SARS-CoV-2 virus as it is a crucial enzyme for viral replication. Herein, we report three plausible allosteric sites on Mpro that can expand structure-based drug discovery efforts for new Mpro inhibitors. To find these sites, we used mixed-solvent molecular dynamics (MixMD) simulations, an efficient computational protocol that finds binding hotspots through mapping the surface of unbound proteins with 5% cosolvents in water. We have used normal mode analysis to support our claim of allosteric control for these sites. Further, we have performed virtual screening against the sites with 361 hits from Mpro screenings available through the National Center for Advancing Translational Sciences (NCATS). We have identified the NCATS inhibitors that bind to the remote sites better than the active site of Mpro, and we propose these molecules may be allosteric regulators of the system. After identifying our sites, new X-ray crystal structures were released that show fragment molecules in the sites we found, supporting the notion that these sites are accurate and druggable.
Collapse
Affiliation(s)
- Debarati DasGupta
- Department of Medicinal Chemistry, University of Michigan, Ann Arbor, Michigan 48109-1065, United States
| | - Wallace K B Chan
- Department of Pharmacology, University of Michigan, Ann Arbor, Michigan 48109-5632, United States
| | - Heather A Carlson
- Department of Medicinal Chemistry, University of Michigan, Ann Arbor, Michigan 48109-1065, United States
| |
Collapse
|
28
|
Goel H, Hazel A, Yu W, Jo S, MacKerell AD. Application of Site-Identification by Ligand Competitive Saturation in Computer-Aided Drug Design. NEW J CHEM 2022; 46:919-932. [PMID: 35210743 PMCID: PMC8863107 DOI: 10.1039/d1nj04028f] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Abstract
Site Identification by Ligand Competitive Saturation (SILCS) is a molecular simulation approach that uses diverse small solutes in aqueous solution to obtain functional group affinity patterns of a protein or other macromolecule. This involves employing a combined Grand Canonical Monte Carlo (GCMC)-molecular dynamics (MD) method to sample the full 3D space of the protein, including deep binding pockets and interior cavities from which functional group free energy maps (FragMaps) are obtained. The information content in the maps, which include contributions from protein flexibilty and both protein and functional group desolvation contributions, can be used in many aspects of the drug discovery process. These include identification of novel ligand binding pockets, including allosteric sites, pharmacophore modeling, prediction of relative protein-ligand binding affinities for database screening and lead optimization efforts, evaluation of protein-protein interactions as well as in the formulation of biologics-based drugs including monoclonal antibodies. The present article summarizes the various tools developed in the context of the SILCS methodology and their utility in computer-aided drug design (CADD) applications, showing how the SILCS toolset can improve the drug-development process on a number of fronts with respect to both accuracy and throughput representing a new avenue of CADD applications.
Collapse
Affiliation(s)
- Himanshu Goel
- Computer Aided Drug Design Center, Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, 20, Penn St. Baltimore, Maryland 21201, United States
| | - Anthony Hazel
- Computer Aided Drug Design Center, Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, 20, Penn St. Baltimore, Maryland 21201, United States
| | - Wenbo Yu
- Computer Aided Drug Design Center, Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, 20, Penn St. Baltimore, Maryland 21201, United States
| | - Sunhwan Jo
- SilcsBio LLC, 1100 Wicomico St. Suite 323, Baltimore, MD, 21230, United States
| | - Alexander D. MacKerell
- Computer Aided Drug Design Center, Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, 20, Penn St. Baltimore, Maryland 21201, United States., SilcsBio LLC, 1100 Wicomico St. Suite 323, Baltimore, MD, 21230, United States.,, Tel: 410-706-7442, Fax: 410-706-5017
| |
Collapse
|
29
|
Kamenik AS, Linker SM, Riniker S. Enhanced sampling without borders: on global biasing functions and how to reweight them. Phys Chem Chem Phys 2022; 24:1225-1236. [PMID: 34935813 PMCID: PMC8768491 DOI: 10.1039/d1cp04809k] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2021] [Accepted: 12/14/2021] [Indexed: 12/17/2022]
Abstract
Molecular dynamics (MD) simulations are a powerful tool to follow the time evolution of biomolecular motions in atomistic resolution. However, the high computational demand of these simulations limits the timescales of motions that can be observed. To resolve this issue, so called enhanced sampling techniques are developed, which extend conventional MD algorithms to speed up the simulation process. Here, we focus on techniques that apply global biasing functions. We provide a broad overview of established enhanced sampling methods and promising new advances. As the ultimate goal is to retrieve unbiased information from biased ensembles, we also discuss benefits and limitations of common reweighting schemes. In addition to concisely summarizing critical assumptions and implications, we highlight the general application opportunities as well as uncertainties of global enhanced sampling.
Collapse
Affiliation(s)
- Anna S Kamenik
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland.
| | - Stephanie M Linker
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland.
| | - Sereina Riniker
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland.
| |
Collapse
|
30
|
Sabanés Zariquiey F, Jacoby E, Vos A, van Vlijmen HWT, Tresadern G, Harvey J. Divide and Conquer. Pocket-Opening Mixed-Solvent Simulations in the Perspective of Docking Virtual Screening Applications for Drug Discovery. J Chem Inf Model 2022; 62:533-543. [PMID: 35041430 DOI: 10.1021/acs.jcim.1c01164] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The existence of a druggable binding pocket is a prerequisite for computational drug-target interaction studies including virtual screening. Retrospective studies have shown that extended sampling methods like Markov State Modeling and mixed-solvent simulations can identify cryptic pockets relevant for drug discovery. Here, we apply a combination of mixed-solvent molecular dynamics (MD) and time-structure independent component analysis (TICA) to four retrospective case studies: NPC2, the CECR2 bromodomain, TEM-1, and MCL-1. We compare previous experimental and computational findings to our results. It is shown that the successful identification of cryptic pockets depends on the system and the cosolvent probes. We used alternative TICA internal features such as the unbiased backbone coordinates or backbone dihedrals versus biased interatomic distances. We found that in the case of NPC2, TEM-1, and MCL-1, the use of unbiased features is able to identify cryptic pockets, although in the case of the CECR2 bromodomain, more specific features are required to properly capture a pocket opening. In the perspective of virtual screening applications, it is shown how docking studies with the parent ligands depend critically on the conformational state of the targets.
Collapse
Affiliation(s)
| | - Edgar Jacoby
- Computational Chemistry, Janssen Research & Development, Turnhoutseweg 30, B-2340 Beerse, Belgium
| | - Ann Vos
- Computational Chemistry, Janssen Research & Development, Turnhoutseweg 30, B-2340 Beerse, Belgium
| | - Herman W T van Vlijmen
- Computational Chemistry, Janssen Research & Development, Turnhoutseweg 30, B-2340 Beerse, Belgium
| | - Gary Tresadern
- Computational Chemistry, Janssen Research & Development, Turnhoutseweg 30, B-2340 Beerse, Belgium
| | - Jeremy Harvey
- Department of Chemistry, KU Leuven, Celestijnenlaan 200F, 3001 Leuven, Belgium
| |
Collapse
|
31
|
Zhang C, Wu J, Chen Q, Tan H, Huang F, Guo J, Zhang X, Yu H, Shi W. Allosteric binding on nuclear receptors: Insights on screening of non-competitive endocrine-disrupting chemicals. ENVIRONMENT INTERNATIONAL 2022; 159:107009. [PMID: 34883459 DOI: 10.1016/j.envint.2021.107009] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/23/2021] [Revised: 11/23/2021] [Accepted: 11/25/2021] [Indexed: 06/13/2023]
Abstract
Endocrine-disrupting chemicals (EDCs) can compete with endogenous hormones and bind to the orthosteric site of nuclear receptors (NRs), affecting normal endocrine system function and causing severe symptoms. Recently, a series of pharmaceuticals and personal care products (PPCPs) have been discovered to bind to the allosteric sites of NRs and induce similar effects. However, it remains unclear how diverse EDCs work in this new way. Therefore, we have systematically summarized the allosteric sites and underlying mechanisms based on existing studies, mainly regarding drugs belonging to the PPCP class. Advanced methods, classified as structural biology, biochemistry and computational simulation, together with their advantages and hurdles for allosteric site recognition and mechanism insight have also been described. Furthermore, we have highlighted two available strategies for virtual screening of numerous EDCs, relying on the structural features of allosteric sites and lead compounds, respectively. We aim to provide reliable theoretical and technical support for a broader view of various allosteric interactions between EDCs and NRs, and to drive high-throughput and accurate screening of potential EDCs with non-competitive effects.
Collapse
Affiliation(s)
- Chi Zhang
- State Key Laboratory of Pollution Control and Resources Reuse, School of the Environment, Nanjing University, Nanjing 210023, Jiangsu, China; Jiangsu Province Ecology and Environment Protection Key Laboratory of Chemical Safety and Health Risk, Nanjing 210023, Jiangsu, China
| | - Jinqiu Wu
- State Key Laboratory of Pollution Control and Resources Reuse, School of the Environment, Nanjing University, Nanjing 210023, Jiangsu, China; Jiangsu Province Ecology and Environment Protection Key Laboratory of Chemical Safety and Health Risk, Nanjing 210023, Jiangsu, China
| | - Qinchang Chen
- State Key Laboratory of Pollution Control and Resources Reuse, School of the Environment, Nanjing University, Nanjing 210023, Jiangsu, China; Jiangsu Province Ecology and Environment Protection Key Laboratory of Chemical Safety and Health Risk, Nanjing 210023, Jiangsu, China
| | - Haoyue Tan
- State Key Laboratory of Pollution Control and Resources Reuse, School of the Environment, Nanjing University, Nanjing 210023, Jiangsu, China; Jiangsu Province Ecology and Environment Protection Key Laboratory of Chemical Safety and Health Risk, Nanjing 210023, Jiangsu, China
| | - Fuyan Huang
- State Key Laboratory of Pollution Control and Resources Reuse, School of the Environment, Nanjing University, Nanjing 210023, Jiangsu, China; Jiangsu Province Ecology and Environment Protection Key Laboratory of Chemical Safety and Health Risk, Nanjing 210023, Jiangsu, China
| | - Jing Guo
- State Key Laboratory of Pollution Control and Resources Reuse, School of the Environment, Nanjing University, Nanjing 210023, Jiangsu, China; Jiangsu Province Ecology and Environment Protection Key Laboratory of Chemical Safety and Health Risk, Nanjing 210023, Jiangsu, China
| | - Xiaowei Zhang
- State Key Laboratory of Pollution Control and Resources Reuse, School of the Environment, Nanjing University, Nanjing 210023, Jiangsu, China; Jiangsu Province Ecology and Environment Protection Key Laboratory of Chemical Safety and Health Risk, Nanjing 210023, Jiangsu, China
| | - Hongxia Yu
- State Key Laboratory of Pollution Control and Resources Reuse, School of the Environment, Nanjing University, Nanjing 210023, Jiangsu, China; Jiangsu Province Ecology and Environment Protection Key Laboratory of Chemical Safety and Health Risk, Nanjing 210023, Jiangsu, China
| | - Wei Shi
- State Key Laboratory of Pollution Control and Resources Reuse, School of the Environment, Nanjing University, Nanjing 210023, Jiangsu, China; Jiangsu Province Ecology and Environment Protection Key Laboratory of Chemical Safety and Health Risk, Nanjing 210023, Jiangsu, China.
| |
Collapse
|
32
|
Rehman AU, Lu S, Khan AA, Khurshid B, Rasheed S, Wadood A, Zhang J. Hidden allosteric sites and De-Novo drug design. Expert Opin Drug Discov 2021; 17:283-295. [PMID: 34933653 DOI: 10.1080/17460441.2022.2017876] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Abstract
INTRODUCTION Hidden allosteric sites are not visible in apo-crystal structures, but they may be visible in holo-structures when a certain ligand binds and maintains the ligand intended conformation. Several computational and experimental techniques have been used to investigate these hidden sites but identifying them remains a challenge. AREAS COVERED This review provides a summary of the many theoretical approaches for predicting hidden allosteric sites in disease-related proteins. Furthermore, promising cases have been thoroughly examined to reveal the hidden allosteric site and its modulator. EXPERT OPINION In the recent past, with the development in scientific techniques and bioinformatics tools, the number of drug targets for complex human diseases has significantly increased but unfortunately most of these targets are undruggable due to several reasons. Alternative strategies such as finding cryptic (hidden) allosteric sites are an attractive approach for exploitation of the discovery of new targets. These hidden sites are difficult to recognize compared to allosteric sites, mainly due to a lack of visibility in the crystal structure. In our opinion, after many years of development, MD simulations are finally becoming successful for obtaining a detailed molecular description of drug-target interaction.
Collapse
Affiliation(s)
- Ashfaq Ur Rehman
- Key Laboratory of Cell Differentiation and Apoptosis of Chinese Ministry of Education, Clinical and Fundamental Research Center, Renji Hospital, Shanghai Jiao-Tong University School of Medicine, Shanghai, China
| | - Shaoyong Lu
- Key Laboratory of Cell Differentiation and Apoptosis of Chinese Ministry of Education, Clinical and Fundamental Research Center, Renji Hospital, Shanghai Jiao-Tong University School of Medicine, Shanghai, China
| | - Abdul Aziz Khan
- Bio-X Institutes, Key Laboratory for the Genetics of Development and Neuropsychiatric Disorders (Ministry of Education), Shanghai Key Laboratory of Psychotic Disorders, and Brain Science and Technology Research Center, Institute of Psychology and Behavioral Sciences, Shanghai Jiao Tong University, Shanghai, China
| | - Beenish Khurshid
- Department of Biochemistry, Abdul Wali Khan University Mardan, Mardan, Pakistan
| | - Salman Rasheed
- National Center for Bioinformatics, Quaid-e-Azam University, Islamabad, Pakistan
| | - Abdul Wadood
- Department of Biochemistry, Abdul Wali Khan University Mardan, Mardan, Pakistan
| | - Jian Zhang
- Key Laboratory of Cell Differentiation and Apoptosis of Chinese Ministry of Education, Clinical and Fundamental Research Center, Renji Hospital, Shanghai Jiao-Tong University School of Medicine, Shanghai, China.,School of Pharmaceutical Sciences, Zhengzhou University, Zhengzhou, China
| |
Collapse
|