Generalized Langevin equation with fractional Gaussian noise: subdiffusion within a single protein molecule

Main role of fractal-like nature of conformational space in subdiffusion in proteins

Protein dynamics involves a myriad of mechanical movements happening at different time and space scales, which make it highly complex. One of the less understood features of protein dynamics is subdiffusivity, defined as sublinear dependence between displacement and time. Here, we use all-atoms molecular dynamics (MD) simulations to directly interrogate an already well-established theory and demonstrate that subdiffusivity arises from the fractal nature of the network of metastable conformations over which the dynamics, thought of as a diffusion process, takes place.

Internal friction as a factor in the anomalous chain length dependence of DNA transcriptional dynamics

Recent experiments by Brückner et al. [Science 380, 1357 (2023)] have observed an anomalous chain length dependence of the time of near approach of widely separated pairs of genomic elements on transcriptionally active chromosomal DNA. In this paper, I suggest that the anomaly may have its roots in internal friction between neighboring segments on the DNA backbone. The basis for this proposal is a model of chain dynamics formulated in terms of a continuum scaled Brownian walk (sBw) of polymerization index N. The sBw is an extension of the simple Brownian walk model widely used in path integral calculations of polymer properties, differing from it in containing an additional parameter H (the Hurst index) that can be tuned to produce varying degrees of correlation between adjacent monomers. A calculation using the sBw of the mean time τc for chain closure predicts-under the Wilemski-Fixman approximation for diffusion-controlled reactions-that at early times, τc varies as the 2/3 power of N, in close agreement with the findings of the Brückner et al. study. Other scaling relations of that study, including those related to the probability of loop formation and the mean square displacements of terminal monomers, are also satisfactorily accounted for by the model.

Work distribution of a colloid in an elongational flow field and under Ornstein-Uhlenbeck noise

The study of thermodynamic properties of microscopic systems, such as a colloid in a fluid, has been of great interest to researchers since the discovery of the fluctuation theorem and associated laws of stochastic thermodynamics. However, most of these studies confine themselves to systems where effective fluctuations acting on the colloid are in the form of delta-correlated Gaussian white noise (GWN). In this study, instead, we look into the work distribution function when a colloid trapped in a harmonic potential moves from one position to another in a fluid medium with an elongational flow field where the effective fluctuations are given by the Ornstein-Uhlenbeck noise, a type of colored noise. We use path integrals to calculate this distribution function and compare and contrast its properties to the case with GWN. We find that the work distribution function turns out to be non-Gaussian as a result of the elongational flow field but continues to obey the fluctuation theorem in both types of noise. Further, we also look into the effects of the various system parameters on the behavior of work fluctuations and find that although the distribution tends to broaden with increasing noise intensity, increased correlation in fluctuations acts to oppose this effect. Additionally, the system is found to consume heat from the surroundings at early times and dissipate it into the media at later times. This study, therefore, is a step towards gaining a better understanding of the thermodynamic properties of colloidal systems under nonlinear complex flows that also display correlated fluctuations.

Analysis of Velocity Autocorrelation Functions from Molecular Dynamics Simulations of a Small Peptide by the Generalized Langevin Equation with a Power-Law Kernel

Internal motions play an essential role in the biological functions of proteins and have been the subject of numerous theoretical and spectroscopic studies. Such complex environments are associated with anomalous diffusion where, in contrast to the classical Brownian motion, the relevant correlation functions have power law decays with time. In this work, we investigate the presence of long memory stochastic processes through the analysis of atomic velocity autocorrelation functions. Analytical expressions of the velocity autocorrelation function spectrum obtained through a Mori-Zwanzig projection approach were shown to be compatible with molecular dynamics simulations of a small helical peptide (8-polyalanine).

Ballistic Brownian Motion of Nanoconfined DNA

Theoretical treatments of polymer dynamics in liquid generally start with the basic assumption that motion at the smallest scale is heavily overdamped; therefore, inertia can be neglected. We report on the Brownian motion of tethered DNA under nanoconfinement, which was analyzed by molecular dynamics simulation and nanoelectrochemistry-based single-electron shuttle experiments. Our results show a transition into the ballistic Brownian motion regime for short DNA in sub-5 nm gaps, with quality coefficients as high as 2 for double-stranded DNA, an effect mainly attributed to a drastic increase in stiffness. The possibility for DNA to enter the underdamped regime could have profound implications on our understanding of the energetics of biomolecular engines such as the replication machinery, which operates in nanocavities that are a few nanometers wide.

Trends in Single-Molecule Total Internal Reflection Fluorescence Imaging and Their Biological Applications with Lab-on-a-Chip Technology

Single-molecule imaging technologies, especially those based on fluorescence, have been developed to probe both the equilibrium and dynamic properties of biomolecules at the single-molecular and quantitative levels. In this review, we provide an overview of the state-of-the-art advancements in single-molecule fluorescence imaging techniques. We systematically explore the advanced implementations of in vitro single-molecule imaging techniques using total internal reflection fluorescence (TIRF) microscopy, which is widely accessible. This includes discussions on sample preparation, passivation techniques, data collection and analysis, and biological applications. Furthermore, we delve into the compatibility of microfluidic technology for single-molecule fluorescence imaging, highlighting its potential benefits and challenges. Finally, we summarize the current challenges and prospects of fluorescence-based single-molecule imaging techniques, paving the way for further advancements in this rapidly evolving field.

Time-fractional Caputo derivative versus other integrodifferential operators in generalized Fokker-Planck and generalized Langevin equations

Fractional diffusion and Fokker-Planck equations are widely used tools to describe anomalous diffusion in a large variety of complex systems. The equivalent formulations in terms of Caputo or Riemann-Liouville fractional derivatives can be derived as continuum limits of continuous-time random walks and are associated with the Mittag-Leffler relaxation of Fourier modes, interpolating between a short-time stretched exponential and a long-time inverse power-law scaling. More recently, a number of other integrodifferential operators have been proposed, including the Caputo-Fabrizio and Atangana-Baleanu forms. Moreover, the conformable derivative has been introduced. We study here the dynamics of the associated generalized Fokker-Planck equations from the perspective of the moments, the time-averaged mean-squared displacements, and the autocovariance functions. We also study generalized Langevin equations based on these generalized operators. The differences between the Fokker-Planck and Langevin equations with different integrodifferential operators are discussed and compared with the dynamic behavior of established models of scaled Brownian motion and fractional Brownian motion. We demonstrate that the integrodifferential operators with exponential and Mittag-Leffler kernels are not suitable to be introduced to Fokker-Planck and Langevin equations for the physically relevant diffusion scenarios discussed in our paper. The conformable and Caputo Langevin equations are unveiled to share similar properties with scaled and fractional Brownian motion, respectively.

