Relationship between sequence determinants of stability for two natural homologous proteins with different folds

In the Cro protein family, an evolutionary change in secondary structure has converted an alpha-helical fold to a mixture of alpha-helix and beta-sheet. P22 Cro and lambda Cro represent the ancestral all-alpha and descendant alpha+beta folds, respectively. The major structural differences between these proteins are at the C-terminal end of the domain (residues 34-56), where two alpha-helices in P22 Cro align with two beta-strands in lambda Cro. We sought to assess the possibility that smooth evolutionary transitions could have converted the all-alpha structure to the alpha+beta structure through sequences that could adopt both folds. First, we used scanning mutagenesis to identify and compare patterns of key stabilizing residues in the C-terminal regions of both P22 Cro and lambda Cro. These patterns exhibited little similarity to each other, with structurally important residues in the two proteins most often occurring at different sequence positions. Second, "hybrid scanning" studies, involving replacement of each wild-type residue in P22 Cro with the aligned wild-type residue in lambda Cro and vice versa, revealed five or six residues in each protein that strongly destabilized the other. These results suggest that key stability determinants for each Cro fold are quite different and that the P22 Cro sequence strongly favors the all-alpha structure while the lambda Cro sequence strongly favors the alpha+beta structure. Nonetheless, we were able to design a "structurally ambivalent" sequence fragment (SASF1), which corresponded to residues 39-56 and simultaneously incorporated most key stabilizing residues for both P22 Cro and lambda Cro. NMR experiments showed SASF1 to stably fold as a beta-hairpin when incorporated into the lambda Cro sequence but as a pair of alpha-helices when incorporated into P22 Cro.     

High Resolution Structural Studies of Cro Repressor Protein and Implications for DNA Recognition

Abstract

Cro repressor is a small dimeric protein that binds to specific sites on the DNA of bacteriophage λ. The structure of Cro has been determined and suggests that the protein binds to its sequence-specific sites with a pair of two-fold related α-helices of the protein located within successive major grooves of the DNA.

From the known three-dimensional structure of the repressor, model building and energy refinement have been used to develop a detailed model for the presumed complex between Cro and DNA. Recognition of specific DNA binding sites appears to occur via multiple hydrogen bonds between amino acid side chains of the protein and base pair atoms exposed within the major groove of DNA. The Cro:DNA model is consistent with the calculated electrostatic potential energy surface of the protein.

From a series of amino acid sequence and gene sequence comparisons, it appears that a number of other DNA-binding proteins have an α-helical DNA-binding region similar to that seen in Cro. The apparent sequence homology includes not only DNA-binding proteins from different bacteriophages, but also gene-regulatory proteins from bacteria and yeast. It has also been found that the conformations of part of the presumed DNA-binding regions of Cro repressor, λ repressor and CAP gene activator proteins are strikingly similar. Taken together, these results strongly suggest that a two-helical structural unit occurs in the DNA-binding region of many proteins that regulate gene expression. However, the results to date do not suggest that there is a simple one-to-one recognition code between amino acids and bases.

Crystals have been obtained of complexes of Cro with six-base-pair and nine-basepair DNA oligomers, and X-ray analysis of these co-crystals is in progress.     

Design of lambda Cro fold: solution structure of a monomeric variant of the de novo protein

One of the classical DNA-binding proteins, bacteriophage lambda Cro, forms a homodimer with a unique fold of alpha-helices and beta-sheets. We have computationally designed an artificial sequence of 60 amino acid residues to stabilize the backbone tertiary structure of the lambda Cro dimer by simulated annealing using knowledge-based structure-sequence compatibility functions. The designed amino acid sequence has 25% identity with that of natural lambda Cro and preserves Phe58, which is important for formation of the stably folded structure of lambda Cro. The designed dimer protein and its monomeric variant, which was redesigned by the insertion of a beta-hairpin sequence at the C-terminal region to prevent dimerization, were synthesized and biochemically characterized to be well folded. The designed protein was monomeric under a wide range of protein concentrations and its solution structure was determined by NMR spectroscopy. The solved structure is similar to that of a monomeric variant of natural lambda Cro with a root-mean-square deviation of the polypeptide backbones at 2.1A and has a well-packed protein core. Thus, our knowledge-based functions provide approximate but essential relationships between amino acid sequences and protein structures, and are useful for finding novel sequences that are foldable into a given target structure.     

