Robustness of DNA looping across multiple cell divisions in individual bacteria.

DNA looping has emerged as a central paradigm of transcriptional regulation, as it is shared across many living systems. One core property of DNA looping-based regulation is its ability to greatly enhance repression or activation of genes with only a few copies of transcriptional regulators. However, this property based on a small number of proteins raises the question of the robustness of such a mechanism with respect to the large intracellular perturbations taking place during growth and division of the cell. Here we address the issue of sensitivity to variations of intracellular parameters of gene regulation by DNA looping. We use the lac system as a prototype to experimentally identify the key features of the robustness of DNA looping in growing Escherichia coli cells. Surprisingly, we observe time intervals of tight repression spanning across division events, which can sometimes exceed 10 generations. Remarkably, the distribution of such long time intervals exhibits memoryless statistics that is mostly insensitive to repressor concentration, cell division events, and the number of distinct loops accessible to the system. By contrast, gene regulation becomes highly sensitive to these perturbations when DNA looping is absent. Using stochastic simulations, we propose that the observed robustness to division emerges from the competition between fast, multiple rebinding events of repressors and slow initiation rate of the RNA polymerase. We argue that fast rebinding events are a direct consequence of DNA looping that ensures robust gene repression across a range of intracellular perturbations.

All-or-none amyloid disassembly via chaperone-triggered fibril unzipping favors clearance of α-synuclein toxic species.

α-synuclein aggregation is present in Parkinson's disease and other neuropathologies. Among the assemblies that populate the amyloid formation process, oligomers and short fibrils are the most cytotoxic. The human Hsc70-based disaggregase system can resolve α-synuclein fibrils, but its ability to target other toxic assemblies has not been studied. Here, we show that this chaperone system preferentially disaggregates toxic oligomers and short fibrils, while its activity against large, less toxic amyloids is severely impaired. Biochemical and kinetic characterization of the disassembly process reveals that this behavior is the result of an all-or-none abrupt solubilization of individual aggregates. High-speed atomic force microscopy explicitly shows that disassembly starts with the destabilization of the tips and rapidly progresses to completion through protofilament unzipping and depolymerization without accumulation of harmful oligomeric intermediates. Our data provide molecular insights into the selective processing of toxic amyloids, which is critical to identify potential therapeutic targets against increasingly prevalent neurodegenerative disorders.

Clearly Detectable, Kinetically Restricted Solid-Solid Phase Transition in cis-Ceramide Monolayers.

Sphingosine [(2 S,3 R,4 E)-2-amino-4-octadecene-1,3-diol] is the most common sphingoid base in mammals. Ceramides are N-acyl sphingosines. Numerous small variations on this canonical structure are known, including the 1-deoxy, the 4,5-dihydro, and many others. However, whenever there is a Δ4 double bond, it adopts the trans (or E) configuration. We synthesized a ceramide containing 4 Z-sphingosine and palmitic acid ( cis-pCer) and studied its behavior in the form of monolayers extended on an air-water interface. cis-pCer acted very differently from the trans isomer in that, upon lateral compression of the monolayer, a solid-solid transition was clearly observed at a mean molecular area ≤44 Å2·molecule-1, whose characteristics depended on the rate of compression. The solid-solid transition, as well as states of domain coexistence, could be imaged by atomic force microscopy and by Brewster-angle microscopy. Atomistic molecular dynamics simulations provided results compatible with the experimentally observed differences between the cis and trans isomers. The data can help in the exploration of other solid-solid transitions in lipids, both in vitro and in vivo, that have gone up to now undetected because of their less obvious change in surface properties along the transition, as compared to cis-pCer.

Literature-based automated reconstruction, expansion, and refinement of the TGF-ß superfamily ligand-receptor network.

