Combinations of the alpha-helix-turn-alpha-helix motif of TetR with respective residues from LacI or 434Cro: DNA recognition, inducer binding, and urea-dependent denaturation

Transcription Factor Based Small-Molecule Sensing with a Rapid Cell Phone Enabled Fluorescent Bead Assay

Recently, allosteric transcription factors (TFs) were identified as a novel class of biorecognition elements for in vitro sensing, whereby an indicator of the differential binding affinity between a TF and its cognate DNA exhibits dose-dependent responsivity to an analyte. Described is a modular bead-based biosensor design that can be applied to such TF-DNA-analyte systems. DNA-functionalized beads enable efficient mixing and spatial separation, while TF-labeled semiconductor quantum dots serve as bright fluorescent indicators of the TF-DNA bound (on bead) and unbound states. The prototype sensor for derivatives of the antibiotic tetracycline exhibits nanomolar sensitivity with visual detection of bead fluorescence. Facile changes to the sensor enable sensor response tuning without necessitating changes to the biomolecular affinities. Assay components self-assemble, and readout by eye or digital camera is possible within 5 minutes of analyte addition, making sensor use facile, rapid, and instrument-free.

Engineering DNA recognition and allosteric response properties of TetR family proteins by using a module-swapping strategy

The development of synthetic biological systems requires modular biomolecular components to flexibly alter response pathways. In previous studies, we have established a module-swapping design principle to engineer allosteric response and DNA recognition properties among regulators in the LacI family, in which the engineered regulators served as effective components for implementing new cellular behavior. Here we introduced this protein engineering strategy to two regulators in the TetR family: TetR (UniProt Accession ID: P04483) and MphR (Q9EVJ6). The TetR DNA-binding module and the MphR ligand-binding module were used to create the TetR-MphR. This resulting hybrid regulator possesses DNA-binding properties of TetR and ligand response properties of MphR, which is able to control gene expression in response to a molecular signal in cells. Furthermore, we studied molecular interactions between the TetR DNA-binding module and MphR ligand-binding module by using mutant analysis. Together, we demonstrated that TetR family regulators contain discrete and functional modules that can be used to build biological components with novel properties. This work highlights the utility of rational design as a means of creating modular parts for cell engineering and introduces new possibilities in rewiring cellular response pathways.

Simplified mechanistic models of gene regulation for analysis and design

Simplified mechanistic models of gene regulation are fundamental to systems biology and essential for synthetic biology. However, conventional simplified models typically have outputs that are not directly measurable and are based on assumptions that do not often hold under experimental conditions. To resolve these issues, we propose a ‘model reduction’ methodology and simplified kinetic models of total mRNA and total protein concentration, which link measurements, models and biochemical mechanisms. The proposed approach is based on assumptions that hold generally and include typical cases in systems and synthetic biology where conventional models do not hold. We use novel assumptions regarding the ‘speed of reactions’, which are required for the methodology to be consistent with experimental data. We also apply the methodology to propose simplified models of gene regulation in the presence of multiple protein binding sites, providing both biological insights and an illustration of the generality of the methodology. Lastly, we show that modelling total protein concentration allows us to address key questions on gene regulation, such as efficiency, burden, competition and modularity.

Explicit analytic equations for multimolecular thermal melting curves

The analysis of thermal melting curves requires the knowledge of equations for the temperature dependence of the relative fraction of folded and unfolded components. To implement these equations as standard tools for curve fitting, they should be as explicit as possible. From the van't Hoff formalism it is known that the equilibrium constant and hence the folded fraction is a function of the absolute temperature, the van't Hoff transition enthalpy, and the melting temperature. The work presented here is devoted to the mathematically self-contained derivation and the listing of explicit equations for the folded fraction as a function of the thermodynamic parameters in the case of arbitrary molecularities. Part of the results are known, others are new. It is in particular shown for the first time that the folded fraction is the composition of a universal function which depends solely on the molecularity and a dimensionless function which is governed by the concrete thermodynamic regime but is independent of the molecularity. The results will prove useful for extracting the thermodynamic parameters from experimental data on the basis of regression analysis. As supporting information, open-source Matlab scripts for the computer implementation of the equations are provided.

Tuning and controlling gene expression noise in synthetic gene networks

