Non-Markovian data-driven modeling of single-cell motility

Trajectories of human breast cancer cells moving on one-dimensional circular tracks are modeled by thenon-Markovian version of the Langevin equation that includes an arbitrary memory function. When averagedover cells, the velocity distribution exhibits spurious non-Gaussian behavior, while single cells are characterizedby Gaussian velocity distributions. Accordingly, the data are described by a linear memory model whichincludes different random walk models that were previously used to account for various aspects of cell motilitysuch as migratory persistence, non-Markovian effects, colored noise, and anomalous diffusion. The memoryfunction is extracted from the trajectory data without restrictions or assumptions, thus making our approachtruly data driven, and is used for unbiased single-cell comparison. The cell memory displays time-delayedsingle-exponential negative friction, which clearly distinguishes cell motion from the simple persistent randomwalk model and suggests a regulatory feedback mechanism that controls cell migration. Based on the extractedmemory function we formulate a generalized exactly solvable cell migration model which indicates thatnegative friction generates cell persistence over long timescales. The nonequilibrium character of cell motionis investigated by mapping the non-Markovian Langevin equation with memory onto a Markovian model thatinvolves a hidden degree of freedom and is equivalent to the underdamped active Ornstein-Uhlenbeck process

Non-Arrhenius barrier crossing dynamics of non-equilibrium non-Markovian systems

The non-equilibrium non-Markovian barrier crossing dynamics of a one-dimensional massive coordinate, described by the non-equilibrium version of the generalized Langevin equation with unequal random and friction relaxation times, is studied by simulations and analytical methods. Within a harmonic approximation, a general formula for the barrier crossing time is derived which agrees favorably with simulations. Non-equilibrium random forces with a relaxation time longer than the friction relaxation time induce non-Arrhenius behavior and dramatically increase the barrier crossing time; within the harmonic theory this corresponds to a reduced effective temperature which also modifies the spatial and velocity distributions

External Potential Modifies Friction of Molecular Solutes in Water

Stokes’s law for the friction of a sphere in water has been argued to work down to molecular scales, provided the effective hydrodynamic radius includes the hydration layer. In interpretations of experiments and in theoretical models, it is tacitly assumed that the solvent friction experienced by a solute does not depend on whether an external confinement potential acts on the solute. Using a novel method to extract the friction memory function from molecular dynamics simulations, we show that the solvent friction of a strongly harmonically confined methane molecule in water increases by 60% compared to its free-solution value, which is caused by an amplification of the slowest component of the memory function. The friction enhancement occurs for potential strengths typical of physical and chemical bonds and is accompanied by a significant slowing-down of the hydration water dynamics. Thus, the solvent friction acting on molecular solutes is not determined by solvent properties and solute-solvent interactions alone but results from the coupling between solute and solvent dynamics and thereby can be tuned by an external potential acting on the solute. This also explains why simulations of positionally constrained solutes do not reproduce free-solution diffusivities. Dynamic scaling arguments suggest similar effects also for macromolecular solutes provided the solution viscosity is sufficiently enhanced

Pair-Reaction Dynamics in Water: Competition of Memory, Potential Shape, and Inertial Effects

When described by a one-dimensional reaction coordinate, pair-reaction rates in a solvent depend, in addition to the potential barrier height and the friction coefficient, on the potential shape, the effective mass, and the friction relaxation spectrum, but a rate theory that accurately accounts for all of these effects does not exist. After a review of classical reaction-rate theories, we show how to extract all parameters of the generalized Langevin equation (GLE) and, in particular, the friction memory function from molecular dynamics (MD) simulations of two prototypical pair reactions in water, the dissociation of NaCl and of two methane molecules. The memory exhibits multiple time scales and, for NaCl, pronounced oscillatory components. Simulations of the GLE by Markovian embedding techniques accurately reproduce the pair-reaction kinetics from MD simulations without any fitting parameters, which confirms the accuracy of the approximative form of the GLE and of the parameter extraction techniques. By modification of the GLE parameters, we investigate the relative importance of memory, mass, and potential shape effects. Neglect of memory slows down NaCl and methane dissociation by roughly a factor of 2; neglect of mass accelerates reactions by a similar factor, and the harmonic approximation of the potential shape gives rise to slight acceleration. This partial error cancellation explains why Kramers’ theory, which neglects memory effects and treats the potential shape in harmonic approximation, describes reaction rates better than more sophisticated theories. In essence, all three effects, friction memory, inertia, and the potential shape nonharmonicity, are important to quantitatively describe pair-reaction kinetics in water

