Human hnRNPA1 reorganizes telomere-bound Replication Protein A.

Human replication protein A (RPA) is a heterotrimeric ssDNA binding protein responsible for many aspects of cellular DNA metabolism. Dynamic interactions of the four RPA DNA binding domains (DBDs) with DNA control replacement of RPA by downstream proteins in various cellular metabolic pathways. RPA plays several important functions at telomeres where it binds to and melts telomeric G-quadruplexes, non-canonical DNA structures formed at the G-rich telomeric ssDNA overhangs. Here, we combine single-molecule total internal reflection fluorescence microscopy (smTIRFM) and mass photometry (MP) with biophysical and biochemical analyses to demonstrate that heterogeneous nuclear ribonucleoprotein A1 (hnRNPA1) specifically remodels RPA bound to telomeric ssDNA by dampening the RPA configurational dynamics and forming a ternary complex. Uniquely, among hnRNPA1 target RNAs, telomeric repeat-containing RNA (TERRA) is selectively capable of releasing hnRNPA1 from the RPA-telomeric DNA complex. We speculate that this telomere specific RPA-DNA-hnRNPA1 complex is an important structure in telomere protection. One Sentence Summary: At the single-stranded ends of human telomeres, the heterogeneous nuclear ribonucleoprotein A1 (hnRNPA1) binds to and modulates conformational dynamics of the ssDNA binding protein RPA forming a ternary complex which is controlled by telomeric repeat-containing RNA (TERRA).

A basic model for the association of ligands with membrane cholesterol: application to cytolysin binding.

Almost all the cholesterol in cellular membranes is associated with phospholipids in simple stoichiometric complexes. This limits the binding of sterol ligands such as filipin and perfringolysin O (PFO) to a small fraction of the total. We offer a simple mathematical model that characterizes this complexity. It posits that the cholesterol accessible to ligands has two forms: active cholesterol, which is that not complexed with phospholipids; and extractable cholesterol, that which ligands can capture competitively from the phospholipid complexes. Simulations based on the model match published data for the association of PFO oligomers with liposomes, plasma membranes, and the isolated endoplasmic reticulum. The model shows how the binding of a probe greatly underestimates cholesterol abundance when its affinity for the sterol is so weak that it competes poorly with the membrane phospholipids. Two examples are the understaining of plasma membranes by filipin and the failure of domain D4 of PFO to label their cytoplasmic leaflets. Conversely, the exaggerated staining of endolysosomes suggests that their cholesterol, being uncomplexed, is readily available. The model is also applicable to the association of cholesterol with intrinsic membrane proteins. For example, it supports the hypothesis that the sharp threshold in the regulation of homeostatic endoplasmic reticulum proteins by cholesterol derives from the cooperativity of their binding to the sterol weakly held by the phospholipids. Thus, the model explicates the complexity inherent in the binding of ligands like PFO and filipin to the small accessible fraction of membrane cholesterol.

A basic model for cell cholesterol homeostasis.

Cells manage their cholesterol by negative feedback using a battery of sterol-responsive proteins. How these activities are coordinated so as to specify the abundance and distribution of the sterol is unclear. We present a simple mathematical model that addresses this question. It assumes that almost all of the cholesterol is associated with phospholipids in stoichiometric complexes. A small fraction of the sterol is uncomplexed and thermodynamically active. It equilibrates among the organelles, setting their sterol level according to the affinity of their phospholipids. The activity of the homeostatic proteins in the cytoplasmic membranes is then set by their fractional saturation with uncomplexed cholesterol in competition with the phospholipids. The high-affinity phospholipids in the plasma membrane (PM) are filled to near stoichiometric equivalence, giving it most of the cell sterol. Notably, the affinity of the phospholipids in the endomembranes (EMs) is lower by orders of magnitude than that of the phospholipids in the PM. Thus, the small amount of sterol in the EMs rests far below stoichiometric capacity. Simulations match a variety of experimental data. The model captures the essence of cell cholesterol homeostasis, makes coherent a diverse set of experimental findings, provides a surprising prediction and suggests new experiments.

KERA: analysis tool for multi-process, multi-state single-molecule data.