Functional expression and affinity selection of single-chain cro by phage display: isolation of novel DNA-binding proteins

A robust selection system affording phage display of the DNA-binding helix-turn-helix protein Cro is presented. The aim of the work was to construct an experimental system allowing for the construction and isolation of Cro-derived protein with new DNA-binding properties. A derivative of the phage lambda Cro repressor, scCro8, in which the protein subunits had been covalently connected via a peptide linker was expressed in fusion with the gene 3 protein of Escherichia coli filamentous phage. The phage-displayed single-chain Cro was shown to retain the DNA binding properties of its wild-type Cro counterpart regarding DNA sequence specificity and binding affinity. A kinetic analysis revealed the rate constant of dissociation of the single-chain Cro-phage/DNA complex to be indistinguishable from that of the free single-chain Cro. Affinity selection using a biotinylated DNA with a target consensus operator sequence allowed for a 3000-fold enrichment of phages displaying single-chain Cro over control phages. The selection was based on entrapment of phage/DNA complexes formed in solution on streptavidin-coated paramagnetic beads. The expression system was subsequently used to isolate variant scCro8 proteins, mutated in their DNA-binding residues, that specifically recognized new, unnatural target DNA ligands.     

N15 Cro and lambda Cro: orthologous DNA-binding domains with completely different but equally effective homodimer interfaces

Bacteriophage Cro proteins bind to target DNA as dimers but do not all dimerize with equal strength, and differ in fold in the region of the dimer interface. We report the structure of the Cro protein from Enterobacteria phage N15 at 1.05 A resolution. The subunit fold contains five alpha-helices and is closely similar to the structure of P22 Cro (1.3 A backbone room mean square difference over 52 residues), but quite different from that of lambda Cro, a structurally diverged member of this family with a mixed alpha-helix/beta-sheet fold. N15 Cro crystallizes as a biological dimer with an extensive interface (1303 A(2) change in accessible surface area per dimer) and also dimerizes in solution with a K(d) of 5.1 +/- 1.5 microM. Its dimerization is much stronger than that of its structural homolog P22 Cro, which does not self-associate detectably in solution. Instead, the level of self-association and interfacial area for N15 Cro is similar to that of lambda Cro, even though these two orthologs do not share the same fold and have dimer interfaces that are qualitatively different in structure. The common Cro ancestor is thought to be an all-helical monomer similar to P22 Cro. We propose that two Cro descendants independently developed stronger dimerization by entirely different mechanisms.     

Complementation of the lytD1 mutation of Escherichia coli by either the cI or cro gene of bacteriophage lambda.

The lytD1 mutant of Escherichia coli exhibits temperature-sensitive growth which is attributed to cellular autolysis at the restrictive temperature. Either of two cloned phage lambda genes, identified as cI and cro, suppressed the lytD1(Ts) lysis phenotype, suggesting that lytD encodes a DNA-binding protein with a DNA-binding specificity similar to that of CI and Cro. LytD may be a repressor of a gene(s) involved in cellular autolysis.     

Secondary structure switching in Cro protein evolution

We report the solution structure of the Cro protein from bacteriophage P22. Comparisons of its sequence and structure to those of lambda Cro strongly suggest an alpha-to-beta secondary structure switching event during Cro evolution. The folds of P22 Cro and lambda Cro share a three alpha helix fragment comprising the N-terminal half of the domain. However, P22 Cro's C terminus folds as two helices, while lambda Cro's folds as a beta hairpin. The all-alpha fold found for P22 Cro appears to be ancestral, since it also occurs in cI proteins, which are anciently duplicated paralogues of Cro. PSI-BLAST and transitive homology analyses strongly suggest that the sequences of P22 Cro and lambda Cro are globally homologous despite encoding different folds. The alpha+beta fold of lambda Cro therefore likely evolved from its all-alpha ancestor by homologous secondary structure switching, rather than by nonhomologous replacement of both sequence and structure.     

Terminase host factor: a histone-like E. coli protein which can bind to the cos region of bacteriophage lambda DNA.

