DNA looping and unlooping by AraC protein

Identification, characterization and classification of prokaryotic nucleoid-associated proteins

Common throughout life is the need to compact and organize the genome. Possible mechanisms involved in this process include supercoiling, phase separation, charge neutralization, macromolecular crowding, and nucleoid-associated proteins (NAPs). NAPs are special in that they can organize the genome at multiple length scales, and thus are often considered as the architects of the genome. NAPs shape the genome by either bending DNA, wrapping DNA, bridging DNA, or forming nucleoprotein filaments on the DNA. In this mini-review, we discuss recent advancements of unique NAPs with differing architectural properties across the tree of life, including NAPs from bacteria, archaea, and viruses. To help the characterization of NAPs from the ever-increasing number of metagenomes, we recommend a set of cheap and simple in vitro biochemical assays that give unambiguous insights into the architectural properties of NAPs. Finally, we highlight and showcase the usefulness of AlphaFold in the characterization of novel NAPs.

Atomic force microscopy measurements and model of DNA bending caused by binding of AraC protein

Atomic force microscopy (AFM) was used to conduct single-molecule imaging of protein/DNA complexes involved in the regulation of the arabinose operon of Escherichia coli. In the presence of arabinose, the transcription regulatory protein AraC binds to a 38 bp region consisting of the araI1 and araI2 half-sites. The domain positioning of full-length AraC, when bound to DNA, was not previously known. In this study, AraC was combined with 302 and 560 bp DNA and arabinose, deposited on a mica substrate, and imaged with AFM in air. High resolution images of 560 bp DNA, where bound protein was visible, showed that AraC induces a bend in the DNA with an angle 60° ± 12° with a median of 55°. These results are consistent with earlier gel electrophoresis measurements that measured the DNA bend angle based on migration rates. By using known domain structures of AraC, geometric constraints, and contacts determined from biochemical experiments, we developed a model of the tertiary and quaternary structure of DNA-bound AraC in the presence of arabinose. The DNA bend angle predicted by the model is in agreement with the measurement values. We discuss the results in view of other regulatory proteins that cause DNA bending and formation of the open complex to initiate transcription.

A Career's Work, the l-Arabinose Operon: How It Functions and How We Learned It

Very few labs have had the good fortune to have been able to focus for more than 50 years on a relatively narrow research topic and to be in a field in which both basic knowledge and the research technology and methods have progressed as rapidly as they have in molecular biology. My research group, first at Brandeis University and then at Johns Hopkins University, has had this opportunity. In this review, therefore, I will describe largely the work from my laboratory that has spanned this period and which was carried out by 40 plus graduate students, several postdoctoral associates, my technician, and me. In addition to presenting the scientific findings or results, I will place many of the topics in scientific context and, because we needed to develop a good many of the experimental methods behind our findings, I will also describe some of these methods and their importance. Also included will be occasional comments on how the research community or my research group functioned. Because a wide variety of approaches were used throughout our work, no ideal organization of this review is apparent. Therefore, I have chosen to use a hybrid structure in which there are six sections. Within each of the sections, experiments and findings will be described roughly in chronological order. Frequent cross references between parts and sections will be made because some findings and experimental approaches could logically have been described in more than one place.

TADfit is a multivariate linear regression model for profiling hierarchical chromatin domains on replicate Hi-C data

Topologically associating domains (TADs) are fundamental building blocks of three dimensional genome, and organized into complex hierarchies. Identifying hierarchical TADs on Hi-C data helps to understand the relationship between genome architectures and gene regulation. Herein we propose TADfit, a multivariate linear regression model for profiling hierarchical chromatin domains, which tries to fit the interaction frequencies in Hi-C contact matrix with and without replicates using all-possible hierarchical TADs, and the significant ones can be determined by the regression coefficients obtained with the help of an online learning solver called Follow-The-Regularized-Leader (FTRL). Beyond the existing methods, TADfit has an ability to handle multiple contact matrix replicates and find partially overlapping TADs on them, which helps to find the comprehensive underlying TADs across replicates from different experiments. The comparative results tell that TADfit has better accuracy and reproducibility, and the hierarchical TADs called by it exhibit a reasonable biological relevance.

The AraC/XylS Protein MxiE and Its Coregulator IpgC Control a Negative Feedback Loop in the Transcriptional Cascade That Regulates Type III Secretion in Shigella flexneri

Members of the AraC family of transcriptional regulators (AFTRs) control the expression of many genes important to cellular processes, including virulence. In Shigella species, the type III secretion system (T3SS), a key determinant for host cell invasion, is regulated by the three-tiered VirF/VirB/MxiE transcriptional cascade. Both VirF and MxiE belong to the AFTRs and are characterized as positive transcriptional regulators. Here, we identify a novel regulatory activity for MxiE and its coregulator IpgC, which manifests as a negative feedback loop in the VirF/VirB/MxiE transcriptional cascade. Our findings show that MxiE and IpgC downregulate the virB promoter and, hence, VirB protein production, thus decreasing VirB-dependent promoter activity at ospD1, one of the nearly 50 VirB-dependent genes. At the virB promoter, regions required for negative MxiE- and IpgC-dependent regulation were mapped and found to be coincident with regions required for positive VirF-dependent regulation. In tandem, negative MxiE- and IpgC-dependent regulation of the virB promoter only occurred in the presence of VirF, suggesting that MxiE and IpgC can function to counter VirF activation of the virB promoter. Lastly, MxiE and IpgC do not downregulate another VirF-activated promoter, icsA, demonstrating that this negative feedback loop targets the virB promoter. Our study provides insight into a mechanism that may reprogram Shigella virulence gene expression following type III secretion and provides the impetus to examine if MxiE and IpgC homologs in other important bacterial pathogens, such as Burkholderia pseudomallei and Salmonella enterica serovars Typhimurium and Typhi, coordinate similar negative feedback loops. IMPORTANCE The large AraC family of transcriptional regulators (AFTRs) control virulence gene expression in many bacterial pathogens. In Shigella species, the AraC/XylS protein MxiE and its coregulator IpgC positively regulate the expression of type III secretion system genes within the three-tiered VirF/VirB/MxiE transcriptional cascade. Our findings suggest a negative feedback loop in the VirF/VirB/MxiE cascade, in which MxiE and IpgC counter VirF-dependent activation of the virB promoter, thus making this the first characterization of negative MxiE- and IpgC-dependent regulation. Our study provides insight into a mechanism that likely reprograms Shigella virulence gene expression following type III secretion, which has implications for other important bacterial pathogens with functional homologs of MxiE and IpgC.

