Institutes
Refine
Document Type
- Preprint (44) (remove)
Language
- English (44)
Has Fulltext
- yes (44)
Is part of the Bibliography
- no (44)
Keywords
- Gross-Neveu model (3)
- inhomogeneous phases (3)
- mean-field (3)
- moat regime (3)
- phase diagram (3)
- stability analysis (3)
- two-point function (3)
- wave-function renormalization (3)
- FOS: Physical sciences (1)
- Nuclear Experiment (nucl-ex) (1)
Institute
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.
DNA binding redistributes activation domain ensemble and accessibility in pioneer factor Sox2
(2023)
More than 1600 human transcription factors orchestrate the transcriptional machinery to control gene expression and cell fate. Their function is conveyed through intrinsically disordered regions (IDRs) containing activation or repression domains but lacking quantitative structural ensemble models prevents their mechanistic decoding. Here we integrate single-molecule FRET and NMR spectroscopy with molecular simulations showing that DNA binding can lead to complex changes in the IDR ensemble and accessibility. The C-terminal IDR of pioneer factor Sox2 is highly disordered but its conformational dynamics are guided by weak and dynamic charge interactions with the folded DNA binding domain. Both DNA and nucleosome binding induce major rearrangements in the IDR ensemble without affecting DNA binding affinity. Remarkably, interdomain interactions are redistributed in complex with DNA leading to variable exposure of two activation domains critical for transcription. Charged intramolecular interactions allowing for dynamic redistributions may be common in transcription factors and necessary for sensitive tuning of structural ensembles.
We study the μ-μ45-T phase diagram of the 2+1-dimensional Gross-Neveu model, where μ denotes the ordinary chemical potential, μ45 the chiral chemical potential and T the temperature. We use the mean-field approximation and two different lattice regularizations with naive chiral fermions. An inhomogeneous phase at finite lattice spacing is found for one of the two regularizations. Our results suggest that there is no inhomogeneous phase in the continuum limit. We show that a chiral chemical potential is equivalent to an isospin chemical potential. Thus, all results presented in this work can also be interpreted in the context of isospin imbalance.
The phase diagram of the (1+1)-dimensional Gross-Neveu model is reanalyzed for (non-)zero chemical potential and (non-)zero temperature within the mean-field approximation. By investigating the momentum dependence of the bosonic two-point function, the well-known second-order phase transition from the Z2 symmetric phase to the so-called inhomogeneous phase is detected. In the latter phase the chiral condensate is periodically varying in space and translational invariance is broken. This work is a proof of concept study that confirms that it is possible to correctly localize second-order phase transition lines between phases without condensation and phases of spatially inhomogeneous condensation via a stability analysis of the homogeneous phase. To complement other works relying on this technique, the stability analysis is explained in detail and its limitations and successes are discussed in context of the Gross-Neveu model. Additionally, we present explicit results for the bosonic wave-function renormalization in the mean-field approximation, which is extracted analytically from the bosonic two-point function. We find regions -- a so-called moat regime -- where the wave function renormalization is negative accompanying the inhomogeneous phase as expected.
The phase diagram of the (1+1)-dimensional Gross-Neveu model is reanalyzed for (non-)zero chemical potential and (non-)zero temperature within the mean-field approximation. By investigating the momentum dependence of the bosonic two-point function, the well-known second-order phase transition from the Z2 symmetric phase to the so-called inhomogeneous phase is detected. In the latter phase the chiral condensate is periodically varying in space and translational invariance is broken. This work is a proof of concept study that confirms that it is possible to correctly localize second-order phase transition lines between phases without condensation and phases of spatially inhomogeneous condensation via a stability analysis of the homogeneous phase. To complement other works relying on this technique, the stability analysis is explained in detail and its limitations and successes are discussed in context of the Gross-Neveu model. Additionally, we present explicit results for the bosonic wave-function renormalization in the mean-field approximation, which is extracted analytically from the bosonic two-point function. We find regions -- a so-called moat regime -- where the wave function renormalization is negative accompanying the inhomogeneous phase as expected.
The phase diagram of the (1+1)-dimensional Gross-Neveu model is reanalyzed for (non-)zero chemical potential and (non-)zero temperature within the mean-field approximation. By investigating the momentum dependence of the bosonic two-point function, the well-known second-order phase transition from the Z2 symmetric phase to the so-called inhomogeneous phase is detected. In the latter phase the chiral condensate is periodically varying in space and translational invariance is broken. This work is a proof of concept study that confirms that it is possible to correctly localize second-order phase transition lines between phases without condensation and phases of spatially inhomogeneous condensation via a stability analysis of the homogeneous phase. To complement other works relying on this technique, the stability analysis is explained in detail and its limitations and successes are discussed in context of the Gross-Neveu model. Additionally, we present explicit results for the bosonic wave-function renormalization in the mean-field approximation, which is extracted analytically from the bosonic two-point function. We find regions -- a so-called moat regime -- where the wave function renormalization is negative accompanying the inhomogeneous phase as expected.
We prove that the projectivized strata of differentials are not contained in pointed Brill-Noether divisors, with only a few exceptions. For a generic element in a stratum of differentials, we show that many of the associated pointed Brill-Noether loci are of expected dimension. We use our results to study the Auel-Haburcak Conjecture: We obtain new non-containments between maximal Brill-Noether loci in Mg. Our results regarding quadratic differentials imply that the quadratic strata in genus 6 are uniruled.
We study the μ-μ45-T phase diagram of the 2+1-dimensional Gross-Neveu model, where μ denotes the ordinary chemical potential, μ45 the chiral chemical potential and T the temperature. We use the mean-field approximation and two different lattice regularizations with naive chiral fermions. An inhomogeneous phase at finite lattice spacing is found for one of the two regularizations. Our results suggest that there is no inhomogeneous phase in the continuum limit. We show that a chiral chemical potential is equivalent to an isospin chemical potential. Thus, all results presented in this work can also be interpreted in the context of isospin imbalance.
We explore the phase structure of the 1+1 dimensional Gross-Neveu model at finite number of fermion flavors using lattice field theory. Besides a chirally symmetric phase and a homogeneously broken phase we find evidence for the existence of an inhomogeneous phase, where the condensate is a spatially oscillating function. Our numerical results include a crude μ-T phase diagram.
In this work, the phase diagram of the 2+1-dimensional Gross-Neveu model is investigated with baryon chemical potential as well as chiral chemical potential in the mean-field approximation. We study the theory using two lattice discretizations, which are both based on naive fermions. An inhomogeneous chiral phase is observed only for one of the two discretizations. Our results suggest that this phase disappears in the continuum limit.
We continue previous investigations of the (inhomogeneous) phase structure of the Gross-Neveu model in a noninteger number of spatial dimensions (1≤d<3) in the limit of an infinite number of fermion species (N→∞) at (non)zero chemical potential μ. In this work, we extend the analysis from zero to nonzero temperature T.
The phase diagram of the Gross-Neveu model in 1≤d<3 spatial dimensions is well known under the assumption of spatially homogeneous condensation with both a symmetry broken and a symmetric phase present for all spatial dimensions. In d=1 one additionally finds an inhomogeneous phase, where the order parameter, the condensate, is varying in space. Similarly, phases of spatially varying condensates are also found in the Gross-Neveu model in d=2 and d=3, as long as the theory is not fully renormalized, i.e., in the presence of a regulator. For d=2, one observes that the inhomogeneous phase vanishes, when the regulator is properly removed (which is not possible for d=3 without introducing additional parameters).
In the present work, we use the stability analysis of the symmetric phase to study the presence (for 1≤d<2) and absence (for 2≤d<3) of these inhomogeneous phases and the related moat regimes in the fully renormalized Gross-Neveu model in the μ,T-plane. We also discuss the relation between "the number of spatial dimensions" and "studying the model with a finite regulator" as well as the possible consequences for the limit d→3.
Inhomogeneous condensation in the Gross-Neveu model in noninteger spatial dimensions 1 ≤ d < 3
(2023)
The Gross-Neveu model in the N→∞ approximation in d=1 spatial dimensions exhibits a chiral inhomogeneous phase (IP), where the chiral condensate has a spatial dependence that spontaneously breaks translational invariance and the Z2 chiral symmetry. This phase is absent in d=2, while in d=3 its existence and extent strongly depends on the regularization and the value of the finite regulator. This work connects these three results smoothly by extending the analysis to non-integer spatial dimensions 1≤d<3, where the model is fully renormalizable. To this end, we adapt the stability analysis, which probes the stability of the homogeneous ground state under inhomogeneous perturbations, to non-integer spatial dimensions. We find that the IP is present for all d<2 and vanishes exactly at d=2. Moreover, we find no instability towards an IP for 2≤d<3, which suggests that the IP in d=3 is solely generated by the presence of a regulator.
We show the absence of an instability of homogeneous (chiral) condensates against spatially inhomogeneous perturbations for various 2+1-dimensional four-fermion and Yukawa models. All models are studied at non-zero baryon chemical potential, while some of them are also subjected to chiral and isospin chemical potential. The considered theories contain up to 16 Lorentz-(pseudo)scalar fermionic interaction channels. We prove the stability of homogeneous condensates by analyzing the bosonic two-point function, which can be expressed in a purely analytical form at zero temperature. Our analysis is presented in a general manner for all of the different discussed models. We argue that the absence of an inhomogeneous chiral phase (where the chiral condensate is spatially non-uniform) follows from this lack of instability. Furthermore, the existence of a moat regime, where the bosonic wave function renormalization is negative, in these models is ruled out.
We show the absence of an instability of homogeneous (chiral) condensates against spatially inhomogeneous perturbations for various 2+1-dimensional four-fermion and Yukawa models. All models are studied at non-zero baryon chemical potential, while some of them are also subjected to chiral and isospin chemical potential. The considered theories contain up to 16 Lorentz-(pseudo)scalar fermionic interaction channels. We prove the stability of homogeneous condensates by analyzing the bosonic two-point function, which can be expressed in a purely analytical form at zero temperature. Our analysis is presented in a general manner for all of the different discussed models. We argue that the absence of an inhomogeneous chiral phase (where the chiral condensate is spatially non-uniform) follows from this lack of instability. Furthermore, the existence of a moat regime, where the bosonic wave function renormalization is negative, in these models is ruled out.
Inhomogeneous condensation in the Gross-Neveu model in non-integer spatial dimensions 1 ≤ d < 3
(2023)
he Gross-Neveu model in the N→∞ approximation in d=1 spatial dimensions exhibits a chiral inhomogeneous phase (IP), where the chiral condensate has a spatial dependence that spontaneously breaks translational invariance and the Z2 chiral symmetry. This phase is absent in d=2, while in d=3 its existence and extent strongly depends on the regularization and the value of the finite regulator. This work connects these three results smoothly by extending the analysis to non-integer spatial dimensions 1≤d<3, where the model is fully renormalizable. To this end, we adapt the stability analysis, which probes the stability of the homogeneous ground state under inhomogeneous perturbations, to non-integer spatial dimensions. We find that the IP is present for all d<2 and vanishes exactly at d=2. Moreover, we find no instability towards an IP for 2≤d<3, which suggests that the IP in d=3 is solely generated by the presence of a regulator.
Vertebrate life depends on renal function to filter excess fluid and remove low-molecular-weight waste products. An essential component of the kidney filtration barrier is the slit diaphragm (SD), a specialized cell-cell junction between podocytes. Although the constituents of the SD are largely known, its molecular organization remains elusive. Here, we use super-resolution correlative light and electron microscopy to quantify a linear rate of reduction in albumin concentration across the filtration barrier under no-flow conditions. Next, we use cryo-electron tomography of vitreous lamellae from high-pressure frozen native glomeruli to analyze the molecular architecture of the SD. The resulting densities resemble a fishnet pattern. Fitting of Nephrin and Neph1, the main constituents of the SD, results in a complex interaction pattern with multiple contact sites between the molecules. Using molecular dynamics simulations, we construct a blueprint of the SD that explains its molecular architecture. Our architectural understanding of the SD reconciles previous findings and provides a mechanistic framework for the development of novel therapies to treat kidney dysfunction.
Vertebrate life depends on renal function to filter excess fluid and remove low-molecular-weight waste products. An essential component of the kidney filtration barrier is the slit diaphragm (SD), a specialized cell-cell junction between podocytes. Although the constituents of the SD are largely known, its molecular organization remains elusive. Here, we use super-resolution correlative light and electron microscopy to quantify a linear rate of reduction in albumin concentration across the filtration barrier. Next, we use cryo-electron tomography of vitreous lamellae from high-pressure frozen native glomeruli to analyze the molecular architecture of the SD. The resulting densities resemble a fishnet pattern. Fitting of Nephrin and Neph1, the main constituents of the SD, results in a complex interaction pattern with multiple contact sites between the molecules. Using molecular dynamics flexible fitting, we construct a blueprint of the SD, where we describe all interactions. Our architectural understanding of the SD reconciles previous findings and provides a mechanistic framework for the development of novel therapies to treat kidney dysfunction.
The most precise measurements to date of the 3ΛH lifetime τ and Λ separation energy BΛ are obtained using the data sample of Pb-Pb collisions at √= 5.02 TeV collected by ALICE at the LHC. The 3ΛH is reconsNN structed via its charged two-body mesonic decay channel (3ΛH→ 3He + π− and the charge-conjugate process). The measured values τ=[253±11 (stat.)±6 (syst.)] ps and BΛ=[102±63 (stat.)±67 (syst.)] keV are compatible with predictions from effective field theories and confirm that the 3ΛH structure is consistent with a weakly-bound system.
The most precise measurements to date of the 3ΛH lifetime τ and Λ separation energy BΛ are obtained using the data sample of Pb-Pb collisions at √sNN = 5.02 TeV collected by ALICE at the LHC. The 3ΛH is reconstructed via its charged two-body mesonic decay channel (3ΛH→ 3He + π− and the charge-conjugate process). The measured values τ=[253±11 (stat.)±6 (syst.)] ps and BΛ=[72±63 (stat.)±36 (syst.)] keV are compatible with predictions from effective field theories and conclusively confirm that the 3ΛH is a weakly-bound system.
The most precise measurements to date of the 3ΛH lifetime τ and Λ separation energy BΛ are obtained using the data sample of Pb-Pb collisions at √sNN = 5.02 TeV collected by ALICE at the LHC. The 3ΛH is reconstructed via its charged two-body mesonic decay channel (3ΛH→ 3He + π− and the charge-conjugate process). The measured values τ=[253±11 (stat.)±6 (syst.)] ps and BΛ=[72±63 (stat.)±35 (syst.)] keV are compatible with predictions from effective field theories and conclusively confirm that the 3ΛH is a weakly-bound system.
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.
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 (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.
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.
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.
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.
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.
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.
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.
Transient receptor potential (TRP) ion channels are among the most well-studied classes of temperature-sensing molecules. Yet, the molecular mechanism and thermodynamic basis for the temperature sensitivity of TRP channels remains to this day poorly understood. One hypothesis is that the temperature-sensing mechanism can simply be described by a difference in heat capacity between the closed and open channel states. While such a two-state model may be simplistic it nonetheless has descriptive value, in the sense that it can be used to to compare overall temperature sensitivity between different channels and mutants. Here, we introduce a mathematical framework based on the two-state model to reliably extract temperature-dependent thermodynamic potentials and heat capacities from measurements of equilibrium constants at different temperatures. Our framework is implemented in an open-source data analysis package that provides a straightforward way to fit both linear and nonlinear van ‘t Hoff plots, thus avoiding some of the previous, potentially erroneous, assumptions when extracting thermodynamic variables from TRP channel electrophysiology data.
Tracking influenza a virus infection in the lung from hematological data with machine learning
(2022)
The tracking of pathogen burden and host responses with minimal-invasive methods during respiratory infections is central for monitoring disease development and guiding treatment decisions. Utilizing a standardized murine model of respiratory Influenza A virus (IAV) infection, we developed and tested different supervised machine learning models to predict viral burden and immune response markers, i.e. cytokines and leukocytes in the lung, from hematological data. We performed independently in vivo infection experiments to acquire extensive data for training and testing purposes of the models. We show here that lung viral load, neutrophil counts, cytokines like IFN-γ and IL-6, and other lung infection markers can be predicted from hematological data. Furthermore, feature analysis of the models shows that blood granulocytes and platelets play a crucial role in prediction and are highly involved in the immune response against IAV. The proposed in silico tools pave the path towards improved tracking and monitoring of influenza infections and possibly other respiratory infections based on minimal-invasively obtained hematological parameters.
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.
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.
Changes in the efficacies of synapses are thought to be the neurobiological basis of learning and memory. The efficacy of a synapse depends on its current number of neurotransmitter receptors. Recent experiments have shown that these receptors are highly dynamic, moving back and forth between synapses on time scales of seconds and minutes. This suggests spontaneous fluctuations in synaptic efficacies and a competition of nearby synapses for available receptors. Here we propose a mathematical model of this competition of synapses for neurotransmitter receptors from a local dendritic pool. Using minimal assumptions, the model produces a fast multiplicative scaling behavior of synapses. Furthermore, the model explains a transient form of heterosynaptic plasticity and predicts that its amount is inversely related to the size of the local receptor pool. Overall, our model reveals logistical tradeoffs during the induction of synaptic plasticity due to the rapid exchange of neurotransmitter receptors between synapses.
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.
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.
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.
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.
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.
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 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.
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.
The electrical and computational properties of neurons in our brains are determined by a rich repertoire of membrane-spanning ion channels and elaborate dendritic trees. However, the precise reason for this inherent complexity remains unknown. Here, we generated large stochastic populations of biophysically realistic hippocampal granule cell models comparing those with all 15 ion channels to their reduced but functional counterparts containing only 5 ion channels. Strikingly, valid parameter combinations in the full models were more frequent and more stable in the face of perturbations to channel expression levels. Scaling up the numbers of ion channels artificially in the reduced models recovered these advantages confirming the key contribution of the actual number of ion channel types. We conclude that the diversity of ion channels gives a neuron greater flexibility and robustness to achieve target excitability.