Kinetic and thermodynamic studies of purine repressor binding to corepressor and operator DNA

Synthetic cell-based materials extract positional information from morphogen gradients

Biomaterials composed of synthetic cells have the potential to adapt and differentiate guided by physicochemical environmental cues. Inspired by biological systems in development, which extract positional information (PI) from morphogen gradients in the presence of uncertainties, we here investigate how well synthetic cells can determine their position within a multicellular structure. To calculate PI, we created and analyzed a large number of synthetic cellular assemblies composed of emulsion droplets connected via lipid bilayer membranes. These droplets contained cell-free feedback gene circuits that responded to gradients of a genetic inducer acting as a morphogen. PI is found to be limited by gene expression noise and affected by the temporal evolution of the morphogen gradient and the cell-free expression system itself. The generation of PI can be rationalized by computational modeling of the system. We scale our approach using three-dimensional printing and demonstrate morphogen-based differentiation in larger tissue-like assemblies.

Biomolecular Assemblies: Moving from Observation to Predictive Design

Biomolecular assembly is a key driving force in nearly all life processes, providing structure, information storage, and communication within cells and at the whole organism level. These assembly processes rely on precise interactions between functional groups on nucleic acids, proteins, carbohydrates, and small molecules, and can be fine-tuned to span a range of time, length, and complexity scales. Recognizing the power of these motifs, researchers have sought to emulate and engineer biomolecular assemblies in the laboratory, with goals ranging from modulating cellular function to the creation of new polymeric materials. In most cases, engineering efforts are inspired or informed by understanding the structure and properties of naturally occurring assemblies, which has in turn fueled the development of predictive models that enable computational design of novel assemblies. This Review will focus on selected examples of protein assemblies, highlighting the story arc from initial discovery of an assembly, through initial engineering attempts, toward the ultimate goal of predictive design. The aim of this Review is to highlight areas where significant progress has been made, as well as to outline remaining challenges, as solving these challenges will be the key that unlocks the full power of biomolecules for advances in technology and medicine.

NF-κB and IRF pathways: cross-regulation on target genes promoter level

Background

The NF-κB and IRF transcription factor families are major players in inflammation and antiviral response and act as two major effectors of the innate immune response (IIR). The regulatory mechanisms of activation of these two pathways and their interactions during the IIR are only partially known.

Results

Our in silico findings report that there is cross-regulation between both pathways at the level of gene transcription regulation, mediated by the presence of binding sites for both factors in promoters of genes essential for these pathways. These findings agree with recent experimental data reporting crosstalk between pathways activated by RIG-I and TLR3 receptors in response to pathogens.

Conclusions

We present an extended crosstalk diagram of the IRF - NF-κB pathways. We conclude that members of the NF-κB family may directly impact regulation of IRF family, while IRF members impact regulation of NF-κB family rather indirectly, via other transcription factors such as AP-1 and SP1.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1511-7) contains supplementary material, which is available to authorized users.

Monte Carlo simulations of single and coupled synthetic molecular motors

We use a minimal model to study the processive motion of coupled synthetic molecular motors along a DNA track and we present data from Monte Carlo (MC) computer simulations based on this model. The model was originally proposed by Bromley et al. [HFSP J. 3, 204 (2009)] for studying the properties of a synthetic protein motor, the "Tumbleweed" (TW), and involves rigid Y-shaped motors diffusively rotating along the track while controlled by a series of periodically injected ligand pulses into the solution. The advantage of the model is that it mimics the mechanical properties of the TW motor in detail. Both the average first passage time which measures the diffusive motion of the motors, and the average dwell time on the track which measures their processivity are investigated by varying the parameters of the model. The latter includes ligand concentration and the range and strength of the binding interaction between motors and the track. In particular, it is of experimental interest to study the dependence of these dynamic time scales of the motors on the ligand concentration. Single rigid TW motors were first studied since no previous MC simulations of these motors have been performed. We first studied single motors for which we found a logarithmic decrease of the average first passage time and a logarithmic increase of the average dwell time with increasing ligand concentration. For two coupled motors, the dependence on ligand concentration is still logarithmic for the average first passage time but becomes linear for the average dwell time. This suggests a much greater stability in the processive motion of coupled motors as compared to single motors in the limit of large ligand concentration. By increasing the number of coupled motors, m, it was found that the average first passage time of the coupled motors only increases slowly with m while the average dwell time increases exponentially with m. Thus the stability of coupled motors on the track can be considerably enhanced by their cooperative motion.

