TrackArt: the user friendly interface for single molecule tracking data analysis and simulation applied to complex diffusion in mica supported lipid bilayers

Membrane lipids drive formation of KRAS4b-RAF1 RBDCRD nanoclusters on the membrane

The oncogene RAS, extensively studied for decades, presents persistent gaps in understanding, hindering the development of effective therapeutic strategies due to a lack of precise details on how RAS initiates MAPK signaling with RAF effector proteins at the plasma membrane. Recent advances in X-ray crystallography, cryo-EM, and super-resolution fluorescence microscopy offer structural and spatial insights, yet the molecular mechanisms involving protein-protein and protein-lipid interactions in RAS-mediated signaling require further characterization. This study utilizes single-molecule experimental techniques, nuclear magnetic resonance spectroscopy, and the computational Machine-Learned Modeling Infrastructure (MuMMI) to examine KRAS4b and RAF1 on a biologically relevant lipid bilayer. MuMMI captures long-timescale events while preserving detailed atomic descriptions, providing testable models for experimental validation. Both in vitro and computational studies reveal that RBDCRD binding alters KRAS lateral diffusion on the lipid bilayer, increasing cluster size and decreasing diffusion. RAS and membrane binding cause hydrophobic residues in the CRD region to penetrate the bilayer, stabilizing complexes through β-strand elongation. These cooperative interactions among lipids, KRAS4b, and RAF1 are proposed as essential for forming nanoclusters, potentially a critical step in MAP kinase signal activation.

Maximum likelihood-based estimation of diffusion coefficient is quick and reliable method for analyzing estradiol actions on surface receptor movements

The rapid effects of estradiol on membrane receptors are in the focus of the estradiol research field, however, the molecular mechanisms of these non-classical estradiol actions are poorly understood. Since the lateral diffusion of membrane receptors is an important indicator of their function, a deeper understanding of the underlying mechanisms of non-classical estradiol actions can be achieved by investigating receptor dynamics. Diffusion coefficient is a crucial and widely used parameter to characterize the movement of receptors in the cell membrane. The aim of this study was to investigate the differences between maximum likelihood-based estimation (MLE) and mean square displacement (MSD) based calculation of diffusion coefficients. In this work we applied both MSD and MLE to calculate diffusion coefficients. Single particle trajectories were extracted from simulation as well as from α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) receptor tracking in live estradiol-treated differentiated PC12 (dPC12) cells. The comparison of the obtained diffusion coefficients revealed the superiority of MLE over the generally used MSD analysis. Our results suggest the use of the MLE of diffusion coefficients because as it has a better performance, especially for large localization errors or slow receptor movements.

Asynchronous Reciprocal Coupling of Martini 2.2 Coarse-Grained and CHARMM36 All-Atom Simulations in an Automated Multiscale Framework

The appeal of multiscale modeling approaches is predicated on the promise of combinatorial synergy. However, this promise can only be realized when distinct scales are combined with reciprocal consistency. Here, we consider multiscale molecular dynamics (MD) simulations that combine the accuracy and macromolecular flexibility accessible to fixed-charge all-atom (AA) representations with the sampling speed accessible to reductive, coarse-grained (CG) representations. AA-to-CG conversions are relatively straightforward because deterministic routines with unique outcomes are achievable. Conversely, CG-to-AA conversions have many solutions due to a surge in the number of degrees of freedom. While automated tools for biomolecular CG-to-AA transformation exist, we find that one popular option, called Backward, is prone to stochastic failure and the AA models that it does generate frequently have compromised protein structure and incorrect stereochemistry. Although these shortcomings can likely be circumvented by human intervention in isolated instances, automated multiscale coupling requires reliable and robust scale conversion. Here, we detail an extension to Multiscale Machine-learned Modeling Infrastructure (MuMMI), including an improved CG-to-AA conversion tool called sinceCG. This tool is reliable (∼98% weakly correlated repeat success rate), automatable (no unrecoverable hangs), and yields AA models that generally preserve protein secondary structure and maintain correct stereochemistry. We describe how the MuMMI framework identifies CG system configurations of interest, converts them to AA representations, and simulates them at the AA scale while on-the-fly analyses provide feedback to update CG parameters. Application to systems containing the peripheral membrane protein RAS and proximal components of RAF kinase on complex eight-component lipid bilayers with ∼1.5 million atoms is discussed in the context of MuMMI.

