Lambda repressor: a model system for understanding protein-DNA interactions and protein stability

Biophysical ambiguities prevent accurate genetic prediction

A goal of biology is to predict how mutations combine to alter phenotypes, fitness and disease. It is often assumed that mutations combine additively or with interactions that can be predicted. Here, we show using simulations that, even for the simple example of the lambda phage transcription factor CI repressing a gene, this assumption is incorrect and that perfect measurements of the effects of mutations on a trait and mechanistic understanding can be insufficient to predict what happens when two mutations are combined. This apparent paradox arises because mutations can have different biophysical effects to cause the same change in a phenotype and the outcome in a double mutant depends upon what these hidden biophysical changes actually are. Pleiotropy and non-monotonic functions further confound prediction of how mutations interact. Accurate prediction of phenotypes and disease will sometimes not be possible unless these biophysical ambiguities can be resolved using additional measurements.

RecA: Regulation and Mechanism of a Molecular Search Engine

Homologous recombination maintains genomic integrity by repairing broken chromosomes. The broken chromosome is partially resected to produce single-stranded DNA (ssDNA) that is used to search for homologous double-stranded DNA (dsDNA). This homology driven 'search and rescue' is catalyzed by a class of DNA strand exchange proteins that are defined in relation to Escherichia coli RecA, which forms a filament on ssDNA. Here, we review the regulation of RecA filament assembly and the mechanism by which RecA quickly and efficiently searches for and identifies a unique homologous sequence among a vast excess of heterologous DNA. Given that RecA is the prototypic DNA strand exchange protein, its behavior affords insight into the actions of eukaryotic RAD51 orthologs and their regulators, BRCA2 and other tumor suppressors.

A self-lysis pathway that enhances the virulence of a pathogenic bacterium

In mammalian cells, programmed cell death (PCD) plays important roles in development, in the removal of damaged cells, and in fighting bacterial infections. Although widespread among multicellular organisms, there are relatively few documented instances of PCD in bacteria. Here we describe a potential PCD pathway in Pseudomonas aeruginosa that enhances the ability of the bacterium to cause disease in a lung infection model. Activation of the system can occur in a subset of cells in response to DNA damage through cleavage of an essential transcription regulator we call AlpR. Cleavage of AlpR triggers a cell lysis program through de-repression of the alpA gene, which encodes a positive regulator that activates expression of the alpBCDE lysis cassette. Although this is lethal to the individual cell in which it occurs, we find it benefits the population as a whole during infection of a mammalian host. Thus, host and pathogen each may use PCD as a survival-promoting strategy. We suggest that activation of the Alp cell lysis pathway is a disease-enhancing response to bacterial DNA damage inflicted by the host immune system.

Fast protein folding kinetics

Fast-folding proteins have been a major focus of computational and experimental study because they are accessible to both techniques: they are small and fast enough to be reasonably simulated with current computational power, but have dynamics slow enough to be observed with specially developed experimental techniques. This coupled study of fast-folding proteins has provided insight into the mechanisms, which allow some proteins to find their native conformation well <1 ms and has uncovered examples of theoretically predicted phenomena such as downhill folding. The study of fast folders also informs our understanding of even 'slow' folding processes: fast folders are small; relatively simple protein domains and the principles that govern their folding also govern the folding of more complex systems. This review summarizes the major theoretical and experimental techniques used to study fast-folding proteins and provides an overview of the major findings of fast-folding research. Finally, we examine the themes that have emerged from studying fast folders and briefly summarize their application to protein folding in general, as well as some work that is left to do.

Insights into the diversity of φRSM phages infecting strains of the phytopathogen Ralstonia solanacearum complex: regulation and evolution

The filamentous φRSM phages (φRSM1 and φRSM3) have integration/excision capabilities in the phytopathogenic bacterium Ralstonia solanacearum. In the present study, we further investigated φRSM-like sequences present in the genomes of R. solanacearum strains belonging to the four major evolutionary lineages (phylotypes I-IV). Based on bioinformatics and comparative genomic analyses, we found that φRSM homologs are highly diverse in R. solanacearum complex strains. We detected an open reading frame (ORF)15 located upstream of the gene for φRSM integrase, which exhibited amino acid sequence similarity to phage repressor proteins. ORF15-encoded protein (a putative repressor) was found to encode a 104-residue polypeptide containing a DNA-binding (helix-turn-helix) domain and was expressed in R. solanacearum lysogenic strains. This suggested that φRSM3-ORF15 might be involved in the establishment and maintenance of a lysogenic state, as well as in phage immunity. Comparison of the putative repressor proteins and their binding sites within φRSM-related prophages provides insights into how these regulatory systems of filamentous phages have evolved and diverged in the R. solanacearum complex. In conclusion, φRSM phages represent a unique group of filamentous phages that are equipped with innate integration/excision (ORF14) and regulatory systems (ORF15).

General properties of transcriptional time series in Escherichia coli

Gene activity is described by the time-series of discrete, stochastic mRNA production events. This transcriptional time-series exhibits intermittent, bursty behavior. One consequence of this temporal intricacy is that gene expression can be tuned by varying different features of the time-series. What schemes for varying the transcriptional time-series are observed in the cell? Are the observed properties of these time-series optimized for cellular function? To address these questions, we characterize mRNA copy-number statistics at single-molecule resolution from multiple Escherichia coli promoters. We find that the degree of burstiness depends only on the gene expression level, while being independent of the details of gene regulation. The observed behavior is explained by the underlying variation in the duration of bursting events. Using information theory, we find that the properties of the transcriptional time series allow the cell to efficiently map the extracellular concentration of inducer molecules to intracellular levels of mRNA and proteins.

Lysogen stability is determined by the frequency of activity bursts from the fate-determining gene

Bacterial lysogeny serves as a simple paradigm for cell differentiation.

We characterize the activity of the fate-determining genes, cI and cro, with single-molecule resolution.

Stability of the lysogenic state is found to depend in a simple manner on the frequency of activity bursts from cI.

The ability of living cells to maintain an inheritable memory of their gene-expression state is key to cellular differentiation. Bacterial lysogeny serves as a simple paradigm for long-term cellular memory. In this study, we address the following question: in the absence of external perturbation, how long will a cell stay in the lysogenic state before spontaneously switching away from that state? We show by direct measurement that lysogen stability exhibits a simple exponential dependence on the frequency of activity bursts from the fate-determining gene, cI. We quantify these gene-activity bursts using single-molecule-resolution mRNA measurements in individual cells, analyzed using a stochastic mathematical model of the gene-network kinetics. The quantitative relation between stability and gene activity is independent of the fine details of gene regulation, suggesting that a quantitative prediction of cell-state stability may also be possible in more complex systems.

