Understanding the cell: Future views of structural biology.

Determining the structure and mechanisms of all individual functional modules of cells at high molecular detail has often been seen as equal to understanding how cells work. Recent technical advances have led to a flush of high-resolution structures of various macromolecular machines, but despite this wealth of detailed information, our understanding of cellular function remains incomplete. Here, we discuss present-day limitations of structural biology and highlight novel technologies that may enable us to analyze molecular functions directly inside cells. We predict that the progression toward structural cell biology will involve a shift toward conceptualizing a 4D virtual reality of cells using digital twins. These will capture cellular segments in a highly enriched molecular detail, include dynamic changes, and facilitate simulations of molecular processes, leading to novel and experimentally testable predictions. Transferring biological questions into algorithms that learn from the existing wealth of data and explore novel solutions may ultimately unveil how cells work.

Disordered regions in the IRE1α ER lumenal domain mediate its stress-induced clustering.

Conserved signaling cascades monitor protein-folding homeostasis to ensure proper cellular function. One of the evolutionary conserved key players is IRE1, which maintains endoplasmic reticulum (ER) homeostasis through the unfolded protein response (UPR). Upon accumulation of misfolded proteins in the ER, IRE1 forms clusters on the ER membrane to initiate UPR signaling. What regulates IRE1 cluster formation is not fully understood. Here, we show that the ER lumenal domain (LD) of human IRE1α forms biomolecular condensates in vitro. IRE1α LD condensates were stabilized both by binding to unfolded polypeptides as well as by tethering to model membranes, suggesting their role in assembling IRE1α into signaling-competent stable clusters. Molecular dynamics simulations indicated that weak multivalent interactions drive IRE1α LD clustering. Mutagenesis experiments identified disordered regions in IRE1α LD to control its clustering in vitro and in cells. Importantly, dysregulated clustering of IRE1α mutants led to defects in IRE1α signaling. Our results revealed that disordered regions in IRE1α LD control its clustering and suggest their role as a common strategy in regulating protein assembly on membranes.

Osh Proteins Control Nanoscale Lipid Organization Necessary for PI(4,5)P2 Synthesis.

The plasma membrane (PM) is composed of a complex lipid mixture that forms heterogeneous membrane environments. Yet, how small-scale lipid organization controls physiological events at the PM remains largely unknown. Here, we show that ORP-related Osh lipid exchange proteins are critical for the synthesis of phosphatidylinositol (4,5)-bisphosphate [PI(4,5)P2], a key regulator of dynamic events at the PM. In real-time assays, we find that unsaturated phosphatidylserine (PS) and sterols, both Osh protein ligands, synergistically stimulate phosphatidylinositol 4-phosphate 5-kinase (PIP5K) activity. Biophysical FRET analyses suggest an unconventional co-distribution of unsaturated PS and phosphatidylinositol 4-phosphate (PI4P) species in sterol-containing membrane bilayers. Moreover, using in vivo imaging approaches and molecular dynamics simulations, we show that Osh protein-mediated unsaturated PI4P and PS membrane lipid organization is sensed by the PIP5K specificity loop. Thus, ORP family members create a nanoscale membrane lipid environment that drives PIP5K activity and PI(4,5)P2 synthesis that ultimately controls global PM organization and dynamics.

Integrated Functions of Membrane Property Sensors and a Hidden Side of the Unfolded Protein Response.

Eukaryotic cells face the challenge of maintaining the complex composition of several coexisting organelles. The molecular mechanisms underlying the homeostasis of subcellular membranes and their adaptation during stress are only now starting to emerge. Here, we discuss three membrane property sensors of the endoplasmic reticulum (ER), namely OPI1, MGA2, and IRE1, each controlling a large cellular program impacting the lipid metabolic network. OPI1 coordinates the production of membrane and storage lipids, MGA2 regulates the production of unsaturated fatty acids required for membrane biogenesis, and IRE1 controls the unfolded protein response (UPR) to adjust ER size, protein folding, and the secretory capacity of the cell. Although these proteins use remarkably distinct sensing mechanisms, they are functionally connected via the ER membrane and cooperate to maintain membrane homeostasis. As a rationalization of the recently described mechanism of UPR activation by lipid bilayer stress, we propose that IRE1 can sense the protein-to-lipid ratio in the ER membrane to ensure a balanced production of membrane proteins and lipids.

