Refine
Year of publication
- 2020 (77) (remove)
Document Type
- Preprint (77) (remove)
Language
- English (77)
Has Fulltext
- yes (77)
Is part of the Bibliography
- no (77)
Institute
- Physik (77) (remove)
Disordered proteins and nucleic acids can condense into droplets that resemble the membraneless organelles observed in living cells. MD simulations offer a unique tool to characterize the molecular interactions governing the formation of these biomolecular condensates, their physico-chemical properties, and the factors controlling their composition and size. However, biopolymer condensation depends sensitively on the balance between different energetic and entropic contributions. Here, we develop a general strategy to fine-tune the potential energy function for molecular dynamics simulations of biopolymer phase separation. We rebalance protein-protein interactions against solvation and entropic contributions to match the excess free energy of transferring proteins between dilute solution and condensate. We illustrate this formalism by simulating liquid droplet formation of the FUS low complexity domain (LCD) with a rebalanced MARTINI model. By scaling the strength of the nonbonded interactions in the coarse-grained MARTINI potential energy function, we map out a phase diagram in the plane of protein concentration and interaction strength. Above a critical scaling factor of αc ≈ 0.6, FUS LCD condensation is observed, where α = 1 and 0 correspond to full and repulsive interactions in the MARTINI model, respectively. For a scaling factor α = 0.65, we recover the experimental densities of the dilute and dense phases, and thus the excess protein transfer free energy into the droplet and the saturation concentration where FUS LCD condenses. In the region of phase separation, we simulate FUS LCD droplets of four different sizes in stable equilibrium with the dilute phase and slabs of condensed FUS LCD for tens of microseconds, and over one millisecond in aggregate. We determine surface tensions in the range of 0.01 to 0.4mN/m from the fluctuations of the droplet shape and from the capillary-wave-like broadening of the interface between the two phases. From the dynamics of the protein end-to-end distance, we estimate shear viscosities from 0.001 to 0.02Pas for the FUS LCD droplets with scaling factors α in the range of 0.625 to 0.75, where we observe liquid droplets. Significant hydration of the interior of the droplets keeps the proteins mobile and the droplets fluid.
The spike (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.
i-TED is an innovative detection system which exploits Compton imaging techniques to achieve a superior signal-to-background ratio in (n,γ) cross-section measurements using time-of-flight technique. This work presents the first experimental validation of the i-TED apparatus for high-resolution time-of-flight experiments and demonstrates for the first time the concept proposed for background rejection. To this aim both 197Au(n,γ) and 56Fe(n,γ) reactions were measured at CERN n\_TOF using an i-TED demonstrator based on only three position-sensitive detectors. Two \cds detectors were also used to benchmark the performance of i-TED. The i-TED prototype built for this study shows a factor of ∼3 higher detection sensitivity than state-of-the-art \cds detectors in the ∼10~keV neutron energy range of astrophysical interest. This paper explores also the perspectives of further enhancement in performance attainable with the final i-TED array consisting of twenty position-sensitive detectors and new analysis methodologies based on Machine-Learning techniques.
Inclusive ψ(2S) production is measured in p-Pb collisions at the centre-of-mass energy per nucleon-nucleon pair sNN−−−√=8.16 TeV, using the ALICE detector at the CERN LHC. The production of ψ(2S) is studied at forward (2.03<ycms<3.53) and backward (−4.46<ycms<−2.96) centre-of-mass rapidity and for transverse momentum pT < 12 GeV/c via the decay to muon pairs. In this paper, we report the integrated as well as the ycms- and pT-differential inclusive production cross sections. Nuclear effects on ψ(2S) production are studied via the determination of the nuclear modification factor that shows a strong suppression at both forward and backward centre-of-mass rapidities. Comparisons with corresponding results for inclusive J/ψ show a similar suppression for the two states at forward rapidity (p-going direction), but a stronger suppression for ψ(2S) at backward rapidity (Pb-going direction). As a function of pT, no clear dependence of the nuclear modification factor is found. The relative size of nuclear effects on ψ(2S) production compared to J/ψ is also studied via the double ratio of production cross sections [σψ(2S)/σJ/ψ]pPb/[σψ(2S)/σJ/ψ]pp between p-Pb and pp collisions. The results are compared with theoretical models that include various effects related to the initial and final state of the collision system and also with previous measurements at sNN−−−√ = 5.02 TeV.
The production of the Λ(1520) baryonic resonance has been measured at midrapidity in inelastic pp collisions at s√ = 7 TeV and in p-Pb collisions at sNN−−−√ = 5.02 TeV for non-single diffractive events and in multiplicity classes. The resonance is reconstructed through its hadronic decay channel Λ(1520) → pK− and the charge conjugate with the ALICE detector. The integrated yields and mean transverse momenta are calculated from the measured transverse momentum distributions in pp and p-Pb collisions. The mean transverse momenta follow mass ordering as previously observed for other hyperons in the same collision systems. A Blast-Wave function constrained by other light hadrons (π, K, K0S, p, Λ) describes the shape of the Λ(1520) transverse momentum distribution up to 3.5 GeV/c in p-Pb collisions. In the framework of this model, this observation suggests that the Λ(1520) resonance participates in the same collective radial flow as other light hadrons. The ratio of the yield of Λ(1520) to the yield of the ground state particle Λ remains constant as a function of charged-particle multiplicity, suggesting that there is no net effect of the hadronic phase in p-Pb collisions on the Λ(1520) yield.
Recurrent cortical network dynamics plays a crucial role for sequential information processing in the brain. While the theoretical framework of reservoir computing provides a conceptual basis for the understanding of recurrent neural computation, it often requires manual adjustments of global network parameters, in particular of the spectral radius of the recurrent synaptic weight matrix. Being a mathematical and relatively complex quantity, the spectral radius is not readily accessible to biological neural networks, which generally adhere to the principle that information about the network state should either be encoded in local intrinsic dynamical quantities (e.g. membrane potentials), or transmitted via synaptic connectivity. We present two synaptic scaling rules for echo state networks that solely rely on locally accessible variables. Both rules work online, in the presence of a continuous stream of input signals. The first rule, termed flow control, is based on a local comparison between the mean squared recurrent membrane potential and the mean squared activity of the neuron itself. It is derived from a global scaling condition on the dynamic flow of neural activities and requires the separability of external and recurrent input currents. We gained further insight into the adaptation dynamics of flow control by using a mean field approximation on the variances of neural activities that allowed us to describe the interplay between network activity and adaptation as a two-dimensional dynamical system. The second rule that we considered, variance control, directly regulates the variance of neural activities by locally scaling the recurrent synaptic weights. The target set point of this homeostatic mechanism is dynamically determined as a function of the variance of the locally measured external input. This functional relation was derived from the same mean-field approach that was used to describe the approximate dynamics of flow control.
The effectiveness of the presented mechanisms was tested numerically using different external input protocols. The network performance after adaptation was evaluated by training the network to perform a time delayed XOR operation on binary sequences. As our main result, we found that flow control can reliably regulate the spectral radius under different input statistics, but precise tuning is negatively affected by interneural correlations. Furthermore, flow control showed a consistent task performance over a wide range of input strengths/variances. Variance control, on the other side, did not yield the desired spectral radii with the same precision. Moreover, task performance was less consistent across different input strengths.
Given the better performance and simpler mathematical form of flow control, we concluded that a local control of the spectral radius via an implicit adaptation scheme is a realistic alternative to approaches using classical “set point” homeostatic feedback controls of neural firing.
Author summary How can a neural network control its recurrent synaptic strengths such that network dynamics are optimal for sequential information processing? An important quantity in this respect, the spectral radius of the recurrent synaptic weight matrix, is a non-local quantity. Therefore, a direct calculation of the spectral radius is not feasible for biological networks. However, we show that there exist a local and biologically plausible adaptation mechanism, flow control, which allows to control the recurrent weight spectral radius while the network is operating under the influence of external inputs. Flow control is based on a theorem of random matrix theory, which is applicable if inter-synaptic correlations are weak. We apply the new adaption rule to echo-state networks having the task to perform a time-delayed XOR operation on random binary input sequences. We find that flow-controlled networks can adapt to a wide range of input strengths while retaining essentially constant task performance.
Recurrent cortical network dynamics plays a crucial role for sequential information processing in the brain. While the theoretical framework of reservoir computing provides a conceptual basis for the understanding of recurrent neural computation, it often requires manual adjustments of global network parameters, in particular of the spectral radius of the recurrent synaptic weight matrix. Being a mathematical and relatively complex quantity, the spectral radius is not readily accessible to biological neural networks, which generally adhere to the principle that information about the network state should either be encoded in local intrinsic dynamical quantities (e.g. membrane potentials), or transmitted via synaptic connectivity. We present two synaptic scaling rules for echo state networks that solely rely on locally accessible variables. Both rules work online, in the presence of a continuous stream of input signals. The first rule, termed flow control, is based on a local comparison between the mean squared recurrent membrane potential and the mean squared activity of the neuron itself. It is derived from a global scaling condition on the dynamic flow of neural activities and requires the separability of external and recurrent input currents. We gained further insight into the adaptation dynamics of flow control by using a mean field approximation on the variances of neural activities that allowed us to describe the interplay between network activity and adaptation as a two-dimensional dynamical system. The second rule that we considered, variance control, directly regulates the variance of neural activities by locally scaling the recurrent synaptic weights. The target set point of this homeostatic mechanism is dynamically determined as a function of the variance of the locally measured external input. This functional relation was derived from the same mean-field approach that was used to describe the approximate dynamics of flow control.
The effectiveness of the presented mechanisms was tested numerically using different external input protocols. The network performance after adaptation was evaluated by training the network to perform a time delayed XOR operation on binary sequences. As our main result, we found that flow control can reliably regulate the spectral radius under different input statistics, but precise tuning is negatively affected by interneural correlations. Furthermore, flow control showed a consistent task performance over a wide range of input strengths/variances. Variance control, on the other side, did not yield the desired spectral radii with the same precision. Moreover, task performance was less consistent across different input strengths.
Given the better performance and simpler mathematical form of flow control, we concluded that a local control of the spectral radius via an implicit adaptation scheme is a realistic alternative to approaches using classical “set point” homeostatic feedback controls of neural firing.
Author summary How can a neural network control its recurrent synaptic strengths such that network dynamics are optimal for sequential information processing? An important quantity in this respect, the spectral radius of the recurrent synaptic weight matrix, is a non-local quantity. Therefore, a direct calculation of the spectral radius is not feasible for biological networks. However, we show that there exist a local and biologically plausible adaptation mechanism, flow control, which allows to control the recurrent weight spectral radius while the network is operating under the influence of external inputs. Flow control is based on a theorem of random matrix theory, which is applicable if inter-synaptic correlations are weak. We apply the new adaption rule to echo-state networks having the task to perform a time-delayed XOR operation on random binary input sequences. We find that flow-controlled networks can adapt to a wide range of input strengths while retaining essentially constant task performance.
Scattering studies with low-energy kaon-proton femtoscopy in proton–proton collisions at the LHC
(2020)
The study of the strength and behaviour of the antikaon-nucleon (K¯¯¯¯N) interaction constitutes one of the key focuses of the strangeness sector in low-energy Quantum Chromodynamics (QCD). In this letter a unique high-precision measurement of the strong interaction between kaons and protons, close and above the kinematic threshold, is presented. The femtoscopic measurements of the correlation function at low pair-frame relative momentum of (K+ p ⊕ K− p¯¯¯) and (K− p ⊕ K+ p¯¯¯) pairs measured in pp collisions at s√ = 5, 7 and 13 TeV are reported. A structure observed around a relative momentum of 58 MeV/c in the measured correlation function of (K− p ⊕ K+ p¯¯¯) with a significance of 4.4. σ constitutes the first experimental evidence for the opening of the (K¯¯¯¯0n⊕K0n¯¯¯) isospin breaking channel due to the mass difference between charged and neutral kaons. The measured correlation functions have been compared to Jülich and Kyoto models in addition to the Coulomb potential. The high-precision data at low relative momenta presented in this work prove femtoscopy to be a powerful complementary tool to scattering experiments and provide new constraints above the K¯¯¯¯N threshold for low-energy QCD chiral models.