The magnitude and colour of noise in genetic negative feedback systems

The comparative ability of transcriptional and small RNA-mediated negative feedback to control fluctuations or ‘noise’ in gene expression remains unexplored. Both autoregulatory mechanisms usually suppress the average (mean) of the protein level and its variability across cells. The variance of the number of proteins per molecule of mean expression is also typically reduced compared with the unregulated system, but is almost never below the value of one. This relative variance often substantially exceeds a recently obtained, theoretical lower limit for biochemical feedback systems. Adding the transcriptional or small RNA-mediated control has different effects. Transcriptional autorepression robustly reduces both the relative variance and persistence (lifetime) of fluctuations. Both benefits combine to reduce noise in downstream gene expression. Autorepression via small RNA can achieve more extreme noise reduction and typically has less effect on the mean expression level. However, it is often more costly to implement and is more sensitive to rate parameters. Theoretical lower limits on the relative variance are known to decrease slowly as a measure of the cost per molecule of mean expression increases. However, the proportional increase in cost to achieve substantial noise suppression can be different away from the optimal frontier—for transcriptional autorepression, it is frequently negligible.

Classes of fast and specific search mechanisms for proteins on DNA

Problems of search and recognition appear over different scales in biological systems. In this review we focus on the challenges posed by interactions between proteins, in particular transcription factors, and DNA and possible mechanisms which allow for fast and selective target location. Initially we argue that DNA-binding proteins can be classified, broadly, into three distinct classes which we illustrate using experimental data. Each class calls for a different search process and we discuss the possible application of different search mechanisms proposed over the years to each class. The main thrust of this review is a new mechanism which is based on barrier discrimination. We introduce the model and analyze in detail its consequences. It is shown that this mechanism applies to all classes of transcription factors and can lead to a fast and specific search. Moreover, it is shown that the mechanism has interesting transient features which allow for stability at the target despite rapid binding and unbinding of the transcription factor from the target.

Bacillus subtilis CodY operators contain overlapping CodY binding sites

CodY is a global transcriptional regulator that is activated by branched-chain amino acids. A palindromic 15-bp sequence motif, AATTTTCNGAAAATT, is associated with CodY DNA binding. A gel mobility shift assay was used to examine the effect of pH on the binding of Bacillus subtilis CodY to the hutPp and ureAp(3) promoters. CodY at pH 6.0 has higher affinity for DNA, more enhanced activation by isoleucine, and a lower propensity for nonspecific DNA binding than CodY at pH 8.0. DNase I footprinting was used to identify the CodY-protected regions in the hutPp and ureAp(3) promoters. The CodY-protected sequences for both promoters were found to contain multiple copies of the 15-bp motif with 6-bp overlaps. Mutational analysis of the hutPp regulatory region revealed that two overlapping sequence motifs were required for CodY-mediated regulation. The presence of overlapping sequence motifs in the regulatory regions of many B. subtilis CodY-regulated genes suggests that CodY binds to native operators that contain overlapping binding sites.

Flexibility in the inducer binding region is crucial for allostery in the Escherichia coli lactose repressor

