Search

In this study, we utilize Protein Residue Networks (PRNs), constructed using Local Spatial Pattern (LSP) alignment, to explore the dynamic behavior of Catabolite Activator Protein (CAP) upon the sequential binding of cAMP. We employed the Degree Centrality of these PRNs to investigate protein dynamics on a sub-nanosecond time scale, hypothesizing that it would reflect changes in CAPs entropy related to its thermal motions. We show that the binding of the first cAMP led to an increase in stability in the Cyclic-Nucleotide Binding Domain A (CNBD-A) and destabilization in CNBD-B, agreeing with previous reports explaining the negative cooperativity of cAMP binding in terms of an entropy-driven allostery. LSP-based PRNs also allow for the study of Betweenness Centrality, another graph-theoretical characteristic of PRNs, providing insights into global residue connectivity within CAP. Using this approach, we were able to correctly identify amino acids that were shown to be critical in mediating allosteric interactions in CAP. The agreement between our studies and previous experimental reports validates our method, particularly with respect to the reliability of Degree Centrality as a proxy for entropy related to protein thermal dynamics. Because LSP-based PRNs can be easily extended to include dynamics of small organic molecules, polynucleotides, or other allosteric proteins, the methods presented here mark a significant advancement in the field, positioning them as vital tools for a fast, cost-effective, and accurate analysis of entropy-driven allostery and identification of allosteric hotspots.

Many microorganisms exhibit nutrient preferences, exemplified by the 'hierarchical' consumption of certain carbon substrates. Here, we systematically investigate under which physiological conditions hierarchical substrate utilization occurs and its mechanisms of implementation. We show utilization hierarchy of Escherichia coli to be ordered by the carbon-uptake flux rather than the identity of the substrates. A detailed study of glycerol uptake finds that it is fully suppressed if the uptake flux of another glycolytic substrate exceeds a threshold, which is set to the influx obtained when grown on glycerol alone. Below this threshold, limited glycerol uptake is 'supplemented' such that the total carbon uptake is maintained at the threshold. This behaviour results from total-flux feedback mediated by cAMP-Crp signalling but also requires inhibition by the regulator fructose 1,6-bisphosphate, which senses the upper-glycolytic flux and ensures that glycerol uptake defers to other glycolytic substrates but not to gluconeogenic ones. A quantitative model reproduces all of the observed utilization patterns, including those of key mutants. The proposed mechanism relies on the differential regulation of uptake enzymes and requires a specific operon organization. This organization is found to be conserved across related species for several uptake systems, suggesting the deployment of similar mechanisms for hierarchical substrate utilization by a spectrum of microorganisms.

Ampk is an energy gatekeeper that responds to decreases in ATP by inhibiting energy-consuming anabolic processes and promoting energy-generating catabolic processes. Recently, we showed that Lkb1, an understudied kinase in B lymphocytes and a major upstream kinase for Ampk, had critical and unexpected roles in activating naïve B cells and in germinal center formation. Therefore, we examined whether Lkb1 activities during B cell activation depend on Ampk and report surprising Ampk activation with in vitro B cell stimulation in the absence of energy stress, coupled to rapid biomass accumulation. Despite Ampk activation and a controlling role for Lkb1 in B cell activation, Ampk knockout did not significantly affect B cell activation, differentiation, nutrient dynamics, gene expression, or humoral immune responses. Instead, Ampk loss specifically repressed the transcriptional expression of IgD and its regulator, Zfp318. Results also reveal that early activation of Ampk by phenformin treatment impairs germinal center formation but does not significantly alter antibody responses. Combined, the data show an unexpectedly specific role for Ampk in the regulation of IgD expression during B cell activation.

