Supramolecular Organization of SARS-CoV and SARS-CoV-2 Virions Revealed by Coarse-Grained Models of Intact Virus Envelopes

Modeling membranes in situ

Molecular dynamics simulations of cellular membranes have come a long way-from simple model lipid bilayers to multicomponent systems capturing the crowded and complex nature of real cell membranes. In this opinionated minireview, we discuss the current challenge to simulate the dynamics of membranes in their native environment, in situ, with the prospect of reaching the level of whole cells and cell organelles using an integrative modeling framework.

Adsorption-Driven Deformation and Footprints of the RBD Proteins in SARS-CoV-2 Variants on Biological and Inanimate Surfaces

Respiratory viruses, carried through airborne microdroplets, frequently adhere to surfaces, including plastics and metals. However, our understanding of the interactions between viruses and materials remains limited, particularly in scenarios involving polarizable surfaces. Here, we investigate the role of the receptor-binding domain (RBD) of the spike protein mutations on the adsorption of SARS-CoV-2 to hydrophobic and hydrophilic surfaces employing molecular simulations. To contextualize our findings, we contrast the interactions on inanimate surfaces with those on native biological interfaces, specifically the angiotensin-converting enzyme 2. Notably, we identify a 2-fold increase in structural deformations for the protein's receptor binding motif (RBM) onto inanimate surfaces, indicative of enhanced shock-absorbing mechanisms. Furthermore, the distribution of adsorbed amino acids (landing footprints) on the inanimate surface reveals a distinct regional asymmetry relative to the biological interface, with roughly half of the adsorbed amino acids arranged in opposite sites. In spite of the H-bonds formed at the hydrophilic substrate, the simulations consistently show a higher number of contacts and interfacial area with the hydrophobic surface, where the wild-type RBD adsorbs more strongly than the Delta or Omicron RBDs. In contrast, the adsorption of Delta and Omicron to hydrophilic surfaces was characterized by a distinctive hopping-pattern. The novel shock-absorbing mechanisms identified in the virus adsorption on inanimate surfaces show the embedded high-deformation capacity of the RBD without losing its secondary structure, which could lead to current experimental strategies in the design of virucidal surfaces.

Computational and experimental validation of phthalocyanine and hypericin as effective SARS-CoV-2 fusion inhibitors

Phthalocyanine and hypericin have been previously identified as possible SARS-CoV-2 Spike glycoprotein fusion inhibitors through a virtual screening procedure. In this paper, atomistic simulations of metal-free phthalocyanines and atomistic and coarse-grained simulations of hypericins, placed around a complete model of the Spike embedded in a viral membrane, allowed to further explore their multi-target inhibitory potential, uncovering their binding to key protein functional regions and their propensity to insert in the membrane. Following computational results, pre-treatment of a pseudovirus expressing the SARS-CoV-2 Spike protein with low compounds concentrations resulted in a strong inhibition of its entry into cells, suggesting the activity of these molecules should involve the direct targeting of the viral envelope surface. The combination of computational and in vitro results hence supports the role of hypericin and phthalocyanine as promising SARS-CoV-2 entry inhibitors, further endorsed by literature reporting the efficacy of these compounds in inhibiting SARS-CoV-2 activity and in treating hospitalized COVID-19 patients.Communicated by Ramaswamy H. Sarma.

Unraveling the binding mechanisms of SARS-CoV-2 variants through molecular simulations

The emergence of SARS-CoV-2 variants like Delta (AY.29) and Omicron (EG.5) poses continued challenges for vaccines and therapeutics. Mutations in the viral spike protein are key in altering infectivity and immune evasion. This study uses computational modeling to investigate the molecular binding mechanisms between spike protein variants and the ACE2 host receptor. Using the MARTNI force field, coarse-grained molecular dynamics (CGMD) simulations and nudged elastic band (NEB) calculations explore spike-ACE2 interactions for the wild type, Delta variant, and Omicron variant. The simulations reveal Omicron has the strongest binding affinity at -128.35 ± 10.91 kcal/mol, followed by Delta and wild type. Key mutations in Delta and Omicron, like Q493R and Q498R, optimize electrostatic contacts, enhancing ACE2 interactions. The wild-type spike has the highest transition state energy barrier at 17.87 kcal/mol, while Delta has the lowest barrier at 9.21 kcal/mol. Despite slightly higher dual barriers, Omicron's increased binding energy lowers its overall barrier to rapidly bind ACE2. These findings provide residue-level insights into mutation effects on SARS-CoV-2 infectivity. The computational modeling elucidates mechanisms underlying spike-ACE2 binding kinetics, aiding the development of vaccines and therapies targeting emerging viral strains.

