A multiscale coarse-grained model of the SARS-CoV-2 virion.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of the COVID-19 pandemic. Computer simulations of complete viral particles can provide theoretical insights into large-scale viral processes including assembly, budding, egress, entry, and fusion. Detailed atomistic simulations are constrained to shorter timescales and require billion-atom simulations for these processes. Here, we report the current status and ongoing development of a largely "bottom-up" coarse-grained (CG) model of the SARS-CoV-2 virion. Data from a combination of cryo-electron microscopy (cryo-EM), x-ray crystallography, and computational predictions were used to build molecular models of structural SARS-CoV-2 proteins, which were then assembled into a complete virion model. We describe how CG molecular interactions can be derived from all-atom simulations, how viral behavior difficult to capture in atomistic simulations can be incorporated into the CG models, and how the CG models can be iteratively improved as new data become publicly available. Our initial CG model and the detailed methods presented are intended to serve as a resource for researchers working on COVID-19 who are interested in performing multiscale simulations of the SARS-CoV-2 virion.

The flexibility of ACE2 in the context of SARS-CoV-2 infection.

The coronavirus disease 2019 (COVID-19) pandemic has swept over the world in the past months, causing significant loss of life and consequences to human health. Although numerous drug and vaccine development efforts are underway, there are many outstanding questions on the mechanism of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) viral association to angiotensin-converting enzyme 2 (ACE2), its main host receptor, and host cell entry. Structural and biophysical studies indicate some degree of flexibility in the viral extracellular spike glycoprotein and at the receptor-binding domain (RBD)-receptor interface, suggesting a role in infection. Here, we perform explicitly solvated, all-atom, molecular dynamics simulations of the glycosylated, full-length, membrane-bound ACE2 receptor in both an apo and spike RBD-bound state to probe the intrinsic dynamics of the ACE2 receptor in the context of the cell surface. A large degree of fluctuation in the full-length structure is observed, indicating hinge bending motions at the linker region connecting the head to the transmembrane helix while still not disrupting the ACE2 homodimer or ACE2-RBD interfaces. This flexibility translates into an ensemble of ACE2 homodimer conformations that could sterically accommodate binding of the spike trimer to more than one ACE2 homodimer and suggests a mechanical contribution of the host receptor toward the large spike conformational changes required for cell fusion. This work presents further structural and functional insights into the role of ACE2 in viral infection that can potentially be exploited for the rational design of effective SARS-CoV-2 therapeutics.

Calcium bridging drives polysaccharide co-adsorption to a proxy sea surface microlayer.

Saccharides comprise a significant mass fraction of organic carbon in sea spray aerosol (SSA), but the mechanisms through which saccharides are transferred from seawater to the ocean surface and eventually into SSA are unclear. It is hypothesized that saccharides cooperatively adsorb to other insoluble organic matter at the air/sea interface, known as the sea surface microlayer (SSML). Using a combination of surface-sensitive infrared reflection-absorption spectroscopy and all-atom molecular dynamics simulations, we demonstrate that the marine-relevant, anionic polysaccharide alginate co-adsorbs to an insoluble palmitic acid monolayer via divalent cationic bridging interactions. Ca2+ induces the greatest extent of alginate co-adsorption to the monolayer, evidenced by the ∼30% increase in surface coverage, whereas Mg2+ only facilitates one-third the extent of co-adsorption at seawater-relevant cation concentrations due to its strong hydration propensity. Na+ cations alone do not facilitate alginate co-adsorption, and palmitic acid protonation hinders the formation of divalent cationic bridges between the palmitate and alginate carboxylate moieties. Alginate co-adsorption is largely confined to the interfacial region beneath the monolayer headgroups, so surface pressure, and thus monolayer surface coverage, only changes the amount of alginate co-adsorption by less than 5%. Our results provide physical and molecular characterization of a potentially significant polysaccharide enrichment mechanism within the SSML.

AI-driven multiscale simulations illuminate mechanisms of SARS-CoV-2 spike dynamics.

We develop a generalizable AI-driven workflow that leverages heterogeneous HPC resources to explore the time-dependent dynamics of molecular systems. We use this workflow to investigate the mechanisms of infectivity of the SARS-CoV-2 spike protein, the main viral infection machinery. Our workflow enables more efficient investigation of spike dynamics in a variety of complex environments, including within a complete SARS-CoV-2 viral envelope simulation, which contains 305 million atoms and shows strong scaling on ORNL Summit using NAMD. We present several novel scientific discoveries, including the elucidation of the spike's full glycan shield, the role of spike glycans in modulating the infectivity of the virus, and the characterization of the flexible interactions between the spike and the human ACE2 receptor. We also demonstrate how AI can accelerate conformational sampling across different systems and pave the way for the future application of such methods to additional studies in SARS-CoV-2 and other molecular systems.