Classification-based motion analysis of single-molecule trajectories using DiffusionLab

Single-particle tracking is a powerful approach to study the motion of individual molecules and particles. It can uncover heterogeneities that are invisible to ensemble techniques, which places it uniquely among techniques to study mass transport. Analysis of the trajectories obtained with single-particle tracking in inorganic porous hosts is often challenging, because trajectories are short and/or motion is heterogeneous. We present the DiffusionLab software package for motion analysis of such challenging data sets. Trajectories are first classified into populations with similar characteristics to which the motion analysis is tailored in a second step. DiffusionLab provides tools to classify trajectories based on the motion type either with machine learning or manually. It also offers quantitative mean squared displacement analysis of the trajectories. The software can compute the diffusion constant for an individual trajectory if it is sufficiently long, or the average diffusion constant for multiple shorter trajectories. We demonstrate the DiffusionLab approach via the analysis of a simulated data set with motion types frequently observed in inorganic porous hosts, such as zeolites. The software package with graphical user interface and its documentation are freely available.

MPTHub: An Open-Source Software for Characterizing the Transport of Particles in Biorelevant Media

The study of particle transport in different environments plays an essential role in understanding interactions with humans and other living organisms. Importantly, obtained data can be directly used for multiple applications in fields such as fundamental biology, toxicology, or medicine. Particle movement in biorelevant media can be readily monitored using microscopy and converted into time-resolved trajectories using freely available tracking software. However, translation into tangible and meaningful parameters is time consuming and not always intuitive. We developed new software—MPTHub—as an open-access, standalone, user-friendly tool for the rapid and reliable analysis of particle trajectories extracted from video microscopy. The software was programmed using Python and allowed to import and analyze trajectory data, as well as to export relevant data such as individual and ensemble time-averaged mean square displacements and effective diffusivity, and anomalous transport exponent. Data processing was reliable, fast (total processing time of less than 10 s), and required minimal memory resources (up to a maximum of around 150 MB in random access memory). Demonstration of software applicability was conducted by studying the transport of different polystyrene nanoparticles (100–200 nm) in mucus surrogates. Overall, MPTHub represents a freely available software tool that can be used even by inexperienced users for studying the transport of particles in biorelevant media.

Recapitulation of cell-like KRAS4b membrane dynamics on complex biomimetic membranes

Understanding the spatiotemporal distribution and dynamics of RAS on the plasma membrane (PM) is the key for elucidating the molecular mechanisms of the RAS signaling pathway. Single particle tracking (SPT) experiments show that in cells, KRAS diffuses in at least three interchanging states on the cellular PM; however, KRAS remains monomeric and always shows homogeneous diffusion on artificial membranes. Here, we show for the first time on a supported lipid bilayer composed of heterogeneous lipid components that we can recapitulate the three-state diffusion of KRAS seen in cells. The use of a biologically relevant eight-lipid system opens a new frontier in the biophysical studies of RAS and other membrane associated proteins on a biomimetic system that recapitulates the complexity of a cellular PM.

•

KRAS4b shows homogeneous diffusion on simple 2-lipids bilayer

•

KRAS4b shows a cell-like, three-state diffusion on a complex 8-lipid bilayer

•

Phase separation in lipids favors the multi-state diffusion of KRAS4b

•

The complex lipid composition favors RAS nanoclustering irrespective of nucleotide state

Taking Particle Tracking into Practice by Novel Software and Screening Approach: Case-Study of Oral Lipid Nanocarriers