Mechanisms of allostery at the viral surface through the eyes of molecular simulation

The outermost surface layer of any virus is formed by either a capsid shell or envelope. Such layers have traditionally been thought of as immovable structures, but it is becoming apparent that they cannot be viewed exclusively as static architectures protecting the viral genome. A limited number of proteins on the virion surface must perform a multitude of functions in order to orchestrate the viral life cycle, and allostery can regulate their structures at multiple levels of organization, spanning individual molecules, protomers, large oligomeric assemblies, or entire viral surfaces. Here, we review recent contributions from the molecular simulation field to viral surface allostery, with a particular focus on the trimeric spike glycoprotein emerging from the coronavirus surface, and the icosahedral flaviviral envelope complex. As emerging viral pathogens continue to pose a global threat, an improved understanding of viral dynamics and allosteric regulation will prove crucial in developing novel therapeutic strategies.

Pragmatic Coarse-Graining of Proteins: Models and Applications

The molecular details involved in the folding, dynamics, organization, and interaction of proteins with other molecules are often difficult to assess by experimental techniques. Consequently, computational models play an ever-increasing role in the field. However, biological processes involving large-scale protein assemblies or long time scale dynamics are still computationally expensive to study in atomistic detail. For these applications, employing coarse-grained (CG) modeling approaches has become a key strategy. In this Review, we provide an overview of what we call pragmatic CG protein models, which are strategies combining, at least in part, a physics-based implementation and a top-down experimental approach to their parametrization. In particular, we focus on CG models in which most protein residues are represented by at least two beads, allowing these models to retain some degree of chemical specificity. A description of the main modern pragmatic protein CG models is provided, including a review of the most recent applications and an outlook on future perspectives in the field.

MicroRNAs Derived from Extracellular Vesicles: Keys to Understanding SARS-CoV-2 Vaccination Response in Cancer Patients?

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) provoked a global pandemic identified as coronavirus disease (COVID-19), with millions of deaths worldwide. However, several important questions regarding its impact on public health remain unanswered, such as the impact of vaccination on vulnerable subpopulations such as cancer patients. Cytokine storm and a sustained inflammatory state are commonly associated with immune cell depletion, being manifested in most immunocompromised individuals. This strong immunosuppression can lead to a dysfunctional antiviral response to natural viral infection and compromised vaccination response. Extracellular vesicles (EVs) are membrane-bound vesicles released from cells that are involved in intercellular communication. EVs carry various molecules including microRNAs that play a crucial role in COVID-19 pathophysiology, influencing cellular responses. This review summarizes the state of the art concerning the role of EV-derived miRNAs in COVID-19 infection and their potential use as prognosis biomarkers for vaccination response in cancer patients.

Understanding Virus Structure and Dynamics through Molecular Simulations

Viral outbreaks remain a serious threat to human and animal populations and motivate the continued development of antiviral drugs and vaccines, which in turn benefits from a detailed understanding of both viral structure and dynamics. While great strides have been made in characterizing these systems experimentally, molecular simulations have proven to be an essential, complementary approach. In this work, we review the contributions of molecular simulations to the understanding of viral structure, functional dynamics, and processes related to the viral life cycle. Approaches ranging from coarse-grained to all-atom representations are discussed, including current efforts at modeling complete viral systems. Overall, this review demonstrates that computational virology plays an essential role in understanding these systems.

Short Peptide Self-Assembly in the Martini Coarse-Grain Force Field Family