Nonergodic Brownian oscillator: High-frequency response

We consider a Brownian oscillator whose coupling to the environment may lead to the formation of a localized normal mode. For lower values of the oscillator's natural frequency ω≤ω_{c}, the localized mode is absent and the unperturbed oscillator reaches thermal equilibrium. For higher values of ω>ω_{c} when the localized mode is formed, the unperturbed oscillator does not thermalize but rather evolves into a nonequilibrium cyclostationary state. We consider the response of such an oscillator to an external periodic force. Despite the coupling to the environment, the oscillator shows the unbounded resonance (with the response linearly increasing with time) when the frequency of the external force coincides with the frequency of the localized mode. An unusual resonance ("quasiresonance") occurs for the oscillator with the critical value of the natural frequency ω=ω_{c}, which separates thermalizing (ergodic) and nonthermalizing (nonergodic) configurations. In that case, the resonance response increases with time sublinearly, which can be interpreted as a resonance between the external force and the incipient localized mode.

Estimating and Assessing Differential Equation Models with Time-Course Data

Ordinary differential equation (ODE) models are widely used to describe chemical or biological processes. This Article considers the estimation and assessment of such models on the basis of time-course data. Due to experimental limitations, time-course data are often noisy, and some components of the system may not be observed. Furthermore, the computational demands of numerical integration have hindered the widespread adoption of time-course analysis using ODEs. To address these challenges, we explore the efficacy of the recently developed MAGI (MAnifold-constrained Gaussian process Inference) method for ODE inference. First, via a range of examples we show that MAGI is capable of inferring the parameters and system trajectories, including unobserved components, with appropriate uncertainty quantification. Second, we illustrate how MAGI can be used to assess and select different ODE models with time-course data based on MAGI's efficient computation of model predictions. Overall, we believe MAGI is a useful method for the analysis of time-course data in the context of ODE models, which bypasses the need for any numerical integration.

Explicit analytical form for memory kernel in the generalized Langevin equation for end-to-end vector of Rouse chains

The generalized Langevin equation (GLE) provides an attractive theoretical framework for investigating the dynamics of conformational fluctuations of polymeric systems. While the memory kernel is a central function in the GLE, explicit analytical forms for this function have been challenging to obtain, even for the simple models of polymer dynamics. Here, we achieve an explicit analytical expression for the memory kernel in the GLE for the end-to-end vector of Rouse chains in the overdamped limit. Our derivation takes advantage of the finding that the dynamics of the end-to-end vector of Rouse chains with both free ends are equivalent to those of Rouse chains with one free end and the other fixed. For the latter model, we first show that the equations of motion of the Rouse modes as well as their statistical properties can be obtained under the boundary conditions where the free end is held fixed temporarily. We then analytically solve the terms associated with intrachain interactions in the GLE. By formally comparing these terms with the GLE based on the Rouse modes, we obtain an explicit expression for the memory kernel, along with analytical forms for the potential field and the random colored noise force. Our analytical memory kernel is confirmed by numerical calculations in the Laplace space and is shown to yield asymptotic behaviors that are consistent with previous studies. Finally, we utilize our analytical result to simulate the cyclization dynamics of Rouse chains and discuss the scaling of the cyclization time with chain length.

Everlasting impact of initial perturbations on first-passage times of non-Markovian random walks

Persistence, defined as the probability that a signal has not reached a threshold up to a given observation time, plays a crucial role in the theory of random processes. Often, persistence decays algebraically with time with non trivial exponents. However, general analytical methods to calculate persistence exponents cannot be applied to the ubiquitous case of non-Markovian systems relaxing transiently after an imposed initial perturbation. Here, we introduce a theoretical framework that enables the non-perturbative determination of persistence exponents of Gaussian non-Markovian processes with non stationary dynamics relaxing to a steady state after an initial perturbation. Two situations are analyzed: either the system is subjected to a temperature quench at initial time, or its past trajectory is assumed to have been observed and thus known. Our theory covers the case of spatial dimension higher than one, opening the way to characterize non-trivial reaction kinetics for complex systems with non-equilibrium initial conditions.

Role of ergodicity, aging, and Gaussianity in resolving the origins of biomolecule subdiffusion

The internal motions of biomolecules are essential to their function. Although biological macromolecules conventionally show subdiffusive dynamics, only recently has subdiffusion been associated with non-ergodicity. These findings have stimulated new questions in biophysics and statistical mechanics. Is non-ergodic subdiffusion a general strategy shared by biomolecules? What underlying mechanisms are responsible for it? Here, we performed extensive molecular dynamics (MD) simulations to characterize the internal dynamics of six different biomolecules, ranging from single or double-stranded DNA, a single domain protein (KRAS), two globular proteins (PGK and SHP2), to an intrinsically disordered protein (SNAP-25). We found that the subdiffusive behavior of these biomolecules falls into two classes. The internal motion of the first three cases is ergodic subdiffusion and can be interpreted by fractional Brownian motion (FBM), while the latter three cases involve non-ergodic subdiffusion and can be modeled by mixed origins of continuous-time random walk (CTRW) and FBM.

Different Stochastic Resonances Induced by Multiplicative Polynomial Trichotomous Noise in a Fractional Order Oscillator with Time Delay and Fractional Gaussian Noise