Ampk is an energy gatekeeper that responds to decreases in ATP by inhibiting energy-consuming anabolic processes and promoting energy-generating catabolic processes. Recently, we showed that Lkb1, an understudied kinase in B lymphocytes and a major upstream kinase for Ampk, had critical and unexpected roles in activating naïve B cells and in germinal center formation. Therefore, we examined whether Lkb1 activities during B cell activation depend on Ampk and report surprising Ampk activation with in vitro B cell stimulation in the absence of energy stress, coupled to rapid biomass accumulation. Despite Ampk activation and a controlling role for Lkb1 in B cell activation, Ampk knockout did not significantly affect B cell activation, differentiation, nutrient dynamics, gene expression, or humoral immune responses. Instead, Ampk loss specifically repressed the transcriptional expression of IgD and its regulator, Zfp318. Results also reveal that early activation of Ampk by phenformin treatment impairs germinal center formation but does not significantly alter antibody responses. Combined, the data show an unexpectedly specific role for Ampk in the regulation of IgD expression during B cell activation.

Two major transcriptional regulators of carbon metabolism in bacteria are Cra and CRP. CRP is considered to be the main mediator of catabolite repression. Unlike for CRP, in vivo DNA binding information of Cra is scarce. Here we generate and integrate ChIP-exo and RNA-seq data to identify 39 binding sites for Cra and 97 regulon genes that are regulated by Cra in Escherichia coli. An integrated metabolic-regulatory network was formed by including experimentally-derived regulatory information and a genome-scale metabolic network reconstruction. Applying analysis methods of systems biology to this integrated network showed that Cra enables optimal bacterial growth on poor carbon sources by redirecting and repressing glycolysis flux, by activating the glyoxylate shunt pathway, and by activating the respiratory pathway. In these regulatory mechanisms, the overriding regulatory activity of Cra over CRP is fundamental. Thus, elucidation of interacting transcriptional regulation of core carbon metabolism in bacteria by two key transcription factors was possible by combining genome-wide experimental measurement and simulation with a genome-scale metabolic model.

Two major transcriptional regulators of carbon metabolism in bacteria are Cra and CRP. CRP is considered to be the main mediator of catabolite repression. Unlike for CRP, in vivo DNA binding information of Cra is scarce. Here we generate and integrate ChIP-exo and RNA-seq data to identify 39 binding sites for Cra and 97 regulon genes that are regulated by Cra in Escherichia coli. An integrated metabolic-regulatory network was formed by including experimentally-derived regulatory information and a genome-scale metabolic network reconstruction. Applying analysis methods of systems biology to this integrated network showed that Cra enables optimal bacterial growth on poor carbon sources by redirecting and repressing glycolysis flux, by activating the glyoxylate shunt pathway, and by activating the respiratory pathway. In these regulatory mechanisms, the overriding regulatory activity of Cra over CRP is fundamental. Thus, elucidation of interacting transcriptional regulation of core carbon metabolism in bacteria by two key transcription factors was possible by combining genome-wide experimental measurement and simulation with a genome-scale metabolic model.

Numerous in vitro studies have yielded a refined picture of the structural and molecular associations between Cyclic-AMP receptor protein (Crp), the DNA motif, and RNA polymerase (RNAP) holoenzyme. In this study, high-resolution ChIP-exonuclease (ChIP-exo) was applied to study Crp binding in vivo and at genome-scale. Surprisingly, Crp was found to provide little to no protection of the DNA motif under activating conditions. Instead, Crp demonstrated binding patterns that closely resembled those generated by σ70. The binding patterns of both Crp and σ70 are indicative of RNAP holoenzyme DNA footprinting profiles associated with stages during transcription initiation that occur post-recruitment. This is marked by a pronounced advancement of the template strand footprint profile to the +20 position relative to the transcription start site and a multimodal distribution on the nontemplate strand. This trend was also observed in the familial transcription factor, Fnr, but full protection of the motif was seen in the repressor ArcA. Given the time-scale of ChIP studies and that the rate-limiting step in transcription initiation is typically post recruitment, we propose a hypothesis where Crp is absent from the DNA motif but remains associated with RNAP holoenzyme post-recruitment during transcription initiation. The release of Crp from the DNA motif may be a result of energetic changes that occur as RNAP holoenzyme traverses the various stable intermediates towards elongation complex formation.

