Machine Learning-Driven Multiscale Modeling: Bridging the Scales with a Next-Generation Simulation Infrastructure
Not Available
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Engineering topics
Nguyen, Kien
Publications and source records attributed to Nguyen, Kien.
Ratchet, swivel, tilt and roll: a complete description of subunit rotation in the ribosome
Protein synthesis by the ribosome requires large-scale rearrangements of the ‘small’ subunit (SSU; ~1 MDa), including inter- and intra-subunit rotational motions. However, with nearly 2000 structures of ribosomes and ribosomal subunits now publicly available, it is exceedingly difficult to design experiments based on analysis of all known rotation states. To overcome this, we developed an approach where the orientation of each SSU head and body is described in terms of three angular coordinates (rotation, tilt and tilt direction) and a single translation. By considering the entire RCSB PDB database, we describe 1208 fully-assembled ribosome complexes and 334 isolated small subunits, which span >50 species. This reveals aspects of subunit rearrangements that are universal, and others that are organism/domain-specific. For example, we show that tilt-like rearrangements of the SSU body (i.e. ‘rolling’) are pervasive in both prokaryotic and eukaryotic (cytosolic and mitochondrial) ribosomes. As another example, domain orientations associated with frameshifting in bacteria are similar to those found in eukaryotic ribosomes. Together, this study establishes a common foundation with which structural, simulation, single-molecule and biochemical efforts can more precisely interrogate the dynamics of this prototypical molecular machine.
59 BASIC BIOLOGICAL SCIENCES↗
Network analysis uncovers the communication structure of SARS-CoV-2 spike protein identifying sites for immunogen design
The COVID-19 pandemic, caused by the SARS-CoV-2 virus, has triggered myriad efforts to understand the structure and dynamics of this complex pathogen. The spike glycoprotein of SARS-CoV-2 is a significant target for immunogens as it is the means by which the virus enters human cells, while simultaneously sporting mutations responsible for immune escape. These functional and escape processes are regulated by complex molecular-level interactions. Our study presents quantitative insights on domain and residue contributions to allosteric communication, immune evasion, and local- and global-level control of functions through the derivation of a weighted graph representation from all-atom MD simulations. Focusing on the ancestral form and the D614G-variant, we provide evidence of the utility of our approach by guiding the selection of a mutation that alters the spike’s stability. Taken together, the network approach serves as a valuable tool to evaluate communication “hot-spots” in proteins to guide design of stable immunogens.
60 APPLIED LIFE SCIENCES↗
Asynchronous Reciprocal Coupling of Martini 2.2 Coarse-Grained and CHARMM36 All-Atom Simulations in an Automated Multiscale Framework
The appeal of multiscale modeling approaches is predicated on the promise of combinatorial synergy. However, this promise can only be realized when distinct scales are combined with reciprocal consistency. Here, we consider multiscale molecular dynamics (MD) simulations that combine the accuracy and macromolecular flexibility accessible to fixed-charge all-atom (AA) representations with the sampling speed accessible to reductive, coarse-grained (CG) representations. AA-to-CG conversions are relatively straightforward because deterministic routines with unique outcomes are achievable. Conversely, CG-to-AA conversions have many solutions due to a surge in the number of degrees of freedom. While automated tools for biomolecular CG-to-AA transformation exist, we find that one popular option, called Backward, is prone to stochastic failure and the AA models that it does generate frequently have compromised protein structure and incorrect stereochemistry. Although these shortcomings can likely be circumvented by human intervention in isolated instances, automated multiscale coupling requires reliable and robust scale conversion. In this work, we detail an extension to Multiscale Machine-learned Modeling Infrastructure (MuMMI), including an improved CG-to-AA conversion tool called sinceCG. This tool is reliable (~98% weakly correlated repeat success rate), automatable (no unrecoverable hangs), and yields AA models that generally preserve protein secondary structure and maintain correct stereochemistry. We describe how the MuMMI framework identifies CG system configurations of interest, converts them to AA representations, and simulates them at the AA scale while on-the-fly analyses provide feedback to update CG parameters. Application to systems containing the peripheral membrane protein RAS and proximal components of RAF kinase on complex eight-component lipid bilayers with ~1.5 million atoms is discussed in the context of MuMMI.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
HIV-1 and SARS-CoV-2: Patterns in the evolution of two pandemic pathogens
Not Available
59 BASIC BIOLOGICAL SCIENCES↗
Exploring the Role of Glycans in the Interaction of SARS-CoV-2 RBD and Human Receptor ACE2
COVID-19 is a highly infectious respiratory disease caused by the novel coronavirus SARS-CoV-2. It has become a global pandemic and its frequent mutations may pose new challenges for vaccine design. During viral infection, the Spike RBD of SARS-CoV-2 binds the human host cell receptor ACE2, enabling the virus to enter the host cell. Both the Spike and ACE2 are densely glycosylated, and it is unclear how distinctive glycan types may modulate the interaction of RBD and ACE2. Detailed understanding of these determinants is key for the development of novel therapeutic strategies. To this end, we perform extensive all-atom simulations of the (i) RBD-ACE2 complex without glycans, (ii) RBD-ACE2 with oligomannose MAN9 glycans in ACE2, and (iii) RBD-ACE2 with complex FA2 glycans in ACE2. These simulations identify the key residues at the RBD-ACE2 interface that form contacts with higher probabilities, thus providing a quantitative evaluation that complements recent structural studies. Notably, we find that this RBD-ACE2 contact signature is not altered by the presence of different glycoforms, suggesting that RBD-ACE2 interaction is robust. Applying our simulated results, we illustrate how the recently prevalent N501Y mutation may alter specific interactions with host ACE2 that facilitate the virus-host binding. Furthermore, our simulations reveal how the glycan on Asn90 of ACE2 can play a distinct role in the binding and unbinding of RBD. Finally, an energetics analysis shows that MAN9 glycans on ACE2 decrease RBD-ACE2 affinity, while FA2 glycans lead to enhanced binding of the complex. Together, our results provide a more comprehensive picture of the detailed interplay between virus and human receptor, which is much needed for the discovery of effective treatments that aim at modulating the physical-chemical properties of this virus.
60 APPLIED LIFE SCIENCES↗
The SARS-CoV-2 Spike variant D614G favors an open conformational state
Symmetrization of contacts in D614G SARS-CoV-2 Spike favors greater probability of infection-capable conformation.
60 APPLIED LIFE SCIENCES↗
The SARS-CoV-2 Spike Variant D614G Favors an Open Conformational State
The COVID-19 pandemic underwent a rapid transition with the emergence of a dominant viral variant (from the “D-form” to the “G-form”) that carried an amino acid substitution D614G in its “Spike” protein. The G-form is more infectious in vitro and associated with increased viral loads. To gain insight into the molecular-level underpinnings of these characteristics, we employed microsecond all-atom simulations. Here we show that changes in the protein energetics favor a higher population of infection-capable (open) states through release of hydrogen bonds of an asymmetry present in the D-form but not the G-form. Thus, the increased infectivity of the G-form is likely due to a higher rate of profitable binding encounters with the host receptor. It is also predicted to be more neutralization sensitive due to enhanced exposure of the receptor binding domain, a key target region for neutralizing antibodies. These results are significant for vaccine design. The Molecular Dynamics Simulations datasets generated in this study, namely that of the soluble form SARS-CoV-2 Spike protein are made available here. Details of the trajectory and supporting files are given in the README.txt file.
59 BASIC BIOLOGICAL SCIENCES↗
Root cause analysis
Explore the source record for details and available documents.
root cause analysis↗