A general investigation on the mechanism of stochastic resonance is reported in a time-delay fractional Langevin system, which endues a nonlinear form multiplicative colored noise and fractional Gaussian noise. In terms of theoretical analysis, both the expressions of output steady-state amplitude and that of the first moment of system response are obtained by utilizing stochastic averaging method, fractional Shapiro and Laplace methods. Due to the presence of trichotomous colored noise, the excitation frequency can induce fresh multimodal Bona fide stochastic resonance, exhibiting much more novel dynamical behaviors than the non-disturbance case. It is verified that multimodal pattern only appears with small noise switching rate and memory damping order. The explicit expressions of critical noise intensity corresponding to the generalized stochastic resonance are given for the first time, by which it is determined that nonlinear form colored noise induces much more of a comprehensive resonant phenomena than the linear form. In the case of slow transfer rate noise, a newfangled phenomenon of double hypersensitive response induced by a variation in noise intensity is discovered and verified for the first time, with the necessary range of parameters for this phenomenon given. In terms of numerical scheme, an efficient and feasible algorithm for generating trichotomous noise is proposed, by which an algorithm based on the Caputo fractional derivative are applied. The numerical results match well with the analytical ones.

Multiscale relaxation dynamics and diffusion of myelin basic protein in solution studied by quasielastic neutron scattering

We report an analysis of high-resolution quasielastic neutron scattering spectra from Myelin Basic Protein (MBP) in solution, comparing the spectra at three different temperatures (283, 303, and 323 K) for a pure D₂O buffer and a mixture of D₂O buffer with 30% of deuterated trifluoroethanol (TFE). Accompanying experiments with dynamic light scattering and Circular Dichroism (CD) spectroscopy have been performed to obtain, respectively, the global diffusion constant and the secondary structure content of the molecule for both buffers as a function of temperature. Modeling the decay of the neutron intermediate scattering function by the Mittag-Leffler relaxation function, ϕ(t) = E_α(-(t/τ)^α) (0 < α < 1), we find that trifluoroethanol slows down the relaxation dynamics of the protein at 283 K and leads to a broader relaxation rate spectrum. This effect vanishes with increasing temperature, and at 323 K, its relaxation dynamics is identical in both solvents. These results are coherent with the data from dynamic light scattering, which show that the hydrodynamic radius of MBP in TFE-enriched solutions does not depend on temperature and is only slightly smaller compared to the pure D₂O buffer, except for 283 K, where it is much reduced. In accordance with these observations, the CD spectra reveal that TFE induces essentially a partial transition from β-strands to α-helices, but only a weak increase in the total secondary structure content, leaving about 50% of the protein unfolded. The results show that MBP is for all temperatures and in both buffers an intrinsically disordered protein and that TFE essentially induces a reduction in its hydrodynamic radius and its relaxation dynamics at low temperatures.

Non-ergodicity of a globular protein extending beyond its functional timescale

Internal motions of folded proteins have been assumed to be ergodic, i.e., that the dynamics of a single protein molecule averaged over a very long time resembles that of an ensemble. Here, by performing single-molecule fluorescence resonance energy transfer (smFRET) experiments and molecular dynamics (MD) simulations of a multi-domain globular protein, cytoplasmic protein-tyrosine phosphatase (SHP2), we demonstrate that the functional inter-domain motion is observationally non-ergodic over the time spans 10⁻¹² to 10⁻⁷ s and 10⁻¹ to 10² s. The difference between observational non-ergodicity and simple non-convergence is discussed. In comparison, a single-strand DNA of similar size behaves ergodically with an energy landscape resembling a one-dimensional linear chain. The observed non-ergodicity results from the hierarchical connectivity of the high-dimensional energy landscape of the protein molecule. As the characteristic time for the protein to conduct its dephosphorylation function is ∼10 s, our findings suggest that, due to the non-ergodicity, individual, seemingly identical protein molecules can be dynamically and functionally different.

Internal motions of folded proteins have been assumed to be ergodic, i.e., that the dynamics of a single protein molecule averaged over a very long time resembles that of an ensemble.

Particle dynamics in viscoelastic media: Effects of non-thermal white noise on barrier crossing rates

The growing interest in the dynamics of self-driven particle motion has brought increased attention to the effects of non-thermal noise on condensed phase diffusion. Thanks to data recently collected by Ferrer et al. on activated dynamics in the presence of memory [Phys. Rev. Lett. 126, 108001 (2021)], some of these effects can now be characterized quantitatively. In the present paper, the data collected by Ferrer et al. are used to calculate the extent to which non-thermal white noise alters the time taken by single micron-sized silica particles in a viscoelastic medium to cross the barrier separating the two wells of an optically created bistable potential. The calculation-based on a generalized version of Kramers's flux-over-population approach-indicates that the added noise causes the barrier crossing rate (compared to the noise-free case) to first increase as a function of the noise strength and then to plateau to a constant value. The precise degree of rate enhancement may depend on how the data from the experiments conducted by Ferrer et al. are used in the flux-over-population approach. As claimed by Ferrer et al., this approach predicts barrier crossing times for the original silica-fluid system that agree almost perfectly with their experimental counterparts. However, this near-perfect agreement between theory and experiment is only achieved if the theoretical crossing times are obtained from the most probable values of a crossing time distribution constructed from the distributions of various parameters in Kramers's rate expression. If the mean values of these parameters are used in the expression instead, as would be commonly done, the theoretical crossing times are found to be as much as 1.5 times higher than the experimental values. However, these times turn out to be consistent with an alternative model of viscoelastic barrier crossing based on a mean first passage time formalism, which also uses mean parameter values in its rate expression. The rate enhancements predicted for barrier crossing under non-thermal noise are based on these mean parameter values and are open to experimental verification.

Subdiffusive-Brownian crossover in membrane proteins: a generalized Langevin equation-based approach

In this work, we propose a generalized Langevin equation-based model to describe the lateral diffusion of a protein in a lipid bilayer. The memory kernel is represented in terms of a viscous (instantaneous) and an elastic (noninstantaneous) component modeled through a Dirac δ function and a three-parameter Mittag-Leffler type function, respectively. By imposing a specific relationship between the parameters of the three-parameter Mittag-Leffler function, the different dynamical regimes-namely ballistic, subdiffusive, and Brownian, as well as the crossover from one regime to another-are retrieved. Within this approach, the transition time from the ballistic to the subdiffusive regime and the spectrum of relaxation times underlying the transition from the subdiffusive to the Brownian regime are given. The reliability of the model is tested by comparing the mean-square displacement derived in the framework of this model and the mean-square displacement of a protein diffusing in a membrane calculated through molecular dynamics simulations.

Classification of particle trajectories in living cells: Machine learning versus statistical testing hypothesis for fractional anomalous diffusion