Numerous in vitro studies have yielded a refined picture of the structural and molecular associations between Cyclic-AMP receptor protein (Crp), the DNA motif, and RNA polymerase (RNAP) holoenzyme. In this study, high-resolution ChIP-exonuclease (ChIP-exo) was applied to study Crp binding in vivo and at genome-scale. Surprisingly, Crp was found to provide little to no protection of the DNA motif under activating conditions. Instead, Crp demonstrated binding patterns that closely resembled those generated by σ70. The binding patterns of both Crp and σ70 are indicative of RNAP holoenzyme DNA footprinting profiles associated with stages during transcription initiation that occur post-recruitment. This is marked by a pronounced advancement of the template strand footprint profile to the +20 position relative to the transcription start site and a multimodal distribution on the nontemplate strand. This trend was also observed in the familial transcription factor, Fnr, but full protection of the motif was seen in the repressor ArcA. Given the time-scale of ChIP studies and that the rate-limiting step in transcription initiation is typically post recruitment, we propose a hypothesis where Crp is absent from the DNA motif but remains associated with RNAP holoenzyme post-recruitment during transcription initiation. The release of Crp from the DNA motif may be a result of energetic changes that occur as RNAP holoenzyme traverses the various stable intermediates towards elongation complex formation.

When bacteria are cultured in medium with multiple carbon substrates, they frequently consume these substrates simultaneously. Building on recent advances in the understanding of metabolic coordination exhibited by Escherichia coli cells through cAMP-Crp signaling, we show that this signaling system responds to the total carbon-uptake flux when substrates are co-utilized and derive a mathematical formula that accurately predicts the resulting growth rate, based only on the growth rates on individual substrates.

When bacteria are cultured in medium with multiple carbon substrates, they frequently consume these substrates simultaneously. Building on recent advances in the understanding of metabolic coordination exhibited by Escherichia coli cells through cAMP-Crp signaling, we show that this signaling system responds to the total carbon-uptake flux when substrates are co-utilized and derive a mathematical formula that accurately predicts the resulting growth rate, based only on the growth rates on individual substrates.

The phosphoenolpyruvate:sugar phosphotransferase system (PTS) has been shown to control transport, cell metabolism and gene expression. We here present results supporting the novel suggestion that in certain instances it also regulates the mutation rate. Directed mutations are defined as mutations that occur at higher frequencies when beneficial than when neutral or detrimental. To date, the occurrence of directed point mutations has not been documented and confirmed, but a few examples of transposon-mediated directed mutations have been reported. Here we focus on the first and best-studied example of directed mutation, which involves the regulation of insertion sequence-5 hopping into a specific site upstream of the glpFK glycerol utilization operon in Escherichia coli. This insertional event specifically activates expression of the glpFK operon, allowing the growth of wild-type cells with glycerol as a carbon source in the presence of nonmetabolizable glucose analogues which normally block glycerol utilization. The sugar-transporting PTS controls this process by regulating levels of cytoplasmic glycerol-3-phosphate and cyclic (c)AMP as established in previous publications. Direct involvement of the glycerol repressor, GlpR, and the cAMP receptor protein, Crp, in the regulation of transposon-mediated directed mutation has been demonstrated.

The phosphoenolpyruvate:sugar phosphotransferase system (PTS) has been shown to control transport, cell metabolism and gene expression. We here present results supporting the novel suggestion that in certain instances it also regulates the mutation rate. Directed mutations are defined as mutations that occur at higher frequencies when beneficial than when neutral or detrimental. To date, the occurrence of directed point mutations has not been documented and confirmed, but a few examples of transposon-mediated directed mutations have been reported. Here we focus on the first and best-studied example of directed mutation, which involves the regulation of insertion sequence-5 hopping into a specific site upstream of the glpFK glycerol utilization operon in Escherichia coli. This insertional event specifically activates expression of the glpFK operon, allowing the growth of wild-type cells with glycerol as a carbon source in the presence of nonmetabolizable glucose analogues which normally block glycerol utilization. The sugar-transporting PTS controls this process by regulating levels of cytoplasmic glycerol-3-phosphate and cyclic (c)AMP as established in previous publications. Direct involvement of the glycerol repressor, GlpR, and the cAMP receptor protein, Crp, in the regulation of transposon-mediated directed mutation has been demonstrated.

