1
|
Pirnia A, Maqdisi R, Mittal S, Sener M, Singharoy A. Perspective on Integrative Simulations of Bioenergetic Domains. J Phys Chem B 2024; 128:3302-3319. [PMID: 38562105 DOI: 10.1021/acs.jpcb.3c07335] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/04/2024]
Abstract
Bioenergetic processes in cells, such as photosynthesis or respiration, integrate many time and length scales, which makes the simulation of energy conversion with a mere single level of theory impossible. Just like the myriad of experimental techniques required to examine each level of organization, an array of overlapping computational techniques is necessary to model energy conversion. Here, a perspective is presented on recent efforts for modeling bioenergetic phenomena with a focus on molecular dynamics simulations and its variants as a primary method. An overview of the various classical, quantum mechanical, enhanced sampling, coarse-grained, Brownian dynamics, and Monte Carlo methods is presented. Example applications discussed include multiscale simulations of membrane-wide electron transport, rate kinetics of ATP turnover from electrochemical gradients, and finally, integrative modeling of the chromatophore, a photosynthetic pseudo-organelle.
Collapse
Affiliation(s)
- Adam Pirnia
- School of Molecular Sciences, Arizona State University, Tempe, Arizona 85287-1004, United States
| | - Ranel Maqdisi
- School of Molecular Sciences, Arizona State University, Tempe, Arizona 85287-1004, United States
| | - Sumit Mittal
- VIT Bhopal University, Sehore 466114, Madhya Pradesh, India
| | - Melih Sener
- School of Molecular Sciences, Arizona State University, Tempe, Arizona 85287-1004, United States
- Beckman Institute, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| | - Abhishek Singharoy
- School of Molecular Sciences, Arizona State University, Tempe, Arizona 85287-1004, United States
| |
Collapse
|
2
|
Asim A. Approaches to Backbone Flexibility in Protein-Protein Docking. Methods Mol Biol 2024; 2780:45-68. [PMID: 38987463 DOI: 10.1007/978-1-0716-3985-6_4] [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: 07/12/2024]
Abstract
Proteins are the fundamental organic macromolecules in living systems that play a key role in a variety of biological functions including immunological detection, intracellular trafficking, and signal transduction. The docking of proteins has greatly advanced during recent decades and has become a crucial complement to experimental methods. Protein-protein docking is a helpful method for simulating protein complexes whose structures have not yet been solved experimentally. This chapter focuses on major search tactics along with various docking programs used in protein-protein docking algorithms, which include: direct search, exhaustive global search, local shape feature matching, randomized search, and broad category of post-docking approaches. As backbone flexibility predictions and interactions in high-resolution protein-protein docking remain important issues in the overall optimization context, we have put forward several methods and solutions used to handle backbone flexibility. In addition, various docking methods that are utilized for flexible backbone docking, including ATTRACT, FlexDock, FLIPDock, HADDOCK, RosettaDock, FiberDock, etc., along with their scoring functions, algorithms, advantages, and limitations are discussed. Moreover, what progress in search technology is expected, including not only the creation of new search algorithms but also the enhancement of existing ones, has been debated. As conformational flexibility is one of the most crucial factors affecting docking success, more work should be put into evaluating the conformational flexibility upon binding for a particular case in addition to developing new algorithms to replace the rigid body docking and scoring approach.
Collapse
Affiliation(s)
- Ayesha Asim
- Department of Synthesis and Chemical Technology of Pharmaceutical Substances with Computer Modeling Laboratory, Faculty of Pharmacy, Medical University of Lublin, Lublin, Poland
| |
Collapse
|
3
|
Induced fit with replica exchange improves protein complex structure prediction. PLoS Comput Biol 2022; 18:e1010124. [PMID: 35658008 PMCID: PMC9200320 DOI: 10.1371/journal.pcbi.1010124] [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: 12/22/2021] [Revised: 06/15/2022] [Accepted: 04/20/2022] [Indexed: 11/19/2022] Open
Abstract
Despite the progress in prediction of protein complexes over the last decade, recent blind protein complex structure prediction challenges revealed limited success rates (less than 20% models with DockQ score > 0.4) on targets that exhibit significant conformational change upon binding. To overcome limitations in capturing backbone motions, we developed a new, aggressive sampling method that incorporates temperature replica exchange Monte Carlo (T-REMC) and conformational sampling techniques within docking protocols in Rosetta. Our method, ReplicaDock 2.0, mimics induced-fit mechanism of protein binding to sample backbone motions across putative interface residues on-the-fly, thereby recapitulating binding-partner induced conformational changes. Furthermore, ReplicaDock 2.0 clocks in at 150-500 CPU hours per target (protein-size dependent); a runtime that is significantly faster than Molecular Dynamics based approaches. For a benchmark set of 88 proteins with moderate to high flexibility (unbound-to-bound iRMSD over 1.2 Å), ReplicaDock 2.0 successfully docks 61% of moderately flexible complexes and 35% of highly flexible complexes. Additionally, we demonstrate that by biasing backbone sampling particularly towards residues comprising flexible loops or hinge domains, highly flexible targets can be predicted to under 2 Å accuracy. This indicates that additional gains are possible when mobile protein segments are known.
Collapse
|
4
|
Afrasiabi F, Dehghanpoor R, Haspel N. Integrating Rigidity Analysis into the Exploration of Protein Conformational Pathways Using RRT* and MC. Molecules 2021; 26:molecules26082329. [PMID: 33923805 PMCID: PMC8073574 DOI: 10.3390/molecules26082329] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2021] [Revised: 04/12/2021] [Accepted: 04/13/2021] [Indexed: 11/16/2022] Open
Abstract
To understand how proteins function on a cellular level, it is of paramount importance to understand their structures and dynamics, including the conformational changes they undergo to carry out their function. For the aforementioned reasons, the study of large conformational changes in proteins has been an interest to researchers for years. However, since some proteins experience rapid and transient conformational changes, it is hard to experimentally capture the intermediate structures. Additionally, computational brute force methods are computationally intractable, which makes it impossible to find these pathways which require a search in a high-dimensional, complex space. In our previous work, we implemented a hybrid algorithm that combines Monte-Carlo (MC) sampling and RRT*, a version of the Rapidly Exploring Random Trees (RRT) robotics-based method, to make the conformational exploration more accurate and efficient, and produce smooth conformational pathways. In this work, we integrated the rigidity analysis of proteins into our algorithm to guide the search to explore flexible regions. We demonstrate that rigidity analysis dramatically reduces the run time and accelerates convergence.
Collapse
|
5
|
Harmalkar A, Gray JJ. Advances to tackle backbone flexibility in protein docking. Curr Opin Struct Biol 2020; 67:178-186. [PMID: 33360497 DOI: 10.1016/j.sbi.2020.11.011] [Citation(s) in RCA: 25] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2020] [Revised: 11/18/2020] [Accepted: 11/25/2020] [Indexed: 12/11/2022]
Abstract
Computational docking methods can provide structural models of protein-protein complexes, but protein backbone flexibility upon association often thwarts accurate predictions. In recent blind challenges, medium or high accuracy models were submitted in less than 20% of the 'difficult' targets (with significant backbone change or uncertainty). Here, we describe recent developments in protein-protein docking and highlight advances that tackle backbone flexibility. In molecular dynamics and Monte Carlo approaches, enhanced sampling techniques have reduced time-scale limitations. Internal coordinate formulations can now capture realistic motions of monomers and complexes using harmonic dynamics. And machine learning approaches adaptively guide docking trajectories or generate novel binding site predictions from deep neural networks trained on protein interfaces. These tools poise the field to break through the longstanding challenge of correctly predicting complex structures with significant conformational change.
Collapse
Affiliation(s)
- Ameya Harmalkar
- Department of Chemical and Biomolecular Engineering, Johns Hopkins University, Baltimore, MD, USA
| | - Jeffrey J Gray
- Department of Chemical and Biomolecular Engineering, Johns Hopkins University, Baltimore, MD, USA; Program in Molecular Biophysics, Institute for Nanobiotechnology, and Center for Computational Biology, Johns Hopkins University, Baltimore, MD, USA.
| |
Collapse
|
6
|
Cirauqui Diaz N, Frezza E, Martin J. Using normal mode analysis on protein structural models. How far can we go on our predictions? Proteins 2020; 89:531-543. [PMID: 33349977 DOI: 10.1002/prot.26037] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/17/2020] [Accepted: 12/12/2020] [Indexed: 01/01/2023]
Abstract
Normal mode analysis (NMA) is a fast and inexpensive approach that is largely used to gain insight into functional protein motions, and more recently to create conformations for further computational studies. However, when the protein structure is unknown, the use of computational models is necessary. Here, we analyze the capacity of NMA in internal coordinate space to predict protein motion, its intrinsic flexibility, and atomic displacements, using protein models instead of native structures, and the possibility to use it for model refinement. Our results show that NMA is quite insensitive to modeling errors, but that calculations are strictly reliable only for very accurate models. Our study also suggests that internal NMA is a more suitable tool for the improvement of structural models, and for integrating them with experimental data or in other computational techniques, such as protein docking or more refined molecular dynamics simulations.
Collapse
Affiliation(s)
- Nuria Cirauqui Diaz
- CNRS, UMR 5086 Molecular Microbiology and Structural Biochemistry, Université de Lyon, Lyon, France
| | - Elisa Frezza
- CiTCoM, CNRS, Université de Paris, Paris, France
| | - Juliette Martin
- CNRS, UMR 5086 Molecular Microbiology and Structural Biochemistry, Université de Lyon, Lyon, France
| |
Collapse
|
7
|
Zaidman D, Prilusky J, London N. PRosettaC: Rosetta Based Modeling of PROTAC Mediated Ternary Complexes. J Chem Inf Model 2020; 60:4894-4903. [PMID: 32976709 PMCID: PMC7592117 DOI: 10.1021/acs.jcim.0c00589] [Citation(s) in RCA: 85] [Impact Index Per Article: 21.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/27/2020] [Indexed: 12/22/2022]
Abstract
Proteolysis-targeting chimeras (PROTACs), which induce degradation by recruitment of an E3 ligase to a target protein, are gaining much interest as a new pharmacological modality. However, designing PROTACs is challenging. Formation of a ternary complex between the protein target, the PROTAC, and the recruited E3 ligase is considered paramount for successful degradation. A structural model of this ternary complex could in principle inform rational PROTAC design. Unfortunately, only a handful of structures are available for such complexes, necessitating tools for their modeling. We developed a combined protocol for the modeling of a ternary complex induced by a given PROTAC. Our protocol alternates between sampling of the protein-protein interaction space and the PROTAC molecule conformational space. Application of this protocol-PRosettaC-to a benchmark of known PROTAC ternary complexes results in near-native predictions, with often atomic accuracy prediction of the protein chains, as well as the PROTAC binding moieties. It allowed the modeling of a CRBN/BTK complex that recapitulated experimental results for a series of PROTACs. PRosettaC generated models may be used to design PROTACs for new targets, as well as improve PROTACs for existing targets, potentially cutting down time and synthesis efforts. To enable wide access to this protocol, we have made it available through a web server (https://prosettac.weizmann.ac.il/).
Collapse
Affiliation(s)
- Daniel Zaidman
- Department
of Organic Chemistry, The Weizmann Institute
of Science, 76100, Rehovot, Israel
| | - Jaime Prilusky
- Life
Sciences Core Facilities, Weizmann Institute
of Science, 76100, Rehovot, Israel
| | - Nir London
- Department
of Organic Chemistry, The Weizmann Institute
of Science, 76100, Rehovot, Israel
| |
Collapse
|
8
|
Rallabandi HR, Ganesan P, Kim YJ. Targeting the C-Terminal Domain Small Phosphatase 1. Life (Basel) 2020; 10:life10050057. [PMID: 32397221 PMCID: PMC7281111 DOI: 10.3390/life10050057] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2020] [Revised: 05/05/2020] [Accepted: 05/07/2020] [Indexed: 12/15/2022] Open
Abstract
The human C-terminal domain small phosphatase 1 (CTDSP1/SCP1) is a protein phosphatase with a conserved catalytic site of DXDXT/V. CTDSP1’s major activity has been identified as dephosphorylation of the 5th Ser residue of the tandem heptad repeat of the RNA polymerase II C-terminal domain (RNAP II CTD). It is also implicated in various pivotal biological activities, such as acting as a driving factor in repressor element 1 (RE-1)-silencing transcription factor (REST) complex, which silences the neuronal genes in non-neuronal cells, G1/S phase transition, and osteoblast differentiation. Recent findings have denoted that negative regulation of CTDSP1 results in suppression of cancer invasion in neuroglioma cells. Several researchers have focused on the development of regulating materials of CTDSP1, due to the significant roles it has in various biological activities. In this review, we focused on this emerging target and explored the biological significance, challenges, and opportunities in targeting CTDSP1 from a drug designing perspective.
Collapse
|
9
|
Sivakumar KC, Haixiao J, Naman CB, Sajeevan TP. Prospects of multitarget drug designing strategies by linking molecular docking and molecular dynamics to explore the protein-ligand recognition process. Drug Dev Res 2020; 81:685-699. [PMID: 32329098 DOI: 10.1002/ddr.21673] [Citation(s) in RCA: 24] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/22/2019] [Revised: 03/24/2020] [Accepted: 04/06/2020] [Indexed: 12/14/2022]
Abstract
The designing of drugs that can simultaneously affect different protein targets is one novel and promising way to treat complex diseases. Multitarget drugs act on multiple protein receptors each implicated in the same disease state, and may be considered to be more beneficial than conventional drug therapies. For example, these drugs can have improved therapeutic potency due to synergistic effects on multiple targets, as well as improved safety and resistance profiles due to the combined regulation of potential primary therapeutic targets and compensatory elements and lower dosage typically required. This review analyzes in-silico methods that facilitate multitarget drug design that facilitate the discovery and development of novel therapeutic agents. Here presented is a summary of the progress in structure-based drug discovery techniques that study the process of molecular recognition of targets and ligands, moving from static molecular docking to improved molecular dynamics approaches in multitarget drug design, and the advantages and limitations of each.
Collapse
Affiliation(s)
- Krishnankutty Chandrika Sivakumar
- National Centre for Aquatic Animal Health, Cochin University of Science and Technology, Kochi, India.,Bioinformatics Facility, Rajiv Gandhi Centre for Biotechnology, Thiruvananthapuram, India
| | - Jin Haixiao
- Li Dak Sum Marine Biopharmaceutical Research Center, Department of Marine Pharmacy, College of Food and Pharmaceutical Sciences, Ningbo University, Ningbo, China
| | - C Benjamin Naman
- Li Dak Sum Marine Biopharmaceutical Research Center, Department of Marine Pharmacy, College of Food and Pharmaceutical Sciences, Ningbo University, Ningbo, China.,Center for Marine Biotechnology and Biomedicine, Scripps Institution of Oceanography, University of California, San Diego, La Jolla, California, USA
| | - T P Sajeevan
- National Centre for Aquatic Animal Health, Cochin University of Science and Technology, Kochi, India
| |
Collapse
|
10
|
Marze NA, Roy Burman SS, Sheffler W, Gray JJ. Efficient flexible backbone protein-protein docking for challenging targets. Bioinformatics 2019; 34:3461-3469. [PMID: 29718115 DOI: 10.1093/bioinformatics/bty355] [Citation(s) in RCA: 100] [Impact Index Per Article: 20.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2017] [Accepted: 04/27/2018] [Indexed: 11/15/2022] Open
Abstract
Motivation Binding-induced conformational changes challenge current computational docking algorithms by exponentially increasing the conformational space to be explored. To restrict this search to relevant space, some computational docking algorithms exploit the inherent flexibility of the protein monomers to simulate conformational selection from pre-generated ensembles. As the ensemble size expands with increased flexibility, these methods struggle with efficiency and high false positive rates. Results Here, we develop and benchmark RosettaDock 4.0, which efficiently samples large conformational ensembles of flexible proteins and docks them using a novel, six-dimensional, coarse-grained score function. A strong discriminative ability allows an eight-fold higher enrichment of near-native candidate structures in the coarse-grained phase compared to RosettaDock 3.2. It adaptively samples 100 conformations each of the ligand and the receptor backbone while increasing computational time by only 20-80%. In local docking of a benchmark set of 88 proteins of varying degrees of flexibility, the expected success rate (defined as cases with ≥50% chance of achieving 3 near-native structures in the 5 top-ranked ones) for blind predictions after resampling is 77% for rigid complexes, 49% for moderately flexible complexes and 31% for highly flexible complexes. These success rates on flexible complexes are a substantial step forward from all existing methods. Additionally, for highly flexible proteins, we demonstrate that when a suitable conformer generation method exists, the method successfully docks the complex. Availability and implementation As a part of the Rosetta software suite, RosettaDock 4.0 is available at https://www.rosettacommons.org to all non-commercial users for free and to commercial users for a fee. Supplementary information Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Nicholas A Marze
- Department of Chemical & Biomolecular Engineering, Johns Hopkins University, Baltimore, MD, USA
| | - Shourya S Roy Burman
- Department of Chemical & Biomolecular Engineering, Johns Hopkins University, Baltimore, MD, USA
| | - William Sheffler
- Department of Biochemistry, University of Washington, Seattle, WA, USA.,Institute for Protein Design, University of Washington, Seattle, WA, USA
| | - Jeffrey J Gray
- Department of Chemical & Biomolecular Engineering, Johns Hopkins University, Baltimore, MD, USA.,Program in Molecular Biophysics, Johns Hopkins University, Baltimore, MD, USA.,Institute for NanoBioTechnology, Johns Hopkins University, Baltimore, MD, USA.,Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University, Baltimore, MD, USA
| |
Collapse
|
11
|
Weiss D, Karpiak J, Huang XP, Sassano MF, Lyu J, Roth BL, Shoichet BK. Selectivity Challenges in Docking Screens for GPCR Targets and Antitargets. J Med Chem 2018; 61:6830-6845. [PMID: 29990431 PMCID: PMC6105036 DOI: 10.1021/acs.jmedchem.8b00718] [Citation(s) in RCA: 25] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/04/2018] [Indexed: 12/14/2022]
Abstract
To investigate large library docking's ability to find molecules with joint activity against on-targets and selectivity versus antitargets, the dopamine D2 and serotonin 5-HT2A receptors were targeted, seeking selectivity against the histamine H1 receptor. In a second campaign, κ-opioid receptor ligands were sought with selectivity versus the μ-opioid receptor. While hit rates ranged from 40% to 63% against the on-targets, they were just as good against the antitargets, even though the molecules were selected for their putative lack of binding to the off-targets. Affinities, too, were often as good or better for the off-targets. Even though it was occasionally possible to find selective molecules, such as a mid-nanomolar D2/5-HT2A ligand with 21-fold selectivity versus the H1 receptor, this was the exception. Whereas false-negatives are tolerable in docking screens against on-targets, they are intolerable against antitargets; addressing this problem may demand new strategies in the field.
Collapse
Affiliation(s)
- Dahlia
R. Weiss
- Department
of Pharmaceutical Chemistry, University
of California—San Francisco, San Francisco, California 94158-2550, United States
| | - Joel Karpiak
- Department
of Pharmaceutical Chemistry, University
of California—San Francisco, San Francisco, California 94158-2550, United States
| | - Xi-Ping Huang
- Department
of Pharmacology and National Institute of Mental Health Psychoactive
Drug Screening Program, School of Medicine, University of North Carolina, Chapel Hill, North Carolina 27599, United States
| | - Maria F. Sassano
- Department
of Pharmacology and National Institute of Mental Health Psychoactive
Drug Screening Program, School of Medicine, University of North Carolina, Chapel Hill, North Carolina 27599, United States
| | - Jiankun Lyu
- Department
of Pharmaceutical Chemistry, University
of California—San Francisco, San Francisco, California 94158-2550, United States
| | - Bryan L. Roth
- Department
of Pharmacology and National Institute of Mental Health Psychoactive
Drug Screening Program, School of Medicine, University of North Carolina, Chapel Hill, North Carolina 27599, United States
| | - Brian K. Shoichet
- Department
of Pharmaceutical Chemistry, University
of California—San Francisco, San Francisco, California 94158-2550, United States
| |
Collapse
|