1
|
Banerjee A, Gosavi S. Potential Self-Peptide Inhibitors of the SARS-CoV-2 Main Protease. J Phys Chem B 2023; 127:855-865. [PMID: 36689738 PMCID: PMC9883841 DOI: 10.1021/acs.jpcb.2c05917] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2022] [Revised: 11/23/2022] [Indexed: 01/24/2023]
Abstract
The SARS-CoV-2 main protease (Mpro) plays an essential role in viral replication, cleaving viral polyproteins into functional proteins. This makes Mpro an important drug target. Mpro consists of an N-terminal catalytic domain and a C-terminal α-helical domain (MproC). Previous studies have shown that peptides derived from a given protein sequence (self-peptides) can affect the folding and, in turn, the function of that protein. Since the SARS-CoV-1 MproC is known to stabilize its Mpro and regulate its function, we hypothesized that SARS-CoV-2 MproC-derived self-peptides may modulate the folding and the function of SARS-CoV-2 Mpro. To test this, we studied the folding of MproC in the presence of various self-peptides using coarse-grained structure-based models and molecular dynamics simulations. In these simulations of MproC and one self-peptide, we found that two self-peptides, the α1-helix and the loop between α4 and α5 (loop4), could replace the equivalent native sequences in the MproC structure. Replacement of either sequence in full-length Mpro should, in principle, be able to perturb Mpro function albeit through different mechanisms. Some general principles for the rational design of self-peptide inhibitors emerge: The simulations show that prefolded self-peptides are more likely to replace native sequences than those which do not possess structure. Additionally, the α1-helix self-peptide is kinetically stable and once inserted rarely exchanges with the native α1-helix, while the loop4 self-peptide is easily replaced by the native loop4, making it less useful for modulating function. In summary, a prefolded α1-derived peptide should be able to inhibit SARS-CoV-2 Mpro function.
Collapse
Affiliation(s)
- Arkadeep Banerjee
- Simons Centre for the Study
of Living Machines, National Centre for Biological Sciences, Tata Institute of Fundamental Research, Bengaluru 560065, India
| | - Shachi Gosavi
- Simons Centre for the Study
of Living Machines, National Centre for Biological Sciences, Tata Institute of Fundamental Research, Bengaluru 560065, India
| |
Collapse
|
2
|
Yadahalli S, Jayanthi LP, Gosavi S. A Method for Assessing the Robustness of Protein Structures by Randomizing Packing Interactions. Front Mol Biosci 2022; 9:849272. [PMID: 35832734 PMCID: PMC9271847 DOI: 10.3389/fmolb.2022.849272] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/05/2022] [Accepted: 04/27/2022] [Indexed: 12/02/2022] Open
Abstract
Many single-domain proteins are not only stable and water-soluble, but they also populate few to no intermediates during folding. This reduces interactions between partially folded proteins, misfolding, and aggregation, and makes the proteins tractable in biotechnological applications. Natural proteins fold thus, not necessarily only because their structures are well-suited for folding, but because their sequences optimize packing and fit their structures well. In contrast, folding experiments on the de novo designed Top7 suggest that it populates several intermediates. Additionally, in de novo protein design, where sequences are designed for natural and new non-natural structures, tens of sequences still need to be tested before success is achieved. Both these issues may be caused by the specific scaffolds used in design, i.e., some protein scaffolds may be more tolerant to packing perturbations and varied sequences. Here, we report a computational method for assessing the response of protein structures to packing perturbations. We then benchmark this method using designed proteins and find that it can identify scaffolds whose folding gets disrupted upon perturbing packing, leading to the population of intermediates. The method can also isolate regions of both natural and designed scaffolds that are sensitive to such perturbations and identify contacts which when present can rescue folding. Overall, this method can be used to identify protein scaffolds that are more amenable to whole protein design as well as to identify protein regions which are sensitive to perturbations and where further mutations should be avoided during protein engineering.
Collapse
|
3
|
Zerihun MB, Pucci F, Schug A. CoCoNet-boosting RNA contact prediction by convolutional neural networks. Nucleic Acids Res 2021; 49:12661-12672. [PMID: 34871451 PMCID: PMC8682773 DOI: 10.1093/nar/gkab1144] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/23/2020] [Revised: 10/27/2021] [Accepted: 11/05/2021] [Indexed: 11/24/2022] Open
Abstract
Co-evolutionary models such as direct coupling analysis (DCA) in combination with machine learning (ML) techniques based on deep neural networks are able to predict accurate protein contact or distance maps. Such information can be used as constraints in structure prediction and massively increase prediction accuracy. Unfortunately, the same ML methods cannot readily be applied to RNA as they rely on large structural datasets only available for proteins. Here, we demonstrate how the available smaller data for RNA can be used to improve prediction of RNA contact maps. We introduce an algorithm called CoCoNet that is based on a combination of a Coevolutionary model and a shallow Convolutional Neural Network. Despite its simplicity and the small number of trained parameters, the method boosts the positive predictive value (PPV) of predicted contacts by about 70% with respect to DCA as tested by cross-validation of about eighty RNA structures. However, the direct inclusion of the CoCoNet contacts in 3D modeling tools does not result in a proportional increase of the 3D RNA structure prediction accuracy. Therefore, we suggest that the field develops, in addition to contact PPV, metrics which estimate the expected impact for 3D structure modeling tools better. CoCoNet is freely available and can be found at https://github.com/KIT-MBS/coconet.
Collapse
Affiliation(s)
- Mehari B Zerihun
- John von Neumann Institute for Computing, Jülich Supercomputing Centre, Forschungszentrum Jülich, 52428 Jülich, Germany.,Steinbuch Centre for Computing, Karlsruhe Institute of Technology, 76344 Eggenstein-Leopoldshafen, Germany
| | - Fabrizio Pucci
- John von Neumann Institute for Computing, Jülich Supercomputing Centre, Forschungszentrum Jülich, 52428 Jülich, Germany.,Computational Biology and Bioinformatics, Université Libre de Bruxelles 1050, Brussels, Belgium
| | - Alexander Schug
- John von Neumann Institute for Computing, Jülich Supercomputing Centre, Forschungszentrum Jülich, 52428 Jülich, Germany.,Faculty of Biology, University of Duisburg-Essen, 45117 Essen, Germany
| |
Collapse
|
4
|
Terse VL, Gosavi S. The Molecular Mechanism of Domain Swapping of the C-Terminal Domain of the SARS-Coronavirus Main Protease. Biophys J 2020; 120:504-516. [PMID: 33359834 PMCID: PMC7837137 DOI: 10.1016/j.bpj.2020.11.2277] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2020] [Revised: 11/09/2020] [Accepted: 11/24/2020] [Indexed: 11/25/2022] Open
Abstract
In three-dimensional domain swapping, two protein monomers exchange a part of their structures to form an intertwined homodimer, whose subunits resemble the monomer. Several viral proteins domain swap to increase their structural complexity or functional avidity. The main protease (Mpro) of the severe acute respiratory syndrome (SARS) coronavirus proteolyzes viral polyproteins and has been a target for anti-SARS drug design. Domain swapping in the α-helical C-terminal domain of Mpro (MproC) locks Mpro into a hyperactive octameric form that is hypothesized to promote the early stages of viral replication. However, in the absence of a complete molecular understanding of the mechanism of domain swapping, investigations into the biological relevance of this octameric Mpro have stalled. Isolated MproC can exist as a monomer or a domain-swapped dimer. Here, we investigate the mechanism of domain swapping of MproC using coarse-grained structure-based models and molecular dynamics simulations. Our simulations recapitulate several experimental features of MproC folding. Further, we find that a contact between a tryptophan in the MproC domain-swapping hinge and an arginine elsewhere forms early during folding, modulates the folding route, and promotes domain swapping to the native structure. An examination of the sequence and the structure of the tryptophan containing hinge loop shows that it has a propensity to form multiple secondary structures and contacts, indicating that it could be stabilized into either the monomer- or dimer-promoting conformations by mutations or ligand binding. Finally, because all residues in the tryptophan loop are identical in SARS-CoV and SARS-CoV-2, mutations that modulate domain swapping may provide insights into the role of octameric Mpro in the early-stage viral replication of both viruses.
Collapse
Affiliation(s)
- Vishram L Terse
- Simons Centre for the Study of Living Machines, National Centre for Biological Sciences, Tata Institute of Fundamental Research, Bangalore, India
| | - Shachi Gosavi
- Simons Centre for the Study of Living Machines, National Centre for Biological Sciences, Tata Institute of Fundamental Research, Bangalore, India.
| |
Collapse
|
5
|
Terse VL, Gosavi S. The Sensitivity of Computational Protein Folding to Contact Map Perturbations: The Case of Ubiquitin Folding and Function. J Phys Chem B 2018; 122:11497-11507. [PMID: 30234303 DOI: 10.1021/acs.jpcb.8b07409] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Ubiquitin is a small model protein, commonly used in protein folding experiments and simulations. We simulated ubiquitin using a well-tested structure-based model coarse-grained to a Cα level (Cα-SBM) and found that the simulated folding route did not agree with the experimentally observed one. Simulating the Cα-SBM with a cutoff contact map, instead of a screened contact map, switched the folding route with the new route matching the experimental route. Thus, the simulated folding of ubiquitin is sensitive to contact map definition. The screened contact map, which is used in folding simulations because it captures protein folding cooperativity, removes contacts in which the atoms in contact are occluded by a third atom and is less sensitive to the value of the cutoff distance in well-packed regions of the protein. In sparsely packed regions, the larger cutoff distance creates bridging contacts between atoms which are separated by voids. Such contacts do not seem to affect the folding of most proteins, including those of the ubiquitin fold. However, the surface of ubiquitin has several protruding functional side chains which naturally create bridging contacts. Together, our results show that subtle structural features of a protein that may not be apparent by mere observation can be identified by comparing folding simulations of SBMs in which these features are differently encoded. When such structural features are preserved for functional reasons, differences in computational folding can be leveraged to identify functional features. Notably, such features are accessible to a gradation of SBMs even in commonly studied proteins such as ubiquitin.
Collapse
Affiliation(s)
- Vishram L Terse
- Simons Centre for the Study of Living Machines , National Centre for Biological Sciences , Tata Institute of Fundamental Research, Bangalore 560065 , India
| | - Shachi Gosavi
- Simons Centre for the Study of Living Machines , National Centre for Biological Sciences , Tata Institute of Fundamental Research, Bangalore 560065 , India
| |
Collapse
|
6
|
Reinartz I, Sinner C, Nettels D, Stucki-Buchli B, Stockmar F, Panek PT, Jacob CR, Nienhaus GU, Schuler B, Schug A. Simulation of FRET dyes allows quantitative comparison against experimental data. J Chem Phys 2018; 148:123321. [DOI: 10.1063/1.5010434] [Citation(s) in RCA: 34] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Affiliation(s)
- Ines Reinartz
- Department of Physics, Karlsruhe Institute of Technology, Wolfgang-Gaede-Str. 1, 76131 Karlsruhe, Germany
- Steinbuch Centre for Computing, Karlsruhe Institute of Technology, Hermann-von-Helmholtz-Platz 1, 76344 Eggenstein-Leopoldshafen, Germany
| | - Claude Sinner
- Department of Physics, Karlsruhe Institute of Technology, Wolfgang-Gaede-Str. 1, 76131 Karlsruhe, Germany
- Steinbuch Centre for Computing, Karlsruhe Institute of Technology, Hermann-von-Helmholtz-Platz 1, 76344 Eggenstein-Leopoldshafen, Germany
| | - Daniel Nettels
- Department of Biochemistry, University of Zurich, Winterthurerstrasse 190, 8057 Zurich, Switzerland
| | - Brigitte Stucki-Buchli
- Department of Biochemistry, University of Zurich, Winterthurerstrasse 190, 8057 Zurich, Switzerland
| | - Florian Stockmar
- Institute of Applied Physics, Karlsruhe Institute of Technology, Wolfgang-Gaede-Str. 1, 76131 Karlsruhe, Germany
| | - Pawel T. Panek
- Institute of Physical and Theoretical Chemistry, TU Braunschweig, Gaußstraße 17, 38106 Braunschweig, Germany
| | - Christoph R. Jacob
- Institute of Physical and Theoretical Chemistry, TU Braunschweig, Gaußstraße 17, 38106 Braunschweig, Germany
| | - Gerd Ulrich Nienhaus
- Institute of Applied Physics, Karlsruhe Institute of Technology, Wolfgang-Gaede-Str. 1, 76131 Karlsruhe, Germany
- HEiKA–Heidelberg Karlsruhe Research Partnership, Karlsruhe Institute of Technology, Hermann-von-Helmholtz-Platz 1, 76344 Eggenstein-Leopoldshafen, Germany
- Institute of Nanotechnology and Institute of Toxicology and Genetics, Karlsruhe Institute of Technology, Hermann-von-Helmholtz-Platz 1, 76344 Eggenstein-Leopoldshafen, Germany
- Department of Physics, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| | - Benjamin Schuler
- Department of Biochemistry, University of Zurich, Winterthurerstrasse 190, 8057 Zurich, Switzerland
- Department of Physics, University of Zurich, Winterthurerstrasse 190, 8057 Zurich, Switzerland
| | - Alexander Schug
- Steinbuch Centre for Computing, Karlsruhe Institute of Technology, Hermann-von-Helmholtz-Platz 1, 76344 Eggenstein-Leopoldshafen, Germany
- John von Neumann Institute for Computing, Jülich Supercomputing Centre, Forschungszentrum Jülich, 52425 Jülich, Germany
| |
Collapse
|
7
|
Yadahalli S, Gosavi S. Packing energetics determine the folding routes of the RNase-H proteins. Phys Chem Chem Phys 2017; 19:9164-9173. [DOI: 10.1039/c6cp08940b] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
The RNase-H proteins show a diverse range of folding routes with structurally distinct folding nuclei.
Collapse
Affiliation(s)
- Shilpa Yadahalli
- National Centre for Biological Sciences
- Tata Institute of Fundamental Research
- Bangalore-560065
- India
- Manipal University
| | - Shachi Gosavi
- National Centre for Biological Sciences
- Tata Institute of Fundamental Research
- Bangalore-560065
- India
| |
Collapse
|