Synthetic gene networks can be used to control gene expression and cellular phenotypes in a variety of applications. In many instances, however, such networks can behave unreliably due to gene expression noise. Accordingly, there is a need to develop systematic means to tune gene expression noise, so that it can be suppressed in some cases and harnessed in others, e.g. in cellular differentiation to create population-wide heterogeneity. Here, we present a method for controlling noise in synthetic eukaryotic gene expression systems, utilizing reduction of noise levels by TATA box mutations and noise propagation in transcriptional cascades. Specifically, we introduce TATA box mutations into promoters driving TetR expression and show that these mutations can be used to effectively tune the noise of a target gene while decoupling it from the mean, with negligible effects on the dynamic range and basal expression. We apply mathematical and computational modeling to explain the experimentally observed effects of TATA box mutations. This work, which highlights some important aspects of noise propagation in gene regulatory cascades, has practical implications for implementing gene expression control in synthetic gene networks.

The induction of folding cooperativity by ligand binding drives the allosteric response of tetracycline repressor

Tetracycline (Tc) repressor (TetR) undergoes an allosteric transition upon interaction with the antibiotic, Tc, that abrogates its ability to specifically bind its operator DNA. In this work, by performing equilibrium protein unfolding experiments on wild-type TetR and mutants displaying altered allosteric responses, we have delineated a model to explain TetR allostery. In the absence of Tc, we show that the DNA-binding domains of this homodimeric protein are relatively flexible and unfold independently of the Tc binding/dimerization (TBD) domains. Once Tc is bound, however, the unfolding of the DNA binding domains becomes coupled to the TBD domains, leading to a large increase in DNA-binding domain stability. Noninducible TetR mutants display considerably less interdomain folding cooperativity upon binding to Tc. We conclude that the thermodynamic coupling of the TetR domains caused by Tc binding and the resulting rigidification of the DNA-binding domains into a conformation that is incompatible with DNA binding are the fundamental factors leading to the allosteric response in TetR. This allosteric mechanism can account for properties of the whole TetR family of repressors and may explain the functioning and evolution of other allosteric systems. Our model contrasts with the prevalent view that TetR populates two distinct conformations and that Tc causes a switch between these defined conformations.

Engineered Tet repressors with recognition specificity for the tetO-4C5G operator variant

We created a new DNA recognition specificity for tetracycline repressor (TetR) binding to the tet operator variant tetO-4C5G containing four bp exchanges compared to tetO. TetR variants created by doped oligonucleotide mutagenesis of residues in the DNA recognition helix yielded several mutants binding to tetO-4C5G. These variants contained exchanges of the amino acids at positions 36, 37, 39 and 42. The two amino acid exchanges in TetR E37A P39K are sufficient for tetO-4C5G specific binding. The E37A mutation increases the affinity of TetR for tetO variants and seems to be essential for binding to modified operator sequences. The Lys39 residue is in a position to directly contact the fourth and fifth bps of tetO thereby creating specificity for tetO-4C5G. Combinations of these mutations with others that lead to a reverse phenotype or altered inducer specificity yielded new TetR mutants with the respective combined activities. Single chain TetR variants were constructed that contain DNA reading heads with two different operator binding specificities. Specific binding of this TetR mutant to the respective mixed tetO-wt/4C5G variants containing one wild type and one double exchange operator half site was only accomplished at a low expression level of TetR variant, while cross-talk with other operator variants were observed at an elevated expression level. This observation emphasizes the importance of the transcription factor expression level for in vivo DNA binding specificity. These new TetR variants can be useful for multigene regulation systems.

Induction of single chain tetracycline repressor requires the binding of two inducers

In this article we report the in vivo and in vitro characterization of single chain tetracycline repressor (scTetR) variants in Escherichia coli. ScTetR is genetically and proteolytically stable and exhibits the same regulatory properties as dimeric TetR in E.coli. Urea-dependent denaturation of scTetR is independent of the protein concentration and follows the two-state model with a monophasic transition. Contrary to dimeric TetR, scTetR allows the construction of scTetR mutants, in which one subunit contains a defective inducer binding site while the other is functional. We have used this approach to establish that scTetR needs occupation of both inducer binding sites for in vivo and in vitro induction. Single mutations causing loss of induction in dimeric TetR lead to non-inducible scTetR when inserted into one half-side. The construction of scTetR H64K S135L S138I (scTetR(i2)) in which one half-side is specific for 4-dedimethylamino-anhydrotetracycline (4-ddma-atc) and the other for tetracycline (tc) leads to a protein which is only inducible by the mixture of tc and 4-ddma-atc. Fluorescence titration of scTetR(i2) with both inducers revealed distinct occupancy with each of these inducers yielding roughly a 1:1 stoichiometry of each inducer per scTetR(i2). The properties of this gain of function mutant clearly demonstrate that scTetR requires the binding of two inducers for induction of transcription.

