The 100 most-cited articles in Biophysical Journal, 2023

Simulations of biomolecules have enormous potential to inform our understanding of biology but require extremely demanding calculations. For over twenty years, the Folding@home distributed computing project has pioneered a massively parallel approach to biomolecular simulation, harnessing the resources of citizen scientists across the globe. Here, we summarize the scientific and technical advances this perspective has enabled. As the project's name implies, the early years of Folding@home focused on driving advances in our understanding of protein folding by developing statistical methods for capturing long-timescale processes and facilitating insight into complex dynamical processes. Success laid a foundation for broadening the scope of Folding@home to address other functionally relevant conformational changes, such as receptor signaling, enzyme dynamics, and ligand binding. Continued algorithmic advances, hardware developments such as GPU-based computing, and the growing scale of Folding@home have enabled the project to focus on new areas where massively parallel sampling can be…

Protein aggregates, formed from the assembly of aberrant, misfolded proteins, are a hallmark of neurodegenerative diseases. Disease-associated aggregates such as mutant Huntingtin polyQ inclusions, are typically enriched in p62/SQSTM1, an oligomeric protein that binds to and sequesters aberrant proteins. p62 has been suggested to sequester proteins through formation of liquid-like biomolecular condensates, but the physical mechanisms by which p62 condensates may regulate pathological protein aggregation remain unclear. Here, we use a light-inducible biomimetic condensate system to show that p62 condensates enhance coarsening of mutant polyQ aggregates through interface-mediated sequestration, which accelerates polyQ accumulation into larger aggregates. However, the resulting large aggregates accumulate polyubiquitinated proteins, which depletes free p62, ultimately suppressing further p62 condensation. This dynamic interplay between interface-mediated coarsening of solid aggregates and downstream consequences on the phase behavior of associated regulatory proteins could contribute…

Abstract The cytoplasm is a complex, crowded, actively-driven environment whose biophysical characteristics modulate critical cellular processes such as cytoskeletal dynamics, phase separation, and stem-cell fate. Little is known about the variance in these cytoplasmic properties. Here, we employed particle-tracking nanorheology on genetically encoded multimeric 40-nm nanoparticles (GEMs) to measure diffusion within the cytoplasm of the fission yeast Schizosaccharomyces pombe . We found that the apparent diffusion coefficients of individual GEM particles varied over a 400-fold range, while the differences in average particle diffusivity among individual cells spanned a 10-fold range. To determine the origin of this heterogeneity, we developed a Doppelgänger Simulation approach that uses stochastic simulations of GEM diffusion that replicate the experimental statistics on a particle-by-particle basis, such that each experimental track and cell had a one-to-one correspondence with their simulated counterpart. These simulations showed that the large intra- and inter-cellular variations…

The RNA-binding protein TDP-43 is associated with mRNA processing and transport from the nucleus to the cytoplasm. TDP-43 localizes in the nucleus as well as accumulating in cytoplasmic condensates such as stress granules. Aggregation and formation of amyloid-like fibrils of cytoplasmic TDP-43 are hallmarks of numerous neurodegenerative diseases, most strikingly present in >90% of amyotrophic lateral sclerosis (ALS) patients. If excessive accumulation of cytoplasmic TDP-43 causes, or is caused by, neurodegeneration is presently not known. In this work, we use molecular dynamics simulations at multiple resolutions to explore TDP-43 self- and cross-interaction dynamics. A full-length molecular model of TDP-43, all 414 amino acids, was constructed from select structures of the protein functional domains (N-terminal domain, and two RNA recognition motifs, RRM1 and RRM2) and modeling of disordered connecting loops and the low complexity glycine-rich C-terminus domain. All-atom CHARMM36m simulations of single TDP-43 proteins served as guides to construct a coarse-grained Martini 3 model…