The TGF-ß pathway transduces a variety of extracellular signals into intracellular responses that control multiple cellular processes, including cell growth, apoptosis, and differentiation. It encompasses 33 ligands that interact with 7 type II receptors and 5 type I receptors at the plasma membrane to potentially form 1,155 ligand-receptor complexes in mammalian cells. Retrieving the information of the complexes that are actually formed from reading the literature might be tedious and prone to missing links. Here, we have developed an automated literature-mining procedure to obtain the interactions of the TGF-ß ligand-receptor network. By querying the Information Hyperlinked over Proteins (iHOP) online service and processing the results, we were able to find pairwise interactions between ligands and receptors that allowed us to build the network automatically from the literature. Comparison with available published review papers indicates that this method is able to automatically reconstruct and expand the TGF-ß superfamily ligand-receptor network. Retrieving and parsing the full text of the manuscripts containing the interactions allowed us to refine the network interactions for specific cell lines.

The unreasonable effectiveness of equilibrium gene regulation through the cell cycle.

Systems like the prototypical lac operon can reliably hold repression of transcription upon DNA replication across cell cycles with just 10 repressor molecules per cell and behave as if they were at equilibrium. The origin of this phenomenology is still an unresolved question. Here, we develop a general theory to analyze strong perturbations in quasi-equilibrium systems and use it to quantify the effects of DNA replication in gene regulation. We find a scaling law linking actual with predicted equilibrium transcription via a single kinetic parameter. We show that even the lac operon functions beyond the physical limits of naive regulation through compensatory mechanisms that suppress non-equilibrium effects. Synthetic systems without adjuvant activators, such as the cAMP receptor protein (CRP), lack this reliability. Our results provide a rationale for the function of CRP, beyond just being a tunable activator, as a mitigator of cell cycle perturbations.

Nanocrystal Assemblies: Current Advances and Open Problems.

We explore the potential of nanocrystals (a term used equivalently to nanoparticles) as building blocks for nanomaterials, and the current advances and open challenges for fundamental science developments and applications. Nanocrystal assemblies are inherently multiscale, and the generation of revolutionary material properties requires a precise understanding of the relationship between structure and function, the former being determined by classical effects and the latter often by quantum effects. With an emphasis on theory and computation, we discuss challenges that hamper current assembly strategies and to what extent nanocrystal assemblies represent thermodynamic equilibrium or kinetically trapped metastable states. We also examine dynamic effects and optimization of assembly protocols. Finally, we discuss promising material functions and examples of their realization with nanocrystal assemblies.

Systems biophysics of gene expression.

Gene expression is a process central to any form of life. It involves multiple temporal and functional scales that extend from specific protein-DNA interactions to the coordinated regulation of multiple genes in response to intracellular and extracellular changes. This diversity in scales poses fundamental challenges to the use of traditional approaches to fully understand even the simplest gene expression systems. Recent advances in computational systems biophysics have provided promising avenues to reliably integrate the molecular detail of biophysical process into the system behavior. Here, we review recent advances in the description of gene regulation as a system of biophysical processes that extend from specific protein-DNA interactions to the combinatorial assembly of nucleoprotein complexes. There is now basic mechanistic understanding on how promoters controlled by multiple, local and distal, DNA binding sites for transcription factors can actively control transcriptional noise, cell-to-cell variability, and other properties of gene regulation, including precision and flexibility of the transcriptional responses.

Control of gene expression by modulated self-assembly.

Numerous transcription factors self-assemble into different order oligomeric species in a way that is actively regulated by the cell. Until now, no general functional role has been identified for this widespread process. Here, we capture the effects of modulated self-assembly in gene expression with a novel quantitative framework. We show that this mechanism provides precision and flexibility, two seemingly antagonistic properties, to the sensing of diverse cellular signals by systems that share common elements present in transcription factors like p53, NF-κB, STATs, Oct and RXR. Applied to the nuclear hormone receptor RXR, this framework accurately reproduces a broad range of classical, previously unexplained, sets of gene expression data and corroborates the existence of a precise functional regime with flexible properties that can be controlled both at a genome-wide scale and at the individual promoter level.

Dynamics-informed deconvolutional neural networks for super-resolution identification of regime changes in epidemiological time series.

The ability to infer the timing and amplitude of perturbations in epidemiological systems from their stochastically spread low-resolution outcomes is crucial for multiple applications. However, the general problem of connecting epidemiological curves with the underlying incidence lacks the highly effective methodology present in other inverse problems, such as super-resolution and dehazing from computer vision. Here, we develop an unsupervised physics-informed convolutional neural network approach in reverse to connect death records with incidence that allows the identification of regime changes at single-day resolution. Applied to COVID-19 data with proper regularization and model-selection criteria, the approach can identify the implementation and removal of lockdowns and other nonpharmaceutical interventions (NPIs) with 0.93-day accuracy over the time span of a year.

