Predictive Approach Identifies Molecular Targets and Interventions to Restore Angiogenesis in Wounds With Delayed Healing

Impaired angiogenesis is a hallmark of wounds with delayed healing, and currently used therapies to restore angiogenesis have limited efficacy. Here, we employ a computational simulation-based approach to identify influential molecular and cellular processes, as well as protein targets, whose modulation may stimulate angiogenesis in wounds. We developed a mathematical model that captures the time courses for platelets, 9 cell types, 29 proteins, and oxygen, which are involved in inflammation, proliferation, and angiogenesis during wound healing. We validated our model using previously published experimental data. By performing global sensitivity analysis on thousands of simulated wound-healing scenarios, we identified six processes (among the 133 modeled in total) whose modulation may improve angiogenesis in wounds. By simulating knockouts of 25 modeled proteins and by simulating different wound-oxygenation levels, we identified four proteins [namely, transforming growth factor (TGF)-β, vascular endothelial growth factor (VEGF), fibroblast growth factor-2 (FGF-2), and angiopoietin-2 (ANG-2)], as well as oxygen, as therapeutic targets for stimulating angiogenesis in wounds. Our modeling results indicated that simultaneous inhibition of TGF-β and supplementation of either FGF-2 or ANG-2 could be more effective in stimulating wound angiogenesis than the modulation of either protein alone. Our findings suggest experimentally testable intervention strategies to restore angiogenesis in wounds with delayed healing

Positive Autoregulation Shapes Response Timing and Intensity in Two-component Signal Transduction Systems

Positive feedback loops are regulatory elements that can modulate expression output, kinetics and noise in genetic circuits. Transcriptional regulators participating in such loops are often expressed from two promoters, one constitutive and one autoregulated. Here, we investigate the interplay of promoter strengths and the intensity of the stimulus activating the transcriptional regulator in defining the output of a positively autoregulated genetic circuit. Using a mathematical model of twocomponent regulatory systems, which are present in all domains of life, we establish that positive feedback strongly affects the steady-state output levels at both low and high levels of stimulus if the constitutive promoter of the regulator is weak. By contrast, the effect of positive feedback is negligible when the constitutive promoter is sufficiently strong, unless the stimulus intensity is very high. Furthermore, we determine that positive feedback can affect both transient and steady state output levels even in the simplest genetic regulatory systems. We tested our modeling predictions by abolishing the positive feedback loop in the two-component regulatory system PhoP/PhoQ of Salmonella enterica, which resulted in diminished induction of PhoP-activated genes

Evolution and Dynamics of Regulatory Architectures Controlling Polymyxin B Resistance in Enteric Bacteria

Complex genetic networks consist of structural modules that determine the levels and timing of a cellular response. While the functional properties of the regulatory architectures that make up these modules have been extensively studied, the evolutionary history of regulatory architectures has remained largely unexplored. Here, we investigate the transition between direct and indirect regulatory pathways governing inducible resistance to the antibiotic polymyxin B in enteric bacteria. We identify a novel regulatory architecture—designated feedforward connector loop—that relies on a regulatory protein that connects signal transduction systems post-translationally, allowing one system to respond to a signal activating another system. The feedforward connector loop is characterized by rapid activation, slow deactivation, and elevated mRNA expression levels in comparison with the direct regulation circuit. Our results suggest that, both functionally and evolutionarily, the feedforward connector loop is the transitional stage between direct transcriptional control and indirect regulation

Detection of distant evolutionary relationships between protein families using theory of sequence profile-profile comparison

<p>Abstract</p> <p>Background</p> <p>Detection of common evolutionary origin (homology) is a primary means of inferring protein structure and function. At present, comparison of protein families represented as sequence profiles is arguably the most effective homology detection strategy. However, finding the best way to represent evolutionary information of a protein sequence family in the profile, to compare profiles and to estimate the biological significance of such comparisons, remains an active area of research.</p> <p>Results</p> <p>Here, we present a new homology detection method based on sequence profile-profile comparison. The method has a number of new features including position-dependent gap penalties and a global score system. Position-dependent gap penalties provide a more biologically relevant way to represent and align protein families as sequence profiles. The global score system enables an analytical solution of the statistical parameters needed to estimate the statistical significance of profile-profile similarities. The new method, together with other state-of-the-art profile-based methods (HHsearch, COMPASS and PSI-BLAST), is benchmarked in all-against-all comparison of a challenging set of SCOP domains that share at most 20% sequence identity. For benchmarking, we use a reference ("gold standard") free model-based evaluation framework. Evaluation results show that at the level of protein domains our method compares favorably to all other tested methods. We also provide examples of the new method outperforming structure-based similarity detection and alignment. The implementation of the new method both as a standalone software package and as a web server is available at <url>http://www.ibt.lt/bioinformatics/coma</url>.</p> <p>Conclusion</p> <p>Due to a number of developments, the new profile-profile comparison method shows an improved ability to match distantly related protein domains. Therefore, the method should be useful for annotation and homology modeling of uncharacterized proteins.</p

Speed, Sensitivity, and Bistability in Auto-activating Signaling Circuits

