1
|
Kaynak BT, Krieger JM, Dudas B, Dahmani ZL, Costa MGS, Balog E, Scott AL, Doruker P, Perahia D, Bahar I. Sampling of Protein Conformational Space Using Hybrid Simulations: A Critical Assessment of Recent Methods. Front Mol Biosci 2022; 9:832847. [PMID: 35187088 PMCID: PMC8855042 DOI: 10.3389/fmolb.2022.832847] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/10/2021] [Accepted: 01/12/2022] [Indexed: 12/17/2022] Open
Abstract
Recent years have seen several hybrid simulation methods for exploring the conformational space of proteins and their complexes or assemblies. These methods often combine fast analytical approaches with computationally expensive full atomic molecular dynamics (MD) simulations with the goal of rapidly sampling large and cooperative conformational changes at full atomic resolution. We present here a systematic comparison of the utility and limits of four such hybrid methods that have been introduced in recent years: MD with excited normal modes (MDeNM), collective modes-driven MD (CoMD), and elastic network model (ENM)-based generation, clustering, and relaxation of conformations (ClustENM) as well as its updated version integrated with MD simulations (ClustENMD). We analyzed the predicted conformational spaces using each of these four hybrid methods, applied to four well-studied proteins, triosephosphate isomerase (TIM), 3-phosphoglycerate kinase (PGK), HIV-1 protease (PR) and HIV-1 reverse transcriptase (RT), which provide extensive ensembles of experimental structures for benchmarking and comparing the methods. We show that a rigorous multi-faceted comparison and multiple metrics are necessary to properly assess the differences between conformational ensembles and provide an optimal protocol for achieving good agreement with experimental data. While all four hybrid methods perform well in general, being especially useful as computationally efficient methods that retain atomic resolution, the systematic analysis of the same systems by these four hybrid methods highlights the strengths and limitations of the methods and provides guidance for parameters and protocols to be adopted in future studies.
Collapse
Affiliation(s)
- Burak T. Kaynak
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, PA, United States
| | - James M. Krieger
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, PA, United States
| | - Balint Dudas
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, PA, United States
- Laboratoire de Biologie et Pharmacologie Appliquée, Ecole Normale Supérieure Paris-Saclay, Gif-sur-Yvette, France
- Department of Biophysics and Radiation Biology, Semmelweis University, Budapest, Hungary
| | - Zakaria L. Dahmani
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, PA, United States
| | - Mauricio G. S. Costa
- Programa de Computação Científica, Vice-Presiden̂cia de Educação, Informação e Comunicação, Fundação Oswaldo Cruz, Rio de Janeiro, Brazil
| | - Erika Balog
- Department of Biophysics and Radiation Biology, Semmelweis University, Budapest, Hungary
| | - Ana Ligia Scott
- Laboratory of Bioinformatics and Computational Biology, Center of Mathematics, Computation and Cognition, Federal University of ABC-UFABC, Santo André, Brazil
| | - Pemra Doruker
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, PA, United States
- *Correspondence: Ivet Bahar, ; David Perahia, ; Pemra Doruker,
| | - David Perahia
- Laboratoire de Biologie et Pharmacologie Appliquée, Ecole Normale Supérieure Paris-Saclay, Gif-sur-Yvette, France
- *Correspondence: Ivet Bahar, ; David Perahia, ; Pemra Doruker,
| | - Ivet Bahar
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, PA, United States
- *Correspondence: Ivet Bahar, ; David Perahia, ; Pemra Doruker,
| |
Collapse
|
2
|
Jernigan RL, Sankar K, Jia K, Faraggi E, Kloczkowski A. Computational Ways to Enhance Protein Inhibitor Design. Front Mol Biosci 2021; 7:607323. [PMID: 33614705 PMCID: PMC7886686 DOI: 10.3389/fmolb.2020.607323] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/16/2020] [Accepted: 12/08/2020] [Indexed: 11/22/2022] Open
Abstract
Two new computational approaches are described to aid in the design of new peptide-based drugs by evaluating ensembles of protein structures from their dynamics and through the assessing of structures using empirical contact potential. These approaches build on the concept that conformational variability can aid in the binding process and, for disordered proteins, can even facilitate the binding of more diverse ligands. This latter consideration indicates that such a design process should be less restrictive so that multiple inhibitors might be effective. The example chosen here focuses on proteins/peptides that bind to hemagglutinin (HA) to block the large-scale conformational change for activation. Variability in the conformations is considered from sets of experimental structures, or as an alternative, from their simple computed dynamics; the set of designe peptides/small proteins from the David Baker lab designed to bind to hemagglutinin, is the large set considered and is assessed with the new empirical contact potentials.
Collapse
Affiliation(s)
- Robert L. Jernigan
- Roy J. Carver Department of Biochemistry, Biophysics and Molecular Biology, Iowa State University, Ames, IA, United States
| | - Kannan Sankar
- Roy J. Carver Department of Biochemistry, Biophysics and Molecular Biology, Iowa State University, Ames, IA, United States
| | - Kejue Jia
- Roy J. Carver Department of Biochemistry, Biophysics and Molecular Biology, Iowa State University, Ames, IA, United States
| | - Eshel Faraggi
- Research and Information Systems, LLC, Indianapolis, IN, United States
- Department of Physics, Indiana University Purdue University Indianapolis, Indianapolis, IN, United States
| | - Andrzej Kloczkowski
- Battelle Center for Mathematical Medicine, Nationwide Children's Hospital, Columbus, OH, United States
- Department of Pediatrics, The Ohio State University, Columbus, OH, United States
| |
Collapse
|
3
|
Cutinho PF, Roy J, Anand A, Cheluvaraj R, Murahari M, Chimatapu HSV. Design of metronidazole derivatives and flavonoids as potential non-nucleoside reverse transcriptase inhibitors using combined ligand- and structure-based approaches. J Biomol Struct Dyn 2019; 38:1626-1648. [DOI: 10.1080/07391102.2019.1614094] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Pretisha Flora Cutinho
- Department of Pharmaceutical Chemistry, Faculty of Pharmacy, M.S. Ramaiah University of Applied Sciences, Bangalore, India
| | - Jaydeep Roy
- Department of Pharmacology, Faculty of Pharmacy, M.S. Ramaiah University of Applied Sciences, Bangalore, India
| | - Avinash Anand
- Department of Pharmaceutical Chemistry, Faculty of Pharmacy, M.S. Ramaiah University of Applied Sciences, Bangalore, India
| | - Ravishankar Cheluvaraj
- Department of Pharmacology, Faculty of Pharmacy, M.S. Ramaiah University of Applied Sciences, Bangalore, India
| | - Manikanta Murahari
- Department of Pharmaceutical Chemistry, Faculty of Pharmacy, M.S. Ramaiah University of Applied Sciences, Bangalore, India
- Pharmacological Modelling & Simulation Centre, M.S. Ramaiah University of Applied Sciences, Bangalore, India
| | - H. S. Venkataramana Chimatapu
- Department of Pharmaceutical Chemistry, Faculty of Pharmacy, M.S. Ramaiah University of Applied Sciences, Bangalore, India
| |
Collapse
|
4
|
Loop dynamics behind the affinity of DARPins towards ERK2: Molecular dynamics simulations (MDs) and elastic network model (ENM). J Mol Liq 2019. [DOI: 10.1016/j.molliq.2018.10.157] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
|
5
|
Zhao ZW, Xie XS, Ge H. Nonequilibrium Relaxation of Conformational Dynamics Facilitates Catalytic Reaction in an Elastic Network Model of T7 DNA Polymerase. J Phys Chem B 2016; 120:2869-77. [PMID: 26918464 DOI: 10.1021/acs.jpcb.5b11002] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Nucleotide-induced conformational closing of the finger domain of DNA polymerase is crucial for its catalytic action during DNA replication. Such large-amplitude molecular motion is often not fully accessible to either direct experimental monitoring or molecular dynamics simulations. However, a coarse-grained model can offer an informative alternative, especially for probing the relationship between conformational dynamics and catalysis. Here we investigate the dynamics of T7 DNA polymerase catalysis using a Langevin-type elastic network model incorporating detailed structural information on the open conformation without the substrate bound. Such a single-parameter model remarkably captures the induced conformational dynamics of DNA polymerase upon dNTP binding, and reveals its close coupling to the advancement toward transition state along the coordinate of the target reaction, which contributes to significant lowering of the activation energy barrier. Furthermore, analysis of stochastic catalytic rates suggests that when the activation energy barrier has already been significantly lowered and nonequilibrium relaxation toward the closed form dominates the catalytic rate, one must appeal to a picture of two-dimensional free energy surface in order to account for the full spectrum of catalytic modes. Our semiquantitative study illustrates the general role of conformational dynamics in achieving transition-state stabilization, and suggests that such an elastic network model, albeit simplified, possesses the potential to furnish significant mechanistic insights into the functioning of a variety of enzymatic systems.
Collapse
Affiliation(s)
- Ziqing W Zhao
- Department of Chemistry and Chemical Biology, Harvard University , Cambridge, Massachusetts 02138, United States.,Graduate Program in Biophysics, Harvard University , Cambridge, Massachusetts 02138, United States
| | - X Sunney Xie
- Department of Chemistry and Chemical Biology, Harvard University , Cambridge, Massachusetts 02138, United States.,Graduate Program in Biophysics, Harvard University , Cambridge, Massachusetts 02138, United States.,Biodynamic Optical Imaging Center (BIOPIC), Peking University , Beijing 100871, P. R. China
| | - Hao Ge
- Biodynamic Optical Imaging Center (BIOPIC), Peking University , Beijing 100871, P. R. China.,Beijing International Center for Mathematical Research (BICMR), Peking University , Beijing 100871, P. R. China
| |
Collapse
|
6
|
Yan A, Wang Y, Kloczkowski A, Jernigan RL. Effects of protein subunits removal on the computed motions of partial 30S structures of the ribosome. J Chem Theory Comput 2008; 4:1757-1767. [PMID: 19771145 DOI: 10.1021/ct800223g] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The Anisotropic Network Model (ANM) is used to study motions of the 30S small ribosomal subunit. The effect of the absence of certain subunits on the motions of the remaining partial structures was investigated by removing one protein, pairs of proteins and selected sets of proteins at a time. Our results show that the removal of some proteins doesn't change the large-scale dynamics of the partial structures, but the removal of certain subunits does cause significant changes in motion of the remaining structure, and these changes can be reverted by the removal of other subunits, which indicates interdependence between motions of various parts of the 30S ribosomal structure. We further found that the subunits showing such interdependence have strong positive correlation of their motions, which indicates that these subunits function as a unit block in the 30S small ribosomal subunit. Dynamically interdependent subunit pairs identified in this paper are consistent with previous experimental observations that suggested dimerization of those subunits.
Collapse
Affiliation(s)
- Aimin Yan
- Laurence H. Baker Center for Bioinformatics and Biological Statistics and Department of Biochemistry, Biophysics and Molecular Biology Iowa State University, Ames, IA 50011
| | | | | | | |
Collapse
|
7
|
Jernigan RL, Kloczkowski A. Packing regularities in biological structures relate to their dynamics. METHODS IN MOLECULAR BIOLOGY (CLIFTON, N.J.) 2006; 350:251-76. [PMID: 16957327 PMCID: PMC2039702 DOI: 10.1385/1-59745-189-4:251] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 03/28/2023]
Abstract
The high packing density inside proteins leads to certain geometric regularities and also is one of the most important contributors to the high extent of cooperativity manifested by proteins in their cohesive domain motions. The orientations between neighboring nonbonded residues in proteins substantially follow the similar geometric regularities, regardless of whether the residues are on the surface or buried, a direct result of hydrophobicity forces. These orientations are relatively fixed and correspond closely to small deformations from those of the face-centered cubic lattice, which is the way in which identical spheres pack at the highest density. Packing density also is related to the extent of conservation of residues, and we show this relationship for residue packing densities by averaging over a large sample or residue packings. There are three regimes: (1) over a broad range of packing densities the relationship between sequence entropy and inverse packing density is nearly linear, (2) over a limited range of low packing densities the sequence entropy is nearly constant, and (3) at extremely low packing densities the sequence entropy is highly variable. These packing results provide important justification for the simple elastic network models that have been shown for a large number of proteins to represent protein dynamics so successfully, even when the models are extremely coarse grained. Elastic network models for polymeric chains are simple and could be combined with these protein elastic networks to represent partially denatured parts of proteins. Finally, we show results of applications of the elastic network model to study the functional motions of the ribosome, based on its known structure. These results indicate expected correlations among its components for the step-wise processing steps in protein synthesis, and suggest ways to use these elastic network models to develop more detailed mechanisms, an important possibility because most experiments yield only static structures.
Collapse
Affiliation(s)
- Robert L Jernigan
- Department of Biochemistry, Biophysics, and Molecular Biology, Laurence H. Baker Center for Bioinformatics and Biological Statistics, Iowa State University, Ames, IA, USA
| | | |
Collapse
|
8
|
Fernández A, Tawfik DS, Berkhout B, Sanders R, Kloczkowski A, Sen T, Jernigan B. Protein promiscuity: drug resistance and native functions--HIV-1 case. J Biomol Struct Dyn 2005; 22:615-24. [PMID: 15842167 DOI: 10.1080/07391102.2005.10531228] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
Abstract
The association of a drug with its target protein has the effect of blocking the protein activity and is termed a promiscuous function to distinguish from the protein's native function (Tawfik and associates, Nat. Genet. 37, 73-6, 2005). Obviously, a protein has not evolved naturally for drug association or drug resistance. Promiscuous protein functions exhibit unique traits of evolutionary adaptability, or evolvability, which is dependent on the induction of novel phenotypic traits by a small number of mutations. These mutations might have small effects on native functions, but large effects on promiscuous function; for example, an evolving protein could become increasingly drug resistant while maintaining its original function. Ariel Fernandez, in his opinion piece, notes that drug-binding "promiscuity" can hardly be dissociated from native functions; a dominant approach to drug discovery is the protein-native-substrate transition-state mimetic strategy. Thus, man-made ligands (e.g. drugs) have been successfully crafted to restrain enzymatic activity by focusing on the very same structural features that determine the native function. Using the successful inhibition of HIV-1 protease as an example, Fernandez illustrates how drug designers have employed naturally evolved features of the protein to suppress its activity. Based on these arguments, he dismisses the notion that drug binding is quintessentially promiscuous, even though in principle, proteins did not evolve to associate with man made ligands. In short, Fernandez argues that there may not be separate protein domains that one could term promiscuous domains. While acknowledging that drugs may bind promiscuously or in a native-like manner a la Fernandez, Tawfik maintains the role of evolutionary adaptation, even when a drug binds native-like. In the case of HIV-1 protease, drugs bind natively, and the initial onset of mutations results in drug resistance in addition to a dramatic decline in enzymatic activity and fitness of the virus. A chain of compensatory mutations follows this, and then the virus becomes fully fit and drug resistant. Ben Berkhout and Rogier Sanders subscribe to the evolution of new protein functions through gene duplication. With two identical protein domains, one domain can be released from a constraint imposed by the original function and it is thus free to move in sequence space toward a new function without loss of the original function. They emphasize that the forced evolution of drug-resistance differs significantly from the spontaneous evolution of an additional protein function. For instance, the latter process could proceed gradually on an evolutionary time scale, whereas the acquisition of drug-resistance is an all or nothing process for a virus, leading to the failure or success of therapy. They find no evidence to the thesis that resistance-mutations appear more rapidly in promiscuous domains than native domains. Berkhout and Sanders illustrate the genetic plasticity of HIV-1 by citing examples in which well-conserved amino acid residues of catalytic domains are forced to mutate under drug-pressure. HIV drug resistance biology is very complex. Instead of a viral protein, a drug can be targeted at a cellular protein. For example, Berkhout and Sanders claim, a drug targeted at the cellular protein CCR5 inhibits the binding of the viral envelope glycoprotein (Env) to CCR5. However, Env mutates so that it binds to the CCR5-drug complex and develops drug resistance. Interestingly, CCR5 has not evolved to bind to Env, but to a series of chemokines. Andrzej Kloczkowski, Taner Sen, and Bob Jernigan point out the importance of protein motions for binding. They believe it is likely that different ligands can bind to the diverse protein conformations sampled in the course of normal protein conformational fluctuations. They have been applying simple elastic network models to extract the motions as normal modes, which yield relatively small numbers of conformations that are useful for developing protein mechanisms; while these are typically small motions, for some proteins they can be quite large in scale. One of the major advantages of the approach is that only relatively small numbers of modes are important contributors to the overall motion -- so the approach provides a way to systematically map out a protein's motions. These models successfully represent the conformational fluctuations manifested in the crystallographic B-factors, and often suggest motions related to protein functional behaviors, such as those observed for reverse transcriptase, where two dominant hinges clearly relate to the processing steps -- one showing anti-correlation between the polymerase and ribonuclease H sites related to the translation and positioning of the nucleic acid chain, and another for opening and closing the polymerase site. Disordered proteins represent a more extreme case where the set of accessible conformations is much larger; thus they could offer up a broader range of possible binding forms. Whether evolution controls the functional motions for proteins remains little studied. Intriguingly, buried in the existing databases of protein-protein interactions may be information that can shed light on the extent of promiscuous binding among proteins themselves. Within these data there are cases where large numbers of diverse proteins have been shown to interact with a single protein; some of these could represent promiscuous protein-protein binding. Uncovering these promiscuous behaviors could be important for comprehending the details of how proteins can bind promiscuously to one another, and can exhibit even greater promiscuity in their binding to small molecules. The evolutionary routes, the dynamics of the target protein, and the many other aspects that need to be addressed while designing a drug that may dodge drug resistance, indicate the complexity and multi-disciplinary nature of the issue of drug resistance.
Collapse
Affiliation(s)
- Ariel Fernández
- Center for Computational Biology and Bioinformatics, Indiana University School of Medicine, Indianapolis, IN 46202, USA.
| | | | | | | | | | | | | |
Collapse
|