1
|
Domingues TS, Coifman R, Haji-Akbari A. Estimating Position-Dependent and Anisotropic Diffusivity Tensors from Molecular Dynamics Trajectories: Existing Methods and Future Outlook. J Chem Theory Comput 2024; 20:4427-4455. [PMID: 38815171 DOI: 10.1021/acs.jctc.4c00148] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/01/2024]
Abstract
Confinement can substantially alter the physicochemical properties of materials by breaking translational isotropy and rendering all physical properties position-dependent. Molecular dynamics (MD) simulations have proven instrumental in characterizing such spatial heterogeneities and probing the impact of confinement on materials' properties. For static properties, this is a straightforward task and can be achieved via simple spatial binning. Such an approach, however, cannot be readily applied to transport coefficients due to lack of natural extensions of autocorrelations used for their calculation in the bulk. The prime example of this challenge is diffusivity, which, in the bulk, can be readily estimated from the particles' mobility statistics, which satisfy the Fokker-Planck equation. Under confinement, however, such statistics will follow the Smoluchowski equation, which lacks a closed-form analytical solution. This brief review explores the rich history of estimating profiles of the diffusivity tensor from MD simulations and discusses various approximate methods and algorithms developed for this purpose. Besides discussing heuristic extensions of bulk methods, we overview more rigorous algorithms, including kernel-based methods, Bayesian approaches, and operator discretization techniques. Additionally, we outline methods based on applying biasing potentials or imposing constraints on tracer particles. Finally, we discuss approaches that estimate diffusivity from mean first passage time or committor probability profiles, a conceptual framework originally developed in the context of collective variable spaces describing rare events in computational chemistry and biology. In summary, this paper offers a concise survey of diverse approaches for estimating diffusivity from MD trajectories, highlighting challenges and opportunities in this area.
Collapse
Affiliation(s)
- Tiago S Domingues
- Department of Chemical and Environmental Engineering, Yale University, New Haven, Connecticut 06520, United States
| | - Ronald Coifman
- Department of Mathematics, Yale University, New Haven, Connecticut 06520, United States
- Department of Computer Science, Yale University, New Haven, Connecticut 06520, United States
| | - Amir Haji-Akbari
- Department of Chemical and Environmental Engineering, Yale University, New Haven, Connecticut 06520, United States
| |
Collapse
|
2
|
Domingues TS, Coifman RR, Haji-Akbari A. Robust Estimation of Position-Dependent Anisotropic Diffusivity Tensors from Molecular Dynamics Trajectories. J Phys Chem B 2023; 127:8644-8659. [PMID: 37757480 DOI: 10.1021/acs.jpcb.3c03581] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/29/2023]
Abstract
Confinement breaks translational and rotational symmetry in materials and makes all physical properties functions of position. Such spatial variations are key to modulating material properties at the nanoscale, and characterizing them accurately is therefore an intense area of research in the molecular simulations community. This is relatively easy to accomplish for basic mechanical observables. Determining spatial profiles of transport properties, such as diffusivity, is, however, much more challenging, as it requires calculating position-dependent autocorrelations of mechanical observables. In our previous paper (Domingues, T.S.; Coifman, R.; Haji-Akbari, A. J. Phys. Chem. B 2023, 127, 5273 10.1021/acs.jpcb.3c00670), we analytically derive and numerically validate a set of filtered covariance estimators (FCEs) for quantifying spatial variations of the diffusivity tensor from stochastic trajectories. In this work, we adapt these estimators to extract diffusivity profiles from MD trajectories and validate them by applying them to a Lennard-Jones fluid within a slit pore. We find our MD-adapted estimator to exhibit the same qualitative features as its stochastic counterpart, as it accurately estimates the lateral diffusivity across the pore while systematically underestimating the normal diffusivity close to hard boundaries. We introduce a conceptually simple and numerically efficient correction scheme based on simulated annealing and diffusion maps to resolve the latter artifact and obtain normal diffusivity profiles that are consistent with the self-part of the van Hove correlation functions. Our findings demonstrate the potential of this MD-adapted estimator in accurately characterizing spatial variations of diffusivity in confined materials.
Collapse
Affiliation(s)
- Tiago S Domingues
- Department of Chemical and Environmental Engineering, Yale University, New Haven, Connecticut 06520, United States
| | - Ronald R Coifman
- Department of Mathematics, Yale University, New Haven, Connecticut 06520, United States
| | - Amir Haji-Akbari
- Department of Chemical and Environmental Engineering, Yale University, New Haven, Connecticut 06520, United States
| |
Collapse
|
3
|
Domingues TS, Coifman RR, Haji-Akbari A. Robust Estimation of Position-Dependent Anisotropic Diffusivity Tensors from Stochastic Trajectories. J Phys Chem B 2023. [PMID: 37261948 DOI: 10.1021/acs.jpcb.3c00670] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/03/2023]
Abstract
Materials under confinement can possess properties that deviate considerably from their bulk counterparts. Indeed, confinement makes all physical properties position-dependent and possibly anisotropic, and characterizing such spatial variations and directionality has been an intense area of focus in experimental and computational studies of confined matter. While this task is fairly straightforward for simple mechanical observables, it is far more daunting for transport properties such as diffusivity that can only be estimated from autocorrelations of mechanical observables. For instance, there are well established methods for estimating diffusivity from experimentally observed or computationally generated trajectories in bulk systems. No rigorous generalizations of such methods, however, exist for confined systems. In this work, we present two filtered covariance estimators for computing anisotropic and position-dependent diffusivity tensors and validate them by applying them to stochastic trajectories generated according to known diffusivity profiles. These estimators can accurately capture spatial variations that span over several orders of magnitude and that assume different functional forms. Our kernel-based approach is also very robust to implementation details such as the localization function and time discretization and performs significantly better than estimators that are solely based on local covariance. Moreover, the kernel function does not have to be localized and can instead belong to a dictionary of orthogonal functions. Therefore, the proposed estimator can be readily used to obtain functional estimates of diffusivity rather than a tabulated collection of pointwise estimates. Nonetheless, the susceptibility of the proposed estimators to time discretization is higher at the immediate vicinity of hard boundaries. We demonstrate this heightened susceptibility to be common among all covariance-based estimators.
Collapse
Affiliation(s)
- Tiago S Domingues
- Department of Chemical and Environmental Engineering, Yale University, New Haven, Connecticut 06520, United States
| | - Ronald R Coifman
- Department of Mathematics, Yale University, New Haven, Connecticut 06520, United States
| | - Amir Haji-Akbari
- Department of Chemical and Environmental Engineering, Yale University, New Haven, Connecticut 06520, United States
| |
Collapse
|
4
|
Xist nucleates local protein gradients to propagate silencing across the X chromosome. Cell 2021; 184:6174-6192.e32. [PMID: 34813726 PMCID: PMC8671326 DOI: 10.1016/j.cell.2021.10.022] [Citation(s) in RCA: 56] [Impact Index Per Article: 18.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2020] [Revised: 07/29/2021] [Accepted: 10/11/2021] [Indexed: 02/08/2023]
Abstract
The lncRNA Xist forms ∼50 diffraction-limited foci to transcriptionally silence one X chromosome. How this small number of RNA foci and interacting proteins regulate a much larger number of X-linked genes is unknown. We show that Xist foci are locally confined, contain ∼2 RNA molecules, and nucleate supramolecular complexes (SMACs) that include many copies of the critical silencing protein SPEN. Aggregation and exchange of SMAC proteins generate local protein gradients that regulate broad, proximal chromatin regions. Partitioning of numerous SPEN molecules into SMACs is mediated by their intrinsically disordered regions and essential for transcriptional repression. Polycomb deposition via SMACs induces chromatin compaction and the increase in SMACs density around genes, which propagates silencing across the X chromosome. Our findings introduce a mechanism for functional nuclear compartmentalization whereby crowding of transcriptional and architectural regulators enables the silencing of many target genes by few RNA molecules.
Collapse
|
5
|
Parreira P, Martins MCL. The biophysics of bacterial infections: Adhesion events in the light of force spectroscopy. Cell Surf 2021; 7:100048. [PMID: 33665520 PMCID: PMC7898176 DOI: 10.1016/j.tcsw.2021.100048] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2018] [Revised: 08/10/2020] [Accepted: 12/03/2020] [Indexed: 02/08/2023] Open
Abstract
Bacterial infections are the most eminent public health challenge of the 21st century. The primary step leading to infection is bacterial adhesion to the surface of host cells or medical devices, which is mediated by a multitude of molecular interactions. At the interface of life sciences and physics, last years advances in atomic force microscopy (AFM)-based force spectroscopy techniques have made possible to measure the forces driving bacteria-cell and bacteria-materials interactions on a single molecule/cell basis (single molecule/cell force spectroscopy). Among the bacteria-(bio)materials surface interactions, the life-threatening infections associated to medical devices involving Staphylococcus aureus and Escherichia coli are the most eminent. On the other hand, Pseudomonas aeruginosa binding to the pulmonary and urinary tract or the Helicobacter pylori binding to the gastric mucosa, are classical examples of bacteria-host cell interactions that end in serious infections. As we approach the end of the antibiotic era, acquisition of a deeper knowledge of the fundamental forces involved in bacteria - host cells/(bio)materials surface adhesion is crucial for the identification of new ligand-binding events and its assessment as novel targets for alternative anti-infective therapies. This article aims to highlight the potential of AFM-based force spectroscopy for new targeted therapies development against bacterial infections in which adhesion plays a pivotal role and does not aim to be an extensive overview on the AFM technical capabilities and theory of single molecule force spectroscopy.
Collapse
Affiliation(s)
- Paula Parreira
- INEB – Instituto de Engenharia Biomédica, Universidade do Porto, Portugal
- i3S – Instituto de Investigação e Inovação em Saúde, Universidade do Porto, Portugal
| | - M. Cristina L. Martins
- INEB – Instituto de Engenharia Biomédica, Universidade do Porto, Portugal
- i3S – Instituto de Investigação e Inovação em Saúde, Universidade do Porto, Portugal
- ICBAS – Instituto de Ciências Biomédicas Abel Salazar, Universidade do Porto, Portugal
| |
Collapse
|
6
|
Serov AS, Laurent F, Floderer C, Perronet K, Favard C, Muriaux D, Westbrook N, Vestergaard CL, Masson JB. Statistical Tests for Force Inference in Heterogeneous Environments. Sci Rep 2020; 10:3783. [PMID: 32123194 PMCID: PMC7052274 DOI: 10.1038/s41598-020-60220-1] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/20/2019] [Accepted: 02/05/2020] [Indexed: 01/22/2023] Open
Abstract
We devise a method to detect and estimate forces in a heterogeneous environment based on experimentally recorded stochastic trajectories. In particular, we focus on systems modeled by the heterogeneous overdamped Langevin equation. Here, the observed drift includes a "spurious” force term when the diffusivity varies in space. We show how Bayesian inference can be leveraged to reliably infer forces by taking into account such spurious forces of unknown amplitude as well as experimental sources of error. The method is based on marginalizing the force posterior over all possible spurious force contributions. The approach is combined with a Bayes factor statistical test for the presence of forces. The performance of our method is investigated analytically, numerically and tested on experimental data sets. The main results are obtained in a closed form allowing for direct exploration of their properties and fast computation. The method is incorporated into TRamWAy, an open-source software platform for automated analysis of biomolecule trajectories.
Collapse
Affiliation(s)
- Alexander S Serov
- Decision and Bayesian Computation, USR 3756 (C3BI/DBC) & Neuroscience department CNRS UMR 3751, Institut Pasteur, CNRS, Paris, France.
| | - François Laurent
- Decision and Bayesian Computation, USR 3756 (C3BI/DBC) & Neuroscience department CNRS UMR 3751, Institut Pasteur, CNRS, Paris, France
| | - Charlotte Floderer
- Infectious Disease Research Institute of Montpellier, CNRS UMR 9004, University of Montpellier, Montpellier, France
| | - Karen Perronet
- Laboratoire Charles Fabry, Université Paris-Saclay, Institut d'Optique Graduate School, CNRS UMR8501, 91127, Palaiseau Cedex, France
| | - Cyril Favard
- Infectious Disease Research Institute of Montpellier, CNRS UMR 9004, University of Montpellier, Montpellier, France
| | - Delphine Muriaux
- Infectious Disease Research Institute of Montpellier, CNRS UMR 9004, University of Montpellier, Montpellier, France
| | - Nathalie Westbrook
- Laboratoire Charles Fabry, Université Paris-Saclay, Institut d'Optique Graduate School, CNRS UMR8501, 91127, Palaiseau Cedex, France
| | - Christian L Vestergaard
- Decision and Bayesian Computation, USR 3756 (C3BI/DBC) & Neuroscience department CNRS UMR 3751, Institut Pasteur, CNRS, Paris, France.
| | - Jean-Baptiste Masson
- Decision and Bayesian Computation, USR 3756 (C3BI/DBC) & Neuroscience department CNRS UMR 3751, Institut Pasteur, CNRS, Paris, France.
| |
Collapse
|
7
|
Laurent F, Floderer C, Favard C, Muriaux D, Masson JB, Vestergaard CL. Mapping spatio-temporal dynamics of single biomolecules in living cells. Phys Biol 2019; 17:015003. [PMID: 31765328 DOI: 10.1088/1478-3975/ab5167] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Abstract
We present a Bayesian framework for inferring spatio-temporal maps of diffusivity and potential fields from recorded trajectories of single molecules inside living cells. The framework naturally lets us regularise the high-dimensional inference problem using prior distributions in order to obtain robust results. To overcome the computational complexity of inferring thousands of map parameters from large single particle tracking datasets, we developed a stochastic optimisation method based on local mini-batches and parsimonious gradient calculation. We quantified the gain in convergence speed on numerical simulations, and we demonstrated for the first time temporal regularisation and aligned values of the inferred potential fields across multiple time segments. As a proof-of-concept, we mapped the dynamics of HIV-1 Gag proteins involved in the formation of virus-like particles (VLPs) on the plasma membrane of live T cells at high spatial and temporal resolutions. We focused on transient aggregation events lasting only on tenth of the time required for full VLP formation. The framework and optimisation methods are implemented in the TRamWAy open-source software platform for analysing single biomolecule dynamics.
Collapse
Affiliation(s)
- François Laurent
- Decision and Bayesian Computation, Department of Computational Biology, Department of Neuroscience, CNRS USR 3756, CNRS UMR 3571, Institut Pasteur, 25 rue du Docteur Roux, Paris, 75015, France. Hub de Bioinformatique et Biostatistique - Département Biologie Computationnelle, Institut Pasteur, USR 3756 CNRS, Paris, France
| | | | | | | | | | | |
Collapse
|
8
|
Floderer C, Masson JB, Boilley E, Georgeault S, Merida P, El Beheiry M, Dahan M, Roingeard P, Sibarita JB, Favard C, Muriaux D. Single molecule localisation microscopy reveals how HIV-1 Gag proteins sense membrane virus assembly sites in living host CD4 T cells. Sci Rep 2018; 8:16283. [PMID: 30389967 PMCID: PMC6214999 DOI: 10.1038/s41598-018-34536-y] [Citation(s) in RCA: 29] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/11/2018] [Accepted: 10/12/2018] [Indexed: 11/09/2022] Open
Abstract
Monitoring virus assembly at the nanoscale in host cells remains a major challenge. Human immunodeficiency virus type 1 (HIV-1) components are addressed to the plasma membrane where they assemble to form spherical particles of 100 nm in diameter. Interestingly, HIV-1 Gag protein expression alone is sufficient to produce virus-like particles (VLPs) that resemble the immature virus. Here, we monitored VLP formation at the plasma membrane of host CD4+ T cells using a newly developed workflow allowing the analysis of long duration recordings of single-molecule Gag protein localisation and movement. Comparison of Gag assembling platforms in CD4+ T cells expressing wild type or assembly-defective Gag mutant proteins showed that VLP formation lasts roughly 15 minutes with an assembly time of 5 minutes. Trapping energy maps, built from membrane associated Gag protein movements, showed that one third of the assembling energy is due to direct Gag capsid-capsid interaction while the remaining two thirds require the nucleocapsid-RNA interactions. Finally, we show that the viral RNA genome does not increase the attraction of Gag at the membrane towards the assembling site but rather acts as a spatiotemporal coordinator of the membrane assembly process.
Collapse
Affiliation(s)
- Charlotte Floderer
- Infectious Disease Research Institute of Montpellier (IRIM), UMR9004 CNRS, University of Montpellier, 1919 route de Mende, 34293, Montpellier, France
| | - Jean-Baptiste Masson
- Decision and Bayesian Computation, UMR 3571 CNRS, Pasteur Institute, Paris, France
| | - Elise Boilley
- Infectious Disease Research Institute of Montpellier (IRIM), UMR9004 CNRS, University of Montpellier, 1919 route de Mende, 34293, Montpellier, France
| | - Sonia Georgeault
- INSERM U966 and IBiSA EM Facility, University of Tours, Tours, France
| | - Peggy Merida
- Infectious Disease Research Institute of Montpellier (IRIM), UMR9004 CNRS, University of Montpellier, 1919 route de Mende, 34293, Montpellier, France
| | - Mohamed El Beheiry
- Light and Optical Control of Cellular Organization, Curie Institute, UMR, 168 CNRS, Paris, France
| | - Maxime Dahan
- Light and Optical Control of Cellular Organization, Curie Institute, UMR, 168 CNRS, Paris, France
| | | | - Jean-Baptiste Sibarita
- Interdisciplinary Institute for Neuroscience, UMR 5297 CNRS, University of Bordeaux, Bordeaux, France
| | - Cyril Favard
- Infectious Disease Research Institute of Montpellier (IRIM), UMR9004 CNRS, University of Montpellier, 1919 route de Mende, 34293, Montpellier, France.
| | - Delphine Muriaux
- Infectious Disease Research Institute of Montpellier (IRIM), UMR9004 CNRS, University of Montpellier, 1919 route de Mende, 34293, Montpellier, France.
| |
Collapse
|
9
|
Chang JC, Liu Y, Chou T. Reconstruction of Cell Focal Adhesions using Physical Constraints and Compressive Regularization. Biophys J 2017; 113:2530-2539. [PMID: 29212006 PMCID: PMC5738527 DOI: 10.1016/j.bpj.2017.09.021] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2017] [Revised: 09/13/2017] [Accepted: 09/19/2017] [Indexed: 11/30/2022] Open
Abstract
We develop a method to reconstruct, from measured displacements of an underlying elastic substrate, the spatially dependent forces that cells or tissues impart on it. Given newly available high-resolution images of substrate displacements, it is desirable to be able to reconstruct small-scale, compactly supported focal adhesions that are often localized and exist only within the footprint of a cell. In addition to the standard quadratic data mismatch terms that define least-squares fitting, we motivate a regularization term in the objective function that penalizes vectorial invariants of the reconstructed surface stress while preserving boundaries. We solve this inverse problem by providing a numerical method for setting up a discretized inverse problem that is solvable by standard convex optimization techniques. By minimizing the objective function subject to a number of important physically motivated constraints, we are able to efficiently reconstruct stress fields with localized structure from simulated and experimental substrate displacements. Our method incorporates the exact solution for a stress tensor accurate to first-order finite differences and motivates the use of distance-based cutoffs for data inclusion and problem sparsification.
Collapse
Affiliation(s)
- Joshua C Chang
- Epidemiology and Biostatistics Section, Rehabilitation Medicine, Clinical Center, The National Institutes of Health, Bethesda, Maryland
| | - Yanli Liu
- Department of Mathematics, UCLA, Los Angeles, California
| | - Tom Chou
- Department of Mathematics, UCLA, Los Angeles, California; Department of Biomathematics, UCLA, Los Angeles, California.
| |
Collapse
|
10
|
Arya G. Models for recovering the energy landscape of conformational transitions from single-molecule pulling experiments. MOLECULAR SIMULATION 2016. [DOI: 10.1080/08927022.2015.1123257] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
|
11
|
Fuhrmann A. Advances in the Theory of Single-Molecule Force Spectroscopy: Bond Potentials and Mobilities. Biophys J 2015; 109:851. [PMID: 26331241 PMCID: PMC4564686 DOI: 10.1016/j.bpj.2015.07.029] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/01/2015] [Revised: 06/22/2015] [Accepted: 07/14/2015] [Indexed: 10/23/2022] Open
Affiliation(s)
- Alexander Fuhrmann
- Bioengineering Department, University of California at San Diego, La Jolla, California.
| |
Collapse
|