Refine
Document Type
- Preprint (25) (remove)
Language
- English (25)
Has Fulltext
- yes (25)
Is part of the Bibliography
- no (25)
TriMem: a parallelized hybrid Monte Carlo software for efficient simulations of lipid membranes
(2022)
Lipid membranes are integral building blocks of living cells and perform a multitude of biological functions. Currently, molecular simulations of cellular-scale membrane structures at atomic resolution are nearly impossible, due to their size, complexity, and the large times-scales required. Instead, elastic membrane models are used to simulate membrane topologies and transitions between them, and to infer their properties and functions. Unfortunately, efficiently parallelized open-source simulation code to do so has been lacking. Here, we present TriMem, a parallel hybrid Monte Carlo simulation engine for triangulated lipid membranes. The kernels are efficiently coded in C++ and wrapped with Python for ease-of-use. The parallel implementation of the energy and gradient calculations and of Monte Carlo flip moves of edges in the triangulated membrane enable us to simulate also large and highly curved sub-cellular structures. For validation, we reproduce phase diagrams of vesicles with varying surface-to-volume ratios and area difference. The software can tackle a range of membrane remodelling processes on sub-cellular and cellular scales. Additionally, extensive documentation make the software accessible to the broad biophysics and computational cell biology communities.
Cryo-electron tomography (cryo-ET) is a powerful method to elucidate subcellular architecture and to structurally analyse biomolecules in situ by subtomogram averaging (STA). Specimen thickness is a key factor affecting cryo-ET data quality. Cells that are too thick for transmission imaging can be thinned by cryo-focused-ion-beam (cryo-FIB) milling. However, optimal specimen thickness for cryo-ET on lamellae has not been systematically investigated. Furthermore, the ions used to ablate material can cause damage in the lamellae, thereby reducing STA resolution. Here, we systematically benchmark the resolution depending on lamella thickness and the depth of the particles within the sample. Up to ca. 180 nm, lamella thickness does not negatively impact resolution. This shows that there is no need to generate very thin lamellae and thickness can be chosen such that it captures major cellular features. Furthermore, we show that gallium-ion-induced damage extends to depths of up to 30 nm from either lamella surface.
A key event in cellular physiology is the decision between membrane biogenesis and fat storage. Phosphatidic acid (PA) is an important lipid intermediate and signaling lipid at the branch point of these pathways and constantly monitored by the transcriptional repressor Opi1 to orchestrate lipid metabolism. Here, we report on the mechanism of membrane recognition by Opi1 and identify an amphipathic helix (AH) for the selective binding to membranes containing PA over phosphatidylserine (PS). The insertion of the AH into the hydrophobic core of the membrane renders Opi1 sensitive to the lipid acyl chain composition as an important factor contributing to the regulation of membrane biogenesis. Based on these findings, we rationally designed the membrane binding properties of Opi1 to control its responsiveness in the physiological context. Using extensive molecular dynamics (MD) simulations, we identified two PA-selective three-finger grips that tightly bind the phosphate headgroup, while interacting less intimately and more transiently with PS. This work establishes lipid headgroup selectivity as a new feature in the family of AH-containing membrane property sensors.
Autophagy is a highly conserved catabolic process through which defective or otherwise harmful cellular components are targeted for degradation via the lysosomal route. Regulatory pathways, involving post-translational modifications such as phosphorylation, play a critical role in controlling this tightly orchestrated process. Here, we demonstrate that TBK1 regulates autophagy by phosphorylating autophagy modifiers LC3C and GABARAP-L2 on surface-exposed serine residues (LC3C S93 and S96; GABARAP-L2 S87 and S88). This phosphorylation event impedes their binding to the processing enzyme ATG4 by destabilizing the complex. Phosphorylated LC3C/GABARAP-L2 cannot be removed from liposomes by ATG4 and are thus protected from ATG4-mediated premature removal from nascent autoph-agosomes. This ensures a steady coat of lipidated LC3C/GABARAP-L2 throughout the early steps in autophagosome formation and aids in maintaining a unidirectional flow of the autophagosome to the lysosome. Taken together, we present a new regulatory mechanism of autophagy, which influences the conjugation and de-conjugation of LC3C and GABARAP-L2 to autophagosomes by TBK1-mediated phosphorylation.
Gasdermin-D (GSDMD) is the ultimate effector of pyroptosis, a form of programmed cell death associated with pathogen invasion and inflammation. After proteolytic cleavage by caspases activated by the inflammasome, the GSDMD N-terminal domain (GSDMDNT) assembles on the inner leaflet of the plasma membrane and induces the formation of large membrane pores. We use atomistic molecular dynamics simulations to study GSDMDNT monomers, oligomers, and rings in an asymmetric plasma membrane mimetic. We identify distinct interaction motifs of GSDMDNT with phosphatidylinositol-4,5-bisphosphate (PI(4,5)P2) and phosphatidylserine (PS) head-groups and describe differential lipid binding between the pore and prepore conformations. Oligomers are stabilized by shared lipid binding sites between neighboring monomers acting akin to double-sided tape. We show that already small GSDMDNT oligomers form stable, water-filled and ion-conducting membrane pores bounded by curled beta-sheets. In large-scale simulations, we resolve the process of pore formation by lipid detachment from GSDMDNT arcs and lipid efflux from partial rings. We find that that high-order GSDMDNT oligomers can crack under the line tension of 86 pN created by an open membrane edge to form the slit pores or closed GSDMDNT rings seen in experiment. Our simulations provide a detailed view of key steps in GSDMDNT-induced plasma membrane pore formation, including sublytic pores that explain nonselective ion flux during early pyroptosis.
The SLC26 family of transporters maintains anion equilibria in all kingdoms of life. The family shares a 7 + 7 transmembrane segments inverted repeat architecture with the SLC4 and SLC23 families, but holds a regulatory STAS domain in addition. While the only experimental SLC26 structure is monomeric, SLC26 proteins form structural and functional dimers in the lipid membrane. Here we resolve the structure of an SLC26 dimer embedded in a lipid membrane and characterize its functional relevance by combining PELDOR distance measurements and biochemical studies with MD simulations and spin-label ensemble refinement. Our structural model reveals a unique interface different from the SLC4 and SLC23 families. The functionally relevant STAS domain exerts a stabilizing effect on regions central in this dimer. Characterization of heterodimers indicates that protomers in the dimer functionally interact. The combined structural and functional data define the framework for a mechanistic understanding of functional cooperativity in SLC26 dimers.
Human feline leukaemia virus subgroup C receptor-related proteins 1 and 2 (FLVCR1 and 2) are major facilitator superfamily transporters from the solute carrier family 49. Dysregulation of these ubiquitous transporters has been linked to various haematological and neurological disorders. While both FLVCRs were initially proposed to hold a physiological function in heme transport, subsequent studies questioned this notion. Here, we used structural, computational and biochemical methods and conclude that these two FLVCRs function as human choline transporters. We present cryo-electron microscopy structures of FLVCRs in different inward- and outward-facing conformations, captured in the apo state or in complex with choline in their translocation pathways. Our findings provide insights into the molecular framework of choline coordination and transport, largely mediated by conserved cation-π interactions, and further illuminate the conformational dynamics of the transport cycle. Moreover, we identified a heme binding site on the protein surface of the FLVCR2 N-domain, and observed that heme actively drives the conformational transitions of the protein. This auxiliary binding site might indicate a potential regulatory role of heme in the FLVCR2 transport mechanisms. Our work resolves the contested substrate specificity of the FLVCRs, and sheds light on the process of maintaining cellular choline homeostasis at the molecular level.
Disordered proteins and nucleic acids can condense into droplets that resemble the membraneless organelles observed in living cells. MD simulations offer a unique tool to characterize the molecular interactions governing the formation of these biomolecular condensates, their physico-chemical properties, and the factors controlling their composition and size. However, biopolymer condensation depends sensitively on the balance between different energetic and entropic contributions. Here, we develop a general strategy to fine-tune the potential energy function for molecular dynamics simulations of biopolymer phase separation. We rebalance protein-protein interactions against solvation and entropic contributions to match the excess free energy of transferring proteins between dilute solution and condensate. We illustrate this formalism by simulating liquid droplet formation of the FUS low complexity domain (LCD) with a rebalanced MARTINI model. By scaling the strength of the nonbonded interactions in the coarse-grained MARTINI potential energy function, we map out a phase diagram in the plane of protein concentration and interaction strength. Above a critical scaling factor of αc ≈ 0.6, FUS LCD condensation is observed, where α = 1 and 0 correspond to full and repulsive interactions in the MARTINI model, respectively. For a scaling factor α = 0.65, we recover the experimental densities of the dilute and dense phases, and thus the excess protein transfer free energy into the droplet and the saturation concentration where FUS LCD condenses. In the region of phase separation, we simulate FUS LCD droplets of four different sizes in stable equilibrium with the dilute phase and slabs of condensed FUS LCD for tens of microseconds, and over one millisecond in aggregate. We determine surface tensions in the range of 0.01 to 0.4mN/m from the fluctuations of the droplet shape and from the capillary-wave-like broadening of the interface between the two phases. From the dynamics of the protein end-to-end distance, we estimate shear viscosities from 0.001 to 0.02Pas for the FUS LCD droplets with scaling factors α in the range of 0.625 to 0.75, where we observe liquid droplets. Significant hydration of the interior of the droplets keeps the proteins mobile and the droplets fluid.
The spike protein of SARS-CoV-2 is a highly flexible membrane receptor that triggers the translocation of the virus into cells by attaching to the human receptors. Like other type I membrane receptors, this protein has several extracellular domains connected by flexible hinges. The presence of these hinges results in high flexibility, which consequently results in challenges in defining the conformation of the protein. Here, We developed a new method to define the conformational space based on a few variables inspired by the robotic field’s methods to determine a robotic arm’s forward kinematics. Using newly performed atomistic molecular dynamics (MD) simulations and publicly available data, we found that the Denavit-Hartenberg (DH) parameters can reliably show the changes in the local conformation. Furthermore, the rotational and translational components of the homogenous transformation matrix constructed based on the DH parameters can identify the changes in the global conformation of the spike and also differentiate between the conformation with a similar position of the spike head, which other types of parameters, such as spherical coordinates, fail to distinguish between such conformations. Finally, the new method will be beneficial for looking at the conformational heterogeneity in all other type I membrane receptors.
Upon infection, human immunodeficiency virus (HIV-1) releases its cone-shaped capsid into the cytoplasm of infected T-cells and macrophages. As its largest known cargo, the capsid enters the nuclear pore complex (NPC), driven by interactions with numerous FG-repeat nucleoporins (FG-Nups). Whether NPCs structurally adapt to capsid passage and whether capsids are modified during passage remains unknown, however. Here, we combined super-resolution and correlative microscopy with cryo electron tomography and molecular simulations to study nuclear entry of HIV-1 capsids in primary human macrophages. We found that cytosolically bound cyclophilin A is stripped off capsids entering the NPC, and the capsid hexagonal lattice remains largely intact inside and beyond the central channel. Strikingly, the NPC scaffold rings frequently crack during capsid passage, consistent with computer simulations indicating the need for NPC widening. The unique cone shape of the HIV-1 capsid facilitates its entry into NPCs and helps to crack their rings.