Refine
Year of publication
Language
- English (59)
Has Fulltext
- yes (59)
Is part of the Bibliography
- no (59) (remove)
Keywords
Institute
- Physik (59) (remove)
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.
We present a method that enables the identification and analysis of conformational Markovian transition states from atomistic or coarse-grained molecular dynamics (MD) trajectories. Our algorithm is presented by using both analytical models and examples from MD simulations of the benchmark system helix-forming peptide Ala5, and of larger, biomedically important systems: the 15-lipoxygenase-2 enzyme (15-LOX-2), the epidermal growth factor receptor (EGFR) protein, and the Mga2 fungal transcription factor. The analysis of 15-LOX-2 uses data generated exclusively from biased umbrella sampling simulations carried out at the hybrid ab initio density functional theory (DFT) quantum mechanics/molecular mechanics (QM/MM) level of theory. In all cases, our method automatically identifies the corresponding transition states and metastable conformations in a variationally optimal way, with the input of a set of relevant coordinates, by accurately reproducing the intrinsic slowest relaxation rate of each system. Our approach offers a general yet easy-to-implement analysis method that provides unique insight into the molecular mechanism and the rare but crucial (i.e., rate-limiting) transition states occurring along conformational transition paths in complex dynamical systems such as molecular trajectories.
The ATP-binding cassette transporter TAPL translocates polypeptides from the cytosol into the lysosomal lumen. TAPL can be divided into two functional units: coreTAPL, active in ATP-dependent peptide translocation, and the N-terminal membrane spanning domain, TMD0, responsible for cellular localization and interaction with the lysosomal associated membrane proteins LAMP-1 and LAMP-2. Although the structure and function of ABC transporters were intensively analyzed in the past, the knowledge about accessory membrane embedded domains is limited. Therefore, we expressed the TMD0 of TAPL via a cell-free expression system and confirmed its correct folding by NMR and interaction studies. In cell as well as cell-free expressed TMD0 forms oligomers, which were assigned as dimers by PELDOR spectroscopy and static light scattering. By NMR spectroscopy of uniformly and selectively isotope labeled TMD0 we performed a complete backbone and partial side chain assignment. Accordingly, TMD0 has a four transmembrane helix topology with a short helical segment in a lysosomal loop. The topology of TMD0 was confirmed by paramagnetic relaxation enhancement with paramagnetic stearic acid as well as by nuclear Overhauser effects with c6-DHPC and cross-peaks with water.
Autophagy is a physiological process for the recycling and degradation of cellular materials. Forming the autophagosome from the phagophore, a cup-shaped double-membrane vesicle, is a critical step in autophagy. The origin of the cup shape of the phagophore is poorly understood. In yeast, fusion of a small number of Atg9-containing vesicles is considered a key step in autophagosome biogenesis, aided by Atg1 complexes (ULK1 in mammals) localized at the preautophagosomal structure (PAS). In particular, the S-shaped Atg17-Atg31-Atg29 subcomplex of Atg1 is critical for phagophore nucleation at the PAS. To study this process, we simulated membrane remodeling processes in the presence and absence of membrane associated Atg17. We show that at least three vesicles need to fuse to induce the phagophore shape, consistent with experimental observations. However, fusion alone is not sufficient. Interactions with 34-nm long, S-shaped Atg17 complexes are required to overcome a substantial kinetic barrier in the transition to the cup-shaped phagophore. Our finding rationalizes the recruitment of Atg17 complexes to the yeast PAS, and their unusual shape. In control simulations without Atg17, with weakly binding Atg17, or with straight instead of S-shaped Atg17, the membrane shape transition did not occur. We confirm the critical role of Atg17-membrane interactions experimentally by showing that mutations of putative membrane interaction sites result in reduction or loss of autophagic activity in yeast. Fusion of a small number of vesicles followed by Atg17-guided membrane shape-remodeling thus emerges as a viable route to phagophore formation.
Cysteine cross-linking in native membranes establishes the transmembrane architecture of Ire1
(2021)
The ER is a key organelle of membrane biogenesis and crucial for the folding of both membrane and secretory proteins. Sensors of the unfolded protein response (UPR) monitor the unfolded protein load in the ER and convey effector functions for maintaining ER homeostasis. Aberrant compositions of the ER membrane, referred to as lipid bilayer stress, are equally potent activators of the UPR. How the distinct signals from lipid bilayer stress and unfolded proteins are processed by the conserved UPR transducer Ire1 remains unknown. Here, we have generated a functional, cysteine-less variant of Ire1 and performed systematic cysteine cross-linking experiments in native membranes to establish its transmembrane architecture in signaling-active clusters. We show that the transmembrane helices of two neighboring Ire1 molecules adopt an X-shaped configuration independent of the primary cause for ER stress. This suggests that different forms of stress converge in a common, signaling-active transmembrane architecture of Ire1.
Complex I couples the free energy released from quinone (Q) reduction to pump protons across the biological membrane in the respiratory chains of mitochondria and many bacteria. The Q reduction site is separated by a large distance from the proton-pumping membrane domain. To address the molecular mechanism of this long-range proton-electron coupling, we perform here full atomistic molecular dynamics simulations, free energy calculations, and continuum electrostatics calculations on complex I from Thermus thermophilus. We show that the dynamics of Q is redox-state-dependent, and that quinol, QH2, moves out of its reduction site and into a site in the Q tunnel that is occupied by a Q analog in a crystal structure of Yarrowia lipolytica. We also identify a second Q-binding site near the opening of the Q tunnel in the membrane domain, where the Q headgroup forms strong interactions with a cluster of aromatic and charged residues, while the Q tail resides in the lipid membrane. We estimate the effective diffusion coefficient of Q in the tunnel, and in turn the characteristic time for Q to reach the active site and for QH2 to escape to the membrane. Our simulations show that Q moves along the Q tunnel in a redox-state-dependent manner, with distinct binding sites formed by conserved residue clusters. The motion of Q to these binding sites is proposed to be coupled to the proton-pumping machinery in complex I.
The structure and flexibility of RNA depends sensitively on the microenvironment. Using pulsed electron-electron double-resonance (PELDOR)/double electron-electron resonance (DEER) spectroscopy combined with advanced labeling techniques, we show that the structure of double-stranded RNA (dsRNA) changes upon internalization into Xenopus lævis oocytes. Compared to dilute solution, the dsRNA A-helix is more compact in cells. We recapitulate this compaction in a densely crowded protein solution. Atomic-resolution molecular dynamics simulations of dsRNA semi-quantitatively capture the compaction, and identify non-specific electrostatic interactions between proteins and dsRNA as a possible driver of this effect.
During infection the SARS-CoV-2 virus fuses its viral envelope with cellular membranes of its human host. The viral spike (S) protein mediates both the initial contact with the host cell and the subsequent membrane fusion. Proteolytic cleavage of S at the S2′ site exposes its fusion peptide (FP) as the new N-terminus. By binding to the host membrane, the FP anchors the virus to the host cell. The reorganization of S2 between virus and host then pulls the two membranes together. Here we use molecular dynamics (MD) simulations to study the two core functions of the SARS-CoV-2 FP: to attach quickly to cellular membranes and to form an anchor strong enough to withstand the mechanical force during membrane fusion. In eight 10 μs long MD simulations of FP in proximity to endosomal and plasma membranes, we find that FP binds spontaneously to the membranes and that binding proceeds predominantly by insertion of two short amphipathic helices into the membrane interface. Connected via a flexible linker, the two helices can bind the membrane independently, yet binding of one promotes the binding of the other by tethering it close to the target membrane. By simulating mechanical pulling forces acting on the C-terminus of the FP, we then show that the bound FP can bear forces up to 250 pN before detaching from the membrane. This detachment force is more than 10-fold higher than an estimate of the force required to pull host and viral membranes together for fusion. We identify a fully conserved disulfide bridge in the FP as a major factor for the high mechanical stability of the FP membrane anchor. We conclude, first, that the sequential binding of two short amphipathic helices allows the SARS-CoV-2 FP to insert quickly into the target membrane, before the virion is swept away after shedding the S1 domain connecting it to the host cell receptor. Second, we conclude that the double attachment and the conserved disulfide bridge establish the strong anchoring required for subsequent membrane fusion. Multiple distinct membrane-anchoring elements ensure high avidity and high mechanical strength of FP–membrane binding.
Proton-powered c-ring rotation in mitochondrial ATP synthase is crucial to convert the transmembrane protonmotive force into torque to drive the synthesis of ATP. Capitalizing on recent cryo-EM structures, we aim at a structural and energetic understanding of how functional directional rotation is achieved. We performed multi-microsecond atomistic simulations to determine the free energy profiles along the c-ring rotation angle before and after the arrival of a new proton. Our results reveal that rotation proceeds by dynamic sliding of the ring over the a-subunit surface, during which interactions with conserved polar residues stabilize distinct intermediates. Ordered water chains line up for a Grotthuss-type proton transfer in one of these intermediates. After proton transfer, a high barrier prevents backward rotation and an overall drop in free energy favors forward rotation, ensuring the directionality of c-ring rotation required for the thermodynamically disfavored ATP synthesis. The essential arginine of the a-subunit stabilizes the rotated configuration through a salt-bridge with the c-ring. Overall, we describe a complete mechanism for the rotation step of the ATP synthase rotor, thereby illuminating a process critical to all life at atomic resolution.
Proton-powered c-ring rotation in mitochondrial ATP synthase is crucial to convert the transmembrane protonmotive force into torque to drive the synthesis of ATP. Capitalizing on recent cryo-EM structures, we aim at a structural and energetic understanding of how functional directional rotation is achieved. We performed multi-microsecond atomistic simulations to determine the free energy profiles along the c-ring rotation angle before and after the arrival of a new proton. Our results reveal that rotation proceeds by dynamic sliding of the ring over the a-subunit surface, during which interactions with conserved polar residues stabilize distinct intermediates. Ordered water chains line up for a Grotthuss-type proton transfer in one of these intermediates. After proton transfer, a high barrier prevents backward rotation and an overall drop in free energy favors forward rotation, ensuring the directionality of c-ring rotation required for the thermodynamically disfavored ATP synthesis. The essential arginine of the a-subunit stabilizes the rotated configuration through a salt-bridge with the c-ring. Overall, we describe a complete mechanism for the rotation step of the ATP synthase rotor, thereby illuminating a process critical to all life at atomic resolution.
Cryo-electron tomography (CryoET) resolves individual macromolecules inside living cells. However, the complex composition and high density of cells challenge the faithful identification of features in tomograms. Here, we capitalize on recent advances in electron tomography and demonstrate that 3D template matching (TM) localizes a wide range of structures inside crowded eukaryotic cells with confidence 10 to 100-fold above the noise level. We establish a TM pipeline with systematically tuned parameters for automated, objective and comprehensive feature identification. High-fidelity and high-confidence localizations of nuclear pore complexes, vaults, ribosomes, proteasomes, lipid membranes and microtubules, and individual subunits, demonstrate that TM is generic. We resolve ~100-kDa proteins, connect the functional states of complexes to their cellular localization, and capture vaults carrying ribosomal cargo in situ. By capturing individual molecular events inside living cells with defined statistical confidence, high-confidence TM greatly speeds up the CryoET workflow and sets the stage for visual proteomics.
Molecular mechanisms of inorganic-phosphate release from the core and barbed end of actin filaments
(2023)
The release of inorganic phosphate (Pi) from actin filaments constitutes a key step in their regulated turnover, which is fundamental to many cellular functions. However, the molecular mechanisms underlying Pi release from both the core and barbed end of actin filaments remain unclear. Here, we combine cryo-EM with molecular dynamics simulations and in vitro reconstitution to demonstrate how actin releases Pi through a ‘molecular backdoor’. While constantly open at the barbed end, the backdoor is predominantly closed in filament-core subunits and only opens transiently through concerted backbone movements and rotameric rearrangements of residues close to the nucleotide binding pocket. This mechanism explains why Pi escapes rapidly from the filament end and yet slowly from internal actin subunits. In an actin variant associated with nemaline myopathy, the backdoor is predominantly open in filament-core subunits, resulting in greatly accelerated Pi release after polymerization and filaments with drastically shortened ADP-Pi caps. This demonstrates that the Pi release rate from F-actin is controlled by steric hindrance through the backdoor rather than by the disruption of the ionic bond between Pi and Mg2+ at the nucleotide-binding site. Our results provide the molecular basis for Pi release from actin and exemplify how a single, disease-linked point mutation distorts the nucleotide state distribution and atomic structure of the actin filament.
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.
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.
Intrinsically disordered regions (IDRs) are essential for membrane receptor regulation but often remain unresolved in structural studies. TRPV4, a member of the TRP vanilloid channel family involved in thermo- and osmosensation, has a large N-terminal IDR of approximately 150 amino acids. With an integrated structural biology approach, we analyze the structural ensemble of the TRPV4 IDR and identify a network of regulatory elements that modulate channel activity in a hierarchical lipid-dependent manner through transient long-range interactions. A highly conserved autoinhibitory patch acts as a master regulator by competing with PIP2 binding to attenuate channel activity. Molecular dynamics simulations show that loss of the interaction between the PIP2-binding site and the membrane reduces the force exerted by the IDR on the structured core of TRPV4. This work demonstrates that IDR structural dynamics are coupled to TRPV4 activity and highlights the importance of IDRs for TRP channel function and regulation.
The spike protein (S) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is required for cell entry and is the primary focus for vaccine development. In this study, we combined cryo–electron tomography, subtomogram averaging, and molecular dynamics simulations to structurally analyze S in situ. Compared with the recombinant S, the viral S was more heavily glycosylated and occurred mostly in the closed prefusion conformation. We show that the stalk domain of S contains three hinges, giving the head unexpected orientational freedom. We propose that the hinges allow S to scan the host cell surface, shielded from antibodies by an extensive glycan coat. The structure of native S contributes to our understanding of SARS-CoV-2 infection and potentially to the development of safe vaccines.
A key event in cellular physiology is the decision between membrane biogenesis and fat storage. Phosphatidic acid (PA) is an important intermediate at the branch point of these pathways and is continuously monitored by the transcriptional repressor Opi1 to orchestrate lipid metabolism. In this study, we report on the mechanism of membrane recognition by Opi1 and identify an amphipathic helix (AH) for selective binding of PA over phosphatidylserine (PS). The insertion of the AH into the membrane core renders Opi1 sensitive to the lipid acyl chain composition and provides a means to adjust membrane biogenesis. By rational design of the AH, we tune the membrane-binding properties of Opi1 and control its responsiveness in vivo. Using extensive molecular dynamics simulations, we identify two PA-selective three-finger grips that tightly bind the PA phosphate headgroup while interacting less intimately with PS. This work establishes lipid headgroup selectivity as a new feature in the family of AH-containing membrane property sensors.
The interaction between the Heat Shock Proteins 70 and 40 is at the core of the ATPase regulation of the chaperone machinery that maintains protein homeostasis. However, the structural details of the interaction remain elusive and contrasting models have been proposed for the transient Hsp70/Hsp40 complexes. Here we combine molecular simulations based on both coarse-grained and atomistic models with coevolutionary sequence analysis to shed light on this problem by focusing on the bacterial DnaK/DnaJ system. The integration of these complementary approaches resulted in a novel structural model that rationalizes previous experimental observations. We identify an evolutionarily conserved interaction surface formed by helix II of the DnaJ J-domain and a structurally contiguous region of DnaK, involving lobe IIA of the nucleotide binding domain, the inter-domain linker, and the β-basket of the substrate binding domain.
The interaction between the Heat Shock Proteins 70 and 40 is at the core of the ATPase regulation of the chaperone machinery that maintains protein homeostasis. However, the structural details of this fundamental interaction are still elusive and contrasting models have been proposed for the transient Hsp70/Hsp40 complexes. Here we combine molecular simulations based on both coarsegrained and atomistic models with co-evolutionary sequence analysis to shed light on this problem by focusing on the bacterial DnaK/DnaJ system. The integration of these complementary approaches resulted into a novel structural model that rationalizes previous experimental observations. We identify an evolutionary-conserved interaction surface formed by helix II of the DnaJ J-domain and a groove on lobe IIA of the DnaK nucleotide binding domain, involving the inter-domain linker.
Famotidine inhibits toll-like receptor 3-mediated inflammatory signaling in SARS-CoV-2 infection
(2021)
Apart from prevention using vaccinations, the management options for COVID-19 remain limited. In retrospective cohort studies, use of famotidine, a specific oral H2 receptor antagonist (antihistamine), has been associated with reduced risk of intubation and death in patients hospitalized with COVID-19. In a case series, nonhospitalized patients with COVID-19 experienced rapid symptom resolution after taking famotidine, but the molecular basis of these observations remains elusive. Here we show using biochemical, cellular, and functional assays that famotidine has no effect on viral replication or viral protease activity. However, famotidine can affect histamine-induced signaling processes in infected Caco2 cells. Specifically, famotidine treatment inhibits histamine-induced expression of Toll-like receptor 3 (TLR3) in SARS-CoV-2 infected cells and can reduce TLR3-dependent signaling processes that culminate in activation of IRF3 and the NF-κB pathway, subsequently controlling antiviral and inflammatory responses. SARS-CoV-2-infected cells treated with famotidine demonstrate reduced expression levels of the inflammatory mediators CCL-2 and IL6, drivers of the cytokine release syndrome that precipitates poor outcome for patients with COVID-19. Given that pharmacokinetic studies indicate that famotidine can reach concentrations in blood that suffice to antagonize histamine H2 receptors expressed in mast cells, neutrophils, and eosinophils, these observations explain how famotidine may contribute to the reduced histamine-induced inflammation and cytokine release, thereby improving the outcome for patients with COVID-19.
Secretins form multimeric channels across the outer membrane of Gram-negative bacteria that mediate the import or export of substrates and/or extrusion of type IV pili. The secretin complex of Thermus thermophilus is an oligomer of the 757-residue PilQ protein, essential for DNA uptake and pilus extrusion. Here, we present the cryo-EM structure of this bifunctional complex at a resolution of ~7 Å using a new reconstruction protocol. Thirteen protomers form a large periplasmic domain of six stacked rings and a secretin domain in the outer membrane. A homology model of the PilQ protein was fitted into the cryo-EM map. A crown-like structure outside the outer membrane capping the secretin was found not to be part of PilQ. Mutations in the secretin domain disrupted the crown and abolished DNA uptake, suggesting a central role of the crown in natural transformation.
Mechanism of the electroneutral sodium/proton antiporter PaNhaP from transition-path shooting
(2019)
Na+/H+ antiporters exchange sodium ions and protons on opposite sides of lipid membranes. The electroneutral Na+/H+ antiporter NhaP from archaea Pyrococcus abyssi (PaNhaP) is a functional homolog of the human Na+/H+ exchanger NHE1, which is an important drug target. Here we resolve the Na+ and H+ transport cycle of PaNhaP by transition-path sampling. The resulting molecular dynamics trajectories of repeated ion transport events proceed without bias force, and overcome the enormous time-scale gap between seconds-scale ion exchange and microseconds simulations. The simulations reveal a hydrophobic gate to the extracellular side that opens and closes in response to the transporter domain motion. Weakening the gate by mutagenesis makes the transporter faster, suggesting that the gate balances competing demands of fidelity and efficiency. Transition-path sampling and a committor-based reaction coordinate optimization identify the essential motions and interactions that realize conformational alternation between the two access states in transporter function.
Abstract
The endoplasmic reticulum (ER) is a key organelle of membrane biogenesis and crucial for the folding of both membrane and secretory proteins. Sensors of the unfolded protein response (UPR) monitor the unfolded protein load in the ER and convey effector functions for maintaining ER homeostasis. Aberrant compositions of the ER membrane, referred to as lipid bilayer stress, are equally potent activators of the UPR. How the distinct signals from lipid bilayer stress and unfolded proteins are processed by the conserved UPR transducer Ire1 remains unknown. Here, we have generated a functional, cysteine-less variant of Ire1 and performed systematic cysteine crosslinking experiments in native membranes to establish its transmembrane architecture in signaling-active clusters. We show that the transmembrane helices of two neighboring Ire1 molecules adopt an X-shaped configuration independent of the primary cause for ER stress. This suggests that different forms of stress converge in a common, signaling-active transmembrane architecture of Ire1.
Summary
The endoplasmic reticulum (ER) is a hotspot of lipid biosynthesis and crucial for the folding of membrane and secretory proteins. The unfolded protein response (UPR) controls the size and folding capacity of the ER. The conserved UPR transducer Ire1 senses both unfolded proteins and aberrant lipid compositions to mount adaptive responses. Using a biochemical assay to study Ire1 in signaling-active clusters, Väth et al. provide evidence that the neighboring transmembrane helices of clustered Ire1 form an ‘X’ irrespectively of the primary cause of ER stress. Hence, different forms of ER stress converge in a common, signaling-active transmembrane architecture of Ire1.
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.
Maximum likelihood estimates of diffusion coefficients from single-particle tracking experiments
(2021)
Single-molecule localization microscopy allows practitioners to locate and track labeled molecules in biological systems. When extracting diffusion coefficients from the resulting trajectories, it is common practice to perform a linear fit on mean-squared-displacement curves. However, this strategy is suboptimal and prone to errors. Recently, it was shown that the increments between the observed positions provide a good estimate for the diffusion coefficient, and their statistics are well-suited for likelihood-based analysis methods. Here, we revisit the problem of extracting diffusion coefficients from single-particle tracking experiments subject to static noise and dynamic motion blur using the principle of maximum likelihood. Taking advantage of an efficient real-space formulation, we extend the model to mixtures of subpopulations differing in their diffusion coefficients, which we estimate with the help of the expectation–maximization algorithm. This formulation naturally leads to a probabilistic assignment of trajectories to subpopulations. We employ the theory to analyze experimental tracking data that cannot be explained with a single diffusion coefficient. We test how well a dataset conforms to the assumptions of a diffusion model and determine the optimal number of subpopulations with the help of a quality factor of known analytical distribution. To facilitate use by practitioners, we provide a fast open-source implementation of the theory for the efficient analysis of multiple trajectories in arbitrary dimensions simultaneously.
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.
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.
More than 75% of surface and secreted proteins are modified by covalent addition of complex sugars through N- and O-glycosylation. Unlike proteins, glycans do not typically adopt specific secondary structures and remain very mobile, influencing protein dynamics and interactions with other molecules. Glycan conformational freedom impairs complete structural elucidation of glycoproteins. Computer simulations may be used to model glycan structure and dynamics. However, such simulations typically require thousands of computing hours on specialized supercomputers, thus limiting routine use. Here, we describe a reductionist method that can be implemented on personal computers to graft ensembles of realistic glycan conformers onto static protein structures in a matter of minutes. Using this open-source pipeline, we reconstructed the full glycan cover of SARS-CoV-2 Spike protein (S-protein) and a human GABAA receptor. Focusing on S-protein, we show that GlycoSHIELD recapitulates key features of extended simulations of the glycosylated protein, including epitope masking, and provides new mechanistic insights on N-glycan impact on protein structural dynamics.
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 physicochemical 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. For a scaling factor α = 0.65, we recover 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–0.4 mN/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.02 Pa s for the FUS-LCD droplets with scaling factors α in the range of 0.625–0.75, where we observe liquid droplets. Significant hydration of the interior of the droplets keeps the proteins mobile and the droplets fluid.
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.
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.
The spike (S) protein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is required for cell entry and is the major focus for vaccine development. We combine cryo electron tomography, subtomogram averaging and molecular dynamics simulations to structurally analyze S in situ. Compared to recombinant S, the viral S is more heavily glycosylated and occurs predominantly in a closed pre-fusion conformation. We show that the stalk domain of S contains three hinges that give the globular domain unexpected orientational freedom. We propose that the hinges allow S to scan the host cell surface, shielded from antibodies by an extensive glycan coat. The structure of native S contributes to our understanding of SARS-CoV-2 infection and the development of safe vaccines. The large scale tomography data set of SARS-CoV-2 used for this study is therefore sufficient to resolve structural features to below 5 Ångstrom, and is publicly available at EMPIAR-10453.
The protein Atg2 has been proposed to form a membrane tether that mediates lipid transfer from the ER to the phagophore in autophagy. However, recent kinetic measurements on the human homolog ATG2A indicated a transport rate of only about one lipid per minute, which would be far too slow to deliver the millions of lipids required to form a phagophore on a physiological time scale. Here, we revisit the analysis of the fluorescence quenching experiments. We develop a detailed kinetic model of the lipid transfer between two membranes bridged by a tether that forms a conduit for lipids. The model provides an excellent fit to the fluorescence experiments, with a lipid transfer rate of about 100 per second and protein. At this rate, Atg2-mediated transfer can supply a significant fraction of the lipids required in autophagosome biogenesis. Our kinetic model is generally applicable to lipid-transfer experiments, in particular to proteins forming organelle contact sites in cells.
The severity of the COVID-19 pandemic, caused by the SARS-CoV-2 coronavirus, calls for the urgent development of a vaccine. The primary immunological target is the SARS-CoV-2 spike (S) protein. S is exposed on the viral surface to mediate viral entry into the host cell. To identify possible antibody binding sites not shielded by glycans, 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 vaccine development. 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 in epitope mapping.
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, the GSDMD N-terminal domain (GSDMDNT) assembles on the inner leaflet of the plasma membrane and induces the formation of 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) headgroups and describe their conformational dependence. Oligomers are stabilized by shared lipid binding sites between neighboring monomers acting akin to double-sided tape. We show that already small GSDMDNT oligomers support stable, water-filled, and ion-conducting membrane pores bounded by curled beta-sheets. In large-scale simulations, we resolve the process of pore formation from GSDMDNT arcs and lipid efflux from partial rings. We find 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 atomic force microscopy experiments. 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.
During infection the SARS-CoV-2 virus fuses its viral envelope with cellular membranes of its human host. Initial contact with the host cell and membrane fusion are both mediated by the viral spike (S) protein. Proteolytic cleavage of S at the S2′ site exposes its 40 amino acid long fusion peptide (FP). Binding of the FP to the host membrane anchors the S2 domain of S in both the viral and the host membrane. The reorganization of S2 then pulls the two membranes together. Here we use molecular dynamics (MD) simulations to study the two core functions of the SARS-CoV-2 FP: to attach quickly to cellular membranes and to form an anchor strong enough to withstand the mechanical force during membrane fusion. In eight 10 μs-long MD simulations of FP in proximity to endosomal and plasma membranes, we find that FP binds spontaneously to the membranes and that binding proceeds predominantly by insertion of two short amphipathic helices into the membrane interface. Connected via a flexible linker, the two helices can bind the membrane independently, yet binding of one promotes the binding of the other by tethering it close to the target membrane. By simulating mechanical pulling forces acting on the C-terminus of the FP we then show that the bound FP can bear forces up to 250 pN before detaching from the membrane. This detachment force is more than ten-fold higher than an estimate of the force required to pull host and viral membranes together for fusion. We identify a fully conserved disulfide bridge in the FP as a major factor for the high mechanical stability of the FP membrane anchor. We conclude, first, that the sequential binding of two short amphipathic helices allows the SARS-CoV-2 FP to insert quickly into the target membrane, before the virion is swept away after shedding the S1 domain connecting it to the host cell receptor. Second, we conclude that the double attachment and the conserved disulfide bridge establish the strong anchoring required for subsequent membrane fusion. Multiple distinct membrane-anchoring elements ensure high avidity and high mechanical strength of FP-membrane binding.
Transition path sampling is a powerful tool in the study of rare events. Shooting trial trajectories from configurations along existing transition paths proved particularly efficient in the sampling of reactive trajectories. However, most shooting attempts tend not to result in transition paths, in particular in cases where the transition dynamics has diffusive character. To overcome the resulting efficiency problem, we developed an algorithm for “shooting from the top.” We first define a shooting range through which all paths have to pass and then shoot off trial trajectories only from within this range. For a well chosen shooting range, nearly every shot is successful, resulting in an accepted transition path. To deal with multiple mechanisms, weighted shooting ranges can be used. To cope with the problem of unsuitably placed shooting ranges, we developed an algorithm that iteratively improves the location of the shooting range. The transition path sampling procedure is illustrated for models of diffusive and Langevin dynamics. The method should be particularly useful in cases where the transition paths are long so that only relatively few shots are possible, yet reasonable order parameters are known.
Ion channel gating is essential for cellular homeostasis and is tightly controlled. In some eukaryotic and most bacterial ligand-gated K+ channels, RCK domains regulate ion fluxes. Until now, a single regulatory mechanism has been proposed for all RCK-regulated channels, involving signal transduction from the RCK domain to the gating area. Here, we present an inactive ADP-bound structure of KtrAB from Vibrio alginolyticus, determined by cryo-electron microscopy, which, combined with EPR spectroscopy and molecular dynamics simulations, uncovers a novel regulatory mechanism for ligand-induced action at a distance. Exchange of activating ATP to inactivating ADP triggers short helical segments in the K+-translocating KtrB dimer to organize into two long helices that penetrate deeply into the regulatory RCK domains, thus connecting nucleotide-binding sites and ion gates. As KtrAB and its homolog TrkAH have been implicated as bacterial pathogenicity factors, the discovery of this functionally relevant inactive conformation may advance structure-guided drug development.
New drugs are urgently needed to combat the global TB epidemic. Targeting simultaneously multiple respiratory enzyme complexes of Mycobacterium tuberculosis is regarded as one of the most effective treatment options to shorten drug administration regimes, and reduce the opportunity for the emergence of drug resistance. During infection and proliferation, the cytochrome bd oxidase plays a crucial role for mycobacterial pathophysiology by maintaining aerobic respiration at limited oxygen concentrations. Here, we present the cryo-EM structure of the cytochrome bd oxidase from M. tuberculosis at 2.5 Å. In conjunction with atomistic molecular dynamics (MD) simulation studies we discovered a previously unknown MK-9-binding site, as well as a unique disulfide bond within the Q-loop domain that defines an inactive conformation of the canonical quinol oxidation site in Actinobacteria. Our detailed insights into the long-sought atomic framework of the cytochrome bd oxidase from M. tuberculosis will form the basis for the design of highly specific drugs to act on this enzyme.
Transport of lipids across membranes is fundamental for diverse biological pathways in cells. Multiple ion-coupled transporters participate in lipid translocation, but their mechanisms remain largely unknown. Major facilitator superfamily (MFS) lipid transporters play central roles in cell wall synthesis, brain development and function, lipids recycling, and cell signaling. Recent structures of MFS lipid transporters revealed overlapping architectural features pointing towards a common mechanism. Here we used cysteine disulfide trapping, molecular dynamics simulations, mutagenesis analysis, and transport assays in vitro and in vivo, to investigate the mechanism of LtaA, a proton-dependent MFS lipid transporter essential for lipoteichoic acids synthesis in the pathogen Staphylococcus aureus. We reveal that LtaA displays asymmetric lateral openings with distinct functional relevance and that cycling through outward- and inward-facing conformations is essential for transport activity. We demonstrate that while the entire amphipathic central cavity of LtaA contributes to lipid binding, its hydrophilic pocket dictates substrate specificity. We propose that LtaA catalyzes lipid translocation by a ‘trap-and-flip’ mechanism that might be shared among MFS lipid transporters.
Transport of lipids across membranes is fundamental for diverse biological pathways in cells. Multiple ion-coupled transporters take part in lipid translocation, but their mechanisms remain largely unknown. Major facilitator superfamily (MFS) lipid transporters play central roles in cell wall synthesis, brain development and function, lipids recycling, and cell signaling. Recent structures of MFS lipid transporters revealed overlapping architectural features pointing towards a common mechanism. Here we used cysteine disulfide trapping, molecular dynamics simulations, mutagenesis analysis, and transport assays in vitro and in vivo, to investigate the mechanism of LtaA, a proton-dependent MFS lipid transporter essential for lipoteichoic acid synthesis in the pathogen Staphylococcus aureus. We reveal that LtaA displays asymmetric lateral openings with distinct functional relevance and that cycling through outward- and inward-facing conformations is essential for transport activity. We demonstrate that while the entire amphipathic central cavity of LtaA contributes to lipid binding, its hydrophilic pocket dictates substrate specificity. We propose that LtaA catalyzes lipid translocation by a ‘trap-and-flip’ mechanism that might be shared among MFS lipid transporters.
Binding of the spike protein of SARS-CoV-2 to the human angiotensin converting enzyme 2 (ACE2) receptor triggers translocation of the virus into cells. Both the ACE2 receptor and the spike protein are heavily glycosylated, including at sites near their binding interface. We built fully glycosylated models of the ACE2 receptor bound to the receptor binding domain (RBD) of the SARS-CoV-2 spike protein. Using atomistic molecular dynamics (MD) simulations, we found that the glycosylation of the human ACE2 receptor contributes substantially to the binding of the virus. Interestingly, the glycans at two glycosylation sites, N90 and N322, have opposite effects on spike protein binding. The glycan at the N90 site partly covers the binding interface of the spike RBD. Therefore, this glycan can interfere with the binding of the spike protein and protect against docking of the virus to the cell. By contrast, the glycan at the N322 site interacts tightly with the RBD of the ACE2-bound spike protein and strengthens the complex. Remarkably, the N322 glycan binds into a conserved region of the spike protein identified previously as a cryptic epitope for a neutralizing antibody. By mapping the glycan binding sites, our MD simulations aid in the targeted development of neutralizing antibodies and SARS-CoV-2 fusion inhibitors.
IHMCIF: an extension of the PDBx/mmCIF data standard for integrative structure determination methods
(2024)
IHMCIF (github.com/ihmwg/IHMCIF) is a data information framework that supports archiving and disseminating macromolecular structures determined by integrative or hybrid modeling (IHM), and making them Findable, Accessible, Interoperable, and Reusable (FAIR). IHMCIF is an extension of the Protein Data Bank Exchange/macromolecular Crystallographic Information Framework (PDBx/mmCIF) that serves as the framework for the Protein Data Bank (PDB) to archive experimentally determined atomic structures of biological macromolecules and their complexes with one another and small molecule ligands (e.g., enzyme cofactors and drugs). IHMCIF serves as the foundational data standard for the PDB-Dev prototype system, developed for archiving and disseminating integrative structures. It utilizes a flexible data representation to describe integrative structures that span multiple spatiotemporal scales and structural states with definitions for restraints from a variety of experimental methods contributing to integrative structural biology. The IHMCIF extension was created with the benefit of considerable community input and recommendations gathered by the Worldwide Protein Data Bank (wwPDB) Task Force for Integrative or Hybrid Methods (wwpdb.org/task/hybrid). Herein, we describe the development of IHMCIF to support evolving methodologies and ongoing advancements in integrative structural biology. Ultimately, IHMCIF will facilitate the unification of PDB-Dev data and tools with the PDB archive so that integrative structures can be archived and disseminated through PDB.
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.
Cells maintain membrane fluidity by regulating lipid saturation, but the molecular mechanisms of this homeoviscous adaptation remain poorly understood. We have reconstituted the core machinery for regulating lipid saturation in baker’s yeast to study its molecular mechanism. By combining molecular dynamics simulations with experiments, we uncover a remarkable sensitivity of the transcriptional regulator Mga2 to the abundance, position, and configuration of double bonds in lipid acyl chains, and provide insights into the molecular rules of membrane adaptation. Our data challenge the prevailing hypothesis that membrane fluidity serves as the measured variable for regulating lipid saturation. Rather, we show that Mga2 senses the molecular lipid-packing density in a defined region of the membrane. Our findings suggest that membrane property sensors have evolved remarkable sensitivities to highly specific aspects of membrane structure and dynamics, thus paving the way toward the development of genetically encoded reporters for such properties in the future.
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/DEER 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 is no prerequisite for dimerization. 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.
In response to pathogen infection, gasdermin (GSDM) proteins form membrane pores that induce a host cell death process called pyroptosis1–3. Studies of human and mouse GSDM pores reveal the functions and architectures of 24–33 protomers assemblies4–9, but the mechanism and evolutionary origin of membrane targeting and GSDM pore formation remain unknown. Here we determine a structure of a bacterial GSDM (bGSDM) pore and define a conserved mechanism of pore assembly. Engineering a panel of bGSDMs for site-specific proteolytic activation, we demonstrate that diverse bGSDMs form distinct pore sizes that range from smaller mammalian-like assemblies to exceptionally large pores containing >50 protomers. We determine a 3.3 Å cryo-EM structure of a Vitiosangium bGSDM in an active slinky-like oligomeric conformation and analyze bGSDM pores in a native lipid environment to create an atomic-level model of a full 52-mer bGSDM pore. Combining our structural analysis with molecular dynamics simulations and cellular assays, our results support a stepwise model of GSDM pore assembly and suggest that a covalently bound palmitoyl can leave a hydrophobic sheath and insert into the membrane before formation of the membrane-spanning β-strand regions. These results reveal the diversity of GSDM pores found in nature and explain the function of an ancient post-translational modification in enabling programmed host cell death.
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.
Highlights
• Sampling the large conformational space of disordered proteins requires extensive molecular dynamics (MD) simulations.
• Fragment assembly complements MD simulations to produce extensive ensembles of disordered proteins with atomic detail.
• Hierarchical chain growth (HCG) ensembles capture key experimental descriptors “out of the box”.
• HCG has revealed local structural characteristics associated with protein dysfunction in neurodegeneration.
Abstract
Disordered proteins and nucleic acids play key roles in cellular function and disease. Here, we review recent advances in the computational exploration of the conformational dynamics of flexible biomolecules. While atomistic molecular dynamics (MD) simulation has seen a lot of improvement in recent years, large-scale computing resources and careful validation are required to simulate full-length disordered biopolymers in solution. As a computationally efficient alternative, hierarchical chain growth (HCG) combines pre-sampled chain fragments in a statistically reproducible manner into ensembles of full-length atomically detailed biomolecular structures. Experimental data can be integrated during and after chain assembly. Applications to the neurodegeneration-linked proteins α-synuclein, tau, and TDP-43, including as condensate, illustrate the use of HCG. We conclude by highlighting the emerging connections to AI-based structural modeling including AlphaFold2.
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.
Binding of the spike protein of SARS-CoV-2 to the human angiotensin-converting enzyme 2 (ACE2) receptor triggers translocation of the virus into cells. Both the ACE2 receptor and the spike protein are heavily glycosylated, including at sites near their binding interface. We built fully glycosylated models of the ACE2 receptor bound to the receptor binding domain (RBD) of the SARS-CoV-2 spike protein. Using atomistic molecular dynamics (MD) simulations, we found that the glycosylation of the human ACE2 receptor contributes substantially to the binding of the virus. Interestingly, the glycans at two glycosylation sites, N90 and N322, have opposite effects on spike protein binding. The glycan at the N90 site partly covers the binding interface of the spike RBD. Therefore, this glycan can interfere with the binding of the spike protein and protect against docking of the virus to the cell. By contrast, the glycan at the N322 site interacts tightly with the RBD of the ACE2-bound spike protein and strengthens the complex. Remarkably, the N322 glycan binds to a conserved region of the spike protein identified previously as a cryptic epitope for a neutralizing antibody. By mapping the glycan binding sites, our MD simulations aid in the targeted development of neutralizing antibodies and SARS-CoV-2 fusion inhibitors.
Abstract
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.
Author summary
The SARS-CoV-2 virus has caused a global health crisis. The spike protein exposed at its surface is key for infection and the primary antibody target. However, spike is covered by highly mobile glycan molecules that could impair antibody binding. To identify accessible epitopes, we performed molecular dynamics simulations of an atomistic model of glycosylated spike embedded in a membrane. By combining extensive simulations with bioinformatics analyses, we recovered known antibody binding sites and identified several epitope candidates as targets for further vaccine development.
Human feline leukemia virus subgroup C receptor-related proteins 1 and 2 (FLVCR1 and 2) are members of the major facilitator superfamily1. Their dysfunction is linked to several clinical disorders, including PCARP, HSAN, and Fowler syndrome2–7. Earlier studies concluded that FLVCR1 may function as a putative heme exporter8–12, while FLVCR2 was suggested to act as a heme importer13, yet conclusive biochemical and detailed molecular evidence remained elusive for the function of both transporters14–17. Here, we show that FLVCR1 and FLVCR2 facilitate the transport of choline and ethanolamine across human plasma membranes, utilizing a concentration-driven substrate translocation process. Through structural and computational analyses, we have identified distinct conformational states of FLVCRs and unraveled the coordination chemistry underlying their substrate interactions. Within the binding pocket of both transporters, we identify fully conserved tryptophan and tyrosine residues holding a central role in the formation of cation-π interactions, essential for choline and ethanolamine selectivity. Our findings not only clarify the mechanisms of choline and ethanolamine transport by FLVCR1 and FLVCR2, enhancing our comprehension of disease-associated mutations that interfere with these vital processes, but also shed light on the conformational dynamics of these MFS-type proteins during the transport cycle.
During the co-translational assembly of protein complexes, a fully synthesized subunit engages with the nascent chain of a newly synthesized interaction partner. Such events are thought to contribute to productive assembly, but their exact physiological relevance remains underexplored. Here, we examine structural motifs contained in nucleoporins for their potential to facilitate co-translational assembly. We experimentally test candidate structural motifs and identify several previously unknown co-translational interactions. We demonstrate by selective ribosome profiling that domain invasion motifs of beta-propellers, coiled-coils, and short linear motifs may act as co-translational assembly domains. Such motifs are often contained in proteins that are members of multiple complexes (moonlighters) and engage with closely related paralogs. Surprisingly, moonlighters and paralogs assemble co-translationally in only some but not all of the relevant biogenesis pathways. Our results highlight the regulatory complexity of assembly pathways.
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.
In response to pathogen infection, gasdermin (GSDM) proteins form membrane pores that induce a host cell death process called pyroptosis1–3. Studies of human and mouse GSDM pores reveal the functions and architectures of 24–33 protomers assemblies4–9, but the mechanism and evolutionary origin of membrane targeting and GSDM pore formation remain unknown. Here we determine a structure of a bacterial GSDM (bGSDM) pore and define a conserved mechanism of pore assembly. Engineering a panel of bGSDMs for site-specific proteolytic activation, we demonstrate that diverse bGSDMs form distinct pore sizes that range from smaller mammalian-like assemblies to exceptionally large pores containing >50 protomers. We determine a 3.3 Å cryo-EM structure of a Vitiosangium bGSDM in an active slinky-like oligomeric conformation and analyze bGSDM pores in a native lipid environment to create an atomic-level model of a full 52-mer bGSDM pore. Combining our structural analysis with molecular dynamics simulations and cellular assays, we define a stepwise model of GSDM pore assembly and demonstrate that pore formation is driven by local unfolding of membrane-spanning β-strand regions and pre-insertion of a covalently bound palmitoyl into the target membrane. These results yield insights into the diversity of GSDM pores found in nature and the function of an ancient post-translational modification in enabling a programmed host cell death process.
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.
Snapshots of acetyl-CoA synthesis, the final step of CO₂ fixation in the Wood-Ljungdahl pathway
(2024)
In the ancient microbial Wood-Ljungdahl pathway, CO2 is fixed in a multi-step process with acetyl-CoA synthesis at the bifunctional carbon monoxide dehydrogenase/acetyl-CoA synthase complex (CODH/ACS). Here, we present catalytic snapshots of the CODH/ACS from the gas-converting acetogen Clostridium autoethanogenum, characterizing the molecular choreography of the overall reaction including electron transfer to the CODH for CO2 reduction, methyl transfer from the corrinoid iron-sulfur protein (CoFeSP) partner to the ACS active site and acetyl-CoA production. Unlike CODH, the multidomain ACS undergoes large conformational changes to form an internal connection to the CODH active site, accommodate the CoFeSP for methyl transfer and protect the reaction intermediates. Altogether, the structures allow us to draw a detailed reaction mechanism of this enzyme crucial for CO2 fixation in anaerobic organisms.
Human feline leukaemia virus subgroup C receptor-related proteins 1 and 2 (FLVCR1 and FLVCR2) are members of the major facilitator superfamily1. Their dysfunction is linked to several clinical disorders, including PCARP, HSAN and Fowler syndrome2,3,4,5,6,7. Earlier studies concluded that FLVCR1 may function as a haem exporter8,9,10,11,12, whereas FLVCR2 was suggested to act as a haem importer13, yet conclusive biochemical and detailed molecular evidence remained elusive for the function of both transporters14,15,16. Here, we show that FLVCR1 and FLVCR2 facilitate the transport of choline and ethanolamine across the plasma membrane, using a concentration-driven substrate translocation process. Through structural and computational analyses, we have identified distinct conformational states of FLVCRs and unravelled the coordination chemistry underlying their substrate interactions. Fully conserved tryptophan and tyrosine residues form the binding pocket of both transporters and confer selectivity for choline and ethanolamine through cation–π interactions. Our findings clarify the mechanisms of choline and ethanolamine transport by FLVCR1 and FLVCR2, enhance our comprehension of disease-associated mutations that interfere with these vital processes and shed light on the conformational dynamics of these major facilitator superfamily proteins during the transport cycle.