1
|
Raisinghani N, Alshahrani M, Gupta G, Xiao S, Tao P, Verkhivker G. Exploring conformational landscapes and binding mechanisms of convergent evolution for the SARS-CoV-2 spike Omicron variant complexes with the ACE2 receptor using AlphaFold2-based structural ensembles and molecular dynamics simulations. Phys Chem Chem Phys 2024; 26:17720-17744. [PMID: 38869513 DOI: 10.1039/d4cp01372g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/14/2024]
Abstract
In this study, we combined AlphaFold-based approaches for atomistic modeling of multiple protein states and microsecond molecular simulations to accurately characterize conformational ensembles evolution and binding mechanisms of convergent evolution for the SARS-CoV-2 spike Omicron variants BA.1, BA.2, BA.2.75, BA.3, BA.4/BA.5 and BQ.1.1. We employed and validated several different adaptations of the AlphaFold methodology for modeling of conformational ensembles including the introduced randomized full sequence scanning for manipulation of sequence variations to systematically explore conformational dynamics of Omicron spike protein complexes with the ACE2 receptor. Microsecond atomistic molecular dynamics (MD) simulations provide a detailed characterization of the conformational landscapes and thermodynamic stability of the Omicron variant complexes. By integrating the predictions of conformational ensembles from different AlphaFold adaptations and applying statistical confidence metrics we can expand characterization of the conformational ensembles and identify functional protein conformations that determine the equilibrium dynamics for the Omicron spike complexes with the ACE2. Conformational ensembles of the Omicron RBD-ACE2 complexes obtained using AlphaFold-based approaches for modeling protein states and MD simulations are employed for accurate comparative prediction of the binding energetics revealing an excellent agreement with the experimental data. In particular, the results demonstrated that AlphaFold-generated extended conformational ensembles can produce accurate binding energies for the Omicron RBD-ACE2 complexes. The results of this study suggested complementarities and potential synergies between AlphaFold predictions of protein conformational ensembles and MD simulations showing that integrating information from both methods can potentially yield a more adequate characterization of the conformational landscapes for the Omicron RBD-ACE2 complexes. This study provides insights in the interplay between conformational dynamics and binding, showing that evolution of Omicron variants through acquisition of convergent mutational sites may leverage conformational adaptability and dynamic couplings between key binding energy hotspots to optimize ACE2 binding affinity and enable immune evasion.
Collapse
Affiliation(s)
- Nishank Raisinghani
- 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.
| | - Mohammed Alshahrani
- 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.
| | - Grace Gupta
- 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.
| | - Sian Xiao
- Department of Chemistry, Center for Research Computing, Center for Drug Discovery, Design, and Delivery (CD4), Southern Methodist University, Dallas, Texas, 75275, USA
| | - Peng Tao
- Department of Chemistry, Center for Research Computing, Center for Drug Discovery, Design, and Delivery (CD4), Southern Methodist University, Dallas, Texas, 75275, USA
| | - 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.
- Department of Biomedical and Pharmaceutical Sciences, Chapman University School of Pharmacy, Irvine, CA 92618, USA
| |
Collapse
|
2
|
Raisinghani N, Alshahrani M, Gupta G, Xiao S, Tao P, Verkhivker G. AlphaFold2 Predictions of Conformational Ensembles and Atomistic Simulations of the SARS-CoV-2 Spike XBB Lineages Reveal Epistatic Couplings between Convergent Mutational Hotspots that Control ACE2 Affinity. J Phys Chem B 2024; 128:4696-4715. [PMID: 38696745 DOI: 10.1021/acs.jpcb.4c01341] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/04/2024]
Abstract
In this study, we combined AlphaFold-based atomistic structural modeling, microsecond molecular simulations, mutational profiling, and network analysis to characterize binding mechanisms of the SARS-CoV-2 spike protein with the host receptor ACE2 for a series of Omicron XBB variants including XBB.1.5, XBB.1.5+L455F, XBB.1.5+F456L, and XBB.1.5+L455F+F456L. AlphaFold-based structural and dynamic modeling of SARS-CoV-2 Spike XBB lineages can accurately predict the experimental structures and characterize conformational ensembles of the spike protein complexes with the ACE2. Microsecond molecular dynamics simulations identified important differences in the conformational landscapes and equilibrium ensembles of the XBB variants, suggesting that combining AlphaFold predictions of multiple conformations with molecular dynamics simulations can provide a complementary approach for the characterization of functional protein states and binding mechanisms. Using the ensemble-based mutational profiling of protein residues and physics-based rigorous calculations of binding affinities, we identified binding energy hotspots and characterized the molecular basis underlying epistatic couplings between convergent mutational hotspots. Consistent with the experiments, the results revealed the mediating role of the Q493 hotspot in the synchronization of epistatic couplings between L455F and F456L mutations, providing a quantitative insight into the energetic determinants underlying binding differences between XBB lineages. We also proposed a network-based perturbation approach for mutational profiling of allosteric communications and uncovered the important relationships between allosteric centers mediating long-range communication and binding hotspots of epistatic couplings. The results of this study support a mechanism in which the binding mechanisms of the XBB variants may be determined by epistatic effects between convergent evolutionary hotspots that control ACE2 binding.
Collapse
Affiliation(s)
- Nishank Raisinghani
- 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
| | - Mohammed Alshahrani
- 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
| | - Grace Gupta
- 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
| | - Sian Xiao
- Department of Chemistry, Center for Research Computing, Center for Drug Discovery, Design, and Delivery (CD4), Southern Methodist University, Dallas, Texas 75275, United States
| | - Peng Tao
- Department of Chemistry, Center for Research Computing, Center for Drug Discovery, Design, and Delivery (CD4), Southern Methodist University, Dallas, Texas 75275, United States
| | - Gennady 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
| |
Collapse
|
3
|
Božič A, Podgornik R. Changes in total charge on spike protein of SARS-CoV-2 in emerging lineages. BIOINFORMATICS ADVANCES 2024; 4:vbae053. [PMID: 38645718 PMCID: PMC11031363 DOI: 10.1093/bioadv/vbae053] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/08/2023] [Revised: 03/13/2024] [Accepted: 04/04/2024] [Indexed: 04/23/2024]
Abstract
Motivation Charged amino acid residues on the spike protein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) have been shown to influence its binding to different cell surface receptors, its non-specific electrostatic interactions with the environment, and its structural stability and conformation. It is therefore important to obtain a good understanding of amino acid mutations that affect the total charge on the spike protein which have arisen across different SARS-CoV-2 lineages during the course of the virus' evolution. Results We analyse the change in the number of ionizable amino acids and the corresponding total charge on the spike proteins of almost 2200 SARS-CoV-2 lineages that have emerged over the span of the pandemic. Our results show that the previously observed trend toward an increase in the positive charge on the spike protein of SARS-CoV-2 variants of concern has essentially stopped with the emergence of the early omicron variants. Furthermore, recently emerged lineages show a greater diversity in terms of their composition of ionizable amino acids. We also demonstrate that the patterns of change in the number of ionizable amino acids on the spike protein are characteristic of related lineages within the broader clade division of the SARS-CoV-2 phylogenetic tree. Due to the ubiquity of electrostatic interactions in the biological environment, our findings are relevant for a broad range of studies dealing with the structural stability of SARS-CoV-2 and its interactions with the environment. Availability and implementation The data underlying the article are available in the Supplementary material.
Collapse
Affiliation(s)
- Anže Božič
- Department of Theoretical Physics, Jožef Stefan Institute, Ljubljana 1000, Slovenia
| | - Rudolf Podgornik
- School of Physical Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
- Kavli Institute for Theoretical Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
- Wenzhou Institute, University of Chinese Academy of Sciences, Wenzhou 325001, China
| |
Collapse
|
4
|
Raisinghani N, Alshahrani M, Gupta G, Xiao S, Tao P, Verkhivker G. Predicting Functional Conformational Ensembles and Binding Mechanisms of Convergent Evolution for SARS-CoV-2 Spike Omicron Variants Using AlphaFold2 Sequence Scanning Adaptations and Molecular Dynamics Simulations. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.04.02.587850. [PMID: 38617283 PMCID: PMC11014522 DOI: 10.1101/2024.04.02.587850] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/16/2024]
Abstract
In this study, we combined AlphaFold-based approaches for atomistic modeling of multiple protein states and microsecond molecular simulations to accurately characterize conformational ensembles and binding mechanisms of convergent evolution for the SARS-CoV-2 Spike Omicron variants BA.1, BA.2, BA.2.75, BA.3, BA.4/BA.5 and BQ.1.1. We employed and validated several different adaptations of the AlphaFold methodology for modeling of conformational ensembles including the introduced randomized full sequence scanning for manipulation of sequence variations to systematically explore conformational dynamics of Omicron Spike protein complexes with the ACE2 receptor. Microsecond atomistic molecular dynamic simulations provide a detailed characterization of the conformational landscapes and thermodynamic stability of the Omicron variant complexes. By integrating the predictions of conformational ensembles from different AlphaFold adaptations and applying statistical confidence metrics we can expand characterization of the conformational ensembles and identify functional protein conformations that determine the equilibrium dynamics for the Omicron Spike complexes with the ACE2. Conformational ensembles of the Omicron RBD-ACE2 complexes obtained using AlphaFold-based approaches for modeling protein states and molecular dynamics simulations are employed for accurate comparative prediction of the binding energetics revealing an excellent agreement with the experimental data. In particular, the results demonstrated that AlphaFold-generated extended conformational ensembles can produce accurate binding energies for the Omicron RBD-ACE2 complexes. The results of this study suggested complementarities and potential synergies between AlphaFold predictions of protein conformational ensembles and molecular dynamics simulations showing that integrating information from both methods can potentially yield a more adequate characterization of the conformational landscapes for the Omicron RBD-ACE2 complexes. This study provides insights in the interplay between conformational dynamics and binding, showing that evolution of Omicron variants through acquisition of convergent mutational sites may leverage conformational adaptability and dynamic couplings between key binding energy hotspots to optimize ACE2 binding affinity and enable immune evasion.
Collapse
|
5
|
Raisinghani N, Alshahrani M, Gupta G, Xiao S, Tao P, Verkhivker G. AlphaFold2-Enabled Atomistic Modeling of Structure, Conformational Ensembles, and Binding Energetics of the SARS-CoV-2 Omicron BA.2.86 Spike Protein with ACE2 Host Receptor and Antibodies: Compensatory Functional Effects of Binding Hotspots in Modulating Mechanisms of Receptor Binding and Immune Escape. J Chem Inf Model 2024; 64:1657-1681. [PMID: 38373700 DOI: 10.1021/acs.jcim.3c01857] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/21/2024]
Abstract
The latest wave of SARS-CoV-2 Omicron variants displayed a growth advantage and increased viral fitness through convergent evolution of functional hotspots that work synchronously to balance fitness requirements for productive receptor binding and efficient immune evasion. In this study, we combined AlphaFold2-based structural modeling approaches with atomistic simulations and mutational profiling of binding energetics and stability for prediction and comprehensive analysis of the structure, dynamics, and binding of the SARS-CoV-2 Omicron BA.2.86 spike variant with ACE2 host receptor and distinct classes of antibodies. We adapted several AlphaFold2 approaches to predict both the structure and conformational ensembles of the Omicron BA.2.86 spike protein in the complex with the host receptor. The results showed that the AlphaFold2-predicted structural ensemble of the BA.2.86 spike protein complex with ACE2 can accurately capture the main conformational states of the Omicron variant. Complementary to AlphaFold2 structural predictions, microsecond molecular dynamics simulations reveal the details of the conformational landscape and produced equilibrium ensembles of the BA.2.86 structures that are used to perform mutational scanning of spike residues and characterize structural stability and binding energy hotspots. The ensemble-based mutational profiling of the receptor binding domain residues in the BA.2 and BA.2.86 spike complexes with ACE2 revealed a group of conserved hydrophobic hotspots and critical variant-specific contributions of the BA.2.86 convergent mutational hotspots R403K, F486P, and R493Q. To examine the immune evasion properties of BA.2.86 in atomistic detail, we performed structure-based mutational profiling of the spike protein binding interfaces with distinct classes of antibodies that displayed significantly reduced neutralization against the BA.2.86 variant. The results revealed the molecular basis of compensatory functional effects of the binding hotspots, showing that BA.2.86 lineage may have evolved to outcompete other Omicron subvariants by improving immune evasion while preserving binding affinity with ACE2 via through a compensatory effect of R493Q and F486P convergent mutational hotspots. This study demonstrated that an integrative approach combining AlphaFold2 predictions with complementary atomistic molecular dynamics simulations and robust ensemble-based mutational profiling of spike residues can enable accurate and comprehensive characterization of structure, dynamics, and binding mechanisms of newly emerging Omicron variants.
Collapse
Affiliation(s)
- Nishank Raisinghani
- 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 of America
| | - Mohammed Alshahrani
- 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 of America
| | - Grace Gupta
- 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 of America
| | - Sian Xiao
- Department of Chemistry, Center for Research Computing, Center for Drug Discovery, Design, and Delivery (CD4), Southern Methodist University, Dallas, Texas 75275, United States of America
| | - Peng Tao
- Department of Chemistry, Center for Research Computing, Center for Drug Discovery, Design, and Delivery (CD4), Southern Methodist University, Dallas, Texas 75275, United States of America
| | - Gennady 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 of America
- Department of Biomedical and Pharmaceutical Sciences, Chapman University School of Pharmacy, Irvine, California 92618, United States of America
| |
Collapse
|
6
|
Hristova SH, Zhivkov AM. Three-Dimensional Structural Stability and Local Electrostatic Potential at Point Mutations in Spike Protein of SARS-CoV-2 Coronavirus. Int J Mol Sci 2024; 25:2174. [PMID: 38396850 PMCID: PMC10889838 DOI: 10.3390/ijms25042174] [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: 01/17/2024] [Revised: 01/30/2024] [Accepted: 02/01/2024] [Indexed: 02/25/2024] Open
Abstract
The contagiousness of SARS-CoV-2 β-coronavirus is determined by the virus-receptor electrostatic association of its positively charged spike (S) protein with the negatively charged angiotensin converting enzyme-2 (ACE2 receptor) of the epithelial cells. If some mutations occur, the electrostatic potential on the surface of the receptor-binding domain (RBD) could be altered, and the S-ACE2 association could become stronger or weaker. The aim of the current research is to investigate whether point mutations can noticeably alter the electrostatic potential on the RBD and the 3D stability of the S1-subunit of the S-protein. For this purpose, 15 mutants with different hydrophilicity and electric charge (positive, negative, or uncharged) of the substituted and substituting amino acid residues, located on the RBD at the S1-ACE2 interface, are selected, and the 3D structure of the S1-subunit is reconstructed on the base of the crystallographic structure of the S-protein of the wild-type strain and the amino acid sequence of the unfolded polypeptide chain of the mutants. Then, the Gibbs free energy of folding, isoelectric point, and pH-dependent surface electrostatic potential of the S1-subunit are computed using programs for protein electrostatics. The results show alterations in the local electrostatic potential in the vicinity of the mutant amino acid residue, which can influence the S-ACE2 association. This approach allows prediction of the relative infectivity, transmissibility, and contagiousness (at equal social immune status) of new SARS-CoV-2 mutants by reconstruction of the 3D structure of the S1-subunit and calculation of the surface electrostatic potential.
Collapse
Affiliation(s)
- Svetlana H. Hristova
- Department of Medical Physics and Biophysics, Medical Faculty, Medical University—Sofia, Zdrave Street 2, 1431 Sofia, Bulgaria;
| | - Alexandar M. Zhivkov
- Scientific Research Center, “St. Kliment Ohridski” Sofia University, 8 Dragan Tsankov Blvd., 1164 Sofia, Bulgaria
| |
Collapse
|
7
|
Raisinghani N, Alshahrani M, Gupta G, Xiao S, Tao P, Verkhivker G. AlphaFold2-Enabled Atomistic Modeling of Epistatic Binding Mechanisms for the SARS-CoV-2 Spike Omicron XBB.1.5, EG.5 and FLip Variants: Convergent Evolution Hotspots Cooperate to Control Stability and Conformational Adaptability in Balancing ACE2 Binding and Antibody Resistance. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.12.11.571185. [PMID: 38168257 PMCID: PMC10760024 DOI: 10.1101/2023.12.11.571185] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/05/2024]
Abstract
In this study, we combined AI-based atomistic structural modeling and microsecond molecular simulations of the SARS-CoV-2 Spike complexes with the host receptor ACE2 for XBB.1.5+L455F, XBB.1.5+F456L(EG.5) and XBB.1.5+L455F/F456L (FLip) lineages to examine the mechanisms underlying the role of convergent evolution hotspots in balancing ACE2 binding and antibody evasion. Using the ensemble-based mutational scanning of the spike protein residues and physics-based rigorous computations of binding affinities, we identified binding energy hotspots and characterized molecular basis underlying epistatic couplings between convergent mutational hotspots. Consistent with the experiments, the results revealed the mediating role of Q493 hotspot in synchronization of epistatic couplings between L455F and F456L mutations providing a quantitative insight into the mechanism underlying differences between XBB lineages. Mutational profiling is combined with network-based model of epistatic couplings showing that the Q493, L455 and F456 sites mediate stable communities at the binding interface with ACE2 and can serve as stable mediators of non-additive couplings. Structure-based mutational analysis of Spike protein binding with the class 1 antibodies quantified the critical role of F456L and F486P mutations in eliciting strong immune evasion response. The results of this analysis support a mechanism in which the emergence of EG.5 and FLip variants may have been dictated by leveraging strong epistatic effects between several convergent revolutionary hotspots that provide synergy between the improved ACE2 binding and broad neutralization resistance. This interpretation is consistent with the notion that functionally balanced substitutions which simultaneously optimize immune evasion and high ACE2 affinity may continue to emerge through lineages with beneficial pair or triplet combinations of RBD mutations involving mediators of epistatic couplings and sites in highly adaptable RBD regions.
Collapse
|
8
|
Raisinghani N, Alshahrani M, Gupta G, Xiao S, Tao P, Verkhivker G. Accurate Characterization of Conformational Ensembles and Binding Mechanisms of the SARS-CoV-2 Omicron BA.2 and BA.2.86 Spike Protein with the Host Receptor and Distinct Classes of Antibodies Using AlphaFold2-Augmented Integrative Computational Modeling. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.11.18.567697. [PMID: 38045395 PMCID: PMC10690158 DOI: 10.1101/2023.11.18.567697] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/05/2023]
Abstract
The latest wave SARS-CoV-2 Omicron variants displayed a growth advantage and the increased viral fitness through convergent evolution of functional hotspots that work synchronously to balance fitness requirements for productive receptor binding and efficient immune evasion. In this study, we combined AlphaFold2-based structural modeling approaches with all-atom MD simulations and mutational profiling of binding energetics and stability for prediction and comprehensive analysis of the structure, dynamics, and binding of the SARS-CoV-2 Omicron BA.2.86 spike variant with ACE2 host receptor and distinct classes of antibodies. We adapted several AlphaFold2 approaches to predict both structure and conformational ensembles of the Omicron BA.2.86 spike protein in the complex with the host receptor. The results showed that AlphaFold2-predicted conformational ensemble of the BA.2.86 spike protein complex can accurately capture the main dynamics signatures obtained from microscond molecular dynamics simulations. The ensemble-based dynamic mutational scanning of the receptor binding domain residues in the BA.2 and BA.2.86 spike complexes with ACE2 dissected the role of the BA.2 and BA.2.86 backgrounds in modulating binding free energy changes revealing a group of conserved hydrophobic hotspots and critical variant-specific contributions of the BA.2.86 mutational sites R403K, F486P and R493Q. To examine immune evasion properties of BA.2.86 in atomistic detail, we performed large scale structure-based mutational profiling of the S protein binding interfaces with distinct classes of antibodies that displayed significantly reduced neutralization against BA.2.86 variant. The results quantified specific function of the BA.2.86 mutations to ensure broad resistance against different classes of RBD antibodies. This study revealed the molecular basis of compensatory functional effects of the binding hotspots, showing that BA.2.86 lineage may have primarily evolved to improve immune escape while modulating binding affinity with ACE2 through cooperative effect of R403K, F486P and R493Q mutations. The study supports a hypothesis that the impact of the increased ACE2 binding affinity on viral fitness is more universal and is mediated through cross-talk between convergent mutational hotspots, while the effect of immune evasion could be more variant-dependent.
Collapse
|
9
|
Zhang S, Wang N, Zhang Q, Guan R, Qu Z, Sun L, Li J. The Rise of Electroactive Materials in Face Masks for Preventing Virus Infections. ACS APPLIED MATERIALS & INTERFACES 2023; 15:48839-48854. [PMID: 37815875 DOI: 10.1021/acsami.3c10465] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/12/2023]
Abstract
Air-transmitted pathogens may cause severe epidemics, posing considerable threats to public health and safety. Wearing a face mask is one of the most effective ways to prevent respiratory virus infection transmission. Especially since the new coronavirus pandemic, electroactive materials have received much attention in antiviral face masks due to their highly efficient antiviral capabilities, flexible structural design, excellent sustainability, and outstanding safety. This review first introduces the mechanism for preventing viral infection or the inactivation of viruses by electroactive materials. Then, the applications of electrostatic-, conductive-, triboelectric-, and microbattery-based materials in face masks are described in detail. Finally, the problems of various electroactive antiviral materials are summarized, and the prospects for their future development directions are discussed. In conclusion, electroactive materials have attracted great attention for antiviral face masks, and this review will provide a reference for materials scientists and engineers in antiviral materials and interfaces.
Collapse
Affiliation(s)
- Shaohua Zhang
- College of Textiles and Clothing, the Affiliated Hospital of Qingdao University, Qingdao University, Qingdao 266071, People's Republic of China
| | - Na Wang
- College of Textiles and Clothing, the Affiliated Hospital of Qingdao University, Qingdao University, Qingdao 266071, People's Republic of China
- Industrial Research Institute of Nonwovens and Technical Textiles, Shandong Center for Engineered Nonwovens, Qingdao 266071, People's Republic of China
| | - Qian Zhang
- Department of Respirology, Qingdao Women and Children's Hospital, Qingdao 266034, People's Republic of China
| | - Renzheng Guan
- College of Textiles and Clothing, the Affiliated Hospital of Qingdao University, Qingdao University, Qingdao 266071, People's Republic of China
| | - Zhenghai Qu
- College of Textiles and Clothing, the Affiliated Hospital of Qingdao University, Qingdao University, Qingdao 266071, People's Republic of China
| | - Lirong Sun
- College of Textiles and Clothing, the Affiliated Hospital of Qingdao University, Qingdao University, Qingdao 266071, People's Republic of China
| | - Jiwei Li
- College of Textiles and Clothing, the Affiliated Hospital of Qingdao University, Qingdao University, Qingdao 266071, People's Republic of China
- Industrial Research Institute of Nonwovens and Technical Textiles, Shandong Center for Engineered Nonwovens, Qingdao 266071, People's Republic of China
| |
Collapse
|
10
|
Alshahrani M, Gupta G, Xiao S, Tao P, Verkhivker G. Comparative Analysis of Conformational Dynamics and Systematic Characterization of Cryptic Pockets in the SARS-CoV-2 Omicron BA.2, BA.2.75 and XBB.1 Spike Complexes with the ACE2 Host Receptor: Confluence of Binding and Structural Plasticity in Mediating Networks of Conserved Allosteric Sites. Viruses 2023; 15:2073. [PMID: 37896850 PMCID: PMC10612107 DOI: 10.3390/v15102073] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2023] [Revised: 10/04/2023] [Accepted: 10/06/2023] [Indexed: 10/29/2023] Open
Abstract
In the current study, we explore coarse-grained simulations and atomistic molecular dynamics together with binding energetics scanning and cryptic pocket detection in a comparative examination of conformational landscapes and systematic characterization of allosteric binding sites in the SARS-CoV-2 Omicron BA.2, BA.2.75 and XBB.1 spike full-length trimer complexes with the host receptor ACE2. Microsecond simulations, Markov state models and mutational scanning of binding energies of the SARS-CoV-2 BA.2 and BA.2.75 receptor binding domain complexes revealed the increased thermodynamic stabilization of the BA.2.75 variant and significant dynamic differences between these Omicron variants. Molecular simulations of the SARS-CoV-2 Omicron spike full-length trimer complexes with the ACE2 receptor complemented atomistic studies and enabled an in-depth analysis of mutational and binding effects on conformational dynamic and functional adaptability of the Omicron variants. Despite considerable structural similarities, Omicron variants BA.2, BA.2.75 and XBB.1 can induce unique conformational dynamic signatures and specific distributions of the conformational states. Using conformational ensembles of the SARS-CoV-2 Omicron spike trimer complexes with ACE2, we conducted a comprehensive cryptic pocket screening to examine the role of Omicron mutations and ACE2 binding on the distribution and functional mechanisms of the emerging allosteric binding sites. This analysis captured all experimentally known allosteric sites and discovered networks of inter-connected and functionally relevant allosteric sites that are governed by variant-sensitive conformational adaptability of the SARS-CoV-2 spike structures. The results detailed how ACE2 binding and Omicron mutations in the BA.2, BA.2.75 and XBB.1 spike complexes modulate the distribution of conserved and druggable allosteric pockets harboring functionally important regions. The results are significant for understanding the functional roles of druggable cryptic pockets that can be used for allostery-mediated therapeutic intervention targeting conformational states of the Omicron variants.
Collapse
Affiliation(s)
- Mohammed Alshahrani
- 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; (M.A.); (G.G.)
| | - Grace Gupta
- 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; (M.A.); (G.G.)
| | - Sian Xiao
- Department of Chemistry, Center for Research Computing, Center for Drug Discovery, Design, and Delivery (CD4), Southern Methodist University, Dallas, TX 75275, USA; (S.X.); (P.T.)
| | - Peng Tao
- Department of Chemistry, Center for Research Computing, Center for Drug Discovery, Design, and Delivery (CD4), Southern Methodist University, Dallas, TX 75275, USA; (S.X.); (P.T.)
| | - 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; (M.A.); (G.G.)
- Department of Biomedical and Pharmaceutical Sciences, Chapman University School of Pharmacy, Irvine, CA 92618, USA
| |
Collapse
|
11
|
Hristova SH, Zhivkov AM. Omicron Coronavirus: pH-Dependent Electrostatic Potential and Energy of Association of Spike Protein to ACE2 Receptor. Viruses 2023; 15:1752. [PMID: 37632094 PMCID: PMC10460073 DOI: 10.3390/v15081752] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2023] [Revised: 08/12/2023] [Accepted: 08/14/2023] [Indexed: 08/27/2023] Open
Abstract
The association of the S-protein of the SARS-CoV-2 beta coronavirus to ACE2 receptors of the human epithelial cells determines its contagiousness and pathogenicity. We computed the pH-dependent electric potential on the surface of the interacting globular proteins and pH-dependent Gibbs free energy at the association of the wild-type strain and the omicron variant. The calculated isoelectric points of the ACE2 receptor (pI 5.4) and the S-protein in trimeric form (pI 7.3, wild type), (pI 7.8, omicron variant), experimentally verified by isoelectric focusing, show that at pH 6-7, the S1-ACE2 association is conditioned by electrostatic attraction of the oppositely charged receptor and viral protein. The comparison of the local electrostatic potentials of the omicron variant and the wild-type strain shows that the point mutations alter the electrostatic potential in a relatively small area on the surface of the receptor-binding domain (RBD) of the S1 subunit. The appearance of seven charge-changing point mutations in RBD (equivalent to three additional positive charges) leads to a stronger S1-ACE2 association at pH 5.5 (typical for the respiratory tract) and a weaker one at pH 7.4 (characteristic of the blood plasma); this reveals the reason for the higher contagiousness but lower pathogenicity of the omicron variant in comparison to the wild-type strain.
Collapse
Affiliation(s)
- Svetlana H. Hristova
- Department of Medical Physics and Biophysics, Medical Faculty, Medical University—Sofia, Zdrave Str. 2, 1431 Sofia, Bulgaria;
| | - Alexandar M. Zhivkov
- Institute of Physical Chemistry, Bulgarian Academy of Sciences, Acad. G. Bonchev Str., bl. 11, 1113 Sofia, Bulgaria
| |
Collapse
|
12
|
Božič A, Podgornik R. Evolutionary changes in the number of dissociable amino acids on spike proteins and nucleoproteins of SARS-CoV-2 variants. Virus Evol 2023; 9:vead040. [PMID: 37583936 PMCID: PMC10424713 DOI: 10.1093/ve/vead040] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2023] [Revised: 06/12/2023] [Accepted: 06/27/2023] [Indexed: 08/17/2023] Open
Abstract
The spike protein of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for target recognition, cellular entry, and endosomal escape of the virus. At the same time, it is the part of the virus that exhibits the greatest sequence variation across the many variants which have emerged during its evolution. Recent studies have indicated that with progressive lineage emergence, the positive charge on the spike protein has been increasing, with certain positively charged amino acids (AAs) improving the binding of the spike protein to cell receptors. We have performed a detailed analysis of dissociable AAs of more than 1400 different SARS-CoV-2 lineages, which confirms these observations while suggesting that this progression has reached a plateau with Omicron and its subvariants and that the positive charge is not increasing further. Analysis of the nucleocapsid protein shows no similar increase in positive charge with novel variants, which further indicates that positive charge of the spike protein is being evolutionarily selected for. Furthermore, comparison with the spike proteins of known coronaviruses shows that already the wild-type SARS-CoV-2 spike protein carries an unusually large amount of positively charged AAs when compared to most other betacoronaviruses. Our study sheds light on the evolutionary changes in the number of dissociable AAs on the spike protein of SARS-CoV-2, complementing existing studies and providing a stepping stone towards a better understanding of the relationship between the spike protein charge and viral infectivity and transmissibility.
Collapse
Affiliation(s)
- Anže Božič
- Department of Theoretical Physics, Jožef Stefan Institute, Jamova 39, Ljubljana SI-1000, Slovenia
| | - Rudolf Podgornik
- School of Physical Sciences, University of Chinese Academy of Sciences, No. 19A Yuquan Road, Shijingshan District, Beijing 100049, China
- Kavli Institute for Theoretical Sciences, University of Chinese Academy of Sciences, No. 3 Nanyitiao, Zhongguancun, Haidian District, Beijing 100049, China
- CAS Key Laboratory of Soft Matter Physics, Institute of Physics, Chinese Academy of Sciences, No. 8 3rd South Street, Zhongguancun, Haidian District, Beijing 100190, China
- Wenzhou Institute of the University of Chinese Academy of Sciences, No.1 Jinlian Road, Wenzhou, Zhejiang 325001, China
- Department of Physics, Faculty of Mathematics and Physics, University of Ljubljana, Jadranska 19, Ljubljana SI-1000, Slovenia
| |
Collapse
|
13
|
Lauster D, Osterrieder K, Haag R, Ballauff M, Herrmann A. Respiratory viruses interacting with cells: the importance of electrostatics. Front Microbiol 2023; 14:1169547. [PMID: 37440888 PMCID: PMC10333706 DOI: 10.3389/fmicb.2023.1169547] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2023] [Accepted: 06/08/2023] [Indexed: 07/15/2023] Open
Abstract
The COVID-19 pandemic has rekindled interest in the molecular mechanisms involved in the early steps of infection of cells by viruses. Compared to SARS-CoV-1 which only caused a relatively small albeit deadly outbreak, SARS-CoV-2 has led to fulminant spread and a full-scale pandemic characterized by efficient virus transmission worldwide within a very short time. Moreover, the mutations the virus acquired over the many months of virus transmission, particularly those seen in the Omicron variant, have turned out to result in an even more transmissible virus. Here, we focus on the early events of virus infection of cells. We review evidence that the first decisive step in this process is the electrostatic interaction of the spike protein with heparan sulfate chains present on the surface of target cells: Patches of cationic amino acids located on the surface of the spike protein can interact intimately with the negatively charged heparan sulfate chains, which results in the binding of the virion to the cell surface. In a second step, the specific interaction of the receptor binding domain (RBD) within the spike with the angiotensin-converting enzyme 2 (ACE2) receptor leads to the uptake of bound virions into the cell. We show that these events can be expressed as a semi-quantitative model by calculating the surface potential of different spike proteins using the Adaptive Poison-Boltzmann-Solver (APBS). This software allows visualization of the positive surface potential caused by the cationic patches, which increased markedly from the original Wuhan strain of SARS-CoV-2 to the Omicron variant. The surface potential thus enhanced leads to a much stronger binding of the Omicron variant as compared to the original wild-type virus. At the same time, data taken from the literature demonstrate that the interaction of the RBD of the spike protein with the ACE2 receptor remains constant within the limits of error. Finally, we briefly digress to other viruses and show the usefulness of these electrostatic processes and calculations for cell-virus interactions more generally.
Collapse
Affiliation(s)
- Daniel Lauster
- Institut für Pharmazie, Biopharmazeutika, Freie Universität Berlin, Berlin, Germany
| | | | - Rainer Haag
- Institut für Chemie und Biochemie, SupraFAB, Freie Universität Berlin, Berlin, Germany
| | - Matthias Ballauff
- Institut für Chemie und Biochemie, SupraFAB, Freie Universität Berlin, Berlin, Germany
| | - Andreas Herrmann
- Institut für Chemie und Biochemie, SupraFAB, Freie Universität Berlin, Berlin, Germany
| |
Collapse
|
14
|
Ciccozzi M, Pascarella S. Two sides of the same coin: the N-terminal and the receptor binding domains of SARS-CoV-2 Spike. Future Virol 2023:10.2217/fvl-2022-0181. [PMID: 36896145 PMCID: PMC9987531 DOI: 10.2217/fvl-2022-0181] [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: 10/03/2022] [Accepted: 01/31/2023] [Indexed: 03/08/2023]
Abstract
The SARS-CoV-2 Spike receptor binding domain and N-terminal domain interact with each other in an intricate mechanism. Mutations modulate the interplay between the Spike and host molecules. This editorial comments on the intricacies of SARS-CoV-2 Spike interactions.
Collapse
Affiliation(s)
- Massimo Ciccozzi
- Medical Statistic & Molecular Epidemiology Unit, University of Biomedical Campus, Rome, Italy
| | - Stefano Pascarella
- Department of Biochemical Sciences 'A Rossi Fanelli', Sapienza Università di Roma, Rome, 00185, Italy
| |
Collapse
|
15
|
Neamtu A, Mocci F, Laaksonen A, Barroso da Silva FL. Towards an optimal monoclonal antibody with higher binding affinity to the receptor-binding domain of SARS-CoV-2 spike proteins from different variants. Colloids Surf B Biointerfaces 2023; 221:112986. [PMID: 36375294 PMCID: PMC9617679 DOI: 10.1016/j.colsurfb.2022.112986] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2022] [Revised: 09/13/2022] [Accepted: 10/27/2022] [Indexed: 11/13/2022]
Abstract
A highly efficient and robust multiple scales in silico protocol, consisting of atomistic Molecular Dynamics (MD), coarse-grain (CG) MD, and constant-pH CG Monte Carlo (MC), has been developed and used to study the binding affinities of selected antigen-binding fragments of the monoclonal antibody (mAbs) CR3022 and several of its here optimized versions against 11 SARS-CoV-2 variants including the wild type. Totally 235,000 mAbs structures were initially generated using the RosettaAntibodyDesign software, resulting in top 10 scored CR3022-like-RBD complexes with critical mutations and compared to the native one, all having the potential to block virus-host cell interaction. Of these 10 finalists, two candidates were further identified in the CG simulations to be the best against all SARS-CoV-2 variants. Surprisingly, all 10 candidates and the native CR3022 exhibited a higher affinity for the Omicron variant despite its highest number of mutations. The multiscale protocol gives us a powerful rational tool to design efficient mAbs. The electrostatic interactions play a crucial role and appear to be controlling the affinity and complex building. Studied mAbs carrying a more negative total net charge show a higher affinity. Structural determinants could be identified in atomistic simulations and their roles are discussed in detail to further hint at a strategy for designing the best RBD binder. Although the SARS-CoV-2 was specifically targeted in this work, our approach is generally suitable for many diseases and viral and bacterial pathogens, leukemia, cancer, multiple sclerosis, rheumatoid, arthritis, lupus, and more.
Collapse
Affiliation(s)
- Andrei Neamtu
- Department of Physiology, "Grigore T. Popa" University of Medicine and Pharmacy of Iasi, Str. Universitatii nr. 16, 700051 Iasi, România; TRANSCEND Centre - Regional Institute of Oncology (IRO) Iasi, Str. General Henri Mathias Berthelot, Nr. 2-4 Iași, România
| | - Francesca Mocci
- University of Cagliari, Department of Chemical and Geological Sciences, Campus Monserrato, SS 554 bivio per Sestu, 09042 Monserrato, Italy
| | - Aatto Laaksonen
- Centre of Advanced Research in Bionanoconjugates and Biopolymers, PetruPoni Institute of Macromolecular Chemistry Aleea Grigore Ghica-Voda, 41 A, 700487 Iasi, Romania; University of Cagliari, Department of Chemical and Geological Sciences, Campus Monserrato, SS 554 bivio per Sestu, 09042 Monserrato, Italy; Department of Materials and Environmental Chemistry, Arrhenius Laboratory, Stockholm University, SE-106 91 Stockholm, Sweden; State Key Laboratory of Materials-Oriented and Chemical Engineering, Nanjing Tech University, Nanjing 210009, PR China; Department of Engineering Sciences and Mathematics, Division of Energy Science, Luleå University of Technology, SE-97187 Luleå, Sweden
| | - Fernando L Barroso da Silva
- Universidade de São Paulo, Departamento de Ciências Biomoleculares, Faculdade de Ciências Farmacêuticas de Ribeirão Preto, Av. café, s/no - campus da USP, BR-14040-903 Ribeirão Preto, SP, Brazil; Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, NC 27695, USA.
| |
Collapse
|
16
|
Lee LYY, Suryadinata R, McCafferty C, Ignjatovic V, Purcell DFJ, Robinson P, Morton CJ, Parker MW, Anderson GP, Monagle P, Subbarao K, Neil JA. Heparin Inhibits SARS-CoV-2 Replication in Human Nasal Epithelial Cells. Viruses 2022; 14:v14122620. [PMID: 36560624 PMCID: PMC9785945 DOI: 10.3390/v14122620] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2022] [Revised: 11/18/2022] [Accepted: 11/18/2022] [Indexed: 11/27/2022] Open
Abstract
SARS-CoV-2 is the causative agent of the COVID-19 pandemic. Vaccination, supported by social and public health measures, has proven efficacious for reducing disease severity and virus spread. However, the emergence of highly transmissible viral variants that escape prior immunity highlights the need for additional mitigation approaches. Heparin binds the SARS-CoV-2 spike protein and can inhibit virus entry and replication in susceptible human cell lines and bronchial epithelial cells. Primary infection predominantly occurs via the nasal epithelium, but the nasal cell biology of SARS-CoV-2 is not well studied. We hypothesized that prophylactic intranasal administration of heparin may provide strain-agnostic protection for household contacts or those in high-risk settings against SARS-CoV-2 infection. Therefore, we investigated the ability of heparin to inhibit SARS-CoV-2 infection and replication in differentiated human nasal epithelial cells and showed that prolonged exposure to heparin inhibits virus infection. Furthermore, we establish a method for PCR detection of SARS-CoV-2 viral genomes in heparin-treated samples that can be adapted for the detection of viruses in clinical studies.
Collapse
Affiliation(s)
- Leo Yi Yang Lee
- Department of Microbiology and Immunology, The University of Melbourne at the Peter Doherty Institute for Infection and Immunity, Melbourne, VIC 3000, Australia
| | - Randy Suryadinata
- Department of Respiratory Medicine, Royal Children’s Hospital, Parkville, VIC 3052, Australia
- Infection and Immunity, Murdoch Children’s Research Institute, Royal Children’s Hospital, Parkville, VIC 3052, Australia
| | - Conor McCafferty
- Department of Paediatrics, The University of Melbourne, Parkville, VIC 3052, Australia
- Haematology, Murdoch Children’s Research Institute, Parkville, VIC 3052, Australia
| | - Vera Ignjatovic
- Department of Paediatrics, The University of Melbourne, Parkville, VIC 3052, Australia
| | - Damian F. J. Purcell
- Department of Microbiology and Immunology, The University of Melbourne at the Peter Doherty Institute for Infection and Immunity, Melbourne, VIC 3000, Australia
| | - Phil Robinson
- Department of Respiratory Medicine, Royal Children’s Hospital, Parkville, VIC 3052, Australia
- Infection and Immunity, Murdoch Children’s Research Institute, Royal Children’s Hospital, Parkville, VIC 3052, Australia
- Department of Paediatrics, The University of Melbourne, Parkville, VIC 3052, Australia
| | - Craig J. Morton
- Department of Biochemistry and Pharmacology, Bio21 Molecular Science and Biotechnology Institute, The University of Melbourne, Parkville, VIC 3010, Australia
| | - Michael W. Parker
- Department of Biochemistry and Pharmacology, Bio21 Molecular Science and Biotechnology Institute, The University of Melbourne, Parkville, VIC 3010, Australia
- St. Vincent’s Institute of Medical Research, Fitzroy, VIC 3065, Australia
| | - Gary P. Anderson
- Lung Health Research Centre, Department of Biochemistry and Pharmacology, The University of Melbourne, Parkville, VIC 3010, Australia
| | - Paul Monagle
- Department of Paediatrics, The University of Melbourne, Parkville, VIC 3052, Australia
- Haematology, Murdoch Children’s Research Institute, Parkville, VIC 3052, Australia
- Department of Haematology, Royal Children’s Hospital, Parkville, VIC 3052, Australia
- Kids Cancer Centre, Sydney Children’s Hospital, Randwick, NSW 2031, Australia
| | - Kanta Subbarao
- Department of Microbiology and Immunology, The University of Melbourne at the Peter Doherty Institute for Infection and Immunity, Melbourne, VIC 3000, Australia
- WHO Collaborating Centre for Reference and Research on Influenza, Melbourne, VIC 3000, Australia
- Correspondence:
| | - Jessica A. Neil
- Department of Microbiology and Immunology, The University of Melbourne at the Peter Doherty Institute for Infection and Immunity, Melbourne, VIC 3000, Australia
| |
Collapse
|