Comment in Biophys J. 2023 Sep 5;122(17):3544-3548.

E-cadherin plays a central role in cell-cell adhesion. The ectodomains of wild-type cadherins form a crystalline-like two-dimensional lattice in cell-cell interfaces mediated by both trans (apposed cell) and cis (same cell) interactions. In addition to these extracellular forces, adhesive strength is further regulated by cytosolic phenomena involving α and β catenin-mediated interactions between cadherin and the actin cytoskeleton. Cell-cell adhesion can be further strengthened under tension through mechanisms that have not been definitively characterized in molecular detail. Here we quantitatively determine the role of the cadherin ectodomain in mechanosensing. To this end, we devise an E-cadherin-coated emulsion system, in which droplet surface tension is balanced by protein binding strength to give rise to stable areas of adhesion. To reach the honeycomb/cohesive limit, an initial emulsion compression by centrifugation facilitates E-cadherin trans binding, whereas a high protein surface concentration enables the cis-enhanced stabilization of the interface. We observe an abrupt concentration dependence on recruitment into adhesions of constant crystalline density, reminiscent of a first-order phase transition. Removing the lateral cis interaction with a "cis mutant" shifts this transition to higher surface densities leading to denser, yet weaker adhesions. In both proteins, the stabilization of progressively larger areas of deformation is consistent with single-molecule experiments that show a force-dependent lifetime enhancement in the cadherin ectodomain, which may be attributed to the "X-dimer" bond.

In the presence of polyvalent cations, long double-stranded DNA (dsDNA) in dilute solution undergoes a single-molecule, first-order, phase transition ("condensation"), a phenomenon that has been documented and analyzed by many years of experimental and theoretical studies. There has been no systematic effort, however, to determine whether long single-stranded RNA (ssRNA) shows an analogous behavior. In this study, using dynamic light scattering, analytical ultracentrifugation, and gel electrophoresis, we examine the effects of increasing polyvalent cation concentrations on the effective size of long ssRNAs ranging from 3000 to 12,000 nucleotides. Our results indicate that ssRNA does not undergo a discontinuous condensation as does dsDNA but rather a "continuous" decrease in size with increasing polyvalent cation concentration. And, instead of the 10-fold decrease in size shown by long dsDNA, we document a 50% decrease, as demonstrated for a range of lengths and sequences of ssRNA.

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.

Chromosomal dynamics plays a central role in a number of critical biological processes, such as transcriptional regulation, genetic recombination, and DNA replication. However, visualization of chromatin is generally limited to live imaging of a few fluorescently labeled chromosomal loci or high-resolution reconstruction of multiple loci from a single time frame. To aid in mapping the underlying chromosomal structure based on parsimonious experimental measurements, we present an exact analytical expression for the evolution of the polymer configuration based on a flexible-polymer model, and we propose an algorithm that tracks the polymer configuration from live images of chromatin marked with several fluorescent marks. Our theory identifies the resolution of microscopy needed to achieve high-accuracy tracking for a given spacing of markers, establishing the statistical confidence in the assignment of genome identity to the visualized marks. We then leverage experimental data of locus-tracking measurements to demonstrate the validity of our modeling approach and to establish a basis for the design of experiments with a desired resolution. Altogether, this work provides a computational approach founded on polymer physics that vastly improves the interpretation of in vivo measurements of biopolymer dynamics.

Genomic stability in proliferating cells critically depends on telomere maintenance by telomerase reverse transcriptase. Here we report the development and proof-of-concept results of a single-molecule approach to monitor the catalytic activity of human telomerase in real time and with single-nucleotide resolution. Using zero-mode waveguides and multicolor FRET, we recorded the processive addition of multiple telomeric repeats to individual DNA primers. Unlike existing biophysical and biochemical tools, the novel approach enables the quantification of nucleotide-binding kinetics before nucleotide incorporation. Moreover, it provides a means to dissect the unique translocation dynamics that telomerase must undergo after synthesis of each hexameric DNA repeat. We observed an unexpectedly prolonged binding dwell time of dGTP in the enzyme active site at the start of each repeat synthesis cycle, suggesting that telomerase translocation is composed of multiple rate-contributing sub-steps that evade classical biochemical analysis.

