Energy landscape of the SARS-CoV-2 reveals extensive conformational heterogeneity

Structural basis of closed groove scrambling by a TMEM16 protein

Activation of Ca2+-dependent TMEM16 scramblases induces phosphatidylserine externalization, a key step in multiple signaling processes. Current models suggest that the TMEM16s scramble lipids by deforming the membrane near a hydrophilic groove and that Ca2+ dependence arises from the different association of lipids with an open or closed groove. However, the molecular rearrangements underlying groove opening and how lipids reorganize outside the closed groove remain unknown. Here we directly visualize how lipids associate at the closed groove of Ca2+-bound fungal nhTMEM16 in nanodiscs using cryo-EM. Functional experiments pinpoint lipid-protein interaction sites critical for closed groove scrambling. Structural and functional analyses suggest groove opening entails the sequential appearance of two π-helical turns in the groove-lining TM6 helix and identify critical rearrangements. Finally, we show that the choice of scaffold protein and lipids affects the conformations of nhTMEM16 and their distribution, highlighting a key role of these factors in cryo-EM structure determination.

Structural basis of closed groove scrambling by a TMEM16 protein

Activation of Ca2+-dependent TMEM16 scramblases induces the externalization of phosphatidylserine, a key molecule in multiple signaling processes. Current models suggest that the TMEM16s scramble lipids by deforming the membrane near a hydrophilic groove, and that Ca2+ dependence arises from the different association of lipids with an open or closed groove. However, the molecular rearrangements involved in groove opening and of how lipids reorganize outside the closed groove remain unknown. Using cryogenic electron microscopy, we directly visualize how lipids associate at the closed groove of Ca2+-bound nhTMEM16 in nanodiscs. Functional experiments pinpoint the lipid-protein interaction sites critical for closed groove scrambling. Structural and functional analyses suggest groove opening entails the sequential appearance of two π-helical turns in the groove-lining TM6 helix and identify critical rearrangements. Finally, we show that the choice of scaffold protein and lipids affects the conformations of nhTMEM16 and their distribution, highlighting a key role of these factors in cryoEM structure determination.

Adaptive Ensemble Refinement of Protein Structures in High Resolution Electron Microscopy Density Maps with Radical Augmented Molecular Dynamics Flexible Fitting

Recent advances in cryo-electron microscopy (cryo-EM) have enabled modeling macromolecular complexes that are essential components of the cellular machinery. The density maps derived from cryo-EM experiments are often integrated with manual, knowledge-driven or artificial intelligence-driven and physics-guided computational methods to build, fit, and refine molecular structures. Going beyond a single stationary-structure determination scheme, it is becoming more common to interpret the experimental data with an ensemble of models that contributes to an average observation. Hence, there is a need to decide on the quality of an ensemble of protein structures on-the-fly while refining them against the density maps. We introduce such an adaptive decision-making scheme during the molecular dynamics flexible fitting (MDFF) of biomolecules. Using RADICAL-Cybertools, the new RADICAL augmented MDFF implementation (R-MDFF) is examined in high-performance computing environments for refinement of two prototypical protein systems, adenylate kinase and carbon monoxide dehydrogenase. For these test cases, use of multiple replicas in flexible fitting with adaptive decision making in R-MDFF improves the overall correlation to the density by 40% relative to the refinements of the brute-force MDFF. The improvements are particularly significant at high, 2-3 Å map resolutions. More importantly, the ensemble model captures key features of biologically relevant molecular dynamics that are inaccessible to a single-model interpretation. Finally, the pipeline is applicable to systems of growing sizes, which is demonstrated using ensemble refinement of capsid proteins from the chimpanzee adenovirus. The overhead for decision making remains low and robust to computing environments. The software is publicly available on GitHub and includes a short user guide to install R-MDFF on different computing environments, from local Linux-based workstations to high-performance computing environments.

Structural basis of closed groove scrambling by a TMEM16 protein

