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.     

A trade between similar but nonequivalent intrasubunit and intersubunit contacts in Cro dimer evolution

The homodimeric lambda Cro protein has a "ball-and-socket" interface that includes insertion of an aromatic side chain, Phe 58, from each subunit into a cavity in the hydrophobic core of the other subunit. This overlap between the subunit core and dimer interface hypothetically explains the strong dimerization and weak monomer stability of lambda Cro in comparison to homologues. According to a model developed here and in a previous study [LeFevre, K. R., and Cordes, M. H. (2003) Proc. Natl. Acad. Sci. U.S.A. 100, 2345-2350], the socket cavity evolved in part by replacement of a buried tryptophan in an ancestral stable monomer with a smaller side chain (Ala 33 in lambda Cro). The resulting core defect was in effect repaired by insertion of a different side chain (Phe 58) from a second subunit, generating the ball and socket. Consistent with such an evolutionary trade between intrasubunit and intersubunit interactions, we showed in the previous study that restoration of the ancestral Trp 33 in lambda Cro stabilized the monomer and reduced the extent of dimerization. Here, we report the solution structure of a stable lambda Cro monomer containing the Ala33Trp mutation, which confirms that the restored tryptophan fulfills its ancestral role as a core side chain, filling part of the socket cavity occupied by Phe 58 in the wild-type dimer. The structure also reveals, however, that the cavity is not completely filled by Trp 33, suggesting that its formation could have involved multiple mutations that reduced side chain volume. We offer suggestive evidence of a role of mutations at a second position.     

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.     

Domains in lambda Cro repressor. A calorimetric study.

Thermodynamic properties of a mutant lambda Cro repressor with Cys replacing Val55 were studied calorimetrically. Formation of the S-S cross-link between neighboring Cys55 residues in this dimeric molecule leads to stabilization of a structure formed by the C-terminal parts of the two polypeptide chains, which behave as a single cooperative domain upon protein denaturation by heating. This composite domain is very stable at neutral pH and disrupts at 110 degrees C. The S-S-cross-linked tryptic fragment (residues 22-66), which includes this C-terminal domain, has similar stability. The N-terminal parts of the polypeptide chains do not form any stable structure when isolated, but in S-S-cross-linked dimer, they form a single cooperative block which melts in an all-or-none way 9 degrees C higher than the un-cross-linked protein. The observed cooperation of the distant N-terminal parts in dimer raises questions regarding lambda Cro repressor structure in solution.     

Solution structure and dynamics of a designed monomeric variant of the lambda Cro repressor.

The solution structure of a monomeric variant of the lambda Cro repressor has been determined by multidimensional NMR. Cro K56[DGEVK] differs from wild-type Cro by the insertion of five amino acids at the center of the dimer interface. 1H and 15N resonances for 70 of the 71 residues have been assigned. Thirty-two structures were calculated by hybrid distance geometry/simulated annealing methods using 463 NOE-distance restraints, 26 hydrogen-bond, and 39 dihedral-angle restraints. The root-mean-square deviation (RMSD) from the average structure for atoms in residues 3-60 is 1.03 +/- 0.44 A for the peptide backbone and 1.6 +/- 0.73 A for all nonhydrogen atoms. The overall structure conforms very well to the original design. Although the five inserted residues form a beta hairpin as expected, this engineered turn as well as other turns in the structure are not well defined by the NMR data. Dynamics studies of backbone amides reveal T1/T2 ratios of residues in the alpha2-alpha3, beta2-beta3, and engineered turn that are reflective of chemical exchange or internal motion. The solution structure and dynamics are discussed in light of the conformational variation that has been observed in other Cro structures, and the importance of flexibility in DNA recognition.     

Folding and assembly of lambda Cro repressor dimers are kinetically limited by proline isomerization

