Living cells constantly remodel the shape of their lipid membranes. In the endo-plasmic reticulum (ER), the reticulon homology domain (RHD) of the reticulophagy regulator 1 (RETR1/FAM134B) forms dense autophagic puncta that are associated with membrane removal by ER-phagy. In molecular dynamics (MD) simulations, we find that FAM134B-RHD spontaneously forms clusters, driven in part by curvature-mediated attraction. At a critical size, the FAM134B-RHD clusters induce the formation of membrane buds. The kinetics of budding depends sensitively on protein concentration and bilayer asymmetry. Our MD simulations shed light on the role of FAM134B-RHD in ER-phagy and show that membrane asymmetry can be used to modulate the kinetics barrier for membrane remodeling.
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.
Living cells constantly remodel the shape of their lipid membranes. In the endoplasmic reticulum (ER), the reticulon homology domain (RHD) of the reticulophagy regulator 1 (RETR1/FAM134B) forms dense autophagic puncta that are associated with membrane removal by ER-phagy. In molecular dynamics (MD) simulations, we find that FAM134B-RHD spontaneously forms clusters, driven in part by curvature-mediated attractions. At a critical size, as in a nucleation process, the FAM134B-RHD clusters induce the formation of membrane buds. The kinetics of budding depends sensitively on protein concentration and bilayer asymmetry. Our MD simulations shed light on the role of FAM134B-RHD in ER-phagy and show that membrane asymmetry can be used to modulate the kinetic barrier for membrane remodeling.
Nuclear pore complexes (NPCs) mediate nucleocytoplasmic transport. Their intricate 120 MDa architecture remains incompletely understood. Here, we report a near-complete structural model of the human NPC scaffold with explicit membrane and in multiple conformational states. We combined AI-based structure prediction with in situ and in cellulo cryo-electron tomography and integrative modeling. We show that linker Nups spatially organize the scaffold within and across subcomplexes to establish the higher-order structure. Microsecond-long molecular dynamics simulations suggest that the scaffold is not required to stabilize the inner and outer nuclear membrane fusion, but rather widens the central pore. Our work exemplifies how AI-based modeling can be integrated with in situ structural biology to understand subcellular architecture across spatial organization levels.
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 remodeling processes at atomic resolution are extremely difficult, due to their size, complexity, and the large times-scales on which these processes occur. Instead, elastic membrane models are used to simulate membrane shapes and transitions between them and to infer their properties and functions. Unfortunately, an 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 large and highly curved membrane structures. For validation, we reproduce phase diagrams of vesicles with varying surface-to-volume ratios and area difference. We also compute the density of states to verify correct Boltzmann sampling. The software can be used to tackle a range of large-scale membrane remodeling processes as a step toward cell-scale simulations. Additionally, extensive documentation make the software accessible to the broad biophysics and computational cell biology communities.