Growth of E. coli on several carbon sources is dependent on the catabolite repressor/activator (Cra) protein although a Cra consensus DNA-binding site is not present in the control regions of the relevant catabolic operons. We show that Cra regulates growth by activating expression of the crp gene. It thereby mediates catabolite repression of catabolic operons by an indirect mechanism.

Growth of E. coli on several carbon sources is dependent on the catabolite repressor/activator (Cra) protein although a Cra consensus DNA-binding site is not present in the control regions of the relevant catabolic operons. We show that Cra regulates growth by activating expression of the crp gene. It thereby mediates catabolite repression of catabolic operons by an indirect mechanism.

Growth of E. coli on several carbon sources is dependent on the catabolite repressor/activator (Cra) protein although a Cra consensus DNA-binding site is not present in the control regions of the relevant catabolic operons. We show that Cra regulates growth by activating expression of the crp gene. It thereby mediates catabolite repression of catabolic operons by an indirect mechanism.

Finding conserved motifs in genomic sequences represents one of essential bioinformatic problems. However, achieving high discovery performance without imposing substantial auxiliary constraints on possible motif features remains a key algorithmic challenge. This work describes BAMBI-a sequential Monte Carlo motif-identification algorithm, which is based on a position weight matrix model that does not require additional constraints and is able to estimate such motif properties as length, logo, number of instances and their locations solely on the basis of primary nucleotide sequence data. Furthermore, should biologically meaningful information about motif attributes be available, BAMBI takes advantage of this knowledge to further refine the discovery results. In practical applications, we show that the proposed approach can be used to find sites of such diverse DNA-binding molecules as the cAMP receptor protein (CRP) and Din-family site-specific serine recombinases. Results obtained by BAMBI in these and other settings demonstrate better statistical performance than any of the four widely-used profile-based motif discovery methods: MEME, BioProspector with BioOptimizer, SeSiMCMC and Motif Sampler as measured by the nucleotide-level correlation coefficient. Additionally, in the case of Din-family recombinase target site discovery, the BAMBI-inferred motif is found to be the only one functionally accurate from the underlying biochemical mechanism standpoint. C++ and Matlab code is available at http://www.ee.columbia.edu/~guido/BAMBI or http://genomics.lbl.gov/BAMBI/.

Finding conserved motifs in genomic sequences represents one of essential bioinformatic problems. However, achieving high discovery performance without imposing substantial auxiliary constraints on possible motif features remains a key algorithmic challenge. This work describes BAMBI-a sequential Monte Carlo motif-identification algorithm, which is based on a position weight matrix model that does not require additional constraints and is able to estimate such motif properties as length, logo, number of instances and their locations solely on the basis of primary nucleotide sequence data. Furthermore, should biologically meaningful information about motif attributes be available, BAMBI takes advantage of this knowledge to further refine the discovery results. In practical applications, we show that the proposed approach can be used to find sites of such diverse DNA-binding molecules as the cAMP receptor protein (CRP) and Din-family site-specific serine recombinases. Results obtained by BAMBI in these and other settings demonstrate better statistical performance than any of the four widely-used profile-based motif discovery methods: MEME, BioProspector with BioOptimizer, SeSiMCMC and Motif Sampler as measured by the nucleotide-level correlation coefficient. Additionally, in the case of Din-family recombinase target site discovery, the BAMBI-inferred motif is found to be the only one functionally accurate from the underlying biochemical mechanism standpoint. C++ and Matlab code is available at http://www.ee.columbia.edu/~guido/BAMBI or http://genomics.lbl.gov/BAMBI/.