Cro binds to operator sites in lambda DNA as a dimer. Dimerization of this small repressor protein is weak, however, and proline residues in the dimer interface suggest that folding and assembly of active repressors may be complex. Cro and selected variants have been studied by circular dichroism and fluorescence. Fluorescent probes include a unique tryptophan residue in the dimer interface and extrinsic resonance energy transfer probes that monitor dimerization. Both folding and unfolding are characterized by two distinct kinetic phases. Fast processes that are complete within the 5-10 ms dead time of stopped flow experiments account for the majority of the change in the CD signal and abrupt changes in both tryptophan fluorescence and energy transfer. The slow phases show all the hallmarks of proline isomerization. The rates of the slow phases are between 0.005 and 0.02 s(-1), are relatively independent of protein and denaturant concentration, display activation energies of 20 kcal/mol, and are accelerated by the peptidyl-prolyl isomerase SlyD. Although CD measurements indicate that more than 70% of the secondary structure is regained in the refolding burst phase, intermolecular fluorescence resonance energy transfer experiments indicate that less than 25% of these subunits are assembled into dimers. Full folding and dimerization requires isomerization of the non-native prolyl isomers over hundreds of seconds.     

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.     

Packing interface energetics in different crystal forms of the λ Cro dimer

Variation among crystal structures of the λ Cro dimer highlights conformational flexibility. The structures range from a wild type closed to a mutant fully open conformation, but it is unclear if each represents a stable solution state or if one may be the result of crystal packing. Here we use molecular dynamics (MD) simulation to investigate the energetics of crystal packing interfaces and the influence of site‐directed mutagenesis on them in order to examine the effect of crystal packing on wild type and mutant Cro dimer conformation. Replica exchange MD of mutant Cro in solution shows that the observed conformational differences between the wild type and mutant protein are not the direct consequence of mutation. Instead, simulation of Cro in different crystal environments reveals that mutation affects the stability of crystal forms. Molecular Mechanics Poisson‐Boltzmann Surface Area binding energy calculations reveal the detailed energetics of packing interfaces. Packing interfaces can have diverse properties in strength, energetic components, and some are stronger than the biological dimer interface. Further analysis shows that mutation can strengthen packing interfaces by as much as ～5 kcal/mol in either crystal environment. Thus, in the case of Cro, mutation provides an additional energetic contribution during crystal formation that may stabilize a fully open higher energy state. Moreover, the effect of mutation in the lattice can extend to packing interfaces not involving mutation sites. Our results provide insight into possible models for the effect of crystallization on Cro conformational dynamics and emphasize careful consideration of protein crystal structures. Proteins 2014; 82:1128–1141. © 2013 Wiley Periodicals, Inc.     

The structural basis for enhanced stability and reduced DNA binding seen in engineered second-generation Cro monomers and dimers

It was previously shown that the Cro repressor from phage lambda, which is a dimer, can be converted into a stable monomer by a five-amino acid insertion. Phe58 is the key residue involved in this transition, switching from interactions which stabilize the dimer to those which stabilize the monomer. Structural studies, however, suggested that Phe58 did not penetrate into the core of the monomer as well as it did into the native dimer. This was strongly supported by the finding that certain core-repacking mutations, including in particular, Phe58-->Trp, increased the stability of the monomer. Unexpectedly, the same substitution also increased the stability of the native dimer. At the same time it decreased the affinity of the dimer for operator DNA. Here we describe the crystal structures of the Cro F58W mutant, both as the monomer and as the dimer. The F58W monomer crystallized in a form different from that of the original monomer. In contrast to that structure, which resembled the DNA-bound form of Cro, the F58W monomer is closer in structure to wild-type (i.e. non-bound) Cro. The F58W dimer also crystallizes in a form different from the native dimer but has a remarkably similar overall structure which tends to confirm the large changes in conformation of Cro on binding DNA. Introduction of Trp58 perturbs the position occupied by the side-chain of Arg38, a DNA-contact residue, providing a structural explanation for the reduction in DNA-binding affinity.The improved thermal stability is seen to be due to the enhanced solvent transfer free energy of Trp58 relative to Phe58, supplemented in the dimer structure, although not the monomer, by a reduction in volume of internal cavities.     