Lactose repressor protein (LacI) utilizes an allosteric mechanism to regulate transcription in Escherichia coli, and the transition between inducer- and operator-bound states has been simulated by targeted molecular dynamics (TMD). The side chains of amino acids 149 and 193 interact and were predicted by TMD simulation to play a critical role in the early stages of the LacI conformational change. D149 contacts IPTG directly, and variations at this site provide the opportunity to dissect its role in inducer binding and signal transduction. Single mutants at D149 or S193 exhibit a minimal change in operator binding, and alterations in inducer binding parallel changes in operator release, indicating normal allosteric response. The observation that the double mutant D149A/S193A exhibits wild-type properties excludes the requirement for inter-residue hydrogen bond formation in the allosteric response. The double mutant D149C/S193C purified from cell extracts shows decreased sensitivity to inducer binding while retaining wild-type binding affinities and kinetic constants for both operator and inducer. By manipulating cysteine oxidation, we show that the more reduced state of D149C/S193C responds to inducer more like the wild-type protein, whereas the more oxidized state displays diminished inducer sensitivity. These features of D149C/S193C indicate that the novel disulfide bond formed in this mutant impedes the allosteric transition, consistent with the role of this region predicted by TMD simulation. Together, these results establish the requirement for flexibility in the spatial relationship between D149 and S193 rather than a specific D149-S193 interaction in the LacI allosteric response to inducer.

Functional consequences of exchanging domains between LacI and PurR are mediated by the intervening linker sequence

Homologue function can be differentiated by changing residues that affect binding sites or long-range interactions. LacI and PurR are two proteins that represent the LacI/GalR family (>500 members) of bacterial transcription regulators. All members have distinct DNA-binding and regulatory domains linked by approximately 18 amino acids. Each homologue has specificity for different DNA and regulatory effector ligands; LacI and PurR also exhibit differences in allosteric communication between DNA and effector binding sites. A comparative study of LacI and PurR suggested that alterations in the interface between the regulatory domain and linker are important for differentiating their functions. Four residues (equivalent to LacI positions 48, 55, 58, and 61) appear particularly important for creating a unique interface and were predicted to be necessary for allosteric regulation. However, nearby residues in the linker interact with DNA ligand. Thus, differences observed in interactions between linker and regulatory domain may be the cause of altered function or an effect of the two proteins binding different DNA ligands. To separate these possibilities, we created a chimeric protein with the LacI DNA-binding domain/linker and the PurR regulatory domain (LLhP). If the interface requires homologue-specific interactions in order to propagate the signal from effector binding, then LLhP repression should not be allosterically regulated by effector binding. Experiments show that LLhP is capable of repression from lacO1 and, contrary to expectation, allosteric response is intact. Further, restoring the potential for PurR-like interactions via substitutions in the LLhP linker tends to diminish repression. These effects are especially pronounced for residues 58 and 61. Clearly, binding affinity of LLhP for the lacO1 DNA site is sensitive to long-range changes in the linker. This result also raises the possibility that mutations at positions 58 and 61 co-evolved with changes in the DNA-binding site. In addition, repression measured in the absence and presence of effector ligand shows that allosteric response increases for several LLhP variants with substitutions at positions 48 and 55. Thus, while side chain variation at these sites does not generally dictate the presence or absence of allostery, the nature of the amino acid can modulate the response to effector.

Innate immune responses: crosstalk of signaling and regulation of gene transcription

Innate immune responses to pathogens such as bacteria and viruses are triggered by recognition of specific structures of invading pathogens called pathogen-associated molecular patterns (PAMPs) by cellular pattern recognition receptors (PRRs) that are located at plasma membrane or inside cells. Stimulation of different PAMPs activates Toll-like receptor (TLR)-dependent and -independent signaling pathways that lead to activation of transcription factors nuclear factor-kappaB (NF-kappaB), interferon regulatory factor 3/7 (IRF3/7) and/or activator protein-1 (AP-1), which collaborate to induce transcription of a large number of downstream genes. This review focuses on the rapid progress that has recently improved our understanding of the crosstalk among the pathways and the precise regulation of transcription of the downstream genes.

Evidence that Bacillus catabolite control protein CcpA interacts with RNA polymerase to inhibit transcription