The success on the design of new oral nanocarriers greatly depends on the identification of the best physicochemical properties that would allow their diffusion across the mucus layer that protects the intestinal epithelium. In this context, particle tracking (PT) has arisen in the pharmaceutical field as an excellent tool to evaluate the diffusion of individual particles across the intestinal mucus. In PT, the trajectories of individual particles are characterized by the mean square displacement (MSD), which is used to calculate the coefficient of diffusion (D) and the anomalous diffusion parameter (α) as MSD=4Dτα. Unfortunately, there is no stablished criteria to evaluate the goodness-of-fit of the experimental data to the mathematical model. This work shows that the commonly used R2 parameter may lead to an overestimation of the diffusion capacity of oral nanocarriers. We propose a screening approach based on a combination of R2 with further statistical parameters. We have analyzed the effect of this approach to study the intestinal mucodiffusion of lipid oral nanocarriers, compared to the conventional screening approach. Last, we have developed software able to perform the whole PT analysis in a time-saving, user-friendly, and rational fashion.

Single-molecule analysis of dynamics and interactions of the SecYEG translocon

Protein translocation and insertion into the bacterial cytoplasmic membrane are the essential processes mediated by the Sec machinery. The core machinery is composed of the membrane-embedded translocon SecYEG that interacts with the secretion-dedicated ATPase SecA and translating ribosomes. Despite the simplicity and the available structural insights on the system, diverse molecular mechanisms and functional dynamics have been proposed. Here, we employ total internal reflection fluorescence microscopy to study the oligomeric state and diffusion of SecYEG translocons in supported lipid bilayers at the single-molecule level. Silane-based coating ensured the mobility of lipids and reconstituted translocons within the bilayer. Brightness analysis suggested that approx. 70% of the translocons were monomeric. The translocons remained in a monomeric form upon ribosome binding, but partial oligomerization occurred in the presence of nucleotide-free SecA. Individual trajectories of SecYEG in the lipid bilayer revealed dynamic heterogeneity of diffusion, as translocons commonly switched between slow and fast mobility modes with corresponding diffusion coefficients of 0.03 and 0.7 µm² ·s^-1 . Interactions with SecA ATPase had a minor effect on the lateral mobility, while bound ribosome:nascent chain complexes substantially hindered the diffusion of single translocons. Notably, the mobility of the translocon:ribosome complexes was not affected by the solvent viscosity or macromolecular crowding modulated by Ficoll PM 70, so it was largely determined by interactions within the lipid bilayer and at the interface. We suggest that the complex mobility of SecYEG arises from the conformational dynamics of the translocon and protein:lipid interactions.

Membrane interactions of the globular domain and the hypervariable region of KRAS4b define its unique diffusion behavior

The RAS proteins are GTP-dependent switches that regulate signaling pathways and are frequently mutated in cancer. RAS proteins concentrate in the plasma membrane via lipid-tethers and hypervariable region side-chain interactions in distinct nano-domains. However, little is known about RAS membrane dynamics and the details of RAS activation of downstream signaling. Here, we characterize RAS in live human and mouse cells using single-molecule-tracking methods and estimate RAS mobility parameters. KRAS4b exhibits confined mobility with three diffusive states distinct from the other RAS isoforms (KRAS4a, NRAS, and HRAS); and although most of the amino acid differences between RAS isoforms lie within the hypervariable region, the additional confinement of KRAS4b is largely determined by the protein’s globular domain. To understand the altered mobility of an oncogenic KRAS4b, we used complementary experimental and molecular dynamics simulation approaches to reveal a detailed mechanism.

TCR-pMHC bond conformation controls TCR ligand discrimination