Surfactant Charge Modulates Structure and Stability of Lipase-Embedded Monolayers at Marine-Relevant Aerosol Surfaces.

Lipases, as well as other enzymes, are present and active within the sea surface microlayer (SSML). Upon bubble bursting, lipases partition into sea spray aerosol (SSA) along with surface-active molecules such as lipids. Lipases are likely to be embedded in the lipid monolayer at the SSA surface and thus have the potential to influence SSA interfacial structure and chemistry. Elucidating the structure of the lipid monolayer at SSA interfaces and how this structure is altered upon interaction with a protein system like lipase is of interest, given the importance of how aerosols interact with sunlight, influence cloud formation, and provide surfaces for chemical reactions. Herein, we report an integrated experimental and computational study of Burkholderia cepacia lipase (BCL) embedded in a lipid monolayer and highlight the important role of electrostatic, rather than hydrophobic, interactions as a driver for monolayer stability. Specifically, we combine Langmuir film experiments and molecular dynamics (MD) simulations to examine the detailed interactions between the zwitterionic dipalmitoylphosphatidylcholine (DPPC) monolayer and BCL. Upon insertion of BCL from the underlying subphase into the lipid monolayer, it is shown that BCL permeates and largely disorders the monolayer while strongly interacting with zwitterionic DPPC molecules, as experimentally observed by Langmuir adsorption curves and infrared reflectance absorbance spectroscopy. Explicitly solvated, all-atom MD is then used to provide insights into inter- and intramolecular interactions that drive these observations, with specific attention to the formation of salt bridges or ionic-bonding interactions. We show that after insertion into the DPPC monolayer, lipase is maintained at high surface pressures and in large BCL concentrations by forming a salt-bridge-stabilized lipase-DPPC complex. In comparison, when embedded in an anionic monolayer at low surface pressures, BCL preferentially forms intramolecular salt bridges, reducing its total favorable interactions with the surfactant and partitioning out of the monolayer shortly after injection. Overall, this study shows that the structure and dynamics of lipase-embedded SSA surfaces vary based on surface charge and pressure and that these variations have the potential to differentially modulate the properties of marine aerosols.

Revealing the Impacts of Chemical Complexity on Submicrometer Sea Spray Aerosol Morphology.

Sea spray aerosol (SSA) ejected through bursting bubbles at the ocean surface is a complex mixture of salts and organic species. Submicrometer SSA particles have long atmospheric lifetimes and play a critical role in the climate system. Composition impacts their ability to form marine clouds, yet their cloud-forming potential is difficult to study due to their small size. Here, we use large-scale molecular dynamics (MD) simulations as a "computational microscope" to provide never-before-seen views of 40 nm model aerosol particles and their molecular morphologies. We investigate how increasing chemical complexity impacts the distribution of organic material throughout individual particles for a range of organic constituents with varying chemical properties. Our simulations show that common organic marine surfactants readily partition between both the surface and interior of the aerosol, indicating that nascent SSA may be more heterogeneous than traditional morphological models suggest. We support our computational observations of SSA surface heterogeneity with Brewster angle microscopy on model interfaces. These observations indicate that increased chemical complexity in submicrometer SSA leads to a reduced surface coverage by marine organics, which may facilitate water uptake in the atmosphere. Our work thus establishes large-scale MD simulations as a novel technique for interrogating aerosols at the single-particle level.

Amino Acids Are Driven to the Interface by Salts and Acidic Environments.

Amino acids (AAs), the building blocks of proteins, are enriched by several orders of magnitude in sea spray aerosols compared to ocean waters. This suggests that AAs may reside at the air-water interface and be highly surface active. Using surface tension measurements, infrared reflection-absorption spectroscopy, and molecular dynamics simulations, we show that AAs are surface active and that salts and low-pH environments are drivers of surface activity. At typical sea spray salt concentrations and pH values, we determine that the surface coverage of hydrophobic AAs increases by approximately 1 order of magnitude. Additionally, divalent cations such as Ca2+ and Mg2+ can further increase AA surface propensity, particularly at neutral pH. Overall, these results indicate that AAs are likely to be found at increased concentrations at the surface of sea spray aerosols, where they can impact the cloud activation properties of the aerosol and enhance peptide formation under certain conditions.

Cation-Driven Lipopolysaccharide Morphological Changes Impact Heterogeneous Reactions of Nitric Acid with Sea Spray Aerosol Particles.