Single-particle tracking (SPT) has become a popular tool to study the intracellular transport of molecules in living cells. Inferring the character of their dynamics is important, because it determines the organization and functions of the cells. For this reason, one of the first steps in the analysis of SPT data is the identification of the diffusion type of the observed particles. The most popular method to identify the class of a trajectory is based on the mean-square displacement (MSD). However, due to its known limitations, several other approaches have been already proposed. With the recent advances in algorithms and the developments of modern hardware, the classification attempts rooted in machine learning (ML) are of particular interest. In this work, we adopt two ML ensemble algorithms, i.e., random forest and gradient boosting, to the problem of trajectory classification. We present a new set of features used to transform the raw trajectories data into input vectors required by the classifiers. The resulting models are then applied to real data for G protein-coupled receptors and G proteins. The classification results are compared to recent statistical methods going beyond MSD.

A near analytic solution of a stochastic immune response model considering variability in virus and T-cell dynamics

Biological processes at the cellular level are stochastic in nature, and the immune response system is no different. Therefore, models that attempt to explain this system need to also incorporate noise or fluctuations that can account for the observed variability. In this work, a stochastic model of the immune response system is presented in terms of the dynamics of T cells and virus particles. Making use of the Green's function and the Wilemski-Fixman approximation, this model is then solved to obtain the analytical expression for the joint probability density function of these variables in the early and late stages of infection. This is then also used to calculate the average level of virus particles in the system. Upon comparing the theoretically predicted average virus levels to those of COVID-19 patients, it is hypothesized that the long-lived dynamics that are characteristics of such viral infections are due to the long range correlations in the temporal fluctuations of the virions. This model, therefore, provides an insight into the effects of noise on viral dynamics.

Microscopic insights into dynamic disorder in the isomerization dynamics of the protein BPTI

Understanding the dynamic disorder behind a process, i.e., the dynamic effect of fluctuations that occur on a timescale slower or comparable with the timescale of the process, is essential for elucidating the dynamics and kinetics of complicated molecular processes in biomolecules and liquids. Despite numerous theoretical studies of single-molecule kinetics, our microscopic understanding of dynamic disorder remains limited. In the present study, we investigate the microscopic aspects of dynamic disorder in the isomerization dynamics of the Cys14-Cys38 disulfide bond in the protein bovine pancreatic trypsin inhibitor, which has been observed by nuclear magnetic resonance. We use a theoretical model with a stochastic transition rate coefficient, which is calculated from the 1-ms-long time molecular dynamics trajectory obtained by Shaw et al. [Science 330, 341-346 (2010)]. The isomerization dynamics are expressed by the transitions between coarse-grained states consisting of internal states, i.e., conformational sub-states. In this description, the rate for the transition from the coarse-grained states is stochastically modulated due to fluctuations between internal states. We examine the survival probability for the conformational transitions from a coarse-grained state using a theoretical model, which is a good approximation to the directly calculated survival probability. The dynamic disorder changes from a slow modulation limit to a fast modulation limit depending on the aspects of the coarse-grained states. Our analysis of the rate modulations behind the survival probability, in relation to the fluctuations between internal states, reveals the microscopic origin of dynamic disorder.

Statistical moments of the interspike intervals for a neuron model driven by trichotomous noise

The influence of a colored three-level input noise (trichotomous noise) on the spike generation of a perfect integrate-and-fire (PIF) model of neurons is studied. Using a first-passage-time formulation, exact expressions for the Laplace transform of the output interspike interval (ISI) density and for the statistical moments of the ISIs (such as the coefficient of variation, the skewness, the serial correlation coefficient, and the Fano factor) are derived. To model the anomalous subdiffusion that can arise from, e.g., the trapping properties of dendritic spines, the model is extended by including a random operational time in the form of an inverse strictly increasing Lévy-type subordinator, and exact formulas for ISI statistics are given for this case as well. Particularly, it is shown that at some parameter regimes, the ISI density exhibits a three-modal structure. The results for the extended model show that the ISI serial correlation coefficient and the Fano factor are nonmonotonic with respect to the input current, which indicates that at an intermediate value of the input current the variability of the output spike trains is minimal. Similarities and differences between the behavior of the presented models and the previously investigated PIF models driven by dichotomous noise are also discussed.

Characterizing stochastic resonance in a triple cavity

Many biological systems possess confined structures, which produce novel influences on the dynamics. Here, stochastic resonance (SR) in a triple cavity that consists of three units and is subjected to noise, periodic force and vertical constance force is studied, by calculating the spectral amplification η numerically. Meanwhile, SR in the given triple cavity and differences from other structures are explored. First, it is found that the cavity parameters can eliminate or regulate the maximum of η and the noise intensity that induces this maximum. Second, compared to a double cavity with similar maximum/minimum widths and distances between two maximum widths as the triple cavity, η in the triple one shows a larger maximum. Next, the conversion of the natural boundary in the pure potential to the reflection boundary in the triple cavity will create the necessity of a vertical force to induce SR and lead to a decrease in the maximum of η. In addition, η monotonically decreases with the increase of the vertical force and frequency of the periodic force, while it presents several trends when increasing the periodic force's amplitude for different noise intensities. Finally, our studies are extended to the impact of fractional Gaussian noise excitations. This article is part of the theme issue 'Vibrational and stochastic resonance in driven nonlinear systems (part 2)'.

Petrov-Galerkin methods for the construction of non-Markovian dynamics preserving nonlocal statistics

A common observation in coarse-graining a molecular system is the non-Markovian behavior, primarily due to the lack of scale separations. This is reflected in the strong memory effect and the non-white noise spectrum, which must be incorporated into a coarse-grained description to correctly predict dynamic properties. To construct a stochastic model that gives rise to the correct non-Markovian dynamics, we propose a Galerkin projection approach, which transforms the exhausting effort of finding an appropriate model to choosing appropriate subspaces in terms of the derivatives of the coarse-grained variables and, at the same time, provides an accurate approximation to the generalized Langevin equation. We introduce the notion of fractional statistics that embodies nonlocal properties. More importantly, we show how to pick subspaces in the Galerkin projection so that those statistics are automatically matched.

Effect of Memory and Inertial Contribution on Transition-Time Distributions: Theory and Simulations

