A dynamic network of transcription in LPS-treated human subjects

<p>Abstract</p> <p>Background</p> <p>Understanding the transcriptional regulatory networks that map out the coordinated dynamic responses of signaling proteins, transcription factors and target genes over time would represent a significant advance in the application of genome wide expression analysis. The primary challenge is monitoring transcription factor activities over time, which is not yet available at the large scale. Instead, there have been several developments to estimate activities computationally. For example, Network Component Analysis (NCA) is an approach that can predict transcription factor activities over time as well as the relative regulatory influence of factors on each target gene.</p> <p>Results</p> <p>In this study, we analyzed a gene expression data set in blood leukocytes from human subjects administered with lipopolysaccharide (LPS), a prototypical inflammatory challenge, in the context of a reconstructed regulatory network including 10 transcription factors, 99 target genes and 149 regulatory interactions. We found that the computationally estimated activities were well correlated to their coordinated action. Furthermore, we found that clustering the genes in the context of regulatory influences greatly facilitated interpretation of the expression data, as clusters of gene expression corresponded to the activity of specific factors or more interestingly, factor combinations which suggest coordinated regulation of gene expression. The resulting clusters were therefore more biologically meaningful, and also led to identification of additional genes under the same regulation.</p> <p>Conclusion</p> <p>Using NCA, we were able to build a network that accounted for between 8–11% genes in the known transcriptional response to LPS in humans. The dynamic network illustrated changes of transcription factor activities and gene expressions as well as interactions of signaling proteins, transcription factors and target genes.</p

Amino-acid site variability among natural and designed proteins

Computational protein design attempts to create protein sequences that fold stably into pre-specified structures. Here we compare alignments of designed proteins to alignments of natural proteins and assess how closely designed sequences recapitulate patterns of sequence variation found in natural protein sequences. We design proteins using RosettaDesign, and we evaluate both fixed-backbone designs and variable-backbone designs with different amounts of backbone flexibility. We find that proteins designed with a fixed backbone tend to underestimate the amount of site variability observed in natural proteins while proteins designed with an intermediate amount of backbone flexibility result in more realistic site variability. Further, the correlation between solvent exposure and site variability in designed proteins is lower than that in natural proteins. This finding suggests that site variability is too uniform across different solvent exposure states (i.e., buried residues are too variable or exposed residues too conserved). When comparing the amino acid frequencies in the designed proteins with those in natural proteins we find that in the designed proteins hydrophobic residues are underrepresented in the core. From these results we conclude that intermediate backbone flexibility during design results in more accurate protein design and that either scoring functions or backbone sampling methods require further improvement to accurately replicate structural constraints on site variability

Mycobacterium intermedium Granulomatous Dermatitis from Hot Tub Exposure

Nontuberculous mycobacteria, which are widespread in the environment, frequently cause opportunistic infections in immunocompromised patients. We report the first case of a patient with chronic granulomatous dermatitis caused by a rarely described organism, Mycobacterium intermedium. The infection was associated with exposure in a home hot tub

Atmospheric sulfur cycling in the southeastern Pacific – longitudinal distribution, vertical profile, and diel variability observed during VOCALS-REx