Optical trapping in biophysics typically uses micron-scale beads made of materials like polystyrene or glass to probe the target of interest. Using smaller beads made of higher-index materials could increase the time resolution of these measurements. We characterized the trapping of nanoscale beads made of diamond and titanium dioxide (TiO2) in a single-beam gradient trap. Calculating theoretical expectations for the trapping stiffness of these beads, we found good agreement with measured values. Trap stiffness was significantly higher for TiO2 beads, owing to notable enhancement from nonlinear optical effects, not previously observed for continuous-wave trapping. Trap stiffness was over 6-fold higher for TiO2 beads than polystyrene beads of similar size at 70 mW laser power. These results suggest that diamond and TiO2 nanobeads can be used to improve time resolution in optical tweezers measurements.

Traction patterns of adherent cells provide important information on their interaction with the environment, cell migration, or tissue patterns and morphogenesis. Traction force microscopy is a method aimed at revealing these traction patterns for adherent cells on engineered substrates with known constitutive elastic properties from deformation information obtained from substrate images. Conventionally, the substrate deformation information is processed by numerical algorithms of varying complexity to give the corresponding traction field via solution of an ill-posed inverse elastic problem. We explore the capabilities of a deep convolutional neural network as a computationally more efficient and robust approach to solve this inversion problem. We develop a general purpose training process based on collections of circular force patches as synthetic training data, which can be subjected to different noise levels for additional robustness. The performance and the robustness of our approach against noise is systematically characterized for synthetic data, artificial cell models, and real cell images, which are subjected to different noise levels. A comparison with state-of-the-art Bayesian Fourier transform traction cytometry reveals the precision, robustness, and speed improvements achieved by our approach, leading to an acceleration of traction force microscopy methods in practical applications.

Ryanodine receptors (RyRs) are Ca2+ release channels, gated by Ca2+ in the cytosol and the sarcoplasmic reticulum lumen. Their regulation is impaired in certain cardiac and muscle diseases. Although a lot of data is available on the luminal Ca2+ regulation of RyR, its interpretation is complicated by the possibility that the divalent ions used to probe the luminal binding sites may contaminate the cytoplasmic sites by crossing the channel pore. In this study, we used Eu3+, an impermeable agonist of Ca2+ binding sites, as a probe to avoid this complication and to gain more specific information about the function of the luminal Ca2+ sensor. Single-channel currents were measured from skeletal muscle and cardiac RyRs (RyR1 and RyR2) using the lipid bilayer technique. We show that RyR2 is activated by the luminal addition of Ca2+, whereas RyR1 is inhibited. These results were qualitatively reproducible using Eu3+. The luminal regulation of RyR1 carrying a mutation associated with malignant hyperthermia was not different from that of the wild-type. RyR1 inhibition by Eu3+ was extremely voltage dependent, whereas RyR2 activation did not depend on the membrane potential. These results suggest that the RyR1 inhibition site is in the membrane's electric field (channel pore), whereas the RyR2 activation site is outside. Using in silico analysis and previous results, we predicted putative Ca2+ binding site sequences. We propose that RyR2 bears an activation site, which is missing in RyR1, but both isoforms share the same inhibitory Ca2+ binding site near the channel gate.

Using all-atom replica-exchange molecular dynamics simulations, we mapped the mechanisms of binding of the nuclear localization signal (NLS) sequence from Venezuelan equine encephalitis virus (VEEV) capsid protein to importin-α (impα) transport protein. Our objective was to identify the VEEV NLS sequence fragment that confers native, experimentally resolved binding to impα as well as to study associated binding energetics and conformational ensembles. The two selected VEEV NLS peptide fragments, KKPK and KKPKKE, show strikingly different binding mechanisms. The minNLS peptide KKPK binds non-natively and nonspecifically by adopting five diverse conformational clusters with low similarity to the x-ray structure 3VE6 of NLS-impα complex. Despite the prevalence of non-native interactions, the minNLS peptide still largely binds to the impα major NLS binding site. In contrast, the coreNLS peptide KKPKKE binds specifically and natively, adopting a largely homogeneous binding ensemble with a dominant, highly native-like conformational cluster. The coreNLS peptide retains most of native binding interactions, including π-cation contacts and a tryptophan cage. While KKPK binding is governed by a complex multistate free energy landscape featuring transitions between multiple binding poses, the coreNLS peptide free energy map is simple, exhibiting a single dominant native-like bound basin. We argue that the origin of the coreNLS peptide binding specificity is several electrostatic interactions formed by the two C-terminal amino acids, Lys10 and Glu11, with impα. The coreNLS sequence is then sufficient for native binding, but none of the amino acids flanking minNLS, including Lys10 and Glu11, are strictly necessary for the native pose. Our analyses indicate that the VEEV coreNLS sequence is virtually unique among human and viral proteins interacting with impα making it a potential target for VEEV-specific inhibitors.