ConspectusPivotal to the success of any computational experiment is the ability to make reliable predictions about the system under study and the time required to yield these results. Biomolecular interactions is one area of research that sits in every camp of resolution vs the time required, from the quantum mechanical level to in vivo studies. At an approximate midpoint, there is coarse-grained molecular dynamics, for which the Martini force fields have become the most widely used, fast enough to simulate the entire membrane of a mitochondrion though lacking atom-specific precision. While many force fields have been parametrized to account for a specific system under study, the Martini force field has aimed at casting a wider net with more generalized bead types that have demonstrated suitability for broad use and reuse in applications from protein-graphene oxide coassembly to polysaccharides interactions.In this Account, the progressive (Martini versions 1 through 3) and peripheral (Sour Martini, constant pH, Martini Straight, Dry Martini, etc.) developmental trajectory of the Martini force field will be analyzed in terms of self-assembling systems with a focus on short (two to three amino acids) peptide self-assembly in aqueous environments. In particular, this will focus on the effects of the Martini solvent model and compare how changes in bead definitions and mapping have effects on different systems. Considerable effort in the development of Martini has been expended to reduce the "stickiness" of amino acids to better simulate proteins in bilayers. We have included in this Account a short study of dipeptide self-assembly in water, using all mainstream Martini force fields, to examine their ability to reproduce this behavior. The three most recently released versions of Martini and variations in their solvents are used to simulate in triplicate all 400 dipeptides of the 20 gene-encoded amino acids. The ability of the force fields to model the self-assembly of the dipeptides in aqueoues environments is determined by the measurement of the aggregation propensity, and additional descriptors are used to gain further insight into the dipeptide aggregates.

Molecular architecture and dynamics of SARS-CoV-2 envelope by integrative modeling

Despite tremendous efforts, the exact structure of SARS-CoV-2 and related betacoronaviruses remains elusive. SARS-CoV-2 envelope is a key structural component of the virion that encapsulates viral RNA. It is composed of three structural proteins, spike, membrane (M), and envelope, which interact with each other and with the lipids acquired from the host membranes. Here, we developed and applied an integrative multi-scale computational approach to model the envelope structure of SARS-CoV-2 with near atomistic detail, focusing on studying the dynamic nature and molecular interactions of its most abundant, but largely understudied, M protein. The molecular dynamics simulations allowed us to test the envelope stability under different configurations and revealed that the M dimers agglomerated into large, filament-like, macromolecular assemblies with distinct molecular patterns. These results are in good agreement with current experimental data, demonstrating a generic and versatile approach to model the structure of a virus de novo.

SARS-CoV-2 Omicron Subvariants Balance Host Cell Membrane, Receptor, and Antibody Docking via an Overlapping Target Site

Variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are emerging rapidly and offer surfaces that are optimized for recognition of host cell membranes while also evading antibodies arising from vaccinations and previous infections. Host cell infection is a multi-step process in which spike heads engage lipid bilayers and one or more angiotensin-converting enzyme 2 (ACE-2) receptors. Here, the membrane binding surfaces of Omicron subvariants are compared using cryo-electron microscopy (cEM) structures of spike trimers from BA.2, BA.2.12.1, BA.2.13, BA.2.75, BA.3, BA.4, and BA.5 viruses. Despite significant differences around mutated sites, they all maintain strong membrane binding propensities that first appeared in BA.1. Both their closed and open states retain elevated membrane docking capacities, although the presence of more closed than open states diminishes opportunities to bind receptors while enhancing membrane engagement. The electrostatic dipoles are generally conserved. However, the BA.2.75 spike dipole is compromised, and its ACE-2 affinity is increased, and BA.3 exhibits the opposite pattern. We propose that balancing the functional imperatives of a stable, readily cleavable spike that engages both lipid bilayers and receptors while avoiding host defenses underlies betacoronavirus evolution. This provides predictive criteria for rationalizing future pandemic waves and COVID-19 transmissibility while illuminating critical sites and strategies for simultaneously combating multiple variants.

Role of cholesterol-recognition motifs in the infectivity of SARS-CoV-2 variants

The presence of linear amino acid motifs with the capacity to recognize the neutral lipid cholesterol, known as Cholesterol Recognition/interaction Amino acid Consensus sequence (CRAC), and its inverse or mirror image, CARC, has recently been reported in the primary sequence of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike S homotrimeric glycoprotein. These motifs also occur in the two other pathogenic coronaviruses, SARS-CoV, and Middle-East respiratory syndrome CoV (MERS-CoV), most conspicuously in the transmembrane domain, the fusion peptide, the amino-terminal domain, and the receptor binding domain of SARS-CoV-2 S protein. Here we analyze the presence of cholesterol-recognition motifs in these key regions of the spike glycoprotein in the pathogenic CoVs. We disclose the inherent pathophysiological implications of the cholesterol motifs in the virus-host cell interactions and variant infectivity.