Summary Bacilluscatabolite control protein (CcpA) mediates carbon catabolite repression (CCR) by controlling expression of catabolite responsive (CR) genes or operons through interaction with catabolite responsive elements (cres) located within or outside of CR promoters. Here, we investigated how CcpA inhibits the transcription of CR promoters in vitro. CcpA has different affinities for different cres, but this does not correlate with its ability to inhibit transcription. In the amyE promoter, which overlaps a CcpA binding site (amyE cre centred at +4.5), CcpA does not prevent RNA polymerase (RNAP) binding to the promoter; it may even interact with RNAP. Inserting non-integral turns of helix (1.5 and 2.5) between the amyE promoter (-10 hexamer) and the amyE cre relieved CCR of amyE expression. In the xyl operon, despite the downstream location of its cre (a major cre centred at +130.5), CcpA blocked transcription initiation, not elongation (roadblock) at the site of the cre. Taken together, our results strongly suggest that CcpA requires interactions with RNAP to inhibit transcription.

Quantitative interdependence of coeffectors, CcpA and cre in carbon catabolite regulation of Bacillus subtilis

The phosphoproteins HPrSerP and CrhP are the main effectors for CcpA-mediated carbon catabolite regulation (CCR) in Bacillus subtilis. Complexes of CcpA with HPrSerP or CrhP regulate genes by binding to the catabolite responsive elements (cre). We present a quantitative analysis of HPrSerP and CrhP interaction with CcpA by surface plasmon resonance (SPR) revealing small and similar equilibrium constants of 4.8 +/- 0.4 microm for HPrSerP-CcpA and 19.1 +/- 2.5 microm for CrhP-CcpA complex dissociation. Forty millimolar fructose-1,6-bisphosphate (FBP) or glucose-6-phosphate (Glc6-P) increases the affinity of HPrSerP to CcpA at least twofold, but have no effect on CrhP-CcpA binding. Saturation of binding of CcpA to cre as studied by fluorescence and SPR is dependent on 50 microm of HPrSerP or > 200 microm CrhP. The rate constants of HPrSerP-CcpA-cre complex formation are k(a) = 3 +/- 1 x 10(6) m(-1).s(-1) and k(d) = 2.0 +/- 0.4 x 10(-3).s(-1), resulting in a K(D) of 0.6 +/- 0.3 nm. FBP and Glc6-P stimulate CcpA-HPrSerP but not CcpA-CrhP binding to cre. Maximal HPrSerP-CcpA-cre complex formation in the presence of 10 mm FBP requires about 10-fold less HPrSerP. These data suggest a specific role for FBP and Glc6-P in enhancing only HPrSerP-mediated CCR.

Design principles for regulator gene expression in a repressible gene circuit

We consider the design of a type of repressible gene circuit that is common in bacteria. In this type of circuit, a regulator protein acts to coordinately repress the expression of effector genes when a signal molecule with which it interacts is present. The regulator protein can also independently influence the expression of its own gene, such that regulator gene expression is repressible (like effector genes), constitutive, or inducible. Thus, a signal-directed change in the activity of the regulator protein can result in one of three patterns of coupled regulator and effector gene expression: direct coupling, in which regulator and effector gene expression change in the same direction; uncoupling, in which regulator gene expression remains constant while effector gene expression changes; or inverse coupling, in which regulator and effector gene expression change in opposite directions. We have investigated the functional consequences of each form of coupling using a mathematical model to compare alternative circuits on the basis of engineering-inspired criteria for functional effectiveness. The results depend on whether the regulator protein acts as a repressor or activator of transcription at the promoters of effector genes. In the case of repressor control of effector gene expression, direct coupling is optimal among the three forms of coupling, whereas in the case of activator control, inverse coupling is optimal. Results also depend on the sensitivity of effector gene expression to changes in the level of a signal molecule; the optimal form of coupling can be physically realized only for circuits with sufficiently small sensitivity. These theoretical results provide a rationale for autoregulation of regulator genes in repressible gene circuits and lead to testable predictions, which we have compared with data available in the literature and electronic databases.