HIV-1 VIF and human APOBEC3G interaction directly observed through molecular specific labeling using a new dual promotor vector

Over the last few decades, protein NMR isotope labeling methods using E. coli based expression have revolutionized the information accessible from biomolecular NMR experiments. Selective labeling of a protein of interest in a multi-protein complex can significantly reduce the number of cross-peaks and allow for study of large protein complexes. However, limitations still remain since some proteins are not stable independently and cannot be separately labeled in either NMR active isotope enriched or unenriched media and reconstituted into a multimeric complex. To overcome this limitation, the LEGO NMR method was previously developed using protein expression plasmids containing T7 or araBAD promoters to separately express proteins in the same E. coli after changing between labeled and unlabeled media. Building on this, we developed a method to label the Human Immunodeficiency Virus type 1 viral infectivity factor (HIV-1 Vif), a monomerically unstable protein, in complex with CBFβ, it's host binding partner. We designed a dual promoter plasmid containing both T7 and araBAD promoters to independently control the expression of HIV-1 Vif in NMR active isotope enriched media and CBFβ in unenriched media. Using this method, we assigned the backbone resonance and directly observed the binding of HIV-1 Vif with APOBEC3G, a host restriction factor to HIV-1.

Co-evolution of drug resistance and broadened substrate recognition in HIV protease variants isolated from an Escherichia coli genetic selection system

A genetic selection system for activity of HIV protease is described that is based on a synthetic substrate constructed as a modified AraC regulatory protein that when cleaved stimulate L-arabinose metabolism in an Escherichia coli araC strain. Growth stimulation on selective plates was shown to depend on active HIV protease and the scissile bond in the substrate. In addition, the growth of cells correlated well with the established cleavage efficiency of the sites in the viral polyprotein, Gag, when these sites were individually introduced into the synthetic substate of the selection system. Plasmids encoding protease variants selected based on stimulation of cell growth in the presence of saquinavir or cleavage of a site not cleaved by wild-type protease, were indistinguishable with respect to both phenotypes. Also, both groups of selected plasmids encoded side chain substitutions known from clinical isolates or displayed different side chain substitutions but at identical positions. One highly frequent side chain substitution, E34V, not regarded as a major drug resistance substitution was found in variants obtained under both selective conditions and is suggested to improve protease processing of the synthetic substrate. This substitution is away from the substrate-binding cavity and together with other substitutions in the selected reading frames supports the previous suggestion of a substrate-binding site extended from the active site binding pocket itself.

Quantifying the regulatory role of individual transcription factors in Escherichia coli

Gene regulation often results from the action of multiple transcription factors (TFs) acting at a promoter, obscuring the individual regulatory effect of each TF on RNA polymerase (RNAP). Here we measure the fundamental regulatory interactions of TFs in E. coli by designing synthetic target genes that isolate individual TFs' regulatory effects. Using a thermodynamic model, each TF's regulatory interactions are decoupled from TF occupancy and interpreted as acting through (de)stabilization of RNAP and (de)acceleration of transcription initiation. We find that the contribution of each mechanism depends on TF identity and binding location; regulation immediately downstream of the promoter is insensitive to TF identity, but the same TFs regulate by distinct mechanisms upstream of the promoter. These two mechanisms are uncoupled and can act coherently, to reinforce the observed regulatory role (activation/repression), or incoherently, wherein the TF regulates two distinct steps with opposing effects.

Structure of the master regulator Rns reveals an inhibitor of enterotoxigenic Escherichia coli virulence regulons

Enteric infections caused by the gram-negative bacteria enterotoxigenic Escherichia coli (ETEC), Vibrio cholerae, Shigella flexneri, and Salmonella enterica are among the most common and affect billions of people each year. These bacteria control expression of virulence factors using a network of transcriptional regulators, some of which are modulated by small molecules as has been shown for ToxT, an AraC family member from V. cholerae. In ETEC the expression of many types of adhesive pili is dependent upon the AraC family member Rns. We present here the 3 Å crystal structure of Rns and show it closely resembles ToxT. Rns crystallized as a dimer via an interface similar to that observed in other dimeric AraC's. Furthermore, the structure of Rns revealed the presence of a ligand, decanoic acid, that inhibits its activity in a manner similar to the fatty acid mediated inhibition observed for ToxT and the S. enterica homologue HilD. Together, these results support our hypothesis that fatty acids regulate virulence controlling AraC family members in a common manner across a number of enteric pathogens. Furthermore, for the first time this work identifies a small molecule capable of inhibiting the ETEC Rns regulon, providing a basis for development of therapeutics against this deadly human pathogen.

OnfD, an AraC-Type Transcriptional Regulator Encoded by Rhizobium tropici CIAT 899 and Involved in Nod Factor Synthesis and Symbiosis

Rhizobium tropici CIAT 899 is a broad-host-range rhizobial strain that establishes symbiotic interactions with legumes and tolerates different environmental stresses such as heat, acidity, or salinity. This rhizobial strain produces a wide variety of symbiotically active nodulation factors (NF) induced not only by the presence of plant-released flavonoids but also under osmotic stress conditions through the LysR-type transcriptional regulators NodD1 (flavonoids) and NodD2 (osmotic stress). However, the activation of NodD2 under high-osmotic-stress conditions remains elusive. Here, we have studied the role of a new AraC-type regulator (named as OnfD) in the symbiotic interaction of R. tropici CIAT 899 with Phaseolus vulgaris and Lotus plants. We determined that OnfD is required under salt stress conditions for the transcriptional activation of the nodulation genes and therefore the synthesis and export of NF, which are required for a successful symbiosis with P. vulgaris Moreover, using bacterial two-hybrid analysis, we demonstrated that the OnfD and NodD2 proteins form homodimers and OnfD/NodD2 form heterodimers, which could be involved in the production of NF in the presence of osmotic stress conditions since both regulators are required for NF synthesis in the presence of salt. A structural model of OnfD is presented and discussed.IMPORTANCE The synthesis and export of rhizobial NF are mediated by a conserved group of LysR-type regulators, the NodD proteins. Here, we have demonstrated that a non-LysR-type regulator, an AraC-type protein, is required for the transcriptional activation of symbiotic genes and for the synthesis of symbiotically active NF under salt stress conditions.