Two structures of a lambda Cro variant highlight dimer flexibility but disfavor major dimer distortions upon specific binding of cognate DNA

Previously reported crystal structures of free and DNA-bound dimers of λ Cro differ strongly (about 4 Å backbone rmsd), suggesting both flexibility of the dimer interface and induced-fit protein structure changes caused by sequence-specific DNA binding. Here, we present two crystal structures, in space groups P3₂21 and C2 at 1.35 and 1.40 Å resolution, respectively, of a variant of λ Cro with three mutations in its recognition helix (Q27P/A29S/K32Q, or PSQ for short). One dimer structure (P3₂21; PSQ form 1) resembles the DNA-bound wild-type Cro dimer (1.0 Å backbone rmsd), while the other (C2; PSQ form 2) resembles neither unbound (3.6 Å) nor bound (2.4 Å) wild-type Cro. Both PSQ form 2 and unbound wild-type dimer crystals have a similar interdimer β-sheet interaction between the β1 strands at the edges of the dimer. In the former, an infinite, open β-structure along one crystal axis results, while in the latter, a closed tetrameric barrel is formed. Neither the DNA-bound wild-type structure nor PSQ form 1 contains these interdimer interactions. We propose that β-sheet superstructures resulting from crystal contact interactions distort Cro dimers from their preferred solution conformation, which actually resembles the DNA-bound structure. These results highlight the remarkable flexibility of λ Cro but also suggest that sequence-specific DNA binding may not induce large changes in the protein structure.     

A role for indels in the evolution of Cro protein folds

Insertions and deletions in protein sequences, or indels, can disrupt structure and may result in changes in protein folds during evolution or in association with alternative splicing. Pfl 6 and Xfaso 1 are two proteins in the Cro family that share a common ancestor but have different folds. Sequence alignments of the two proteins show two gaps, one at the N terminus, where the sequence of Xfaso 1 is two residues shorter, and one near the center of the sequence, where the sequence of Pfl 6 is five residues shorter. To test the potential importance of indels in Cro protein evolution, we generated hybrid variants of Pfl 6 and Xfaso 1 with indels in one or both regions, chosen according to several plausible sequence alignments. All but one deletion variant completely unfolded both proteins, showing that a longer N‐terminal sequence was critical for Pfl 6 folding and a longer central region sequence was critical for Xfaso 1 folding. By contrast, Xfaso 1 tolerated a longer N‐terminal sequence with little destabilization, and Pfl 6 tolerated central region insertions, albeit with substantial effects on thermal stability and some perturbation of the surrounding structure. None of the mutations appeared to convert one stable fold into the other. On the basis of this two‐protein comparison, short insertion and deletion mutations probably played a role in evolutionary fold change in the Cro family, but were also not the only factors. Proteins 2013; 81:1988–1996. © 2013 Wiley Periodicals, Inc.     

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.     

Crystal structure of an engineered Cro monomer bound nonspecifically to DNA: possible implications for nonspecific binding by the wild-type protein.

The structure has been determined at 3.0 A resolution of a complex of engineered monomeric Cro repressor with a seven-base pair DNA fragment. Although the sequence of the DNA corresponds to the consensus half-operator that is recognized by each subunit of the wild-type Cro dimer, the complex that is formed in the crystals by the isolated monomer appears to correspond to a sequence-independent mode of association. The overall orientation of the protein relative to the DNA is markedly different from that observed for Cro dimer bound to a consensus operator. The recognition helix is rotated 48 degrees further out of the major groove, while the turn region of the helix-turn-helix remains in contact with the DNA backbone. All of the direct base-specific interactions seen in the wild-type Cro-operator complex are lost. Virtually all of the ionic interactions with the DNA backbone, however, are maintained, as is the subset of contacts between the DNA backbone and a channel on the protein surface. Overall, 25% less surface area is buried at the protein DNA interface than for half of the wild-type Cro-operator complex, and the contacts are more ionic in character due to a reduction of hydrogen bonding and van der Waals interactions. Based on this crystal structure, model building was used to develop a possible model for the sequence-nonspecific interaction of the wild-type Cro dimer with DNA. In the sequence-specific complex, the DNA is bent, the protein dimer undergoes a large hinge-bending motion relative to the uncomplexed form, and the complex is twofold symmetric. In contrast, in the proposed nonspecific complex the DNA is straight, the protein retains a conformation similar to the apo form, and the complex lacks twofold symmetry. The model is consistent with thermodynamic, chemical, and mutagenic studies, and suggests that hinge bending of the Cro dimer may be critical in permitting the transition from the binding of protein at generic sites on the DNA to binding at high affinity operator sites.     