Noise in transcription negative feedback loops: simulation and experimental analysis

Negative feedback loops have been invoked as a way to control and decrease transcriptional noise. Here, we have built three circuits to test the effect of negative feedback loops on transcriptional noise of an autoregulated gene encoding a transcription factor (TF) and a downstream gene (DG), regulated by this TF. Experimental analysis shows that self-repression decreases noise compared to expression from a non-regulated promoter. Interestingly enough, we find that noise minimization by negative feedback loop is optimal within a range of repression strength. Repression values outside this range result in noise increase producing a U-shaped behaviour. This behaviour is the result of external noise probably arising from plasmid fluctuations as shown by simulation of the network. Regarding the target gene of a self-repressed TF (sTF), we find a strong decrease of noise when repression by the sTF is strong and a higher degree of noise anti-correlation between sTF and its target. Simulations of the circuits indicate that the main source of noise in these circuits could come from plasmid variation and therefore that negative feedback loops play an important role in suppressing both external and internal noise. An important observation is that DG expression without negative feedback exhibits bimodality at intermediate TF repression values. This bimodal behaviour seems to be the result of external noise as it can only be found in those simulations that include plasmid variation.

Two-way interdomain signal transduction in tetracycline repressor

The transcription of genes encoding resistance to the antibiotic, tetracycline (Tc), is repressed by tetracycline repressor (TetR), which is a homodimeric alpha-helical protein possessing a small N-terminal DNA binding domain (DNB domain) and a larger C-terminal domain (TBD domain). Binding of Tc to the TBD domain induces a subtle conformational change in the DNB domain that leads to abrogation of its DNA-binding activity, and induction of Tc resistance. While most previous studies on TetR have focused on the effects of Tc-binding on the DNB domain conformation, here we have investigated the role of the DNB domain in modulating Tc binding. We have discovered that a TBD domain construct entirely lacking the DNB domain displays a drastic reduction in Tc-binding affinity even though the DNB domain is far from the Tc-binding site. In the context of full-length TetR, highly destabilizing amino acid substitutions in the DNB domain cause reductions in Tc-binding activity. Strikingly, the DNB domains of these mutants, which are completely unfolded in the absence of Tc, are induced to fold when Tc is bound. These results demonstrate that there is a previously unrecognized two-way interdomain signaling mechanism in TetR whereby the DNB domain is required for maximal Tc-binding by the TBD domain, and Tc-binding in the TBD domain leads to stabilization of the DNB domain. Furthermore, our work suggests that detailed thermodynamic and kinetic studies on mutant forms of other allosteric proteins may also reveal surprising and previously undetected modes of interdomain communication.

Pharmacologic transgene control systems for gene therapy

Pharmacologic transgene-expression dosing is considered essential for future gene therapy scenarios. Genetic interventions require precise transcription or translation fine-tuning of therapeutic transgenes to enable their titration into the therapeutic window, to adapt them to daily changing dosing regimes of the patient, to integrate them seamlessly into the patient's transcriptome orchestra, and to terminate their expression after successful therapy. In recent years, decisive progress has been achieved in designing high-precision trigger-inducible mammalian transgene control modalities responsive to clinically licensed and inert heterologous molecules or to endogenous physiologic signals. Availability of a portfolio of compatible transcription control systems has enabled assembly of higher-order control circuitries providing simultaneous or independent control of several transgenes and the design of (semi-)synthetic gene networks, which emulate digital expression switches, regulatory transcription cascades, epigenetic expression imprinting, and cellular transcription memories. This review provides an overview of cutting-edge developments in transgene control systems, of the design of synthetic gene networks, and of the delivery of such systems for the prototype treatment of prominent human disease phenotypes.

Sequence correlations between Cro recognition helices and cognate O(R) consensus half-sites suggest conserved rules of protein-DNA recognition