Mechanical stresses generated at the cell-cell level and cell-substrate level have been suggested to be important in a host of physiological and pathological processes. However, the influence various chemical compounds have on the mechanical stresses mentioned above is poorly understood, hindering the discovery of novel therapeutics, and representing a barrier in the field. To overcome this barrier, we implemented two approaches: 1) monolayer boundary predictor and 2) discretized window predictor utilizing either stepwise linear regression or quadratic support vector machine machine learning model to predict the dose-dependent response of tractions and intercellular stresses to chemical perturbation. We used experimental traction and intercellular stress data gathered from samples subject to 0.2 or 2 μg/mL drug concentrations along with cell morphological properties extracted from the bright-field images as predictors to train our model. To demonstrate the predictive capability of our machine learning models, we predicted tractions and intercellular stresses in response to 0 and 1 μg/mL drug concentrations which were not utilized in the training sets. Results revealed the discretized window predictor trained just with four samples (292 images) to best predict both intercellular stresses and tractions using the quadratic support vector machine and stepwise linear regression models, respectively, for the unseen sample images.

Adherent filopodia are elongated finger-like membrane protrusions, extending from the edges of diverse cell types and participating in cell adhesion, spreading, migration, and environmental sensing. The formation and elongation of filopodia are driven by the polymerization of parallel actin filaments, comprising the filopodia cytoskeletal core. Here, we report that adherent filopodia, formed during the spreading of cultured cells on galectin-8-coated substrates, tend to change the direction of their extension in a chiral fashion, acquiring a left-bent shape. Cryoelectron tomography examination indicated that turning of the filopodia tip to the left is accompanied by the displacement of the actin core bundle to the right of the filopodia midline. Reduction of the adhesion to galectin-8 by treatment with thiodigalactoside abolished this filopodia chirality. By modulating the expression of a variety of actin-associated filopodia proteins, we identified myosin-X and formin DAAM1 as major filopodia chirality promoting factors. Formin mDia1, actin filament elongation factor VASP, and actin filament cross-linker fascin were also shown to be involved. Thus, the simple actin cytoskeleton of filopodia, together with a small number of associated proteins are sufficient to drive a complex navigation process, manifested by the development of left-right asymmetry in these cellular protrusions.

Cell migration is a complex phenomenon. Not only do different cells migrate in different default modes, but the same cell can also change its migration mode to adapt to different terrains. This complexity has riddled cell biologists and biophysicists for decades in that, despite the development of many powerful tools over the past 30 years, how cells move is still being actively investigated. This is because we have yet to fully understand the mystery of cell migration plasticity, particularly the reciprocal relation between force generation and migration mode transition. Herein we explore the future directions, in terms of measurement platforms and imaging-based techniques, to facilitate the undertaking of elucidating the relation between force generation machinery and migration mode transition. By briefly reviewing the evolution of the platforms and techniques developed in the past, we propose the desirable features to be added to achieve high measurement accuracy and improved temporal and spatial resolution, permitting us to unveil the mystery of cell migration plasticity.

Fluorescence recovery after photobleaching (FRAP) has emerged as one of the most widely utilized techniques to quantify binding and diffusion kinetics of biomolecules in biophysics. Since its inception in the mid-1970s, FRAP has been used to address an enormous array of questions including the characteristic features of lipid rafts, how cells regulate the viscosity of their cytoplasm, and the dynamics of biomolecules inside condensates formed by liquid-liquid phase separation. In this perspective, I briefly summarize the history of the field and discuss why FRAP has proven to be so incredibly versatile and popular. Next, I provide an overview of the extensive body of knowledge that has emerged on best practices for quantitative FRAP data analysis, followed by some recent examples of biological lessons learned using this powerful approach. Finally, I touch on new directions and opportunities for biophysicists to contribute to the continued development of this still-relevant research tool.