Transition paths refer to the time taken by molecules to cross a barrier separating two molecular conformations. In this work, we study how memory, as well as inertial contribution in the dynamics along a reaction coordinate, can affect the distribution of the transition-path time. We use a simple model of dynamics governed by a generalized Langevin equation with a power-law memory along with the inertial term, which was neglected in previous studies, where memory effects were explored only in the overdamped limit. We derive an approximate expression for the transit-time distribution and discuss our results for the short- and long-time limits and also compare it with known results in the high friction (overdamped) limit as well as in the Markovian limit. We have developed a numerical algorithm to test our theoretical results against extensive numerical simulations.

Inference of dynamic systems from noisy and sparse data via manifold-constrained Gaussian processes

Parameter estimation for nonlinear dynamic system models, represented by ordinary differential equations (ODEs), using noisy and sparse data, is a vital task in many fields. We propose a fast and accurate method, manifold-constrained Gaussian process inference (MAGI), for this task. MAGI uses a Gaussian process model over time series data, explicitly conditioned on the manifold constraint that derivatives of the Gaussian process must satisfy the ODE system. By doing so, we completely bypass the need for numerical integration and achieve substantial savings in computational time. MAGI is also suitable for inference with unobserved system components, which often occur in real experiments. MAGI is distinct from existing approaches as we provide a principled statistical construction under a Bayesian framework, which incorporates the ODE system through the manifold constraint. We demonstrate the accuracy and speed of MAGI using realistic examples based on physical experiments.

Heterogeneous diffusion processes and nonergodicity with Gaussian colored noise in layered diffusivity landscapes

Heterogeneous diffusion processes (HDPs) with space-dependent diffusion coefficients D(x) are found in a number of real-world systems, such as for diffusion of macromolecules or submicron tracers in biological cells. Here, we examine HDPs in quenched-disorder systems with Gaussian colored noise (GCN) characterized by a diffusion coefficient with a power-law dependence on the particle position and with a spatially random scaling exponent. Typically, D(x) is considered to be centerd at the origin and the entire x axis is characterized by a single scaling exponent α. In this work we consider a spatially random scenario: in periodic intervals ("layers") in space D(x) is centerd to the midpoint of each interval. In each interval the scaling exponent α is randomly chosen from a Gaussian distribution. The effects of the variation of the scaling exponents, the periodicity of the domains ("layer thickness") of the diffusion coefficient in this stratified system, and the correlation time of the GCN are analyzed numerically in detail. We discuss the regimes of superdiffusion, subdiffusion, and normal diffusion realisable in this system. We observe and quantify the domains where nonergodic and non-Gaussian behaviors emerge in this system. Our results provide new insights into the understanding of weak ergodicity breaking for HDPs driven by colored noise, with potential applications in quenched layered systems, typical model systems for diffusion in biological cells and tissues, as well as for diffusion in geophysical systems.

Impact of Feature Choice on Machine Learning Classification of Fractional Anomalous Diffusion

The growing interest in machine learning methods has raised the need for a careful study of their application to the experimental single-particle tracking data. In this paper, we present the differences in the classification of the fractional anomalous diffusion trajectories that arise from the selection of the features used in random forest and gradient boosting algorithms. Comparing two recently used sets of human-engineered attributes with a new one, which was tailor-made for the problem, we show the importance of a thoughtful choice of the features and parameters. We also analyse the influence of alterations of synthetic training data set on the classification results. The trained classifiers are tested on real trajectories of G proteins and their receptors on a plasma membrane.

Origin of subdiffusions in proteins: Insight from peptide systems

Subdiffusive kinetics are popular in proteins and peptides as observed in experiments and simulations. For protein systems with diverse interactions, are there multiple mechanisms to produce the common subdiffusion behavior? To approach this problem, long trajectories of two model peptides are simulated to study the mechanism of subdiffusion and the relations with their interactions. The free-energy profiles and the subdiffusive kinetics are observed for these two peptides. A hierarchical plateau analysis is employed to extract the features of the landscape from the mean square of displacement. The mechanism of subdiffusions can be postulated by comparing the exponents by simulations with those based on various models. The results indicate that the mechanisms of these two peptides are different and are related to the characteristics of their energy landscapes. The subdiffusion of the flexible peptide is mainly caused by depth distribution of traps on the energy landscape, while the subdiffusion of the helical peptide is attributed to the fractal topology of local minima on the landscape. The emergence of these different mechanisms reflects different kinetic scenarios in peptide systems though the peptides behave in a similar way of diffusion. To confirm these ideas, the transition networks between various conformations of these peptides are generated. Based on the network description, the controlled kinetics based only on the topology of the networks are calculated and compared with the results based on simulations. For the flexible peptide, the feature of controlled diffusion is distinct from that of simulation, and for the helical peptide, two kinds of kinetics have a similar exponent of subdiffusion. These results further exemplify the importance of the landscape topology in the kinetics of structural proteins and the effect of depth distribution of traps for the subdiffusion of disordered peptides.

Look at Tempered Subdiffusion in a Conjugate Map: Desire for the Confinement

The Laplace distribution of random processes was observed in numerous situations that include glasses, colloidal suspensions, live cells, and firm growth. Its origin is not so trivial as in the case of Gaussian distribution, supported by the central limit theorem. Sums of Laplace distributed random variables are not Laplace distributed. We discovered a new mechanism leading to the Laplace distribution of observable values. This mechanism changes the contribution ratio between a jump and a continuous parts of random processes. Our concept uses properties of Bernstein functions and subordinators connected with them.

Getting around the cell: physical transport in the intracellular world

Eukaryotic cells face the challenging task of transporting a variety of particles through the complex intracellular milieu in order to deliver, distribute, and mix the many components that support cell function. In this review, we explore the biological objectives and physical mechanisms of intracellular transport. Our focus is on cytoplasmic and intra-organelle transport at the whole-cell scale. We outline several key biological functions that depend on physically transporting components across the cell, including the delivery of secreted proteins, support of cell growth and repair, propagation of intracellular signals, establishment of organelle contacts, and spatial organization of metabolic gradients. We then review the three primary physical modes of transport in eukaryotic cells: diffusive motion, motor-driven transport, and advection by cytoplasmic flow. For each mechanism, we identify the main factors that determine speed and directionality. We also highlight the efficiency of each transport mode in fulfilling various key objectives of transport, such as particle mixing, directed delivery, and rapid target search. Taken together, the interplay of diffusion, molecular motors, and flows supports the intracellular transport needs that underlie a broad variety of biological phenomena.