Effect of Disulfide Bridge on the Binding of SARS-CoV-2 Fusion Peptide to Cell Membrane: A Coarse-Grained Study

In this paper, we present the parameterization of the CAVS coarse-grained (CG) force field for 20 amino acids, and our CG simulations show that the CAVS force field could accurately predict the amino acid tendency of the secondary structure. Then, we used the CAVS force field to investigate the binding of a severe acute respiratory syndrome-associated coronavirus fusion peptide (SARS-CoV-2 FP) to a phospholipid bilayer: a long FP (FP-L) containing 40 amino acids and a short FP (FP-S) containing 26 amino acids. Our CAVS CG simulations displayed that the binding affinity of the FP-L to the bilayer is higher than that of the FP-S. We found that the FP-L interacted more strongly with membrane cholesterol than the FP-S, which should be attributed to the stable helical structure of the FP-L at the C-terminus. In addition, we discovered that the FP-S had one major and two minor membrane-bound states, in agreement with previous all-atom molecular dynamics (MD) studies. However, we found that both the C-terminal and N-terminal amino acid residues of the FP-L can strongly interact with the bilayer membrane. Furthermore, we found that the disulfide bond formed between Cys840 and Cys851 stabilized the helices of the FP-L at the C-terminus, enhancing the interaction between the FP-L and the bilayer membrane. Our work indicates that the stable helical structure is crucial for binding the SARS-CoV-2 FP to cell membranes. In particular, the helical stability of FP should have a significant influence on the FP-membrane binding.

COGRIMEN: Coarse-Grained Method for Modeling of Membrane Proteins in Implicit Environments

The presented methodology is based on coarse-grained representation of biomolecules in implicit environments and is designed for the molecular dynamics simulations of membrane proteins and their complexes. The membrane proteins are not only found in the cell membrane but also in all membranous compartments of the cell: Golgi apparatus, mitochondria, endosomes and lysosomes, and they usually form large complexes. To investigate such systems the methodology is proposed based on two independent approaches combining the coarse-grained MARTINI model for proteins and the effective energy function to mimic the water/membrane environments. The latter is based on the implicit environment developed for all-atom simulations in the IMM1 method. The force field solvation parameters for COGRIMEN were initially calculated from IMM1 all-atom parameters and then optimized using Genetic Algorithms. The new methodology was tested on membrane proteins, their complexes and oligomers. COGRIMEN method is implemented as a patch for NAMD program and can be useful for fast and brief studies of large membrane protein complexes.

Phytochemical Discrimination, Biological Activity and Molecular Docking of Water-Soluble Inhibitors from Saussurea costus Herb against Main Protease of SARS-CoV-2

Siddha medicine is one of the oldest medical systems in the world and is believed to have originated more than 10,000 years ago and is prevalent across ancient Tamil land. It is undeniable that inhibitor preferences rise with increasing solubility in water due to the considerations pertaining to the bioavailability and the ease of which unabsorbed residues can be disposed of. In this study, we showed the phytochemical discrimination of Saussurea costus extracted with water at room temperature as a green extraction procedure. A total of 48 compounds were identified using gas chromatography-mass spectrometry (GC-MS). The fatty acids had a high phytochemical abundance at 73.8%, followed by tannins at 8.2%, carbohydrates at 6.9%, terpenoids at 4.3%, carboxylic acids at 2.5%, hydrocarbons at 2.4%, phenolic compounds at 0.2%, and sterols at 1.5%. Of these compounds, 22 were docked on the active side and on the catalytic dyad of His41 and Cys145 of the main protease of SARS-CoV-2 (M^pro). Eight active inhibitors were carbohydrates, five were fatty acids, three were terpenoids, two were carboxylic acids, one was a tannin, one was a phenolic compound, and one was a sterol. The best inhibitors were 4,8,13-Cyclotetradecatriene-1,3-diol, 1,5,9-trimethyl-12-(1-methylethyl), Andrographolide, and delta.4-Androstene-3.beta.,17.beta.-diol, with a binding affinity that ranged from −6.1 kcal/mol to −6.5 kcal/mol. The inhibitory effect of Saussurea costus of SARS-CoV-2 entry into the cell was studied using a pseudovirus with Spike proteins from the D614G variant and the VOC variants Gamma and Delta. Based on the viral cycle of SARS-CoV-2, our results suggest that the Saussurea costus aqueous extract has no virucidal effect and inhibits the virus in the events after cell entry. Furthermore, the biological activity of the aqueous extract was investigated against HSV-1 virus and two bacterial strains, namely Staphylococcus aureus ATCC BAA 1026 and Escherichia coli ATCC 9637. According to this study, an enormous number of water-soluble inhibitors were identified from Saussurea costus against the M^pro, and this is unprecedented as far as we know.