Multi-landmark alignment of genomic signals reveals conserved expression patterns across transcription start sites.

The prevalent one-dimensional alignment of genomic signals to a reference landmark is a cornerstone of current methods to study transcription and its DNA-dependent processes but it is prone to mask potential relations among multiple DNA elements. We developed a systematic approach to align genomic signals to multiple locations simultaneously by expanding the dimensionality of the genomic-coordinate space. We analyzed transcription in human and uncovered a complex dependence on the relative position of neighboring transcription start sites (TSSs) that is consistently conserved among cell types. The dependence ranges from enhancement to suppression of transcription depending on the relative distances to the TSSs, their intragenic position, and the transcriptional activity of the gene. Our results reveal a conserved hierarchy of alternative TSS usage within a previously unrecognized level of genomic organization and provide a general methodology to analyze complex functional relationships among multiple types of DNA elements.

Trafficking coordinate description of intracellular transport control of signaling networks.

Many cellular networks rely on the regulated transport of their components to transduce extracellular information into precise intracellular signals. The dynamics of these networks is typically described in terms of compartmentalized chemical reactions. There are many important situations, however, in which the properties of the compartments change continuously in a way that cannot naturally be described by chemical reactions. Here, we develop an approach based on transport along a trafficking coordinate to precisely describe these processes and we apply it explicitly to the TGF-ß signal transduction network, which plays a fundamental role in many diseases and cellular processes. The results of this newly introduced approach accurately capture the distinct TGF-ß signaling dynamics of cells with and without cancerous backgrounds and provide an avenue to predict the effects of chemical perturbations in a way that closely recapitulates the observed cellular behavior.

CplexA: a Mathematica package to study macromolecular-assembly control of gene expression.

SUMMARY: Macromolecular assembly coordinates essential cellular processes, such as gene regulation and signal transduction. A major challenge for conventional computational methods to study these processes is tackling the exponential increase of the number of configurational states with the number of components. CplexA is a Mathematica package that uses functional programming to efficiently compute probabilities and average properties over such exponentially large number of states from the energetics of the interactions. The package is particularly suited to study gene expression at complex promoters controlled by multiple, local and distal, DNA binding sites for transcription factors. AVAILABILITY: CplexA is freely available together with documentation at http://sourceforge.net/projects/cplexa/.

Computational design of single-stranded DNA hairpin aptamers immobilized on a biosensor substrate.

Aptamer interactions with a surface of attachment are central to the design and performance of aptamer-based biosensors. We have developed a computational modeling approach to study different system designs-including different aptamer-attachment ends, aptamer surface densities, aptamer orientations, and solvent solutions-and applied it to an anti MUC1 aptamer tethered to a silica biosensor substrate. Amongst all the system designs explored, we found that attaching the anti MUC1 aptamer through the 5' terminal end, in a high surface density configuration, and solvated in a 0.8 M NaCl solution provided the best exposure of the aptamer MUC1 binding regions and resulted in the least amount of aptamer backbone fluctuations. Many of the other designs led to non-functional systems, with the aptamer collapsing onto the surface. The computational approach we have developed and the resulting analysis techniques can be employed for the rational design of aptamer-based biosensors and provide a valuable tool for improving biosensor performance and repeatability.

Ascertaining the initiation of epidemic resurgences: an application to the COVID-19 second surges in Europe and the Northeast United States.

Assessing a potential resurgence of an epidemic outbreak with certainty is as important as it is challenging. The low number of infectious individuals after a long regression, and the randomness associated with it, makes it difficult to ascertain whether the infectious population is growing or just fluctuating. We have developed an approach to compute confidence intervals for the switching time from decay to growth and to compute the corresponding multiple-location aggregated quantities over a region to increase the precision of the determination. We estimated the aggregate prevalence over time for Europe and the northeast United States to characterize the COVID-19 second surge in these regions during year 2020. We find a starting date as early as 3 July (95% confidence interval (CI): 1-6 July) for Europe and 19 August (95% CI: 16-23 August) for the northeast United States; subsequent infectious populations that, as of 31 December, have always increased or remained stagnant; and the resurgences being the collective effect of each overall region with no location, either country or state, dominating the regional dynamics by itself.