The O(R) regions from several lambdoid bacteriophages contain the three regulatory sites O(R)1, O(R)2 and O(R)3, to which the Cro and CI proteins can bind. These sites show imperfect dyad symmetry, have similar sequences, and generally lie on the same face of the DNA double helix. We have developed a computational method, which analyzes the O(R) regions of additional phages and predicts the location of these three sites. After tuning the method to predict known O(R) sites accurately, we used it to predict unknown sites, and ultimately compiled a database of 32 known and predicted O(R) binding site sets. We then identified sequences of the recognition helices (RH) for the cognate Cro proteins through manual inspection of multiple sequence alignments. Comparison of Cro RH and consensus O(R) half-site sequences revealed strong one-to-one correlations between two amino acids at each of three RH positions and two bases at each of three half-site positions (H1-->2, H3-->5 and H6-->6). In each of these three cases, one of the two amino acid/base-pairings corresponds to a contact observed in the crystal structure of a lambda Cro/consensus operator complex. The alternate amino acid/base combinations were rationalized using structural models. We suggest that the pairs of amino acid residues act as binary switches that efficiently modulate specificity for different consensus half-site variants during evolution. The observation of structurally reasonable amino acid-to-base correlations suggests that Cro proteins share some common rules of recognition despite their functional and structural diversity.

Tet repressor residues indirectly recognizing anhydrotetracycline

Two tetracycline repressor (TetR) sequence variants sharing 63% identical amino acids were investigated in terms of their recognition specificity for tetracycline and anhydrotetracycline. Thermodynamic complex stabilities determined by urea-dependent unfolding reveal that tetracycline stabilizes both variants to a similar extent but that anhydrotetracycline discriminates between them significantly. Isofunctional TetR hybrid proteins of these sequence variants were constructed and their denaturation profiles identified residues 57 and 61 as the complex stability determinant. Association kinetics reveal different recognition of these TetR variants by anhydrotetracycline, but the binding constants indicate similar stabilization. The identified residues connect to an internal water network, which suggests that the discrepancy in the observed thermodynamics may be caused by an entropy effect. Exchange of these interacting residues between the two TetR variants appears to influence the flexibility of this water organization, demonstrating the importance of buried, structural water molecules for ligand recognition and protein function. Therefore, this structural module seems to be a key requisite for the plasticity of the multiple ligand binding protein TetR.

Folding and association of an extremely stable dimeric protein from Sulfolobus islandicus

ORF56 is a plasmid-encoded protein from Sulfolobus islandicus, which probably controls the copy number of the pRN1 plasmid by binding to its own promotor. The protein showed an extremely high stability in denaturant, heat, and pH-induced unfolding transitions, which can be well described by a two-state reaction between native dimers and unfolded monomers. The homodimeric character of native ORF56 was confirmed by analytical ultracentrifugation. Far-UV circular dichroism and fluorescence spectroscopy gave superimposable denaturant-induced unfolding transitions and the midpoints of both heat as well as denaturant-induced unfolding depend on the protein concentration supporting the two-state model. This model was confirmed by GdmSCN-induced unfolding monitored by heteronuclear 2D NMR spectroscopy. Chemical denaturation was accomplished by GdmCl and GdmSCN, revealing a Gibbs free energy of stabilization of -85.1 kJ/mol at 25 degrees C. Thermal unfolding was possible only above 1 M GdmCl, which shifted the melting temperature (t(m)) below the boiling point of water. Linear extrapolation of t(m) to 0 M GdmCl yielded a t(m) of 107.5 degrees C (5 microM monomer concentration). Additionally, ORF56 remains natively structured over a remarkable pH range from pH 2 to pH 12. Folding kinetics were followed by far-UV CD and fluorescence after either stopped-flow or manual mixing. All kinetic traces showed only a single phase and the two probes revealed coincident folding rates (k(f), k(u)), indicating the absence of intermediates. Apparent first-order refolding rates depend linearly on the protein concentration, whereas the unfolding rates do not. Both lnk(f) and lnk(u) depend linearly on the GdmCl concentration. Together, folding and association of homodimeric ORF56 are concurrent events. In the absence of denaturant ORF56 refolds fast (7.0 x 10(7)M(-1)s(-1)) and unfolds extremely slowly (5.7 year(-1)). Therefore, high stability is coupled to a slow unfolding rate, which is often observed for proteins of extremophilic organisms.

Two mutations in the tetracycline repressor change the inducer anhydrotetracycline to a corepressor