Butane dihedral angle dynamics in water is dominated by internal friction

The dihedral dynamics of butane in water is known to be rather insensitive to the water viscosity; possible explanations for this involve inertial effects or Kramers’ turnover, the finite memory time of friction, and the presence of so-called internal friction. To disentangle these factors, we introduce a method to directly extract the friction memory function from unconstrained simulations in the presence of an arbitrary free-energy landscape. By analysis of the dihedral friction in butane for varying water viscosity, we demonstrate the existence of an internal friction contribution that does not scale linearly with water viscosity. At normal water viscosity, the internal friction turns out to be eight times larger than the solvent friction and thus completely dominates the effective friction. By comparison with simulations of a constrained butane molecule that has the dihedral as the only degree of freedom, we show that internal friction comes from the six additional degrees of freedom in unconstrained butane that are orthogonal to the dihedral angle reaction coordinate. While the insensitivity of butane’s dihedral dynamics to water viscosity is solely due to the presence of internal friction, inertial effects nevertheless crucially influence the resultant transition rates. In contrast, non-Markovian effects due to the finite memory time are present but do not significantly influence the dihedral barrier-crossing rate of butane. These results not only settle the character of dihedral dynamics in small solvated molecular systems such as butane, they also have important implications for the folding of polymers and proteins

Water Structuring Induces Nonuniversal Hydration Repulsion between Polar Surfaces: Quantitative Comparison between Molecular Simulations, Theory, and Experiments

Polar surfaces in water typically repel each other at close separations, even if they are charge-neutral. This so-called hydration repulsion balances the van der Waals attraction and gives rise to a stable nanometric water layer between the polar surfaces. The resulting hydration water layer is crucial for the properties of concentrated suspensions of lipid membranes and hydrophilic particles in biology and technology, but its origin is unclear. It has been suggested that surface-induced molecular water structuring is responsible for the hydration repulsion, but a quantitative proof of this water-structuring hypothesis is missing. To gain an understanding of the mechanism causing hydration repulsion, we perform molecular simulations of different planar polar surfaces in water. Our simulated hydration forces between phospholipid bilayers agree perfectly with experiments, validating the simulation model and methods. For the comparison with theory, it is important to split the simulated total surface interaction force into a direct contribution from surface–surface molecular interactions and an indirect water-mediated contribution. We find the indirect hydration force and the structural water-ordering profiles from the simulations to be in perfect agreement with the predictions from theoretical models that account for the surface-induced water ordering, which strongly supports the water-structuring hypothesis for the hydration force. However, the comparison between the simulations for polar surfaces with different headgroup architectures reveals significantly different decay lengths of the indirect water-mediated hydration-force, which for laterally homogeneous water structuring would imply different bulk-water properties. We conclude that laterally inhomogeneous water ordering, induced by laterally inhomogeneous surface structures, shapes the hydration repulsion between polar surfaces in a decisive manner. Thus, the indirect water-mediated part of the hydration repulsion is caused by surface-induced water structuring but is surface-specific and thus nonuniversal

Proton-transfer spectroscopy beyond the normal-mode scenario

A stochastic theory is developed to predict the spectral signature of proton-transfer processes and is applied to infrared spectra computed from ab initio molecular-dynamics simulations of a single H5O2+ cation. By constraining the oxygen atoms to a fixed distance, this system serves as a tunable model for general proton-transfer processes with variable barrier height. Three spectral contributions at distinct frequencies are identified and analytically predicted: the quasi-harmonic motion around the most probable configuration, amenable to normal-mode analysis, the contribution due to transfer paths when the proton moves over the barrier, and a shoulder for low frequencies stemming from the stochastic transfer-waiting-time distribution; the latter two contributions are not captured by normal-mode analysis but exclusively reported on the proton-transfer kinetics. In accordance with reaction rate theory, the transfer-waiting-contribution frequency depends inversely exponentially on the barrier height, whereas the transfer-path-contribution frequency is rather insensitive to the barrier height