Reliably quantifying the evolving worldwide dynamic state of the COVID-19 outbreak from death records, clinical parametrization, and demographic data.

The dynamic characterization of the COVID-19 outbreak is critical to implement effective actions for its control and eradication but the information available at a global scale is not sufficiently reliable to be used directly. Here, we develop a quantitative approach to reliably quantify its temporal evolution and controllability through the integration of multiple data sources, including death records, clinical parametrization of the disease, and demographic data, and we explicitly apply it to countries worldwide, covering 97.4% of the human population, and to states within the United States (US). The validation of the approach shows that it can accurately reproduce the available prevalence data and that it can precisely infer the timing of nonpharmaceutical interventions. The results of the analysis identified general patterns of recession, stabilization, and resurgence. The diversity of dynamic behaviors of the outbreak across countries is paralleled by those of states and territories in the US, converging to remarkably similar global states in both cases. Our results offer precise insights into the dynamics of the outbreak and an efficient avenue for the estimation of the prevalence rates over time.

Ab initio thermodynamic modeling of distal multisite transcription regulation.

Transcription regulation typically involves the binding of proteins over long distances on multiple DNA sites that are brought close to each other by the formation of DNA loops. The inherent complexity of assembling regulatory complexes on looped DNA challenges the understanding of even the simplest genetic systems, including the prototypical lac operon. Here we implement a scalable approach based on thermodynamic molecular properties to model ab initio systems regulated through multiple DNA sites with looping. We show that this approach applied to the lac operon accurately predicts the system behavior for a wide range of cellular conditions, which include the transcription rate over five orders of magnitude as a function of the repressor concentration for wild type and all seven combinations of deletions of three operators, as well as the observed induction curves for cells with and without active catabolite activator protein. Our results provide new insights into the detailed functioning of the lac operon and reveal an efficient avenue to incorporate the required underlying molecular complexity into fully predictive models of gene regulation.

DNA looping in gene regulation: from the assembly of macromolecular complexes to the control of transcriptional noise.

The formation of DNA loops by the binding of proteins and protein complexes at distal DNA sites plays a central role in many cellular processes, such as transcription, recombination and replication. Important thermodynamic concepts underlie the assembly of macromolecular complexes on looped DNA. The effects that this process has on the properties of gene regulation extend beyond the traditional view of DNA looping as a mechanism to increase the affinity of regulatory molecules for their cognate sites. Recent developments indicate that DNA looping can also lead to the suppression of cell-to-cell variability, the control of transcriptional noise, and the activation of cooperative interactions on demand.

DNA looping: the consequences and its control.

The formation of DNA loops by proteins and protein complexes is ubiquitous to many fundamental cellular processes, including transcription, recombination and replication. Recently, advances have been made in understanding the properties of DNA looping in its natural context and how they propagate to cellular behavior through gene regulation. The result of connecting the molecular properties of DNA looping with cellular physiology measurements indicates that looping of DNA in vivo is much more complex and easier than predicted from current models, and reveals a wealth of previously unappreciated details.

Stochastic dynamics of macromolecular-assembly networks.

The formation and regulation of macromolecular complexes provides the backbone of most cellular processes, including gene regulation and signal transduction. The inherent complexity of assembling macromolecular structures makes current computational methods strongly limited for understanding how the physical interactions between cellular components give rise to systemic properties of cells. Here, we present a stochastic approach to study the dynamics of networks formed by macromolecular complexes in terms of the molecular interactions of their components. Exploiting key thermodynamic concepts, this approach makes it possible to both estimate reaction rates and incorporate the resulting assembly dynamics into the stochastic kinetics of cellular networks. As prototype systems, we consider the lac operon and phage lambda induction switches, which rely on the formation of DNA loops by proteins and on the integration of these protein-DNA complexes into intracellular networks. This cross-scale approach offers an effective starting point to move forward from network diagrams, such as those of protein-protein and DNA-protein interaction networks, to the actual dynamics of cellular processes.