Lipopolysaccharides (LPS) in sea spray aerosol (SSA) particles have recently been shown to undergo heterogeneous reactions with HNO3 in the atmosphere. Here, we integrate theory and experiment to further investigate how the most abundant sea salt cations, Na+, Mg2+, and Ca2+, impact HNO3 reactions with LPS-containing SSA particles. Aerosol reaction flow tube studies show that heterogeneous reactions of SSA particles with divalent cation (Mg2+ and Ca2+) and LPS signatures were less reactive with HNO3 than those dominated by monovalent cations (Na+). All-atom molecular dynamics simulations of model LPS aggregates suggest that divalent cations cross-link the oligosaccharide chains to increase molecular aggregation and rigidity, which changes the particle phase and morphology, decreases water diffusion, and consequently decreases the reactive uptake of HNO3. This study provides new insight into how complex chemical interactions between ocean-derived salts and biogenic organic species can impact the heterogeneous reactivity of SSA particles.

Mesoscale All-Atom Influenza Virus Simulations Suggest New Substrate Binding Mechanism.

Influenza virus circulates in human, avian, and swine hosts, causing seasonal epidemic and occasional pandemic outbreaks. Influenza neuraminidase, a viral surface glycoprotein, has two sialic acid binding sites. The catalytic (primary) site, which also binds inhibitors such as oseltamivir carboxylate, is responsible for cleaving the sialic acid linkages that bind viral progeny to the host cell. In contrast, the functional annotation of the secondary site remains unclear. Here, we better characterize these two sites through the development of an all-atom, explicitly solvated, and experimentally based integrative model of the pandemic influenza A H1N1 2009 viral envelope, containing ∼160 million atoms and spanning ∼115 nm in diameter. Molecular dynamics simulations of this crowded subcellular environment, coupled with Markov state model theory, provide a novel framework for studying realistic molecular systems at the mesoscale and allow us to quantify the kinetics of the neuraminidase 150-loop transition between the open and closed states. An analysis of chloride ion occupancy along the neuraminidase surface implies a potential new role for the neuraminidase secondary site, wherein the terminal sialic acid residues of the linkages may bind before transfer to the primary site where enzymatic cleavage occurs. Altogether, our work breaks new ground for molecular simulation in terms of size, complexity, and methodological analyses of the components. It also provides fundamental insights into the understanding of substrate recognition processes for this vital influenza drug target, suggesting a new strategy for the development of anti-influenza therapeutics.

A Multiscale Coarse-grained Model of the SARS-CoV-2 Virion.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of the COVID-19 pandemic. Computer simulations of complete viral particles can provide theoretical insights into large-scale viral processes including assembly, budding, egress, entry, and fusion. Detailed atomistic simulations, however, are constrained to shorter timescales and require billion-atom simulations for these processes. Here, we report the current status and on-going development of a largely "bottom-up" coarse-grained (CG) model of the SARS-CoV-2 virion. Structural data from a combination of cryo-electron microscopy (cryo-EM), x-ray crystallography, and computational predictions were used to build molecular models of structural SARS-CoV-2 proteins, which were then assembled into a complete virion model. We describe how CG molecular interactions can be derived from all-atom simulations, how viral behavior difficult to capture in atomistic simulations can be incorporated into the CG models, and how the CG models can be iteratively improved as new data becomes publicly available. Our initial CG model and the detailed methods presented are intended to serve as a resource for researchers working on COVID-19 who are interested in performing multiscale simulations of the SARS-CoV-2 virion. SIGNIFICANCE STATEMENT: This study reports the construction of a molecular model for the SARS-CoV-2 virion and details our multiscale approach towards model refinement. The resulting model and methods can be applied to and enable the simulation of SARS-CoV-2 virions.

Insights into the behavior of nonanoic acid and its conjugate base at the air/water interface through a combined experimental and theoretical approach.

The partitioning of medium-chain fatty acid surfactants such as nonanoic acid (NA) between the bulk phase and the air/water interface is of interest to a number of fields including marine and atmospheric chemistry. However, questions remain about the behavior of these molecules, the contributions of various relevant chemical equilibria, and the impact of pH, salt and bulk surfactant concentrations. In this study, the surface adsorption of nonanoic acid and its conjugate base is quantitatively investigated at various pH values, surfactant concentrations and the presence of salts. Surface concentrations of protonated and deprotonated species are dictated by surface-bulk equilibria which can be calculated from thermodynamic considerations. Notably we conclude that the surface dissociation constant of soluble surfactants cannot be directly obtained from these experimental measurements, however, we show that molecular dynamics (MD) simulation methods, such as free energy perturbation (FEP), can be used to calculate the surface acid dissociation constant relative to that in the bulk. These simulations show that nonanoic acid is less acidic at the surface compared to in the bulk solution with a pK a shift of 1.1 ± 0.6, yielding a predicted surface pK a of 5.9 ± 0.6. A thermodynamic cycle for nonanoic acid and its conjugate base between the air/water interface and the bulk phase can therefore be established. Furthermore, the effect of salts, namely NaCl, on the surface activity of protonated and deprotonated forms of nonanoic acid is also examined. Interestingly, salts cause both a decrease in the bulk pK a of nonanoic acid and a stabilization of both the protonated and deprotonated forms at the surface. Overall, these results suggest that the deprotonated medium-chain fatty acids under ocean conditions can also be present within the sea surface microlayer (SSML) present at the ocean/atmosphere interface due to the stabilization effect of the salts in the ocean. This allows the transfer of these species into sea spray aerosols (SSAs). More generally, we present a framework with which the behavior of partially soluble species at the air/water interface can be predicted from surface adsorption models and the surface pK a can be predicted from MD simulations.