Effects of transient subordinators on the firing statistics of a neuron model driven by dichotomous noise

The behavior of a stochastic perfect integrate-and-fire (PIF) model of neurons is considered. The effect of temporally correlated random activity of synaptic inputs is modeled as a combination of an asymmetric dichotomous noise and a random operation time in the form of an inverse strictly increasing Lévy-type subordinator. Using a first-passage-time formulation, we find exact expressions for the output interspike interval (ISI) statistics. Particularly, it is shown that at some parameter regimes the ISI density exhibits a multimodal structure. Moreover, it is demonstrated that the coefficient of variation, the serial correlation coefficient, and the Fano factor display a nonmonotonic dependence on the mean input current μ, i.e., the ISI's regularity is maximized at an intermediate value of μ. The features of spike statistics, analytically revealed in our study, are compared with previously obtained results for a perfect integrate-and-fire neuron model driven by dichotomous noise (without subordination).

Mechanical Unfolding of Single Polyubiquitin Molecules Reveals Evidence of Dynamic Disorder

Mechanical unfolding of single polyubiquitin molecules subjected to a constant stretching force showed nonexponentiality in the measured probability density of unfolding (waiting time distribution) and the survival probability of the folded state during the course of the measurements. These observations explored the relevance of disorder present in the system under study with implications for a static disorder approach to rationalize the experimental results. Here, an approach for dynamic disorder is presented based on Zwanzig's fluctuating bottleneck (FB) model, in which the rate of the reaction is controlled by the passage through the cross-sectional area of the bottleneck. The radius of the latter undergoes stochastic fluctuations that in turn is described in terms of the end-to-end distance fluctuations of the Rouse-like dynamics using a non-Markovian generalized Langevin equation with a memory kernel and Gaussian colored noise. Our results are comprised of analytical expressions for the survival probability and waiting time distribution, which show excellent agreement with the experimental data throughout the range of the applied forces. In addition, by fitting the survival probabilities at different stretching forces, we quantify two system parameters, namely, the average free energy ΔG _av and the average distance to the transition state Δx _av, both perfectly recovered the experimental estimates. These agreements validate the present model of polymer dynamics, which captures the very essence of dynamic disorder in single-molecule pulling experiments.

Generalization of Langevin Dynamics from Spatio-Temporal Dressed Dynamics Perspective

Understanding multiscale dynamics characteristics has been the holy grail in a broad range of scientific disciplines from physics, chemistry, to biology, and beyond. The seminal Langevin equation successfully unravels remarkable details of Brownian motion dynamics involving stochastic collisions from the fluidic medium. However, extraordinary phenomena beyond the scope of Langevin dynamics were observed to exhibit a series of multiscale dynamic features in recent years. Here, an explicit spatio-temporal coupled kernel is developed to provide the microscopic account for delicate dynamic coupling between particle and medium in the dressed dynamics perspective. This methodology is applied to investigate an aqueous solvation shell model with an explicit spatial boundary to illustrate the significance of spatial and temporal coupling, leading to a general temporal profile of dressed dynamics over 12 orders of magnitude in time In particular, the time-resolved viscosity is formulated to address the remarkable enhancement of viscosity measured by the liquid cell electron microscopy. The understanding of dressed dynamics processes will be greatly enriched by further studies on the detailed dynamics that manifests the microscopic inhomogeneity of medium.

Unraveling Dynamic Transitions in Time-Resolved Biomolecular Motions by A Dressed Diffusion Model

Recent experimental data reveal the complexity of diffusion dynamics beyond the scope of classical Brownian dynamics. The particles exhibit diverse diffusive motions from the anomalous toward classical diffusion over a wide range of temporal scales. Here a dressed diffusion model is developed to account for non-Brownian phenomena. By coupling the particle dynamics with a local field, the dressed diffusion model generalizes the Langevin equation through coupled damping kernels and generates the salient feature of time-dependent diffusion dynamics reported in the experimental measurements of biomolecules. The dressed diffusion model provides one quantitative aspect for future endeavors in this rapid-growing field.

Transition path dynamics in the binding of intrinsically disordered proteins: A simulation study

Association of proteins and other biopolymers is a ubiquitous process in living systems. Recent single-molecule measurements probe the dynamics of association in unprecedented detail by measuring the properties of association transition paths, i.e., short segments of molecular trajectories between the time the proteins are close enough to interact and the formation of the final complex. Interpretation of such measurements requires adequate models for describing the dynamics of experimental observables. In an effort to develop such models, here we report a simulation study of the association dynamics of two oppositely charged, disordered polymers. We mimic experimental measurements by monitoring intermonomer distances, which we treat as "experimental reaction coordinates." While the dynamics of the distance between the centers of mass of the molecules is found to be memoryless and diffusive, the dynamics of the experimental reaction coordinates displays significant memory and can be described by a generalized Langevin equation with a memory kernel. We compute the most commonly measured property of transition paths, the distribution of the transition path time, and show that, despite the non-Markovianity of the underlying dynamics, it is well approximated as one-dimensional diffusion in the potential of mean force provided that an apparent value of the diffusion coefficient is used. This apparent value is intermediate between the slow (low frequency) and fast (high frequency) limits of the memory kernel. We have further studied how the mean transition path time depends on the ionic strength and found only weak dependence despite strong electrostatic attraction between the polymers.

Non-Markovian harmonic oscillator across a magnetic field and time-dependent force fields

We study the non-Markovian Brownian motion of an electrically charged harmonic oscillator through the action of both a constant magnetic field and time-dependent force fields. The generalized Langevin equation with a friction memory kernel is used to derive the generalized phase-space Fokker-Planck equation for the harmonic oscillator in the absence and in the presence of time-dependent force fields. To achieve our goal, the characteristic function method is applied to obtain, in an accurate way, the theoretical description of the problem. We explicitly calculate the correlation and cross-correlation functions for the position and velocity vectors. We show that the relevant physics behind the theory is contained in the generalized diffusion coefficient, which accounts for the natural coupling between both the harmonic oscillator and magnetic field. Our theoretical results are compared with those previously reported in the literature.