Organogenesis arises from the collective arrangement of cells into progressively 3D-shaped tissue. The acquisition of a correctly shaped organ is then the result of a complex interplay between molecular cues, responsible for differentiation and patterning, and the mechanical properties of the system, which generate the necessary forces that drive correct shape emergence. Nowadays, technological advances in the fields of microscopy, molecular biology, and computer science are making it possible to see and record such complex interactions in incredible, unforeseen detail within the global context of the developing embryo. A quantitative and interdisciplinary perspective of developmental biology becomes then necessary for a comprehensive understanding of morphogenesis. Here, we provide a roadmap to quantify the events that lead to morphogenesis from imaging to image analysis, quantification, and modeling, focusing on the discrete cellular and tissue shape changes, as well as their mechanical properties.

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.

Research on the locomotion of single cells on hard, flat surfaces brought insight into the mechanisms of leading-edge protrusion, spatially graded adhesion, front-rear coordination, and how intracellular and traction forces are harnessed to execute various maneuvers. Here, we highlight how, by studying a variety of cell types, shapes, and movements, Ken Jacobson and his collaborators made several discoveries that triggered the mechanistic understanding of cell motility. We then review the recent advancements and current perspectives in this field.

Mitochondria adapt to changing cellular environments, stress stimuli, and metabolic demands through dramatic morphological remodeling of their shape, and thus function. Such mitochondrial dynamics is often dependent on cytoskeletal filament interactions. However, the precise organization of these filamentous assemblies remains speculative. Here, we apply cryogenic electron tomography to directly image the nanoscale architecture of the cytoskeletal-membrane interactions involved in mitochondrial dynamics in response to damage. We induced mitochondrial damage via membrane depolarization, a cellular stress associated with mitochondrial fragmentation and mitophagy. We find that, in response to acute membrane depolarization, mammalian mitochondria predominantly organize into tubular morphology that abundantly displays constrictions. We observe long bundles of both unbranched actin and septin filaments enriched at these constrictions. We also observed septin-microtubule interactions at these sites and elsewhere, suggesting that these two filaments guide each other in the cytosolic space. Together, our results provide empirical parameters for the architecture of mitochondrial constriction factors to validate/refine existing models and inform the development of new ones.

Fetal lung fibroblasts contribute dynamic infrastructure for the developing lung. These cells undergo dynamic mechanical transitions, including cyclic stretch and spreading, which are integral to lung growth in utero. We investigated the role of the nuclear envelope protein emerin in cellular responses to these dynamic mechanical transitions. In contrast to control cells, which briskly realigned their nuclei, actin cytoskeleton, and extracellular matrices in response to cyclic stretch, fibroblasts that were acutely downregulated for emerin showed incomplete reorientation of both nuclei and actin cytoskeleton. Emerin-downregulated fibroblasts were also aberrantly circular in contrast to the spindle-shaped controls and exhibited an altered pattern of filamentous actin organization that was disconnected from the nucleus. Emerin knockdown was also associated with reduced myosin light chain phosphorylation during cell spreading. Interestingly, emerin-downregulated fibroblasts also demonstrated reduced fibronectin fibrillogenesis and production. These findings indicate that nuclear-cytoskeletal coupling serves a role in the dynamic regulation of cytoskeletal structure and function and may also impact the transmission of traction force to the extracellular matrix microenvironment.

Imaging two or more fluorescent biosensors in the same living cell can reveal the spatiotemporal coordination of protein activities. However, using multiple Förster resonance energy transfer (FRET) biosensors together is challenging due to toxicity and the need for orthogonal fluorophores. Here we generate a biosensor component that binds selectively to the activated conformation of three different proteins. This enabled multiplexed FRET with fewer fluorophores, and reduced toxicity. We generated this MultiBinder (MB) reagent for the GTPases RhoA, Rac1, and Cdc42 by combining portions of the downstream effector proteins Pak1 and Rhotekin. Using FRET between mCherry on the MB and YPet or mAmetrine on two target proteins, the activities of any pair of GTPases could be distinguished. The MB was used to image Rac1 and RhoA together with a third, dye-based biosensor for Cdc42. Quantifying effects of biosensor combinations on the frequency, duration, and velocity of cell protrusions and retractions demonstrated reduced toxicity. Multiplexed imaging revealed signaling hierarchies between the three proteins at the cell edge where they regulate motility.