The fusion of lipid membranes progresses through a series of hemifusion intermediates with two significant energy barriers related to the formation of stalk and fusion pore, respectively. These energy barriers determine the speed and success rate of many critical biological processes, including the fusion of highly curved membranes, for example synaptic vesicles and enveloped viruses. Here we use continuum elastic theory of lipid monolayers to determine the relationship between membrane shape and energy barriers to fusion. We find that the stalk formation energy decreases with curvature by up to 31 k B T in a 20-nm-radius vesicle compared with planar membranes and by up to 8 k B T in the fusion of highly curved, long, tubular membranes. In contrast, the fusion pore formation energy barrier shows a more complicated behavior. Immediately after stalk expansion to the hemifusion diaphragm, the fusion pore formation energy barrier is low (15–25 k B T) due to lipid stretching in the distal monolayers and increased tension in highly curved vesicles. Therefore, the opening of the fusion…

The actin cortex is a complex cytoskeletal machinery that drives and responds to changes in cell shape. It must generate or adapt to plasma membrane curvature to facilitate diverse functions such as cell division, migration, and phagocytosis. Due to the complex molecular makeup of the actin cortex, it remains unclear whether actin networks are inherently able to sense and generate membrane curvature, or whether they rely on their diverse binding partners to accomplish this. Here, we show that curvature sensing is an inherent capability of branched actin networks nucleated by Arp2/3 and VCA. We develop a robust method to encapsulate actin inside giant unilamellar vesicles (GUVs) and assemble an actin cortex at the inner surface of the GUV membrane. We show that actin forms a uniform and thin cortical layer when present at high concentration and distinct patches associated with negative membrane curvature at low concentration. Serendipitously, we find that the GUV production method also produces dumbbell-shaped GUVs, which we explain using mathematical modeling in terms of membrane…

Biomolecular condensates, thought to form via liquid–liquid phase separation of intracellular mixtures, are multicomponent systems that can include diverse types of proteins and RNAs. RNA is a critical modulator of RNA-protein condensate stability, as it induces an RNA-concentration dependent reentrant phase transition—increasing stability at low RNA concentrations and decreasing it at high concentrations. Beyond concentration, RNAs inside condensates can be heterogeneous in length, sequence, and structure. Here, we use multiscale simulations to understanding how different RNA parameters interact with one another to modulate the properties of RNA-protein condensates. To do so, we perform residue/nucleotide-resolution coarse-grained Molecular Dynamics simulations of multicomponent RNA-protein condensates containing RNAs of different lengths and concentrations, and either FUS or PR 25 proteins. Our simulations reveal that RNA length regulates the reentrant phase behaviour of RNA-protein condensates: increasing RNA length sensitively rises the maximum value that the critical…

Cell science has made significant progress by focusing on understanding individual cellular processes through reductionist approaches. However, the sheer volume of knowledge collected presents challenges in integrating this information across different scales of space and time to comprehend cellular behaviors, as well as making the data and methods more accessible for the community to tackle complex biological questions. This perspective proposes the creation of next-generation virtual cells, which are dynamic 3D models that integrate information from diverse sources, including simulations, biophysical models, image-based models, and evidence-based knowledge graphs. These virtual cells would provide statistically accurate and holistic views of real cells, bridging the gap between theoretical concepts and experimental data, and facilitating productive new collaborations among researchers across related fields.

disappeared in ordinary metabolically active cells, implying that living cells fluidize their internal states such that they approach the critical jamming point.

The snapshot distribution of mRNA counts per cell can be measured using single-molecule fluorescence in situ hybridization or single-cell RNA sequencing. These distributions are often fit to the steady-state distribution of the two-state telegraph model to estimate the three transcriptional parameters for a gene of interest: mRNA synthesis rate, the switching on rate (the on state being the active transcriptional state), and the switching off rate. This model assumes no extrinsic noise, i.e., parameters do not vary between cells, and thus estimated parameters are to be understood as approximating the average values in a population. The accuracy of this approximation is currently unclear. Here, we develop a theory that explains the size and sign of estimation bias when inferring parameters from single-cell data using the standard telegraph model. We find specific bias signatures depending on the source of extrinsic noise (which parameter is most variable across cells) and the mode of transcriptional activity. If gene expression is not bursty then the population averages of all three…