A major unanswered question is how a TCR discriminates between foreign and self-peptides presented on the APC surface. Here, we used in situ fluorescence resonance energy transfer (FRET) to measure the distances of single TCR–pMHC bonds and the conformations of individual TCR–CD3ζ receptors at the membranes of live primary T cells. We found that a TCR discriminates between closely related peptides by forming single TCR–pMHC bonds with different conformations, and the most potent pMHC forms the shortest bond. The bond conformation is an intrinsic property that is independent of the binding affinity and kinetics, TCR microcluster formation, and CD4 binding. The bond conformation dictates the degree of CD3ζ dissociation from the inner leaflet of the plasma membrane via a positive calcium signaling feedback loop to precisely control the accessibility of CD3ζ ITAMs for phosphorylation. Our data revealed the mechanism by which a TCR deciphers the structural differences among peptides via the TCR–pMHC bond conformation.

Cell morphology and nucleoid dynamics in dividing Deinococcus radiodurans

Our knowledge of bacterial nucleoids originates mostly from studies of rod- or crescent-shaped bacteria. Here we reveal that Deinococcus radiodurans, a relatively large spherical bacterium with a multipartite genome, constitutes a valuable system for the study of the nucleoid in cocci. Using advanced microscopy, we show that D. radiodurans undergoes coordinated morphological changes at both the cellular and nucleoid level as it progresses through its cell cycle. The nucleoid is highly condensed, but also surprisingly dynamic, adopting multiple configurations and presenting an unusual arrangement in which oriC loci are radially distributed around clustered ter sites maintained at the cell centre. Single-particle tracking and fluorescence recovery after photobleaching studies of the histone-like HU protein suggest that its loose binding to DNA may contribute to this remarkable plasticity. These findings demonstrate that nucleoid organization is complex and tightly coupled to cell cycle progression in this organism.

Diaphanous-1 affects the nanoscale clustering and lateral diffusion of receptor for advanced glycation endproducts (RAGE)

The interactions between the cytoplasmic protein diaphanous-1 (Diaph1) and the receptor for advanced glycation endproducts (RAGE) drive the negative consequences of RAGE signaling in several disease processes. Reported in this work is how Diaph1 affects the nanoscale clustering and diffusion of RAGE measured using super-resolution stochastic optical reconstruction microscopy (STORM) and single particle tracking (SPT). Altering the Diaph1 binding site has a different impact on RAGE diffusion compared to when Diaph1 expression is reduced in HEK293 cells. In cells with reduced Diaph1 expression (RAGE-Diaph1^-/-), the average RAGE diffusion coefficient is increased by 35%. RAGE diffusion is known to be influenced by the dynamics of the actin cytoskeleton. Actin labeling shows that a reduced Diaph1 expression leads to cells with reduced filopodia density and length. In contrast, when two RAGE amino acids that interact with Diaph1 are mutated (RAGE_RQ/AA), the average RAGE diffusion coefficient is decreased by 16%. Since RAGE diffusion is slowed when the interaction between Diaph1 and RAGE is disrupted, the interaction of the two proteins results in faster RAGE diffusion. In both RAGE_RQ/AA and RAGE-Diaph1^-/- cells the number and size of RAGE clusters are decreased compared to cells expressing RAGE and native concentrations of Diaph1. This work shows that Diaph1 has a role in affecting RAGE clusters and diffusion.

Bacterial cell wall nanoimaging by autoblinking microscopy

Spurious blinking fluorescent spots are often seen in bacteria during single-molecule localization microscopy experiments. Although this 'autoblinking' phenomenon is widespread, its origin remains unclear. In Deinococcus strains, we observed particularly strong autoblinking at the periphery of the bacteria, facilitating its comprehensive characterization. A systematic evaluation of the contributions of different components of the sample environment to autoblinking levels and the in-depth analysis of the photophysical properties of autoblinking molecules indicate that the phenomenon results from transient binding of fluorophores originating mostly from the growth medium to the bacterial cell wall, which produces single-molecule fluorescence through a Point Accumulation for Imaging in Nanoscale Topography (PAINT) mechanism. Our data suggest that the autoblinking molecules preferentially bind to the plasma membrane of bacterial cells. Autoblinking microscopy was used to acquire nanoscale images of live, unlabeled D. radiodurans and could be combined with PALM imaging of PAmCherry-labeled bacteria in two-color experiments. Autoblinking-based super-resolved images provided insight into the formation of septa in dividing bacteria and revealed heterogeneities in the distribution and dynamics of autoblinking molecules within the cell wall.