Beyond Shielding: The Roles of Glycans in the SARS-CoV-2 Spike Protein.

The ongoing COVID-19 pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has resulted in more than 28,000,000 infections and 900,000 deaths worldwide to date. Antibody development efforts mainly revolve around the extensively glycosylated SARS-CoV-2 spike (S) protein, which mediates host cell entry by binding to the angiotensin-converting enzyme 2 (ACE2). Similar to many other viral fusion proteins, the SARS-CoV-2 spike utilizes a glycan shield to thwart the host immune response. Here, we built a full-length model of the glycosylated SARS-CoV-2 S protein, both in the open and closed states, augmenting the available structural and biological data. Multiple microsecond-long, all-atom molecular dynamics simulations were used to provide an atomistic perspective on the roles of glycans and on the protein structure and dynamics. We reveal an essential structural role of N-glycans at sites N165 and N234 in modulating the conformational dynamics of the spike's receptor binding domain (RBD), which is responsible for ACE2 recognition. This finding is corroborated by biolayer interferometry experiments, which show that deletion of these glycans through N165A and N234A mutations significantly reduces binding to ACE2 as a result of the RBD conformational shift toward the "down" state. Additionally, end-to-end accessibility analyses outline a complete overview of the vulnerabilities of the glycan shield of the SARS-CoV-2 S protein, which may be exploited in the therapeutic efforts targeting this molecular machine. Overall, this work presents hitherto unseen functional and structural insights into the SARS-CoV-2 S protein and its glycan coat, providing a strategy to control the conformational plasticity of the RBD that could be harnessed for vaccine development.

Beyond Shielding: The Roles of Glycans in SARS-CoV-2 Spike Protein.

The ongoing COVID-19 pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has resulted in more than 15,000,000 infections and 600,000 deaths worldwide to date. Antibody development efforts mainly revolve around the extensively glycosylated SARS-CoV-2 spike (S) protein, which mediates the host cell entry by binding to the angiotensin-converting enzyme 2 (ACE2). Similar to many other viruses, the SARS-CoV-2 spike utilizes a glycan shield to thwart the host immune response. Here, we built a full-length model of glycosylated SARS-CoV-2 S protein, both in the open and closed states, augmenting the available structural and biological data. Multiple microsecond-long, all-atom molecular dynamics simulations were used to provide an atomistic perspective on the roles of glycans, and the protein structure and dynamics. We reveal an essential structural role of N-glycans at sites N165 and N234 in modulating the conformational dynamics of the spike's receptor binding domain (RBD), which is responsible for ACE2 recognition. This finding is corroborated by biolayer interferometry experiments, which show that deletion of these glycans through N165A and N234A mutations significantly reduces binding to ACE2 as a result of the RBD conformational shift towards the "down" state. Additionally, end-to-end accessibility analyses outline a complete overview of the vulnerabilities of the glycan shield of SARS-CoV-2 S protein, which may be exploited by therapeutic efforts targeting this molecular machine. Overall, this work presents hitherto unseen functional and structural insights into the SARS-CoV-2 S protein and its glycan coat, providing a strategy to control the conformational plasticity of the RBD that could be harnessed for vaccine development.

Sea Spray Aerosol Structure and Composition Using Cryogenic Transmission Electron Microscopy.

The composition and surface properties of atmospheric aerosol particles largely control their impact on climate by affecting their ability to uptake water, react heterogeneously, and nucleate ice in clouds. However, in the vacuum of a conventional electron microscope, the native surface and internal structure often undergo physicochemical rearrangement resulting in surfaces that are quite different from their atmospheric configurations. Herein, we report the development of cryogenic transmission electron microscopy where laboratory generated sea spray aerosol particles are flash frozen in their native state with iterative and controlled thermal and/or pressure exposures and then probed by electron microscopy. This unique approach allows for the detection of not only mixed salts, but also soft materials including whole hydrated bacteria, diatoms, virus particles, marine vesicles, as well as gel networks within hydrated salt droplets-all of which will have distinct biological, chemical, and physical processes. We anticipate this method will open up a new avenue of analysis for aerosol particles, not only for ocean-derived aerosols, but for those produced from other sources where there is interest in the transfer of organic or biological species from the biosphere to the atmosphere.