To characterize the mechanisms governing the diffusion of particles in biological scenarios, it is essential to accurately determine their diffusive properties. To do so, we propose a machine-learning method to characterize diffusion processes with time-dependent properties at the experimental time resolution. Our approach operates at the single-trajectory level predicting the properties of interest, such as the diffusion coefficient or the anomalous diffusion exponent, at every time step of the trajectory. In this way, changes in the diffusive properties occurring along the trajectory emerge naturally in the prediction and thus allow the characterization without any prior knowledge or assumption about the system. We first benchmark the method on synthetic trajectories simulated under several conditions. We show that our approach can successfully characterize both abrupt and continuous changes in the diffusion coefficient or the anomalous diffusion exponent. Finally, we leverage the method to analyze experiments of single-molecule diffusion of two membrane proteins in living cells:…

The distribution of red blood cells (RBCs) in the microcirculation determines the oxygen delivery and solute transport to tissues. This process relies on the partitioning of RBCs at successive bifurcations throughout the microvascular network, and it has been known since the last century that RBCs partition disproportionately to the fractional blood flow rate, therefore leading to heterogeneity of the hematocrit (i.e., volume fraction of RBCs in blood) in microvessels. Usually, downstream of a microvascular bifurcation, the vessel branch with a higher fraction of blood flow receives an even higher fraction of RBC flux. However, both temporal and time-average deviations from this phase-separation law have been observed in recent studies. Here, we quantify how the microscopic behavior of RBC lingering (i.e., RBCs temporarily residing near the bifurcation apex with diminished velocity) influences their partitioning, through combined in vivo experiments and in silico simulations. We developed an approach to quantify the cell lingering at highly confined capillary-level bifurcations and…

ABSTRACT Upon nutrient limitation, budding yeast of Saccharomyces cerevisiae shift from fast growth (the log stage) to quiescence (the stationary stage). This shift is accompanied by liquid-liquid phase separation in the membrane of the vacuole, an endosomal organelle. Recent work indicates that the resulting micron-scale domains in vacuole membranes enable yeast to survive periods of stress. An outstanding question is which molecular changes might cause this membrane phase separation. Here, we conduct lipidomics of vacuole membranes in both the log and stationary stages. Isolation of pure vacuole membranes is challenging in the stationary stage, when lipid droplets are in close contact with vacuoles. Immuno-isolation has previously been shown to successfully purify log-stage vacuole membranes with high organelle specificity, but it was not previously possible to immuno-isolate stationary stage vacuole membranes. Here, we develop Mam3 as a bait protein for vacuole immuno-isolation, and demonstrate low contamination by non-vacuolar membranes. We find that stationary stage vacuole…

Biological condensates are known to retain a large fraction of water to remain in a liquid and reversible state. Local solvation contributions from water hydrating hydrophilic and hydrophobic protein surfaces were proposed to play a prominent role for the formation of condensates through liquid-liquid phase separation (LLPS). However, although the total free energy is accessible by calorimetry, the partial solvent contributions to the free energy changes upon LLPS remained experimentally inaccessible so far. Here, we show that the recently developed THz calorimetry approach allows to quantify local hydration enthalpy and entropy changes upon LLPS of α-elastin in real time, directly from experimental THz spectroscopy data. We find that hydrophobic solvation dominates the entropic solvation term, whereas hydrophilic solvation mainly contributes to the enthalpy. Both terms are in the order of hundreds of kJ/mol, which is more than one order of magnitude larger than the total free energy changes at play during LLPS. However, since we show that entropy/enthalpy mostly compensates, a…

The physical stability of peptide-based drugs is of great interest to the pharmaceutical industry. Glucagon-like peptide 1 (GLP-1) is a 31-amino acid peptide hormone, the analogs of which are frequently used in the treatment of type 2 diabetes. We investigated the physical stability of GLP-1 and its C-terminal amide derivative, GLP-1-Am, both of which aggregate into amyloid fibrils. While off-pathway oligomers have been proposed to explain the unusual aggregation kinetics observed previously for GLP-1 under specific conditions, these oligomers have not been studied in any detail. Such states are important as they may represent potential sources of cytotoxicity and immunogenicity. Here, we identified and isolated stable, low-molecular-weight oligomers of GLP-1 and GLP-1-Am, using size-exclusion chromatography. Under the conditions studied, isolated oligomers were shown to be resistant to fibrillation or dissociation. These oligomers contain between two and five polypeptide chains and they have a highly disordered structure as indicated by a variety of spectroscopic techniques. They…