Terminase Host Factor (THF), an E. coli protein capable of fulfilling the host factor requirement for in vitro bacteriophage lambda terminase activity, displays properties characteristic of the prokaryotic type II DNA-binding or "histone-like" proteins. It is a 22 K basic, heat- and acid-stable protein which binds non-specifically to various DNAs. Conditions can be established, however, where THF binds preferentially to the cohesive end site (cos) of lambda DNA forming several distinct complexes as visualized by band retardation in polyacrylamide gels. DNase I footprinting reveals that THF can protect several regions of the top strand on the right side (+) of cos but does not bind as well to the left side (-). The binding regions are separated either by unprotected or by DNase I- hypersensitive bases. Under the conditions used in these experiments, DNA which does not contain cos lambda sequences does not show this pattern of protection. Several repeated motifs in the cos lambda nucleotide sequence may represent a consensus sequence for THF interaction. THF may be similar to other "histone-like" proteins which display both non-specific and selective DNA-binding capacities.     

Folding kinetics of phage 434 Cro protein

Folding kinetics for phage 434 Cro protein are examined and compared with those reported for lambda(6-85), the N-terminal domain of the repressor of phage lambda. The two proteins have similar all-helical structures consisting of five helices but different stabilities. In contrast to lambda(6-85), sharp and distinct aromatic (1)H NMR signals without exchange broadening characterize the native and urea-denatured 434 Cro forms at equilibrium at 20 degrees C, indicating slow interconversion on the NMR time scale. Stopped-flow fluorescence data using the single 434 Cro tryptophan indicate strongly urea-dependent refolding rates and smaller urea dependencies of the unfolding rates, suggesting a native-like transition state ensemble. Refolding rates are slower and unfolding rates considerably faster at pH 4 than at pH 6. This accounts for the lower stability of 434 Cro at pH 4 and suggests the existence of pH-dependent, possibly salt bridge interactions that are more stabilizing at pH 6. At <2 M urea, decreased folding amplitudes and nonlinear urea dependencies that are apparent at pH 6 indicate deviation from two-state behavior and suggest the formation of an early folding intermediate. The folding behavior of 434 Cro and why it folds 2 orders of magnitude slower than lambda(6-85) are rationalized in terms of the lower intrinsic helix stabilities and putative charge interactions in 434 Cro.     

Conservation and diversity in the immunity regions of wild phages with the immunity specificity of phage lambda

The gene regulatory circuitry of phage lambda is among the best-understood circuits. Much of the circuitry centres around the immunity region, which includes genes for two repressors, CI and Cro, and their cis-acting sites. Related phages, termed lambdoid phages, have different immunity regions, but similar regulatory circuitry and genome organization to that of lambda, and show a mosaic organization, arising by recombination between lambdoid phages. We sequenced the immunity regions of several wild phages with the immunity specificity of lambda, both to determine whether natural variation exists in regulation, and to analyse conservation and variability in a region rich in well-studied regulatory elements. CI, Cro and their cis-acting sites are almost identical to those in lambda, implying that regulatory mechanisms controlled by the immunity region are conserved. A segment adjacent to one of the operator regions is also conserved, and may be a novel regulatory element. In most isolates, different alleles of two regulatory proteins (N and CII) flank the immunity region; possibly the lysis-lysogeny decision is more variable among isolates. Extensive mosaicism was observed for several elements flanking the immunity region. Very short sequence elements or microhomologies were also identified. Our findings suggest mechanisms by which fine-scale mosaicism arises.     

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.     

Folding propensities of synthetic peptide fragments covering the entire sequence of phage 434 Cro protein.