To generate forces that drive migration of a eukaryotic cell, arrays of actin filaments (F-actin) are assembled at the cell's leading membrane edge. To maintain cell propulsion and respond to dynamic external cues, actin filaments must be disassembled to regenerate the actin monomers (G-actin), and transport of G-actin from sites of disassembly back to the leading edge completes the treadmilling cycle and limits the flux of F-actin assembly. Whether or not molecular diffusion is sufficient for G-actin transport has been a long-standing topic of debate, in part because the dynamic nature of cell motility and migration hinders the estimation of transport parameters. In this work, we applied an experimental system in which cells adopt an approximately constant and symmetrical shape; they cannot migrate but exhibit fast, steady treadmilling in the thin region protruding from the cell. Using fluorescence recovery after photobleaching, we quantified the relative concentrations and corresponding fluxes of F- and G-actin in this system. In conjunction with mathematical modeling, constrained by measured features of each region of interest, this approach revealed that diffusion alone cannot account for the transport of G-actin to the leading edge. Although G-actin diffusion and vectorial transport might vary with position in the protruding region, good agreement with the fluorescence recovery after photobleaching measurements was achieved by a model with constant G-actin diffusivity ∼2 μm2/s and anterograde G-actin velocity less than 1 μm/s.

Dectin-1A is a C-type lectin innate immunoreceptor that recognizes β-(1,3;1,6)-glucan, a structural component of Candida species cell walls. β-Glucans can adopt solution structures ranging from random coil to insoluble fiber due to tertiary (helical) and quaternary structure. Fungal β-glucans of medium and high molecular weight are highly structured, but low molecular weight glucan is much less structured. Despite similar affinity for Dectin-1, the ability of glucans to induce Dectin-1A-mediated signaling correlates with degree of structure. Glucan denaturation experiments showed that glucan structure determines agonistic potential, but not receptor binding affinity. We explored the impact of glucan structure on molecular aggregation of Dectin-1A. Stimulation with glucan signaling decreased Dectin-1A diffusion coefficient. Fluorescence measurements provided direct evidence of ligation-induced Dectin-1A aggregation, which positively correlated with increasing glucan structure content. In contrast, Dectin-1A is predominantly in a low aggregation state in resting cells. Molecular aggregates formed during interaction with highly structured, agonistic glucans did not exceed relatively small (<15 nm) clusters of a few engaged receptors. Finally, we observed increased molecular aggregation of Dectin-1A at fungal particle contact sites in a manner that positively correlated with the degree of exposed glucan on the particle surface. These results indicate that Dectin-1A senses the solution conformation of β-glucans through their varying ability to drive receptor dimer/oligomer formation and activation of membrane proximal signaling events.

Upon vascular injury, platelets form a hemostatic plug by binding to the subendothelium and to each other. Platelet-to-matrix binding is initially mediated by von Willebrand factor (VWF) and platelet-to-platelet binding is mediated mainly by fibrinogen and VWF. After binding, the actin cytoskeleton of a platelet drives its contraction, generating traction forces that are important to the cessation of bleeding. Our understanding of the relationship between adhesive environment, F-actin morphology, and traction forces is limited. Here, we examined F-actin morphology of platelets attached to surfaces coated with fibrinogen and VWF. We identified distinct F-actin patterns induced by these protein coatings and found that these patterns were identifiable into three classifications via machine learning: solid, nodular, and hollow. We observed that traction forces for platelets were significantly higher on VWF than on fibrinogen coatings and these forces varied by F-actin pattern. In addition, we analyzed the F-actin orientation in platelets and noted that their filaments were more circumferential when on fibrinogen coatings and having a hollow F-actin pattern, while they were more radial on VWF and having a solid F-actin pattern. Finally, we noted that subcellular localization of traction forces corresponded to protein coating and F-actin pattern: VWF-bound, solid platelets had higher forces at their central region while fibrinogen-bound, hollow platelets had higher forces at their periphery. These distinct F-actin patterns on fibrinogen and VWF and their differences in F-actin orientation, force magnitude, and force localization could have implications in hemostasis, thrombus architecture, and venous versus arterial thrombosis.