Non-Markovian modeling of protein folding.

We extract the folding free energy landscape and the time-dependent friction function, the two ingredients of the generalized Langevin equation (GLE), from explicit-water molecular dynamics (MD) simulations of the α-helix forming polypeptide [Formula: see text] for a one-dimensional reaction coordinate based on the sum of the native H-bond distances. Folding and unfolding times from numerical integration of the GLE agree accurately with MD results, which demonstrate the robustness of our GLE-based non-Markovian model. In contrast, Markovian models do not accurately describe the peptide kinetics and in particular, cannot reproduce the folding and unfolding kinetics simultaneously, even if a spatially dependent friction profile is used. Analysis of the GLE demonstrates that memory effects in the friction significantly speed up peptide folding and unfolding kinetics, as predicted by the Grote-Hynes theory, and are the cause of anomalous diffusion in configuration space. Our methods are applicable to any reaction coordinate and in principle, also to experimental trajectories from single-molecule experiments. Our results demonstrate that a consistent description of protein-folding dynamics must account for memory friction effects

Integrative determination of the atomic structure of mutant huntingtin exon 1 fibrils from Huntington's disease

Neurodegeneration in Huntington's disease (HD) is accompanied by the aggregation of fragments of the mutant huntingtin protein, a biomarker of disease progression. A particular pathogenic role has been attributed to the aggregation-prone huntingtin exon 1 (HttEx1) fragment, whose polyglutamine (polyQ) segment is expanded. Unlike amyloid fibrils from Parkinson's and Alzheimer's diseases, the atomic-level structure of HttEx1 fibrils has remained unknown, limiting diagnostic and treatment efforts. We present and analyze the structure of fibrils formed by polyQ peptides and polyQ-expanded HttEx1. Atomic-resolution perspectives are enabled by an integrative analysis and unrestrained all-atom molecular dynamics (MD) simulations incorporating experimental data from electron microscopy (EM), solid-state NMR, and other techniques. Visualizing the HttEx1 subdomains in atomic detail helps explaining the biological properties of these protein aggregates, as well as paves the way for targeting them for detection and degradation.</p

Water structuring induces nonuniversal hydration repulsion between polar surfaces : quantitative comparison between molecular simulations, theory, and experiments

Polar surfaces in water typically repel each other at close separations, even if they are charge-neutral. This so-called hydration repulsion balances the van der Waals attraction and gives rise to a stable nanometric water layer between the polar surfaces. The resulting hydration water layer is crucial for the properties of concentrated suspensions of lipid membranes and hydrophilic particles in biology and technology, but its origin is unclear. It has been suggested that surface-induced molecular water structuring is responsible for the hydration repulsion, but a quantitative proof of this water-structuring hypothesis is missing. To gain an understanding of the mechanism causing hydration repulsion, we perform molecular simulations of different planar polar surfaces in water. Our simulated hydration forces between phospholipid bilayers agree perfectly with experiments, validating the simulation model and methods. For the comparison with theory, it is important to split the simulated total surface interaction force into a direct contribution from surface-surface molecular interactions and an indirect water-mediated contribution. We find the indirect hydration force and the structural water-ordering profiles from the simulations to be in perfect agreement with the predictions from theoretical models that account for the surface-induced water ordering, which strongly supports the water-structuring hypothesis for the hydration force. However, the comparison between the simulations for polar surfaces with different headgroup architectures reveals significantly different decay lengths of the indirect water-mediated hydration-force, which for laterally homogeneous water structuring would imply different bulk-water properties. We conclude that laterally inhomogeneous water ordering, induced by laterally inhomogeneous surface structures, shapes the hydration repulsion between polar surfaces in a decisive manner. Thus, the indirect water-mediated part of the hydration repulsion is caused by surface-induced water structuring but is surface-specific and thus nonuniversal.Deutsche ForschungsgemeinschaftStuttgart Center for Simulation Science (SimTech)Slovenian Research and Innovation Agency ARI