A comparative three-dimensional model of the carboxy-terminal domain of the lambda repressor and its use to build intact repressor tetramer models bound to adjacent operator sites

A model for residues 93-236 of the lambda repressor (1gfx) was predicted, based on the UmuD(') crystal structure, as part of four intact repressor molecules bound to two adjacent operator sites. The structure of region 136-230 in 1gfx was found to be nearly identical to the independently determined crystal structure of the 132-236 fragment, 1f39, released later by the PDB. Later, two more tetrameric models of the lambda repressor tetramer bound to two adjacent operator sites were constructed by us; in one of these, 1j5g, the N-domain and C-domain coordinates and hence monomer-monomer and dimer-dimer interactions are almost the same as in 1gfx, but the structure of the linker region is partly based on the linker region of the LexA dimer in 1jhe; in the other, 1lwq, the crystalline tetramer for region 140-236 has been coopted from the crystal structure deposited in 1kca, the operator DNA and N-domain coordinates of which are same as those in 1gfx and 1j5g, but the linker region is partly based on the LexA dimer structures 1jhe and 1jhh. Monomer-monomer interactions at the same operator site are stabilized by exposed hydrophobic side chains in beta-strands while cooperative interactions are mostly confined to beta(6) and some adjacent residues in both 1gfx and 1j5g. Mutational data, existence of a twofold axis relating two C-domains within a dimer, and minimization of DNA distortion between adjacent operator sites allow us to roughly position the C-domain with respect to the N-domain for both 1gfx and 1j5g. The study correlates these models with functional, biochemical, biophysical, and immunological data on the repressor in the literature. The oligomerization mode observed in the crystal structure of 132-236 may not exist in the intact repressor bound to the operator since it is shown to contradict several published biochemical data on the intact repressor.

Non-standard insulin design: structure-activity relationships at the periphery of the insulin receptor

The design of insulin analogues has emphasized stabilization or destabilization of structural elements according to established principles of protein folding. To this end, solvent-exposed side-chains extrinsic to the receptor-binding surface provide convenient sites of modification. An example is provided by an unfavorable helical C-cap (Thr(A8)) whose substitution by favorable amino acids (His(A8) or Arg(A8)) has yielded analogues of improved stability. Remarkably, these analogues also exhibit enhanced activity, suggesting that activity may correlate with stability. Here, we test this hypothesis by substitution of diaminobutyric acid (Dab(A8)), like threonine an amino acid of low helical propensity. The crystal structure of Dab(A8)-insulin is similar to those of native insulin and the related analogue Lys(A8)-insulin. Although no more stable than native insulin, the non-standard analogue is twice as active. Stability and affinity can therefore be uncoupled. To investigate alternative mechanisms by which A8 substitutions enhance activity, multiple substitutions were introduced. Surprisingly, diverse aliphatic, aromatic and polar side-chains enhance receptor binding and biological activity. Because no relationship is observed between activity and helical propensity, we propose that local interactions between the A8 side-chain and an edge of the hormone-receptor interface modulate affinity. Dab(A8)-insulin illustrates the utility of non-standard amino acids in hypothesis-driven protein design.

Activities of monomeric insulin analogs at position A8 are uncorrelated with their thermodynamic stabilities

Previous studies have demonstrated that the potency and thermodynamic stability of human insulin are enhanced in concert by substitution of Thr(A8) by arginine or histidine. These surface substitutions stabilize the N-terminal alpha-helix of the A chain, a key element of hormone-receptor recognition. Does enhanced stability necessarily imply enhanced activity? Here, we test by structure-based mutagenesis the relationship between the stability and activity of the hormone. To circumvent confounding effects of insulin self-association, A chain analogs were combined with a variant B chain (Asp(B10), Lys(B28), and Pro(B29) (DKP)) to create a monomeric template. Five analogs were obtained by chain combination; disulfide pairing proceeded in each case with native yield. CD and (1)H NMR spectra of the DKP analogs are essentially identical to those of DKP-insulin, indicating a correspondence of structures. Receptor binding affinities were determined by competitive displacement of (125)I-insulin from human placental membranes. Thermodynamic stabilities were measured by CD titration; unfolding was monitored as a function of guanidine concentration. In this broader collection of analogs receptor binding affinities are uncorrelated with stability. We suggest that receptor binding affinities of A8 analogs reflect local features of the hormone-receptor interface rather than the stability of the free hormone or the intrinsic C-capping propensity of the A8 side chain.

A chimeric activator of transcription that uses two DNA-binding domains to make simultaneous contact with pairs of recognition sites

Many well-known transcriptional regulatory proteins are composed of at least two independently folding domains and, typically, only one of these is a DNA-binding domain. However, some transcriptional regulators have been described that have more than one DNA-binding domain. Regulators with a single DNA-binding domain often bind co-operatively to the DNA in homotypic or heterotypic combinations, and two or more DNA-binding domains of a single regulatory protein can also bind co-operatively to suitably positioned recognition sequences. Here, we examine the behaviour of a chimeric activator of transcription with two different DNA-binding domains, that of the bacteriophage lambda cI protein and that of the Escherichia coli cyclic AMP receptor protein. We show that these two DNA-binding moieties, when present in the same molecule, can bind co-operatively to a pair of cognate recognition sites located upstream of a test promoter, thereby permitting the chimera to function as a particularly strong activator of transcription from this promoter. Our results show how such a bivalent DNA-binding protein can be used to regulate transcription differentially from promoters that bear either one or both recognition sites.

Dual regulatory control of a particle maturation function of bacteriophage P1

A unique arrangement of promoter elements was found upstream of the bacteriophage P1 particle maturation gene (mat). A P1-specific late-promoter sequence with conserved elements located at positions -22 and -10 was expected from the function of the gene in phage morphogenesis. In addition to a late-promoter sequence, a -35 element and an operator sequence for the major repressor protein, C1, were found. The -35 and -10 elements constituted an active Escherichia coli sigma(70) consensus promoter, which was converted into a P1-regulated early promoter by the superimposition of a C1 operator. This combination of early- and late-promoter elements regulates and fine-tunes the expression of the particle maturation gene. During lysogenic growth the gene is turned off by P1 immunity functions. Upon induction of lytic growth, the expression of mat starts simultaneously with the expression of other C1-regulated P1 early functions. However, while most of the latter functions are downregulated during late stages of lytic growth the expression of mat continues throughout the entire lytic growth cycle of bacteriophage P1. Thus, the maturation function has a head start on the structural components of the phage particle.

Mechanism for a transcriptional activator that works at the isomerization step

Transcriptional activators in prokaryotes have been shown to stimulate different steps in the initiation process including the initial binding of RNA polymerase (RNAP) to the promoter and a postbinding step known as the isomerization step. Evidence suggests that activators that affect initial binding can work by a cooperative binding mechanism by making energetically favorable contacts with RNAP, but the mechanism by which activators affect the isomerization step is unclear. A well-studied example of an activator that normally exerts its effect exclusively on the isomerization step is the bacteriophage lambda cI protein (lambdacI), which has been shown genetically to interact with the C-terminal region of the final sigma(70) subunit of RNAP. We show here that the interaction between lambdacI and final sigma can stimulate transcription even when the relevant portion of final sigma is transplanted to another subunit of RNAP. This activation depends on the ability of lambdacI to stabilize the binding of the transplanted final sigma moiety to an ectopic -35 element. Based on these and previous findings, we discuss a simple model that explains how an activator's ability to stabilize the binding of an RNAP subdomain to the DNA can account for its effect on either the initial binding of RNAP to a promoter or the isomerization step.

Protein-DNA recognition complexes: conservation of structure and binding energy in the transition state

This paper considers how enzymes that catalyze reactions at specific DNA sites have been engineered to overcome the problem of competitive inhibition by excess nonspecific binding sites on DNA. The formation of a specific protein-DNA recognition complex is discussed from both structural and thermodynamic perspectives, and contrasted with formation of nonspecific complexes. Evidence (from EcoRI and BamHI endonucleases) is presented that a wide variety of perturbations of the DNA substrate alter binding free energy but do not affect the free energy of activation for the chemical step; that is, many energetic factors contribute equally to the recognition complex and the transition-state complex. This implies that the specific recognition complex bears a close resemblance to the transition-state complex, such that very tight binding to the recognition site on the DNA substrate does not inhibit catalysis, but instead provides energy that is efficiently utilized along the path to the transition state. It is suggested that this view can be usefully extended to "noncatalytic" site-specific DNA-binding proteins like transcriptional activators and general transcription factors.

One nanosecond molecular dynamics simulation of the N-terminal domain of the lambda repressor protein

We have carried out molecular dynamics simulation of the N-terminal domain of the lambda repressor protein in a surrounding environment including explicit waters and ions. We observe two apparent dynamics substates in the nanosecond protein simulation, the transition occurring around 500 ps. The existence of these two apparent substates results from a high flexibility of the arm in each monomer, a relative flexibility of both arms with respect to each other, and a relative displacement of the recognition helices from 30 to 40 A of interhelical distance. Many amino acid residues, including those involved in DNA recognition, undergo a simultaneous transition in their side-chain conformations, consistent with the relationship between side-chain conformation and secondary structural elements, as observed in protein crystal structures. This result suggests plausible conformational changes experienced by the protein upon DNA binding. On the whole, the non-consensus monomer appears to be more flexible than its consensus counterpart.

The genetic switch for the regulatory pathway of Lactobacillus plantarum phage (phi)g1e: characterization of the promoter P(L), the repressor gene cpg, and the cpg-encoded protein Cpg in Escherichia coli

The structural and functional features of the approximately 530 bp P(L)/Gb5-Gb6-cpg-Gb7 region (P(L) overlaps Gb5) for the lysogenic pathway of L. plantarum phage (phi)gle were investigated using the cat gene of E. coli plasmid pKK232-8 as a reporter. In E. coli XL1-Blue, a recombinant plasmid pKPL2 (cat under P(L)/Gb5-Gb6) exhibited distinct CAT activity, whereas the activity of pKPLCP1 (cat under P(L)/Gb5-Gb6-cpg) was only marginal. When pKPL2 was coexistent with a compatible derivative of plasmid pACYC177 carrying P(L)/Gb5-Gb6-cpg, the CAT activity was declined to the level of pKPLCP1. On the other hand, the cpg-encoded protein Cpg was overproduced in E. coli under P(T7). The molecular mass of the purified Cpg (14.5 kDa on a SDS gel) corresponded well with that (15.1 kDa) predicted from the DNA sequence. Gel-shift and footprinting assays demonstrated that Cpg selectively binds to about 25 bp bases centered around the GATAC-box (from 1 to 7). Moreover, protein crosslinking experiments using glutaraldehyde showed that Cpg most likely functions as a dimeric form. Thus, the present results indicate that Cpg probably represses P(L) through binding to the operator GATAC-box(es), and the P(L)/cpg region might participate in the lysogenic pathway.

Electrostatic contributions to the binding free energy of the lambdacI repressor to DNA

A model based on the nonlinear Poisson-Boltzmann (NLPB) equation is used to study the electrostatic contribution to the binding free energy of the lambdacI repressor to its operator DNA. In particular, we use the Poisson-Boltzmann model to calculate the pKa shift of individual ionizable amino acids upon binding. We find that three residues on each monomer, Glu34, Glu83, and the amino terminus, have significant changes in their pKa and titrate between pH 4 and 9. This information is then used to calculate the pH dependence of the binding free energy. We find that the calculated pH dependence of binding accurately reproduces the available experimental data over a range of physiological pH values. The NLPB equation is then used to develop an overall picture of the electrostatics of the lambdacI repressor-operator interaction. We find that long-range Coulombic forces associated with the highly charged nucleic acid provide a strong driving force for the interaction of the protein with the DNA. These favorable electrostatic interactions are opposed, however, by unfavorable changes in the solvation of both the protein and the DNA upon binding. Specifically, the formation of a protein-DNA complex removes both charged and polar groups at the binding interface from solvent while it displaces salt from around the nucleic acid. As a result, the electrostatic contribution to the lambdacI repressor-operator interaction opposes binding by approximately 73 kcal/mol at physiological salt concentrations and neutral pH. A variety of entropic terms also oppose binding. The major force driving the binding process appears to be release of interfacial water from the protein and DNA surfaces upon complexation and, possibly, enhanced packing interactions between the protein and DNA in the interface. When the various nonelectrostatic terms are described with simple models that have been applied previously to other binding processes, a general picture of protein/DNA association emerges in which binding is driven by the nonpolar interactions, whereas specificity results from electrostatic interactions that weaken binding but are necessary components of any protein/DNA complex.

Amino acid-amino acid contacts at the cooperativity interface of the bacteriophage lambda and P22 repressors

The bacteriophage lambda repressor and its relatives bind cooperatively to adjacent as well as artificially separated operator sites. This cooperativity is mediated by a protein-protein interaction between the DNA-bound dimers. Here we use a genetic approach to identify two pairs of amino acids that interact at the dimer-dimer interface. One of these pairs is nonconserved in the aligned sequences of the lambda and P22 repressors; we show that a lambda repressor variant bearing the P22 residues at these two positions interacts specifically with the P22 repressor. The other pair consists of a conserved ion pair; we reverse the charges at these two positions and demonstrate that, whereas the individual substitutions abolish the interaction of the DNA-bound dimers, these changes in combination restore the interaction of both lambdacI and P22c2 dimers.

Conversion of the omega subunit of Escherichia coli RNA polymerase into a transcriptional activator or an activation target

Evidence obtained in both eukaryotes and prokaryotes indicates that arbitrary contacts between DNA-bound proteins and components of the transcriptional machinery can activate transcription. Here we demonstrate that the Escherichia coli omega protein, which copurifies with RNA polymerase, can function as a transcriptional activator when linked covalently to a DNA-binding protein. We show further that omega can function as an activation target when this covalent linkage is replaced by a pair of interacting polypeptides fused to the DNA-binding protein and to omega, respectively. Our findings imply that the omega protein is associated with RNA polymerase holoenzyme in vivo, and provide support for the hypothesis that contact between a DNA-bound protein and any component of E. coli RNA polymerase can activate transcription.

Changes in the periplasmic linker and in the expression level affect the activity of ToxR and lambda-ToxR fusion proteins in Escherichia coli

In order to assess the potentiality of Vibrio cholerae ToxR protein and of bacteriophage lambda repressor as indicators of the dimerization of periplasmic proteins in Escherichia coli, we have constructed a series of plasmids encoding transmembrane fusion proteins. The amino-terminal part, containing the DNA binding domain of either ToxR or lambda repressor, is located in the cytoplasm and acts as reporter for dimerization. As models of periplasmic proteins we have used alkaline phosphatase (a dimer) and beta-lactamase (a monomer). Both the expression level and the distance between the transmembrane segment and the periplasmic protein substantially affect the activity of the reporter domains.

Activation of prokaryotic transcription through arbitrary protein-protein contacts

Many transcriptional activators in prokaryotes are known to bind near a promoter and contact RNA polymerase, but it is not clear whether a protein-protein contact between an activator and RNA polymerase is enough to activate gene transcription. Here we show that contact between a DNA-bound protein and a heterologous protein domain fused to RNA polymerase can elicit transcriptional activation; moreover, the strength of this engineered protein-protein interaction determines the amount of gene activation. Our results indicate that an arbitrary interaction between a DNA-bound protein and RNA polymerase can activate transcription. We also find that when the DNA-bound 'activator' makes contact with two different components of the polymerase, the effect of these two interactions on transcription is synergistic.

An aromatic stacking interaction between subunits helps mediate DNA sequence specificity: operator site discrimination by phage lambda cI repressor

Sequence specific DNA binding by regulatory proteins provides the basis for regulation of initiation of transcription. A great deal of progress has been made toward understanding sequence specific recognition by individual protein subunits. An additional level of control that needs to be understood is that due to coupling between the subunits of oligomeric regulatory proteins. An example is the bacteriophage lambda cI repressor, a dimeric protein that regulates the lysogenic to lytic genetic switch of the phage. Two levels of specificity are critical to this regulation. First, like all transcriptional regulators, dimers distinguish operator from nonspecific DNA. Direct readout of the DNA sequence by the recognition helix is considered the well understood mechanism for this. However, differential affinity for O(R)1, O(R)2 and O(R)3 is equally critical to the switch because it mediates opposing regulation of divergent promoters. Site specificity at this second level is less well understood. Conformational adaptation by both the repressor and the different operators appears to be important. To evaluate how subunit-subunit interactions are involved in this process, we investigated the effects on both dimer stability and operator binding of amino acid substitutions at the contacts between the symmetrically related helices-5 in the dimer interface. Substitutions for Tyr88 alter dimer stability and greatly perturb differential operator affinity, but generally do not affect operator versus non-operator specificity. The pattern of these effects suggests that the geometry of the face-to-face aromatic stacking interaction between symmetrically related Tyr88 in each subunit, a group in the dimer interface but far removed from the DNA binding interface, plays a critical role in operator discrimination. Conformational changes in the tertiary structure of the subunits appears to be involved. By contrast, the significant effect of I84S substitution is to greatly decrease affinity for all three operators. Presumably, the altered packing of the dimer interface causes a quarternary structural change that moves the two helix-turn-helix motifs out of register with successive DNA major grooves.

Genome structure of the Lactobacillus temperate phage phi g1e: the whole genome sequence and the putative promoter/repressor system

The complete genome sequence of a Lactobacillus temperate phage phi g1e was established. The double-stranded DNA is composed of 42,259 bp, and encodes for sixty-two possible open reading frames (ORF) as well as several potential regulatory sequences. Based on comparative analysis with other related proteins of the Lactobacillus and Lactococcus phages as well as the Escherichia coli phages (such as lambda), functions were putatively assigned to several phi g1e ORFs: cng and cpg (encoding for repressors), hel (helicase), ntp (NTPase), and several ORFs (e.g., minor capsid proteins). An about 1000-bp DNA region of phi g1e containing cpg and cng was inferred to function as a promoter/repressor system for the phi g1e lysogenic and lytic pathway.

Contribution of a tyrosine side chain to ribonuclease A catalysis and stability

An intricate architecture of covalent bonds and noncovalent interactions appear to position the side chain of Lys 41 properly within the active site of bovine pancreatic ribonuclease A (RNase A). One of these interactions arises from Tyr 97, which is conserved in all 41 RNase A homologues of known sequence. Tyr 97 has a solvent-inaccessible side chain that donates a hydrogen bond to the main-chain oxygen of Lys 41. Here, the role of Tyr 97 was examined by replacing Tyr 97 with a phenylalanine, alanine, or glycine residue. All three mutant proteins have diminished catalytic activity, with the value of Kcat being perturbed more significantly than that of Km. The free energies with which Y97F, Y97A, and Y97G RNase A bind to the rate-limiting transition state during the cleavage of poly(cytidylic acid) are diminished by 0.74, 3.3, and 3.8 kcal/mol, respectively. These results show that even though Tyr 97 is remote from the active site, its side chain contributes to catalysis. The role of Tyr 97 in the thermal stability of RNase A is large. The conformational free energies of native Y97F, Y97A, and Y97G RNase A are decreased by 3.54, 12.0, and 11.7 kcal/mol, respectively. The unusually large decrease in stability caused by the Tyr-->Phe mutation could result from a decrease in the barrier to isomerization of the Lys 41-Pro 42 peptide bond.

Identification of prophage genes expressed in lysogens of the Lactococcus lactis bacteriophage BK5-T

Bacteriophage BK5-T is a small isometric-headed temperate phage that infects Lactococcus lactis subsp. cremoris. Northern (RNA) analysis of mRNA produced by lysogenic strains containing BK5-T prophage revealed four major BK5-T transcripts that are 0.8, 1.3, 1.8, and 1.8 kb in size and enabled a transcription map of the prophage genome to be prepared. The position and size of each transcript corresponded closely to the position and size of open reading frames predicted from the nucleotide sequence of BK5-T. Analysis of the transcripts suggested that one of them was derived from the gene encoding the BK5-T integrase and another was from the gene encoding the BK5-T homolog of the lambda cI repressor. Computer analysis of the nucleotide sequence upstream of the BK5-T cI homolog predicted the presence of a pair of divergent promoters and three inverted repeat sequences, features characteristic of temperature-phage immunity regions. By analogy with lambda, the three inverted repeat sequences could be binding sites for cI or Cro homologs and the two divergent promoters could initiate transcription through the BK5-T equivalents of cI and cro.

Submillisecond folding of monomeric lambda repressor

The folding kinetics of a truncated form of the N-terminal domain of phage lambda repressor [lambda 6-85] has been investigated by using the technique of dynamic NMR. lambda 6-85 has been shown previously to fold in a purely two-state fashion. This allows the determination of folding and unfolding rates from simulation of the exchange-broadened aromatic resonances of Tyr-22. The folding kinetics were determined over a range of 1.35 to 3.14 M urea. The urea dependence of both folding and unfolding rate constants is exponential, suggesting that the rate-determining step is invariant at the urea concentrations studied. The folding and unfolding rates extrapolated to 0 M urea at 37 degrees C are 3600 +/- 400 s-1 and 27 +/- 6 s-1, respectively. The observed lambda 6-85 folding rate constant exceeds that of other fast-folding globular proteins by a factor of 14-54. The urea dependence of the folding and unfolding rate constants suggests that the transition state of the rate-determining step is considerably more exposed to solvent than previously studied protein-folding transition states. The surprising rapidity of lambda 6-85 folding and unfolding may be the consequence of its all-helical secondary structure. These kinetic results clearly demonstrate that all of the fundamental events of protein folding can occur on the submillisecond time scale.

Structure and stability of monomeric lambda repressor: NMR evidence for two-state folding

The absence of equilibrium intermediates in protein folding reactions (i.e., two-state folding) simplifies thermodynamic and kinetic analyses but is difficult to prove rigorously. We demonstrate a sensitive method for detecting partially folded species based on using proton chemical shifts as local probes of structure. The coincidence of denaturation curves for probes throughout the molecule is a particularly stringent test for two-state folding. In this study we investigate a new form of the N-terminal domain of bacteriophage lambda repressor consisting of residues 6-85 (lambda 6-85) using nuclear magnetic resonance (NMR) and circular dichroism (CD). This truncated version lacks the residues required for dimerization and is monomeric under the conditions used for NMR. Heteronuclear NMR was used to assign the 1H, 15N, and backbone 13C resonances. The secondary and tertiary structure of lambda 6-85 is very similar to that reported for the crystal structure of the DNA-bound 1-92 fragment [Beamer, L. J., and Pabo, C. O. (1992) J. Mol. Biol. 227, 177-196], as judged by analysis of chemical shifts, amide hydrogen exchange, amide-alpha coupling constants, and nuclear Overhauser enhancements. Thermal and urea denaturation studies were conducted using the chemical shifts of the four aromatic side chains as local probes and the CD signal at 222 nm as a global probe. Plots of the fraction denatured versus denaturant concentration obtained from these studies are identical for all probes under all conditions studied. This observation provides strong evidence for two-state folding, indicating that there are no populated intermediates in the folding of lambda 6-85.

The interactions of Escherichia coli trp repressor with tryptophan and with an operator oligonucleotide. NMR studies using selectively 15N-labelled protein

The effects of the binding of the corepressor L-tryptophan and an operator oligonucleotide to Escherichia coli trp repressor have been studied, using selective 15N labelling to permit observation of the backbone amide resonances of 50 of the 107 residues of the protein monomer. Repressor molecules selectively labelled in turn with [15N]alanine, [15N]glutamate, [15N]isoleucine, [15N]leucine and [15N]methionine were prepared by isolating them from prototrophic E. coli cells grown in media containing a mixture of unlabelled and the appropriate 15N-enriched amino acids. Analysis of the heteronuclear correlation spectra of the labelled repressors shows the value of selective labelling in resolving the crosspeaks of, for example, the 19 leucine and 12 glutamate residues. All 50 residues studied show measurable changes in amide 1H and/or 15N chemical shift on the binding of tryptophan and/or the operator oligonucleotide, showing clearly that ligand binding has effects which are transmitted throughout almost the whole protein. Large chemical shift changes on ligand binding are seen in residues in the tryptophan binding site and in the 'helix-turn-helix' DNA-binding domain, but also in residues in helices C and F remote from the ligand binding sites. On operator binding there is selective broadening of the signals of residues in the N-terminal region of the protein and in the DNA-binding domain, perhaps reflecting a conformational equilibrium.

Interaction of RNA hairpins with the human U1A N-terminal RNA binding domain

The isolated 102 amino acid N-terminal RNA binding domain (RBD) of the human U1A protein specifically interacts with a short RNA hairpin containing the U1 snRNA stem/loop II sequence. This recognition is nucleotide-specific, for substitutions of critical nucleotides in the RNA loop decrease binding affinity up to 10(6)-fold, as measured by nitrocellulose filter binding experiments. The magnitude of the loss of binding free energy with single-nucleotide substitution in the conserved GCA sequence suggests that the interaction between the RBD and RNA occurs through a number of interdependent specific contacts in the complex. 13C and 15N NMR experiments, using isotopically-labeled RNA together with unlabeled protein, show that the chemical shifts of many protons from the bound RNA are substantially different from those of the free RNA, especially in the loop region of the hairpin. All these data suggest that there is a conformational change in the RNA upon formation of the RBD-RNA complex.

In vitro selection of zinc fingers with altered DNA-binding specificity

We have used random mutagenesis and phage display to alter the DNA-binding specificity of Zif268, a transcription factor that contains three zinc finger domains. Four residues in the helix of finger 1 of Zif268 that potentially mediate DNA binding were identified from an X-ray structure of the Zif268-DNA complex. A library was constructed in which these residues were randomly mutated and the Zif268 variants were fused to a truncated version of the gene III coat protein on the surface of M13 filamentous phage particles. The phage displayed the mutant proteins in a monovalent fashion and were sorted by repeated binding and elution from affinity matrices containing different DNA sequences. When the matrix contained the natural nine base pair operator sequence 5'-GCG-TGG-GCG-3', native-like zinc fingers were isolated. New finger 1 variants were found by sorting with two different operators in which the singly modified triplets, GTG and TCG, replaced the native finger 1 triplet, GCG. Overall, the selected finger 1 variants contained a preponderance of polar residues at the four sites. Interestingly, the net charge of the four residues in any selected finger never derived more that one unit from neutrality despite the fact that about half the variants contained three or four charged residues over the four sites. Measurements of the dissociation constants for two of these purified finger 1 variants by gel-shift assay showed their specificities to vary over a 10-fold range, with the greatest affinity being for the DNA binding site for which they were sorted.(ABSTRACT TRUNCATED AT 250 WORDS)

Specificity determinants for the interaction of lambda repressor and P22 repressor dimers

The related phage lambda and phage P22 repressors each bind cooperatively to adjacent and separated operator sites, an interaction that involves a pair of repressor dimers. The specificities of these interactions differ: Each dimer interacts with its own type but not with dimers of the heterologous repressor. The two repressors exhibit significant amino acid sequence homology in their carboxy-terminal domains, which are responsible for both dimer formation and the dimer-dimer interaction. Here, we identify a collection of amino acid substitutions that disrupt the protein-protein interaction of DNA-bound lambda repressor dimers and show that several of these substitutions have the same effect when introduced at the corresponding positions of P22 repressor. We use this information to construct a variant of the lambda repressor bearing only six non-wild-type amino acids that has a switched specificity; that is, it binds cooperatively with P22 repressor, but not with wild-type lambda repressor. These results identify a series of residues that determine the specificities of the two interactions.

The crystal structure of a mutant protein with altered but improved hydrophobic core packing

The dense packing observed in protein interiors appears to be crucial for stabilizing the native structure--even subtle internal substitutions are usually destabilizing. Thus, steric complementarity of core residues is thought to be an important criterion for "inverse folding" predictive methods, which judge whether a newly determined sequence is consistent with any known folds. A major problem in the development of useful core packing evaluation algorithms, however, is that there are occasional mutations that are predicted to disrupt native packing but that yield an equally or more stable protein. We have solved the crystal structure of such a variant of lambda repressor, which, despite having three larger core substitutions, is more stable than the wild type. The structure reveals that the protein accommodates the potentially disruptive residues with shifts in its alpha-helical arrangement. The variant is apparently more stable because its packing is improved--the core has a higher packing density and little geometric strain. These rearrangements, however, cause repositioning of functional residues, which result in reduced DNA binding activity. By comparing these results with the predictions of two core packing algorithms, it is clear that the protein possesses a relatively high degree of main-chain flexibility that must be accounted for in order to predict the full spectrum of compatible core sequences. This study also shows how, in protein evolution, a particular set of core residue identities might be selected not because they provide optimal stability but because they provide sufficient stability in addition to the precise structure required for optimal activity.

Separation of the complex DNA binding domain of EBNA-1 into DNA recognition and dimerization subdomains of novel structure

EBNA-1 is essential for replication of the latent episomal form of the Epstein-Barr virus genome and is involved in regulation of viral latency promoters. EBNA-1 activity is mediated through direct DNA binding. The DNA binding and dimerization functions of EBNA-1 have previously been located to a carboxy-terminal domain, amino acids (aa) 459 to 607. To identify and define the subdomains for these two functions, we created an extensive series of deletions and point mutations in an EBNA-1 (aa 408 to 641) background. The ability of the EBNA-1 mutants to heterodimerize with a wild-type EBNA-1 (aa 459 to 641) Immunoprecipitation assays with a monoclonal antibody, EBNA.OT1x, that recognizes EBNA-1 (aa 408 to 641) but not EBNA-1 (aa 459 to 641). These experiments revealed that mutations affecting dimerization occurred over two separate regions, aa 501 to 532 and aa 554 to 598. DNA binding was tested in mobility shift assays against a panel of oligonucleotide-binding sites. Dimerization was a prerequisite for DNA binding. The DNA recognition domain was localized to a separate region, aa 459 to 487, upstream of the dimerization domain. EBNA-1 variants carrying substitutions at aa 467 and 468 and at aa 477 gave a pattern of binding to mutant oligonucleotide probes that implicates these particular amino acids in DNA recognition. EBNA-1 appears to utilize novel mechanisms for both DNA recognition and dimerization since neither domain conforms to previously described structural motifs.

A mutation outside the two zinc fingers of ADR1 can suppress defects in either finger

A second-site mutation that restored DNA binding to ADR1 mutants altered at different positions in the two zinc fingers was identified. This mutation (called IS1) was a conservative change of arginine 91 to lysine in a region amino terminal to the two zinc fingers and known from previous experiments to be necessary for DNA binding. IS1 increased binding to the UAS1 sequence two- to sevenfold for various ADR1 mutants and twofold for wild-type ADR1. The change of arginine 91 to glycine decreased binding twofold, suggesting that this arginine is involved in DNA binding in the wild-type protein. The increase in binding by IS1 did not involve protein-protein interactions between the two ADR1 monomers, nor did it require the presence of the sequences flanking UAS1. However, the effect of IS1 was influenced by the sequence of the first finger, suggesting that interactions between the region amino terminal to the fingers and the fingers themselves could exist. A model for the role of the amino-terminal region based on these results and sequence homologies with other DNA-binding motifs is proposed.

A mutation outside the two zinc fingers of ADR1 can suppress defects in either finger

A second-site mutation that restored DNA binding to ADR1 mutants altered at different positions in the two zinc fingers was identified. This mutation (called IS1) was a conservative change of arginine 91 to lysine in a region amino terminal to the two zinc fingers and known from previous experiments to be necessary for DNA binding. IS1 increased binding to the UAS1 sequence two- to sevenfold for various ADR1 mutants and twofold for wild-type ADR1. The change of arginine 91 to glycine decreased binding twofold, suggesting that this arginine is involved in DNA binding in the wild-type protein. The increase in binding by IS1 did not involve protein-protein interactions between the two ADR1 monomers, nor did it require the presence of the sequences flanking UAS1. However, the effect of IS1 was influenced by the sequence of the first finger, suggesting that interactions between the region amino terminal to the fingers and the fingers themselves could exist. A model for the role of the amino-terminal region based on these results and sequence homologies with other DNA-binding motifs is proposed.

Refined 1.8 A crystal structure of the lambda repressor-operator complex

The crystal structure of the lambda repressor-operator complex has been refined to an R-factor of 18.9% at 1.8 A resolution. This refinement, using data collected at low temperature, has revealed the structure of the N-terminal arm and shows that the interactions of repressor with the two halves of the pseudo-symmetric operator site are significantly different. The two halves of the complex are most similar near the outer edge of the operator site (in a region where the lambda and 434 repressors make similar contacts), but they become increasingly different toward the center of the operator. There are striking differences near the center of the site where it appears that the arm makes significant contacts to only one half of the DNA site. This suggested a new way of aligning the operator sites in phage lambda. The high resolution structure confirms many of the previously noted features of the complex, but also reveals a number of new protein-DNA contacts. It also gives a better view of the extensive H-bonding networks that couple contacts made by different residues and different regions of the protein, and reveals important new details about the helix-turn-helix (HTH) region, and the positions of many water molecules in the complex.

Contacts between Tet repressor and tet operator revealed by new recognition specificities of single amino acid replacement mutants

We have analyzed the DNA binding properties of Tet-repressor mutants with single amino acid residue replacements at eight positions within the alpha-helix-turn-alpha-helix DNA-binding motif. A saturation mutagenesis of Gln38, Pro39, Thr40, Tyr42, Trp43 and His44 in the second alpha-helix was performed; in addition, several substitutions of Thr27 and Arg28 in the first alpha-helix were constructed. The abilities of these mutant repressors to bind a set of 16 operator variants were determined and revealed 23 new binding specificities. All repressor mutants with DNA-binding activity were inducible by tetracycline, while mutants lacking binding activity were trans-dominant over the wild-type. All mutant proteins were present at the same intracellular steady-state concentrations as the wild-type. These results suggest the structural integrity of the mutant repressors. On the basis of the new recognition specificities, five contacts between a repressor monomer and each operator half-site and the chemical nature of these repressor-operator interactions are proposed. We suggest that Arg28 contacts guanine of the G.C base-pair at operator position 2 with two H-bonds, Gln38 binds adenine of the A.T base-pair at position 3 with two H-bonds, and the methyl group of Thr40 participates in a van der Waals' contact with cytosine of the G.C base-pair at position 6 of tet operator. A previously unrecognized type of interaction is proposed for Pro39, which inserts its side-chain between the methyl groups of the thymines of T.A and A.T base-pairs at positions 4 and 5. Computer modeling of these proposed contacts reveals that they are possible using the canonical structures of the helix-turn-helix motif and B-DNA. These contacts suggest an inverse orientation of the Tet repressor helix-turn-helix with respect to the operator center as compared with non-inducible repressor-operator complexes, and are supported by similar contacts of other repressor-operator complexes.

How Lac repressor finds lac operator in vitro

Filter-binding and gel mobility shift assays were used to analyse the kinetics of the interaction of Lac repressor with lac operator. A comparison of the two techniques reveals that filter-binding assays with tetrameric Lac repressor have often been misinterpreted. It has been assumed that all complexes of Lac repressor and lac operator DNA bind with equal affinity to nitrocellulose filters. This assumption is wrong. Sandwich or loop complexes where two lac operators bind to one tetrameric Lac repressor are not or are only badly retained on nitrocellulose filters under normal conditions. Taking this into account, dimeric and tetrameric Lac repressor do not show any DNA-length dependence of their association and dissociation rate constants when they bind to DNA fragments smaller than 2455 base-pairs carrying a single symmetric ideal lac operator. A ninefold increased association rate to ideal lac operator on lambda DNA is observed for tetrameric but not dimeric Lac repressor. It is presumably due to intersegment transfer involving lac operator-like sequences.

Binding of EBNA-1 to DNA creates a protease-resistant domain that encompasses the DNA recognition and dimerization functions

The Epstein-Barr virus nuclear antigen EBNA-1 is essential for replication of the viral DNA during latency. EBNA-1 binds as a dimer to palindromic recognition sequences within the plasmid origin of replication, ori-P. In this study, proteinase K susceptibility has been used to further characterize the DNA-binding domain of EBNA-1. Limited protease digestion of EBNA-1 (amino acids 408 to 641) generated a smaller DNA-binding species that had a degree of inherent protease resistance. When EBNA-1 was preincubated with a specific DNA probe, the protease resistance of the smaller binding species increased 100-fold, suggesting that the conformation of EBNA-1 changes on binding. The protease-resistant species comprised an 18-kDa polypeptide that was further cleaved at high levels of protease to 11- and 5.4-kDa products. A model of the proposed protease-resistant domain structure is presented. Constructions carrying serial, internal deletions across the 18-kDa domain were created. Each of the deletions perturbed dimerization ability and abolished DNA binding. These studies suggest that the DNA-binding and dimerization motifs of EBNA-1 lie within a conformationally discrete domain whose overall integrity is necessary for EBNA-1-DNA interaction.

In vivo analysis of the helix-turn-helix motif of the fushi tarazu homeo domain of Drosophila melanogaster

We report a systematic mutational analysis of the helix-turn-helix motif (HTH) of the fushi tarazu (ftz) homeo domain (HD) of Drosophila. We started out by testing the function of chimeric ftz proteins containing either a part of the Sex combs reduced (Scr) or the muscle segment homeobox (msh) HDs. By complementation tests in transgenic flies, cotransfection assays in cultured Drosophila cells and in vitro DNA-binding assays, we have found that the ftz activity is retained in the ftz-Scr chimera but is lost in the ftz-msh chimera, which is defective in binding to an Antennapedia (Antp)-class target site. Further studies with a series of back-mutants of the ftz-msh chimera have revealed that a set of class-specific DNA backbone-contacting residues in the HTH, particularly Arg-28 and Arg-43, are required for efficient target site recognition and, hence, full ftz activity both in vitro and in vivo.

Glucocorticoid receptor DNA-binding specificity is increased by the organization of DNA in nucleosomes

A DNA fragment containing glucocorticoid receptor binding sites in the mouse mammary tumor virus promoter was reconstituted in vitro with histones to form nucleosome cores, which become positioned on the DNA fragment in a sequence-specific manner. Glucocorticoid receptor binding to specific DNA sequences was analyzed by quantitative DNase I footprinting. The receptor interacted with surprisingly high affinity with one of the binding sites in the reconstituted promoter, although it was reduced by a factor of approximately 2 compared with the same site in protein-free DNA. By contrast, the affinity for random genomic nucleosomal sites was drastically reduced compared with histone-free DNA. Thus, reconstituting the promoter in vitro resulted in a 60- to 70-fold increase in binding specificity. Such an increase in selective binding may help to explain the ability of glucocorticoid receptor to effectively locate its target sites in chromatin.

Structural basis for the nucleic acid binding cooperativity of bacteriophage T4 gene 32 protein: the (Lys/Arg)3(Ser/Thr)2 (LAST) motif

To identify the functional residues of the N-terminal B region of bacteriophage T4 gene 32 protein involved in its cooperative binding to single-stranded nucleic acids, a process dependent on homotypic protein-protein interaction, we have studied the interaction of the protein with synthetic peptides containing different portions of this domain. Gel-permeation chromatography showed that a 6-acryloyl-2-dimethylaminonaphthalene (acrylodan)-labeled fluorescent peptide corresponding to the first 17 residues of gene 32 protein formed a complex with whole protein. The fluorescence was blue-shifted 14 nm upon interaction with intact protein, and somewhat less so (7-11 nm) with cleavage products of the protein lacking B domains. The intrinsic tryptophan fluorescence of whole and truncated protein was quenched by this peptide and by the nonderivatized peptide. The peptide bound tightly to truncated protein at both 0.015 and 0.44 M Na+, with a stoichiometry of 1:1. Similar tryptophan quenching or acrylodan blue shifts were obtained with peptides corresponding to residues 1-9 and 3-8, but not residues 1-4, 5-9, or 5-17, indicating that the essential amino acids are contained within positions 3-8, Lys-Arg-Lys-Ser-Thr-Ala. Several other DNA binding proteins contain a LAST motif with documented involvement of these residues in nucleic acid interaction. The amino acid and coding sequence of residues 110-114, a region proposed to be involved in nucleic acid binding, is virtually identical to that of residues 3-7. Based on these observations, we have formulated a model for the cooperative interactions of gene 32 protein with single-stranded nucleic acids.

Analysis and modulation of protein stability

Numerous site-directed mutagenesis experiments have provided new insights into the stabilizing role of the individual forces and interactions within a globular protein molecule. Some useful guidelines and procedures are now available for producing genetically more stable proteins. Examples are the introduction of disulfide bonds, ion-binding sites, salt bridges, hydrophobic residues or hydrogen bonds, and the improvement of hydrophobic packing or alpha-helix propensity. Moreover, it is now clearly recognized that thermophilic (and, in general, extremophilic) bacteria produce highly stable proteins and enzymes of practical interest.

The role of internal packing interactions in determining the structure and stability of a protein

Cassette mutagenesis has been used to investigate how internal packing interactions help to specify a protein's three-dimensional structure and stability. Three interacting residues in the hydrophobic core of the N-terminal domain of lambda repressor were randomized combinatorially. The randomization was restricted to the five amino acids Val, Leu, Ile, Met and Phe, thereby generating a sterically diverse set of core sequences composed solely of hydrophobic residues. We have isolated 78 of the 125 possible sequences generated by this randomization. Approximately 70% of the isolated sequences show some level of biological activity, and thus still carry sufficient information to encode the basic structure of lambda repressor. An assay based on the temperature dependence of activity in vivo has been used to estimate the relative activities and thermal stabilities of the set of mutants. In addition, nine mutants have been purified and their stabilities and DNA binding activities characterized in vitro. Of the 56 active sequences, only two, in addition to the wild-type, maintain the wild-type level of stability and activity. All three of these proteins satisfy stringent requirements for specifically shaped residues at each position. All of the remaining active sequences have reduced stabilities and/or reduced DNA binding affinities. These and previous results suggest that there are two levels of structural information encoded in core residues. At the first level, the basic structural information appears to reside largely in the hydrophobic character of these residues. The majority of sequences that simply maintain hydrophobicity at core positions are able to adopt the overall lambda repressor fold and maintain moderate stability. At the second, more detailed level, specific steric features of these residues and their packing interactions clearly act as important determinants of the protein's precise structure and stability. These results imply that many of the basic structural features of a protein could be predicted from relatively simple, degenerate sequence patterns.

A syntactic representation of units of genetic information--a syntax of units of genetic information

The experimental study of the molecular mechanisms involved in the regulation of operons and other units of genetic information (UGIs) has not been accompanied by a parallel effort in integrative approaches. Following a recently obtained mathematical justification in the search of a grammatical theory of regulation, in this paper it is shown that a syntactic representation of prokaryotic simple UGIs can help the search for general rules governing their organization and regulation. Molecular categories like promoter, operator, structural gene, etc, are the elements for this level of analysis. Based on diverse types of evidence available in the literature, a principle which establishes a strictly successive representation of UGIs regulated at the initiation of transcription is derived. This linear array of categories is not enough for proposing general rules, which can be obtained by identifying groups of these categories as clusters or syntactic categories. It is shown that the notion of syntactic categories is implicitly used in the classical definition of an operon. Based on a hierarchy of biological restrictions for the construction of UGIs, a grammatical principle that defines a hierarchical relation among activator regions, promoters and operators is proposed. These proposals are integrated in a grammar that accounts for simple positively and negatively regulated UGIs.