Nanoparticle transport phenomena in confined flows

Nanoparticles submerged in confined flow fields occur in several technological applications involving heat and mass transfer in nanoscale systems. Describing the transport with nanoparticles in confined flows poses additional challenges due to the coupling between the thermal effects and fluid forces. Here, we focus on the relevant literature related to Brownian motion, hydrodynamic interactions and transport associated with nanoparticles in confined flows. We review the literature on the several techniques that are based on the principles of non-equilibrium statistical mechanics and computational fluid dynamics in order to simultaneously preserve the fluctuation-dissipation relationship and the prevailing hydrodynamic correlations. Through a review of select examples, we discuss the treatments of the temporal dynamics from the colloidal scales to the molecular scales pertaining to nanoscale fluid dynamics and heat transfer. As evident from this review, there, indeed has been little progress made in regard to the accurate modeling of heat transport in nanofluids flowing in confined geometries such as tubes. Therefore the associated mechanisms with such processes remain unexplained. This review has revealed that the information available in open literature on the transport properties of nanofluids is often contradictory and confusing. It has been very difficult to draw definitive conclusions. The quality of work reported on this topic is non-uniform. A significant portion of this review pertains to the treatment of the fluid dynamic aspects of the nanoparticle transport problem. By simultaneously treating the energy transport in ways discussed in this review as related to momentum transport, the ultimate goal of understanding nanoscale heat transport in confined flows may be achieved.

Transport and escape in a deformable channel driven by fractional Gaussian noise

Fractional Gaussian noise (FGN) with the Hurst exponent H is an important tool to model various phenomena in biophysical systems, like subdiffusion in a single protein molecule. Considering that there also exists a confined structure which can be modeled as a channel in these systems, transport and escape driven by FGN in a deformable channel are investigated in this paper. By calculating the mean velocity, and the mean first passage time (MFPT) for crossing the nearest bottleneck and the probability distribution of the final position, effects of FGN and channel structure on the system dynamics are illustrated. Our results indicate that FGN has a complex influence mechanism under different combinations of H and the noise intensity. For a persistence case (H>0.5), the mean velocity decreases but MFPT increases with the increase of the noise intensity and H. While for an antipersistence case (H<0.5), when H is small, the relationships among the mean velocity, MFPT and the noise intensity are exactly the opposite to persistence cases. When H has a large value, the mean velocity tends to first decrease and then increase. Moreover, effects of the bottleneck and channel asymmetry are investigated. It is shown that a small H and a large channel width can lead to a large mean velocity and fast crossing. Besides, a channel asymmetry can affect the system dynamics by inducing asymmetric structure and adjusting the width of bottleneck. However, the effect of the bottleneck is the main factor. Therefore, a combination of channel with wide bottleneck and FGN in an antipersistence regime is the optimal choice to promote the transport and escape. These results provide a basis for the explanation of molecular activity in living organisms and the design of particle mixture separators.

Passage through a sub-diffusing geometrical bottleneck

The usual Kramers theory of reaction rates in condensed media predict the rate to have an inverse dependence on the viscosity of the medium, η. However, experiments on ligand binding to proteins, performed long ago, showed the rate to have η^-ν dependence, with ν in the range of 0.4-0.8. Zwanzig [J. Chem. Phys. 97, 3587 (1992)] suggested a model in which the ligand has to pass through a fluctuating opening to reach the binding site. This fluctuating gate model predicted the rate to be proportional to η^-1/2. More recently, experiments performed by Xie et al. [Phys. Rev. Lett. 93, 180603 (2004)] showed that the distance between two groups in a protein undergoes not normal diffusion, but subdiffusion. Hence, in this paper, we suggest and solve a generalization of the Zwanzig model, viz., passage through an opening, whose size undergoes subdiffusion. Our solution shows that the rate is proportional to η^-ν with ν in the range of 0.5-1, and hence, the subdiffusion model can explain the experimental observations.

Fractional electron transfer kinetics and a quantum breaking of ergodicity

The dissipative curve-crossing problem provides a paradigm for electron-transfer (ET) processes in condensed media. It establishes the simplest conceptual test bed to study the influence of the medium's dynamics on ET kinetics both on the ensemble level, and on the level of single particles. Single electron description is particularly important for nanoscaled systems like proteins, or molecular wires. Especially insightful is this framework in the semiclassical limit, where the environment can be treated classically, and an exact analytical treatment becomes feasible. Slow medium's dynamics is capable of enslaving ET and bringing it on the ensemble level from a quantum regime of nonadiabatic tunneling to the classical adiabatic regime, where electrons follow the nuclei rearrangements. This classical adiabatic textbook picture contradicts, however, in a very spectacular fashion to the statistics of single electron transitions, even in the Debye, memoryless media, also named Ohmic in the parlance of the famed spin-boson model. On the single particle level, ET always remains quantum, and this was named a quantum breaking of ergodicity in the adiabatic ET regime. What happens in the case of subdiffusive, fractional, or sub-Ohmic medium's dynamics, which is featured by power-law decaying dynamical memory effects typical, e.g., for protein macromolecules, and other viscoelastic media? Such a memory is vividly manifested by anomalous Cole-Cole dielectric response in such media. We address this question based both on accurate numerics and analytical theory. The ensemble theory remarkably agrees with the numerical dynamics of electronic populations, revealing a power-law relaxation tail even in a profoundly nonadiabatic electron transfer regime. In other words, ET in such media should typically display fractional kinetics. However, a profound difference with the numerically accurate results occurs for the distribution of residence times in the electronic states, both on the ensemble level and the level of single trajectories. Ergodicity is broken dynamically even in a more spectacular way than in the memoryless case. Our results question the applicability of all the existing and widely accepted ensemble theories of electron transfer in fractional, sub-Ohmic environments, on the level of single molecules, and provide a real challenge to face, both for theorists and experimentalists.

Modeling Biological Oscillations: Integration of Short Reaction Pauses into a Stationary Model of a Negative Feedback Loop Generates Sustained Long Oscillations

