1
|
Evans D, Sheraz S, Lau A. SARS-CoV-2 3CLPro Dihedral Angles Reveal Allosteric Signaling. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.05.22.595309. [PMID: 38826232 PMCID: PMC11142162 DOI: 10.1101/2024.05.22.595309] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/04/2024]
Abstract
In allosteric proteins, identifying the pathways that signals take from allosteric ligand-binding sites to enzyme active sites or binding pockets and interfaces remains challenging. This avenue of research is motivated by the goals of understanding particular macromolecular systems of interest and creating general methods for their study. An especially important protein that is the subject of many investigations in allostery is the SARS-CoV-2 main protease (Mpro), which is necessary for coronaviral replication. It is both an attractive drug target and, due to intense interest in it for the development of pharmaceutical compounds, a gauge of the state-of-the-art approaches in studying protein inhibition. Here we develop a computational method for characterizing protein allostery and use it to study Mpro. We propose a role of the protein's C-terminal tail in allosteric modulation and warn of unintuitive traps that can plague studies of the role of protein dihedrals angles in transmitting allosteric signals.
Collapse
Affiliation(s)
- Daniel Evans
- Department of Biophysics and Biophysical Chemistry, Johns Hopkins University School of Medicine, Baltimore, Maryland 21205, United States
| | - Samreen Sheraz
- Department of Biophysics and Biophysical Chemistry, Johns Hopkins University School of Medicine, Baltimore, Maryland 21205, United States
| | - Albert Lau
- Department of Biophysics and Biophysical Chemistry, Johns Hopkins University School of Medicine, Baltimore, Maryland 21205, United States
| |
Collapse
|
2
|
Ose NJ, Campitelli P, Modi T, Kazan IC, Kumar S, Ozkan SB. Some mechanistic underpinnings of molecular adaptations of SARS-COV-2 spike protein by integrating candidate adaptive polymorphisms with protein dynamics. eLife 2024; 12:RP92063. [PMID: 38713502 PMCID: PMC11076047 DOI: 10.7554/elife.92063] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/08/2024] Open
Abstract
We integrate evolutionary predictions based on the neutral theory of molecular evolution with protein dynamics to generate mechanistic insight into the molecular adaptations of the SARS-COV-2 spike (S) protein. With this approach, we first identified candidate adaptive polymorphisms (CAPs) of the SARS-CoV-2 S protein and assessed the impact of these CAPs through dynamics analysis. Not only have we found that CAPs frequently overlap with well-known functional sites, but also, using several different dynamics-based metrics, we reveal the critical allosteric interplay between SARS-CoV-2 CAPs and the S protein binding sites with the human ACE2 (hACE2) protein. CAPs interact far differently with the hACE2 binding site residues in the open conformation of the S protein compared to the closed form. In particular, the CAP sites control the dynamics of binding residues in the open state, suggesting an allosteric control of hACE2 binding. We also explored the characteristic mutations of different SARS-CoV-2 strains to find dynamic hallmarks and potential effects of future mutations. Our analyses reveal that Delta strain-specific variants have non-additive (i.e., epistatic) interactions with CAP sites, whereas the less pathogenic Omicron strains have mostly additive mutations. Finally, our dynamics-based analysis suggests that the novel mutations observed in the Omicron strain epistatically interact with the CAP sites to help escape antibody binding.
Collapse
Affiliation(s)
- Nicholas James Ose
- Department of Physics and Center for Biological Physics, Arizona State UniversityTempeUnited States
| | - Paul Campitelli
- Department of Physics and Center for Biological Physics, Arizona State UniversityTempeUnited States
| | - Tushar Modi
- Department of Physics and Center for Biological Physics, Arizona State UniversityTempeUnited States
| | - I Can Kazan
- Department of Physics and Center for Biological Physics, Arizona State UniversityTempeUnited States
| | - Sudhir Kumar
- Institute for Genomics and Evolutionary Medicine, Temple UniversityPhiladelphiaUnited States
- Department of Biology, Temple UniversityPhiladelphiaUnited States
- Center for Genomic Medicine Research, King Abdulaziz UniversityJeddahSaudi Arabia
| | - Sefika Banu Ozkan
- Department of Physics and Center for Biological Physics, Arizona State UniversityTempeUnited States
| |
Collapse
|
3
|
Bhat ZA, Khan MM, Rehman A, Iqbal J, Sanjeev BS, Madhumalar A. MD simulations indicate Omicron P132H of SARS-CoV-2 M pro is a potential allosteric mutant involved in modulating the dynamics of catalytic site entry loop. Int J Biol Macromol 2024; 262:130077. [PMID: 38346625 DOI: 10.1016/j.ijbiomac.2024.130077] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/18/2023] [Revised: 02/07/2024] [Accepted: 02/07/2024] [Indexed: 02/20/2024]
Abstract
The SARS-CoV-2 main protease Mpro, essential for viral replication is an important drug target. It plays a critical role in processing viral polyproteins necessary for viral replication assembly. One of the predominant SARS-CoV-2 Mpro mutations of Omicron variant is Pro132His. Structurally, this mutation site is located ∼22 Å away from the catalytic site. The solved crystal structure of this mutant in complex with inhibitors as well as its reported catalytic efficiency did not show any difference with respect to the wild type. Thus, the mutation was concluded to be non-allosteric. Based on microsecond long MD simulation of the Pro132His mutant and wild type, we show that Pro132His mutation affects the conformational equilibrium with more population of conformational substates having open catalytic site, modulated by the dynamics of the catalytic site entry loop, implying the allosteric nature of this mutation. The structural analysis indicates that rearrangement of hydrogen bonds between His132 and adjacent residues enhances the dynamics of the linker, which in turn is augmented by the inherent dynamic flexibility of the catalytic pocket entry site due to the presence of charged residues. The altered dynamics leading to loss of secondary structures corroborate well with the reported compromised thermal stability.
Collapse
Affiliation(s)
- Zahoor Ahmad Bhat
- Multidisciplinary Centre for Advanced Research and Studies, Jamia Millia Islamia, New Delhi 110025, India
| | - Mohd Muzammil Khan
- Multidisciplinary Centre for Advanced Research and Studies, Jamia Millia Islamia, New Delhi 110025, India
| | - Ayyub Rehman
- Multidisciplinary Centre for Advanced Research and Studies, Jamia Millia Islamia, New Delhi 110025, India
| | - Jawed Iqbal
- Multidisciplinary Centre for Advanced Research and Studies, Jamia Millia Islamia, New Delhi 110025, India
| | - B S Sanjeev
- Department of Applied Sciences, Indian Institute of Information Technology, Prayagraj -211012, India
| | - Arumugam Madhumalar
- Multidisciplinary Centre for Advanced Research and Studies, Jamia Millia Islamia, New Delhi 110025, India.
| |
Collapse
|
4
|
Ose NJ, Campitelli P, Modi T, Can Kazan I, Kumar S, Banu Ozkan S. Some mechanistic underpinnings of molecular adaptations of SARS-COV-2 spike protein by integrating candidate adaptive polymorphisms with protein dynamics. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2023.09.14.557827. [PMID: 37745560 PMCID: PMC10515954 DOI: 10.1101/2023.09.14.557827] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 09/26/2023]
Abstract
We integrate evolutionary predictions based on the neutral theory of molecular evolution with protein dynamics to generate mechanistic insight into the molecular adaptations of the SARS-COV-2 Spike (S) protein. With this approach, we first identified Candidate Adaptive Polymorphisms (CAPs) of the SARS-CoV-2 Spike protein and assessed the impact of these CAPs through dynamics analysis. Not only have we found that CAPs frequently overlap with well-known functional sites, but also, using several different dynamics-based metrics, we reveal the critical allosteric interplay between SARS-CoV-2 CAPs and the S protein binding sites with the human ACE2 (hACE2) protein. CAPs interact far differently with the hACE2 binding site residues in the open conformation of the S protein compared to the closed form. In particular, the CAP sites control the dynamics of binding residues in the open state, suggesting an allosteric control of hACE2 binding. We also explored the characteristic mutations of different SARS-CoV-2 strains to find dynamic hallmarks and potential effects of future mutations. Our analyses reveal that Delta strain-specific variants have non-additive (i.e., epistatic) interactions with CAP sites, whereas the less pathogenic Omicron strains have mostly additive mutations. Finally, our dynamics-based analysis suggests that the novel mutations observed in the Omicron strain epistatically interact with the CAP sites to help escape antibody binding.
Collapse
Affiliation(s)
- Nicholas J. Ose
- Department of Physics and Center for Biological Physics, Arizona State University, Tempe, Arizona, United States of America
| | - Paul Campitelli
- Department of Physics and Center for Biological Physics, Arizona State University, Tempe, Arizona, United States of America
| | - Tushar Modi
- Department of Physics and Center for Biological Physics, Arizona State University, Tempe, Arizona, United States of America
| | - I. Can Kazan
- Department of Physics and Center for Biological Physics, Arizona State University, Tempe, Arizona, United States of America
| | - Sudhir Kumar
- Institute for Genomics and Evolutionary Medicine, Temple University, Philadelphia, Pennsylvania, United States of America
- Department of Biology, Temple University, Philadelphia, Pennsylvania, United States of America
- Center for Genomic Medicine Research, King Abdulaziz University, Jeddah, Saudi Arabia
| | - S. Banu Ozkan
- Department of Physics and Center for Biological Physics, Arizona State University, Tempe, Arizona, United States of America
| |
Collapse
|
5
|
Ose NJ, Campitelli P, Patel R, Kumar S, Ozkan SB. Protein dynamics provide mechanistic insights about epistasis among common missense polymorphisms. Biophys J 2023; 122:2938-2947. [PMID: 36726312 PMCID: PMC10398253 DOI: 10.1016/j.bpj.2023.01.037] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2022] [Revised: 12/20/2022] [Accepted: 01/26/2023] [Indexed: 02/03/2023] Open
Abstract
Sequencing of the protein coding genome has revealed many different missense mutations of human proteins and different population frequencies of corresponding haplotypes, which consist of different sets of those mutations. Here, we present evidence for pairwise intramolecular epistasis (i.e., nonadditive interactions) between many such mutations through an analysis of protein dynamics. We suggest that functional compensation for conserving protein dynamics is a likely evolutionary mechanism that maintains high-frequency mutations that are individually nonneutral but epistatically compensating within proteins. This analysis is the first of its type to look at human proteins with specific high population frequency mutations and examine the relationship between mutations that make up that observed high-frequency protein haplotype. Importantly, protein dynamics revealed a separation between high and low frequency haplotypes within a target protein cytochrome P450 2A7, with the high-frequency haplotypes showing behavior closer to the wild-type protein. Common protein haplotypes containing two mutations display dynamic compensation in which one mutation can correct for the dynamic effects of the other. We also utilize a dynamics-based metric, EpiScore, that evaluates the epistatic interactions and allows us to see dynamic compensation within many other proteins.
Collapse
Affiliation(s)
- Nicholas J Ose
- Department of Physics and Center for Biological Physics, Arizona State University, Tempe, Arizona
| | - Paul Campitelli
- Department of Physics and Center for Biological Physics, Arizona State University, Tempe, Arizona
| | - Ravi Patel
- Institute for Genomics and Evolutionary Medicine, Temple University, Philadelphia, Pennsylvania; Department of Biology, Temple University, Philadelphia, Pennsylvania
| | - Sudhir Kumar
- Institute for Genomics and Evolutionary Medicine, Temple University, Philadelphia, Pennsylvania; Department of Biology, Temple University, Philadelphia, Pennsylvania; Center for Genomic Medicine Research, King Abdulaziz University, Jeddah, Saudi Arabia.
| | - S Banu Ozkan
- Department of Physics and Center for Biological Physics, Arizona State University, Tempe, Arizona.
| |
Collapse
|
6
|
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
|
7
|
Essa AF, El-Hawary SS, Kubacy TM, El-Din A M El-Khrisy E, El-Desoky AH, Elshamy AI, Younis IY. Integration of LC/MS, NMR and Molecular Docking for Profiling of Bioactive Diterpenes from Euphorbia mauritanica L. with in Vitro Anti-SARS-CoV-2 Activity. Chem Biodivers 2023; 20:e202200918. [PMID: 36602020 DOI: 10.1002/cbdv.202200918] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/28/2022] [Revised: 01/04/2023] [Accepted: 01/04/2023] [Indexed: 01/06/2023]
Abstract
In spite of tremendous efforts exerted in the management of COVID-19, the absence of specific treatments and the prevalence of delayed and long-term complications termed post-COVID syndrome still urged all concerned researchers to develop a potent inhibitor of SARS-Cov-2. The hydromethanolic extracts of different parts of E. mauritanica were in vitro screened for anti-SARS-Cov-2 activity. Then, using an integrated strategy of LC/MS/MS, molecular networking and NMR, the chemical profile of the active extract was determined. To determine the optimum target for these compounds, docking experiments of the active extract's identified compounds were conducted at several viral targets. The leaves extract showed the best inhibitory effect with IC50 8.231±0.04 μg/ml. The jatrophane diterpenes were provisionally annotated as the primary metabolites of the bioactive leaves extract based on multiplex of LC/MS/MS, molecular network, and NMR. In silico studies revealed the potentiality of the compounds in the most active extract to 3CLpro, where compound 20 showed the best binding affinity. Further attention should be paid to the isolation of various jatrophane diterpenes from Euphorbia and evaluating their effects on SARS-Cov-2 and its molecular targets.
Collapse
Affiliation(s)
- Ahmed F Essa
- Chemistry of Natural Compounds Department, National Research Center, 33 El Bohouth St., Dokki, Giza, 12622, Egypt
| | - Seham S El-Hawary
- Pharmacognosy Department, Faculty of Pharmacy, Cairo University, Cairo, 12613, Egypt
| | - Tahia M Kubacy
- Chemistry of Natural Compounds Department, National Research Center, 33 El Bohouth St., Dokki, Giza, 12622, Egypt
| | - Ezz El-Din A M El-Khrisy
- Chemistry of Natural Compounds Department, National Research Center, 33 El Bohouth St., Dokki, Giza, 12622, Egypt
| | - Ahmed H El-Desoky
- Pharmacognosy Department, National Research Center, 33 El Bohouth St., Dokki, Giza, 12622, Egypt
| | - Abdelsamed I Elshamy
- Chemistry of Natural Compounds Department, National Research Center, 33 El Bohouth St., Dokki, Giza, 12622, Egypt
| | - Inas Y Younis
- Pharmacognosy Department, Faculty of Pharmacy, Cairo University, Cairo, 12613, Egypt
| |
Collapse
|
8
|
Wordom update 2: A user-friendly program for the analysis of molecular structures and conformational ensembles. Comput Struct Biotechnol J 2023; 21:1390-1402. [PMID: 36817953 PMCID: PMC9929209 DOI: 10.1016/j.csbj.2023.01.026] [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: 12/14/2022] [Revised: 01/19/2023] [Accepted: 01/20/2023] [Indexed: 01/29/2023] Open
Abstract
We present the second update of Wordom, a user-friendly and efficient program for manipulation and analysis of conformational ensembles from molecular simulations. The actual update expands some of the existing modules and adds 21 new modules to the update 1 published in 2011. The new adds can be divided into three sets that: 1) analyze atomic fluctuations and structural communication; 2) explore ion-channel conformational dynamics and ionic translocation; and 3) compute geometrical indices of structural deformation. Set 1 serves to compute correlations of motions, find geometrically stable domains, identify a dynamically invariant core, find changes in domain-domain separation and mutual orientation, perform wavelet analysis of large-scale simulations, process the output of principal component analysis of atomic fluctuations, perform functional mode analysis, infer regions of mechanical rigidity, analyze overall fluctuations, and perform the perturbation response scanning. Set 2 includes modules specific for ion channels, which serve to monitor the pore radius as well as water or ion fluxes, and measure functional collective motions like receptor twisting or tilting angles. Finally, set 3 includes tools to monitor structural deformations by computing angles, perimeter, area, volume, β-sheet curvature, radial distribution function, and center of mass. The ring perception module is also included, helpful to monitor supramolecular self-assemblies. This update places Wordom among the most suitable, complete, user-friendly, and efficient software for the analysis of biomolecular simulations. The source code of Wordom and the relative documentation are available under the GNU general public license at http://wordom.sf.net.
Collapse
|
9
|
Yildirim A, Tekpinar M. Building Quantitative Bridges between Dynamics and Sequences of SARS-CoV-2 Main Protease and a Diverse Set of Thirty-Two Proteins. J Chem Inf Model 2023; 63:9-19. [PMID: 36513349 DOI: 10.1021/acs.jcim.2c01206] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
Proteases are major drug targets for many viral diseases. However, mutations can render several antiprotease drugs inefficient rapidly even though these mutations may not alter protein structures significantly. Understanding relations between quickly mutating residues, protease structures, and the dynamics of the proteases is crucial for designing potent drugs. Due to this reason, we studied relations between the evolutionary information on residues in the amino acid sequences and protein dynamics for SARS-CoV-2 main protease. More precisely, we analyzed three dynamical quantities (Schlitter entropy, root-mean-square fluctuations, and dynamical flexibility index) and their relation to the amino acid conservation extracted from multiple sequence alignments of the main protease. We showed that a quantifiable similarity can be built between a sequence-based quantity called Jensen-Shannon conservation and those three dynamical quantities. We validated this similarity for a diverse set of 32 different proteins, other than the SARS-CoV-2 main protease. We believe that establishing these kinds of quantitative bridges will have larger implications for all viral proteases as well as all proteins.
Collapse
Affiliation(s)
- Ahmet Yildirim
- Department of Biology, Siirt University, 56100Siirt, Turkey
| | - Mustafa Tekpinar
- CNRS, IBPS, Laboratoire de Biologie Computationnelle et Quantitative - UMR 7238, Sorbonne University, 75005Paris, France
| |
Collapse
|