Activation of the Unfolded Protein Response by Lipid Bilayer Stress.

The unfolded protein response (UPR) is a conserved homeostatic program that is activated by misfolded proteins in the lumen of the endoplasmic reticulum (ER). Recently, it became evident that aberrant lipid compositions of the ER membrane, referred to as lipid bilayer stress, are equally potent in activating the UPR. The underlying molecular mechanism, however, remained unclear. We show that the most conserved transducer of ER stress, Ire1, uses an amphipathic helix (AH) to sense membrane aberrancies and control UPR activity. In vivo and in vitro experiments, together with molecular dynamics (MD) simulations, identify the physicochemical properties of the membrane environment that control Ire1 oligomerization. This work establishes the molecular mechanism of UPR activation by lipid bilayer stress.

Antibody accessibility determines location of spike surface mutations in SARS-CoV-2 variants.

The steady emergence of SARS-CoV-2 variants gives us a real-time view of the interplay between viral evolution and the host immune defense. The spike protein of SARS-CoV-2 is the primary target of antibodies. Here, we show that steric accessibility to antibodies provides a strong predictor of mutation activity in the spike protein of SARS-CoV-2 variants, including Omicron. We introduce an antibody accessibility score (AAS) that accounts for the steric shielding effect of glycans at the surface of spike. We find that high values of the AAS correlate strongly with the sites of mutations in the spike proteins of newly emerging SARS-CoV-2 variants. We use the AAS to assess the escapability of variant spike proteins, i.e., their ability to escape antibody-based immune responses. The high calculated escapability of the Omicron variant BA.5 with respect to both wild-type (WT) vaccination and BA.1 infection is consistent with its rapid spread despite high rates of vaccination and prior infection with earlier variants. We calculated the AAS from structural and molecular dynamics simulation data that were available early in the pandemic, in the spring of 2020. The AAS thus allows us to prospectively assess the ability of variant spike proteins to escape antibody-based immune responses and to pinpoint regions of expected mutation activity in future variants.

A Eukaryotic Sensor for Membrane Lipid Saturation.

Maintaining a fluid bilayer is essential for cell signaling and survival. Lipid saturation is a key factor determining lipid packing and membrane fluidity, and it must be tightly controlled to guarantee organelle function and identity. A dedicated eukaryotic mechanism of lipid saturation sensing, however, remains elusive. Here we show that Mga2, a transcription factor conserved among fungi, acts as a lipid-packing sensor in the ER membrane to control the production of unsaturated fatty acids. Systematic mutagenesis, molecular dynamics simulations, and electron paramagnetic resonance spectroscopy identify a pivotal role of the oligomeric transmembrane helix (TMH) of Mga2 for intra-membrane sensing, and they show that the lipid environment controls the proteolytic activation of Mga2 by stabilizing alternative rotational orientations of the TMH region. This work establishes a eukaryotic strategy of lipid saturation sensing that differs significantly from the analogous bacterial mechanism relying on hydrophobic thickness.

Computational epitope map of SARS-CoV-2 spike protein.