Adherent cells use actomyosin contractility to generate mechanical force and to sense the physical properties of their environment, with dramatic consequences for migration, division, differentiation, and fate. However, the organization of the actomyosin system within cells is highly variable, with its assembly and function being controlled by small GTPases from the Rho family. To understand better how activation of these regulators translates into cell-scale force generation in the context of different physical environments, here we combine recent advances in non-neuronal optogenetics with micropatterning and traction force microscopy on soft elastic substrates. We find that, after whole-cell RhoA activation by the CRY2/CIBN optogenetic system with a short pulse of 100 ms, single cells contract on a minute timescale in proportion to their original traction force, before returning to their original tension setpoint with near perfect precision, on a longer timescale of several minutes. To decouple the biochemical and mechanical elements of this response, we introduce a mathematical…

Computer simulation techniques form a versatile tool, a computational microscope, for exploring biological processes. This tool has been particularly effective in exploring different features of biological membranes. In recent years, thanks to elegant multiscale simulation schemes, some fundamental limitations of investigations by distinct simulation techniques have been resolved. As a result, we are now capable of exploring processes spanning multiple scales beyond the capacity of any single technique. In this perspective, we argue that mesoscale simulations require more attention and must be further developed to fill evident gaps in a quest toward simulating and modeling living cell membranes.

Multiphasic architectures are found ubiquitously in biomolecular condensates and are thought to have important implications for the organization of multiple chemical reactions within the same compartment. Many of these multiphasic condensates contain RNA in addition to proteins. Here, we investigate the importance of different interactions in multiphasic condensates comprising two different proteins and RNA using computer simulations with a residue-resolution coarse-grained model of proteins and RNA. We find that in multilayered condensates containing RNA in both phases, protein-RNA interactions dominate, with aromatic residues and arginine forming the key stabilizing interactions. The total aromatic and arginine content of the two proteins must be appreciably different for distinct phases to form, and we show that this difference increases as the system is driven toward greater multiphasicity. Using the trends observed in the different interaction energies of this system, we demonstrate that we can also construct multilayered condensates with RNA preferentially concentrated in one…

Abstract When multivesicular endosomes (MVEs) fuse with the plasma membrane, exosomes are released into the extracellular space where they can affect other cells. Whether exosomes regulate cells nearby or further away depends on whether they remain attached to the secreting cell membrane. The regulation and kinetics of exosome secretion are not well characterized, but probes for directly imaging single MVE fusion events have allowed for visualization of the fusion and release process. In particular, the design of an exosome marker with a pH sensitive dye in the middle of the tetraspanin protein CD63 has facilitated studies of individual MVE fusion events. Using TIRF microscopy, single MVE fusion events were measured in A549 cells held at 23-37°C and events were identified using an automated detection algorithm. Stable docking precedes fusion almost all of the time and a decrease in temperature was accompanied by decrease in the rate of content loss and a decrease in the frequency of fusion events. The loss of CD63-pHluorin fluorescence was measured at fusion sites and fit with a…

NADH and NADPH play key roles in the regulation of metabolism. Their endogenous fluorescence is sensitive to enzyme binding, allowing changes in cellular metabolic state to be determined using fluorescence lifetime imaging microscopy (FLIM). However, to fully uncover the underlying biochemistry, the relationships between their fluorescence and binding dynamics require greater understanding. Here we accomplish this through time- and polarization-resolved fluorescence and polarized two-photon absorption measurements. Two lifetimes result from binding of both NADH to lactate dehydrogenase and NADPH to isocitrate dehydrogenase. The composite fluorescence anisotropy indicates the shorter (1.3-1.6 ns) decay component to be accompanied by local motion of the nicotinamide ring, pointing to attachment solely via the adenine moiety. For the longer lifetime (3.2-4.4 ns), the nicotinamide conformational freedom is found to be fully restricted. As full and partial nicotinamide binding are recognized steps in dehydrogenase catalysis, our results unify photophysical, structural, and functional…