Pathogenic Mechanism and Multi-omics Analysis of Oral Manifestations in COVID-19

Coronavirus disease 2019 (COVID-19) is a respiratory infectious disease that seriously threatens human life. The clinical manifestations of severe COVID-19 include acute respiratory distress syndrome and multiple organ failure. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causal agent of COVID-19, spreads through contaminated droplets. SARS-CoV-2 particles have been detected in the saliva of COVID-19 patients, implying that the virus can infect and damage the oral cavity. The oral manifestations of COVID-19 include xerostomia and gustatory dysfunction. Numerous studies showed that the four structural proteins of SARS-CoV-2 are its potential pathogenic factors, especially the S protein, which binds to human ACE2 receptors facilitating the entry of the virus into the host cells. Usually, upon entry into the host cell, a pathogen triggers the host’s immune response. However, a mount of multi-omics and immunological analyses revealed that COVID-19 is caused by immune dysregulation. A decrease in the number and phenotypes of immune cells, IFN-1 production and excessive release of certain cytokines have also been reported. In conclusion, this review summarizes the oral manifestations of COVID-19 and multi-omics analysis of SARS-CoV-2 infection.

Electrostatic Map of the SARS-CoV-2 Virion Specifies Binding Sites of the Antiviral Cationic Photosensitizer

Electrostatics is an important part of virus life. Understanding the detailed distribution of charges over the surface of a virus is important to predict its interactions with host cells, antibodies, drugs, and different materials. Using a coarse-grained model of the entire viral envelope developed by D. Korkin and S.-J. Marrink’s scientific groups, we created an electrostatic map of the external surface of SARS-CoV-2 and found a highly heterogeneous distribution of the electrostatic potential field of the viral envelope. Numerous negative patches originate mainly from negatively charged lipid domains in the viral membrane and negatively charged areas on the “stalks” of the spike (S) proteins. Membrane (M) and envelope (E) proteins with the total positive charge tend to colocalize with the negatively charged lipids. In the E protein pentamer exposed to the outer surface, negatively charged glutamate residues and surrounding lipids form a negative electrostatic potential ring around the channel entrance. We simulated the interaction of the antiviral octacationic photosensitizer octakis(cholinyl)zinc phthalocyanine with the surface structures of the entire model virion using the Brownian dynamics computational method implemented in ProKSim software (version r661). All mentioned negatively charged envelope components attracted the photosensitizer molecules and are thus potential targets for reactive oxygen generated in photosensitized reactions.

Multifaceted membrane binding head of the SARS-CoV-2 spike protein

The SARS-CoV-2 spike protein presents a surface with enormous membrane binding potential to host tissues and organelles of infected cells. Its exposed trimeric head binds not only the angiotensin-converting enzyme 2 (ACE2), but also host phospholipids which are missing from all existing structures. Hence, the membrane interaction surfaces that mediate viral fusion, entry, assembly and egress remain unclear. Here the spike:membrane docking sites are identified based on membrane optimal docking area (MODA) analysis of 3D structures of spike proteins in closed and open conformations at endocytic and neutral pH levels as well as ligand complexes. This reveals multiple membrane binding sites in the closed spike head that together prefer convex membranes and are modulated by pH, fatty acids and post-translational modifications including glycosylation. The exposure of the various membrane interaction sites adjusts upon domain repositioning within the trimer, allowing formation of intermediate bilayer complexes that lead to the prefusion state while also enabling ACE2 receptor recognition. In contrast, all antibodies that target the spike head would block the membrane docking process that precedes ACE2 recognition. Together this illuminates the engagements of the spike protein with plasma, endocytic, ER or exocytic vesicle membranes that help to drive the cycle of viral infection, and offers novel sites for intervention.