Evolution of Regulated Transcription

The genomes of all organisms abound with various cis-regulatory elements, which control gene activity. Transcriptional enhancers are a key group of such elements in eukaryotes and are DNA regions that form physical contacts with gene promoters and precisely orchestrate gene expression programs. Here, we follow gradual evolution of this regulatory system and discuss its features in different organisms. In eubacteria, an enhancer-like element is often a single regulatory element, is usually proximal to the core promoter, and is occupied by one or a few activators. Activation of gene expression in archaea is accompanied by the recruitment of an activator to several enhancer-like sites in the upstream promoter region. In eukaryotes, activation of expression is accompanied by the recruitment of activators to multiple enhancers, which may be distant from the core promoter, and the activators act through coactivators. The role of the general DNA architecture in transcription control increases in evolution. As a whole, it can be seen that enhancers of multicellular eukaryotes evolved from the corresponding prototypic enhancer-like regulatory elements with the gradually increasing genome size of organisms.

Synthetic circuits reveal how mechanisms of gene regulatory networks constrain evolution

Phenotypic variation is the raw material of adaptive Darwinian evolution. The phenotypic variation found in organismal development is biased towards certain phenotypes, but the molecular mechanisms behind such biases are still poorly understood. Gene regulatory networks have been proposed as one cause of constrained phenotypic variation. However, most pertinent evidence is theoretical rather than experimental. Here, we study evolutionary biases in two synthetic gene regulatory circuits expressed in Escherichia coli that produce a gene expression stripe-a pivotal pattern in embryonic development. The two parental circuits produce the same phenotype, but create it through different regulatory mechanisms. We show that mutations cause distinct novel phenotypes in the two networks and use a combination of experimental measurements, mathematical modelling and DNA sequencing to understand why mutations bring forth only some but not other novel gene expression phenotypes. Our results reveal that the regulatory mechanisms of networks restrict the possible phenotypic variation upon mutation. Consequently, seemingly equivalent networks can indeed be distinct in how they constrain the outcome of further evolution.

A novel DMAPP-responding genetic circuit sensor for high-throughput screening and evolving isoprene synthase

High-throughput screening is a popular tool for collating biological data which would otherwise require the use of excessive resources. In this study, an artificial genetic circuit sensor responding to dimethylallyl diphosphate (DMAPP) was constructed based on a modified L-arabinose operon for high-throughput screening and isoprene synthase (ispS) evolution in Escherichia coli (E. coli). As a first step, the DNA sequence of the L-arabinose ligand-binding domain (LBD) was replaced with an ispS gene to enable the AraC operon responding to DMAPP, which is the substrate of the IspS enzyme. Then, an enhanced GFP (eGFP) was also introduced as a reporter for pBAD promoter. The expression level of the reporter was monitored using either of the two tools: flow cytometer (FCM) and microplate reader. Sequentially, we observed that a high DMAPP concentration led to low eGFP fluorescence, and the overexpression of ispS gene, which consumes DMAPP, resulted in a high eGFP expression. These results demonstrated that the artificial genetic circuit sensor responded directly to the intracellular concentration of DMAPP, and the expression of IspS enzyme could be positively correlated to the expression level of eGFP. Finally, we identified two IspS mutants with different activities from an ispS gene library and further validated the screening method.

Protein-induced DNA linking number change by sequence-specific DNA binding proteins and its biological effects

Sequence-specific DNA-binding proteins play essential roles in many fundamental biological events such as DNA replication, recombination, and transcription. One common feature of sequence-specific DNA-binding proteins is to introduce structural changes to their DNA recognition sites including DNA-bending and DNA linking number change (ΔLk). In this article, I review recent progress in studying protein-induced ΔLk by several sequence-specific DNA-binding proteins, such as E. coli cAMP receptor protein (CRP) and lactose repressor (LacI). It was demonstrated recently that protein-induced ΔLk is an intrinsic property for sequence-specific DNA-binding proteins and does not correlate to protein-induced other structural changes, such as DNA bending. For instance, although CRP bends its DNA recognition site by 90°, it was not able to introduce a ΔLk to it. However, LacI was able to simultaneously bend and introduce a ΔLk to its DNA binding sites. Intriguingly, LacI also constrained superhelicity within LacI-lac O1 complexes if (-) supercoiled DNA templates were provided. I also discuss how protein-induced ΔLk help sequence-specific DNA-binding proteins regulate their biological functions. For example, it was shown recently that LacI utilizes the constrained superhelicity (ΔLk) in LacI-lac O1 complexes and serves as a topological barrier to constrain free, unconstrained (-) supercoils within the 401-bp DNA loop. These constrained (-) supercoils enhance LacI's binding affinity and therefore the repression of the lac promoter. Other biological functions include how DNA replication initiators λ O and DnaA use the induced ΔLk to open/melt bacterial DNA replication origins.

Physical, chemical, and metabolic state sensors expand the synthetic biology toolbox for Synechocystis sp. PCC 6803

Many under-developed organisms possess important traits that can boost the effectiveness and sustainability of microbial biotechnology. Photoautotrophic cyanobacteria can utilize the energy captured from light to fix carbon dioxide for their metabolic needs while living in environments not suited for growing crops. Various value-added compounds have been produced by cyanobacteria in the laboratory; yet, the products' titers and yields are often not industrially relevant and lag behind what have been accomplished in heterotrophic microbes. Genetic tools for biological process control are needed to take advantage of cyanobacteria's beneficial qualities, as tool development also lags behind what has been created in common heterotrophic hosts. To address this problem, we developed a suite of sensors that regulate transcription in the model cyanobacterium Synechocystis sp. PCC 6803 in response to metabolically relevant signals, including light and the cell's nitrogen status, and a family of sensors that respond to the inexpensive chemical, l-arabinose. Increasing the number of available tools enables more complex and precise control of gene expression. Expanding the synthetic biology toolbox for this cyanobacterium also improves our ability to utilize this important under-developed organism in biotechnology. Biotechnol. Bioeng. 2017;114: 1561-1569. © 2017 Wiley Periodicals, Inc.