The primary immunological target of COVID-19 vaccines is the SARS-CoV-2 spike (S) protein. S is exposed on the viral surface and mediates viral entry into the host cell. To identify possible antibody binding sites, we performed multi-microsecond molecular dynamics simulations of a 4.1 million atom system containing a patch of viral membrane with four full-length, fully glycosylated and palmitoylated S proteins. By mapping steric accessibility, structural rigidity, sequence conservation, and generic antibody binding signatures, we recover known epitopes on S and reveal promising epitope candidates for structure-based vaccine design. We find that the extensive and inherently flexible glycan coat shields a surface area larger than expected from static structures, highlighting the importance of structural dynamics. The protective glycan shield and the high flexibility of its hinges give the stalk overall low epitope scores. Our computational epitope-mapping procedure is general and should thus prove useful for other viral envelope proteins whose structures have been characterized.

Information-theoretical measures identify accurate low-resolution representations of protein configurational space.

The steadily growing computational power employed to perform molecular dynamics simulations of biological macromolecules represents at the same time an immense opportunity and a formidable challenge. In fact, large amounts of data are produced, from which useful, synthetic, and intelligible information has to be extracted to make the crucial step from knowing to understanding. Here we tackled the problem of coarsening the conformational space sampled by proteins in the course of molecular dynamics simulations. We applied different schemes to cluster the frames of a dataset of protein simulations; we then employed an information-theoretical framework, based on the notion of resolution and relevance, to gauge how well the various clustering methods accomplish this simplification of the configurational space. Our approach allowed us to identify the level of resolution that optimally balances simplicity and informativeness; furthermore, we found that the most physically accurate clustering procedures are those that induce an ultrametric structure of the low-resolution space, consistently with the hypothesis that the protein conformational landscape has a self-similar organisation. The proposed strategy is general and its applicability extends beyond that of computational biophysics, making it a valuable tool to extract useful information from large datasets.

Intrinsic map dynamics exploration for uncharted effective free-energy landscapes.

We describe and implement a computer-assisted approach for accelerating the exploration of uncharted effective free-energy surfaces (FESs). More generally, the aim is the extraction of coarse-grained, macroscopic information from stochastic or atomistic simulations, such as molecular dynamics (MD). The approach functionally links the MD simulator with nonlinear manifold learning techniques. The added value comes from biasing the simulator toward unexplored phase-space regions by exploiting the smoothness of the gradually revealed intrinsic low-dimensional geometry of the FES.

Molecular free energy profiles from force spectroscopy experiments by inversion of observed committors.

In single-molecule force spectroscopy experiments, a biomolecule is attached to a force probe via polymer linkers and the total extension of the molecule plus apparatus is monitored as a function of time. In a typical unfolding experiment at constant force, the total extension jumps between two values that correspond to the folded and unfolded states of the molecule. For several biomolecular systems, the committor, which is the probability to fold starting from a given extension, has been used to extract the molecular activation barrier (a technique known as "committor inversion"). In this work, we study the influence of the force probe, which is much larger than the molecule being measured, on the activation barrier obtained by committor inversion. We use a two-dimensional framework in which the diffusion coefficient of the molecule and of the pulling device can differ. We systematically study the free energy profile along the total extension obtained from the committor by numerically solving the Onsager equation and using Brownian dynamics simulations. We analyze the dependence of the extracted barrier on the linker stiffness, molecular barrier height, and diffusion anisotropy and, thus, establish the range of validity of committor inversion. Along the way, we showcase the committor of 2-dimensional diffusive models and illustrate how it is affected by barrier asymmetry and diffusion anisotropy.

Dominant folding pathways of a WW domain.

We investigate the folding mechanism of the WW domain Fip35 using a realistic atomistic force field by applying the Dominant Reaction Pathways approach. We find evidence for the existence of two folding pathways, which differ by the order of formation of the two hairpins. This result is consistent with the analysis of the experimental data on the folding kinetics of WW domains and with the results obtained from large-scale molecular dynamics simulations of this system. Free-energy calculations performed in two coarse-grained models support the robustness of our results and suggest that the qualitative structure of the dominant paths are mostly shaped by the native interactions. Computing a folding trajectory in atomistic detail only required about one hour on 48 Central Processing Units. The gain in computational efficiency opens the door to a systematic investigation of the folding pathways of a large number of globular proteins.