Islets of Langerhans operate as multicellular networks in which several hundred β cells work in synchrony to produce secretory pulses of insulin, a hormone crucial for controlling metabolic homeostasis. Their collective rhythmic activity is facilitated by gap junctional coupling and affected by their functional heterogeneity, but the details of this robust and coordinated behavior are still not fully understood. Recent advances in multicellular imaging and optogenetic and photopharmacological strategies, as well as in network science, have led to the discovery of specialized β cell subpopulations that were suggested to critically determine the collective dynamics in the islets. In particular hubs, i.e., β cells with many functional connections, are believed to significantly enhance communication capacities of the intercellular network and facilitate an efficient spreading of intercellular Ca 2+ waves, whereas wave-initiator cells trigger intercellular signals in their cohorts. Here, we determined Ca 2+ signaling characteristics of these two β cell subpopulations and the relationship…

In biology, release of Ca 2+ ions in the cytosol is essential to trigger or control many cell functions. Calcium signaling acutely depends on lipid membrane permeability to Ca 2+ . For proper understanding of membrane permeability to Ca 2+ , both membrane hydration and the structure of the hydrophobic core must be taken into account. Here, we vary the hydrophobic core of bilayer membranes and observe different types of behavior in high-throughput wide-field second harmonic imaging. Ca 2+ translocation is observed through mono-unsaturated (DOPC:DOPA) membranes, reduced upon the addition of cholesterol, and completely inhibited for branched (DPhPC:DPhPA) and poly-unsaturated (SLPC:SLPA) lipid membranes. We propose, using molecular dynamics simulations, that ion transport occurs through ion-induced transient pores, which requires nonequilibrium membrane restructuring. This results in different rates at different locations and suggests that the hydrophobic structure of lipids plays a much more sophisticated regulating role than previously thought.

The innate immune system provides a crucial first line of defense against invading pathogens attacking the body. As the only member of the human cathelicidin family, the antimicrobial peptide LL-37 has been shown to have antiviral, antifungal, and antibacterial properties. In complexation with nucleic acids, LL-37 is suggested to maintain its beneficial health effects while also acting as a condensation agent for the nucleic acid. Complexes formed by LL-37 and nucleic acids have been shown to be immunostimulatory with a positive impact on the human innate immune system. However, some studies also suggest that in some circumstances, LL-37/nucleic acid complexes may be a contributing factor to autoimmune disorders such as psoriasis and systemic lupus erythematosus. This review provides a comprehensive discussion of research highlighting the beneficial health effects of LL-37/nucleic acid complexes, as well as discussing observed detrimental effects. We will emphasize why it is important to investigate and elucidate structural characteristics, such as condensation patterns of nucleic…

The heat shock protein 90 (Hsp90) is a molecular chaperone, which plays a key role in eukaryotic protein homeostasis. Co-chaperones assist Hsp90 in client maturation and in regulating essential cellular processes such as cell survival, signal transduction, gene regulation, hormone signaling, and neurodegeneration. Aha1 (activator of Hsp90 ATPase) is a unique co-chaperone known to stimulate the ATP hydrolysis of Hsp90, but the mechanism of their interaction is still unclear. In this report, we show that one or two Aha1 molecules can bind to one Hsp90 dimer and that the binding stoichiometry affects Hsp90's conformation, kinetics, ATPase activity, and stability. In particular, a coordination of two Aha1 molecules can be seen in stimulating the ATPase activity of Hsp90 and the unfolding of the middle domain, whereas the conformational equilibrium and kinetics are hardly affected by the stoichiometry of bound Aha1. Altogether, we show a regulation mechanism through the stoichiometry of Aha1 going far beyond a regulation of Hsp90's conformation.