Pseudomonas aeruginosa AmrZ Binds to Four Sites in the algD Promoter, Inducing DNA-AmrZ Complex Formation and Transcriptional Activation

During late stages of cystic fibrosis pulmonary infections, Pseudomonas aeruginosa often overproduces the exopolysaccharide alginate, protecting the bacterial community from host immunity and antimicrobials. The transcription of the alginate biosynthesis operon is under tight control by a number of factors, including AmrZ, the focus of this study. Interestingly, multiple transcription factors interact with the far-upstream region of this promoter (PalgD), in which one AmrZ binding site has been identified previously. The mechanisms of AmrZ binding and subsequent activation remain unclear and require more-detailed investigation. In this study, in-depth examinations elucidated four AmrZ binding sites, and their disruption eliminated AmrZ binding and promoter activation. Furthermore, our in vitro fluorescence resonance energy transfer experiments suggest that AmrZ holds together multiple binding sites in PalgD and thereafter induces the formation of higher-order DNA-AmrZ complexes. To determine the importance of interactions between those AmrZ oligomers in the cell, a DNA phasing experiment was performed. PalgD transcription was significantly impaired when the relative phase between AmrZ binding sites was reversed (5 bp), while a full-DNA-turn insertion (10 bp) restored promoter activity. Taken together, the investigations presented here provide a deeper mechanistic understanding of AmrZ-mediated binding to PalgD IMPORTANCE: Overproduction of the exopolysaccharide alginate provides protection to Pseudomonas aeruginosa against antimicrobial treatments and is associated with chronic P. aeruginosa infections in the lungs of cystic fibrosis patients. In this study, we combined a variety of microbiological, genetic, biochemical, and biophysical approaches to investigate the activation of the alginate biosynthesis operon promoter by a key transcription factor named AmrZ. This study has provided important new information on the mechanism of activation of this extremely complex promoter.

Protein-induced DNA linking number change by sequence-specific DNA binding proteins and its biological effects

Sequence-specific DNA-binding proteins play essential roles in many fundamental biological events such as DNA replication, recombination, and transcription. One common feature of sequence-specific DNA-binding proteins is to introduce structural changes to their DNA recognition sites including DNA-bending and DNA linking number change (ΔLk). In this article, I review recent progress in studying protein-induced ΔLk by several sequence-specific DNA-binding proteins, such as E. coli cAMP receptor protein (CRP) and lactose repressor (LacI). It was demonstrated recently that protein-induced ΔLk is an intrinsic property for sequence-specific DNA-binding proteins and does not correlate to protein-induced other structural changes, such as DNA bending. For instance, although CRP bends its DNA recognition site by 90°, it was not able to introduce a ΔLk to it. However, LacI was able to simultaneously bend and introduce a ΔLk to its DNA binding sites. Intriguingly, LacI also constrained superhelicity within LacI-lac O1 complexes if (-) supercoiled DNA templates were provided. I also discuss how protein-induced ΔLk help sequence-specific DNA-binding proteins regulate their biological functions. For example, it was shown recently that LacI utilizes the constrained superhelicity (ΔLk) in LacI-lac O1 complexes and serves as a topological barrier to constrain free, unconstrained (-) supercoils within the 401-bp DNA loop. These constrained (-) supercoils enhance LacI's binding affinity and therefore the repression of the lac promoter. Other biological functions include how DNA replication initiators λ O and DnaA use the induced ΔLk to open/melt bacterial DNA replication origins.

A large family of anti-activators accompanying XylS/AraC family regulatory proteins

AraC Negative Regulators (ANR) suppress virulence genes by directly down‐regulating AraC/XylS members in Gram‐negative bacteria. In this study, we sought to investigate the distribution and molecular mechanisms of regulatory function for ANRs among different bacterial pathogens. We identified more than 200 ANRs distributed in diverse clinically important gram negative pathogens, including Vibrio spp., Salmonella spp., Shigella spp., Yersinia spp., Citrobacter spp., enterotoxigenic (ETEC) and enteroaggregative E. coli (EAEC), and members of the Pasteurellaceae. By employing a bacterial two hybrid system, pull down assays and surface plasmon resonance (SPR) analysis, we demonstrate that Aar (AggR‐activated regulator), a prototype member of the ANR family in EAEC, binds with high affinity to the central linker domain of AraC‐like member AggR. ANR‐AggR binding disrupted AggR dimerization and prevented AggR‐DNA binding. ANR homologs of Vibrio cholerae, Citrobacter rodentium, Salmonella enterica and ETEC were capable of complementing Aar activity by repressing aggR expression in EAEC strain 042. ANR homologs of ETEC and Vibrio cholerae bound to AggR as well as to other members of the AraC family, including Rns and ToxT. The predicted proteins of all ANR members exhibit three highly conserved predicted α‐helices. Site‐directed mutagenesis studies suggest that at least predicted α‐helices 2 and 3 are required for Aar activity. In sum, our data strongly suggest that members of the novel ANR family act by directly binding to their cognate AraC partners.

Mutual regulation causes co-entrainment between a synthetic oscillator and the bacterial cell cycle

The correct functioning of cells requires the orchestration of multiple cellular processes, many of which are inherently dynamical. The conditions under which these dynamical processes entrain each other remain unclear. Here we use synthetic biology to address this question in the case of concurrent cellular oscillations. Specifically, we study at the single-cell level the interaction between the cell division cycle and a robust synthetic gene oscillator in Escherichia coli. Our results suggest that cell division is able to partially entrain the synthetic oscillations under normal growth conditions, by driving the periodic replication of the genes involved in the oscillator. Coupling the synthetic oscillations back into the cell cycle via the expression of a key regulator of chromosome replication increases the synchronization between the two periodic processes. A simple computational model allows us to confirm this effect.

Epistatic Interactions in the Arabinose Cis-Regulatory Element