Sustained oscillations are frequently observed in biological systems consisting of a negative feedback loop, but a mathematical model with two ordinary differential equations (ODE) that has a negative feedback loop structure fails to produce sustained oscillations. Only when a time delay is introduced into the system by expanding to a three-ODE model, transforming to a two-delay differential equations (DDE) model, or introducing a bistable trigger do stable oscillations present themselves. In this study, we propose another mechanism for producing sustained oscillations based on periodic reaction pauses of chemical reactions in a negative feedback system. We model the oscillatory system behavior by allowing the coefficients in the two-ODE model to be periodic functions of time-called pulsate functions-to account for reactions with go-stop pulses. We find that replacing coefficients in the two-ODE system with pulsate functions with microscale (several seconds) pauses can produce stable system-wide oscillations that have periods of approximately 1 to several hours long. We also compare our two-ODE and three-ODE models with the two-DDE, three-ODE, and three-DDE models without the pulsate functions. Our numerical experiments suggest that sustained long oscillations in biological systems with a negative feedback loop may be an intrinsic property arising from the slow diffusion-based pulsate behavior of biochemical reactions.

Cyclization and Relaxation Dynamics of Finite-Length Collapsed Self-Avoiding Polymers

We study the cyclization and relaxation dynamics of ideal as well as interacting polymers as a function of chain length N. For the cyclization time τ_{cyc} of ideal chains we recover the known scaling τ_{cyc}∼N^{2} for different backbone models, for a self-avoiding slightly collapsed chain we obtain from Langevin simulations and scaling theory a modified scaling τ_{cyc}∼N^{5/3}. The cyclization and relaxation dynamics of a finite-length collapsed chain scale differently; this unexpected dynamic multiscale behavior is rationalized by the crossover between swollen and collapsed chain behavior.

Simulation Best Practices for Lipid Membranes [Article v1.0]

We establish a reliable and robust standardization of settings for practical molecular dynamics (MD) simulations of pure and mixed (single- and multi-component) lipid bilayer membranes. In lipid membranes research, particle-based molecular simulations are a powerful tool alongside continuum theory, lipidomics, and model, in vitro, and in vivo experiments. Molecular simulations can provide precise and reproducible spatiotemporal (atomic- and femtosecond-level) information about membrane structure, mechanics, thermodynamics, kinetics, and dynamics. Yet the simulation of lipid membranes can be a daunting task, given the uniqueness of lipid membranes relative to conventional liquid-liquid and solid-liquid interfaces, the immense and complex thermodynamic and statistical mechanical theory, the diversity of multiscale lipid models, limitations of modern computing power, the difficulty and ambiguity of simulation controls, finite size effects, competitive continuum simulation alternatives, and the desired application, including vesicle experiments and biological membranes. These issues can complicate an essential understanding of the field of lipid membranes, and create major bottlenecks to simulation advancement. In this article, we clarify these issues and present a consistent, thorough, and user-friendly framework for the design of state-of-the-art lipid membrane MD simulations. We hope to allow early-career researchers to quickly overcome common obstacles in the field of lipid membranes and reach maximal impact in their simulations.

Adhesion of a polymer-grafted nanoparticle to cells explored using generalized Langevin dynamics

We model a polymer-grafted stealth nanoparticle (SNP) as a composite system consisting of a spherical core coated with a porous polymeric brush with end-ligands. Adjacent to target cells, the near-wall hydrodynamics, thermal fluctuations, and thermodynamic adhesive interactions simultaneously impact the transient motion of the SNP. Employing both the Langevin framework for the effective hard sphere dynamics and the coupled generalized Langevin framework for the nanoparticle-polymer dynamics, we comprehensively investigate the velocity and position temporal relaxations of the SNP in the absence and presence of end-to-end distance fluctuations for the tethered polymer. We demonstrate that polymer structural relaxations substantially impact the SNP adhesive dynamics, especially when the grafted polymer is more flexible. Moreover, a long-time tail with t-3/2 scaling due to polymer chain-length fluctuations is observed in the velocity autocorrelation for a bound SNP. Finally, the thermodynamic effects of membrane morphology on SNP adhesion are explored by modifying the membrane-mediated binding potential of mean force.

Non-Markovian intracellular transport with sub-diffusion and run-length dependent detachment rate

Intracellular transport of organelles is fundamental to cell function and health. The mounting evidence suggests that this transport is in fact anomalous. However, the reasons for the anomaly is still under debate. We examined experimental trajectories of organelles inside a living cell and propose a mathematical model that describes the previously reported transition from sub-diffusive to super-diffusive motion. In order to explain super-diffusive behaviour at long times, we introduce non-Markovian detachment kinetics of the cargo: the rate of detachment is inversely proportional to the time since the last attachment. Recently, we observed the non-Markovian detachment rate experimentally in eukaryotic cells. Here we further discuss different scenarios of how this effective non-Markovian detachment rate could arise. The non-Markovian model is successful in simultaneously describing the time averaged variance (the time averaged mean squared displacement corrected for directed motion), the mean first passage time of trajectories and the multiple peaks observed in the distributions of cargo velocities. We argue that non-Markovian kinetics could be biologically beneficial compared to the Markovian kinetics commonly used for modelling, by increasing the average distance the cargoes travel when a microtubule is blocked by other filaments. In turn, sub-diffusion allows cargoes to reach neighbouring filaments with higher probability, which promotes active motion along the microtubules.

Fluid heterogeneity detection based on the asymptotic distribution of the time-averaged mean squared displacement in single particle tracking experiments

A tracer particle is called anomalously diffusive if its mean squared displacement grows approximately as σ 2 t α as a function of time t for some constant σ 2, where the diffusion exponent satisfies α ≠ 1. In this article, we use recent results on the asymptotic distribution of the time-averaged mean squared displacement [20] to construct statistical tests for detecting physical heterogeneity in viscoelastic fluid samples starting from one or multiple observed anomalously diffusive paths. The methods are asymptotically valid for the range 0 < α < 3/2 and involve a mathematical characterization of time-averaged mean squared displacement bias and the effect of correlated disturbance errors. The assumptions on particle motion cover a broad family of fractional Gaussian processes, including fractional Brownian motion and many fractional instances of the generalized Langevin equation framework. We apply the proposed methods in experimental data from treated P. aeruginosa biofilms generated by the collaboration of the Hill and Schoenfisch Labs at UNC-Chapel Hill.