1
|
Muenker TM, Knotz G, Krüger M, Betz T. Accessing activity and viscoelastic properties of artificial and living systems from passive measurement. NATURE MATERIALS 2024; 23:1283-1291. [PMID: 39085417 DOI: 10.1038/s41563-024-01957-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/29/2022] [Accepted: 06/26/2024] [Indexed: 08/02/2024]
Abstract
Living systems are complex dynamic entities that operate far from thermodynamic equilibrium. Their active, non-equilibrium behaviour requires energy to drive cellular organization and dynamics. Unfortunately, most statistical mechanics approaches are not valid in non-equilibrium situations, forcing researchers to use intricate and often invasive methods to study living processes. Here we experimentally demonstrate that an observable termed mean back relaxation quantifies the active mechanics of living cells from passively observed particle trajectories. The mean back relaxation represents the average trajectory of a particle after a recent motion and is calculated from three-point probabilities. We show that this parameter allows the detection of broken detailed balance in confined systems. We experimentally observe that it provides access to the non-equilibrium generating energy and viscoelastic properties of artificial bulk materials and living cells. These findings suggest that the mean back relaxation can function as a marker of non-equilibrium dynamics and is a non-invasive avenue to determine viscoelastic material properties from passive measurements.
Collapse
Affiliation(s)
- Till M Muenker
- Third Institute of Physics, Georg August Universität Göttingen, Göttingen, Germany
| | - Gabriel Knotz
- Institute of Theoretical Physics, Georg August Universität Göttingen, Göttingen, Germany
| | - Matthias Krüger
- Institute of Theoretical Physics, Georg August Universität Göttingen, Göttingen, Germany.
| | - Timo Betz
- Third Institute of Physics, Georg August Universität Göttingen, Göttingen, Germany.
- Cluster of Excellence 'Multiscale Bioimaging: From Molecular Machines to Networks of Excitable Cells' (MBExC), Georg August Universität Göttingen, Göttingen, Germany.
| |
Collapse
|
2
|
Netz RR. Derivation of the nonequilibrium generalized Langevin equation from a time-dependent many-body Hamiltonian. Phys Rev E 2024; 110:014123. [PMID: 39160956 DOI: 10.1103/physreve.110.014123] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/01/2023] [Accepted: 06/20/2024] [Indexed: 08/21/2024]
Abstract
It has become standard practice to describe systems that remain far from equilibrium even in their steady state by Langevin equations with colored noise which is chosen independently from the friction contribution. Since these Langevin equations are typically not derived from first-principle Hamiltonian dynamics, it is not clear whether they correspond to physically realizable scenarios. By exact Mori projection in phase space we derive the nonequilibrium generalized Langevin equation (GLE) for an arbitrary phase-space dependent observable A from a generic many-body Hamiltonian with a time-dependent external force h(t) acting on the same observable A. This is the same Hamiltonian from which the standard fluctuation-dissipation theorem is derived, which reflects the generality of our approach. The observable A could, for example, be the position of an atom, of a molecule or of a macroscopic object, the distance between two such entities or a more complex phase-space function such as the reaction coordinate of a chemical reaction or of the folding of a protein. The Hamiltonian could, for example, describe a fluid, a solid, a viscoelastic medium, or even a turbulent inhomogeneous environment. The GLE, which is a closed-form equation of motion for the observable A, is obtained in explicit form to all orders in h(t) and without restrictions on the type of many-body Hamiltonian or the observable A. If the dynamics of the observable A corresponds to a Gaussian process, the resultant GLE has a similar form as the equilibrium Mori GLE, and in particular the friction memory kernel is given by the two-point autocorrelation function of the sum of the complementary and the external force h(t). This is a nontrivial and useful result, as many observables that characterize nonequilibrium systems display Gaussian statistics. For non-Gaussian nonequilibrium observables correction terms appear in the GLE and in the relation between the force autocorrelation and the friction memory kernel, which are explicitly given in terms of cubic correlation functions of A. Interpreting the external force h(t) as a stochastic process, we derive nonequilibrium corrections to the fluctuation-dissipation theorem and present methods to extract all GLE parameters from experimental or simulation time-series data, thus making our nonequilibrium GLE a practical tool to study and model general nonequilibrium systems.
Collapse
|
3
|
Klimek A, Mondal D, Block S, Sharma P, Netz RR. Data-driven classification of individual cells by their non-Markovian motion. Biophys J 2024; 123:1173-1183. [PMID: 38515300 PMCID: PMC11140416 DOI: 10.1016/j.bpj.2024.03.023] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2024] [Revised: 03/11/2024] [Accepted: 03/18/2024] [Indexed: 03/23/2024] Open
Abstract
We present a method to differentiate organisms solely by their motion based on the generalized Langevin equation (GLE) and use it to distinguish two different swimming modes of strongly confined unicellular microalgae Chlamydomonas reinhardtii. The GLE is a general model for active or passive motion of organisms and particles that can be derived from a time-dependent general many-body Hamiltonian and in particular includes non-Markovian effects (i.e., the trajectory memory of its past). We extract all GLE parameters from individual cell trajectories and perform an unbiased cluster analysis to group them into different classes. For the specific cell population employed in the experiments, the GLE-based assignment into the two different swimming modes works perfectly, as checked by control experiments. The classification and sorting of single cells and organisms is important in different areas; our method, which is based on motion trajectories, offers wide-ranging applications in biology and medicine.
Collapse
Affiliation(s)
- Anton Klimek
- Fachbereich Physik, Freie Universität Berlin, Berlin, Germany
| | - Debasmita Mondal
- Department of Physics, Indian Institute of Science, Bangalore, India; James Franck Institute, University of Chicago, Chicago, Illinois
| | - Stephan Block
- Institut für Chemie und Biochemie, Freie Universität Berlin, Berlin, Germany
| | - Prerna Sharma
- Department of Physics, Indian Institute of Science, Bangalore, India; Department of Bioengineering, Indian Institute of Science, Bangalore, India
| | - Roland R Netz
- Fachbereich Physik, Freie Universität Berlin, Berlin, Germany.
| |
Collapse
|
4
|
Klages R. Cell migration: Beyond Brownian motion. Biophys J 2024; 123:1167-1169. [PMID: 38637988 PMCID: PMC11140459 DOI: 10.1016/j.bpj.2024.04.013] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2024] [Revised: 04/07/2024] [Accepted: 04/12/2024] [Indexed: 04/20/2024] Open
Affiliation(s)
- Rainer Klages
- Centre for Complex Systems, School of Mathematical Sciences, Queen Mary University of London, London, United Kingdom; London Mathematical Laboratory, London, United Kingdom.
| |
Collapse
|
5
|
Dalton BA, Kiefer H, Netz RR. The role of memory-dependent friction and solvent viscosity in isomerization kinetics in viscogenic media. Nat Commun 2024; 15:3761. [PMID: 38704367 PMCID: PMC11069540 DOI: 10.1038/s41467-024-48016-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2023] [Accepted: 04/18/2024] [Indexed: 05/06/2024] Open
Abstract
Molecular isomerization kinetics in liquid solvent depends on a complex interplay between the solvent friction acting on the molecule, internal dissipation effects (also known as internal friction), the viscosity of the solvent, and the dihedral free energy profile. Due to the absence of accurate techniques to directly evaluate isomerization friction, it has not been possible to explore these relationships in full. By combining extensive molecular dynamics simulations with friction memory-kernel extraction techniques we consider a variety of small, isomerising molecules under a range of different viscogenic conditions and directly evaluate the viscosity dependence of the friction acting on a rotating dihedral. We reveal that the influence of different viscogenic media on isomerization kinetics can be dramatically different, even when measured at the same viscosity. This is due to the dynamic solute-solvent coupling, mediated by time-dependent friction memory kernels. We also show that deviations from the linear dependence of isomerization rates on solvent viscosity, which are often simply attributed to internal friction effects, are due to the simultaneous violation of two fundamental relationships: the Stokes-Einstein relation and the overdamped Kramers prediction for the barrier-crossing rate, both of which require explicit knowledge of friction.
Collapse
Affiliation(s)
| | - Henrik Kiefer
- Freie Universität Berlin, Fachbereich Physik, Berlin, Germany
| | - Roland R Netz
- Freie Universität Berlin, Fachbereich Physik, Berlin, Germany.
| |
Collapse
|
6
|
Tepper L, Dalton B, Netz RR. Accurate Memory Kernel Extraction from Discretized Time-Series Data. J Chem Theory Comput 2024; 20:3061-3068. [PMID: 38603471 PMCID: PMC11044577 DOI: 10.1021/acs.jctc.3c01289] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2023] [Revised: 03/25/2024] [Accepted: 04/01/2024] [Indexed: 04/13/2024]
Abstract
Memory effects emerge as a fundamental consequence of dimensionality reduction when low-dimensional observables are used to describe the dynamics of complex many-body systems. In the context of molecular dynamics (MD) data analysis, accounting for memory effects using the framework of the generalized Langevin equation (GLE) has proven efficient, accurate, and insightful, particularly when working with high-resolution time series data. However, in experimental systems, high-resolution data are often unavailable, raising questions about the impact of the data resolution on the estimated GLE parameters. This study demonstrates that direct memory extraction from time series data remains accurate when the discretization time is below the memory time. To obtain memory functions reliably, even when the discretization time exceeds the memory time, we introduce a Gaussian Process Optimization (GPO) scheme. This scheme minimizes the deviation of discretized two-point correlation functions between time series data and GLE simulations and is able to estimate accurate memory kernels as long as the discretization time stays below the longest time scale in the data, typically the barrier crossing time.
Collapse
Affiliation(s)
- Lucas Tepper
- Department of Physics, Freie
Universität Berlin, 14195 Berlin, Germany
| | - Benjamin Dalton
- Department of Physics, Freie
Universität Berlin, 14195 Berlin, Germany
| | - Roland R. Netz
- Department of Physics, Freie
Universität Berlin, 14195 Berlin, Germany
| |
Collapse
|
7
|
Heyn JCJ, Rädler JO, Falcke M. Mesenchymal cell migration on one-dimensional micropatterns. Front Cell Dev Biol 2024; 12:1352279. [PMID: 38694822 PMCID: PMC11062138 DOI: 10.3389/fcell.2024.1352279] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2023] [Accepted: 03/29/2024] [Indexed: 05/04/2024] Open
Abstract
Quantitative studies of mesenchymal cell motion are important to elucidate cytoskeleton function and mechanisms of cell migration. To this end, confinement of cell motion to one dimension (1D) significantly simplifies the problem of cell shape in experimental and theoretical investigations. Here we review 1D migration assays employing micro-fabricated lanes and reflect on the advantages of such platforms. Data are analyzed using biophysical models of cell migration that reproduce the rich scenario of morphodynamic behavior found in 1D. We describe basic model assumptions and model behavior. It appears that mechanical models explain the occurrence of universal relations conserved across different cell lines such as the adhesion-velocity relation and the universal correlation between speed and persistence (UCSP). We highlight the unique opportunity of reproducible and standardized 1D assays to validate theory based on statistical measures from large data of trajectories and discuss the potential of experimental settings embedding controlled perturbations to probe response in migratory behavior.
Collapse
Affiliation(s)
- Johannes C. J. Heyn
- Fakultät für Physik, Ludwig-Maximilians-Universität München (LMU), Munich, Germany
| | - Joachim O. Rädler
- Fakultät für Physik, Ludwig-Maximilians-Universität München (LMU), Munich, Germany
| | - Martin Falcke
- Max Delbrück Center for Molecular Medicine in the Helmholtz Association, Berlin, Germany
- Department of Physics, Humboldt University, Berlin, Germany
| |
Collapse
|
8
|
Shea J, Jung G, Schmid F. Force renormalization for probes immersed in an active bath. SOFT MATTER 2024; 20:1767-1785. [PMID: 38305056 DOI: 10.1039/d3sm01387a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/03/2024]
Abstract
Langevin equations or generalized Langevin equations (GLEs) are popular models for describing the motion of a particle in a fluid medium in an effective manner. Here we examine particles immersed in an inherently nonequilibrium fluid, i.e., an active bath, which are subject to an external force. Specifically, we consider two types of forces that are highly relevant for microrheological studies: A harmonic, trapping force and a constant, "drag" force. We study such systems by molecular simulations and use the simulation data to extract an effective GLE description. We find that within this description, in an active bath, the external force in the GLE is not equal to the physical external force, but rather a renormalized external force, which can be significantly smaller. The effect cannot be attributed to the mere temperature renormalization, which is also observed.
Collapse
Affiliation(s)
- Jeanine Shea
- Institut für Theoretische Physik, Technische Universität Berlin, Hardenbergstr. 36, 10623 Berlin, Germany.
| | - Gerhard Jung
- Laboratoire Charles Coulomb (L2C), Université de Montpellier, CNRS, 34095 Montpellier, France
| | - Friederike Schmid
- Institut für Physik, Johannes Gutenberg-Universität Mainz, 55099 Mainz, Germany.
| |
Collapse
|
9
|
Zhdanov VP. Kinetics of cancer metastasis. Biosystems 2024; 235:105098. [PMID: 38056592 DOI: 10.1016/j.biosystems.2023.105098] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/03/2023] [Revised: 11/14/2023] [Accepted: 12/02/2023] [Indexed: 12/08/2023]
Abstract
The formation of metastases during cancer is now considered to be induced by migrating metastatic stem cells (MetSCs) in preexisting niches or niches induced by MetSCs or tumor-derived exosomes (TDEs). I propose and compare two simplest generic models describing these two scenarios. The number of tumors is predicted (i) to increase exponentially in the case of preexisting niches and (ii) to diverge during a finite time interval in the case of induced niches. The latter prediction is novel and of interest because rapid collapse in the end of a finite time interval is a well-known feature of the cancer metastasis. Two advanced models describing the two scenarios of cancer metastasis have been scrutinized as well. These models clarify the likely role of various specific factors in the metastasis. In particular, the equations derived in the framework of the advanced model with preexisting niches have been solved analytically allowing (i) to clarify the factors determining the duration of the period from the initiation of the primary tumor to the phase when the metastases start to dominate, (ii) to estimate the number of metastases in the end of this period, and (iii) to explains why the use of chemotherapy typically results in the improvement of the patient state only for a relatively short period. The equations derived in the framework of the advanced model with induced niches have no analytical solution, and their analysis merits additional attention.
Collapse
Affiliation(s)
- Vladimir P Zhdanov
- Boreskov Institute of Catalysis, Russian Academy of Sciences, Novosibirsk, Russia.
| |
Collapse
|
10
|
Sayer T, Montoya-Castillo A. Compact and complete description of non-Markovian dynamics. J Chem Phys 2023; 158:014105. [PMID: 36610963 DOI: 10.1063/5.0132614] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023] Open
Abstract
Generalized master equations provide a theoretically rigorous framework to capture the dynamics of processes ranging from energy harvesting in plants and photovoltaic devices to qubit decoherence in quantum technologies and even protein folding. At their center is the concept of memory. The explicit time-nonlocal description of memory is both protracted and elaborate. When physical intuition is at a premium, one would desire a more compact, yet complete, description. Here, we demonstrate how and when the time-convolutionless formalism constitutes such a description. In particular, by focusing on the dissipative dynamics of the spin-boson and Frenkel exciton models, we show how to: easily construct the time-local generator from reference reduced dynamics, elucidate the dependence of its existence on the system parameters and the choice of reduced observables, identify the physical origin of its apparent divergences, and offer analysis tools to diagnose their severity and circumvent their deleterious effects. We demonstrate that, when applicable, the time-local approach requires as little information as the more commonly used time-nonlocal scheme, with the important advantages of providing a more compact description, greater algorithmic simplicity, and physical interpretability. We conclude by introducing the discrete-time analog and a straightforward protocol to employ it in cases where the reference dynamics have limited resolution. The insights we present here offer the potential for extending the reach of dynamical methods, reducing both their cost and conceptual complexity.
Collapse
Affiliation(s)
- Thomas Sayer
- Department of Chemistry, University of Colorado Boulder, Boulder, Colorado 80309, USA
| | | |
Collapse
|
11
|
Xue C, Huang Y, Zheng X, Hu G. Hopping Behavior Mediates the Anomalous Confined Diffusion of Nanoparticles in Porous Hydrogels. J Phys Chem Lett 2022; 13:10612-10620. [PMID: 36350083 DOI: 10.1021/acs.jpclett.2c02733] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/16/2023]
Abstract
Diffusion is an essential means of mass transport in porous materials such as hydrogels, which are appealing in various biomedical applications. Herein, we investigate the diffusive motion of nanoparticles (NPs) in porous hydrogels to provide a microscopic view of confined diffusion. Based on the mean square displacement from particle tracking experiments, we elucidate the anomalous diffusion dynamics of the embedded NPs and reveal the heterogeneous pore structures in hydrogels. The results demonstrate that diffusive NPs can intermittently escape from single pores through void connective pathways and exhibit non-Gaussian displacement probability distribution. We simulate this scenario using the Monte Carlo method and clarify the existence of hopping events in porous diffusion. The resultant anomalous diffusion can be fully depicted by combining the hopping mechanism and the hydrodynamic effect. Our results highlight the hopping behavior through the connective pathways and establish a hybrid model to predict NP transport in porous environments.
Collapse
Affiliation(s)
- Chundong Xue
- State Key Laboratory of Nonlinear Mechanics, Beijing Key Laboratory of Engineered Construction and Mechanobiology, Institute of Mechanics, Chinese Academy of Sciences, Beijing100190, China
- School of Optoelectronic Engineering and Instrumentation Science, Dalian University of Technology, Dalian116024, China
| | - Yirong Huang
- State Key Laboratory of Nonlinear Mechanics, Beijing Key Laboratory of Engineered Construction and Mechanobiology, Institute of Mechanics, Chinese Academy of Sciences, Beijing100190, China
| | - Xu Zheng
- State Key Laboratory of Nonlinear Mechanics, Beijing Key Laboratory of Engineered Construction and Mechanobiology, Institute of Mechanics, Chinese Academy of Sciences, Beijing100190, China
| | - Guoqing Hu
- Department of Engineering Mechanics, Zhejiang University, Hangzhou310027, China
| |
Collapse
|
12
|
Shaebani MR, Piel M, Lautenschläger F. Distinct speed and direction memories of migrating dendritic cells diversify their search strategies. Biophys J 2022; 121:4099-4108. [PMID: 36181271 PMCID: PMC9675022 DOI: 10.1016/j.bpj.2022.09.033] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/21/2022] [Revised: 08/10/2022] [Accepted: 09/26/2022] [Indexed: 11/19/2022] Open
Abstract
Migrating cells exhibit various motility patterns, resulting from different migration mechanisms, cell properties, or cell-environment interactions. The complexity of cell dynamics is reflected, e.g., in the diversity of the observed forms of velocity autocorrelation function-which has been widely served as a measure of diffusivity and spreading. By analyzing the dynamics of migrating dendritic cells in vitro, we disentangle the contributions of direction θ and speed v to the velocity autocorrelation. We find that the ability of cells to maintain their speed or direction of motion is unequal, reflected in different temporal decays of speed and direction autocorrelation functions, ACv(t)∼t-1.2 and ACθ(t)∼t-0.5, respectively. The larger power-law exponent of ACv(t) indicates that the cells lose their speed memory considerably faster than the direction memory. Using numerical simulations, we investigate the influence of ACθ and ACv as well as the direction-speed cross correlation Cθ-v on the search time of a persistent random walker to find a randomly located target in confinement. Although ACθ and Cθ-v play the major roles, we find that the speed autocorrelation ACv can be also tuned to minimize the search time. Adopting an optimal ACv can reduce the search time even up to 10% compared with uncorrelated spontaneous speeds. Our results suggest that migrating cells can improve their search efficiency, especially in crowded environments, through the directional or speed persistence or the speed-direction correlation.
Collapse
Affiliation(s)
- M Reza Shaebani
- Department of Theoretical Physics, Saarland University, Saarbrücken, Germany; Centre for Biophysics, Saarland University, Saarbrücken, Germany.
| | - Matthieu Piel
- Institut Curie and Institut Pierre Gilles de Gennes, PSL Research University, CNRS, UMR 144, Paris, France
| | - Franziska Lautenschläger
- Centre for Biophysics, Saarland University, Saarbrücken, Germany; Department of Experimental Physics, Saarland University, Saarbrücken, Germany
| |
Collapse
|
13
|
Sinha S, Li X, Das R, Thirumalai D. Mechanical feedback controls the emergence of dynamical memory in growing tissue monolayers. J Chem Phys 2022; 156:245101. [DOI: 10.1063/5.0087815] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The growth of a tissue, which depends on cell–cell interactions and biologically relevant processes such as cell division and apoptosis, is regulated by a mechanical feedback mechanism. We account for these effects in a minimal two-dimensional model in order to investigate the consequences of mechanical feedback, which is controlled by a critical pressure, p c. A cell can only grow and divide if its pressure, due to interaction with its neighbors, is less than p c. Because temperature is not a relevant variable, the cell dynamics is driven by self-generated active forces (SGAFs) that arise due to cell division. We show that even in the absence of intercellular interactions, cells undergo diffusive behavior. The SGAF-driven diffusion is indistinguishable from the well-known dynamics of a free Brownian particle at a fixed finite temperature. When intercellular interactions are taken into account, we find persistent temporal correlations in the force–force autocorrelation function (FAF) that extends over a timescale of several cell division times. The time-dependence of the FAF reveals memory effects, which increases as p c increases. The observed non-Markovian effects emerge due to the interplay of cell division and mechanical feedback and are inherently a non-equilibrium phenomenon.
Collapse
Affiliation(s)
- Sumit Sinha
- Department of Physics, University of Texas at Austin, Austin, Texas 78712, USA
| | - Xin Li
- Department of Chemistry, University of Texas at Austin, Austin, Texas 78712, USA
| | - Rajsekhar Das
- Department of Chemistry, University of Texas at Austin, Austin, Texas 78712, USA
| | - D. Thirumalai
- Department of Chemistry, University of Texas at Austin, Austin, Texas 78712, USA
| |
Collapse
|
14
|
Rickert CA, Lieleg O. Machine learning approaches for biomolecular, biophysical, and biomaterials research. BIOPHYSICS REVIEWS 2022; 3:021306. [PMID: 38505413 PMCID: PMC10914139 DOI: 10.1063/5.0082179] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/13/2021] [Accepted: 05/12/2022] [Indexed: 03/21/2024]
Abstract
A fluent conversation with a virtual assistant, person-tailored news feeds, and deep-fake images created within seconds-all those things that have been unthinkable for a long time are now a part of our everyday lives. What these examples have in common is that they are realized by different means of machine learning (ML), a technology that has fundamentally changed many aspects of the modern world. The possibility to process enormous amount of data in multi-hierarchical, digital constructs has paved the way not only for creating intelligent systems but also for obtaining surprising new insight into many scientific problems. However, in the different areas of biosciences, which typically rely heavily on the collection of time-consuming experimental data, applying ML methods is a bit more challenging: Here, difficulties can arise from small datasets and the inherent, broad variability, and complexity associated with studying biological objects and phenomena. In this Review, we give an overview of commonly used ML algorithms (which are often referred to as "machines") and learning strategies as well as their applications in different bio-disciplines such as molecular biology, drug development, biophysics, and biomaterials science. We highlight how selected research questions from those fields were successfully translated into machine readable formats, discuss typical problems that can arise in this context, and provide an overview of how to resolve those encountered difficulties.
Collapse
|
15
|
Vroylandt H, Goudenège L, Monmarché P, Pietrucci F, Rotenberg B. Likelihood-based non-Markovian models from molecular dynamics. Proc Natl Acad Sci U S A 2022; 119:e2117586119. [PMID: 35320038 PMCID: PMC9060509 DOI: 10.1073/pnas.2117586119] [Citation(s) in RCA: 18] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2021] [Accepted: 02/16/2022] [Indexed: 01/09/2023] Open
Abstract
SignificanceThe analysis of complex systems with many degrees of freedom generally involves the definition of low-dimensional collective variables more amenable to physical understanding. Their dynamics can be modeled by generalized Langevin equations, whose coefficients have to be estimated from simulations of the initial high-dimensional system. These equations feature a memory kernel describing the mutual influence of the low-dimensional variables and their environment. We introduce and implement an approach where the generalized Langevin equation is designed to maximize the statistical likelihood of the observed data. This provides an efficient way to generate reduced models to study dynamical properties of complex processes such as chemical reactions in solution, conformational changes in biomolecules, or phase transitions in condensed matter systems.
Collapse
Affiliation(s)
- Hadrien Vroylandt
- Institut des Sciences du Calcul et des Données, Sorbonne Université, F-75005 Paris, France
| | - Ludovic Goudenège
- CNRS, FR 3487, Fédération de Mathématiques de CentraleSupélec, CentraleSupélec, Université Paris-Saclay, 91190 Gif-sur-Yvette, France
| | - Pierre Monmarché
- Laboratoire Jacques-Louis Lions, Sorbonne Université, F-75005 Paris, France
- Laboratoire de Chimie Théorique, Sorbonne Université, F-75005 Paris, France
| | - Fabio Pietrucci
- Muséum National d’Histoire Naturelle, UMR CNRS 7590, Institut de Minéralogie, de Physique des Matériaux et de Cosmochimie, Sorbonne Université, F-75005 Paris, France
| | - Benjamin Rotenberg
- Physicochimie des Électrolytes et Nanosystèmes Interfaciaux, Sorbonne Université, CNRS, F-75005 Paris, France
| |
Collapse
|
16
|
Jung G. Non-Markovian systems out of equilibrium: exact results for two routes of coarse graining. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2022; 34:204004. [PMID: 35180708 DOI: 10.1088/1361-648x/ac56a7] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/25/2021] [Accepted: 02/18/2022] [Indexed: 06/14/2023]
Abstract
Generalized Langevin equations (GLEs) can be systematically derived via dimensional reduction from high-dimensional microscopic systems. For linear models the derivation can either be based on projection operator techniques such as the Mori-Zwanzig (MZ) formalism or by 'integrating out' the bath degrees of freedom. Based on exact analytical results we show that both routes can lead to fundamentally different GLEs and that the origin of these differences is based inherently on the non-equilibrium nature of the microscopic stochastic model. The most important conceptional difference between the two routes is that the MZ result intrinsically fulfills the generalized second fluctuation-dissipation theorem while the integration result can lead to its violation. We supplement our theoretical findings with numerical and simulation results for two popular non-equilibrium systems: time-delayed feedback control and the active Ornstein-Uhlenbeck process.
Collapse
Affiliation(s)
- Gerhard Jung
- Department of Chemical Engineering, Kyoto University, Japan
- Laboratoire Charles Coulomb (L2C), Université de Montpellier, CNRS, 34095 Montpellier, France
| |
Collapse
|
17
|
Disentangling cadherin-mediated cell-cell interactions in collective cancer cell migration. Biophys J 2022; 121:44-60. [PMID: 34890578 PMCID: PMC8758422 DOI: 10.1016/j.bpj.2021.12.006] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2021] [Revised: 10/30/2021] [Accepted: 12/06/2021] [Indexed: 01/07/2023] Open
Abstract
Cell dispersion from a confined area is fundamental in a number of biological processes, including cancer metastasis. To date, a quantitative understanding of the interplay of single-cell motility, cell proliferation, and intercellular contacts remains elusive. In particular, the role of E- and N-cadherin junctions, central components of intercellular contacts, is still controversial. Combining theoretical modeling with in vitro observations, we investigate the collective spreading behavior of colonies of human cancer cells (T24). The spreading of these colonies is driven by stochastic single-cell migration with frequent transient cell-cell contacts. We find that inhibition of E- and N-cadherin junctions decreases colony spreading and average spreading velocities, without affecting the strength of correlations in spreading velocities of neighboring cells. Based on a biophysical simulation model for cell migration, we show that the behavioral changes upon disruption of these junctions can be explained by reduced repulsive excluded volume interactions between cells. This suggests that in cancer cell migration, cadherin-based intercellular contacts sharpen cell boundaries leading to repulsive rather than cohesive interactions between cells, thereby promoting efficient cell spreading during collective migration.
Collapse
|
18
|
Abstract
Protein-folding kinetics is often described as Markovian (i.e., memoryless) diffusion in a one-dimensional free energy landscape, governed by an instantaneous friction coefficient that is fitted to reproduce experimental or simulated folding times. For the α-helix forming polypeptide alanine9 and a specific reaction coordinate that consists of the summed native hydrogen-bond lengths, we demonstrate that the friction extracted from molecular dynamics simulations exhibits significant memory with a decay time that is in the nanosecond range and thus, of the same order as the folding and unfolding times. Our non-Markovian modeling not only reproduces the molecular dynamics simulations accurately but also demonstrates that memory friction effects lead to anomalous and drastically accelerated protein kinetics. 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 alanine9 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.
Collapse
|
19
|
Abstract
Creating artificial macromolecular transport systems that can support the movement of molecules along defined routes is a key goal of nanotechnology. Here, we report the bottom-up construction of a macromolecular transport system in which molecular pistons diffusively move through micrometer-long, hollow filaments. The pistons can cover micrometer distances in fractions of seconds. We build the system using multi-layer DNA origami and analyze the structures of the components using transmission electron microscopy. We study the motion of the pistons along the tubes using single-molecule fluorescence microscopy and perform Langevin simulations to reveal details of the free energy surface that directs the motions of the pistons. The tubular transport system achieves diffusivities and displacement ranges known from natural molecular motors and realizes mobility improvements over five orders of magnitude compared to previous artificial random walker designs. Electric fields can also be employed to actively pull the pistons along the filaments, thereby realizing a nanoscale electric rail system. Our system presents a platform for artificial motors that move autonomously driven by chemical fuels and for performing nanotribology studies, and it could form a basis for future molecular transportation networks. DNA origami can be used to control the movement of nanoscale assemblies. Here the authors construct multiple-micrometer-long hollow DNA filaments through which DNA pistons move with micrometer-per-second speeds.
Collapse
|
20
|
Jung G, Schmid F. Fluctuation-dissipation relations far from equilibrium: a case study. SOFT MATTER 2021; 17:6413-6425. [PMID: 34132298 PMCID: PMC8262459 DOI: 10.1039/d1sm00521a] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 04/07/2021] [Accepted: 06/06/2021] [Indexed: 06/12/2023]
Abstract
Fluctuation-dissipation relations or "theorems" (FDTs) are fundamental for statistical physics and can be rigorously derived for equilibrium systems. Their applicability to non-equilibrium systems is, however, debated. Here, we simulate an active microrheology experiment, in which a spherical colloid is pulled with a constant external force through a fluid, creating near-equilibrium and far-from-equilibrium systems. We characterize the structural and dynamical properties of these systems, and reconstruct an effective generalized Langevin equation (GLE) for the colloid dynamics. Specifically, we test the validity of two FDTs: The first FDT relates the non-equilibrium response of a system to equilibrium correlation functions, and the second FDT relates the memory friction kernel in the GLE to the stochastic force. We find that the validity of the first FDT depends strongly on the strength of the external driving: it is fulfilled close to equilibrium and breaks down far from it. In contrast, we observe that the second FDT is always fulfilled. We provide a mathematical argument why this generally holds for memory kernels reconstructed from a deterministic Volterra equation for correlation functions, even for non-stationary non-equilibrium systems. Motivated by the Mori-Zwanzig formalism, we therefore suggest to impose an orthogonality constraint on the stochastic force, which is in fact equivalent to the validity of this Volterra equation. Such GLEs automatically satisfy the second FDT and are unique, which is desirable when using GLEs for coarse-grained modeling.
Collapse
Affiliation(s)
- Gerhard Jung
- Institut für Theoretische Physik, Universität Innsbruck, Technikerstraße 21A, A-6020 Innsbruck, Austria.
| | - Friederike Schmid
- Institut für Physik, Johannes Gutenberg-Universität Mainz, 55099 Mainz, Germany.
| |
Collapse
|
21
|
d’Alessandro J, Barbier--Chebbah A, Cellerin V, Benichou O, Mège RM, Voituriez R, Ladoux B. Cell migration guided by long-lived spatial memory. Nat Commun 2021; 12:4118. [PMID: 34226542 PMCID: PMC8257581 DOI: 10.1038/s41467-021-24249-8] [Citation(s) in RCA: 36] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2021] [Accepted: 06/08/2021] [Indexed: 02/06/2023] Open
Abstract
Living cells actively migrate in their environment to perform key biological functions-from unicellular organisms looking for food to single cells such as fibroblasts, leukocytes or cancer cells that can shape, patrol or invade tissues. Cell migration results from complex intracellular processes that enable cell self-propulsion, and has been shown to also integrate various chemical or physical extracellular signals. While it is established that cells can modify their environment by depositing biochemical signals or mechanically remodelling the extracellular matrix, the impact of such self-induced environmental perturbations on cell trajectories at various scales remains unexplored. Here, we show that cells can retrieve their path: by confining motile cells on 1D and 2D micropatterned surfaces, we demonstrate that they leave long-lived physicochemical footprints along their way, which determine their future path. On this basis, we argue that cell trajectories belong to the general class of self-interacting random walks, and show that self-interactions can rule large scale exploration by inducing long-lived ageing, subdiffusion and anomalous first-passage statistics. Altogether, our joint experimental and theoretical approach points to a generic coupling between motile cells and their environment, which endows cells with a spatial memory of their path and can dramatically change their space exploration.
Collapse
Affiliation(s)
- Joseph d’Alessandro
- grid.508487.60000 0004 7885 7602Université de Paris, CNRS, Institut Jacques Monod, Paris, F-75006 France
| | - Alex Barbier--Chebbah
- grid.462844.80000 0001 2308 1657Laboratoire de Physique Théorique de la Matière Condensée, CNRS/Sorbonne Université, Paris, France
| | - Victor Cellerin
- grid.508487.60000 0004 7885 7602Université de Paris, CNRS, Institut Jacques Monod, Paris, F-75006 France
| | - Olivier Benichou
- grid.462844.80000 0001 2308 1657Laboratoire de Physique Théorique de la Matière Condensée, CNRS/Sorbonne Université, Paris, France
| | - René Marc Mège
- grid.508487.60000 0004 7885 7602Université de Paris, CNRS, Institut Jacques Monod, Paris, F-75006 France
| | - Raphaël Voituriez
- grid.462844.80000 0001 2308 1657Laboratoire Jean Perrin and Laboratoire de Physique Théorique de la Matière Condensée, CNRS/Sorbonne Université, Paris, France
| | - Benoît Ladoux
- grid.508487.60000 0004 7885 7602Université de Paris, CNRS, Institut Jacques Monod, Paris, F-75006 France
| |
Collapse
|
22
|
Klippenstein V, Tripathy M, Jung G, Schmid F, van der Vegt NFA. Introducing Memory in Coarse-Grained Molecular Simulations. J Phys Chem B 2021; 125:4931-4954. [PMID: 33982567 PMCID: PMC8154603 DOI: 10.1021/acs.jpcb.1c01120] [Citation(s) in RCA: 40] [Impact Index Per Article: 13.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Abstract
Preserving the correct dynamics at the coarse-grained (CG) level is a pressing problem in the development of systematic CG models in soft matter simulation. Starting from the seminal idea of simple time-scale mapping, there have been many efforts over the years toward establishing a meticulous connection between the CG and fine-grained (FG) dynamics based on fundamental statistical mechanics approaches. One of the most successful attempts in this context has been the development of CG models based on the Mori-Zwanzig (MZ) theory, where the resulting equation of motion has the form of a generalized Langevin equation (GLE) and closely preserves the underlying FG dynamics. In this Review, we describe some of the recent studies in this regard. We focus on the construction and simulation of dynamically consistent systematic CG models based on the GLE, both in the simple Markovian limit and the non-Markovian case. Some recent studies of physical effects of memory are also discussed. The Review is aimed at summarizing recent developments in the field while highlighting the major challenges and possible future directions.
Collapse
Affiliation(s)
- Viktor Klippenstein
- Eduard-Zintl-Institut für Anorganische und Physikalische Chemie, Technische Universität Darmstadt, 64287 Darmstadt, Germany
| | - Madhusmita Tripathy
- Eduard-Zintl-Institut für Anorganische und Physikalische Chemie, Technische Universität Darmstadt, 64287 Darmstadt, Germany
| | - Gerhard Jung
- Institut für Theoretische Physik, Universität Innsbruck, Technikerstraße 21 A, A-6020 Innsbruck, Austria
| | - Friederike Schmid
- Institut für Physik, Johannes Gutenberg-Universität Mainz, Staudingerweg 9, 55128 Mainz, Germany
| | - Nico F A van der Vegt
- Eduard-Zintl-Institut für Anorganische und Physikalische Chemie, Technische Universität Darmstadt, 64287 Darmstadt, Germany
| |
Collapse
|
23
|
Transient surface hydration impacts biogeography and intercellular interactions of non-motile bacteria. Appl Environ Microbiol 2021; 87:AEM.03067-20. [PMID: 33579687 PMCID: PMC8091113 DOI: 10.1128/aem.03067-20] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022] Open
Abstract
There are many hydrated surface niches that are neither static nor continuously flowing that are colonized by microbes such as bacteria. Such periodic hydrodynamic regimes are distinct from aquatic systems where microbial dissemination is reasonably predicted by assuming continuous flow or static systems where motile microbes largely control their own fate. Here we show how non-motile bacteria exhibit rapid, dispersive bursts of movement over surfaces using transient confluent hydration from the environment, which we term "surface hydrodispersion" where cells traverse thousands of cell lengths within minutes. The fraction of the population disseminated by surface hydrodispersion is small-on order of 1 cell per million. Thus, surface hydrodispersion can promote isolated distribution of single cells, which is unlike other characterized active and passive surface motilities. We describe this translocation using a continuous time random walk modeling approach and find in computational simulations that transient fluid accumulation, dilution, and gravitational pull are the contributing factors. Surface hydrodispersion, consistent with advection, is unlike simple colony expansion as it dramatically alters spatial relationships, shown here with Staphylococcus aureus, which becomes increasingly virulent when isolated from Corynebacterium striatum Surface hydrodispersion of non-motile bacteria exploiting transient fluid availability and gravity is a mechanism that can result in sporadic and sudden shifts in microbial community behavior. To better understand how this movement can impact biogeography on the millimeter scale, this work describes a system for study of primary factors behind this movement as well as a stochastic model describing this dispersal.Importance: Understanding the dynamics within microbiome communities is a challenge. Knowledge of phylogeny and spatial arrangement has led to increased understanding of numerous polymicrobial communities yet, these snapshots do not convey the dynamics of populations over time. The actual biogeography of any microbiome controls the potential interactions, governing any possible antagonistic or synergistic behavior. Accordingly, a shift in biogeography can enable new behavior. Little is known about the movement mechanisms of "non-motile" microbes. Here we characterize a universal means of movement we term hydrodispersion where non-motile bacteria are transported thousands of cell lengths in minutes. We show that only a small fraction of the population is translocated by hydrodispersion and describe this movement further using a random-walk mathematical model approach in silico We demonstrate the importance of hydrodispersion by showing that Staphylococcus aureus can separate from a coculture inoculation with Corynebacterium striatum thus permitting transition to a more virulent state.
Collapse
|
24
|
Mitterwallner BG, Lavacchi L, Netz RR. Negative friction memory induces persistent motion. THE EUROPEAN PHYSICAL JOURNAL. E, SOFT MATTER 2020; 43:67. [PMID: 33099707 DOI: 10.1140/epje/i2020-11992-5] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/25/2020] [Accepted: 09/30/2020] [Indexed: 06/11/2023]
Abstract
We investigate the mean-square displacement (MSD) for random motion governed by the generalized Langevin equation for memory functions that contain two different time scales: In the first model, the memory kernel consists of a delta peak and a single-exponential and in the second model of the sum of two exponentials. In particular, we investigate the scenario where the long-time exponential kernel contribution is negative. The competition between positive and negative friction memory contributions produces an enhanced transient persistent regime in the MSD, which is relevant for biological motility and active matter systems.
Collapse
Affiliation(s)
| | - Laura Lavacchi
- Fachbereich Physik, Freie Universität Berlin, 14195, Berlin, Germany
| | - Roland R Netz
- Fachbereich Physik, Freie Universität Berlin, 14195, Berlin, Germany
| |
Collapse
|