Terence Hwa

On schemes of combinatorial transcription logic

Cells receive a wide variety of cellular and environmental signals, which are often processed combinatorially to generate specific genetic responses. Here we explore theoretically the potentials and limitations of combinatorial signal integration at the level of cis-regulatory transcription control. Our analysis suggests that many complex transcription-control functions of the type encountered in higher eukaryotes are already implementable within the much simpler bacterial transcription system. Using a quantitative model of bacterial transcription and invoking only specific protein–DNA interaction and weak glue-like interaction between regulatory proteins, we show explicit schemes to implement regulatory logic functions of increasing complexity by appropriately selecting the strengths and arranging the relative positions of the relevant protein-binding DNA sequences in the cis-regulatory region. The architectures that emerge are naturally modular and evolvable. Our results suggest that the transcription regulatory apparatus is a “programmable” computing machine, belonging formally to the class of Boltzmann machines. Crucial to our results is the ability to regulate gene expression at a distance. In bacteria, this can be achieved for isolated genes via DNA looping controlled by the dimerization of DNA-bound proteins. However, if adopted extensively in the genome, long-distance interaction can cause unintentional intergenic cross talk, a detrimental side effect difficult to overcome by the known bacterial transcription-regulation systems. This may be a key factor limiting the genome-wide adoption of complex transcription control in bacteria. Implications of our findings for combinatorial transcription control in eukaryotes are discussed.

Interdependence of cell growth and gene expression: origins and consequences.

Theory of Growth Control Although quantitative studies of growth in bacterial cultures have been made for over 50 years, the relationship between cell proliferation and gene expression has not been clear. Scott et al. (p. 1099; see the Perspective by Lerman and Palsson) have revealed that mass per cell exponentially increased with linear increases in growth rate and that ribosome abundance increased linearly with growth rate depending on the rate of translation. Hence, the systems properties of the biological processes involved in growth can be derived without any molecular understanding of their basis and can be used to establish fundamental properties for the design of biotechnological procedures. Simple mathematical models describe the relationship between bacterial replication, cellular resources, and protein expression. In bacteria, the rate of cell proliferation and the level of gene expression are intimately intertwined. Elucidating these relations is important both for understanding the physiological functions of endogenous genetic circuits and for designing robust synthetic systems. We describe a phenomenological study that reveals intrinsic constraints governing the allocation of resources toward protein synthesis and other aspects of cell growth. A theory incorporating these constraints can accurately predict how cell proliferation and gene expression affect one another, quantitatively accounting for the effect of translation-inhibiting antibiotics on gene expression and the effect of gratuitous protein expression on cell growth. The use of such empirical relations, analogous to phenomenological laws, may facilitate our understanding and manipulation of complex biological systems before underlying regulatory circuits are elucidated.

Direct-coupling analysis of residue coevolution captures native contacts across many protein families

The similarity in the three-dimensional structures of homologous proteins imposes strong constraints on their sequence variability. It has long been suggested that the resulting correlations among amino acid compositions at different sequence positions can be exploited to infer spatial contacts within the tertiary protein structure. Crucial to this inference is the ability to disentangle direct and indirect correlations, as accomplished by the recently introduced direct-coupling analysis (DCA). Here we develop a computationally efficient implementation of DCA, which allows us to evaluate the accuracy of contact prediction by DCA for a large number of protein domains, based purely on sequence information. DCA is shown to yield a large number of correctly predicted contacts, recapitulating the global structure of the contact map for the majority of the protein domains examined. Furthermore, our analysis captures clear signals beyond intradomain residue contacts, arising, e.g., from alternative protein conformations, ligand-mediated residue couplings, and interdomain interactions in protein oligomers. Our findings suggest that contacts predicted by DCA can be used as a reliable guide to facilitate computational predictions of alternative protein conformations, protein complex formation, and even the de novo prediction of protein domain structures, contingent on the existence of a large number of homologous sequences which are being rapidly made available due to advances in genome sequencing.

Identification of direct residue contacts in protein–protein interaction by message passing

Understanding the molecular determinants of specificity in protein–protein interaction is an outstanding challenge of postgenome biology. The availability of large protein databases generated from sequences of hundreds of bacterial genomes enables various statistical approaches to this problem. In this context covariance-based methods have been used to identify correlation between amino acid positions in interacting proteins. However, these methods have an important shortcoming, in that they cannot distinguish between directly and indirectly correlated residues. We developed a method that combines covariance analysis with global inference analysis, adopted from use in statistical physics. Applied to a set of >2,500 representatives of the bacterial two-component signal transduction system, the combination of covariance with global inference successfully and robustly identified residue pairs that are proximal in space without resorting to ad hoc tuning parameters, both for heterointeractions between sensor kinase (SK) and response regulator (RR) proteins and for homointeractions between RR proteins. The spectacular success of this approach illustrates the effectiveness of the global inference approach in identifying direct interaction based on sequence information alone. We expect this method to be applicable soon to interaction surfaces between proteins present in only 1 copy per genome as the number of sequenced genomes continues to expand. Use of this method could significantly increase the potential targets for therapeutic intervention, shed light on the mechanism of protein–protein interaction, and establish the foundation for the accurate prediction of interacting protein partners.

Growth Rate-Dependent Global Effects on Gene Expression in Bacteria