Cells employ a myriad of signaling circuits to detect environmental signals and drive specific gene expression responses. A common motif in these circuits is inducible auto-activation: a transcription factor that activates its own transcription upon activation by a ligand or by post-transcriptional modification. Examples range from the two-component signaling systems in bacteria and plants to the genetic circuits of animal viruses such as HIV. We here present a theoretical study of such circuits, based on analytical calculations, numerical computations, and simulation. Our results reveal several surprising characteristics. They show that auto-activation can drastically enhance the sensitivity of the circuit's response to input signals: even without molecular cooperativity, an ultra-sensitive threshold response can be obtained. However, the increased sensitivity comes at a cost: auto-activation tends to severely slow down the speed of induction, a stochastic effect that was strongly underestimated by earlier deterministic models. This slow-induction effect again requires no molecular cooperativity and is intimately related to the bimodality recently observed in non-cooperative auto-activation circuits. These phenomena pose strong constraints on the use of auto-activation in signaling networks. To achieve both a high sensitivity and a rapid induction, an inducible auto-activation circuit is predicted to acquire low cooperativity and low fold-induction. Examples from Escherichia coli's two-component signaling systems support these predictions

Two Component Systems: Physiological Effect of a Third Component

Signal transduction systems mediate the response and adaptation of organisms to environmental changes. In prokaryotes, this signal transduction is often done through Two Component Systems (TCS). These TCS are phosphotransfer protein cascades, and in their prototypical form they are composed by a kinase that senses the environmental signals (SK) and by a response regulator (RR) that regulates the cellular response. This basic motif can be modified by the addition of a third protein that interacts either with the SK or the RR in a way that could change the dynamic response of the TCS module. In this work we aim at understanding the effect of such an additional protein (which we call “third component”) on the functional properties of a prototypical TCS. To do so we build mathematical models of TCS with alternative designs for their interaction with that third component. These mathematical models are analyzed in order to identify the differences in dynamic behavior inherent to each design, with respect to functionally relevant properties such as sensitivity to changes in either the parameter values or the molecular concentrations, temporal responsiveness, possibility of multiple steady states, or stochastic fluctuations in the system. The differences are then correlated to the physiological requirements that impinge on the functioning of the TCS. This analysis sheds light on both, the dynamic behavior of synthetically designed TCS, and the conditions under which natural selection might favor each of the designs. We find that a third component that modulates SK activity increases the parameter space where a bistable response of the TCS module to signals is possible, if SK is monofunctional, but decreases it when the SK is bifunctional. The presence of a third component that modulates RR activity decreases the parameter space where a bistable response of the TCS module to signals is possible

A Probabilistic Model of Local Sequence Alignment That Simplifies Statistical Significance Estimation

Sequence database searches require accurate estimation of the statistical significance of scores. Optimal local sequence alignment scores follow Gumbel distributions, but determining an important parameter of the distribution (λ) requires time-consuming computational simulation. Moreover, optimal alignment scores are less powerful than probabilistic scores that integrate over alignment uncertainty (“Forward” scores), but the expected distribution of Forward scores remains unknown. Here, I conjecture that both expected score distributions have simple, predictable forms when full probabilistic modeling methods are used. For a probabilistic model of local sequence alignment, optimal alignment bit scores (“Viterbi” scores) are Gumbel-distributed with constant λ = log 2, and the high scoring tail of Forward scores is exponential with the same constant λ. Simulation studies support these conjectures over a wide range of profile/sequence comparisons, using 9,318 profile-hidden Markov models from the Pfam database. This enables efficient and accurate determination of expectation values (E-values) for both Viterbi and Forward scores for probabilistic local alignments

Information Routing Driven by Background Chatter in a Signaling Network

Living systems are capable of processing multiple sources of information simultaneously. This is true even at the cellular level, where not only coexisting signals stimulate the cell, but also the presence of fluctuating conditions is significant. When information is received by a cell signaling network via one specific input, the existence of other stimuli can provide a background activity –or chatter– that may affect signal transmission through the network and, therefore, the response of the cell. Here we study the modulation of information processing by chatter in the signaling network of a human cell, specifically, in a Boolean model of the signal transduction network of a fibroblast. We observe that the level of external chatter shapes the response of the system to information carrying signals in a nontrivial manner, modulates the activity levels of the network outputs, and effectively determines the paths of information flow. Our results show that the interactions and node dynamics, far from being random, confer versatility to the signaling network and allow transitions between different information-processing scenarios

Attenuation of the Sensing Capabilities of PhoQ in Transition to Obligate Insect–Bacterial Association

Sodalis glossinidius, a maternally inherited endosymbiont of the tsetse fly, maintains genes encoding homologues of the PhoP-PhoQ two-component regulatory system. This two-component system has been extensively studied in facultative bacterial pathogens and is known to serve as an environmental magnesium sensor and a regulator of key virulence determinants. In the current study, we show that the inactivation of the response regulator, phoP, renders S. glossinidius sensitive to insect derived cationic antimicrobial peptides (AMPs). The resulting mutant strain displays reduced expression of genes involved in the structural modification of lipid A that facilitates resistance to AMPs. In addition, the inactivation of phoP alters the expression of type-III secretion system (TTSS) genes encoded within three distinct chromosomal regions, indicating that PhoP-PhoQ also serves as a master regulator of TTSS gene expression. In the absence of phoP, S. glossinidius is unable to superinfect either its natural tsetse fly host or a closely related hippoboscid louse fly. Furthermore, we show that the S. glossinidius PhoQ sensor kinase has undergone functional adaptations that result in a substantially diminished ability to sense ancestral signals. The loss of PhoQ's sensory capability is predicted to represent a novel adaptation to the static symbiotic lifestyle, allowing S. glossinidius to constitutively express genes that facilitate resistance to host derived AMPs