Molecular machines within cells dynamically assemble, disassemble and reorganize. Molecular interactions between their components can be observed at the single-molecule level and quantified using colocalization single-molecule spectroscopy, in which individual labeled molecules are seen transiently associating with a surface-tethered partner, or other total internal reflection fluorescence microscopy approaches in which the interactions elicit changes in fluorescence in the labeled surface-tethered partner. When multiple interacting partners can form ternary, quaternary and higher order complexes, the types of spatial and temporal organization of these complexes can be deduced from the order of appearance and reorganization of the components. Time evolution of complex architectures can be followed by changes in the fluorescence behavior in multiple channels. Here, we describe the kinetic event resolving algorithm (KERA), a software tool for organizing and sorting the discretized fluorescent trajectories from a range of single-molecule experiments. KERA organizes the data in groups by transition patterns, and displays exhaustive dwell time data for each interaction sequence. Enumerating and quantifying sequences of molecular interactions provides important information regarding the underlying mechanism of the assembly, dynamics and architecture of the macromolecular complexes. We demonstrate KERA's utility by analyzing conformational dynamics of two DNA binding proteins: replication protein A and xeroderma pigmentosum complementation group D helicase.

Distinguishing synaptic vesicle precursor navigation of microtubule ends with a single rate constant model.

Axonal motor driven cargo utilizes the microtubule cytoskeleton in order to direct cargo, such as synaptic vesicle precursors (SVP), to where they are needed. This transport requires vesicles to travel up to microns in distance. It has recently been observed that finite microtubule lengths can act as roadblocks inhibiting SVP and increasing the time required for transport. SVPs reach the end of a microtubule and pause until they can navigate to a neighboring microtubule in order to continue transport. The mechanism(s) by which axonal SVPs navigate the end of a microtubule in order to continue mobility is unknown. In this manuscript we model experimentally observed vesicle pausing at microtubule ends in C. elegans. We show that a single rate-constant model reproduces the time SVPs pause at MT-ends. This model is based on the time an SVP must detach from its current microtubule and re-attach to a neighboring microtubule. We show that vesicle pause times are different for anterograde and retrograde motion, suggesting that vesicles utilize different proteins at plus and minus end sites. Last, we show that vesicles do not likely utilize a tug-of-war like mechanism and reverse direction in order to navigate microtubule ends.

Effects of Intercalating Molecules on the Polymer Properties of DNA.

Atomic force microscopy (AFM) enables determination of physical properties from single DNA molecules. Insertion of aromatic molecules into the structure of DNA results in morphological changes. However, the accompanying changes to elastic properties due to this insertion are not fully understood. AFM was used to examine the morphological effects of intercalator binding and report changes in the elastic properties of intrinsically straight DNA molecules. The persistence length and polymer extension were characterized in the presence of three intercalating molecules: ethidium bromide and the less well studied chloroquine and acridine. It was found that all three intercalators significantly increased the bending persistence length. In addition, an analysis of the normal bending modes of the static molecules corroborated these results. This approach of measuring binding effects of intercalators on DNA physical properties using a model system of intrinsically straight DNA is applicable to other DNA binding ligands and other modes of DNA interaction.

Dynamic motions of molecular motors in the actin cytoskeleton.

During intracellular transport, cellular cargos, such as organelles, vesicles, and proteins, are transported within cells. Intracellular transport plays an important role in diverse cellular functions. Molecular motors walking on the cytoskeleton facilitate active intracellular transport, which is more efficient than diffusion-based passive transport. Active transport driven by kinesin and dynein walking on microtubules has been studied well during recent decades. However, mechanisms of active transport occurring in disorganized actin networks via myosin motors remain elusive. To provide physiologically relevant insights, we probed motions of myosin motors in actin networks under various conditions using our well-established computational model that rigorously accounts for the mechanical and dynamical behaviors of the actin cytoskeleton. We demonstrated that myosin motions can be confined due to three different reasons in the absence of F-actin turnover. We verified mechanisms of motor stalling using in vitro reconstituted actomyosin networks. We also found that with F-actin turnover, motors consistently move for a long time without significant confinement. Our study sheds light on the importance of F-actin turnover for effective active transport in the actin cytoskeleton.

Anomalous intracellular transport phases depend on cytoskeletal network features.

