Raúl Guantes

Dynamical Principles of Two-Component Genetic Oscillators

Genetic oscillators based on the interaction of a small set of molecular components have been shown to be involved in the regulation of the cell cycle, the circadian rhythms, or the response of several signaling pathways. Uncovering the functional properties of such oscillators then becomes important for the understanding of these cellular processes and for the characterization of fundamental properties of more complex clocks. Here, we show how the dynamics of a minimal two-component oscillator is drastically affected by its genetic implementation. We consider a repressor and activator element combined in a simple logical motif. While activation is always exerted at the transcriptional level, repression is alternatively operating at the transcriptional (Design I) or post-translational (Design II) level. These designs display differences on basic oscillatory features and on their behavior with respect to molecular noise or entrainment by periodic signals. In particular, Design I induces oscillations with large activator amplitudes and arbitrarily small frequencies, and acts as an “integrator” of external stimuli, while Design II shows emergence of oscillations with finite, and less variable, frequencies and smaller amplitudes, and detects better frequency-encoded signals (“resonator”). Similar types of stimulus response are observed in neurons, and thus this work enables us to connect very different biological contexts. These dynamical principles are relevant for the characterization of the physiological roles of simple oscillator motifs, the understanding of core machineries of complex clocks, and the bio-engineering of synthetic oscillatory circuits.

Global variability in gene expression and alternative splicing is modulated by mitochondrial content

Noise in gene expression is a main determinant of phenotypic variability. Increasing experimental evidence suggests that genome-wide cellular constraints largely contribute to the heterogeneity observed in gene products. It is still unclear, however, which global factors affect gene expression noise and to what extent. Since eukaryotic gene expression is an energy demanding process, differences in the energy budget of each cell could determine gene expression differences. Here, we quantify the contribution of mitochondrial variability (a natural source of ATP variation) to global variability in gene expression. We find that changes in mitochondrial content can account for ∼50% of the variability observed in protein levels. This is the combined result of the effect of mitochondria dosage on transcription and translation apparatus content and activities. Moreover, we find that mitochondrial levels have a large impact on alternative splicing, thus modulating both the abundance and type of mRNAs. A simple mathematical model in which mitochondrial content simultaneously affects transcription rate and splicing site choice can explain the alternative splicing data. The results of this study show that mitochondrial content (and/or probably function) influences mRNA abundance, translation, and alternative splicing, which ultimately affects cellular phenotype.

Sodium pumps adapt spike bursting to stimulus statistics

Pump activity is a homeostatic mechanism that maintains ionic gradients. Here we examined whether the slow reduction in excitability induced by sodium-pump activity that has been seen in many neuronal types is also involved in neuronal coding. We took intracellular recordings from a spike-bursting sensory neuron in the leech Hirudo medicinalis in response to naturalistic tactile stimuli with different statistical distributions. We show that regulation of excitability by sodium pumps is necessary for the neuron to make different responses depending on the statistical context of the stimuli. In particular, sodium-pump activity allowed spike-burst sizes and rates to code not for stimulus values per se, but for their ratio with the standard deviation of the stimulus distribution. Modeling further showed that sodium pumps can be a general mechanism of adaptation to statistics on the time scale of 1 min. These results implicate the ubiquitous pump activity in the adaptation of neural codes to statistics.

A Tunable Coarse-Grained Model for Ligand-Receptor Interaction

Cell-surface receptors are the most common target for therapeutic drugs. The design and optimization of next generation synthetic drugs require a detailed understanding of the interaction with their corresponding receptors. Mathematical approximations to study ligand-receptor systems based on reaction kinetics strongly simplify the spatial constraints of the interaction, while full atomistic ligand-receptor models do not allow for a statistical many-particle analysis, due to their high computational requirements. Here we present a generic coarse-grained model for ligand-receptor systems that accounts for the essential spatial characteristics of the interaction, while allowing statistical analysis. The model captures the main features of ligand-receptor kinetics, such as diffusion dependence of affinity and dissociation rates. Our model is used to characterize chimeric compounds, designed to take advantage of the receptor over-expression phenotype of certain diseases to selectively target unhealthy cells. Molecular dynamics simulations of chimeric ligands are used to study how selectivity can be optimized based on receptor abundance, ligand-receptor affinity and length of the linker between both ligand subunits. Overall, this coarse-grained model is a useful approximation in the study of systems with complex ligand-receptor interactions or spatial constraints.