Activation of Ca2+-dependent TMEM16 scramblases induces the externalization of phosphatidylserine, a key molecule in multiple signaling processes. Current models suggest that the TMEM16s scramble lipids by deforming the membrane near a hydrophilic groove, and that Ca2+ dependence arises from the different association of lipids with an open or closed groove. However, the molecular rearrangements involved in groove opening and of how lipids reorganize outside the closed groove remain unknown. Using cryogenic electron microscopy, we directly visualize how lipids associate at the closed groove of Ca2+-bound nhTMEM16 in nanodiscs. Functional experiments pinpoint the lipid-protein interaction sites critical for closed groove scrambling. Structural and functional analyses suggest groove opening entails the sequential appearance of two π-helical turns in the groove-lining TM6 helix and identify critical rearrangements. Finally, we show that the choice of scaffold protein and lipids affects the conformations of nhTMEM16 and their distribution, highlighting a key role of these factors in cryoEM structure determination.

Conformational heterogeneity and probability distributions from single-particle cryo-electron microscopy

Single-particle cryo-electron microscopy (cryo-EM) is a technique that takes projection images of biomolecules frozen at cryogenic temperatures. A major advantage of this technique is its ability to image single biomolecules in heterogeneous conformations. While this poses a challenge for data analysis, recent algorithmic advances have enabled the recovery of heterogeneous conformations from the noisy imaging data. Here, we review methods for the reconstruction and heterogeneity analysis of cryo-EM images, ranging from linear-transformation-based methods to nonlinear deep generative models. We overview the dimensionality-reduction techniques used in heterogeneous 3D reconstruction methods and specify what information each method can infer from the data. Then, we review the methods that use cryo-EM images to estimate probability distributions over conformations in reduced subspaces or predefined by atomistic simulations. We conclude with the ongoing challenges for the cryo-EM community.

Delineating the impact of pathogenic mutations on the conformational dynamics of HDL's vital protein ApoA1: a combined computational and molecular dynamic simulation approach

Apolipoprotein A1 (ApoA1), is the important component of high-density lipoproteins (HDL), that has key role in HDL biogenesis, cholesterol trafficking, and reverse cholesterol transport (RCT). Non-synonymous Single Nucleotide Polymorphisms (nsSNPs) in ApoA1 have been linked to cardiovascular diseases and amyloidosis as they alter the protein's native structure and function. Therefore in this study, we attempted to understand the molecular pathogenicity profile of nsSNPs of ApoA1 using various computational approaches. We used state-of-the-art computational methods to thoroughly investigate the 295 ApoA1 nsSNPs at sequence and structural levels. Seven nsSNPs (L13R, L84R, L84P, L99P, R173P, L187P, and L238P) out of 295 were classified as the most deleterious and destabilizing. In order to estimate the effect of such destabilizing mutations on the protein conformation, all-atom molecular dynamics simulations (MDS) of ApoA1 wild-type (WT), L99P and R173P for 100 ns, was carried out using GROMACS 5.0.1 package. The MD simulation investigation revealed significant structural alterations in L99P and R173P. In addition, they had changed principal component analysis and electrostatic surface potential, decreased structural compactness, and intramolecular hydrogen bonds, which supported the rationale underpinning ApoA1 dysfunction with such mutations. This work sheds light on ApoA1 dysfunction due to single amino acid alterations, and offers new insight into the molecular basis of ApoA1-related diseases progression.Communicated by Ramaswamy H. Sarma.

SARS-CoV-2 spike opening dynamics and energetics reveal the individual roles of glycans and their collective impact

The trimeric spike (S) glycoprotein, which protrudes from the SARS-CoV-2 viral envelope, binds to human ACE2, initiated by at least one protomer's receptor binding domain (RBD) switching from a "down" (closed) to an "up" (open) state. Here, we used large-scale molecular dynamics simulations and two-dimensional replica exchange umbrella sampling calculations with more than a thousand windows and an aggregate total of 160 μs of simulation to investigate this transition with and without glycans. We find that the glycosylated spike has a higher barrier to opening and also energetically favors the down state over the up state. Analysis of the S-protein opening pathway reveals that glycans at N165 and N122 interfere with hydrogen bonds between the RBD and the N-terminal domain in the up state, while glycans at N165 and N343 can stabilize both the down and up states. Finally, we estimate how epitope exposure for several known antibodies changes along the opening path. We find that the BD-368-2 antibody's epitope is continuously exposed, explaining its high efficacy.