We report for the first time the in vitro characterization of a reverse tetracycline repressor (revTetR). The dimeric wild-type repressor (TetR) binds to tet operator tetO in the absence of the inducer anhydrotetracycline (atc) to confer tight repression. We have isolated the revTetR G96E L205S mutant, which, contrary to TetR, binds tetO only in the presence of atc. This reverse acting mutant was overproduced and purified. Effector and DNA binding properties were analyzed by EMSA and quantified by fluorescence titration and surface plasmon resonance. The association constant K(A) of revTetR for binding of [atcMg](+) is approximately 10(8) M(-1), four orders of magnitude lower than that of TetR. The affinity of TetR for tetO is 5.6 +/- 2 x 10(9) M(-1) and that for revTetR in the presence of atc is 1 +/- 0.2 x 10(8) M(-1). Both induced forms, the atc-bound TetR and the free revTetR, have the same low affinity of 4 +/- 1 x 10(5) M(-1) for DNA. Therefore, atc does not act as a dimerization agent for revTetR. We discuss the structural differences between TetR and revTetR potentially underlying this reversal of activity.

A protein contortionist: core mutations of GB1 that induce dimerization and domain swapping

Immunoglobulin-binding domain B1 of streptococcal protein G (GB1), a small (56 residues), stable, single-domain protein, is one of the most extensively used model systems in the area of protein folding and design. Recently, NMR and X-ray structures of a quintuple GB1 core mutant (L5V/A26F/F30V/Y33F/A34F) that showed an unexpected, intertwined tetrameric architecture were determined. Here, we report the NMR structure of another mutant, derived from the tetramer by reverting the single amino acid position F26 back to the wild-type sequence A26. The structure reveals a domain-swapped dimer that involves exchange of the second beta-hairpin. The resulting overall structure comprises an eight-stranded beta-sheet whose concave side is covered by two alpha helices. The dimer dissociates into a partially folded, monomeric species with a dissociation constant of 93(+/-10)microM.

A differential scanning calorimetry study of tetracycline repressor

Tetracycline repressor (TetR), which constitutes the most common mechanism of bacterial resistance to an antibiotic, is a homodimeric protein composed of two identical subunits, each of which contains a domain possessing a helix-turn-helix motif and a domain responsible for binding tetracycline. Binding of tetracycline in the protein pocket is accompanied by conformational changes in TetR, which abolish the specific interaction between the protein and DNA. Differential scanning calorimetry (DSC) and CD measurements, performed at pH 8.0, were used to observe the thermal denaturation of TetR in the absence and presence of tetracycline. The DSC results show that, in the absence of tetracycline, the thermally induced transitions of TetR can be described as an irreversible process, strongly dependent on scan rate and indicating that the protein denaturation is under kinetic control described by the simple kinetic scheme: N(2)--->D(2), where k is a first-order kinetic constant, N is the native state, and D is the denatured state. On the other hand, analysis of the scan rate effect on the transitions of TetR in the presence of tetracycline shows that thermal unfolding of the protein can be described by the two-state model: N(2)<--->U(2)--->D. In the proposed model, TetR in the presence of tetracycline undergoes co-operative unfolding, characterized by an enthalpy change (DeltaH(cal) = 1067 kJ x mol(-1)) and an entropy change (DeltaS = 3.1 kJ x mol(-1)).

Equilibrium unfolding of dimeric human prostatic acid phosphatase involves an inactive monomeric intermediate