Quantification of the influence of HPrSer46P on CcpA-cre interaction

Carbon catabolite repression (CCR) of the Bacillus megateriumxyl operon is dependent on the catabolite responsive element cre, the catabolite control protein (CcpA) and the histidine-containing phosphocarrier protein phosphorylated at the serine 46 residue (HPrSer46P). The latter is formed in the presence of glucose and mediates CCR via CcpA. We present evidence for the presence of HPrSer46P in a ternary complex with CcpA and cre. We also demonstrate increased stability of this complex compared to the CcpA-cre complex by electrophoretic mobility shift analysis (EMSA). This stabilization by HPrSer46P is the same for the xyl cre and an improved cre. Thus, HPrSer46P is a co-repressor for CcpA. In addition, surface plasmon resonance (SPR) experiments yielded binding constants of CcpA and the CcpA-HPrSer46P complex with cre. HPrSer46P stimulated CcpA binding to cre 50-fold. The binding constant is 4.9(+/- 0.5) x 10(6) M(-1). Non-phosphorylated HPr did not affect the complex formation between CcpA and cre. Previously proposed effects by glucose-6-phosphate, fructose-1,6-diphosphate and NADP on CcpA-cre or CcpA-HPrSer46P-cre formation were not found in EMSA and SPR experiments.

Fine-tuning function: correlation of hinge domain interactions with functional distinctions between LacI and PurR

LacI and PurR are highly homologous proteins. Their functional units are homodimers, with an N-terminal DNA binding domain that comprises the helix-turn-helix (HTH), N-linker, and hinge regions from both monomers. Hinge structural changes are known to occur upon DNA dissociation but are difficult to monitor experimentally. The initial steps of hinge unfolding were therefore examined using molecular dynamics simulations, utilizing a truncated, chimeric protein comprising the LacI HTH/N-linker and PurR hinge. A terminal Gly-Cys-Gly was added to allow "dimerization" through disulfide bond formation. Simulations indicate that differences in LacI and PurR hinge primary sequence affect the quaternary structure of the hinge x hinge' interface. However, these alternate hinge orientations would be sterically restricted by the core domain. These results prompted detailed comparison of recently available DNA-bound structures for LacI and truncated LacI(1-62) with the PurR structure. Examination revealed that different N-linker and hinge contacts to the core domain of the partner monomer (which binds effector molecule) affect the juxtapositions of the HTH, N-linker, and hinge regions in the DNA binding domain. In addition, the two full-length repressors exhibit significant differences in the interactions between the core and the C-linker connection to the DNA binding domain. Both linkers and the hinge have been implicated in the allosteric response of these repressors. Intriguingly, one functional difference between these two proteins is that they exhibit opposite allosteric response to effector. Simulations and observed structural distinctions are correlated with mutational analysis and sequence information from the LacI/GalR family to formulate a mechanism for fine-tuning individual repressor function.

Real-time kinetic analysis of hCG-monoclonal antibody interaction using radiolabeled hCG probe: presence of two forms of Ag-mAb complex as revealed by solid phase dissociation studies

Real-time kinetics of ligand-ligate interaction has predominantly been studied by either fluorescence or surface plasmon resonance based methods. Almost all such studies are based on association between the ligand and the ligate. This paper reports our analysis of dissociation data of monoclonal antibody-antigen (hCG) system using radio-iodinated hCG as a probe and nitrocellulose as a solid support to immobilize mAb. The data was analyzed quantitatively for a one-step and a two-step model. The data fits well into the two-step model. We also found that a fraction of what is bound is non-dissociable (tight-binding portion (TBP)). The TBP was neither an artifact of immobilization nor does it interfere with analysis. It was present when the reaction was carried out in homogeneous solution in liquid phase. The rate constants obtained from the two methods were comparable. The work reported here shows that real-time kinetics of other ligand-ligate interaction can be studied using nitrocellulose as a solid support.