Plasticity facilitates sustainable growth in the commons

Matteo Cavaliere and Juan F. Poyatos (2013). Journal of the Royal Society Interface 10(81):20121006.

In the commons, communities whose growth depends on public good, individuals often rely on surprisingly simple strategies, or heuristics, to decide whether to contribute to the shared resource (at risk of exploitation by free-riders). Although this appears a limitation, we show here how four heuristics lead to sustainable growth when coupled to specific ecological constraints. The two simplest ones—contribute permanently or switch stochastically between contributing or not—are first shown to bring sustainability when the public good efficiently promotes growth. If efficiency declines and the commons is structured in small groups, the most effective strategy resides in contributing only when a majority of individuals are also contributors. In contrast, when group size becomes large, the most effective behaviour follows a minimal-effort rule: contribute only when it is strictly necessary. Both plastic strategies are observed in natural scenarios across scales that present them as relevant social motifs for the sustainable management of public goods.

Protein Disorder in Arabidopsis thaliana - Implications for Plant Environmental Adaptation

Natalia Pietrosemoli, Juan A. García-Martín, Roberto Solano and Florencio Pazos (2013). PLoS One 8(2):e55524.

Intrinsically Disordered Proteins/Regions (IDPs/IDRs) lack a fixed 3D structure in their native (functional) form. These highly mobile polypeptides are generally involved in functions associated to transient protein interactions, since their conformational flexibility makes them very suitable for this kind of molecular recognition. For this reason, IDPs are involved in key cellular process such as cell-cycle control and signaling, where this type of molecular interactions are frequent. Due to their involvement in important cellular processes, the repertory of IDPs in a given proteome provides insight into important biological aspects of the corresponding organism. In general, high IDP content has been associated to complexity, both at the organism and cellular process levels.

Since IDPs/IDRs are frequently involved in transient interactions, introducing disorder in a given system is a way of increasing its "wiring" (the number of interactions/connections) and consequently the complexity. The proteome of A. thaliana contains less IDPs than that of Human. Nevertheless, some particular processes, mainly related to response to external stimulus contain more disorder in the plant than in Human. We propose that since plants cannot escape from environmental conditions as animals do, they use disorder as a simple and fast mechanism, independent of transcriptional control, for introducing versatility in the interaction networks underlying these biological processes so that they can quickly adapt and respond to challenging environmental conditions.

Balancing noise and plasticity in eukaryotic gene expression

Djordje Bajić and Juan F Poyatos (2012). BMC Genomics 13:343.

BACKGROUND:Coupling the control of expression stochasticity (noise) to the ability of expression change (plasticity) can alter gene function and influence adaptation. A number of factors, such as transcription re-initiation, strong chromatin regulation or genome neighboring organization, underlie this coupling. However, these factors do not necessarily combine in equivalent ways and strengths in all genes. Can we identify then alternative architectures that modulate in distinct ways the linkage of noise and plasticity?

RESULTS:Here we first show that strong chromatin regulation, commonly viewed as a source of coupling, can lead to plasticity without noise. The nature of this regulation is relevant too, with plastic but noiseless genes being subjected to general activators whereas plastic and noisy genes experience more specific repression. Contrarily, in genes exhibiting poor transcriptional control, it is translational efficiency what separates noise from plasticity, a pattern related to transcript length. This additionally implies that genome neighboring organization -as modifier- appears only effective in highly plastic genes. In this class, we confirm bidirectional promoters (bipromoters) as a configuration capable to reduce coupling by abating noise but also reveal an important trade-off, since bipromoters also decrease plasticity. This presents ultimately a paradox between intergenic distances and modulation, with short intergenic distances both associated and disassociated to noise at different plasticity levels.

CONCLUSIONS:Balancing the coupling among different types of expression variability appears as a potential shaping force of genome regulation and organization. This is reflected in the use of different control strategies at genes with different sets of functional constraints.

Interactive calculation and visualization of function-related conservation patterns in multiple sequence alignments and structures

Thilo Muth, Juan A. García-Martín, Antonio Rausell, David Juan, Alfonso Valencia and Florencio Pazos (2012). Bioinformatics 28(4):584-6.

JDet is a multiplatform software for the interactive calculation and visualization of function-related conservation patterns in multiple sequence alignments and structures. It contains the set of tools and features we consider critical for the daily work with this kind of data, and that previously were disseminated in different packages and web servers. The package includes two of our recently developed programs for extracting this kind of information from protein alignments. JDet can be extended via plugins to implement other methods or parse their results.

Towards the prediction of protein interaction partners using physical docking

Mark N. Wass, Gloria Fuentes, Carles Pons, Florencio Pazos and Alfonso Valencia. (2011). Mol Syst Biol. 7:469.