Single-molecule studies of flavivirus envelope dynamics: Experiment and computation

Flaviviruses are simple enveloped viruses exhibiting complex structural and functional heterogeneities. Decades of research have provided crucial basic insights, antiviral medication and moderately successful gene therapy trials. The most infectious particle is, however, not always the most abundant one in a population, questioning the utility of classic ensemble-averaging virology approaches. Indeed, viral replication is often not particularly efficient, prone to errors or containing parallel routes. Here, we review different single-molecule sensitive fluorescence methods that are employed to investigate flaviviruses. In particular, we review how (i) time-resolved Förster resonance energy transfer (trFRET) was applied to probe dengue envelope conformations; (ii) FRET-fluorescence correlation spectroscopy to investigate dengue envelope intrinsic dynamics and (iii) single particle tracking to follow the path of dengue viruses in cells. We also discuss how such methods may be supported by molecular dynamics (MD) simulations over a range of spatio-temporal scales, to provide complementary data on the structure and dynamics of flaviviral systems. We describe recent improvements in multiscale MD approaches that allowed the simulation of dengue particle envelopes in near-atomic resolution. We hope this review is an incentive for setting up and applying similar single-molecule studies and combine them with MD simulations to investigate structural dynamics of entire flavivirus particles over the nanosecond-to-millisecond time-scale and follow viruses during infection in cells over milliseconds to minutes.

Interactions and Translational Dynamics of Phosphatidylinositol Bisphosphate (PIP2) Lipids in Asymmetric Lipid Bilayers

Phosphatidylinositol phosphate (PIP) lipids are critical to many cell signaling pathways, in part by acting as molecular beacons that recruit peripheral membrane proteins to specific locations within the plasma membrane. Understanding the biophysics of PIP-protein interactions is critical to developing a chemically detailed model of cell communication. Resolving such interactions is challenging, even in model membrane systems, because of the difficulty in preparing PIP-containing membranes with high fluidity and integrity. Here we report on a simple, vesicle-based protocol for preparing asymmetric supported lipid bilayers in which fluorescent PIP lipid analogues are found only on the top leaflet of the supported membrane facing the bulk solution. With this asymmetric distribution of lipids between the leaflets, the fluorescent signal from the PIP lipid analogue reports directly on interactions between the peripheral molecules and the top leaflet of the membrane. Asymmetric PIP-containing bilayers are an ideal platform to investigate the interaction of PIP with peripheral membrane proteins using fluorescence-based imaging approaches. We demonstrate their usefulness here with a combined fluorescence correlation spectroscopy and single particle tracking study of the interaction between PIP2 lipids and a polycationic polymer, quaternized polyvinylpyridine (QPVP). With this approach we are able to quantify the microscopic features of the mobility coupling between PIP2 lipids and polybasic QPVP. With single particle tracking we observe individual PIP2 lipids switch from Brownian to intermittent motion as they become transiently trapped by QPVP.

In Situ Mapping of the Mechanical Properties of Biofilms by Particle-tracking Microrheology

Bacterial cells are able to form surface-attached biofilm communities known as biofilms by encasing themselves in extracellular polymeric substances (EPS). The EPS serves as a physical and protective scaffold that houses the bacterial cells and consists of a variety of materials that includes proteins, exopolysaccharides and DNA. The composition of the EPS may change, which remodels the mechanic properties of the biofilm to further develop or support alternative biofilm structures, such as streamers, as a response to environmental cues. Despite this, there are little quantitative descriptions on how EPS components contribute to the mechanical properties and function of biofilms. Rheology, the study of the flow of matter, is of particular relevance to biofilms as many biofilms grow in flow conditions and are constantly exposed to shear stress. It also provides measurement and insight on the spreading of the biofilm on a surface. Here, particle-tracking microrheology is used to examine the viscoelasticity and effective crosslinking roles of different matrix components in various parts of the biofilm during development. This approach allows researchers to measure mechanic properties of biofilms at the micro-scale, which might provide useful information for controlling and engineering biofilms.