A mathematical model for the rational design of chimeric ligands in selective drug therapies

Chimeric drugs with selective potential toward specific cell types constitute one of the most promising forefronts of modern Pharmacology. We present a mathematical model to test and optimize these synthetic constructs, as an alternative to conventional empirical design. We take as a case study a chimeric construct composed of epidermal growth factor (EGF) linked to different mutants of interferon (IFN). Our model quantitatively reproduces all the experimental results, illustrating how chimeras using mutants of IFN with reduced affinity exhibit enhanced selectivity against cell overexpressing EGF receptor. We also investigate how chimeric selectivity can be improved based on the balance between affinity rates, receptor abundance, activity of ligand subunits, and linker length between subunits. The simplicity and generality of the model facilitate a straightforward application to other chimeric constructs, providing a quantitative systematic design and optimization of these selective drugs against certain cell‐based diseases, such as Alzheimers and cancer.

Mitochondria and the non‐genetic origins of cell‐to‐cell variability: More is different

Gene expression activity is heterogeneous in a population of isogenic cells. Identifying the molecular basis of this variability will improve our understanding of phenomena like tumor resistance to drugs, virus infection, or cell fate choice. The complexity of the molecular steps and machines involved in transcription and translation could introduce sources of randomness at many levels, but a common constraint to most of these processes is its energy dependence. In eukaryotic cells, most of this energy is provided by mitochondria. A clonal population of cells may show a large variability in the number and functionality of mitochondria. Here, we discuss how differences in the mitochondrial content of each cell contribute to heterogeneity in gene products. Changes in the amount of mitochondria can also entail drastic alterations of a cells gene expression program, which ultimately leads to phenotypic diversity.

Positive regulatory dynamics by a small noncoding RNA: speeding up responses under temperature stress

Recent discoveries of noncoding regulatory RNAs have led to further understanding of the elements controlling genetic expression. In E. coli, most of those ncRNAs for which functional knowledge is available were shown to be dependent on the Hfq RNA chaperone and to act as inhibitors of translation by base pairing with their mRNA target. Nevertheless, there are also some examples where the sRNA plays a role of a translational activator, structurally enhancing ribosome binding to mRNA. In this work, we seek to understand the dynamics of DsrA-based positive regulation of rpoS mRNA, encoding the σ(S) RNA polymerase subunit, and to understand how it helps to mitigate environmental stress in bacteria. Our analysis is based on the first absolute quantification of the copy number of both the sRNA and of its corresponding mRNA in combination with mathematical models for post-transcriptional regulation. We show that on average, DsrA is present at a ratio of 3 to 24 copies per cell, while an rpoS transcript is present at a level of 1 to 4 copies per cell, both levels increasing when temperature is decreased. Our analysis supports the idea that temperature dependency of DsrA degradation is not a crucial condition for the attainment of observed DsrA steady levels, but highlights that this may have a marked influence on the dynamics of the regulation, notably to speed up the time of recovery to normal RNA levels after ending the stress signal. Further, our analysis also reveals how reversibility of RNA complex formation and σ(S)-regulated degradation act to reduce intrinsic noise in σ(S) induction. Taking into account the importance of this master regulator, which allows E. coli as well as other important pathogens to survive their environment, the present work contributes to complete the panel of multiple signals used to regulate bacterial transcription.

Trade-offs and noise tolerance in signal detection by genetic circuits