The phage 434 Cro protein, the N-terminal domain of its repressor (R1-69) and that of phage lambda (lambda6-85) constitute a group of small, monomeric, single-domain folding units consisting of five helices with striking structural similarity. The intrinsic helix stabilities in lambda6-85 have been correlated to its rapid folding behavior, and a residual hydrophobic cluster found in R1-69 in 7 M urea has been proposed as a folding initiation site. To understand the early events in the folding of 434 Cro, and for comparison with R1-69 and lambda6-85, we examined the conformational behavior of five peptides covering the entire 434 Cro sequence in water, 40% (by volume) TFE/water, and 7 M urea solutions using CD and NMR. Each peptide corresponds to a helix and adjacent residues as identified in the native 434 Cro NMR and crystal structures. All are soluble and monomeric in the solution conditions examined except for the peptide corresponding to the 434 Cro helix 4, which has low water solubility. Helix formation is observed for the 434 Cro helix 1 and helix 2 peptides in water, for all the peptides in 40% TFE and for none in 7 M urea. NMR data indicate that the helix limits in the peptides are similar to those in the native protein helices. The number of side-chain NOEs in water and TFE correlates with the helix content, and essentially none are observed in 7 M urea for any peptide, except that for helix 5, where a hydrophobic cluster may be present. The low intrinsic folding propensities of the five helices could account for the observed stability and folding behavior of 434 Cro and is, at least qualitatively, in accord with the results of the recently described diffusion-collision model incorporating intrinsic helix propensities.     

Intersubunit disulfide-bonded lambda-Cro protein

Site-directed mutagenesis has been employed to substitute cysteine for valine at position 55, which is located on the dimer interface of the Cro protein of bacteriophage lambda. It has been found that the Cys55 Cro protein (Cro VC55) spontaneously forms a stable disulfide-bonded dimer in the absence of a reducing agent. UV-CD and NMR data showed that the mutant protein retains the conformation of the wild Cro protein and has acquired significant heat-stability. However, its specific DNA-binding activity is reduced several times compared with that of the wild Cro. Photochemically induced dynamic nuclear polarization (CIDNP) spectra demonstrated that a conformational change of Cro VC55 did not take place upon the formation of a complex with OR3, in contrast to the case of the wild Cro. These data suggest that the induced fitting, like loosening, of the two subunits of the wild Cro dimer contributes to the enhancement of its affinity to its operator DNA, which results in a specific interaction between Cro and OR3.     

Many gene-regulatory proteins appear to have a similar α-helical fold that binds DNA and evolved from a common precursor

Summary Amino acid and DNA sequence comparisons suggest that many sequence-specific DNA-binding proteins have in common and homologous region of about 22 amino acids. This region corresponds to two consecutive α-helices that occur in bot Cro and cI repressor proteins of bacteriophage λ and in catabolite gene activator protein ofEscherichia coli and are presumed to interact with DNA. The results obtained here suggest that this α-helical DNA-binding fold occurs in many proteins that regulate gene expression. It also appears that this DNA-binding unit evolved from a common evolutionary precursor.     

Phage λ Cro Protein and Ci Repressor Use Two Different Patterns of Specific Protein-DNA Interactions to Achieve Sequence Specificity in Vivo

By assaying the binding of wild-type Cro to a set of 40 mutant lambda operators in vivo, we have determined that the 14 outermost base pairs of the 17 base pair, consensus lambda operator are critical for Cro binding. Cro protein recognizes 4 base pairs in a lambda operator half-site in different ways than cI repressor. The sequence determinants of Cro binding at these critical positions in vivo are nearly perfectly consistent with the model proposed by W. F. ANDERSON, D. H. OHLENDORF, Y. TAKEDA and B. W. MATTHEWS and modified by Y. TAKEDA, A. SARAI and V. M. RIVERA for the specific interactions between Cro and its operator, and explain the relative order of affinities of the six natural lambda operators for Cro. Our data call into question the idea that lambda repressor and Cro protein recognize the consensus lambda operator by nearly identical patterns of specific interactions.     

Histones H3 and H2a are homologous to the lambda repressor and cro proteins in 22 residue segments implicated in DNA binding

The histones H3 and H2a from calf thymus are homologous to the repressor and cro repressor proteins of bacteriophage lambda in a 22-residue segment that has been implicated by mutational and model-building studies in DNA binding. In the lambda proteins this segment is folded into a helix-turn-helix unit of supersecondary structure, and we propose that the homologous regions in the histones possess the same fold. Homology was quantified with a unified procedure based on criteria of identity of key residues, primary structural homology and similarity of secondary structural potential. It has previously been shown that a set of other prokaryotic DNA-binding proteins have primary structural homology with the two lambda proteins. Homologies detected between the histones H4 and H2b and members of this set suggest that these histones also contain the putative DNA-binding fold.