To facilitate rapid changes in morphology without endangering cell integrity, each cell possesses a substantial amount of cell surface excess (CSE) that can be promptly deployed to cover cell extensions. CSE can be stored in different types of small surface projections such as filopodia, microvilli, and ridges, with rounded bleb-like projections being the most common and rapidly achieved form of storage. We demonstrate that, similar to rounded cells in 2D culture, rounded cells in 3D collagen contain large amounts of CSE and use it to cover developing protrusions. Upon retraction of a protrusion, the CSE this produces is stored over the cell body similar to the CSE produced by cell rounding. We present high-resolution imaging of F-actin and microtubules (MTs) for different cell lines in a 3D environment and demonstrate the correlated changes between CSE and protrusion dynamics. To coordinate CSE storage and release with protrusion formation and motility, we expect cells to have specific mechanisms for regulating CSE, and we hypothesize that MTs play a substantial role in this mechanism by reducing cell surface dynamics and stabilizing CSE. We also suggest that different effects of MT depolymerization on cell motility, such as inhibiting mesenchymal motility and enhancing amoeboid, can be explained by this role of MTs in CSE regulation.

Membrane cholesterol-rich domains have been shown to be important for regulating a range of membrane protein activities. Low-density lipoprotein receptor (LDLR)-mediated internalization of cholesterol-rich LDL particles is tightly regulated by feedback mechanisms involving intracellular sterol sensors. Since LDLR plays a role in maintaining cellular cholesterol homeostasis, we explore the role that membrane domains may have in regulating LDLR activity. We expressed a fluorescent LDLR-mEGFP construct in HEK293T cells and imaged the unligated receptor or bound to an LDL/DiI fluorescent ligand using total internal reflection fluorescence microscopy. We studied the receptor's spatiotemporal dynamics using fluorescence fluctuation analysis methods. Image cross correlation spectroscopy reveals a lower LDL-to-LDLR binding fraction when membrane cholesterol concentrations are augmented using cholesterol esterase, and a higher binding fraction when the cells are treated with methyl-β-cyclodextrin) to lower membrane cholesterol. This suggests that LDLR's ability to metabolize LDL particles is negatively correlated to membrane cholesterol concentrations. We then tested if a change in activity is accompanied by a change in membrane localization. Image mean-square displacement analysis reveals that unligated LDLR-mEGFP and ligated LDLR-mEGFP/LDL-DiI constructs are transiently confined on the cell membrane, and the size of their confinement domains increases with augmented cholesterol concentrations. Receptor diffusion within the domains and their domain-escape probabilities decrease upon treatment with methyl-β-cyclodextrin, consistent with a change in receptor populations to more confined domains, likely clathrin-coated pits. We propose a feedback model to account for regulation of LDLR within the cell membrane: when membrane cholesterol concentrations are high, LDLR is sequestered in cholesterol-rich domains. These LDLR populations are attenuated in their efficacy to bind and internalize LDL. However, when membrane cholesterol levels drop, LDL has a higher binding affinity to its receptor and the LDLR transits to nascent clathrin-coated domains, where it diffuses at a slower rate while awaiting internalization.

Fluorescence redistribution after photobleaching is a commonly used method to understand the dynamic behavior of molecules within cells. Analytic solutions have been developed for specific, well-defined models of dynamic behavior in idealized geometries, but these solutions are inaccurate in complex geometries or when complex binding and diffusion behaviors exist. We demonstrate the use of numerical reaction-diffusion simulations using the Virtual Cell software platform to model fluorescence redistribution after photobleaching experiments. Multiple simulations employing parameter scans and varying bleaching locations and sizes can help to bracket diffusion coefficients and kinetic rate constants in complex image-based geometries. This approach is applied to problems in membrane surface diffusion as well as diffusion and binding in cytosolic volumes in complex cell geometries. In addition, we model diffusion and binding within phase-separated biomolecular condensates (liquid droplets). These are modeled as spherical low-affinity binding domains that also define a high viscosity medium for exchange of the free fluorescently labeled ligand with the external cytosol.