Bacterial gene expression depends not only on specific regulatory mechanisms, but also on bacterial growth, because important global parameters such as the abundance of RNA polymerases and ribosomes are all growth-rate dependent. Understanding of these global effects is necessary for a quantitative understanding of gene regulation and for the design of synthetic genetic circuits. We find that the observed growth-rate dependence of constitutive gene expression can be explained by a simple model using the measured growth-rate dependence of the relevant cellular parameters. More complex growth dependencies for genetic circuits involving activators, repressors, and feedback control were analyzed and verified experimentally with synthetic circuits. Additional results suggest a feedback mechanism mediated by general growth-dependent effects that does not require explicit gene regulation if the expressed protein affects cell growth. This mechanism can lead to growth bistability and promote the acquisition of important physiological functions such as antibiotic resistance and tolerance (persistence).

Physical constraints and functional characteristics of transcription factor–DNA interaction

We study theoretical “design principles” for transcription factor (TF)–DNA interaction in bacteria, focusing particularly on the statistical interaction of the TFs with the genomic background (i.e., the genome without the target sites). We introduce and motivate the concept of programmability, i.e., the ability to set the threshold concentration for TF binding over a wide range merely by mutating the binding sequence of a target site. This functional demand, together with physical constraints arising from the thermodynamics and kinetics of TF–DNA interaction, leads us to a narrow range of “optimal” interaction parameters. We find that this parameter set agrees well with experimental data for the interaction parameters of a few exemplary prokaryotic TFs, which indicates that TF–DNA interaction is indeed programmable. We suggest further experiments to test whether this is a general feature for a large class of TFs.

Combinatorial transcriptional control of the lactose operon of Escherichia coli

The goal of systems biology is to understand the behavior of the whole in terms of knowledge of the parts. This is hard to achieve in many cases due to the difficulty of characterizing the many constituents involved in a biological system and their complex web of interactions. The lac promoter of Escherichia coli offers the possibility of confronting “system-level” properties of transcriptional regulation with the known biochemistry of the molecular constituents and their mutual interactions. Such confrontations can reveal previously unknown constituents and interactions, as well as offer insight into how the components work together as a whole. Here we study the combinatorial control of the lac promoter by the regulators Lac repressor (LacR) and cAMP-receptor protein (CRP). A previous in vivo study [Setty Y, Mayo AE, Surette MG, Alon U (2003) Proc Natl Acad Sci USA 100:7702–7707] found gross disagreement between the observed promoter activities and the expected behavior based on the known molecular mechanisms. We repeated the study by identifying and removing several extraneous factors that significantly modulated the expression of the lac promoter. Through quantitative, systematic characterization of promoter activity for a number of key mutants and guided by the thermodynamic model of transcriptional regulation, we were able to account for the combinatorial control of the lac promoter quantitatively, in terms of a cooperative interaction between CRP and LacR-mediated DNA looping. Specifically, our analysis indicates that the sensitivity of the inducer response results from LacR-mediated DNA looping, which is significantly enhanced by CRP.

Coordination of bacterial proteome with metabolism by cyclic AMP signalling

The cyclic AMP (cAMP)-dependent catabolite repression effect in Escherichia coli is among the most intensely studied regulatory processes in biology. However, the physiological function(s) of cAMP signalling and its molecular triggers remain elusive. Here we use a quantitative physiological approach to show that cAMP signalling tightly coordinates the expression of catabolic proteins with biosynthetic and ribosomal proteins, in accordance with the cellular metabolic needs during exponential growth. The expression of carbon catabolic genes increased linearly with decreasing growth rates upon limitation of carbon influx, but decreased linearly with decreasing growth rate upon limitation of nitrogen or sulphur influx. In contrast, the expression of biosynthetic genes showed the opposite linear growth-rate dependence as the catabolic genes. A coarse-grained mathematical model provides a quantitative framework for understanding and predicting gene expression responses to catabolic and anabolic limitations. A scheme of integral feedback control featuring the inhibition of cAMP signalling by metabolic precursors is proposed and validated. These results reveal a key physiological role of cAMP-dependent catabolite repression: to ensure that proteomic resources are spent on distinct metabolic sectors as needed in different nutrient environments. Our findings underscore the power of quantitative physiology in unravelling the underlying functions of complex molecular signalling networks.

Sequential Establishment of Stripe Patterns in an Expanding Cell Population

A synthetic circuit implementing density-controlled bacterial motility autonomously produces a tunable stripe pattern. Periodic stripe patterns are ubiquitous in living organisms, yet the underlying developmental processes are complex and difficult to disentangle. We describe a synthetic genetic circuit that couples cell density and motility. This system enabled programmed Escherichia coli cells to form periodic stripes of high and low cell densities sequentially and autonomously. Theoretical and experimental analyses reveal that the spatial structure arises from a recurrent aggregation process at the front of the continuously expanding cell population. The number of stripes formed could be tuned by modulating the basal expression of a single gene. The results establish motility control as a simple route to establishing recurrent structures without requiring an extrinsic pacemaker.

Genomics-aided structure prediction

We introduce a theoretical framework that exploits the ever-increasing genomic sequence information for protein structure prediction. Structure-based models are modified to incorporate constraints by a large number of non-local contacts estimated from direct coupling analysis (DCA) of co-evolving genomic sequences. A simple hybrid method, called DCA-fold, integrating DCA contacts with an accurate knowledge of local information (e.g., the local secondary structure) is sufficient to fold proteins in the range of 1–3 Å resolution.