Dimethylsulfide (DMS) emitted from the ocean is a biogenic precursor gas for sulfur dioxide (SO&lt;sub&gt;2&lt;/sub&gt;) and non-sea-salt sulfate aerosols (SO&lt;sub&gt;4&lt;/sub&gt;&lt;sup&gt;2&amp;minus;&lt;/sup&gt;). During the VAMOS-Ocean-Cloud-Atmosphere-Land Study Regional Experiment (VOCALS-REx) in 2008, multiple instrumented platforms were deployed in the Southeastern Pacific (SEP) off the coast of Chile and Peru to study the linkage between aerosols and stratocumulus clouds. We present here observations from the NOAA Ship &lt;i&gt;Ronald H. Brown&lt;/i&gt; and the NSF/NCAR C-130 aircraft along ~20° S from the coast (70° W) to a remote marine atmosphere (85° W). While SO&lt;sub&gt;4&lt;/sub&gt;&lt;sup&gt;2&amp;minus;&lt;/sup&gt; and SO&lt;sub&gt;2&lt;/sub&gt; concentrations were distinctly elevated above background levels in the coastal marine boundary layer (MBL) due to anthropogenic influence (~800 and 80 pptv, respectively), their concentrations rapidly decreased west of 78° W (~100 and 25 pptv). In the remote region, entrainment from the free troposphere (FT) increased MBL SO&lt;sub&gt;2&lt;/sub&gt; burden at a rate of 0.05 &amp;plusmn; 0.02 μmoles m&lt;sup&gt;&amp;minus;2&lt;/sup&gt; day&lt;sup&gt;&amp;minus;1&lt;/sup&gt; and diluted MBL SO&lt;sub&gt;4&lt;/sub&gt;&lt;sup&gt;2&lt;/sup&gt; burden at a rate of 0.5 &amp;plusmn; 0.3 μmoles m&lt;sup&gt;&amp;minus;2&lt;/sup&gt; day&lt;sup&gt;&amp;minus;1&lt;/sup&gt;, while the sea-to-air DMS flux (3.8 &amp;plusmn; 0.4 μmoles m&lt;sup&gt;&amp;minus;2&lt;/sup&gt; day&lt;sup&gt;&amp;minus;1&lt;/sup&gt;) remained the predominant source of sulfur mass to the MBL. In-cloud oxidation was found to be the most important mechanism for SO&lt;sub&gt;2&lt;/sub&gt; removal and in situ SO&lt;sub&gt;4&lt;/sub&gt;&lt;sup&gt;2&amp;minus;&lt;/sup&gt; production. Surface SO&lt;sub&gt;4&lt;/sub&gt;&lt;sup&gt;2&amp;minus;&lt;/sup&gt; concentration in the remote MBL displayed pronounced diel variability, increasing rapidly in the first few hours after sunset and decaying for the rest of the day. We theorize that the increase in SO&lt;sub&gt;4&lt;/sub&gt;&lt;sup&gt;2&amp;minus;&lt;/sup&gt; was due to nighttime recoupling of the MBL that mixed down cloud-processed air, while decoupling and sporadic precipitation scavenging were responsible for the daytime decline in SO&lt;sub&gt;4&lt;/sub&gt;&lt;sup&gt;2&amp;minus;&lt;/sup&gt;

Incorporation of enzyme concentrations into FBA and identification of optimal metabolic pathways

<p>Abstract</p> <p>Background</p> <p>In the present article, we propose a method for determining optimal metabolic pathways in terms of the level of concentration of the enzymes catalyzing various reactions in the entire metabolic network. The method, first of all, generates data on reaction fluxes in a pathway based on steady state condition. A set of constraints is formulated incorporating weighting coefficients corresponding to concentration of enzymes catalyzing reactions in the pathway. Finally, the rate of yield of the target metabolite, starting with a given substrate, is maximized in order to identify an optimal pathway through these weighting coefficients.</p> <p>Results</p> <p>The effectiveness of the present method is demonstrated on two synthetic systems existing in the literature, two pentose phosphate, two glycolytic pathways, core carbon metabolism and a large network of carotenoid biosynthesis pathway of various organisms belonging to different phylogeny. A comparative study with the existing extreme pathway analysis also forms a part of this investigation. Biological relevance and validation of the results are provided. Finally, the impact of the method on metabolic engineering is explained with a few examples.</p> <p>Conclusions</p> <p>The method may be viewed as determining an optimal set of enzymes that is required to get an optimal metabolic pathway. Although it is a simple one, it has been able to identify a carotenoid biosynthesis pathway and the optimal pathway of core carbon metabolic network that is closer to some earlier investigations than that obtained by the extreme pathway analysis. Moreover, the present method has identified correctly optimal pathways for pentose phosphate and glycolytic pathways. It has been mentioned using some examples how the method can suitably be used in the context of metabolic engineering.</p

Prospecting environmental mycobacteria: combined molecular approaches reveal unprecedented diversity