We studied the mechanical leaflet coupling of prototypic mammalian plasma membranes using neutron spin-echo spectroscopy. In particular, we examined a series of asymmetric phospholipid vesicles with phosphatidylcholine and sphingomyelin enriched in the outer leaflet and inner leaflets composed of phosphatidylethanolamine/phosphatidylserine mixtures. The bending rigidities of most asymmetric membranes were anomalously high, exceeding even those of symmetric membranes formed from their cognate leaflets. Only asymmetric vesicles with outer leaflets enriched in sphingolipid displayed bending rigidities in conformity with these symmetric controls. We performed complementary small-angle neutron and x-ray experiments on the same vesicles to examine possible links to structural coupling mechanisms, which would show up in corresponding changes in membrane thickness. In addition, we estimated differential stress between leaflets originating either from a mismatch of their lateral areas or spontaneous curvatures. However, no correlation with asymmetry-induced membrane stiffening was observed.…

Cell division during early embryogenesis is linked to key morphogenic events such as embryo symmetry breaking and tissue patterning. It is thought that the physical surrounding of cells together with cell intrinsic cues act as a mechanical "mold," guiding cell division to ensure these events are robust. To quantify how cell division is affected by the mechanical and geometrical environment, we present a novel computational mechanical model of cytokinesis, the final phase of cell division. Simulations with the model reproduced experimentally observed furrow dynamics and describe the volume ratio of daughter cells in asymmetric cell divisions, based on the position and orientation of the mitotic spindle. For dividing cells in geometrically confined environments, we show how the orientation of confinement relative to the division axis modulates the volume ratio in asymmetric cell division. Further, we quantified how cortex viscosity and surface tension determine the shape of a dividing cell and govern bubble-instabilities in asymmetric cell division. Finally, we simulated the formation…

Imaging fluorescence correlation spectroscopy (FCS) is a powerful tool to extract information on molecular mobilities, actions, and interactions in live cells, tissues, and organisms. Nevertheless, several limitations restrict its applicability. First, FCS is data hungry, requiring 50,000 frames at 1-ms time resolution to obtain accurate parameter estimates. Second, the data size makes evaluation slow. Third, as FCS evaluation is model dependent, data evaluation is significantly slowed unless analytic models are available. Here, we introduce two convolutional neural networks-FCSNet and ImFCSNet-for correlation and intensity trace analysis, respectively. FCSNet robustly predicts parameters in 2D and 3D live samples. ImFCSNet reduces the amount of data required for accurate parameter retrieval by at least one order of magnitude and makes correct estimates even in moderately defocused samples. Both convolutional neural networks are trained on simulated data, are model agnostic, and allow autonomous, real-time evaluation of imaging FCS measurements.

Formation of the immunological synapse (IS) is a key event during initiation of an adaptive immune response to a specific antigen. During this process, a T cell and an antigen presenting cell form a stable contact that allows the T cell to integrate both internal and external stimuli in order to decide whether to activate. The threshold for T cell activation depends on the strength and frequency of the calcium (Ca 2+ ) signaling induced by antigen recognition, and it must be tightly regulated to avoid undesired harm to healthy cells. Potassium (K + ) channels are recruited to the IS to maintain the negative membrane potential required to sustain Ca 2+ entry. However, the precise localization of K + channels within the IS remains unknown. Here, we visualized the dynamic subsynaptic distribution of Kv1.3, the main voltage-gated potassium channel in human T cells. Upon T cell receptor engagement, Kv1.3 polarized toward the synaptic cleft and diffused throughout the F-actin rich distal compartment of the synaptic interface—an effect enhanced by CD2-CD58 corolla formation. As the synapse…