The prediction of protein interaction partners (which protein interacts with which one) and the prediction of the structure of a complex between two proteins known to interact have so far been regarded as two separated issues. Physical docking programs are intended for the second goal: predict the arrangement of two protein chains of known 3D structure when they interact. These programs produce a list of complexes (docking models), representing different poses/arrangements of the two protein chains, sorted by an score representing, in the simplest case, the geometric complementarity between the two chains in these arrangements. In general, these programs have a high degree of error and the native (real) complex is usually not among the top scoring ones in this list. Additional information and constraints are usually needed to filter these lists in order to locate the native complex. This relatively poor performance, together with the intrinsically high computational requirements of these programs, precluded their usage in the search for interaction partners to date.

We have shown that, in spite of this low accuracy of the basic docking methods, in a significant number of cases they can be used to locate interaction partners within large sets of alternative structures with diverse characteristics. The distribution of raw scores of the list of complexes is shifted to better scores for the docking experiment involving a real interaction partner, than for those involving likely non-interacting chains. Paradoxically, this happens even in cases where the docking experiment involving the truly interacting chain fails in locating the real complex

A possible explanation for this observation could be related to the "energy-funnel" theory of protein interaction. According with this theory, not only the real complex between two interacting proteins has a favourable energy but there are many alternative complexes with favourable energies which help in "directing" the interaction to the real one. It is difficult to think that only the real complex has a favourable energy and that consequently the interaction process has to rely purely on random encounters until it is found, taking into account the millions of alternatives in which two proteins can face each other. Other conformations "around" the real one might also be favoured in a way that once the two chains find any of them, their way to the real complex is facilitated. That would explain the fact that, even if the real complex is not found by the docking program, these complexes are enough for shifting the energy distribution.

The balance of weak and strong interactions in genetic networks

Poyatos, J. F. (2011). PLoS ONE, 6, e14598.

Genetic interactions are being quantitatively characterized in a comprehensive way in several model organisms. These data are then globally represented in terms of genetic networks. How are interaction strengths distributed in these networks? And what type of functional organization of the underlying genomic systems is revealed by such distribution patterns?

This work demonstrates a non-random organization of interaction strengths in genetic networks of Caenorhabditis elegans and Saccharomyces. cerevisiae, with implications for the structure of genetic buffering, the absence of strong genetic cascades, the functional role of genetic hubs and the phenotypic/evolutionary features of the genes involved. This non-random organization, common to other complex networks, could ultimately reflect how genetic variation is eventually influencing the phenotype.

MBRole: enrichment analysis of metabolomic data

Chagoyen M, Pazos F. (2011). Bioinformatics, 27(5):730-1.

MBRole is a web server for the functional interpretation of metabolomics results. MBRole performs overrepresentation (enrichment) analysis of annotations for a set of compounds. The annotations include biological and chemical information obtained from different public databases (like KEGG, PubChem, HMDB...).

pBAM1: an all-synthetic genetic tool for analysis and construction of complex bacterial phenotypes

Martinez-Garcia E, Calles B, Arevalo-Rodriguez M, de Lorenzo V (2011). BMC Microbiol 11: 38.

This paper reports the design and performance of two all-synthetic mini-transposon vectors, pBAM1 and pBAM1-GFP. The power of the both systemsis shown in the paper with a couple of examples.

pBAM1 plasmid map. Functional elements of the plasmid include relevant restriction sites, antibiotic markers (Ap, ampicillin, Km, kanamycin), transposase (tnpA), origin of replication (R6K), the origin of transfer region (oriT), mosaic element O (ME-O), and mosaic element I (ME-I), are shown.

Regulatory exaptation of the catabolite repression protein (Crp)-cAMP system in Pseudomonas putida

Milanesio P, Arce-Rodriguez A, Munoz A, Calles B, de Lorenzo V (2011) Environ Microbiol 13: 324-339.

In this paper, we described the function of the Crp-cAMP system of P. putida. We found that both components of the system, i.e. the genes crp and cyaA (coding for the cAMP receptor protein and adenylate cyclase, respectively) produce fully functional proteins that are not involved in the carbon catabolite repression fenomena, like the E. coli counterparts. In contrast, this system is controlling cell surface related functions, like resistance to antibiotics, toxic compounds and utilization of dipeptides as C and N source. We propose that this dissimilar functional scope of Crp-cAMP system in bacteria is a case of a molecular regulatory exaptation mechanism, i.e. the process through which a property evolved for a particular function is co-opted for a new use.

Fructose 1-Phosphate is the preferred effector of the metabolic regulator Cra of Pseudomonas putida

Chavarria M, Santiago C, Platero R, Krell T, Casasnovas JM, et al. (2011). J Biol Chem.