Changes in gene expression are an important mode of evolution; however, the proximate mechanism of these changes is poorly understood. In particular, little is known about the effects of mutations within cis binding sites for transcription factors, or the nature of epistatic interactions between these mutations. Here, we tested the effects of single and double mutants in two cis binding sites involved in the transcriptional regulation of the Escherichia coli araBAD operon, a component of arabinose metabolism, using a synthetic system. This system decouples transcriptional control from any posttranslational effects on fitness, allowing a precise estimate of the effect of single and double mutations, and hence epistasis, on gene expression. We found that epistatic interactions between mutations in the araBAD cis-regulatory element are common, and that the predominant form of epistasis is negative. The magnitude of the interactions depended on whether the mutations are located in the same or in different operator sites. Importantly, these epistatic interactions were dependent on the presence of arabinose, a native inducer of the araBAD operon in vivo, with some interactions changing in sign (e.g., from negative to positive) in its presence. This study thus reveals that mutations in even relatively simple cis-regulatory elements interact in complex ways such that selection on the level of gene expression in one environment might perturb regulation in the other environment in an unpredictable and uncorrelated manner.

Impact of Residual Inducer on Titratable Expression Systems

Inducible expression systems are widely employed for the titratable control of gene expression, yet molecules inadvertently present in the growth medium or synthesized by the host cells can alter the response profile of some of these systems. Here, we explored the quantitative impact of these residual inducers on the apparent response properties of inducible systems. Using a simple mathematical model, we found that the presence of residual inducer shrinks the apparent dynamic range and causes the apparent Hill coefficient to converge to one. We also found that activating systems were more sensitive than repressing systems to the presence of residual inducer and the response parameters were most heavily dependent on the original Hill coefficient. Experimental interrogation of common titratable systems based on an L-arabinose inducible promoter or a thiamine pyrophosphate-repressing riboswitch in Escherichia coli confirmed the predicted trends. We finally found that residual inducer had a distinct effect on "all-or-none" systems, which exhibited increased sensitivity to the added inducer until becoming fully induced. Our findings indicate that residual inducer or repressor alters the quantitative response properties of titratable systems, impacting their utility for scientific discovery and pathway engineering.

Mechanism of regulation of 'chromosome kissing' induced by Fob1 and its physiological significance

Protein-mediated "chromosome kissing" between two DNA sites in trans (or in cis) is known to facilitate three-dimensional control of gene expression and DNA replication. However, the mechanisms of regulation of the long-range interactions are unknown. Here, we show that the replication terminator protein Fob1 of Saccharomyces cerevisiae promoted chromosome kissing that initiated rDNA recombination and controlled the replicative life span (RLS). Oligomerization of Fob1 caused synaptic (kissing) interactions between pairs of terminator (Ter) sites that initiated recombination in rDNA. Fob1 oligomerization and Ter-Ter kissing were regulated by intramolecular inhibitory interactions between the C-terminal domain (C-Fob1) and the N-terminal domain (N-Fob1). Phosphomimetic substitutions of specific residues of C-Fob1 counteracted the inhibitory interaction. A mutation in either N-Fob1 that blocked Fob1 oligomerization or C-Fob1 that blocked its phosphorylation antagonized chromosome kissing and recombination and enhanced the RLS. The results provide novel insights into a mechanism of regulation of Fob1-mediated chromosome kissing.

Ellis Englesberg and the discovery of positive control in gene regulation

Based on his work with the Escherichia coli l-arabinose operon, Ellis Englesberg proposed in 1965 that the regulatory gene araC was an "activator gene" required for positive control of the ara operon. This challenged the widely held belief in a universal mechanism of negative regulation proposed earlier by Jacob and Monod. For years, Englesberg's model was met with deep skepticism. Despite much frustration with complex ad hoc explanations used to challenge his model, Englesberg persisted until the evidence for positive control in ara and other systems became overwhelming. Englesberg's pioneering work enriched the original operon model and had a lasting impact in opening new and exciting ways of thinking about transcriptional regulation.

Domain characterization of Bacillus subtilis GabR, a pyridoxal 5'-phosphate-dependent transcriptional regulator

Bacillus subtilis GabR is a transcriptional regulator consisting of a helix-turn-helix N-terminal DNA-binding domain, a pyridoxal 5'-phosphate (PLP)-binding C-terminal domain that has a structure homologous to aminotransferases, and a linker of 29 amino acid residues. In the presence of γ-aminobutyrate (GABA), GabR activates the transcription of gabT and gabD, which encode GABA aminotransferase and succinate semialdehyde dehydrogenase, respectively. We expressed N-terminal and C-terminal domain fragments (named N'-GabR and C'-GabR) in Escherichia coli cells, and obtained N'-GabR as a soluble monomer and C'-GabR as a soluble dimer. Spectroscopic studies suggested that C'-GabR contains PLP and binds to d-Ala, β-Ala, d-Asn and d-Gln, as well as GABA, although the intact GabR binds only to GABA. N'-GabR does not bind to the DNA fragment containing the GabR-binding sequence regardless of the presence or absence of C'-GabR. A fusion protein consisting of N'-GabR and 2-aminoadipate aminotransferase of Thermus thermophilus bound to the DNA fragment. These results suggested that each domain of GabR could be an independent folding unit. The C-terminal domain provides the N-terminal domain with DNA-binding ability via dimerization. The N-terminal domain controls the ligand specificity of the C-terminal domain. Connection by the linker is indispensable for the mutual interaction of the domains.

Redox-dependent DNA distortion in a SoxR protein-promoter complex studied using fluorescent probes