Background: Environmental mycobacteria (EM) include species commonly found in various terrestrial and aquatic environments, encompassing animal and human pathogens in addition to saprophytes. Approximately 150 EM species can be separated into fast and slow growers based on sequence and copy number differences of their 16S rRNA genes. Cultivation methods are not appropriate for diversity studies; few studies have investigated EM diversity in soil despite their importance as potential reservoirs of pathogens and their hypothesized role in masking or blocking M. bovis BCG vaccine. Methods: We report here the development, optimization and validation of molecular assays targeting the 16S rRNA gene to assess diversity and prevalence of fast and slow growing EM in representative soils from semi tropical and temperate areas. New primer sets were designed also to target uniquely slow growing mycobacteria and used with PCR-DGGE, tag-encoded Titanium amplicon pyrosequencing and quantitative PCR. Results: PCR-DGGE and pyrosequencing provided a consensus of EM diversity; for example, a high abundance of pyrosequencing reads and DGGE bands corresponded to M. moriokaense, M. colombiense and M. riyadhense. As expected pyrosequencing provided more comprehensive information; additional prevalent species included M. chlorophenolicum, M. neglectum, M. gordonae, M. aemonae. Prevalence of the total Mycobacterium genus in the soil samples ranged from 2.3×107 to 2.7×108 gene targets g−1; slow growers prevalence from 2.9×105 to 1.2×107 cells g−1. Conclusions: This combined molecular approach enabled an unprecedented qualitative and quantitative assessment of EM across soil samples. Good concordance was found between methods and the bioinformatics analysis was validated by random resampling. Sequences from most pathogenic groups associated with slow growth were identified in extenso in all soils tested with a specific assay, allowing to unmask them from the Mycobacterium whole genus, in which, as minority members, they would have remained undetected

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

An analytic and systematic framework for estimating metabolic flux ratios from 13C tracer experiments

<p>Abstract</p> <p>Background</p> <p>Metabolic fluxes provide invaluable insight on the integrated response of a cell to environmental stimuli or genetic modifications. Current computational methods for estimating the metabolic fluxes from ¹³<it>C </it>isotopomer measurement data rely either on manual derivation of analytic equations constraining the fluxes or on the numerical solution of a highly nonlinear system of isotopomer balance equations. In the first approach, analytic equations have to be tediously derived for each organism, substrate or labelling pattern, while in the second approach, the global nature of an optimum solution is difficult to prove and comprehensive measurements of external fluxes to augment the ¹³<it>C </it>isotopomer data are typically needed.</p> <p>Results</p> <p>We present a novel analytic framework for estimating metabolic flux ratios in the cell from ¹³<it>C </it>isotopomer measurement data. In the presented framework, equation systems constraining the fluxes are derived automatically from the model of the metabolism of an organism. The framework is designed to be applicable with all metabolic network topologies, ¹³<it>C </it>isotopomer measurement techniques, substrates and substrate labelling patterns.</p> <p>By analyzing nuclear magnetic resonance (NMR) and mass spectrometry (MS) measurement data obtained from the experiments on glucose with the model micro-organisms <it>Bacillus subtilis </it>and <it>Saccharomyces cerevisiae </it>we show that our framework is able to automatically produce the flux ratios discovered so far by the domain experts with tedious manual analysis. Furthermore, we show by <it>in silico </it>calculability analysis that our framework can rapidly produce flux ratio equations – as well as predict when the flux ratios are unobtainable by linear means – also for substrates not related to glucose.</p> <p>Conclusion</p> <p>The core of ¹³<it>C </it>metabolic flux analysis framework introduced in this article constitutes of flow and independence analysis of metabolic fragments and techniques for manipulating isotopomer measurements with vector space techniques. These methods facilitate efficient, analytic computation of the ratios between the fluxes of pathways that converge to a common junction metabolite. The framework can been seen as a generalization and formalization of existing tradition for computing metabolic flux ratios where equations constraining flux ratios are manually derived, usually without explicitly showing the formal proofs of the validity of the equations.</p

Tension and Robustness in Multitasking Cellular Networks

Cellular networks multitask by exhibiting distinct, context-dependent dynamics. However, network states (parameters) that generate a particular dynamic are often sub-optimal for others, defining a source of “tension” between them. Though multitasking is pervasive, it is not clear where tension arises, what consequences it has, and how it is resolved. We developed a generic computational framework to examine the source and consequences of tension between pairs of dynamics exhibited by the well-studied RB-E2F switch regulating cell cycle entry. We found that tension arose from task-dependent shifts in parameters associated with network modules. Although parameter sets common to distinct dynamics did exist, tension reduced both their accessibility and resilience to perturbation, indicating a trade-off between “one-size-fits-all” solutions and robustness. With high tension, robustness can be preserved by dynamic shifting of modules, enabling the network to toggle between tasks, and by increasing network complexity, in this case by gene duplication. We propose that tension is a general constraint on the architecture and operation of multitasking biological networks. To this end, our work provides a framework to quantify the extent of tension between any network dynamics and how it affects network robustness. Such analysis would suggest new ways to interfere with network elements to elucidate the design principles of cellular networks