The catabolite repressor/activator (Cra) protein is a global sensor and regulator of carbon fluxes through the central metabolic pathways of Gram-negative bacteria. Our results showed that Cra regulates the PTSFru system in P. putida by binding strongly as a dimer to a single palindromic sequence within the cra/fruB intergenic region and that only F1P can lift such down-regulation.

The interplay of the EIIA(Ntr) component of the nitrogen-related phosphotransferase system (PTS(Ntr)) of Pseudomonas putida with pyruvate dehydrogenase.

Pfluger-Grau K, Chavarria M, de Lorenzo V (2011) Biochim Biophys Acta.

Both genetic and biochemical evidence revealed that the non-phosphorylated form of PtsN down-regulates the Pyruvate dehydrogenase activity by means of direct protein–protein interactions with AceE.This result suggests that EIIANtr takes part in the node of C metabolism that checks the flux of carbon from carbohydrates into the Krebs cycle.

Monitoring biodegradative enzymes with nanobodies raised in Camelus dromedarius with mixtures of catabolic proteins.

Zafra O, Fraile S, Gutierrez C, Haro A, Paez-Espino AD, et al. (2011) Environ Microbiol.

Generation of nanobody libraries workflow. The sketch summarized the main steps of the process. African camels (C. dromedarius) are inoculated with a mixture of proteins and following a 6 week period, RNA from lymphocytes is extracted, retrotrancribed to DNA and the sequences corresponding to the VHH domains amplified with suitable primers for capture in a phage display vector. This library is then separately subject to two rounds of panning on plates with the individual antigens of the original inoculation cocktail. Those M13-VHH clones that retained a strong binding activity, as detected with a phage-ELISA assay were kept for further analysis.

Cooperative amino acid changes shift the response of the sigma(54) -dependent regulator XylR from natural m-xylene towards xenobiotic 2,4-dinitrotoluene.

de Las Heras A, de Lorenzo V (2011). Mol Microbiol 79: 1248-1259.

The A domain of the transcriptional regulator XylR of Pseudomonas putida can be evolved in the Laboratory to recreate the itinerary that goes from one exisiting effector specificity towards a new one. The figure shows amino acids that participate in such changes in effector specificity

Genes that move the window of viability of life: lessons from bacteria thriving at the cold extreme: mesophiles can be turned into extremophiles by substituting essential genes.

de Lorenzo V (2011). Bioessays 33: 38-42.

The Biotechnology of the future will increasingly depend on combining biologica activities mined from the environment (eg metagenomes) with genomic chasses engineered with the tools of Synthetic Biology.

Linking genes to microbial growth kinetics-An integrated biochemical systems engineering approach.

Koutinas M, Kiparissides A, Silva-Rocha R, Lam MC, Martins Dos Santos VA, et al. (2011) Metab Eng.

In this paper, we created a complete model of the TOL regulatory network, including both regulatory and metabolic processes, using set ordinary differential equations. Literature mining together with specifically designed experiments allows the determination of the missing kinetic parameters necessary to feed the model. Simulations of the systems allow the determination of critical regulatory steps with major contribution to the network operation. Also, predictive simulation (using in silico perturbation of the network) were performed and the results were nicely reproduced in validation experiments. We propose a general integrative strategy for the modeling of hybrid regulatory-metabolic network applicable to biodegradation process.

Improving the prediction of Pseudomonas putida mt-2 growth kinetics with the use of a gene expression regulation model of the TOL plasmid

Koutinas M, Kiparissides A, Lam MC, Silva-Rocha R, Godinho M, de Lorenzo V, Martins Dos Santos VA, Pistikopoulos EN, Mantalaris A. (2011) Biochemical Engineering Journal.

Although Boolean models are useful when few kinetic constants are available to a given network, they cannot provide high resolution information about the circuit operation. Such information is desired when regulatory/metabolic optimization is the final goal of the work, which is often the case in the field of bioengineering. In the paper, we use a model based on ordinary differential equations to try to link the dynamics of gene activation in the TOL network to biomass production and subtract consumption, final outcome of the process. The model focus only in the interplay between the two regulators of the system (XylR and XylS) and demonstrated that the outcome of the system can be successfully explained by analyze the dynamics of gene activation at this layer of the network.

The logicome of environmental bacteria: merging catabolic and regulatory events with Boolean formalisms

Silva-Rocha R, Tamames J, Dos Santos VM, de Lorenzo V. (2011) Environ Microbiol.

In this work we address the usage of Boolean formalism for the modeling of complex regulatory networks in bacteria. Using the pathway for degradation of toluene/m-xylene of Pseudomonas putida mt-2 as an example, we explore the potential of such formalism as a convenient tool to inspect network behavior when little kinetic information is available for the system. We argue that logic gate-driven hypothesis generation could significantly impulse the field of biodegradation.


Systems Biology Program @CNB-CSIC

c/ Darwin, 3 - 28049 Madrid - Spain