Intracellular transport in eukaryotic cells consists of phases of passive, diffusion-based transport and active, motor-driven transport along filaments that make up the cell's cytoskeleton. The interplay between superdiffusive transport along cytoskeletal filaments and the anomalous nature of subdiffusion in the bulk can lead to novel effects in transport behavior at the cellular scale. Here we develop a computational model of the process with cargo being ballistically transported along explicitly modeled cytoskeletal filament networks and passively transported in the cytoplasm by a subdiffusive continuous-time random walk (CTRW). We show that, over a physiologically relevant range of filament lengths and numbers, the network introduces a filament-length sensitive superdiffusive phase at early times which crosses over to a phase where the CTRW is dominant and produces subdiffusion at late times. We apply our approach to the problem of insulin secretion from cells and show that the superdiffusive phase introduced by the filament network manifests as a peak in the secretion at early times followed by an extended sustained release phase that is dominated by the CTRW process at late times. Our results are consistent with in vivo observations of insulin transport in healthy cells and shed light on the potential for the cell to tune functionally important transport phases by altering its cytoskeletal network.

Dynamics and selective remodeling of the DNA-binding domains of RPA.

Replication protein A (RPA) coordinates important DNA metabolic events by stabilizing single-stranded DNA (ssDNA) intermediates, activating the DNA-damage response and handing off ssDNA to the appropriate downstream players. Six DNA-binding domains (DBDs) in RPA promote high-affinity binding to ssDNA yet also allow RPA displacement by lower affinity proteins. We generated fluorescent versions of Saccharomyces cerevisiae RPA and visualized the conformational dynamics of individual DBDs in the context of the full-length protein. We show that both DBD-A and DBD-D rapidly bind to and dissociate from ssDNA while RPA remains bound to ssDNA. The recombination mediator protein Rad52 selectively modulates the dynamics of DBD-D. These findings reveal how RPA-interacting proteins with lower ssDNA binding affinities can access the occluded ssDNA and remodel individual DBDs to replace RPA.

Entropy production rate is maximized in non-contractile actomyosin.

The actin cytoskeleton is an active semi-flexible polymer network whose non-equilibrium properties coordinate both stable and contractile behaviors to maintain or change cell shape. While myosin motors drive the actin cytoskeleton out-of-equilibrium, the role of myosin-driven active stresses in the accumulation and dissipation of mechanical energy is unclear. To investigate this, we synthesize an actomyosin material in vitro whose active stress content can tune the network from stable to contractile. Each increment in activity determines a characteristic spectrum of actin filament fluctuations which is used to calculate the total mechanical work and the production of entropy in the material. We find that the balance of work and entropy does not increase monotonically and the entropy production rate is maximized in the non-contractile, stable state of actomyosin. Our study provides evidence that the origins of entropy production and activity-dependent dissipation relate to disorder in the molecular interactions between actin and myosin.

Improved statistical methods enable greater sensitivity in rhythm detection for genome-wide data.

Robust methods for identifying patterns of expression in genome-wide data are important for generating hypotheses regarding gene function. To this end, several analytic methods have been developed for detecting periodic patterns. We improve one such method, JTK_CYCLE, by explicitly calculating the null distribution such that it accounts for multiple hypothesis testing and by including non-sinusoidal reference waveforms. We term this method empirical JTK_CYCLE with asymmetry search, and we compare its performance to JTK_CYCLE with Bonferroni and Benjamini-Hochberg multiple hypothesis testing correction, as well as to five other methods: cyclohedron test, address reduction, stable persistence, ANOVA, and F24. We find that ANOVA, F24, and JTK_CYCLE consistently outperform the other three methods when data are limited and noisy; empirical JTK_CYCLE with asymmetry search gives the greatest sensitivity while controlling for the false discovery rate. Our analysis also provides insight into experimental design and we find that, for a fixed number of samples, better sensitivity and specificity are achieved with higher numbers of replicates than with higher sampling density. Application of the methods to detecting circadian rhythms in a metadataset of microarrays that quantify time-dependent gene expression in whole heads of Drosophila melanogaster reveals annotations that are enriched among genes with highly asymmetric waveforms. These include a wide range of oxidation reduction and metabolic genes, as well as genes with transcripts that have multiple splice forms.

Distribution of directional change as a signature of complex dynamics.

Analyses of random walks traditionally use the mean square displacement (MSD) as an order parameter characterizing dynamics. We show that the distribution of relative angles of motion between successive time intervals of random walks in two or more dimensions provides information about stochastic processes beyond the MSD. We illustrate the behavior of this measure for common models and apply it to experimental particle tracking data. For a colloidal system, the distribution of relative angles reports sensitively on caging as the density varies. For transport mediated by molecular motors on filament networks in vitro and in vivo, we discover self-similar properties that cannot be described by existing models and discuss possible scenarios that can lead to the elucidated statistical features.

Stability and stoichiometry of bilayer phospholipid-cholesterol complexes: relationship to cellular sterol distribution and homeostasis.

