1
|
Dziadek ŁJ, Sieradzan AK, Czaplewski C, Zalewski M, Banaś F, Toczek M, Nisterenko W, Grudinin S, Liwo A, Giełdoń A. Assessment of Four Theoretical Approaches to Predict Protein Flexibility in the Crystal Phase and Solution. J Chem Theory Comput 2024. [PMID: 39171852 DOI: 10.1021/acs.jctc.4c00754] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/23/2024]
Abstract
In this paper, we evaluated the ability of four coarse-grained methods to predict protein flexible regions with potential biological importance, UNRES-flex, UNRES-DSSP-flex (based on the united residue model of polypeptide chains without and with secondary structure restraints, respectively), CABS-flex (based on the C-α, C-β, and side chain model), and nonlinear rigid block normal mode analysis (NOLB) with a set of 100 protein structures determined by NMR spectroscopy or X-ray crystallography, with all secondary structure types. End regions with high fluctuations were excluded from analysis. The Pearson and Spearman correlation coefficients were used to quantify the conformity between the calculated and experimental fluctuation profiles, the latter determined from NMR ensembles and X-ray B-factors, respectively. For X-ray structures (corresponding to proteins in a crowded environment), NOLB resulted in the best agreement between the predicted and experimental fluctuation profiles, while for NMR structures (corresponding to proteins in solution), the ranking of performance is CABS-flex > UNRES-DSSP-flex > UNRES-flex > NOLB; however, CABS-flex sometimes exaggerated the extent of small fluctuations, as opposed to UNRES-DSSP-flex.
Collapse
Affiliation(s)
- Ł J Dziadek
- Faculty of Chemistry, University of Gdansk, ul. Wita-Stwosza 63, 80-308 Gdańsk, Poland
| | - A K Sieradzan
- Faculty of Chemistry, University of Gdansk, ul. Wita-Stwosza 63, 80-308 Gdańsk, Poland
| | - C Czaplewski
- Faculty of Chemistry, University of Gdansk, ul. Wita-Stwosza 63, 80-308 Gdańsk, Poland
- School of Computational Sciences, Korea Institute for Advanced Study, 85 Hoegiro, Dongdaemun-gu, Seoul 02455, Republic of Korea
| | - M Zalewski
- Faculty of Chemistry, University of Gdansk, ul. Wita-Stwosza 63, 80-308 Gdańsk, Poland
| | - F Banaś
- Faculty of Chemistry, University of Gdansk, ul. Wita-Stwosza 63, 80-308 Gdańsk, Poland
| | - M Toczek
- Faculty of Chemistry, University of Gdansk, ul. Wita-Stwosza 63, 80-308 Gdańsk, Poland
| | - W Nisterenko
- Faculty of Chemistry, University of Gdansk, ul. Wita-Stwosza 63, 80-308 Gdańsk, Poland
| | - S Grudinin
- LJK, University Grenoble Alpes, CNRS, Grenoble INP, F-38000 Grenoble, France
| | - A Liwo
- Faculty of Chemistry, University of Gdansk, ul. Wita-Stwosza 63, 80-308 Gdańsk, Poland
| | - A Giełdoń
- Faculty of Chemistry, University of Gdansk, ul. Wita-Stwosza 63, 80-308 Gdańsk, Poland
| |
Collapse
|
2
|
Naveed M, Mahmood S, Aziz T, Azeem A, Hussain I, Waseem M, Ali A, Alharbi M, Alshammari A, Alasmari AF. Designing a novel chimeric multi-epitope vaccine subunit against Staphylococcus argenteus through artificial intelligence approach integrating pan-genome analysis, in vitro identification, and immunogenicity profiling. J Biomol Struct Dyn 2023:1-16. [PMID: 37695632 DOI: 10.1080/07391102.2023.2256881] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/05/2023] [Accepted: 09/04/2023] [Indexed: 09/12/2023]
Abstract
Staphylococcus argenteus is a newly identified pathogen that causes respiratory tract infections, skin infections, such as cellulitis, abscesses, and impetigo, and currently, there is no licensed vaccine available against it. To develop a vaccine against S. argenteus, a bacterial pan-genome analysis was applied to identify potential vaccine candidates. A total of 4908 core proteins were retrieved and utilized for identifying four proteins, including SG38 Panton-Valentine leukocidin LukS-PV protein, SG62 staphylococcal enterotoxin type A protein, SG39 enterotoxin B protein, and SG43 enterotoxin type C3 protein as potential vaccine candidates. Epitopes were predicted for these proteins using different types of B and T-cell epitope prediction tools, and only those with a non-toxic profile, antigenic, non-allergenic, and immunogenic were selected. The selected epitopes were linked to each other to form a multi-epitope vaccine construct, which was further linked to the PADRE sequence (AKFVAAWTLKAAA) and 50s ribosomal L7/L12 protein to enhance the vaccine's antigenicity. The three-dimensional structure of the vaccine construct was assessed to determine its binding affinity with key Toll-like receptor 9 (TLR-9) and Toll-like receptor 5 (TLR-5) immune cell receptors. Our findings demonstrate that the vaccine exhibits favorable binding interactions with these immune cell receptors, indicating its potential efficacy. Molecular dynamic simulations further confirmed the accessibility of vaccine epitopes to the host immune system, substantiating its ability to elicit protective immune responses. Taken together, this study highlights the promising candidacy of the modeled vaccine construct for future in vivo and in vitro experimental investigations.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Muhammad Naveed
- Department of Biotechnology, Faculty of Life Sciences, University of Central Punjab, Lahore, Pakistan
| | - Sarmad Mahmood
- Department of Biotechnology, Faculty of Life Sciences, University of Central Punjab, Lahore, Pakistan
| | - Tariq Aziz
- Laboratory of Animal Health, Food Hygiene and Quality, Department of Agriculture, University of Ioannina, Arta, Greece
| | - Arooj Azeem
- Department of Biotechnology, Faculty of Life Sciences, University of Central Punjab, Lahore, Pakistan
| | - Ibrar Hussain
- Punjab Health Facility Management Company, Chakwal, Pakistan
| | - Muhammad Waseem
- Department of Biotechnology, Faculty of Life Sciences, University of Central Punjab, Lahore, Pakistan
| | - Ayaz Ali
- Department of Biotechnology, University of Malakand, Chakdara, Pakistan
| | - Metab Alharbi
- Department of Pharmacology and Toxicology, College of Pharmacy, King Saud University, Riyadh, Saudi Arabia
| | - Abdulrahman Alshammari
- Department of Pharmacology and Toxicology, College of Pharmacy, King Saud University, Riyadh, Saudi Arabia
| | - Abdullah F Alasmari
- Department of Pharmacology and Toxicology, College of Pharmacy, King Saud University, Riyadh, Saudi Arabia
| |
Collapse
|
3
|
Chkioua L, Amri Y, Sahli C, Rhouma FB, Chehida AB, Tebib N, Messaoud T, Abdennebi HB, Laradi S. Identification of mutations that causes glucose-6-phosphate transporter defect in tunisian patients with glycogenosis type 1b. Diabetol Metab Syndr 2023; 15:86. [PMID: 37118808 PMCID: PMC10142411 DOI: 10.1186/s13098-023-01065-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 12/23/2022] [Accepted: 04/18/2023] [Indexed: 04/30/2023] Open
Abstract
BACKGROUND Glycogen storage disease type 1b (GSD1b) is an autosomal recessive lysosomal storage disease caused by defective glucose-6-phosphate transporter encoded by SLC37A4 leading to the accumulation of glycogen in various tissues. The high rate of consanguineous marriages in Tunisian population provides an ideal environment to facilitate the identification of homozygous pathogenic mutations. We aimed to determine the clinical and genetic profiles of patients with GSD1b to evaluate SLC37A4 mutations spectrum in Tunisian patients. METHODS All exons and flanking intron regions of SLC37A4 gene were screened by direct sequencing to identify mutations and polymorphisms in three unrelated families with GSD1b. Bioinformatics tools were then used to predict the impacts of identified mutations on the structure and function of protein in order to propose a function-structure relationship of the G6PT1 protein. RESULTS Three patients (MT, MB and SI) in Families I, II and III who had the severe phenotype were homoallelic for the two identified mutations: p.R300H (famillies I, II) and p.W393X (Family III), respectively. One of the alterations was a missense mutation p.R300H of exon 6 in SLC37A4 gene. The analysis of the protein structure flexibility upon p.R300H mutation using DynaMut tool and CABS-flex 2.0 server showed that the reported mutation increase the molecule flexibility of in the cytosol region and would probably lead to significant conformational changes. CONCLUSION This is the first Tunisian report of SLC37A4 mutations identified in Tunisia causing the glycogenosis type Ib disease. Bioinformatics analysis allowed us to establish an approximate structure-function relationship for the G6PT1 protein, thereby providing better genotype/phenotype correlation knowledge.
Collapse
Affiliation(s)
- Latifa Chkioua
- Research Laboratory of Human Genome and Multifactorial Diseases, Faculty of Pharmacy, University of Monastir, Monastir, Tunisia.
- Faculty of Pharmacy, University of Monastir, Street Avicenne, Monastir, 5000, Tunisia.
| | - Yessine Amri
- Biochemistry Laboratory (LR 00SP03), Bechir Hamza Children's Hospital, Tunis, Tunisia
- Department of Educational Sciences, Higher Institute of Applied Studies in Humanity, University of Jendouba, Le Kef, Tunis, Tunisia
| | - Chayma Sahli
- Biochemistry Laboratory (LR 00SP03), Bechir Hamza Children's Hospital, Tunis, Tunisia
| | - Ferdawes Ben Rhouma
- Research Laboratory of Human Genome and Multifactorial Diseases, Faculty of Pharmacy, University of Monastir, Monastir, Tunisia
| | - Amel Ben Chehida
- Pediatrics Department, La Rabta Hospital, Tunis, Tunisia
- Research Laboratory: LR12SPO2 Investigation and Management of Inherited Metabolic Diseases, Faculty of Medicine of Tunis, Tunis El Manar University, Tunis, Tunisia
| | - Neji Tebib
- Pediatrics Department, La Rabta Hospital, Tunis, Tunisia
| | - Taieb Messaoud
- Biochemistry Laboratory (LR 00SP03), Bechir Hamza Children's Hospital, Tunis, Tunisia
| | - Hassen Ben Abdennebi
- Research Laboratory of Human Genome and Multifactorial Diseases, Faculty of Pharmacy, University of Monastir, Monastir, Tunisia
| | - Sandrine Laradi
- The Auvergne-Rhône-Alpes Regional Branch of the French National Blood System EFS/GIMAP, EA 3064, Saint Etienne, 42100, France
| |
Collapse
|
4
|
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
|
5
|
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
|
6
|
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
|
7
|
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
|
8
|
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]
|
9
|
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
|
10
|
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
|