Guanidine hydrochloride (GdnHCl)-induced unfolding of human prostatic acid phosphatase (hPAP), a homodimer of 50 kDa subunit molecular weight, was investigated with activity measurements, size exclusion HPLC, tryptophan fluorescence, 1-anilinonaphtalene-8-sulfonate (ANS) binding and reactivity with 2-(4'-maleimidoanilino)naphthalene-6-sulfonate (MIANS). Equilibrium analysis was performed to shed light on the role of dimerization in the folding and stability of the catalytically active oligomeric protein. Unfolding was reversible, as verified by activity measurements and tryptophan fluorescence. The noncoincidence of the unfolding curves obtained by different techniques suggests the occurrence of a multiphasic process. The reaction of hPAP inactivation is accompanied by dissociation of the dimer into two monomers. The midpoint of this transition is at 0.65 M GdnHCl with 4.24+/-0.12 kcalmol(-1) free energy change. Binding of ANS to the inactive phosphatase monomer, especially remarkable in the region from 0.8 to 1.25M GdnHCl, suggests that the hydrophobic probe indicates exposition of the intersubunit hydrophobic surface and a loosening of the monomer's tertiary structure. Strong fluorescence of thiol group derivatives, the products of their reaction with MIANS, appears in a limited range of GdnHCl concentrations (1.2-1.6M). This shows that in the relaxed structure of the intermediate, the reagent is allowed to penetrate into the hydrophobic environment of the partially hidden thiol groups. The equilibrium unfolding reaction of hPAP, as monitored by tryptophan fluorescence, does not depend on the protein concentration and displays a single transition curve with a midpoint at 1.7 M GdnHCl and value of DeltaG(unf)(H(2)O)=3.38+/-0.08 kcalmol(-1) per monomer, a result implying that this transition is related to the conformational change of the earlier dissociated and already inactive subunit of the protein.

Characterization of interactions between the transcriptional repressor PhlF and its binding site at the phlA promoter in Pseudomonas fluorescens F113

The phlACBD genes responsible for the biosynthesis of the antifungal metabolite 2,4-diacetylphloroglucinol (PHL) by the biocontrol strain Pseudomonas fluorescens F113 are regulated at the transcriptional level by the pathway-specific repressor PhlF. Strong evidence suggests that this regulation occurs mainly in the early logarithmic phase of growth. First, the expression of the phlF gene is relatively high between 3 and 13 h of growth and relatively low thereafter, with the phlACBD operon following an opposite expression profile. Second, the kinetics of PHL biosynthesis are specifically altered in the logarithmic phase in a P. fluorescens F113 phlF mutant. The phlA-phlF intergenic region presents a complex organization in that phlACBD is transcribed from a sigma(70) RNA polymerase-dependent promoter that is likely to overlap the promoter of the divergently transcribed phlF gene. The repression by PhlF is due to its interaction with an inverted repeated sequence, phO, located downstream of the phlA transcriptional start site. Cross-linking experiments indicate that PhlF can dimerize in solution, and thus PhlF may bind phO as a dimer or higher-order complex. Furthermore, it is now demonstrated that certain regulators of PHL synthesis act by modulating PhlF binding to phO. PHL, which has previously been shown to be an autoinducer of PHL biosynthesis, interacts with PhlF to destabilize the PhlF-phO complex. Conversely, the PhlF-phO complex is stabilized by the presence of salicylate, which has been shown to be an inhibitor of phlA expression.

Mechanism of Tet repressor induction by tetracyclines: length compensates for sequence in the alpha8-alpha9 loop

Natural Tet repressor (TetR) variants are alpha-helical proteins bearing a large loop between helices 8 and 9, which is variable in sequence and length. We have deleted this loop consisting of 14 amino acid residues in TetR(D) and rebuilt it stepwise with up to 42 alanine residues. All except the mutant with the longest alanine loop show wild-type repression, but none is inducible with tetracycline. This demonstrates the importance of the alpha8-alpha9 loop and its amino acid sequence for induction. The induction efficiencies increase with loop length, when the more tightly binding inducer anhydrotetracycline is used. The largest increase of inducibility was observed for TetR mutants with loop lengths between eight and 17 alanine residues. Since loop residues Asp/Glu157 and Arg158 are conserved in the natural TetR sequence variants, we constructed a mutant in which all other residues of the loop were replaced by alanine. This mutant exhibits increased anhydrotetracycline induction compared to the corresponding alanine variant. Thus, these residues are important for induction. Binding constants for the anhydrotetracycline-TetR interaction are below the detection level of 10(5) M(-1) for the mutant with a loop of two alanine residues and increase sharply until a loop size of ten residues is reached. TetR variants with longer loops have similar anhydrotetracycline-binding constants, ranging between 2.6 x 10(9) M(-1) and 8.0 x 10(9) M(-1), about 500-fold lower than wild-type TetR. The increase of the affinity occurs at shorter loop lengths than that of inducibility. We conclude that the induction defect of the polyalanine variants arises from two increments: (i) the loop must have a minimal length-to allow efficient inducer binding; (ii) the loop must structurally participate in the conformational change associated with induction.

The thermodynamic stability of the proteins of the ccd plasmid addiction system

The two opponents, toxin (CcdB, LetB or LetD, protein G, LynB) and antidote (CcdA, LetA, protein H, LynA), in the plasmid addiction system ccd of the F plasmid were studied by different biophysical methods. The thermodynamic stability was measured at different temperatures combining denaturant and thermally induced unfolding. It was found that both proteins denature in a two-state equilibrium (native dimer versus unfolded monomer) and that CcdA has a significantly lower thermodynamic stability. Using a numerical model, which was developed earlier by us, and on the basis of the determined thermodynamic parameters the concentration dependence of the denaturation transition temperature was obtained for both proteins. This concentration dependence may be of physiological significance, as the concentration of both ccd addiction proteins cannot exceed a certain limit because their expression is controlled by autoregulation. The influence of DNA on the thermal stability of the two proteins was probed. It was found that cognate DNA increases the melting temperature of CcdA. In the presence of non-specific DNA the thermal stability was not changed. The melting temperature of CcdB was not influenced by the applied double-stranded oligonucleotides, neither cognate nor unspecific.

Engineering stability in gene networks by autoregulation

The genetic and biochemical networks which underlie such things as homeostasis in metabolism and the developmental programs of living cells, must withstand considerable variations and random perturbations of biochemical parameters. These occur as transient changes in, for example, transcription, translation, and RNA and protein degradation. The intensity and duration of these perturbations differ between cells in a population. The unique state of cells, and thus the diversity in a population, is owing to the different environmental stimuli the individual cells experience and the inherent stochastic nature of biochemical processes (for example, refs 5 and 6). It has been proposed, but not demonstrated, that autoregulatory, negative feedback loops in gene circuits provide stability, thereby limiting the range over which the concentrations of network components fluctuate. Here we have designed and constructed simple gene circuits consisting of a regulator and transcriptional repressor modules in Escherichia coli and we show the gain of stability produced by negative feedback.

The ionization of a buried glutamic acid is thermodynamically linked to the stability of Leishmania mexicana triose phosphate isomerase

The amino acid sequence of Leishmania mexicana triose phosphate isomerase is unique in having at position 65 a glutamic acid instead of a glutamine. The stability properties of LmTIM and the E65Q mutant were investigated by pH and guanidinium chloride-induced unfolding. The crystal structure of E65Q was determined. Three important observations were made: (a) there are no structural rearrangements as the result of the substitution; (b) the mutant is more stable than the wild-type; and (c) the stability of the wild-type enzyme shows strong pH dependence, which can be attributed to the ionization of Glu65. Burying of the Glu65 side chain in the uncharged environment of the dimer interface results in a shift in pKa of more than 3 units. The pH-dependent decrease in overall stability is due to weakening of the monomer-monomer interactions (in the dimer). The E65Q substitution causes an increase in stability as the result of the formation of an additional hydrogen bond in each subunit (DeltaDeltaG degrees of 2 kcal.mol-1 per monomer) and the elimination of a charged group in the dimer interface (DeltaDeltaG degrees of at least 9 kcal.mol-1 per dimer). The computated shift in pKa and the stability of the dimer calculated from the charge distribution in the protein structure agree closely with the experimental results. The guanidinium chloride dependence of the unfolding constant was smaller than expected from studies involving monomeric model proteins. No intermediates could be identified in the unfolding equilibrium by combining fluorescence and CD measurements. Study of a stable monomeric triose phosphate isomerase variant confirmed that the phenomenon persists in the monomer.

Solvent-exposed residues in the Tet repressor (TetR) four-helix bundle contribute to subunit recognition and dimer stability

Dimerization specificity of Tet repressor (TetR) can be altered by changes in the core of the four-helix bundle that mediates protein-protein recognition. We demonstrate here that the affinity of subunit interaction depends also on the solvent-exposed residues at positions 128 and 179'-184', which interact across the dimerization surface. TetR(B) and (D), two naturally occurring sequence variants, differ at position 128 with respect to the monomer-monomer distances in the crystal structures and the charge of the amino acids, being glutamate in TetR(B) and arginine in TetR(D). In vivo analysis of chimeric TetR(B/D) variants revealed that the single E128R exchange does not alter the dimerization specificity of TetR(B) to the one of TetR(D). When combined with specificity mutations in alpha10, it is, however, able to increase dimerization efficiency of the TetR(B/D) chimera with TetR(D). A loss of contact analysis revealed a positive interaction between Arg-128 and residues located at positions 179'-184' of the second monomer. We constructed a hyperstable TetR(B) variant by replacing residues 128 and 179-184 by the respective TetR(D) sequence. These results establish that in addition to a region in the hydrophobic core residues at the solvent-exposed periphery of the dimerization surface participate in protein-protein recognition in the TetR four-helix bundle.

Thermodynamics and kinetics of unfolding of the thermostable trimeric adenylate kinase from the archaeon Sulfolobus acidocaldarius

The thermal stability of adenylate kinase from the thermoacidophilic archaeon Sulfolobus acidocaldarius was characterized comprehensively using denaturant-induced unfolding, differential scanning calorimetry, circular dichroism spectroscopy, and enzymological inactivation studies. The thermally induced unfolding of the protein is irreversible due to aggregation, whereas the unfolding induced by guanidinium chloride is reversible. The protein is known to be a homotrimer in its native state and we established that it unfolds upon dissociation in the case of denaturant unfolding. We measured the thermodynamic stability of the protein in a temperature range from 5 to 70 degrees C using denaturant unfolding. The protein has a maximum of stability (intrinsic free energy) of 31 kcal/mol-trimer (130 kJ/mol-trimer) at 32 degrees C (based on the linear extrapolation model). The heat capacity change upon unfolding DeltaCp and the m-value were considered to be constant in this temperature range and calculated to be 2.86 kcal/mol-trimer (11.9 kJ/mol-trimer) and 5.67 kcal/mol-trimer M (23.7 kJ/mol-trimer M), respectively. The influence of trimerization on thermodynamic stability was investigated. The several interrelated aspects of thermal stability such as unfolding kinetics, the temperature-dependence of the free energy, and the concentration and temperature-dependencies of the fraction of denatured protein are described quantitatively. The properties of the Gibbs-Helmholtz function of the adenylate kinase from S. acidocaldarius, in particular, and of oligomeric proteins, in general terms, are discussed and compared with the properties of the analogous function for monomeric proteins. Moreover, we discuss methodological aspects: we obtained the analytical expression of the denaturant-unfolding isotherm for homotrimeric proteins; we include a formula Appendix containing the derivations of the expressions used.

Stepwise selection of TetR variants recognizing tet operator 4C with high affinity and specificity

The TetR PQ39 mutant exhibits a new recognition specificity for the tetO-4C operator, but the affinity is not sufficiently high for use in vivo. A stepwise selection of additional mutations by cassette mutagenesis with randomization of residues in the TetR alpha-helix-turn-alpha-helix motif (HTH) yielded mutant TetR EA37PQ39YM42 showing a similar affinity and increased specificity for tetO-4C as wild-type TetR for tetO. A set of mutants obtained by that approach revealed that the fourth residue of the HTH (Leu41), which points towards the core of the DNA binding domain in TetR, alters the recognition of base-pair 4, e.g. the mutant TetR LV41YM42 exhibits a new recognition specificity for tetO-4G. A small residue at the last position in the turn of the HTH increases the affinity and specificity of DNA binding of TetR mutants containing the PQ39 exchange. Thus, cooperation between residues at positions 37, 39, 41 and 42 in the HTH of TetR is necessary to optimize recognition of base-pair 4. We conclude that creating a new DNA recognition specificity in the HTH of TetR with high affinity for the tetO-4C operator variant requires exchanges altering flexibility and/or adjustment of the recognition alpha-helix to the target DNA in addition to the contacting residue.

Determinants of protein-protein recognition by four helix bundles: changing the dimerization specificity of Tet repressor

Homo- and heterodimerization is essential for the activity of many proteins, particularly transcription factors. One widely distributed structural motif for protein recognition is the four helix bundle. To understand the molecular details determining specificity of subunit recognition in a dimer formed by a four helix bundle, we investigated Tet repressor (TetR) sequence variants TetR(B) and TetR(D), which do not form heterodimers. We used molecular modeling to identify residues with the potential to determine recognition of subunits. Directed mutagenesis of these residues in TetR(B) by the TetR(D) sequence resulted in chimeric TetR(B/D) repressors with new subunit recognition specificities. The single LS192 exchange in TetR(B/D)192 in the center of the helix bundle leads to a relaxed specificity since this variant dimerizes with TetR(B) and (D). To construct a variant with a new specificity it was not sufficient to mutate the contacting residue, F197, in the other subunit. Instead, it was necessary to exchange two more residues in the vicinity of F197 and S192. The resulting TetR(B/D)188, 192,193,197 forms dimers with TetR(D) but not with TetR(B), indicating that four amino acid exchanges are sufficient to change subunit recognition. These results establish that targeted alterations in the structural complementarity of protein-protein interaction surfaces can be used to construct new recognition specificities. However, it is not sufficient to adjust the complementary residues since the surrounding amino acids contribute essentially to protein-protein recognition.