Folding pathways of a knotted protein with a realistic atomistic force field.

We report on atomistic simulation of the folding of a natively-knotted protein, MJ0366, based on a realistic force field. To the best of our knowledge this is the first reported effort where a realistic force field is used to investigate the folding pathways of a protein with complex native topology. By using the dominant-reaction pathway scheme we collected about 30 successful folding trajectories for the 82-amino acid long trefoil-knotted protein. Despite the dissimilarity of their initial unfolded configuration, these trajectories reach the natively-knotted state through a remarkably similar succession of steps. In particular it is found that knotting occurs essentially through a threading mechanism, involving the passage of the C-terminal through an open region created by the formation of the native [Formula: see text]-sheet at an earlier stage. The dominance of the knotting by threading mechanism is not observed in MJ0366 folding simulations using simplified, native-centric models. This points to a previously underappreciated role of concerted amino acid interactions, including non-native ones, in aiding the appropriate order of contact formation to achieve knotting.

Making lipids very unhappy to discover how they bind to proteins.

Membrane lipid composition is maintained by conserved lipid transfer proteins, but computational approaches to study their lipid-binding mechanisms are limiting. Srinivasan et al. (https://doi.org/10.1083/jcb.202312055) develop a clever molecular dynamics simulations assay to accurately model lipid-binding poses in lipid transfer proteins.

Sampling a Rare Protein Transition Using Quantum Annealing.

Simulating spontaneous structural rearrangements in macromolecules with classical molecular dynamics is an outstanding challenge. Conventional supercomputers can access time intervals of up to tens of µs, while many key events occur on exponentially longer time scales. Path sampling techniques have the advantage of focusing the computational power on barrier-crossing trajectories, but generating uncorrelated transition paths that explore diverse conformational regions remains a problem. We employ a hybrid path-sampling paradigm that addresses this issue by generating trial transition paths using a quantum annealing (QA) machine. We first employ a classical computer to perform an uncharted exploration of the conformational space. The data set generated in this exploration is then postprocessed using a path integral-based method to yield a coarse-grained network representation of the reactive kinetics. By resorting to a quantum annealer, quantum superposition can be exploited to encode all of the transition pathways in the initial quantum state, thus potentially solving the path exploration problem. Furthermore, each QA cycle yields a completely uncorrelated trial trajectory. We previously validated this scheme on a prototypically simple transition, which could be extensively characterized on a desktop computer. Here, we scale up in complexity and perform an all-atom simulation of a protein conformational transition that occurs on the millisecond time scale, obtaining results that match those of the Anton special-purpose supercomputer. Despite limitations due to the available quantum annealers, our study highlights how realistic biomolecular simulations provide potentially impactful new ground for applying, testing, and advancing quantum technologies.

ATG3 proteins possess a unique amphipathic α-helix essential for the Atg8/LC3 lipidation reaction.

In our recent paper, we uncovered that ATG3 exhibits a large degree of structural dynamics on autophagic membranes to efficiently carry out LC3 lipidation. ATG3 proteins possess an amphipathic α-helix (AH) identified by a small number of bulky and hydrophobic residues. This biophysical fingerprint allows for transient membrane association of ATG3 and facilitates its enzymatic reaction. This study will pave the way for a structural and mechanistic understanding of how membrane association of ATG proteins is orchestrated during autophagosome formation.

Membrane association of the ATG8 conjugation machinery emerges as a key regulatory feature for autophagosome biogenesis.

Autophagy is a highly conserved intracellular pathway that is essential for survival in all eukaryotes. In healthy cells, autophagy is used to remove damaged intracellular components, which can be as simple as unfolded proteins or as complex as whole mitochondria. Once the damaged component is captured, the autophagosome engulfs it and closes, isolating the content from the cytoplasm. The autophagosome then fuses with the late endosome and/or lysosome to deliver its content to the lysosome for degradation. Formation of the autophagosome, sequestration or capture of content, and closure all require the ATG proteins, which constitute the essential core autophagy protein machinery. This brief 'nutshell' will highlight recent data revealing the importance of small membrane-associated domains in the ATG proteins. In particular, recent findings from two parallel studies reveal the unexpected key role of α-helical structures in the ATG8 conjugation machinery and ATG8s. These studies illustrate how unique membrane association modules can control the formation of autophagosomes.

IGF2BP1 phosphorylation in the disordered linkers regulates ribonucleoprotein condensate formation and RNA metabolism.

The insulin-like growth factor 2 mRNA binding protein 1 (IGF2BP1) is a conserved RNA-binding protein that regulates RNA stability, localization and translation. IGF2BP1 is part of various ribonucleoprotein (RNP) condensates. However, the mechanism that regulates its assembly into condensates remains unknown. By using proteomics, we demonstrate that phosphorylation of IGF2BP1 at S181 in a disordered linker is regulated in a stress-dependent manner. Phosphomimetic mutations in two disordered linkers, S181E and Y396E, modulate RNP condensate formation by IGF2BP1 without impacting its binding affinity for RNA. Intriguingly, the S181E mutant, which lies in linker 1, impairs IGF2BP1 condensate formation in vitro and in cells, whereas a Y396E mutant in the second linker increases condensate size and dynamics. Structural approaches show that the first linker binds RNAs nonspecifically through its RGG/RG motif, an interaction weakened in the S181E mutant. Notably, linker 2 interacts with IGF2BP1's folded domains and these interactions are partially impaired in the Y396E mutant. Importantly, the phosphomimetic mutants impact IGF2BP1's interaction with RNAs and remodel the transcriptome in cells. Our data reveal how phosphorylation modulates low-affinity interaction networks in disordered linkers to regulate RNP condensate formation and RNA metabolism.

Molecular Free Energies, Rates, and Mechanisms from Data-Efficient Path Sampling Simulations.

Molecular dynamics is a powerful tool for studying the thermodynamics and kinetics of complex molecular events. However, these simulations can rarely sample the required time scales in practice. Transition path sampling overcomes this limitation by collecting unbiased trajectories and capturing the relevant events. Moreover, the integration of machine learning can boost the sampling while simultaneously learning a quantitative representation of the mechanism. Still, the resulting trajectories are by construction non-Boltzmann-distributed, preventing the calculation of free energies and rates. We developed an algorithm to approximate the equilibrium path ensemble from machine-learning-guided path sampling data. At the same time, our algorithm provides efficient sampling, mechanism, free energy, and rates of rare molecular events at a very moderate computational cost. We tested the method on the folding of the mini-protein chignolin. Our algorithm is straightforward and data-efficient, opening the door to applications in many challenging molecular systems.

Unique amphipathic α helix drives membrane insertion and enzymatic activity of ATG3.

Autophagosome biogenesis requires a localized perturbation of lipid membrane dynamics and a unique protein-lipid conjugate. Autophagy-related (ATG) proteins catalyze this biogenesis on cellular membranes, but the underlying molecular mechanism remains unclear. Focusing on the final step of the protein-lipid conjugation reaction, the ATG8/LC3 lipidation, we show how the membrane association of the conjugation machinery is organized and fine-tuned at the atomistic level. Amphipathic α helices in ATG3 proteins (AHATG3) have low hydrophobicity and contain less bulky residues. Molecular dynamics simulations reveal that AHATG3 regulates the dynamics and accessibility of the thioester bond of the ATG3~LC3 conjugate to lipids, enabling the covalent lipidation of LC3. Live-cell imaging shows that the transient membrane association of ATG3 with autophagic membranes is governed by the less bulky-hydrophobic feature of AHATG3. The unique properties of AHATG3 facilitate protein-lipid bilayer association, leading to the remodeling of the lipid bilayer required for the formation of autophagosomes.