Transcription factor (TF) proteins regulate gene expression by binding to specific sites on the genome. In the facilitated diffusion model, an optimized search process is achieved by the TF alternating between 3D diffusion in the bulk and 1D diffusion along DNA. While undergoing 1D diffusion, the protein can switch from a search mode for fast diffusion along nonspecific DNA to a recognition mode for stable binding to specific DNA. It was recently noticed that, for a small TF domain protein, reorientations on DNA happen between the nonspecific and specific DNA binding. We here conducted all-atom molecular dynamics simulations with steering forces to reveal the protein-DNA binding free energetics, confirming that the search and recognition modes are distinguished primarily by protein orientations on the DNA. As the binding free energy difference between the specific and nonspecific DNA system slightly deviates from that being estimated directly from dissociation constants on 15-bp DNA constructs, we hypothesize that the discrepancy can come from DNA sequences flanking the 6-bp central…

Cells exert forces on mechanically compliant environments to sense stiffness, migrate, and remodel tissue. Cells can sense environmental stiffness via myosin-generated pulling forces acting on F-actin, which is in turn mechanically coupled to the environment via adhesive proteins, akin to a clutch in a drivetrain. In this "motor-clutch" framework, the force transmitted depends on the complex interplay of motor, clutch, and environmental properties. Previous mean-field analysis of the motor-clutch model identified the conditions for optimal stiffness for maximal force transmission via a dimensionless number that combines motor-clutch parameters. However, in this and other previous mean-field analyses, the motor-clutch system is assumed to have balanced motors and clutches and did not consider force-dependent clutch reinforcement and catch bond behavior. Here, we generalize the motor-clutch analytical framework to include imbalanced motor-clutch regimes, with clutch reinforcement and catch bonding, and investigate optimality with respect to all parameters. We found that traction force…

Long noncoding RNAs (lncRNAs) perform several important functions in cells including cis-regulation of transcription. Barring a few specific cases, the mechanisms underlying transcriptional regulation by lncRNAs remain poorly understood. Transcriptional proteins can form condensates via phase separation at protein-binding loci (BL) on the genome (e.g., enhancers and promoters). lncRNA-coding genes are present at loci in close genomic proximity of these BL and these RNAs can interact with transcriptional proteins via attractive heterotypic interactions mediated by their net charge. Motivated by these observations, we propose that lncRNAs can dynamically regulate transcription in cis via charge-based heterotypic interactions with transcriptional proteins in condensates. To study the consequences of this mechanism, we developed and studied a dynamical phase-field model. We find that proximal lncRNAs can promote condensate formation at the BL. Vicinally localized lncRNA can migrate to the BL to attract more protein because of favorable interaction free energies. However, increasing the…

Cell motility on flat substrates exhibits coexisting steady and oscillatory morphodynamics, the biphasic adhesion-velocity relation, and the universal correlation between speed and persistence (UCSP) as simultaneous observations common to many cell types. Their universality and concurrency suggest a unifying mechanism causing all three of them. Stick-slip models for cells on one-dimensional lanes suggest multistability to arise from the nonlinear friction of retrograde flow. This study suggests a mechanical mechanism controlled by integrin signaling on the basis of a biophysical model and analysis of trajectories of MDA-MB-231 cells on fibronectin lanes, which additionally explains the constitutive relations. The experiments exhibit cells with steady or oscillatory morphodynamics and either spread or moving with spontaneous transitions between the dynamic regimes, spread and moving, and spontaneous direction reversals. Our biophysical model is based on the force balance at the protrusion edge, the noisy clutch of retrograde flow, and a response function of friction and membrane drag…

Abstract Highly homologous ubiquitin-binding shuttle proteins UBQLN1, UBQLN2 and UBQLN4 differ in both their specific protein quality control functions and their propensities to localize to stress-induced condensates, cellular aggregates and aggresomes. We previously showed that UBQLN2 phase separates in vitro , and that the phase separation propensities of UBQLN2 deletion constructs correlate with their ability to form condensates in cells. Here, we demonstrated that full-length UBQLN1, UBQLN2 and UBQLN4 exhibit distinct phase behaviors in vitro . Strikingly, UBQLN4 phase separates at a much lower saturation concentration than UBQLN1. However, neither UBQLN1 nor UBQLN4 phase separates with a strong temperature dependence, unlike UBQLN2. We determined that the temperature-dependent phase behavior of UBQLN2 stems from its unique proline-rich (Pxx) region, which is absent in the other UBQLNs. We found that the short N-terminal disordered regions of UBQLN1, UBQLN2 and UBQLN4 inhibit UBQLN phase separation via electrostatics interactions. Charge variants of the N-terminal regions…