Insights into dimerization and four-helix bundle formation found by dissection of the dimer interface of the GrpE protein from Escherichia coli

The GrpE heat shock protein from Escherichia coli has a homodimeric structure. The dimer interface encompasses two long alpha-helices at the NH(2)-terminal end from each monomer (forming a "tail"), which lead into a small four-helix bundle from which each monomer contributes two short sequential alpha-helices in an antiparallel topological arrangement. We have created a number of different deletion mutants of GrpE that have portions of the dimer interface to investigate requirements for dimerization and to study four-helix bundle formation. Using chemical crosslinking and analytical ultracentrifugation techniques to probe for multimeric states, we find that a mutant containing only the long alpha-helical tail portion (GrpE1-88) is unable to form a dimer, most likely due to a decrease in alpha-helical content as determined by circular dichroism spectroscopy, thus one reason for a dimeric structure for the GrpE protein is to support the tail region. Mutants containing both of the short alpha-helices (GrpE1-138 and GrpE88-197) are able to form a dimer and presumably the four-helix bundle at the dimer interface. These two mutants have equilibrium constants for the monomer-dimer equilibrium that are very similar to the full-length protein suggesting that the tail region does not contribute significantly to the stability of the dimer. Interestingly, one mutant that contains just one of the short alpha-helices (GrpE1-112) exists as a tetrameric species, which presumably is forming a four-helix bundle structure. A proposed model is discussed for this mutant and its relevance for factors influencing four-helix bundle formation.     

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.     

Crystal structure of N-domain of FKBP22 from Shewanella sp. SIB1: dimer dissociation by disruption of Val-Leu knot

FK506‐binding protein 22 (FKBP22) from the psychrotophic bacterium Shewanella sp. SIB1 (SIB1 FKBP22) is a homodimeric protein with peptidyl prolyl cis‐trans isomerase (PPIase) activity. Each monomer consists of the N‐terminal domain responsible for dimerization and C‐terminal catalytic domain. To reveal interactions at the dimer interface of SIB1 FKBP22, the crystal structure of the N‐domain of SIB1 FKBP22 (SN‐FKBP22, residues 1‐68) was determined at 1.9 Å resolution. SN‐FKBP22 forms a dimer, in which each monomer consists of three helices (α1, α2, and α3N). In the dimer, two monomers have head‐to‐head interactions, in which residues 8–64 of one monomer form tight interface with the corresponding residues of the other. The interface is featured by the presence of a Val‐Leu knot, in which Val37 and Leu41 of one monomer interact with Val41 and Leu37 of the other, respectively. To examine whether SIB1 FKBP22 is dissociated into the monomers by disruption of this knot, the mutant protein V37R/L41R‐FKBP22, in which Val37 and Leu41 of SIB1 FKBP22 are simultaneously replaced by Arg, was constructed and biochemically characterized. This mutant protein was indistinguishable from the SIB1 FKBP22 derivative lacking the N‐domain in oligomeric state, far‐UV CD spectrum, thermal denaturation curve, PPIase activity, and binding ability to a folding intermediate of protein, suggesting that the N‐domain of V37R/L41R‐FKBP22 is disordered. We propose that a Val‐Leu knot at the dimer interface of SIB1 FKBP22 is important for dimerization and dimerization is required for folding of the N‐domain.