Genetic circuits can implement elaborated tasks of amplitude or frequency signal detection. What type of constraints could circuits experience in the performance of these tasks, and how are they affected by molecular noise? Here, we consider a simple detection process–a signal acting on a two-component module–to analyze these issues. We show that the presence of a feedback interaction in the detection module imposes a trade-off on amplitude and frequency detection, whose intensity depends on feedback strength. A direct interaction between the signal and the output species, in a type of feed-forward loop architecture, greatly modifies these trade-offs. Indeed, we observe that coherent feed-forward loops can act simultaneously as good frequency and amplitude noise-tolerant detectors. Alternatively, incoherent feed-forward loop structures can work as high-pass filters improving high frequency detection, and reaching noise tolerance by means of noise filtering. Analysis of experimental data from several specific coherent and incoherent feed-forward loops shows that these properties can be realized in a natural context. Overall, our results emphasize the limits imposed by circuit structure on its characteristic stimulus response, the functional plasticity of coherent feed-forward loops, and the seemingly paradoxical advantage of improving signal detection with noisy circuit components.

Mitochondrial levels determine variability in cell death by modulating apoptotic gene expression

Fractional killing is the main cause of tumour resistance to chemotherapy. This phenomenon is observed even in genetically identical cancer cells in homogeneous microenvironments. To understand this variable resistance, here we investigate the individual responses to TRAIL in a clonal population of HeLa cells using live-cell microscopy and computational modelling. We show that the cellular mitochondrial content determines the apoptotic fate and modulates the time to death, cells with higher mitochondrial content are more prone to die. We find that all apoptotic protein levels are modulated by the mitochondrial content. Modelling the apoptotic network, we demonstrate that these correlations, and especially the differential control of anti- and pro-apoptotic protein pairs, confer mitochondria a powerful discriminatory capacity of apoptotic fate. We find a similar correlation between the mitochondria and apoptotic proteins in colon cancer biopsies. Our results reveal a different role of mitochondria in apoptosis as the global regulator of apoptotic protein expression.It is unclear what causes variation in cell death in response to chemotherapy. Here, the authors show that cellular mitochondrial content modulates apoptotic protein levels, which in turn regulates response to agents such as TRAIL.

Validation of the Hirst-Type Spore Trap for Simultaneous Monitoring of Prokaryotic and Eukaryotic Biodiversities in Urban Air Samples by Next-Generation Sequencing

ABSTRACT Pollen, fungi, and bacteria are the main microscopic biological entities present in outdoor air, causing allergy symptoms and disease transmission and having a significant role in atmosphere dynamics. Despite their relevance, a method for monitoring simultaneously these biological particles in metropolitan environments has not yet been developed. Here, we assessed the use of the Hirst-type spore trap to characterize the global airborne biota by high-throughput DNA sequencing, selecting regions of the 16S rRNA gene and internal transcribed spacer for the taxonomic assignment. We showed that aerobiological communities are well represented by this approach. The operational taxonomic units (OTUs) of two traps working synchronically compiled >87% of the total relative abundance for bacterial diversity collected in each sampler, >89% for fungi, and >97% for pollen. We found a good correspondence between traditional characterization by microscopy and genetic identification, obtaining more-accurate taxonomic assignments and detecting a greater diversity using the latter. We also demonstrated that DNA sequencing accurately detects differences in biodiversity between samples. We concluded that high-throughput DNA sequencing applied to aerobiological samples obtained with Hirst spore traps provides reliable results and can be easily implemented for monitoring prokaryotic and eukaryotic entities present in the air of urban areas. IMPORTANCE Detection, monitoring, and characterization of the wide diversity of biological entities present in the air are difficult tasks that require time and expertise in different disciplines. We have evaluated the use of the Hirst spore trap (an instrument broadly employed in aerobiological studies) to detect and identify these organisms by DNA-based analyses. Our results showed a consistent collection of DNA and a good concordance with traditional methods for identification, suggesting that these devices can be used as a tool for continuous monitoring of the airborne biodiversity, improving taxonomic resolution and characterization together. They are also suitable for acquiring novel DNA amplicon-based information in order to gain a better understanding of the biological particles present in a scarcely known environment such as the air.