The [2Fe-2S] transcriptional factor SoxR, a member of the MerR family, is regulated by the reversible oxidation and reduction of [2Fe-2S] clusters and functions as a sensor of oxidative stress in Escherichia coli. In the oxidized state, distortion of the target DNA promoter region initiates transcription by RNA polymerase, thereby activating transcription. The inactive reduced state of the protein has remained uncharacterized. Here, we directly observed redox-dependent conformational changes in the promoter DNA by site-specifically replacing selected adenine (A) and cytosine (C) bases in the promoter oligonucleotide with the fluorescent probes 2-aminopurine (2Ap) and pyrrolocytosine (pyrrolo-dC), respectively. Reduction of the [2Fe-2S] cluster in the SoxR-DNA complex dramatically weakened the fluorescence intensity of the 2Ap moieties incorporated into the central part of the DNA. In contrast, the fluorescence of 2Ap moieties incorporated at A in other regions and the fluorescence of pyrrolo-dC moieties in the central region of the DNA (C3 and C3') were only slightly decreased by the reduction. These results strongly suggest that the redox change causes a large conformational change within a region confined to the central A-T base pairs in the promoter region of the DNA.

Virulence regulons of enterotoxigenic Escherichia coli

Enterotoxigenic Escherichia coli is frequently associated with travelers' diarrhea and is a leading cause of infant and childhood mortality in developing countries. Disease is dependent upon the orchestrated expression of enterotoxins, flexible adhesive pili, and other virulence factors. Both the heat-labile (LT) and heat-stable (ST-H) enterotoxins are regulated at the level of transcription by cAMP-receptor protein which represses the expression of LT while activating expression of ST-H. The expression of many different serotypes of adhesive pili is regulated by Rns, a member of the AraC family that represents a subgroup of conserved virulence regulators from several enteric pathogens. These Rns-like regulators recognize similar DNA binding sites, and a compiled sequence logo suggests they may bind DNA through both major and minor groove interactions. These regulators are also tempting targets for novel therapeutics because they play pivotal roles during infection. To that end, high-throughput screens have begun to identify compounds that inhibit the activity of these regulators, predominately by interfering with DNA binding.

Simple whole-cell biodetection and bioremediation of heavy metals based on an engineered lead-specific operon

A lead-specific binding protein, PbrR, and promoter pbr from the lead resistance operon, pbr, of Cupriavidus metallidurans CH34 was incorporated into E. coli in conjunction with an engineered downstream RFP (red fluorescence protein), which allowed for highly sensitive and selective whole-cell detection of lead ions. The subsequent display of PbrR on the E. coli cell surface permitted selective adsorption of lead ions from solution containing various heavy metal ions. The surface-engineered E. coli bacteria effectively protected Arabidopsis thaliana seed germination from the toxicity of lead ions at high concentrations. Engineering the E. coli bacteria harboring these lead-specific elements from the pbr operon may potentially be a valuable general strategy for biodetection and bioremediation of toxic heavy metal ions in the environment.

High-level expression of functionally active dengue-2 non-structural antigen 1 production in Escherichia coli

Detection of nonstructural protein (NS1) is an important diagnostic marker during acute phase of dengue infection. Not only for diagnostic purpose, the protein had important role in vaccine design as well, as a candidate for studying virus assembly and maturation. Various researchers employed different expression systems and strategies for recombinant NS1 protein production. Attempts to express NS1 protein in prokaryotic and yeast expression system result in formation of insoluble protein which needs to undergo refolding to attain native structural and functional forms. Here, we report the production of soluble NS1 protein in E. coli by using appropriate vector and employing suitable culture conditions to maximize protein production. Proteins were purified using metal affinity chromatography. SDS-PAGE and western blot analysis reveal the native structure of NS1 protein. Solid phase ELISA using the recombinantly expressed antigen with positive and negative dengue samples showed that the expressed protein retains its antigenic and immunological properties. To our knowledge, this is the first report on the successful production of functionally active recombinant dengue-2 NS1 protein production without undergoing any in vitro posttranslational modification process.

Understanding the basis of a class of paradoxical mutations in AraC through simulations

Most mutations at position 15 in the N-terminal arm of the regulatory protein AraC leave the protein incapable of responding to arabinose and inducing the proteins required for arabinose catabolism. Mutations at other positions of the arm do not have this behavior. Simple energetic analysis of the interactions between the arm and bound arabinose do not explain the uninducibility of AraC with mutations at position 15. Extensive molecular dynamics (MD) simulations, carried out largely on the Open Science Grid, were done of the wild-type protein with and without bound arabinose and of all possible mutations at position 15, many of which were constructed and measured for this work. Good correlation was found for deviation of arm position during the simulations and inducibility as measured in vivo of the same mutant proteins. Analysis of the MD trajectories revealed that preservation of the shape of the arm is critical to inducibility. To maintain the correct shape of the arm, the strengths of three interactions observed to be strong in simulations of the wild-type AraC protein need to be preserved. These interactions are between arabinose and residue 15, arabinose and residues 8-9, and residue 13 and residue 15. The latter interaction is notable because residues L9, Y13, F15, W95, and Y97 form a hydrophobic cluster which needs to be preserved for retention of the correct shape.

LacI-DNA-IPTG loops: equilibria among conformations by single-molecule FRET

The E. coli Lac repressor (LacI) tetramer binds simultaneously to a promoter-proximal DNA binding site (operator) and an auxiliary operator, resulting in a DNA loop, which increases repression efficiency. Induction of the lac operon by allolactose reduces the affinity of LacI for DNA, but induction does not completely prevent looping in vivo. Our previous work on the conformations of LacI loops used a hyperstable model DNA construct, 9C14, that contains a sequence directed bend flanked by operators. Single-molecule fluorescence resonance energy transfer (SM-FRET) on a dual fluorophore-labeled LacI-9C14 loop showed that it adopts a single, stable, high-FRET V-shaped LacI conformation. Ligand-induced changes in loop geometry can affect loop stability, and the current work assesses loop population distributions for LacI-9C14 complexes containing the synthetic inducer IPTG. SM-FRET confirms that the high-FRET LacI-9C14 loop is only partially destabilized by saturating IPTG. LacI titration experiments and FRET fluctuation analysis suggest that the addition of IPTG induces loop conformational dynamics and re-equilibration between loop population distributions that include a mixture of looped states that do not exhibit high-efficiency FRET. The results show that repression by looping even at saturating IPTG should be considered in models for regulation of the operon. We propose that persistent DNA loops near the operator function biologically to accelerate rerepression upon exhaustion of inducer.

Structures of the Escherichia coli transcription activator and regulator of diauxie, XylR: an AraC DNA-binding family member with a LacI/GalR ligand-binding domain

Escherichia coli can rapidly switch to the metabolism of l-arabinose and d-xylose in the absence of its preferred carbon source, glucose, in a process called carbon catabolite repression. Transcription of the genes required for l-arabinose and d-xylose consumption is regulated by the sugar-responsive transcription factors, AraC and XylR. E. coli represents a promising candidate for biofuel production through the metabolism of hemicellulose, which is composed of d-xylose and l-arabinose. Understanding the l-arabinose/d-xylose regulatory network is key for such biocatalyst development. Unlike AraC, which is a well-studied protein, little is known about XylR. To gain insight into XylR function, we performed biochemical and structural studies. XylR contains a C-terminal AraC-like domain. However, its N-terminal d-xylose-binding domain contains a periplasmic-binding protein (PBP) fold with structural homology to LacI/GalR transcription regulators. Like LacI/GalR proteins, the XylR PBP domain mediates dimerization. However, unlike LacI/GalR proteins, which dimerize in a parallel, side-to-side manner, XylR PBP dimers are antiparallel. Strikingly, d-xylose binding to this domain results in a helix to strand transition at the dimer interface that reorients both DNA-binding domains, allowing them to bind and loop distant operator sites. Thus, the combined data reveal the ligand-induced activation mechanism of a new family of DNA-binding proteins.

Multiscale modeling of double-helical DNA and RNA: a unification through Lie groups

Several different mechanical models of double-helical nucleic-acid structures that have been presented in the literature are reviewed here together with a new analysis method that provides a reconciliation between these disparate models. In all cases, terminology and basic results from the theory of Lie groups are used to describe rigid-body motions in a coordinate-free way, and when necessary, coordinates are introduced in a way in which simple equations result. We consider double-helical DNAs and RNAs which, in their unstressed referential state, have backbones that are either straight, slightly precurved, or bent by the action of a protein or other bound molecule. At the coarsest level, we consider worm-like chains with anisotropic bending stiffness. Then, we show how bi-rod models converge to this for sufficiently long filament lengths. At a finer level, we examine elastic networks of rigid bases and show how these relate to the coarser models. Finally, we show how results from molecular dynamics simulation at full atomic resolution (which is the finest scale considered here) and AFM experimental measurements (which is at the coarsest scale) relate to these models.

Development and application of an arabinose-inducible expression system by facilitating inducer uptake in Corynebacterium glutamicum

Corynebacterium glutamicum is currently used for the industrial production of a variety of biological materials. Many available inducible expression systems in this species use lac-derived promoters from Escherichia coli that exhibit much lower levels of inducible expression and leaky basal expression. We developed an arabinose-inducible expression system that contains the L-arabinose regulator AraC, the P(BAD) promoter from the araBAD operon, and the L-arabinose transporter AraE, all of which are derived from E. coli. The level of inducible P(BAD)-based expression could be modulated over a wide concentration range from 0.001 to 0.4% L-arabinose. This system tightly controlled the expression of the uracil phosphoribosyltransferase without leaky expression. When the gene encoding green fluorescent protein (GFP) was under the control of P(BAD) promoter, flow cytometry analysis showed that GFP was expressed in a highly homogeneous profile throughout the cell population. In contrast to the case in E. coli, P(BAD) induction was not significantly affected in the presence of different carbon sources in C. glutamicum, which makes it useful in fermentation applications. We used this system to regulate the expression of the odhI gene from C. glutamicum, which encodes an inhibitor of α-oxoglutarate dehydrogenase, resulting in high levels of glutamate production (up to 13.7 mM) under biotin nonlimiting conditions. This system provides an efficient tool available for molecular biology and metabolic engineering of C. glutamicum.

An E. coli cell-free expression toolbox: application to synthetic gene circuits and artificial cells

Cell-free protein synthesis is becoming a powerful technique to construct and to study complex informational processes in vitro. Engineering synthetic gene circuits in a test tube, however, is seriously limited by the transcription repertoire of modern cell-free systems, composed of only a few bacteriophage regulatory elements. Here, we report the construction and the phenomenological characterization of synthetic gene circuits engineered with a cell-free expression toolbox that works with the seven E. coli sigma factors. The E. coli endogenous holoenzyme E(70) is used as the primary transcription machinery. Elementary circuit motifs, such as multiple stage cascades, AND gate and negative feedback loops are constructed with the six other sigma factors, two bacteriophage RNA polymerases, and a set of repressors. The circuit dynamics reveal the importance of the global mRNA turnover rate and of passive competition-induced transcriptional regulation. Cell-free reactions can be carried out over long periods of time with a small-scale dialysis reactor or in phospholipid vesicles, an artificial cell system. This toolbox is a unique platform to study complex transcription/translation-based biochemical systems in vitro.

Modeling synthetic gene oscillators

Genetic oscillators have long held the fascination of experimental and theoretical synthetic biologists alike. From an experimental standpoint, the creation of synthetic gene oscillators represents a yardstick by which our ability to engineer synthetic gene circuits can be measured. For theorists, synthetic gene oscillators are a playground in which to test mathematical models for the dynamics of gene regulation. Historically, mathematical models of synthetic gene circuits have varied greatly. Often, the differences are determined by the level of biological detail included within each model, or which approximation scheme is used. In this review, we examine, in detail, how mathematical models of synthetic gene oscillators are derived and the biological processes that affect the dynamics of gene regulation.

ChIP-Seq and the complexity of bacterial transcriptional regulation

Transcription factors (TFs) play a central role in regulating gene expression in all bacteria. Yet, until recently, studies of TF binding were limited to a small number of factors at a few genomic locations. Chromatin immunoprecipitation followed by sequencing enables mapping of binding sites for TFs in a global and high-throughput fashion. The NIAID funded TB systems biology project http://www.broadinstitute.org/annotation/tbsysbio/home.html aims to map the binding sites for every transcription factor in the genome of Mycobacterium tuberculosis (MTB), the causative agent of human TB. ChIP-Seq data already released through TBDB.org have provided new insight into the mechanisms of TB pathogenesis. But in addition, data from MTB are beginning to challenge many simplifying assumptions associated with gene regulation in all bacteria. In this chapter, we review the global aspects of TF binding in MTB and discuss the implications of these data for our understanding of bacterial gene regulation. We begin by reviewing the canonical model of bacterial transcriptional regulation using the lac operon as the standard paradigm. We then review the use of ChIP-Seq to map the binding sites of DNA-binding proteins and the application of this method to mapping TF binding sites in MTB. Finally, we discuss two aspects of the binding discovered by ChIP-Seq that were unexpected given the canonical model: the substantial binding outside the proximal promoter region and the large number of weak binding sites.

Resistance to topoisomerase cleavage complex induced lethality in Escherichia coli via titration of transcription regulators PurR and FNR

BACKGROUND

Accumulation of gyrase cleavage complex in Escherichia coli from the action of quinolone antibiotics induces an oxidative damage cell death pathway. The oxidative cell death pathway has also been shown to be involved in the lethality following accumulation of cleavage complex formed by bacterial topoisomerase I with mutations that result in defective DNA religation.

METHODS

A high copy number plasmid clone spanning the upp-purMN region was isolated from screening of an E. coli genomic library and analyzed for conferring increased survival rates following accumulation of mutant topoisomerase I proteins as well as treatment with the gyrase inhibitor norfloxacin.

RESULTS

Analysis of the intergenic region upstream of purM demonstrated a novel mechanism of resistance to the covalent protein-DNA cleavage complex through titration of the cellular transcription regulators FNR and PurR responsible for oxygen sensing and repression of purine nucleotide synthesis respectively. Addition of adenine to defined growth medium had similar protective effect for survival following accumulation of topoisomerase cleavage complex, suggesting that increase in purine level can protect against cell death.

CONCLUSIONS

Perturbation of the global regulator FNR and PurR functions as well as increase in purine nucleotide availability could affect the oxidative damage cell death pathway initiated by topoisomerase cleavage complex.

Active role of the interdomain linker of AraC

Mutations in the interdomain linker of the gene for the AraC regulatory protein of Escherichia coli that severely interfere with the protein's ability either to repress or to activate transcription have been found. These mutations have relatively small effects on the dimerization domain's ability to bind arabinose or to dimerize the protein or on the DNA-binding domain's affinity for a single DNA half-site. The linker mutations, however, dramatically change the affinity of AraC for binding to two direct-repeat DNA half-sites. Less dramatically, the induction-deficient linker variants also display altered DNA sequence selectivity. These results show that changing the sequence of the interdomain linker can profoundly affect the dimerization domain-DNA-binding domain interactions in AraC. The smaller effects on the functions of the individual domains could be the direct result of the linker alterations but more likely are the indirect result of the altered dimerization domain-DNA-binding domain interactions. In summary, the linker does not simply function as a passive and flexible connector between the domains of AraC but, instead, is more directly involved in the protein's dimerization domain-DNA-binding domain interactions.

Negative auto-regulation increases the input dynamic-range of the arabinose system of Escherichia coli

Background

Gene regulation networks are made of recurring regulatory patterns, called network motifs. One of the most common network motifs is negative auto-regulation, in which a transcription factor represses its own production. Negative auto-regulation has several potential functions: it can shorten the response time (time to reach halfway to steady-state), stabilize expression against noise, and linearize the gene's input-output response curve. This latter function of negative auto-regulation, which increases the range of input signals over which downstream genes respond, has been studied by theory and synthetic gene circuits. Here we ask whether negative auto-regulation preserves this function also in the context of a natural system, where it is embedded within many additional interactions. To address this, we studied the negative auto-regulation motif in the arabinose utilization system of Escherichia coli, in which negative auto-regulation is part of a complex regulatory network.

Results

We find that when negative auto-regulation is disrupted by placing the regulator araC under constitutive expression, the input dynamic range of the arabinose system is reduced by 10-fold. The apparent Hill coefficient of the induction curve changes from about n = 1 with negative auto-regulation, to about n = 2 when it is disrupted. We present a mathematical model that describes how negative auto-regulation can increase input dynamic-range, by coupling the transcription factor protein level to the input signal.

Conclusions

Here we demonstrate that the negative auto-regulation motif in the native arabinose system of Escherichia coli increases the range of arabinose signals over which the system can respond. In this way, negative auto-regulation may help to increase the input dynamic-range while maintaining the specificity of cooperative regulatory systems. This function may contribute to explaining the common occurrence of negative auto-regulation in biological systems.

Cellular heterogeneity and live cell arrays

In the past decade, the tendency to move from a global, one-size-fits-all treatment philosophy to personalized medicine is based, in part, on the nuanced differences and sub-classifications of disease states. Our knowledge of these varied states stems from not only the ability to diagnose, classify, and perform experiments on cell populations as a whole, but also from new technologies that allow interrogation of cell populations at the individual cell level. Such departures from conventional thinking are driven by the recognition that clonal cell populations have numerous activities that manifest as significant levels of non-genetic heterogeneity. Clonal populations by definition originate from a single genetic origin so are regarded as having a high level of homogeneity as compared to genetically distinct cell populations. However, analysis at the single cell level has revealed a different phenomenon; cells and organisms require an inherent level of non-genetic heterogeneity to function properly, and in some cases, to survive. The growing understanding of this occurrence has lead to the development of methods to monitor, analyze, and better characterize the heterogeneity in cell populations. Following the trend of DNA- and protein microarrays, platforms capable of simultaneously monitoring each cell in a population have been developed. These cellular microarray platforms and other related formats allow for continuous monitoring of single live cells and simultaneously generate individual cell and average population data that are more descriptive and information-rich than traditional bulk methods. These technological advances have helped develop a better understanding of the intricacies associated with biological processes and afforded greater insight into complex biological systems. The associated instruments, techniques, and reagents now allow for highly multiplexed analyses, which enable multiple cellular activities, processes, or pathways to be monitored simultaneously. This critical review will discuss the paradigm shift associated with cellular heterogeneity, speak to the key developments that have lead to our better understanding of systems biology, and detail the future directions of the discipline (281 references).

Understanding apparent DNA flexibility enhancement by HU and HMGB architectural proteins

Understanding and predicting the mechanical properties of protein/DNA complexes are challenging problems in biophysics. Certain architectural proteins bind DNA without sequence specificity and strongly distort the double helix. These proteins rapidly bind and unbind, seemingly enhancing the flexibility of DNA as measured by cyclization kinetics. The ability of architectural proteins to overcome DNA stiffness has important biological consequences, but the detailed mechanism of apparent DNA flexibility enhancement by these proteins has not been clear. Here, we apply a novel Monte Carlo approach that incorporates the precise effects of protein on DNA structure to interpret new experimental data for the bacterial histone-like HU protein and two eukaryotic high-mobility group class B (HMGB) proteins binding to ∼200-bp DNA molecules. These data (experimental measurement of protein-induced increase in DNA cyclization) are compared with simulated cyclization propensities to deduce the global structure and binding characteristics of the closed protein/DNA assemblies. The simulations account for all observed (chain length and concentration dependent) effects of protein on DNA behavior, including how the experimental cyclization maxima, observed at DNA lengths that are not an integral helical repeat, reflect the deformation of DNA by the architectural proteins and how random DNA binding by different proteins enhances DNA cyclization to different levels. This combination of experiment and simulation provides a powerful new approach to resolve a long-standing problem in the biophysics of protein/DNA interactions.