Increased Microbial Butanol Tolerance by Exogenous Membrane Insertion Molecules

Butanol is an ideal biofuel, although poor titers lead to high recovery costs by distillation. Fluidization of microbial membranes by butanol is one of the major factors limiting titers in butanol-producing bioprocesses. Starting with the hypothesis that certain membrane insertion molecules would stabilize the lipid bilayer in the presence of butanol, we applied a combination of in vivo and in vitro techniques within an in silico framework to describe a new approach to achieve solvent tolerance in bacteria. Single-molecule tracking of a model supported bilayer showed that COE1-5C, a five-ringed oligo-polyphenylenevinylene conjugated oligoelectrolyte (COE), reduced the diffusion rate of phospholipids in a microbially derived lipid bilayer to a greater extent than three-ringed and four-ringed COEs. Furthermore, COE1-5C treatment increased the specific growth rate of E. coli K12 relative to a control at inhibitory butanol concentrations. Consequently, to confer butanol tolerance to microbes by exogenous means is complementary to genetic modification of strains in industrial bioprocesses, extends the physiological range of microbes to match favorable bioprocess conditions, and is amenable with complex and undefined microbial consortia for biobutanol production. Molecular dynamics simulations indicated that the π-conjugated aromatic backbone of COE1-5C likely acts as a hydrophobic tether for glycerophospholipid acyl chains by enhancing bilayer integrity in the presence of high butanol concentrations, which thereby counters membrane fluidization. COE1-5C-mitigated E. coli K12 membrane depolarization by butanol is consistent with the hypothesis that improved growth rates in the presence of butanol are a consequence of improved bilayer stability.

Fluorescence methods for analysis of interactions between Ca(2+) signaling, lysosomes, and endoplasmic reticulum

The endoplasmic reticulum (ER) is both the major source of intracellular Ca(2+) for cell signaling and the organelle that forms the most extensive contacts with the plasma membrane and other organelles. Lysosomes fulfill important roles in degrading cellular materials and in cholesterol handling, but they also contribute to Ca(2+) signaling by both releasing and sequestering Ca(2+). Interactions between ER and other Ca(2+)-transporting membranes, notably mitochondria and the plasma membrane, often occur at sites where the two membranes are closely apposed, allowing local Ca(2+) signaling between them. These interactions are often facilitated by scaffold proteins. Recent evidence suggests similar local interactions between ER and lysosomes. We describe simple fluorescence-based methods that allow the interplay between Ca(2+) signals, the ER, and lysosomes to be examined.

Preparation of mica supported lipid bilayers for high resolution optical microscopy imaging

Supported lipid bilayers (SLBs) are widely used as a model for studying membrane properties (phase separation, clustering, dynamics) and its interaction with other compounds, such as drugs or peptides. However SLB characteristics differ depending on the support used. Commonly used techniques for SLB imaging and measurements are single molecule fluorescence microscopy, FCS and atomic force microscopy (AFM). Because most optical imaging studies are carried out on a glass support, while AFM requires an extremely flat surface (generally mica), results from these techniques cannot be compared directly, since the charge and smoothness properties of these materials strongly influence diffusion. Unfortunately, the high level of manual dexterity required for the cutting and gluing thin slices of mica to the glass slide presents a hurdle to routine use of mica for SLB preparation. Although this would be the method of choice, such prepared mica surfaces often end up being uneven (wavy) and difficult to image, especially with small working distance, high numerical aperture lenses. Here we present a simple and reproducible method for preparing thin, flat mica surfaces for lipid vesicle deposition and SLB preparation. Additionally, our custom made chamber requires only very small volumes of vesicles for SLB formation. The overall procedure results in the efficient, simple and inexpensive production of high quality lipid bilayer surfaces that are directly comparable to those used in AFM studies.