1
|
Badaczewska-Dawid A, Wróblewski K, Kurcinski M, Kmiecik S. Structure prediction of linear and cyclic peptides using CABS-flex. Brief Bioinform 2024; 25:bbae003. [PMID: 38305457 PMCID: PMC10836054 DOI: 10.1093/bib/bbae003] [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: 06/30/2023] [Revised: 12/08/2023] [Accepted: 12/28/2023] [Indexed: 02/03/2024] Open
Abstract
The structural modeling of peptides can be a useful aid in the discovery of new drugs and a deeper understanding of the molecular mechanisms of life. Here we present a novel multiscale protocol for the structure prediction of linear and cyclic peptides. The protocol combines two main stages: coarse-grained simulations using the CABS-flex standalone package and an all-atom reconstruction-optimization process using the Modeller program. We evaluated the protocol on a set of linear peptides and two sets of cyclic peptides, with cyclization through the backbone and disulfide bonds. A comparison with other state-of-the-art tools (APPTEST, PEP-FOLD, ESMFold and AlphaFold implementation in ColabFold) shows that for most cases, AlphaFold offers the highest resolution. However, CABS-flex is competitive, particularly when it comes to short linear peptides. As demonstrated, the protocol performance can be further improved by combination with the residue-residue contact prediction method or more efficient scoring. The protocol is included in the CABS-flex standalone package along with online documentation to aid users in predicting the structure of peptides and mini-proteins.
Collapse
Affiliation(s)
| | - Karol Wróblewski
- Biological and Chemical Research Center, Faculty of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| | - Mateusz Kurcinski
- Biological and Chemical Research Center, Faculty of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| | - Sebastian Kmiecik
- Biological and Chemical Research Center, Faculty of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| |
Collapse
|
2
|
Durairaj J, de Ridder D, van Dijk AD. Beyond sequence: Structure-based machine learning. Comput Struct Biotechnol J 2022; 21:630-643. [PMID: 36659927 PMCID: PMC9826903 DOI: 10.1016/j.csbj.2022.12.039] [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: 09/26/2022] [Revised: 12/21/2022] [Accepted: 12/21/2022] [Indexed: 12/31/2022] Open
Abstract
Recent breakthroughs in protein structure prediction demarcate the start of a new era in structural bioinformatics. Combined with various advances in experimental structure determination and the uninterrupted pace at which new structures are published, this promises an age in which protein structure information is as prevalent and ubiquitous as sequence. Machine learning in protein bioinformatics has been dominated by sequence-based methods, but this is now changing to make use of the deluge of rich structural information as input. Machine learning methods making use of structures are scattered across literature and cover a number of different applications and scopes; while some try to address questions and tasks within a single protein family, others aim to capture characteristics across all available proteins. In this review, we look at the variety of structure-based machine learning approaches, how structures can be used as input, and typical applications of these approaches in protein biology. We also discuss current challenges and opportunities in this all-important and increasingly popular field.
Collapse
Affiliation(s)
- Janani Durairaj
- Biozentrum, University of Basel, Basel, Switzerland
- Bioinformatics Group, Department of Plant Sciences, Wageningen University and Research, Wageningen, the Netherlands
| | - Dick de Ridder
- Bioinformatics Group, Department of Plant Sciences, Wageningen University and Research, Wageningen, the Netherlands
| | - Aalt D.J. van Dijk
- Bioinformatics Group, Department of Plant Sciences, Wageningen University and Research, Wageningen, the Netherlands
| |
Collapse
|
3
|
Verkhivker GM, Agajanian S, Oztas D, Gupta G. Computational analysis of protein stability and allosteric interaction networks in distinct conformational forms of the SARS-CoV-2 spike D614G mutant: reconciling functional mechanisms through allosteric model of spike regulation. J Biomol Struct Dyn 2022; 40:9724-9741. [PMID: 34060425 DOI: 10.1080/07391102.2021.1933594] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
In this study, we used an integrative computational approach to examine molecular mechanisms underlying functional effects of the D614G mutation by exploring atomistic modeling of the SARS-CoV-2 spike proteins as allosteric regulatory machines. We combined coarse-grained simulations, protein stability and dynamic fluctuation communication analysis with network-based community analysis to examine structures of the native and mutant SARS-CoV-2 spike proteins in different functional states. Through distance fluctuations communication analysis, we probed stability and allosteric communication propensities of protein residues in the native and mutant SARS-CoV-2 spike proteins, providing evidence that the D614G mutation can enhance long-range signaling of the allosteric spike engine. By combining functional dynamics analysis and ensemble-based alanine scanning of the SARS-CoV-2 spike proteins we found that the D614G mutation can improve stability of the spike protein in both closed and open forms, but shifting thermodynamic preferences towards the open mutant form. Our results revealed that the D614G mutation can promote the increased number of stable communities and allosteric hub centers in the open form by reorganizing and enhancing the stability of the S1-S2 inter-domain interactions and restricting mobility of the S1 regions. This study provides atomistic-based view of allosteric communications in the SARS-CoV-2 spike proteins, suggesting that the D614G mutation can exert its primary effect through allosterically induced changes on stability and communications in the residue interaction networks.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Gennady M Verkhivker
- Keck Center for Science and Engineering, Schmid College of Science and Technology, Chapman University, Orange, CA, USA.,Depatment of Biomedical and Pharmaceutical Sciences, Chapman University School of Pharmacy, Irvine, CA, USA
| | - Steve Agajanian
- Keck Center for Science and Engineering, Schmid College of Science and Technology, Chapman University, Orange, CA, USA
| | - Deniz Oztas
- Keck Center for Science and Engineering, Schmid College of Science and Technology, Chapman University, Orange, CA, USA
| | - Grace Gupta
- Keck Center for Science and Engineering, Schmid College of Science and Technology, Chapman University, Orange, CA, USA
| |
Collapse
|
4
|
Ślusarz R, Lubecka EA, Czaplewski C, Liwo A. Improvements and new functionalities of UNRES server for coarse-grained modeling of protein structure, dynamics, and interactions. Front Mol Biosci 2022; 9:1071428. [PMID: 36589235 PMCID: PMC9794589 DOI: 10.3389/fmolb.2022.1071428] [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/16/2022] [Accepted: 11/29/2022] [Indexed: 12/15/2022] Open
Abstract
In this paper we report the improvements and extensions of the UNRES server (https://unres-server.chem.ug.edu.pl) for physics-based simulations with the coarse-grained UNRES model of polypeptide chains. The improvements include the replacement of the old code with the recently optimized one and adding the recent scale-consistent variant of the UNRES force field, which performs better in the modeling of proteins with the β and the α+β structures. The scope of applications of the package was extended to data-assisted simulations with restraints from nuclear magnetic resonance (NMR) and chemical crosslink mass-spectroscopy (XL-MS) measurements. NMR restraints can be input in the NMR Exchange Format (NEF), which has become a standard. Ambiguous NMR restraints are handled without expert intervention owing to a specially designed penalty function. The server can be used to run smaller jobs directly or to prepare input data to run larger production jobs by using standalone installations of UNRES.
Collapse
Affiliation(s)
- Rafał Ślusarz
- Faculty of Chemistry, University of Gdańsk, Fahrenheit Union of Universities in Gdańsk, Gdańsk, Poland
| | - Emilia A. Lubecka
- Faculty of Electronics, Telecommunication and Informatics, Gdańsk University of Technology, Fahrenheit Union of Universities in Gdańsk, Gdańsk, Poland
| | - Cezary Czaplewski
- Faculty of Chemistry, University of Gdańsk, Fahrenheit Union of Universities in Gdańsk, Gdańsk, Poland
| | - Adam Liwo
- Faculty of Chemistry, University of Gdańsk, Fahrenheit Union of Universities in Gdańsk, Gdańsk, Poland,*Correspondence: Adam Liwo,
| |
Collapse
|
5
|
Colberg M, Schofield J. Configurational entropy, transition rates, and optimal interactions for rapid folding in coarse-grained model proteins. J Chem Phys 2022; 157:125101. [PMID: 36182418 DOI: 10.1063/5.0098612] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Under certain conditions, the dynamics of coarse-grained models of solvated proteins can be described using a Markov state model, which tracks the evolution of populations of configurations. The transition rates among states that appear in the Markov model can be determined by computing the relative entropy of states and their mean first passage times. In this paper, we present an adaptive method to evaluate the configurational entropy and the mean first passage times for linear chain models with discontinuous potentials. The approach is based on event-driven dynamical sampling in a massively parallel architecture. Using the fact that the transition rate matrix can be calculated for any choice of interaction energies at any temperature, it is demonstrated how each state's energy can be chosen such that the average time to transition between any two states is minimized. The methods are used to analyze the optimization of the folding process of two protein systems: the crambin protein and a model with frustration and misfolding. It is shown that the folding pathways for both systems are comprised of two regimes: first, the rapid establishment of local bonds, followed by the subsequent formation of more distant contacts. The state energies that lead to the most rapid folding encourage multiple pathways, and they either penalize folding pathways through kinetic traps by raising the energies of trapping states or establish an escape route from the trapping states by lowering free energy barriers to other states that rapidly reach the native state.
Collapse
Affiliation(s)
- Margarita Colberg
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto, Toronto, Ontario M5S 3H6, Canada
| | - Jeremy Schofield
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto, Toronto, Ontario M5S 3H6, Canada
| |
Collapse
|
6
|
Tomasiak L, Karch R, Schreiner W. Conformational flexibility of a free and TCR-bound pMHC-I protein investigated by long-term molecular dynamics simulations. BMC Immunol 2022; 23:36. [PMID: 35902791 PMCID: PMC9335952 DOI: 10.1186/s12865-022-00510-7] [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: 07/04/2022] [Accepted: 07/13/2022] [Indexed: 11/10/2022] Open
Abstract
Background Major histocompatibility complexes (MHCs) play a crucial role in the cell-mediated adaptive immune response as they present antigenic peptides (p) which are recognized by host T cells through a complex formation of the T cell receptor (TCR) with pMHC. In the present study, we report on changes in conformational flexibility within a pMHC molecule upon TCR binding by looking at molecular dynamics (MD) simulations of the free and the TCR-bound pMHC-I protein of the LC13-HLA-B*44:05-pEEYLQAFTY complex. Results We performed long-term MD simulations with a total simulation time of 8 µs, employing 10 independent 400 ns replicas for the free and the TCR-bound pMHC system. Upon TCR ligation, we observed a reduced dynamic flexibility in the central residues of the peptide and the MHC α1-helix, altered occurrences of hydrogen bonds between the peptide and the MHC, a reduced conformational entropy of the peptide-binding groove, as well as a decreased solvent accessible surface area. Conclusions In summary, our results from 8 µs MD simulations indicate a restricted conformational space of the MHC peptide-binding groove upon TCR ligation and suggest a minimum simulation time of approximately 100 ns for biomolecules of comparable complexity to draw meaningful conclusions. Given the relatively long total simulation time, our results contribute to a more detailed view on conformational flexibility properties of the investigated free and TCR-bound pMHC-I system. Supplementary Information The online version contains supplementary material available at 10.1186/s12865-022-00510-7.
Collapse
Affiliation(s)
- Lisa Tomasiak
- Institute of Biosimulation and Bioinformatics, Center for Medical Statistics, Informatics and Intelligent Systems, Medical University of Vienna, Spitalgasse 23, 1090, Vienna, Austria
| | - Rudolf Karch
- Institute of Biosimulation and Bioinformatics, Center for Medical Statistics, Informatics and Intelligent Systems, Medical University of Vienna, Spitalgasse 23, 1090, Vienna, Austria.
| | - Wolfgang Schreiner
- Institute of Biosimulation and Bioinformatics, Center for Medical Statistics, Informatics and Intelligent Systems, Medical University of Vienna, Spitalgasse 23, 1090, Vienna, Austria
| |
Collapse
|
7
|
Verkhivker GM. Conformational Dynamics and Mechanisms of Client Protein Integration into the Hsp90 Chaperone Controlled by Allosteric Interactions of Regulatory Switches: Perturbation-Based Network Approach for Mutational Profiling of the Hsp90 Binding and Allostery. J Phys Chem B 2022; 126:5421-5442. [PMID: 35853093 DOI: 10.1021/acs.jpcb.2c03464] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Understanding the allosteric mechanisms of the Hsp90 chaperone interactions with cochaperones and client protein clientele is fundamental to dissect activation and regulation of many proteins. In this work, atomistic simulations are combined with perturbation-based approaches and dynamic network modeling for a comparative mutational profiling of the Hsp90 binding and allosteric interaction networks in the three Hsp90 maturation complexes with FKBP51 and P23 cochaperones and the glucocorticoid receptor (GR) client. The conformational dynamics signatures of the Hsp90 complexes and dynamics fluctuation analysis revealed how the intrinsic plasticity of the Hsp90 dimer can be modulated by cochaperones and client proteins to stabilize the closed dimer state required at the maturation stage of the ATPase cycle. In silico deep mutational scanning of the protein residues characterized the hot spots of protein stability and binding affinity in the Hsp90 complexes, showing that binding hot spots may often coincide with the regulatory centers that modulate dynamic allostery in the Hsp90 dimer. We introduce a perturbation-based network approach for mutational scanning of allosteric residue potentials and characterize allosteric switch clusters that control mechanism of cochaperone-dependent client recognition and remodeling by the Hsp90 chaperone. The results revealed a conserved network of allosteric switches in the Hsp90 complexes that allow cochaperones and GR protein to become integrated into the Hsp90 system by anchoring to the conformational switch points in the functional Hsp90 regions. This study suggests that the Hsp90 binding and allostery may operate under a regulatory mechanism in which activation or repression of the Hsp90 activity can be pre-encoded in the allosterically regulated Hsp90 dimer motions. By binding directly to the conformational switch centers on the Hsp90, cochaperones and interacting proteins can efficiently modulate the allosteric interactions and long-range communications required for client remodeling and activation.
Collapse
Affiliation(s)
- Gennady M Verkhivker
- Keck Center for Science and Engineering, Schmid College of Science and Technology, Chapman University, 1 University Drive, Orange, California 92866, United States
- Depatment of Biomedical and Pharmaceutical Sciences, Chapman University School of Pharmacy, Irvine, California 92618, United States
| |
Collapse
|
8
|
In Search of a Dynamical Vocabulary: A Pipeline to Construct a Basis of Shared Traits in Large-Scale Motions of Proteins. APPLIED SCIENCES-BASEL 2022. [DOI: 10.3390/app12147157] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
The paradigmatic sequence–structure–dynamics–function relation in proteins is currently well established in the scientific community; in particular, a large effort has been made to probe the first connection, indeed providing convincing evidence of its strength and rationalizing it in a quantitative and general framework. In contrast, however, the role of dynamics as a link between structure and function has eluded a similarly clear-cut verification and description. In this work, we propose a pipeline aimed at building a basis for the quantitative characterization of the large-scale dynamics of a set of proteins, starting from the sole knowledge of their native structures. The method hinges on a dynamics-based clusterization, which allows a straightforward comparison with structural and functional protein classifications. The resulting basis set, obtained through the application to a group of related proteins, is shown to reproduce the salient large-scale dynamical features of the dataset. Most interestingly, the basis set is shown to encode the fluctuation patterns of homologous proteins not belonging to the initial dataset, thus highlighting the general applicability of the pipeline used to build it.
Collapse
|
9
|
Nguyen HTL, Huang DM. Systematic bottom-up molecular coarse-graining via force and torque matching using anisotropic particles. J Chem Phys 2022; 156:184118. [DOI: 10.1063/5.0085006] [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/14/2022] Open
Abstract
We derive a systematic and general method for parametrizing coarse-grained molecular models consisting of anisotropic particles from fine-grained (e.g. all-atom) models for condensed-phase molecular dynamics simulations. The method, which we call anisotropic force-matching coarse-graining (AFM-CG), is based on rigorous statistical mechanical principles, enforcing consistency between the coarse-grained and fine-grained phase-space distributions to derive equations for the coarse-grained forces, torques, masses, and moments of inertia in terms of properties of a condensed-phase fine-grained system. We verify the accuracy and efficiency of the method by coarse-graining liquid-state systems of two different anisotropic organic molecules, benzene and perylene, and show that the parametrized coarse-grained models more accurately describe properties of these systems than previous anisotropic coarse-grained models parametrized using other methods that do not account for finite-temperature and many-body effects on the condensed-phase coarse-grained interactions. The AFM-CG method will be useful for developing accurate and efficient dynamical simulation models of condensed-phase systems of molecules consisting of large, rigid, anisotropic fragments, such as liquid crystals, organic semiconductors, and nucleic acids.
Collapse
|
10
|
Verkhivker GM, Agajanian S, Kassab R, Krishnan K. Landscape-Based Protein Stability Analysis and Network Modeling of Multiple Conformational States of the SARS-CoV-2 Spike D614G Mutant: Conformational Plasticity and Frustration-Induced Allostery as Energetic Drivers of Highly Transmissible Spike Variants. J Chem Inf Model 2022; 62:1956-1978. [PMID: 35377633 DOI: 10.1021/acs.jcim.2c00124] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
The structural and functional studies of the SARS-CoV-2 spike protein variants revealed an important role of the D614G mutation that is shared across many variants of concern (VOCs), suggesting the effect of this mutation on the enhanced virus infectivity and transmissibility. The recent structural and biophysical studies provided important evidence about multiple conformational substates of the D614G spike protein. The development of a plausible mechanistic model that can explain the experimental observations from a more unified thermodynamic perspective is an important objective of the current work. In this study, we employed efficient and accurate coarse-grained simulations of multiple structural substates of the D614G spike trimers together with the ensemble-based mutational frustration analysis to characterize the dynamics signatures of the conformational landscapes. By combining the local frustration profiling of the conformational states with residue-based mutational scanning of protein stability and network analysis of allosteric interactions and communications, we determine the patterns of mutational sensitivity in the functional regions and sites of variants. We found that the D614G mutation may induce a considerable conformational adaptability of the open states in the SARS-CoV-2 spike protein without compromising the folding stability and integrity of the spike protein. The results suggest that the D614G mutant may employ a hinge-shift mechanism in which the dynamic couplings between the site of mutation and the interprotomer hinge modulate the interdomain interactions, global mobility change, and the increased stability of the open form. This study proposes that mutation-induced modulation of the conformational flexibility and energetic frustration at the interprotomer interfaces may serve as an efficient mechanism for allosteric regulation of the SARS-CoV-2 spike proteins.
Collapse
Affiliation(s)
- Gennady M Verkhivker
- Keck Center for Science and Engineering, Graduate Program in Computational and Data Sciences, Schmid College of Science and Technology, Chapman University, Orange, California 92866, United States.,Department of Biomedical and Pharmaceutical Sciences, Chapman University School of Pharmacy, Irvine, California 92618, United States
| | - Steve Agajanian
- Keck Center for Science and Engineering, Graduate Program in Computational and Data Sciences, Schmid College of Science and Technology, Chapman University, Orange, California 92866, United States
| | - Ryan Kassab
- Keck Center for Science and Engineering, Graduate Program in Computational and Data Sciences, Schmid College of Science and Technology, Chapman University, Orange, California 92866, United States
| | - Keerthi Krishnan
- Keck Center for Science and Engineering, Graduate Program in Computational and Data Sciences, Schmid College of Science and Technology, Chapman University, Orange, California 92866, United States
| |
Collapse
|
11
|
Allosteric Determinants of the SARS-CoV-2 Spike Protein Binding with Nanobodies: Examining Mechanisms of Mutational Escape and Sensitivity of the Omicron Variant. Int J Mol Sci 2022; 23:ijms23042172. [PMID: 35216287 PMCID: PMC8877688 DOI: 10.3390/ijms23042172] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2021] [Revised: 02/13/2022] [Accepted: 02/14/2022] [Indexed: 02/04/2023] Open
Abstract
Structural and biochemical studies have recently revealed a range of rationally engineered nanobodies with efficient neutralizing capacity against the SARS-CoV-2 virus and resilience against mutational escape. In this study, we performed a comprehensive computational analysis of the SARS-CoV-2 spike trimer complexes with single nanobodies Nb6, VHH E, and complex with VHH E/VHH V nanobody combination. We combined coarse-grained and all-atom molecular simulations and collective dynamics analysis with binding free energy scanning, perturbation-response scanning, and network centrality analysis to examine mechanisms of nanobody-induced allosteric modulation and cooperativity in the SARS-CoV-2 spike trimer complexes with these nanobodies. By quantifying energetic and allosteric determinants of the SARS-CoV-2 spike protein binding with nanobodies, we also examined nanobody-induced modulation of escaping mutations and the effect of the Omicron variant on nanobody binding. The mutational scanning analysis supported the notion that E484A mutation can have a significant detrimental effect on nanobody binding and result in Omicron-induced escape from nanobody neutralization. Our findings showed that SARS-CoV-2 spike protein might exploit the plasticity of specific allosteric hotspots to generate escape mutants that alter response to binding without compromising activity. The network analysis supported these findings showing that VHH E/VHH V nanobody binding can induce long-range couplings between the cryptic binding epitope and ACE2-binding site through a broader ensemble of communication paths that is less dependent on specific mediating centers and therefore may be less sensitive to mutational perturbations of functional residues. The results suggest that binding affinity and long-range communications of the SARS-CoV-2 complexes with nanobodies can be determined by structurally stable regulatory centers and conformationally adaptable hotspots that are allosterically coupled and collectively control resilience to mutational escape.
Collapse
|
12
|
Exploring Mechanisms of Allosteric Regulation and Communication Switching in the Multiprotein Regulatory Complexes of the Hsp90 Chaperone with Cochaperones and Client Proteins : Atomistic Insights from Integrative Biophysical Modeling and Network Analysis of Conformational Landscapes. J Mol Biol 2022; 434:167506. [DOI: 10.1016/j.jmb.2022.167506] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2022] [Revised: 02/14/2022] [Accepted: 02/15/2022] [Indexed: 12/16/2022]
|
13
|
Verkhivker G. Conformational Flexibility and Local Frustration in the Functional States of the SARS-CoV-2 Spike B.1.1.7 and B.1.351 Variants: Mutation-Induced Allosteric Modulation Mechanism of Functional Dynamics and Protein Stability. Int J Mol Sci 2022; 23:ijms23031646. [PMID: 35163572 PMCID: PMC8836237 DOI: 10.3390/ijms23031646] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/16/2021] [Revised: 01/22/2022] [Accepted: 01/29/2022] [Indexed: 02/01/2023] Open
Abstract
Structural and functional studies of the SARS-CoV-2 spike proteins have recently determined distinct functional states of the B.1.1.7 and B.1.351 spike variants, providing a molecular framework for understanding the mechanisms that link the effect of mutations with the enhanced virus infectivity and transmissibility. A detailed dynamic and energetic analysis of these variants was undertaken in the present work to quantify the effects of different mutations on functional conformational changes and stability of the SARS-CoV-2 spike protein. We employed the efficient and accurate coarse-grained (CG) simulations of multiple functional states of the D614G mutant, B.1.1.7 and B.1.351 spike variants to characterize conformational dynamics of the SARS-CoV-2 spike proteins and identify dynamic signatures of the functional regions that regulate transitions between the closed and open forms. By combining molecular simulations with full atomistic reconstruction of the trajectories and the ensemble-based mutational frustration analysis, we characterized how the intrinsic flexibility of specific spike regions can control functional conformational changes required for binding with the host-cell receptor. Using the residue-based mutational scanning of protein stability, we determined protein stability hotspots and identified potential energetic drivers favoring the receptor-accessible open spike states for the B.1.1.7 and B.1.351 spike variants. The results suggested that modulation of the energetic frustration at the inter-protomer interfaces can serve as a mechanism for allosteric couplings between mutational sites and the inter-protomer hinges of functional motions. The proposed mechanism of mutation-induced energetic frustration may result in greater adaptability and the emergence of multiple conformational states in the open form. This study suggested that SARS-CoV-2 B.1.1.7 and B.1.351 variants may leverage the intrinsic plasticity of functional regions in the spike protein for mutation-induced modulation of protein dynamics and allosteric regulation to control binding with the host cell receptor.
Collapse
Affiliation(s)
- Gennady Verkhivker
- Keck Center for Science and Engineering, Graduate Program in Computational and Data Sciences, Schmid College of Science and Technology, Chapman University, Orange, CA 92866, USA; ; Tel.: +17-14-516-4586
- Department of Biomedical and Pharmaceutical Sciences, Chapman University School of Pharmacy, Irvine, CA 92618, USA
| |
Collapse
|
14
|
Kuriata A, Badaczewska-Dawid AE, Pujols J, Ventura S, Kmiecik S. Protocols for Rational Design of Protein Solubility and Aggregation Properties Using Aggrescan3D Standalone. Methods Mol Biol 2022; 2340:17-40. [PMID: 35167068 DOI: 10.1007/978-1-0716-1546-1_2] [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: 06/14/2023]
Abstract
Protein aggregation is a major hurdle in the development and manufacturing of protein-based therapeutics. Development of aggregation-resistant and stable protein variants can be guided by rational redesign using computational tools. Here, we describe the architecture and functionalities of the Aggrescan3D (A3D) standalone package for the rational design of protein solubility and aggregation properties based on three-dimensional protein structures. We present the case studies of the three therapeutic proteins, including antibodies, exploring the practical use of the A3D standalone tool. The case studies demonstrate that protein solubility can be easily improved by the A3D prediction of non-destabilizing amino acid mutations at the protein surfaces.
Collapse
Affiliation(s)
- Aleksander Kuriata
- Faculty of Chemistry, Biological and Chemical Research Centre, University of Warsaw, Warsaw, Poland
| | | | - Jordi Pujols
- Institut de Biotecnologia i de Biomedicina (IBB), Universitat Autònoma de Barcelona, Bellaterra, Spain
- Departament de Bioquímica i Biologia Molecular, Universitat Autònoma de Barcelona, Bellaterra, Spain
| | - Salvador Ventura
- Institut de Biotecnologia i de Biomedicina (IBB), Universitat Autònoma de Barcelona, Bellaterra, Spain
- Departament de Bioquímica i Biologia Molecular, Universitat Autònoma de Barcelona, Bellaterra, Spain
| | - Sebastian Kmiecik
- Biological and Chemical Research Centre, University of Warsaw, Warsaw, Poland.
| |
Collapse
|
15
|
Fritz ZR, Schloss RS, Yarmush ML, Williams LJ. HSymM-guided engineering of the immunodominant p53 transactivation domain putative peptide antigen for improved binding to its anti-p53 monoclonal antibody. Bioorg Med Chem Lett 2021; 51:128341. [PMID: 34454062 PMCID: PMC8526406 DOI: 10.1016/j.bmcl.2021.128341] [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: 06/01/2021] [Revised: 08/17/2021] [Accepted: 08/22/2021] [Indexed: 02/07/2023]
Abstract
A novel engineering strategy to improve autoantibody detection with peptide fragments derived from the parent antigen is presented. The model system studied was the binding of the putative p53 TAD peptide antigen (residues 46-55) to its cognate anti-p53 antibody, ab28. Each engineered peptide contained the full decapeptide epitope and differed only in the flanking regions. Since minimal structural information was available to guide the design, a simple epitope:paratope binding model was applied. The Hidden Symmetry Model, which we recently reported, was used to guide peptide design and estimate per-residue contributions to interaction free energy as a function of added C- and N-terminal flanking peptides. Twenty-four peptide constructs were designed, synthesized, and assessed for binding affinity to ab28 by surface plasmon resonance, and a subset of these peptides were evaluated in a simulated immunoassay for limit of detection. Many peptides exhibited over 200-fold enhancements in binding affinity and improved limits of detection. The epitope was reevaluated and is proposed to be the undecapeptide corresponding to residues 45-55. HSymM calculated binding free energy and experimental data were found to be in good agreement (R2 > 0.75).
Collapse
Affiliation(s)
- Zachary R Fritz
- Department of Biomedical Engineering, Rutgers University, Piscataway, NJ 08854, United States
| | - Rene S Schloss
- Department of Biomedical Engineering, Rutgers University, Piscataway, NJ 08854, United States
| | - Martin L Yarmush
- Department of Biomedical Engineering, Rutgers University, Piscataway, NJ 08854, United States
| | - Lawrence J Williams
- Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, United States.
| |
Collapse
|
16
|
Liwo A, Czaplewski C, Sieradzan AK, Lipska AG, Samsonov SA, Murarka RK. Theory and Practice of Coarse-Grained Molecular Dynamics of Biologically Important Systems. Biomolecules 2021; 11:1347. [PMID: 34572559 PMCID: PMC8465211 DOI: 10.3390/biom11091347] [Citation(s) in RCA: 25] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2021] [Revised: 09/03/2021] [Accepted: 09/09/2021] [Indexed: 12/16/2022] Open
Abstract
Molecular dynamics with coarse-grained models is nowadays extensively used to simulate biomolecular systems at large time and size scales, compared to those accessible to all-atom molecular dynamics. In this review article, we describe the physical basis of coarse-grained molecular dynamics, the coarse-grained force fields, the equations of motion and the respective numerical integration algorithms, and selected practical applications of coarse-grained molecular dynamics. We demonstrate that the motion of coarse-grained sites is governed by the potential of mean force and the friction and stochastic forces, resulting from integrating out the secondary degrees of freedom. Consequently, Langevin dynamics is a natural means of describing the motion of a system at the coarse-grained level and the potential of mean force is the physical basis of the coarse-grained force fields. Moreover, the choice of coarse-grained variables and the fact that coarse-grained sites often do not have spherical symmetry implies a non-diagonal inertia tensor. We describe selected coarse-grained models used in molecular dynamics simulations, including the most popular MARTINI model developed by Marrink's group and the UNICORN model of biological macromolecules developed in our laboratory. We conclude by discussing examples of the application of coarse-grained molecular dynamics to study biologically important processes.
Collapse
Affiliation(s)
- Adam Liwo
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland; (C.C.); (A.K.S.); (A.G.L.); (S.A.S.)
| | - Cezary Czaplewski
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland; (C.C.); (A.K.S.); (A.G.L.); (S.A.S.)
| | - Adam K. Sieradzan
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland; (C.C.); (A.K.S.); (A.G.L.); (S.A.S.)
| | - Agnieszka G. Lipska
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland; (C.C.); (A.K.S.); (A.G.L.); (S.A.S.)
| | - Sergey A. Samsonov
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland; (C.C.); (A.K.S.); (A.G.L.); (S.A.S.)
| | - Rajesh K. Murarka
- Department of Chemistry, Indian Institute of Science Education and Research Bhopal, Bhopal Bypass Road, Bhopal 462066, MP, India;
| |
Collapse
|
17
|
Kurcinski M, Kmiecik S, Zalewski M, Kolinski A. Protein-Protein Docking with Large-Scale Backbone Flexibility Using Coarse-Grained Monte-Carlo Simulations. Int J Mol Sci 2021; 22:ijms22147341. [PMID: 34298961 PMCID: PMC8306105 DOI: 10.3390/ijms22147341] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/19/2021] [Revised: 07/03/2021] [Accepted: 07/04/2021] [Indexed: 12/21/2022] Open
Abstract
Most of the protein–protein docking methods treat proteins as almost rigid objects. Only the side-chains flexibility is usually taken into account. The few approaches enabling docking with a flexible backbone typically work in two steps, in which the search for protein–protein orientations and structure flexibility are simulated separately. In this work, we propose a new straightforward approach for docking sampling. It consists of a single simulation step during which a protein undergoes large-scale backbone rearrangements, rotations, and translations. Simultaneously, the other protein exhibits small backbone fluctuations. Such extensive sampling was possible using the CABS coarse-grained protein model and Replica Exchange Monte Carlo dynamics at a reasonable computational cost. In our proof-of-concept simulations of 62 protein–protein complexes, we obtained acceptable quality models for a significant number of cases.
Collapse
|
18
|
Molecular Analysis of Streptomycin Resistance Genes in Clinical Strains of Mycobacterium tuberculosis and Biocomputational Analysis of the MtGidB L101F Variant. Antibiotics (Basel) 2021; 10:antibiotics10070807. [PMID: 34356728 PMCID: PMC8300841 DOI: 10.3390/antibiotics10070807] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2021] [Revised: 06/29/2021] [Accepted: 06/30/2021] [Indexed: 12/30/2022] Open
Abstract
Globally, tuberculosis (TB) remains a prevalent threat to public health. In 2019, TB affected 10 million people and caused 1.4 million deaths. The major challenge for controlling this infectious disease is the emergence and spread of drug-resistant Mycobacterium tuberculosis, the causative agent of TB. The antibiotic streptomycin is not a current first-line anti-TB drug. However, WHO recommends its use in patients infected with a streptomycin-sensitive strain. Several mutations in the M. tuberculosisrpsL, rrs and gidB genes have proved association with streptomycin resistance. In this study, we performed a molecular analysis of these genes in clinical isolates to determine the prevalence of known or novel mutations. Here, we describe the genetic analysis outcome. Furthermore, a biocomputational analysis of the MtGidB L101F variant, the product of a novel mutation detected in gidB during molecular analysis, is also reported as a theoretical approach to study the apparent genotype-phenotype association.
Collapse
|
19
|
Malik RM, Fazal S, Kamal MA. Computational Analysis of Domains Vulnerable to HPV-16 E6 Oncoprotein and Corresponding Hot Spot Residues. Protein Pept Lett 2021; 28:414-425. [PMID: 32703126 DOI: 10.2174/0929866527666200722134801] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/29/2020] [Revised: 06/19/2020] [Accepted: 06/28/2020] [Indexed: 02/08/2023]
Abstract
BACKGROUND Human Papilloma Virus (HPV) is the primary cause of cancers in cervix, head and neck regions. Oncoprotein E6 of HPV-16, after infecting human body, alters host protein- protein interaction networks. E6 interacts with several proteins, causing the infection to progress into cervical cancer. The molecular basis for these interactions is the presence of short linear peptide motifs on E6 identical to those on human proteins. METHODS Motifs of LXXLL and E/DLLL/V-G after identification on E6, were analyzed for their dynamic fluctuations by use of elastic network models. Correlation analysis of amino acid residues of E6 was also performed in specific regions of motifs. RESULTS Arginine, Leucine, Glutamine, Threonine and Glutamic acid have been identified as hot spot residues of E6 which can subsequently provide a platform for drug designing and understanding of pathogenesis of cervical cancer. These amino acids play a significant role in stabilizing interactions with host proteins, ultimately causing infections and cancers. CONCLUSION Our study validates the role of linear binding motifs of E6 of HPV in interacting with these proteins as an important event in the propagation of HPV in human cells and its transformation into cervical cancer. The study further predicts the domains of protein kinase and armadillo as part of the regions involved in the interaction of E6AP, Paxillin and TNF R1, with viral E6.
Collapse
Affiliation(s)
| | - Sahar Fazal
- Capital University of Science and Technology, Islamabad, Pakistan
| | - Mohammad Amjad Kamal
- West China School of Nursing / Institutes for Systems Genetics, Frontiers Science Center for Disease-related Molecular Network, West China Hospital, Sichuan University, Chengdu 610041, Sichuan, China
| |
Collapse
|
20
|
Verkhivker GM, Agajanian S, Oztas DY, Gupta G. Comparative Perturbation-Based Modeling of the SARS-CoV-2 Spike Protein Binding with Host Receptor and Neutralizing Antibodies: Structurally Adaptable Allosteric Communication Hotspots Define Spike Sites Targeted by Global Circulating Mutations. Biochemistry 2021; 60:1459-1484. [PMID: 33900725 PMCID: PMC8098775 DOI: 10.1021/acs.biochem.1c00139] [Citation(s) in RCA: 41] [Impact Index Per Article: 13.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/20/2021] [Revised: 04/12/2021] [Indexed: 12/11/2022]
Abstract
In this study, we used an integrative computational approach to examine molecular mechanisms and determine functional signatures underlying the role of functional residues in the SARS-CoV-2 spike protein that are targeted by novel mutational variants and antibody-escaping mutations. Atomistic simulations and functional dynamics analysis are combined with alanine scanning and mutational sensitivity profiling of the SARS-CoV-2 spike protein complexes with the ACE2 host receptor and the REGN-COV2 antibody cocktail(REG10987+REG10933). Using alanine scanning and mutational sensitivity analysis, we have shown that K417, E484, and N501 residues correspond to key interacting centers with a significant degree of structural and energetic plasticity that allow mutants in these positions to afford the improved binding affinity with ACE2. Through perturbation-based network modeling and community analysis of the SARS-CoV-2 spike protein complexes with ACE2, we demonstrate that E406, N439, K417, and N501 residues serve as effector centers of allosteric interactions and anchor major intermolecular communities that mediate long-range communication in the complexes. The results provide support to a model according to which mutational variants and antibody-escaping mutations constrained by the requirements for host receptor binding and preservation of stability may preferentially select structurally plastic and energetically adaptable allosteric centers to differentially modulate collective motions and allosteric interactions in the complexes with the ACE2 enzyme and REGN-COV2 antibody combination. This study suggests that the SARS-CoV-2 spike protein may function as a versatile and functionally adaptable allosteric machine that exploits the plasticity of allosteric regulatory centers to fine-tune response to antibody binding without compromising the activity of the spike protein.
Collapse
Affiliation(s)
- Gennady M. Verkhivker
- Keck Center for Science and Engineering, Schmid
College of Science and Technology, Chapman University, One
University Drive, Orange, California 92866, United States
- Depatment of Biomedical and Pharmaceutical Sciences,
Chapman University School of Pharmacy, Irvine, California
92618, United States
| | - Steve Agajanian
- Keck Center for Science and Engineering, Schmid
College of Science and Technology, Chapman University, One
University Drive, Orange, California 92866, United States
| | - Deniz Yazar Oztas
- Keck Center for Science and Engineering, Schmid
College of Science and Technology, Chapman University, One
University Drive, Orange, California 92866, United States
| | - Grace Gupta
- Keck Center for Science and Engineering, Schmid
College of Science and Technology, Chapman University, One
University Drive, Orange, California 92866, United States
| |
Collapse
|
21
|
Verkhivker GM, Di Paola L. Integrated Biophysical Modeling of the SARS-CoV-2 Spike Protein Binding and Allosteric Interactions with Antibodies. J Phys Chem B 2021; 125:4596-4619. [PMID: 33929853 PMCID: PMC8098774 DOI: 10.1021/acs.jpcb.1c00395] [Citation(s) in RCA: 39] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2021] [Revised: 04/16/2021] [Indexed: 02/07/2023]
Abstract
Structural and biochemical studies of the severe acute respiratory syndrome (SARS)-CoV-2 spike glycoproteins and complexes with highly potent antibodies have revealed multiple conformation-dependent epitopes highlighting conformational plasticity of spike proteins and capacity for eliciting specific binding and broad neutralization responses. In this study, we used coevolutionary analysis, molecular simulations, and perturbation-based hierarchical network modeling of the SARS-CoV-2 spike protein complexes with a panel of antibodies targeting distinct epitopes to explore molecular mechanisms underlying binding-induced modulation of dynamics and allosteric signaling in the spike proteins. Through coevolutionary analysis of the SARS-CoV-2 spike proteins, we identified highly coevolving hotspots and functional clusters that enable a functional cross-talk between distant allosteric regions in the SARS-CoV-2 spike complexes with antibodies. Coarse-grained and all-atom molecular dynamics simulations combined with mutational sensitivity mapping and perturbation-based profiling of the SARS-CoV-2 receptor-binding domain (RBD) complexes with CR3022 and CB6 antibodies enabled a detailed validation of the proposed approach and an extensive quantitative comparison with the experimental structural and deep mutagenesis scanning data. By combining in silico mutational scanning, perturbation-based modeling, and network analysis of the SARS-CoV-2 spike trimer complexes with H014, S309, S2M11, and S2E12 antibodies, we demonstrated that antibodies can incur specific and functionally relevant changes by modulating allosteric propensities and collective dynamics of the SARS-CoV-2 spike proteins. The results provide a novel insight into regulatory mechanisms of SARS-CoV-2 S proteins showing that antibody-escaping mutations can preferentially target structurally adaptable energy hotspots and allosteric effector centers that control functional movements and allosteric communication in the complexes.
Collapse
Affiliation(s)
- Gennady M. Verkhivker
- Keck Center for Science and Engineering, Schmid
College of Science and Technology, Chapman University, One
University Drive, Orange, California 92866, United States
- Department of Biomedical and Pharmaceutical Sciences,
Chapman University School of Pharmacy, Irvine, California
92618, United States
| | - Luisa Di Paola
- Unit of Chemical-Physics Fundamentals in Chemical
Engineering, Department of Engineering, Università Campus Bio-Medico
di Roma, via Álvaro del Portillo 21, 00128 Rome,
Italy
| |
Collapse
|
22
|
Malik R, Fazal S. Insights into the Dynamic Fluctuations of the Protein HPV16 E1 and Identification of Motifs by Using Elastic Network Modeling. Protein Pept Lett 2021; 28:1061-1070. [PMID: 33858307 DOI: 10.2174/0929866528666210415114858] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/26/2020] [Revised: 02/14/2021] [Accepted: 02/18/2021] [Indexed: 11/22/2022]
Abstract
BACKGROUND Cancers of cervix, head and neck regions have been found to be associated with Human Papilloma Virus (HPV) infection. E1 protein makes an important papillomavirus replication factor. Among the ORFs of papillomaviruses, the most conserved sequence is that of the E1 ORF. It is the viral helicase with being a member of class of ATPases associated with diverse cellular activities (AAA+) helicases. The interactions of E1 with human DNA and proteins occurs in the presence of short linear peptide motifs on E1 identical to those on human proteins. METHODS Different Motifs were identified on HPV16 E1 by using ELMs. Elastic network models were generated by using 3D structures of E1. Their dynamic fluctuations were analyzed on the basis of B factors, correlation analysis and deformation energies. RESULTS 3 motifs were identified on E1 which can interact with Cdk and Cyclin domains of human proteins. 11 motifs identified on E1 have their CDs of Pkinase on human proteins. LIG_MYND_2 has been identified as involved in stabilizing interaction of E1 with Hsp40 and Hsp70. These motifs and amino acids comprising these motifs play a major role in maintaining interactions with human proteins, ultimately causing infections leading to cancers. CONCLUSION Our study identified various motifs on E1 which interact with specific counter domains found in human proteins, already reported having the interactions with E1. We also validated the involvement of these specific motifs containing regions of E1 by modeling elastic networks of E1. These motif involving interactions could be used as drug targets.
Collapse
Affiliation(s)
- Rabbiah Malik
- Capital University of Science and Technology, Islamabad. Pakistan
| | - Sahar Fazal
- Capital University of Science and Technology, Islamabad. Pakistan
| |
Collapse
|
23
|
Badaczewska-Dawid AE, Kolinski A, Kmiecik S. Protocols for Fast Simulations of Protein Structure Flexibility Using CABS-Flex and SURPASS. Methods Mol Biol 2021; 2165:337-353. [PMID: 32621235 DOI: 10.1007/978-1-0716-0708-4_20] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
Conformational flexibility of protein structures can play an important role in protein function. The flexibility is often studied using computational methods since experimental characterization can be difficult. Depending on protein system size, computational tools may require large computational resources or significant simplifications in the modeled systems to speed up calculations. In this work, we present the protocols for efficient simulations of flexibility of folded protein structures that use coarse-grained simulation tools of different resolutions: medium, represented by CABS-flex, and low, represented by SUPRASS. We test the protocols using a set of 140 globular proteins and compare the results with structure fluctuations observed in MD simulations, ENM modeling, and NMR ensembles. As demonstrated, CABS-flex predictions show high correlation to experimental and MD simulation data, while SURPASS is less accurate but promising in terms of future developments.
Collapse
Affiliation(s)
- Aleksandra E Badaczewska-Dawid
- Faculty of Chemistry, Biological and Chemical Research Center, University of Warsaw, Warsaw, Poland.,Department of Chemistry, Iowa State University, Ames, IA, USA
| | - Andrzej Kolinski
- Faculty of Chemistry, Biological and Chemical Research Center, University of Warsaw, Warsaw, Poland
| | - Sebastian Kmiecik
- Faculty of Chemistry, Biological and Chemical Research Center, University of Warsaw, Warsaw, Poland.
| |
Collapse
|
24
|
Verkhivker GM, Di Paola L. Dynamic Network Modeling of Allosteric Interactions and Communication Pathways in the SARS-CoV-2 Spike Trimer Mutants: Differential Modulation of Conformational Landscapes and Signal Transmission via Cascades of Regulatory Switches. J Phys Chem B 2021; 125:850-873. [PMID: 33448856 PMCID: PMC7839160 DOI: 10.1021/acs.jpcb.0c10637] [Citation(s) in RCA: 44] [Impact Index Per Article: 14.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/26/2020] [Revised: 01/08/2021] [Indexed: 12/13/2022]
Abstract
The rapidly growing body of structural and biochemical studies of the SARS-CoV-2 spike glycoprotein has revealed a variety of distinct functional states with radically different arrangements of the receptor-binding domain, highlighting a remarkable function-driven conformational plasticity and adaptability of the spike proteins. In this study, we examined molecular mechanisms underlying conformational and dynamic changes in the SARS-CoV-2 spike mutant trimers through the lens of dynamic analysis of allosteric interaction networks and atomistic modeling of signal transmission. Using an integrated approach that combined coarse-grained molecular simulations, protein stability analysis, and perturbation-based modeling of residue interaction networks, we examined how mutations in the regulatory regions of the SARS-CoV-2 spike protein can differentially affect dynamics and allosteric signaling in distinct functional states. The results of this study revealed key functional regions and regulatory centers that govern collective dynamics, allosteric interactions, and control signal transmission in the SARS-CoV-2 spike proteins. We found that the experimentally confirmed regulatory hotspots that dictate dynamic switching between conformational states of the SARS-CoV-2 spike protein correspond to the key hinge sites and global mediating centers of the allosteric interaction networks. The results of this study provide a novel insight into allosteric regulatory mechanisms of SARS-CoV-2 spike proteins showing that mutations at the key regulatory positions can differentially modulate distribution of states and determine topography of signal communication pathways operating through state-specific cascades of control switch points. This analysis provides a plausible strategy for allosteric probing of the conformational equilibrium and therapeutic intervention by targeting specific hotspots of allosteric interactions and communications in the SARS-CoV-2 spike proteins.
Collapse
Affiliation(s)
- Gennady M. Verkhivker
- Keck
Center for Science and Engineering, Schmid College of Science and
Technology, Chapman University, One University Drive, Orange, California 92866, United States
- Department
of Biomedical and Pharmaceutical Sciences, Chapman University School of Pharmacy, Irvine, California 92618, United States
| | - Luisa Di Paola
- Unit
of Chemical-Physics Fundamentals in Chemical Engineering, Department
of Engineering, Università Campus
Bio-Medico di Roma, via
Álvaro del Portillo 21, 00128 Rome, Italy
| |
Collapse
|
25
|
Gong W, Liu Y, Zhao Y, Wang S, Han Z, Li C. Equally Weighted Multiscale Elastic Network Model and Its Comparison with Traditional and Parameter-Free Models. J Chem Inf Model 2021; 61:921-937. [PMID: 33496590 DOI: 10.1021/acs.jcim.0c01178] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Dynamical properties of proteins play an essential role in their function exertion. The elastic network model (ENM) is an effective and efficient tool in characterizing the intrinsic dynamical properties encoded in biomacromolecule structures. The Gaussian network model (GNM) and anisotropic network model (ANM) are the two often-used ENM models. Here, we introduce an equally weighted multiscale ENM (equally weighted mENM) based on the original mENM (denoted as mENM), in which fitting weights of Kirchhoff/Hessian matrixes in mENM are removed since they neglect the details of pairwise interactions. Then, we perform its comparison with the mENM, traditional ENM, and parameter-free ENM (pfENM) in reproducing dynamical properties for the six representative proteins whose molecular dynamics (MD) trajectories are available in http://mmb.pcb.ub.es/MoDEL/. In the results, for B-factor prediction, mENM performs best, while the equally weighted mENM performs also well, better than the traditional ENM and pfENM models. As to the dynamical cross-correlation map calculation, mENM performs worst, while the results produced from the equally weighted mENM and pfENM models are close to those from MD trajectories with the latter a little better than the former. Furthermore, encouragingly, the equally weighted mANM displays the best performance in capturing the functional motional modes, followed by pfANM and traditional ANM models, while the mANM fails in all the cases. This work is helpful for strengthening the understanding of the elastic network model and provides a valuable guide for researchers to utilize the model to explore protein dynamics.
Collapse
Affiliation(s)
- Weikang Gong
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing 100124, China.,Beijing International Science and Technology Cooperation Base for Intelligent Physiological Measurement and Clinical Transformation, Beijing University of Technology, Beijing 100124, China
| | - Yang Liu
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing 100124, China.,Beijing International Science and Technology Cooperation Base for Intelligent Physiological Measurement and Clinical Transformation, Beijing University of Technology, Beijing 100124, China
| | - Yanpeng Zhao
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing 100124, China.,Beijing International Science and Technology Cooperation Base for Intelligent Physiological Measurement and Clinical Transformation, Beijing University of Technology, Beijing 100124, China
| | - Shihao Wang
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing 100124, China.,Beijing International Science and Technology Cooperation Base for Intelligent Physiological Measurement and Clinical Transformation, Beijing University of Technology, Beijing 100124, China
| | - Zhongjie Han
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing 100124, China.,Beijing International Science and Technology Cooperation Base for Intelligent Physiological Measurement and Clinical Transformation, Beijing University of Technology, Beijing 100124, China
| | - Chunhua Li
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing 100124, China.,Beijing International Science and Technology Cooperation Base for Intelligent Physiological Measurement and Clinical Transformation, Beijing University of Technology, Beijing 100124, China
| |
Collapse
|
26
|
Molecular Cloning, Purification and Characterization of Mce1R of Mycobacterium tuberculosis. Mol Biotechnol 2021; 63:200-220. [PMID: 33423211 DOI: 10.1007/s12033-020-00293-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 12/18/2020] [Indexed: 10/22/2022]
Abstract
The mce1 operon of Mycobacterium tuberculosis, important for lipid metabolism/transport, host cell invasion, modulation of host immune response and pathogenicity, is under the transcriptional control of Mce1R. Hence characterizing Mce1R is an important step for novel anti-tuberculosis drug discovery. The present study reports functional and in silico characterization of Mce1R. In this work, we have computationally modeled the structure of Mce1R and have validated the structure by computational and experimental methods. Mce1R has been shown to harbor the canonical VanR-like structure with a flexible N-terminal 'arm', carrying conserved positively charged residues, most likely involved in the operator DNA binding. The mce1R gene has been cloned, expressed, purified and its DNA-binding activity has been measured in vitro. The Kd value for Mce1R-operator DNA interaction has been determined to be 0.35 ± 0.02 µM which implies that Mce1R binds to DNA with moderate affinity compared to the other FCD family of regulators. So far, this is the first report for measuring the DNA-binding affinity of any VanR-type protein. Despite significant sequence similarity at the N-terminal domain, the wHTH motif of Mce1R exhibits poor conservancy of amino acid residues, critical for DNA-binding, thus results in moderate DNA-binding affinity. The N-terminal DNA-binding domain is structurally dynamic while the C-terminal domain showed significant stability and such profile of structural dynamics is most likely to be preserved in the structural orthologs of Mce1R. In addition to this, a cavity has been detected in the C-terminal domain of Mce1R which contains a few conserved residues. Comparison with other FCD family of regulators suggests that most of the conserved residues might be critical for binding to specific ligand. The max pKd value and drug score for the cavity are estimated to be 9.04 and 109 respectively suggesting that the cavity represents a suitable target site for novel anti-tuberculosis drug discovery approaches.
Collapse
|
27
|
Verkhivker GM. Molecular Simulations and Network Modeling Reveal an Allosteric Signaling in the SARS-CoV-2 Spike Proteins. J Proteome Res 2020; 19:4587-4608. [PMID: 33006900 PMCID: PMC7640983 DOI: 10.1021/acs.jproteome.0c00654] [Citation(s) in RCA: 46] [Impact Index Per Article: 11.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2020] [Indexed: 12/13/2022]
Abstract
The development of computational strategies for the quantitative characterization of the functional mechanisms of SARS-CoV-2 spike proteins is of paramount importance in efforts to accelerate the discovery of novel therapeutic agents and vaccines combating the COVID-19 pandemic. Structural and biophysical studies have recently characterized the conformational landscapes of the SARS-CoV-2 spike glycoproteins in the prefusion form, revealing a spectrum of stable and more dynamic states. By employing molecular simulations and network modeling approaches, this study systematically examined functional dynamics and identified the regulatory centers of allosteric interactions for distinct functional states of the wild-type and mutant variants of the SARS-CoV-2 prefusion spike trimer. This study presents evidence that the SARS-CoV-2 spike protein can function as an allosteric regulatory engine that fluctuates between dynamically distinct functional states. Perturbation-based modeling of the interaction networks revealed a key role of the cross-talk between the effector hotspots in the receptor binding domain and the fusion peptide proximal region of the SARS-CoV-2 spike protein. The results have shown that the allosteric hotspots of the interaction networks in the SARS-CoV-2 spike protein can control the dynamic switching between functional conformational states that are associated with virus entry to the host receptor. This study offers a useful and novel perspective on the underlying mechanisms of the SARS-CoV-2 spike protein through the lens of allosteric signaling as a regulatory apparatus of virus transmission that could open up opportunities for targeted allosteric drug discovery against SARS-CoV-2 proteins and contribute to the rapid response to the current and potential future pandemic scenarios.
Collapse
Affiliation(s)
- Gennady M. Verkhivker
- Graduate
Program in Computational and Data Sciences, Keck Center for Science
and Engineering, Schmid College of Science and Technology, Chapman University, One University Drive, Orange, California 92866, United States
- Department
of Biomedical and Pharmaceutical Sciences, Chapman University School of Pharmacy, Irvine, California 92618, United States
| |
Collapse
|
28
|
Verkhivker G. Coevolution, Dynamics and Allostery Conspire in Shaping Cooperative Binding and Signal Transmission of the SARS-CoV-2 Spike Protein with Human Angiotensin-Converting Enzyme 2. Int J Mol Sci 2020; 21:ijms21218268. [PMID: 33158276 PMCID: PMC7672574 DOI: 10.3390/ijms21218268] [Citation(s) in RCA: 25] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/18/2020] [Revised: 11/02/2020] [Accepted: 11/03/2020] [Indexed: 02/07/2023] Open
Abstract
Binding to the host receptor is a critical initial step for the coronavirus SARS-CoV-2 spike protein to enter into target cells and trigger virus transmission. A detailed dynamic and energetic view of the binding mechanisms underlying virus entry is not fully understood and the consensus around the molecular origins behind binding preferences of SARS-CoV-2 for binding with the angiotensin-converting enzyme 2 (ACE2) host receptor is yet to be established. In this work, we performed a comprehensive computational investigation in which sequence analysis and modeling of coevolutionary networks are combined with atomistic molecular simulations and comparative binding free energy analysis of the SARS-CoV and SARS-CoV-2 spike protein receptor binding domains with the ACE2 host receptor. Different from other computational studies, we systematically examine the molecular and energetic determinants of the binding mechanisms between SARS-CoV-2 and ACE2 proteins through the lens of coevolution, conformational dynamics, and allosteric interactions that conspire to drive binding interactions and signal transmission. Conformational dynamics analysis revealed the important differences in mobility of the binding interfaces for the SARS-CoV-2 spike protein that are not confined to several binding hotspots, but instead are broadly distributed across many interface residues. Through coevolutionary network analysis and dynamics-based alanine scanning, we established linkages between the binding energy hotspots and potential regulators and carriers of signal communication in the virus-host receptor complexes. The results of this study detailed a binding mechanism in which the energetics of the SARS-CoV-2 association with ACE2 may be determined by cumulative changes of a number of residues distributed across the entire binding interface. The central findings of this study are consistent with structural and biochemical data and highlight drug discovery challenges of inhibiting large and adaptive protein-protein interfaces responsible for virus entry and infection transmission.
Collapse
Affiliation(s)
- Gennady Verkhivker
- Graduate Program in Computational and Data Sciences, Schmid College of Science and Technology, Chapman University, Orange, CA 92866, USA; ; Tel.: +1-714-516-4586
- Department of Biomedical and Pharmaceutical Sciences, Chapman University School of Pharmacy, Irvine, CA 92618, USA
| |
Collapse
|
29
|
Amamuddy OS, Verkhivker GM, Bishop ÖT. Impact of Early Pandemic Stage Mutations on Molecular Dynamics of SARS-CoV-2 M pro. J Chem Inf Model 2020; 60:5080-5102. [PMID: 32853525 PMCID: PMC7496595 DOI: 10.1021/acs.jcim.0c00634] [Citation(s) in RCA: 45] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/05/2020] [Indexed: 12/15/2022]
Abstract
A new coronavirus (SARS-CoV-2) is a global threat to world health and economy. Its dimeric main protease (Mpro), which is required for the proteolytic cleavage of viral precursor proteins, is a good candidate for drug development owing to its conservation and the absence of a human homolog. Improving our understanding of Mpro behavior can accelerate the discovery of effective therapies to reduce mortality. All-atom molecular dynamics (MD) simulations (100 ns) of 50 mutant Mpro dimers obtained from filtered sequences from the GISAID database were analyzed using root-mean-square deviation, root-mean-square fluctuation, Rg, averaged betweenness centrality, and geometry calculations. The results showed that SARS-CoV-2 Mpro essentially behaves in a similar manner to its SAR-CoV homolog. However, we report the following new findings from the variants: (1) Residues GLY15, VAL157, and PRO184 have mutated more than once in SARS CoV-2; (2) the D48E variant has lead to a novel "TSEEMLN"" loop at the binding pocket; (3) inactive apo Mpro does not show signs of dissociation in 100 ns MD; (4) a non-canonical pose for PHE140 widens the substrate binding surface; (5) dual allosteric pockets coinciding with various stabilizing and functional components of the substrate binding pocket were found to display correlated compaction dynamics; (6) high betweenness centrality values for residues 17 and 128 in all Mpro samples suggest their high importance in dimer stability-one such consequence has been observed for the M17I mutation whereby one of the N-fingers was highly unstable. (7) Independent coarse-grained Monte Carlo simulations suggest a relationship between the rigidity/mutability and enzymatic function. Our entire approach combining database preparation, variant retrieval, homology modeling, dynamic residue network (DRN), relevant conformation retrieval from 1-D kernel density estimates from reaction coordinates to other existing approaches of structural analysis, and data visualization within the coronaviral Mpro is also novel and is applicable to other coronaviral proteins.
Collapse
Affiliation(s)
- Olivier Sheik Amamuddy
- Research Unit in Bioinformatics, Department of Microbiology and Biochemistry, Rhodes University, Grahamstown 6140, South Africa
| | - Gennady M. Verkhivker
- Graduate Program in Computational and Data Sciences, Schmid College of Science and Technology, Chapman University, Orange, CA 92866, USA
- Department of Biomedical and Pharmaceutical Sciences, Chapman University School of Pharmacy, Irvine, CA 92618, USA
- Department of Pharmacology, Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California San Diego, 9500 Gilman Drive, La Jolla, CA 92093, USA
| | - Özlem Tastan Bishop
- Research Unit in Bioinformatics, Department of Microbiology and Biochemistry, Rhodes University, Grahamstown 6140, South Africa
| |
Collapse
|
30
|
Wang S, Gong W, Deng X, Liu Y, Li C. Exploring the dynamics of RNA molecules with multiscale Gaussian network model. Chem Phys 2020. [DOI: 10.1016/j.chemphys.2020.110820] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
|
31
|
Saldaño TE, Freixas VM, Tosatto SCE, Parisi G, Fernandez-Alberti S. Exploring Conformational Space with Thermal Fluctuations Obtained by Normal-Mode Analysis. J Chem Inf Model 2020; 60:3068-3080. [PMID: 32216314 DOI: 10.1021/acs.jcim.9b01136] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
Proteins in their native states can be represented as ensembles of conformers in dynamical equilibrium. Thermal fluctuations are responsible for transitions between these conformers. Normal-modes analysis (NMA) using elastic network models (ENMs) provides an efficient procedure to explore global dynamics of proteins commonly associated with conformational transitions. In the present work, we present an iterative approach to explore protein conformational spaces by introducing structural distortions according to their equilibrium dynamics at room temperature. The approach can be used either to perform unbiased explorations of conformational space or to explore guided pathways connecting two different conformations, e.g., apo and holo forms. In order to test its performance, four proteins with different magnitudes of structural distortions upon ligand binding have been tested. In all cases, the conformational selection model has been confirmed and the conformational space between apo and holo forms has been encompassed. Different strategies have been tested that impact on the efficiency either to achieve a desired conformational change or to achieve a balanced exploration of the protein conformational multiplicity.
Collapse
Affiliation(s)
- Tadeo E Saldaño
- Universidad Nacional de Quilmes/CONICET, Roque Saenz Peña 352, B1876BXD Bernal, Argentina
| | - Victor M Freixas
- Universidad Nacional de Quilmes/CONICET, Roque Saenz Peña 352, B1876BXD Bernal, Argentina
| | - Silvio C E Tosatto
- Department of Biomedical Sciences, University of Padova, Viale G. Colombo 3, 5131 Padova, Italy
| | - Gustavo Parisi
- Universidad Nacional de Quilmes/CONICET, Roque Saenz Peña 352, B1876BXD Bernal, Argentina
| | | |
Collapse
|
32
|
Stetz G, Astl L, Verkhivker GM. Exploring Mechanisms of Communication Switching in the Hsp90-Cdc37 Regulatory Complexes with Client Kinases through Allosteric Coupling of Phosphorylation Sites: Perturbation-Based Modeling and Hierarchical Community Analysis of Residue Interaction Networks. J Chem Theory Comput 2020; 16:4706-4725. [PMID: 32492340 DOI: 10.1021/acs.jctc.0c00280] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Understanding molecular principles underlying chaperone-based modulation of kinase client activity is critically important to dissect functions and activation mechanisms of many oncogenic proteins. The recent experimental studies have suggested that phosphorylation sites in the Hsp90 and Cdc37 proteins can serve as conformational communication switches of chaperone regulation and kinase interactions. However, a mechanism of allosteric coupling between phosphorylation sites in the Hsp90 and Cdc37 during client binding is poorly understood, and the molecular signatures underpinning specific roles of phosphorylation sites in the Hsp90 regulation remain unknown. In this work, we employed a combination of evolutionary analysis, coarse-grained molecular simulations together with perturbation-based network modeling and scanning of the unbound and bound Hsp90 and Cdc37 structures to quantify allosteric effects of phosphorylation sites and identify unique signatures that are characteristic for communication switches of kinase-specific client binding. By using network-based metrics of the dynamic intercommunity bridgeness and community centrality, we characterize specific signatures of phosphorylation switches involved in allosteric regulation. Through perturbation-based analysis of the dynamic residue interaction networks, we show that mutations of kinase-specific phosphorylation switches can induce long-range effects and lead to a global rewiring of the allosteric network and signal transmission in the Hsp90-Cdc37-kinase complex. We determine a specific group of phosphorylation sites in the Hsp90 where mutations may have a strong detrimental effect on allosteric interaction network, providing insight into the mechanism of phosphorylation-induced communication switching. The results demonstrate that kinase-specific phosphorylation switches of communications in the Hsp90 may be partly predisposed for their regulatory role based on preexisting allosteric propensities.
Collapse
Affiliation(s)
- Gabrielle Stetz
- Graduate Program in Computational and Data Sciences, Keck Center for Science and Engineering, Schmid College of Science and Technology, Chapman University, One University Drive, Orange, California 92866, United States
| | - Lindy Astl
- Graduate Program in Computational and Data Sciences, Keck Center for Science and Engineering, Schmid College of Science and Technology, Chapman University, One University Drive, Orange, California 92866, United States
| | - Gennady M Verkhivker
- Graduate Program in Computational and Data Sciences, Keck Center for Science and Engineering, Schmid College of Science and Technology, Chapman University, One University Drive, Orange, California 92866, United States.,Department of Biomedical and Pharmaceutical Sciences, Chapman University School of Pharmacy, Irvine, California 92618, United States
| |
Collapse
|
33
|
Badaczewska-Dawid AE, Kmiecik S, Koliński M. Docking of peptides to GPCRs using a combination of CABS-dock with FlexPepDock refinement. Brief Bioinform 2020; 22:5855394. [PMID: 32520310 PMCID: PMC8138832 DOI: 10.1093/bib/bbaa109] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2020] [Revised: 05/04/2020] [Accepted: 05/06/2020] [Indexed: 12/19/2022] Open
Abstract
The structural description of peptide ligands bound to G protein-coupled receptors (GPCRs) is important for the discovery of new drugs and deeper understanding of the molecular mechanisms of life. Here we describe a three-stage protocol for the molecular docking of peptides to GPCRs using a set of different programs: (1) CABS-dock for docking fully flexible peptides; (2) PD2 method for the reconstruction of atomistic structures from C-alpha traces provided by CABS-dock and (3) Rosetta FlexPepDock for the refinement of protein–peptide complex structures and model scoring. We evaluated the proposed protocol on the set of seven different GPCR–peptide complexes (including one containing a cyclic peptide), for which crystallographic structures are available. We show that CABS-dock produces high resolution models in the sets of top-scored models. These sets of models, after reconstruction to all-atom representation, can be further improved by Rosetta high-resolution refinement and/or minimization, leading in most of the cases to sub-Angstrom accuracy in terms of interface root-mean-square-deviation measure.
Collapse
Affiliation(s)
| | | | - Michał Koliński
- Corresponding author: Michał Koliński, Bioinformatics Laboratory, Mossakowski Medical Research Centre, Polish Academy of Sciences, 5 Pawińskiego St, 02-106 Warsaw, Poland. Tel: (+48) 22 849 93 58; Fax: (+48) 22 668 55 32; E-mail:
| |
Collapse
|
34
|
Drug design by machine-trained elastic networks: predicting Ser/Thr-protein kinase inhibitors' activities. Mol Divers 2020; 25:899-909. [PMID: 32222890 DOI: 10.1007/s11030-020-10074-6] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/21/2020] [Accepted: 03/11/2020] [Indexed: 12/23/2022]
Abstract
An elastic network model (ENM) represents a molecule as a matrix of pairwise atomic interactions. Rich in coded information, ENMs are hereby proposed as a novel tool for the prediction of the activity of series of molecules, with widely different chemical structures, but a common biological activity. The new approach is developed and tested using a set of 183 inhibitors of serine/threonine-protein kinase enzyme (Plk3) which is an enzyme implicated in the regulation of cell cycle and tumorigenesis. The elastic network (EN) predictive model is found to exhibit high accuracy and speed compared to descriptor-based machine-trained modeling. EN modeling appears to be a highly promising new tool for the high demands of industrial applications such as drug and material design.
Collapse
|
35
|
Sheik Amamuddy O, Veldman W, Manyumwa C, Khairallah A, Agajanian S, Oluyemi O, Verkhivker GM, Tastan Bishop Ö. Integrated Computational Approaches and Tools forAllosteric Drug Discovery. Int J Mol Sci 2020; 21:E847. [PMID: 32013012 PMCID: PMC7036869 DOI: 10.3390/ijms21030847] [Citation(s) in RCA: 48] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/25/2019] [Revised: 01/20/2020] [Accepted: 01/21/2020] [Indexed: 12/16/2022] Open
Abstract
Understanding molecular mechanisms underlying the complexity of allosteric regulationin proteins has attracted considerable attention in drug discovery due to the benefits and versatilityof allosteric modulators in providing desirable selectivity against protein targets while minimizingtoxicity and other side effects. The proliferation of novel computational approaches for predictingligand-protein interactions and binding using dynamic and network-centric perspectives has ledto new insights into allosteric mechanisms and facilitated computer-based discovery of allostericdrugs. Although no absolute method of experimental and in silico allosteric drug/site discoveryexists, current methods are still being improved. As such, the critical analysis and integration ofestablished approaches into robust, reproducible, and customizable computational pipelines withexperimental feedback could make allosteric drug discovery more efficient and reliable. In this article,we review computational approaches for allosteric drug discovery and discuss how these tools can beutilized to develop consensus workflows for in silico identification of allosteric sites and modulatorswith some applications to pathogen resistance and precision medicine. The emerging realization thatallosteric modulators can exploit distinct regulatory mechanisms and can provide access to targetedmodulation of protein activities could open opportunities for probing biological processes and insilico design of drug combinations with improved therapeutic indices and a broad range of activities.
Collapse
Affiliation(s)
- Olivier Sheik Amamuddy
- Research Unit in Bioinformatics (RUBi), Department of Biochemistry and Microbiology, Rhodes University, Grahamstown 6140, South Africa; (O.S.A.); (W.V.); (C.M.); (A.K.)
| | - Wayde Veldman
- Research Unit in Bioinformatics (RUBi), Department of Biochemistry and Microbiology, Rhodes University, Grahamstown 6140, South Africa; (O.S.A.); (W.V.); (C.M.); (A.K.)
| | - Colleen Manyumwa
- Research Unit in Bioinformatics (RUBi), Department of Biochemistry and Microbiology, Rhodes University, Grahamstown 6140, South Africa; (O.S.A.); (W.V.); (C.M.); (A.K.)
| | - Afrah Khairallah
- Research Unit in Bioinformatics (RUBi), Department of Biochemistry and Microbiology, Rhodes University, Grahamstown 6140, South Africa; (O.S.A.); (W.V.); (C.M.); (A.K.)
| | - Steve Agajanian
- Graduate Program in Computational and Data Sciences, Keck Center for Science and Engineering, Schmid College of Science and Technology, Chapman University, One University Drive, Orange, CA 92866, USA; (S.A.); (O.O.)
| | - Odeyemi Oluyemi
- Graduate Program in Computational and Data Sciences, Keck Center for Science and Engineering, Schmid College of Science and Technology, Chapman University, One University Drive, Orange, CA 92866, USA; (S.A.); (O.O.)
| | - Gennady M. Verkhivker
- Graduate Program in Computational and Data Sciences, Keck Center for Science and Engineering, Schmid College of Science and Technology, Chapman University, One University Drive, Orange, CA 92866, USA; (S.A.); (O.O.)
- Department of Biomedical and Pharmaceutical Sciences, Chapman University School of Pharmacy, Irvine, CA 92618, USA
| | - Özlem Tastan Bishop
- Research Unit in Bioinformatics (RUBi), Department of Biochemistry and Microbiology, Rhodes University, Grahamstown 6140, South Africa; (O.S.A.); (W.V.); (C.M.); (A.K.)
| |
Collapse
|
36
|
Badaczewska-Dawid AE, Kolinski A, Kmiecik S. Computational reconstruction of atomistic protein structures from coarse-grained models. Comput Struct Biotechnol J 2019; 18:162-176. [PMID: 31969975 PMCID: PMC6961067 DOI: 10.1016/j.csbj.2019.12.007] [Citation(s) in RCA: 34] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/10/2019] [Accepted: 12/10/2019] [Indexed: 01/02/2023] Open
Abstract
Three-dimensional protein structures, whether determined experimentally or theoretically, are often too low resolution. In this mini-review, we outline the computational methods for protein structure reconstruction from incomplete coarse-grained to all atomistic models. Typical reconstruction schemes can be divided into four major steps. Usually, the first step is reconstruction of the protein backbone chain starting from the C-alpha trace. This is followed by side-chains rebuilding based on protein backbone geometry. Subsequently, hydrogen atoms can be reconstructed. Finally, the resulting all-atom models may require structure optimization. Many methods are available to perform each of these tasks. We discuss the available tools and their potential applications in integrative modeling pipelines that can transfer coarse-grained information from computational predictions, or experiment, to all atomistic structures.
Collapse
Affiliation(s)
| | | | - Sebastian Kmiecik
- Faculty of Chemistry, Biological and Chemical Research Center, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| |
Collapse
|
37
|
Establishing Computational Approaches Towards Identifying Malarial Allosteric Modulators: A Case Study of Plasmodium falciparum Hsp70s. Int J Mol Sci 2019; 20:ijms20225574. [PMID: 31717270 PMCID: PMC6887781 DOI: 10.3390/ijms20225574] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2019] [Revised: 10/24/2019] [Accepted: 10/27/2019] [Indexed: 02/07/2023] Open
Abstract
Combating malaria is almost a never-ending battle, as Plasmodium parasites develop resistance to the drugs used against them, as observed recently in artemisinin-based combination therapies. The main concern now is if the resistant parasite strains spread from Southeast Asia to Africa, the continent hosting most malaria cases. To prevent catastrophic results, we need to find non-conventional approaches. Allosteric drug targeting sites and modulators might be a new hope for malarial treatments. Heat shock proteins (HSPs) are potential malarial drug targets and have complex allosteric control mechanisms. Yet, studies on designing allosteric modulators against them are limited. Here, we identified allosteric modulators (SANC190 and SANC651) against P. falciparum Hsp70-1 and Hsp70-x, affecting the conformational dynamics of the proteins, delicately balanced by the endogenous ligands. Previously, we established a pipeline to identify allosteric sites and modulators. This study also further investigated alternative approaches to speed up the process by comparing all atom molecular dynamics simulations and dynamic residue network analysis with the coarse-grained (CG) versions of the calculations. Betweenness centrality (BC) profiles for PfHsp70-1 and PfHsp70-x derived from CG simulations not only revealed similar trends but also pointed to the same functional regions and specific residues corresponding to BC profile peaks.
Collapse
|
38
|
Geng H, Chen F, Ye J, Jiang F. Applications of Molecular Dynamics Simulation in Structure Prediction of Peptides and Proteins. Comput Struct Biotechnol J 2019; 17:1162-1170. [PMID: 31462972 PMCID: PMC6709365 DOI: 10.1016/j.csbj.2019.07.010] [Citation(s) in RCA: 43] [Impact Index Per Article: 8.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/02/2019] [Revised: 07/07/2019] [Accepted: 07/23/2019] [Indexed: 12/21/2022] Open
Abstract
Compared with rapid accumulation of protein sequences from high-throughput DNA sequencing, obtaining experimental 3D structures of proteins is still much more difficult, making protein structure prediction (PSP) potentially very useful. Currently, a vast majority of PSP efforts are based on data mining of known sequences, structures and their relationships (informatics-based). However, if closely related template is not available, these methods are usually much less reliable than experiments. They may also be problematic in predicting the structures of naturally occurring or designed peptides. On the other hand, physics-based methods including molecular dynamics (MD) can utilize our understanding of detailed atomic interactions determining biomolecular structures. In this mini-review, we show that all-atom MD can predict structures of cyclic peptides and other peptide foldamers with accuracy similar to experiments. Then, some notable successes in reproducing experimental 3D structures of small proteins through MD simulations (some with replica-exchange) of the folding were summarized. We also describe advancements of MD-based refinement of structure models, and the integration of limited experimental or bioinformatics data into MD-based structure modeling.
Collapse
Affiliation(s)
- Hao Geng
- Lab of Computational Chemistry and Drug Design, State Key Laboratory of Chemical Oncogenomics, Peking University Shenzhen Graduate School, Shenzhen 518055, China
| | - Fangfang Chen
- Guangdong and Shenzhen Key Laboratory of Male Reproductive Medicine and Genetics, Peking University Shenzhen Hospital, Shenzhen PKU-HKUST Medical Center, Shenzhen 518036, China
| | - Jing Ye
- Guangdong and Shenzhen Key Laboratory of Male Reproductive Medicine and Genetics, Peking University Shenzhen Hospital, Shenzhen PKU-HKUST Medical Center, Shenzhen 518036, China
| | - Fan Jiang
- Lab of Computational Chemistry and Drug Design, State Key Laboratory of Chemical Oncogenomics, Peking University Shenzhen Graduate School, Shenzhen 518055, China
- NanoAI Biotech Co.,Ltd., Silicon Valley Compound, Longhua District, Shenzhen 518109, China
- Corresponding author at: Lab of Computational Chemistry and Drug Design, State Key Laboratory of Chemical Oncogenomics, Peking University Shenzhen Graduate School, Shenzhen 518055, China.
| |
Collapse
|
39
|
Karch R, Stocsits C, Ilieva N, Schreiner W. Intramolecular Domain Movements of Free and Bound pMHC and TCR Proteins: A Molecular Dynamics Simulation Study. Cells 2019; 8:cells8070720. [PMID: 31337065 PMCID: PMC6678086 DOI: 10.3390/cells8070720] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2019] [Revised: 07/02/2019] [Accepted: 07/12/2019] [Indexed: 12/11/2022] Open
Abstract
The interaction of antigenic peptides (p) and major histocompatibility complexes (pMHC) with T-cell receptors (TCR) is one of the most important steps during the immune response. Here we present a molecular dynamics simulation study of bound and unbound TCR and pMHC proteins of the LC13-HLA-B*44:05-pEEYLQAFTY complex to monitor differences in relative orientations and movements of domains between bound and unbound states of TCR-pMHC. We generated local coordinate systems for MHC α1- and MHC α2-helices and the variable T-cell receptor regions TCR Vα and TCR Vβ and monitored changes in the distances and mutual orientations of these domains. In comparison to unbound states, we found decreased inter-domain movements in the simulations of bound states. Moreover, increased conformational flexibility was observed for the MHC α2-helix, the peptide, and for the complementary determining regions of the TCR in TCR-unbound states as compared to TCR-bound states.
Collapse
Affiliation(s)
- Rudolf Karch
- Section of Biosimulation and Bioinformatics, Center for Medical Statistics, Informatics and Intelligent Systems (CeMSIIS), Medical University of Vienna, Spitalgasse 23, A-1090 Vienna, Austria
| | - Claudia Stocsits
- Section of Biosimulation and Bioinformatics, Center for Medical Statistics, Informatics and Intelligent Systems (CeMSIIS), Medical University of Vienna, Spitalgasse 23, A-1090 Vienna, Austria
| | - Nevena Ilieva
- Institute of Information and Communication Technologies (IICT), Bulgarian Academy of Sciences, Acad. G. Bonchev Str., Block 25A, 1113 Sofia, Bulgaria
- CERN-TH, Esplanade des Particules 1, 1211 Geneva, Switzerland
| | - Wolfgang Schreiner
- Section of Biosimulation and Bioinformatics, Center for Medical Statistics, Informatics and Intelligent Systems (CeMSIIS), Medical University of Vienna, Spitalgasse 23, A-1090 Vienna, Austria.
| |
Collapse
|
40
|
de Oliveira AA, Faustino J, de Lima ME, Menezes R, Nunes KP. Unveiling the Interplay between the TLR4/MD2 Complex and HSP70 in the Human Cardiovascular System: A Computational Approach. Int J Mol Sci 2019; 20:E3121. [PMID: 31247943 PMCID: PMC6651210 DOI: 10.3390/ijms20133121] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/02/2019] [Revised: 06/19/2019] [Accepted: 06/19/2019] [Indexed: 02/06/2023] Open
Abstract
While precise mechanisms underlying cardiovascular diseases (CVDs) are still not fully understood, previous studies suggest that the innate immune system, through Toll-like receptor 4 (TLR4), plays a crucial part in the pathways leading to these diseases, mainly because of its interplay with endogenous molecules. The Heat-shock protein 70 family (HSP70-70kDa) is of particular interest in cardiovascular tissues as it may have dual effects when interacting with TLR4 pathways. Although the hypothesis of the HSP70 family members acting as TLR4 ligands is becoming widely accepted, to date no co-crystal structure of this complex is available and it is still unknown whether this process requires the co-adaptor MD2. In this study, we aimed at investigating the interplay between the TLR4/MD2 complex and HSP70 family members in the human cardiovascular system through transcriptomic data analysis and at proposing a putative interaction model between these proteins. We report compelling evidence of correlated expression levels between TLR4 and MD2 with HSP70 cognate family members, especially in heart tissue. In our molecular docking simulations, we found that HSP70 in the ATP-bound state presents a better docking score towards the TLR4/MD2 complex compared to the ADP-bound state (-22.60 vs. -10.29 kcal/mol, respectively). Additionally, we show via a proximity ligation assay for HSP70 and TLR4, that cells stimulated with ATP have higher formation of fluorescent spots and that MD2 might be required for the complexation of these proteins. The insights provided by our computational approach are potential scaffolds for future in vivo studies investigating the interplay between the TLR4/MD2 complex and HSP70 family members in the cardiovascular system.
Collapse
Affiliation(s)
- Amanda Almeida de Oliveira
- Department of Biomedical and Chemical Engineering and Sciences, Florida Institute of Technology, Melbourne, FL 32901, USA
| | - Josemar Faustino
- Department of Computer Engineering and Sciences, Florida Institute of Technology, Melbourne, FL 32901, USA
| | - Maria Elena de Lima
- Grupo Santa Casa de Belo Horizonte, Programa de Pós-graduação em Ciências da Saúde, Biomedicina e Medicina, Ensino e Pesquisa da Santa Casa de Belo Horizonte, Belo Horizonte, MG 30150-240, Brazil
| | - Ronaldo Menezes
- Department of Computer Science, University of Exeter, Exeter EX4 4PY, UK
| | - Kenia Pedrosa Nunes
- Department of Biomedical and Chemical Engineering and Sciences, Florida Institute of Technology, Melbourne, FL 32901, USA.
| |
Collapse
|
41
|
Sanbonmatsu KY. Large-scale simulations of nucleoprotein complexes: ribosomes, nucleosomes, chromatin, chromosomes and CRISPR. Curr Opin Struct Biol 2019; 55:104-113. [PMID: 31125796 DOI: 10.1016/j.sbi.2019.03.004] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2019] [Accepted: 03/01/2019] [Indexed: 12/11/2022]
Abstract
Recent advances in biotechnology such as Hi-C, CRISPR/Cas9 and ribosome display have placed nucleoprotein complexes at center stage. Understanding the structural dynamics of these complexes aids in optimizing protocols and interpreting data for these new technologies. The integration of simulation and experiment has helped advance mechanistic understanding of these systems. Coarse-grained simulations, reduced-description models, and explicit solvent molecular dynamics simulations yield useful complementary perspectives on nucleoprotein complex structural dynamics. When combined with Hi-C, cryo-EM, and single molecule measurements, these simulations integrate disparate forms of experimental data into a coherent mechanism.
Collapse
|
42
|
Ciemny MP, Badaczewska-Dawid AE, Pikuzinska M, Kolinski A, Kmiecik S. Modeling of Disordered Protein Structures Using Monte Carlo Simulations and Knowledge-Based Statistical Force Fields. Int J Mol Sci 2019; 20:E606. [PMID: 30708941 PMCID: PMC6386871 DOI: 10.3390/ijms20030606] [Citation(s) in RCA: 33] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2018] [Revised: 01/23/2019] [Accepted: 01/29/2019] [Indexed: 12/20/2022] Open
Abstract
The description of protein disordered states is important for understanding protein folding mechanisms and their functions. In this short review, we briefly describe a simulation approach to modeling protein interactions, which involve disordered peptide partners or intrinsically disordered protein regions, and unfolded states of globular proteins. It is based on the CABS coarse-grained protein model that uses a Monte Carlo (MC) sampling scheme and a knowledge-based statistical force field. We review several case studies showing that description of protein disordered states resulting from CABS simulations is consistent with experimental data. The case studies comprise investigations of protein⁻peptide binding and protein folding processes. The CABS model has been recently made available as the simulation engine of multiscale modeling tools enabling studies of protein⁻peptide docking and protein flexibility. Those tools offer customization of the modeling process, driving the conformational search using distance restraints, reconstruction of selected models to all-atom resolution, and simulation of large protein systems in a reasonable computational time. Therefore, CABS can be combined in integrative modeling pipelines incorporating experimental data and other modeling tools of various resolution.
Collapse
Affiliation(s)
- Maciej Pawel Ciemny
- Faculty of Chemistry, Biological and Chemical Research Center, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland.
- Faculty of Physics, University of Warsaw, Pasteura 5, 02-093 Warsaw, Poland.
| | | | - Monika Pikuzinska
- Faculty of Chemistry, Biological and Chemical Research Center, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland.
| | - Andrzej Kolinski
- Faculty of Chemistry, Biological and Chemical Research Center, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland.
| | - Sebastian Kmiecik
- Faculty of Chemistry, Biological and Chemical Research Center, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland.
| |
Collapse
|