Embryonic development requires the establishment of cell polarity to enable cell fate segregation and tissue morphogenesis. This process is regulated by Par complex proteins, which partition into polarized membrane domains and direct downstream polarized cell behaviors. The kinase aPKC (along with its cofactor Par6) is a key member of this network and can be recruited to the plasma membrane by either the small GTPase Cdc42 or the scaffolding protein Par3. Although in vitro interactions among these proteins are well established, much is still unknown about the complexes they form during development. Here, to enable the study of membrane-associated complexes ex vivo, we used a maleic acid copolymer to rapidly isolate membrane proteins from single C. elegans zygotes into lipid nanodiscs. We show that native lipid nanodisc formation enables detection of endogenous complexes involving Cdc42, which are undetectable when cells are lysed in detergent. We found that Cdc42 interacts more strongly with aPKC/Par6 during polarity maintenance than polarity establishment, two developmental stages…

During the HIV-1 assembly process, the Gag polyprotein multimerizes at the producer cell plasma membrane, resulting in the formation of spherical immature virus particles. Gag-gRNA interactions play a crucial role in the multimerization process, which is yet to be fully understood. We have performed large-scale all-atom molecular dynamics simulations of membrane-bound full-length Gag dimer, hexamer, and 18-mer. The inter-domain dynamic correlation of Gag, quantified by the heterogeneous elastic network model (hENM) applied to the simulated trajectories, is observed to be altered by implicit gRNA binding, as well as the multimerization state of the Gag. The lateral dynamics of our simulated membrane-bound Gag proteins, with and without gRNA binding, agree with prior experimental data and help to validate our simulation models and methods. The gRNA binding is observed to impact mainly the SP1 domain of the 18-mer and the MA-CA linker domain of the hexamer. In the absence of gRNA binding, the independent dynamical motion of the NC domain results in a collapsed state of the dimeric Gag.…

Lysolipids such as lauroyl, myristoyl, and palmitoyl lysophosphatidylcholine (LPC) insert into the outer leaflet of liposomes but do not flip to the inner leaflet over many hours. This way, they create asymmetry stress between the intrinsic areas of the two leaflets. We have studied how this stress is relaxed with particular emphasis on the budding and fission of small (diameter 20-30 nm) daughter vesicles (DVs). Asymmetric flow field-flow fractionation was utilized to quantify the extent of budding from large unilamellar vesicles after exposure to LPC. Budding starts at a low threshold of the order of 2 mol% LPC in the outer (and ≈0 mol% LPC in the inner) leaflet. We see reason to assume that the fractional fluorescence intensity from DVs is a good approximation for the fraction of membrane lipid, POPC, transferred into DVs. Accordingly, budding starts with a "budding power" of ≈6 POPC molecules budding off per LPC added, corresponding to a more than 10-fold accumulation of LPC in the outer leaflet of DVs to ≈24 mol%. As long as budding is possible, little strain is built up in the…

The endoplasmic reticulum (ER) is a dynamic network of interconnected sheets and tubules that orchestrates the distribution of lipids, ions, and proteins throughout the cell. The impact of its complex, dynamic morphology on its function as an intracellular transport hub remains poorly understood. To elucidate the functional consequences of ER network structure and dynamics, we quantify how the heterogeneity of the peripheral ER in COS7 cells affects diffusive protein transport. In vivo imaging of photoactivated ER membrane proteins demonstrates their nonuniform spreading to adjacent regions, in a manner consistent with simulations of diffusing particles on extracted network structures. Using a minimal network model to represent tubule rearrangements, we demonstrate that ER network dynamics are sufficiently slow to have little effect on diffusive protein transport. Furthermore, stochastic simulations reveal a novel consequence of ER network heterogeneity: the existence of "hot spots" where sparse diffusive reactants are more likely to find one another. ER exit sites, specialized…