Is cholesterol distributed among intracellular compartments by passive equilibration down its chemical gradient? If so, its distribution should reflect the relative cholesterol affinity of the constituent membrane phospholipids as well as their capacity for association with the sterol. We examined this issue by analyzing the reactivity to cholesterol oxidase of large unilamellar vesicles (LUVs) containing phospholipids and varied levels of cholesterol. The rates of cholesterol oxidation differed among the various phospholipid environments by roughly 4 orders of magnitude. Furthermore, accessibility to the enzyme increased by orders of magnitude at cholesterol thresholds that suggested cholesterol:phospholipid association ratios of 1:1, 2:3, or 1:2 (moles:moles). The accessibility of cholesterol above these thresholds was still constrained by its particular phospholipid environment. One phospholipid, 1-stearoyl-2-oleoyl-sn-glycero-3-phosphatidylserine, exhibited no threshold. The analysis suggested values for the stoichiometries of the putative cholesterol-phospholipid complexes, their relative stabilities, and the fractions of bilayer cholesterol not in complexes at the threshold equivalence points. Predictably, the saturated phosphorylcholine species had the lowest apparent stoichiometric ratios and the strongest associations with cholesterol. These results are in general agreement with the equilibrium distribution of cholesterol between the various LUVs and methyl-ß-cyclodextrin. In addition, the behavior of the cholesterol in intact human red blood cells matched predictions made from LUVs of the corresponding composition. These results support a passive mechanism for the intracellular distribution of cholesterol that can provide a signal for its homeostatic regulation.

Intracellular transport of insulin granules is a subordinated random walk.

We quantitatively analyzed particle tracking data on insulin granules expressing fluorescent fusion proteins in MIN6 cells to better understand the motions contributing to intracellular transport and, more generally, the means for characterizing systems far from equilibrium. Care was taken to ensure that the statistics reflected intrinsic features of the individual granules rather than details of the measurement and overall cell state. We find anomalous diffusion. Interpreting such data conventionally requires assuming that a process is either ergodic with particles working against fluctuating obstacles (fractional brownian motion) or nonergodic with a broad distribution of dwell times for traps (continuous-time random walk). However, we find that statistical tests based on these two models give conflicting results. We resolve this issue by introducing a subordinated scheme in which particles in cages with random dwell times undergo correlated motions owing to interactions with a fluctuating environment. We relate this picture to the underlying microtubule structure by imaging in the presence of vinblastine. Our results provide a simple physical picture for how diverse pools of insulin granules and, in turn, biphasic secretion could arise.

Short-time evolution in the adaptive immune system.

We exploit a simple model to numerically and analytically investigate the effect of enforcing a time constraint for achieving a system-wide goal during an evolutionary dynamics. This situation is relevant to finding antibody specificities in the adaptive immune response as well as to artificial situations in which an evolutionary dynamics is used to generate a desired capability in a limited number of generations. When the likelihood of finding the target phenotype is low, we find that the optimal mutation rate can exceed the error threshold, in contrast to conventional evolutionary dynamics. We also show how a logarithmic correction to the usual inverse scaling of population size with mutation rate arises. Implications for natural and artificial evolutionary situations are discussed.

Entrainment of a driven oscillator as a dynamical phase transition.

Large deviation theory has emerged as the natural language with which to study transitions between dynamical phases, such as the glass transition. Here we use this approach to show that the entrainment of an oscillator to an external driving force can be described as a phase transition between two regions in a joint space-time representation. Specifically, we numerically obtain exact solutions of the large deviation function for a discrete, finite model of an oscillator under a periodically varying external force. We find that the first derivative of the expectation value of the current diverges in the limit of large system size. For weak forcing, perturbation theory allows us to relate the observed frequency of the oscillator to the spectrum of eigenvalues for the unforced system. We find that in the entrainment region where the oscillator exhibits the frequency of the driving force, there is a strong coupling between the ground state and excited states.

Modeling bistable cell-fate choices in the Drosophila eye: qualitative and quantitative perspectives.

A major goal of developmental biology is to understand the molecular mechanisms whereby genetic signaling networks establish and maintain distinct cell types within multicellular organisms. Here, we review cell-fate decisions in the developing eye of Drosophila melanogaster and the experimental results that have revealed the topology of the underlying signaling circuitries. We then propose that switch-like network motifs based on positive feedback play a central role in cell-fate choice, and discuss how mathematical modeling can be used to understand and predict the bistable or multistable behavior of such networks.