Early after-depolarizations (EADs) are action potential (AP) repolarization abnormalities that can trigger lethal arrhythmias. Simulations using biophysically-detailed cardiac myocyte models can reveal how model parameters influence the probability of these cellular arrhythmias, however such analyses can pose a huge computational burden. We have previously developed a highly simplified approach in which logistic regression models (LRMs) map parameters of complex cell models to the probability of ectopic beats (EBs). Here, we extend this approach to predict the probability of early after-depolarizations (P(EAD)) as a mechanistic metric of arrhythmic risk. We use the LRM to investigate how changes in parameters of the slow-activating delayed rectifier current (IKs) affect P(EAD) for 17 different Long QT syndrome type 1 (LQTS1) mutations. In this LQTS1 clinical arrhythmic risk prediction task, we compared P(EAD) for these 17 mutations with two other recently published model-based arrhythmia risk metrics (action potential morphology metric across populations of myocyte models and transmural repolarization prolongation based on a 1-D tissue-level model). These model-based risk metrics yield similar prediction performance however each fails to stratify clinical risk for a significant number of the 17 studied LQTS1 mutations. Nevertheless, an interpretable ensemble model using multivariate linear regression built by combining all of these model-based risk metrics successfully predicts the clinical risk of 17 mutations. These results illustrate the potential of computational approaches in arrhythmia risk prediction.

Piperine is the principal alkaloid present in black pepper and is well-known for its diverse pharmacological effects, including inhibition of different ion channels. Large conductance Ca2+-activated K+ channels (BK) are widely expressed across several tissues and play a vital role in many physiological functions. In this study, we have investigated the pharmacological effects of piperine on various BK channel subunit compositions (BKα, BKαβ1,4, BKαγ1,3) expressed in HEK293T cells. Piperine in zero Ca2+ reversibly inhibited currents from the pore-forming BKα channels in a dose-dependent manner with an IC50 of 4.8μM. Elevating the internal Ca2+ concentration from 0μM to 100μM significantly attenuated the inhibitory effects of piperine on BKα channels. The mutation G311S in the pore domain failed to alter the modulatory effects of piperine, whereas deletion of the entire cytoplasmic domain from BKα channels ablated its inhibitory effects. Addition of either BKβ1 or β4 regulatory subunits did not alter the efficacy of piperine on BKα channels. Interestingly, co-expression of either BKγ1 or BKγ3 subunits greatly diminished the ability of piperine to inhibit BKα channels. Our findings demonstrate that piperine is a potent natural modulator of BKα/BKαβ1,4 subunits but not ΒΚαγ1,3 subunits. The mechanism of piperine modulation appeared to be allosteric and differs from that of other BK pore blockers (paxilline, peptide toxins and quaternary ammonium compounds). Together, our results unravel the potential of piperine to inhibit BK channels, providing a new tool to explore mechanisms underlying the effects of regulatory subunits.

The β-secretase, BACE1, and the α-secretase, ADAM10, are known to competitively cleave amyloid precursor protein (APP) in the amyloid cascades of Alzheimer's disease. Cleavage of APP by BACE1 produces a 99-residue C-terminal peptide (APP-C99) that is subsequently cleaved by γ-secretase to form amyloid-β (Aβ) protein, whereas cleavage of APP by ADAM10 is nonamyloidogenic. It has been speculated that ADAM10/APP and BACE1/APP interactions are regulated by colocalization within and outside of liquid-ordered membrane domains; however, the mechanism of this regulation and the character of the proteins' transmembrane domains are not well understood. In this work, we have developed and characterized minimal congener sequences for the transmembrane domains of ADAM10 and BACE1 using a multiscale modeling approach combining both temperature replica exchange and conventional molecular dynamics simulations based on the coarse-grained Martini2.2 and all-atom CHARMM36 force fields. Our results show that membrane composition impacts the character of the transmembrane domains of BACE1 and ADAM10, adding credence to the speculation that membrane domains are involved in the etiology of Alzheimer's disease.