Structure of the cro repressor from bacteriophage lambda and its interaction with DNA

Structure of peptide sex pheromone receptor PrgX and PrgX/pheromone complexes and regulation of conjugation in Enterococcus faecalis

Many bacterial activities, including expression of virulence factors, horizontal genetic transfer, and production of antibiotics, are controlled by intercellular signaling using small molecules. To date, understanding of the molecular mechanisms of peptide-mediated cell-cell signaling has been limited by a dearth of published information about the molecular structures of the signaling components. Here, we present the molecular structure of PrgX, a DNA- and peptide-binding protein that regulates expression of the conjugative transfer genes of the Enterococcus faecalis plasmid pCF10 in response to an intercellular peptide pheromone signal. Comparison of the structures of PrgX and the PrgX/pheromone complex suggests that pheromone binding destabilizes PrgX tetramers, opening a 70-bp pCF10 DNA loop required for conjugation repression.

Design of λ 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.

Exo-Taq-Based Detection of DNA-Binding Protein for Homogeneous and Microarray Format

The study of DNA-protein interactions is of great importance to understand basic cellular processes such as transcription, replication and recombination. In this research, we developed a novel detection system for DNA-binding proteins (DBPs) involving the exonuclease (Exo) III and Taq DNA polymerase reactions. The system consists of three steps, as follows: the target DBP in the sample solution is incubated with probe DNA, and the probe is digested with Exo III and then extended with Taq using fluorescent dye-labeled dUTP as a substrate. The DBP protects the probe from digestion by Exo III. Therefore, only the DBP-bound probe allows the following extension. We examined this system using the lambda phage Cro repressor in a homogeneous format. The fluorescence image after gel electrophoresis showed a specific band. We also found that this system could be applied to the rapid and efficient detection of DBPs in stem and loop ds-DNA array formats. These results suggest that our method is useful as a new tool for analyzing DNA-protein interactions.

Detecting DNA-binding helix-turn-helix structural motifs using sequence and structure information

In this work, we analyse the potential for using structural knowledge to improve the detection of the DNA-binding helix–turn–helix (HTH) motif from sequence. Starting from a set of DNA-binding protein structures that include a functional HTH motif and have no apparent sequence similarity to each other, two different libraries of hidden Markov models (HMMs) were built. One library included sequence models of whole DNA-binding domains, which incorporate the HTH motif, the second library included shorter models of ‘partial’ domains, representing only the fraction of the domain that corresponds to the functionally relevant HTH motif itself. The libraries were scanned against a dataset of protein sequences, some containing the HTH motifs, others not. HMM predictions were compared with the results obtained from a previously published structure-based method and subsequently combined with it. The combined method proved more effective than either of the single-featured approaches, showing that information carried by motif sequences and motif structures are to some extent complementary and can successfully be used together for the detection of DNA-binding HTHs in proteins of unknown function.

Equilibrium unfolding of dimeric and engineered monomeric forms of lambda Cro (F58W) repressor and the effect of added salts: evidence for the formation of folded monomer induced by sodium perchlorate

The equilibrium unfolding transitions of Cro repressor variants, dimeric variant Cro F58W and monomer Cro K56[DGEVK]F58W, have been studied by urea and guanidine hydrochloride to probe the folding mechanism. The unfolding transitions of a dimeric variant are well described by a two state process involving native dimer and unfolded monomer with a free energy of unfolding, DeltaG(0,un)(0), of approximately 10-11 kcal/mol. The midpoint of transition curves is dependent on total protein concentration and DeltaG(0,un)(0) is independent of protein concentration, as expected for this model. Unfolding of Cro monomer is well described by the standard two state model. The stability of both forms of protein increases in the presence of salt but decreases with the decrease in pH. Because of the suggested importance of a N2<-->2F dimerization process in DNA binding, we have also studied the effect of sodium perchlorate, containing the chaotropic perchlorate anion, on the conformational transition of Cro dimer by CD, fluorescence and NMR (in addition to urea and guanidine hydrochloride) in an attempt both to characterize the thermodynamics of the process and to identify conditions that lead to an increase in the population of the folded monomers. Data suggest that sodium perchlorate stabilizes the protein at low concentration (<1.5 M) and destabilizes the protein at higher perchlorate concentration with the formation of a "significantly folded" monomer. The tryptophan residue in the "significantly folded" monomer induced by perchlorate is more exposed to the solvent than in native dimer.

Process of protein transport by the type III secretion system

The type III secretion system (TTSS) of gram-negative bacteria is responsible for delivering bacterial proteins, termed effectors, from the bacterial cytosol directly into the interior of host cells. The TTSS is expressed predominantly by pathogenic bacteria and is usually used to introduce deleterious effectors into host cells. While biochemical activities of effectors vary widely, the TTSS apparatus used to deliver these effectors is conserved and shows functional complementarity for secretion and translocation. This review focuses on proteins that constitute the TTSS apparatus and on mechanisms that guide effectors to the TTSS apparatus for transport. The TTSS apparatus includes predicted integral inner membrane proteins that are conserved widely across TTSSs and in the basal body of the bacterial flagellum. It also includes proteins that are specific to the TTSS and contribute to ring-like structures in the inner membrane and includes secretin family members that form ring-like structures in the outer membrane. Most prominently situated on these coaxial, membrane-embedded rings is a needle-like or pilus-like structure that is implicated as a conduit for effector translocation into host cells. A short region of mRNA sequence or protein sequence in effectors acts as a signal sequence, directing proteins for transport through the TTSS. Additionally, a number of effectors require the action of specific TTSS chaperones for efficient and physiologically meaningful translocation into host cells. Numerous models explaining how effectors are transported into host cells have been proposed, but understanding of this process is incomplete and this topic remains an active area of inquiry.

Integron integrase binds to bulged hairpin DNA

Gene cassettes are short, monogenic DNA elements that translocate between integrons through site-specific excision and integration. These events require that an integron-coded tyrosine recombinase forms a reactive complex with two sites, at least one of which belongs to the attC class. An attC site can be divided into two pairs of short repeats flanking a palindromic central region. The nucleotide sequence of attC among different cassettes varies extensively, implying that the site contains a structural recognition determinant with low sequence constraints. Oligonucleotides representing many different sequence modifications in either strand of the site were examined for integrase binding by using an electrophoresis mobility shift assay. The inner repeats, a central triplet and two single-nucleotide asymmetries in the site had the strongest influence on binding strength and strand choice. Our data show that the recombinase binds to a bulged hairpin in attC and that the hairpin distortion due to bulging could define the appropriate orientation of the otherwise symmetrical site. This is consistent with the strong bias for binding of recombinase to the bottom-strand oligonucleotides in vitro. Moreover, it was observed that the mobility-shifted complexes persisted under protein-denaturing assay conditions, indicating that a covalent link is indeed formed between integrase and DNA. Upon substitution of the presumed DNA-attacking residue, Y312, with a phenylalanine, DNA binding remained but there was no covalent linkage.

Comparison of backbone dynamics of monomeric and domain-swapped stefin A

Three-dimensional domain swapping has been observed in increasing number of proteins and has been implicated in the initial stages of protein aggregation, including that of the cystatins. Stefin A folds as a monomer under native conditions, while under some denaturing conditions domain-swapped dimer is formed. We have determined the backbone dynamics of the monomeric and domain-swapped dimeric forms of stefin A by (15)N relaxation using a model-free approach. The overall correlation times of the molecules were determined to be 4.6 +/- 0.1 ns and 9.2 +/- 0.2 ns for the monomer and the dimer, respectively. In the monomer, decreased order parameters indicate an increased mobility for the N-terminal trunk, the first and the second binding loops. At the opposite side of the molecule, the loop connecting the alpha-helix with strand B, the beginning of strand B and the loop connecting strands C and D show increased localized mobility. In the domain-swapped dimer, a distinctive feature of the structure is the concatenation of strands B and C into a single long beta-strand. The newly formed linker region between strands B and C, which substitutes for the first binding loop in the monomer, has order parameters typical for the remainder of the beta-strands. Thus, the interaction between subunits that occurs on domain-swapping has consequences for the dynamics of the protein at long-range from the site of conformational change, where an increased rigidity in the newly formed linker region is accompanied by an increased mobility of loops remote from that site.

Three-dimensional structure of MecI. Molecular basis for transcriptional regulation of staphylococcal methicillin resistance

Methicillin-resistant Staphylococcus aureus is the main cause of nosocomial and community-onset infections that affect millions of people worldwide. Some methicillin-resistant Staphylococcus aureus infections have become essentially untreatable by beta-lactams because of acquired molecular machineries enabling antibiotic resistance. Evasion from methicillin challenge is mainly achieved by the synthesis of a penicillin-binding protein of low affinity for antibiotics, MecA, that replaces regular penicillin-binding proteins in cell wall turnover when these have been inactivated by antibiotics. MecA synthesis is regulated by a signal transduction system consisting of the sensor/transducer MecR1 and the 14-kDa transcriptional repressor MecI (also known as methicillin repressor) that constitutively blocks mecA transcription. The three-dimensional structure of MecI reveals a dimer of two independent winged helix domains, each of which binds a palindromic DNA-operator half site, and two intimately intertwining dimerization domains of novel spiral staircase architecture, held together by a hydrophobic core. Limited proteolytic cleavage by cognate MecR1 within the dimerization domains results in loss of dimer interaction surface, dissociation, and repressor release, which triggers MecA synthesis. Structural information on components of the MecA regulatory pathway, in particular on methicillin repressor, the ultimate transcriptional trigger of mecA-encoded methicillin resistance, is expected to lead to the development of new antimicrobial drugs.

Design of peptide that recognizes double-stranded DNA

A novel molecular tool for double-stranded (ds) DNA detection using synthetic peptide is described. The peptide was designed based on the DNA binding domain of the lambda phage CRO repressor (CRO). The designed peptides contain helix-turn-helix (HTH), which is DNA binding motif. A cyclic peptide and a mutant peptide based on CRO were also designed, and the resulting affinity for dsDNA was increased. Furthermore, native amino acids of the peptide were replaced with arginine to increase the affinity for dsDNA. The affinity of these peptides for DNA binding was assessed by surface plasmon resonance (SPR) technique.

Nuclear transcription factors in the hippocampus

In the mammalian hippocampus, there is a trisynaptic loop that has been often referred to in studies on learning and memory mechanisms and their physiological correlate, the long-term potentiation (LTP). The three sets of synapses are formed by the fibers of perforant pathway terminating on granule cells and by the mossy fibers and Schaeffer collaterals making connections with the pyramidal cells. Each of the three types of synapses can develop LTP. LTP is accompanied by changes in gene expression and it is the nuclear transcription, involving specific transcription factors, that is the starting point for the series of biological amplifications and consolidations both necessary for such sustained changes. The transcription factors are proteins that control gene expression, development and functional formation in every eukaryotic cell. Two categories of transcription factors have been defined to date: general factors that comprise at least 20 proteins to form multiple preinitiation complex at the TATA box (TATA rich sequence) or regulatory factors that bind to promoter or enhancer regions at specific sites on the DNA close to, or distant from, the TATA box. Transcription factors have been divided into five different major classes according to unique protein motifs. These include basic domain, zinc-finger, helix-turn-helix, beta-Scaffold factors with minor groove contacts and other transcription factors not specifically classified. Much evidence has been accumulating in favor of the participation of several transcription factors in the consolidation of memory in the mammalian hippocampus following a spatial memory task. It is, therefore, of great importance that the involvement of transcription factors in de novo protein synthesis relevant to the synaptic mechanisms that mediate the formation of long-term memory should be summarized and discussed. No specific correlation between transduction of extracellular signals and expression of nuclear transcription factors, however, has been demonstrated to date.

3D domain swapping: as domains continue to swap

Three-dimensional (3D) domain swapping creates a bond between two or more protein molecules as they exchange their identical domains. Since the term '3D domain swapping' was first used to describe the dimeric structure of diphtheria toxin, the database of domain-swapped proteins has greatly expanded. Analyses of the now about 40 structurally characterized cases of domain-swapped proteins reveal that most swapped domains are at either the N or C terminus and that the swapped domains are diverse in their primary and secondary structures. In addition to tabulating domain-swapped proteins, we describe in detail several examples of 3D domain swapping which show the swapping of more than one domain in a protein, the structural evidence for 3D domain swapping in amyloid proteins, and the flexibility of hinge loops. We also discuss the physiological relevance of 3D domain swapping and a possible mechanism for 3D domain swapping. The present state of knowledge leads us to suggest that 3D domain swapping can occur under appropriate conditions in any protein with an unconstrained terminus. As domains continue to swap, this review attempts not only a summary of the known domain-swapped proteins, but also a framework for understanding future findings of 3D domain swapping.

Is there a code for protein-DNA recognition? Probab(ilistical)ly

Transcriptional regulation of all genes is initiated by the specific binding of regulatory proteins called transcription factors to specific sites on DNA called promoter regions. Transcription factors employ a variety of mechanisms to recognise their DNA target sites. In the last few decades, attempts have been made to describe these mechanisms by general sets of rules and associated models. We give an overview of these models, starting with a historical review of the somewhat controversial issue of a "recognition code" governing protein-DNA interaction. We then present a probabilistic framework in which advantages and disadvantages of various models can be discussed. Finally, we conclude that simplifying assumptions about additivity of interactions are sufficiently justified in many situations (and can be suitably extended in other situations) to allow a unifying concept of a "probabilistic code" for protein-DNA recognition to be defined.

Single-site mutations induce 3D domain swapping in the B1 domain of protein L from Peptostreptococcus magnus

BACKGROUND

Thermodynamic and kinetic studies of the Protein L B1 domain (Ppl) suggest a folding pathway in which, during the folding transition, the first beta hairpin is formed while the second beta hairpin and the alpha helix are largely unstructured. The same mutations in the two beta turns have opposite effects on the folding and unfolding rates. Three of the four residues composing the second beta turn in Ppl have consecutive positive phi angles, indicating strain in the second beta turn.

RESULTS

We have determined the crystal structures of the beta turn mutants G55A, K54G, and G15A, as well as a core mutant, V49A, in order to investigate how backbone strain affects the overall structure of Ppl. Perturbation of the hydrophobic interactions at the closed interface by the V49A mutation triggered the domain swapping of the C-terminal beta strand that relieved the strain in the second beta turn. Interestingly, the asymmetric unit of V49A contains two monomers and one domain-swapped dimer. The G55A mutation escalated the strain in the second beta turn, and this increased strain shifted the equilibrium toward the domain-swapped dimer. The K54G structure revealed that the increased stability is due to the reduction of strain in the second beta turn, while the G15A structure showed that increased strain alone is insufficient to trigger domain swapping.

CONCLUSIONS

Domain swapping in Ppl is determined by the balance of two opposing components of the free energy. One is the strain in the second beta turn that favors the dimer, and the other is the entropic cost of dimer formation that favors the monomer. A single-site mutation can disrupt this balance and trigger domain swapping.

Recognition of specific DNA sequences

Proteins that recognize specific DNA sequences play a central role in the regulation of transcription. The tremendous increase in structural information on protein-DNA complexes has uncovered a remarkable structural diversity in DNA binding folds, while at the same time revealing common themes in binding to target sites in the genome.

A verotoxin 1 B subunit-lambda CRO chimeric protein specifically binds both DNA and globotriaosylceramide (Gb(3)) to effect nuclear targeting of exogenous DNA in Gb(3) positive cells

Inefficient nuclear incorporation of foreign DNA remains a critical roadblock in the development of effective nonviral gene delivery systems. DNA delivered by traditional protocols remains within endosomal/lysosomal vesicles, or is rapidly degraded in the cytoplasm. Verotoxin I (VT), an AB(5) subunit toxin produced by enterohaemorrhagic Escherichia coli, binds to the cell surface glycolipid, globotriaosylceramide (Gb(3)) and is internalized into preendosomes. VT is then retrograde transported to the Golgi, endoplasmic reticulum (ER), and nucleus of highly VT-sensitive cells. We have utilized this nuclear targeting of VT to design a unique delivery system which transports exogenous DNA via vesicular traffic to the nucleus. The nontoxic VT binding subunit (VTB) was fused to the lambda Cro DNA-binding repressor, generating a 14-kDa VTB-Cro chimera. VTB-Cro binds specifically via the Cro domain to a 25-bp DNA fragment containing the consensus Cro operator. VTB-Cro demonstrates simultaneous specific binding to Gb(3). Treatment of Vero cells with fluorescent-labeled Cro operator DNA in the presence of VTB-Cro, results in DNA internalization to the Golgi, ER, and nucleus, whereas fluorescent DNA alone is incorporated poorly and randomly within the cytoplasm. VTB-Cro mediated nuclear DNA transport is prevented by brefeldin A, consistent with Golgi/ER intracellular routing. Pretreatment with filipin had no effect, indicating that caveoli are not involved. This novel VTB-Cro shuttle protein may find practical applications in the fields of intracellular targeting, gene delivery, and gene therapy.

Designed sequence-specific minor groove ligands

In the past decade, a general design for sequence-specific minor groove ligands has evolved, based on the natural products distamycin and netropsin. By utilizing a basic set of design rules for connecting pyrrole, imidazole, and hydroxypyrrole modules, new ligands can be prepared to target almost any sequence of interest with both high affinity and specificity. In this review we present the design rules with a brief history of how they evolved. The structural basis for sequence-specific recognition is explained, together with developments that allow linking of recognition modules that enable targeting of long DNA sequences. Examples of the affinity and specificity that can be achieved with a number of variations on the basic design are given. Recently these molecules have been used to compete with proteins both in vitro and in vivo, and a brief description of the experimental results are given.

Coupled energetics of lambda cro repressor self-assembly and site-specific DNA operator binding II: cooperative interactions of cro dimers

The bacteriophage lambda relies on interactions of the cI and cro repressors which self assemble and bind the two operators (O(R) and O(L)) of the phage genome to control the lysogenic to lytic switch. While the self assembly and O(R) binding of cI have been investigated in detail, a more complete understanding of gene regulation by phage lambda also requires detailed knowledge of the role of cro repressor as it dimerizes and binds at O(R) sites. Since dimerization and operator binding are coupled processes, a full elucidation of the regulatory energetics in this system requires that the equilibrium constants for dimerization and cooperative binding be determined. The dimerization constant for cro has been measured as a prelude to these binding studies. Here, the energetics of cro binding to O(R) are evaluated using quantitative DNaseI footprint titration techniques. Binding data for wild-type and modified O(R) site combinations have been simultaneously analyzed in concert with the dimerization energetics to obtain both the intrinsic and cooperative DNA binding energies for cro with the three O(R) sites. Binding of cro dimers is strongest to O(R)3, then O(R)1 and lastly, O(R)2. Adjacently bound repressors exhibit positive cooperativity ranging from -0.6 to -1.0 kcal/mol. Implications of these, newly resolved, energetics are discussed in the framework of a dynamic model for gene regulation. This characterization of the DNA-binding properties of cro repressor establishes the foundation on which the system can be explored for other, more complex, regulatory elements such as cI-cro cooperativity.

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 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.

DNA-induced conformational changes in bacteriophage 434 repressor

Although bacteriophage 434 repressor binds to its specific DNA sites only as a dimer, formation of the dimers in solution occurs at concentrations three orders of magnitude higher than those needed to bind the 434 operator DNA. Our results suggest that both specific and non-specific DNA induce conformational changes in repressor that lead to formation of repressor dimers. The repressor conformational changes induced by DNA occur at concentrations much lower than those needed for binding of repressor, suggesting that the alternative conformations of repressor persist even if the protein is not in direct contact with DNA. Hence, DNA acts in a "catalytic" fashion to induce a steady-state amount of an alternative repressor conformation that has an enhanced affinity for its specific binding site. These findings suggest that the repressor conformer induced by non-specific DNA is the form of the repressor that is optimized for searching for DNA binding sites along non-specific DNA. Upon finding a binding site, the repressor protein undergoes an additional conformational change that allows it to "lock-on" to its specific site.

Two-dimensional IR correlation spectroscopy: sequential events in the unfolding process of the lambda cro-V55C repressor protein

A question often posed in protein folding/unfolding studies is whether the process is fully cooperative or whether it contains sequential elements. To address this question, one needs tools capable of resolving different events. It seems that, at least in certain cases, two-dimensional (2D) IR correlation spectroscopy can provide answers to this question. To illustrate this point, we have turned to the Cro-V55C dimer of the lambda Cro repressor, a protein known to undergo thermal unfolding in two discrete steps through a stable equilibrium intermediate. The secondary structure of this intermediate is compatible with that of a partially unfolded protein and involves a reorganization of the N terminus, whereas the antiparallel beta-ribbon formed by the C-terminal part of each subunit remains largely intact. To establish whether the unfolding process involves sequential events, we have performed a 2D correlation analysis of IR spectra recorded over the temperature range of 20-95 degrees C. The 2D IR correlation analysis indeed provides evidence for a sequential formation of the stable intermediate, which is created in three (closely related) steps. A first step entails the unfolding of the short N-terminal beta-strand, followed by the unfolding of the alpha-helices in a second step, and the third step comprises the reorganization of the remaining beta-sheet and of some unordered segments in the protein. The complete unfolding of the stable intermediate at higher temperatures also undergoes sequential events that ultimately end with the breaking of the H bonds between the two beta-strands at the dimer interface.

A2 cro, the lysogenic cycle repressor, specifically binds to the genetic switch region of Lactobacillus casei bacteriophage A2

Lysogenic induction of temperate bacteriophage A2 of Lactobacillus casei is controlled by the action of its cI and cro products at the phage operator region. Three 20-bp inverted repeated DNA segments (subsites O1, O2, and O3) and the two divergent (PL and PR) promoters were mapped within the 153-bp operator region. The A2-encoded Cro product is shown to be the functional homolog of lambda Cro. The binding of Cro to the three operator subsites is noncooperative and yields two discrete protein-DNA complexes of retarded migration in mobility shift assays. The Kapp value for the Cro-PL-PR DNA complex was estimated to be 6 nM. Cro shows a slightly higher affinity for O3 than for O1 and O2 subsites. The O3 subsite overlaps the -35 hexamer of the PL promoter, which directs cI expression. A Cro mutant protein, devoid of the last 12 residues (Cro*), allowed the assignment of the DNA-binding domain to the NH2 end of Cro. The C end enhances its affinity for the DNA and probably stabilizes bending induced by Cro.

Multiprotein-DNA complexes in transcriptional regulation

Transcription in eukaryotes is frequently regulated by a mechanism termed combinatorial control, whereby several different proteins must bind DNA in concert to achieve appropriate regulation of the downstream gene. X-ray crystallographic studies of multiprotein complexes bound to DNA have been carried out to investigate the molecular determinants of complex assembly and DNA binding. This work has provided important insights into the specific protein-protein and protein-DNA interactions that govern the assembly of multiprotein regulatory complexes. The results of these studies are reviewed here, and the general insights into the mechanism of combinatorial gene regulation are discussed.

Dimerisation mutants of Lac repressor. I. A monomeric mutant, L251A, that binds Lac operator DNA as a dimer

Dimer formation between monomers of the Escherichia coli Lac repressor is substantially specificed by the interactions between three alpha-helices in each monomer which form a hydrophobic interface. As a first step in analysing the specificity of this interaction, we examined the mutant L251A. LacR bearing this mutation in a background lacking the C-terminal heptad repeats is completely incapable of forming dimers in solution, with a dimer-monomer equilibrium dissociation constant, or Kd, higher than 10(-5)M. This correlates with a 200-fold decrease in its ability to repress the lac operon in vivo compared to dimeric LacR. Surprisingly, the mutant is still capable of forming dimers upon binding to short operator DNA in vitro. Analysis of the kinetic parameters of binding of the mutant to operator DNA reveals a 2000 to 3000-fold increase in the equilibrium dissociation constant (Kd) of the mutant-DNA complex in comparison to dimeric LacR-operator complexes, with the change almost entirely due to a greater than 1000-fold decrease in association rate. The dissociation rate varies only by a factor of about two, in comparison to dimeric LacR. This change reflects a kinetic pathway in which dimer formation, in solution or on DNA, is the rate-limiting step. These findings have implications for the specificity and stability of the protein-protein interface in question.

Intermolecular contacts between the lambda-Cro repressor and the operator DNA characterized by nuclear magnetic resonance spectroscopy

The specific interaction between lambda phage Cro repressor and the DNA fragment bearing the consensus sequence of operators has been studied using nuclear magnetic resonance (NMR). Using both 15N- and 13C/15N- labeled lambda-Cro in complex with unlabeled DNA, chemical shift assignments of the lambda-Cro-DNA complex were obtained using heteronuclear NMR experiments. Inter-molecular contacts between the protein and DNA were identified using heteronuclear filtered NOESY experiments. The inter-molecular contacts were supplemented with intra-protein and intra-DNA NOE constraints to dock lambda-Cro to the bent B-form DNA using restrained molecular dynamics. The structure of one of the subunits of lambda-Cro in the complex is essentially the same as that of the unbound form. In the complex, inter-molecular NOEs were observed between the "helix-turn-helix" region comprising the alpha2 and alpha3 helices of the lambda-Cro protein and the major groove of the DNA. The methyl group of Thr17 forms a hydrophobic contact with the methyl group of the thymine at base pair 1 in the DNA, and Val25 and Ala29 make hydrophobic contacts with the methyl group of the thymine at base pair 5. The presence and the absence of these contacts can explain the difference in the affinity of lambda-Cro to several variants of the operator sequence.

The structure of plasmid-encoded transcriptional repressor CopG unliganded and bound to its operator

The structure of the 45 amino acid transcriptional repressor, CopG, has been solved unliganded and bound to its target operator DNA. The protein, encoded by the promiscuous streptococcal plasmid pMV158, is involved in the control of plasmid copy number. The structure of this protein repressor, which is the shortest reported to date and the first isolated from a plasmid, has a homodimeric ribbon-helix-helix arrangement. It is the prototype for a family of homologous plasmid repressors. CopG cooperatively associates, completely protecting several turns on one face of the double helix in both directions from a 13-bp pseudosymmetric primary DNA recognition element. In the complex structure, one protein tetramer binds at one face of a 19-bp oligonucleotide, containing the pseudosymmetric element, with two beta-ribbons inserted into the major groove. The DNA is bent 60 degrees by compression of both major and minor grooves. The protein dimer displays topological similarity to Arc and MetJ repressors. Nevertheless, the functional tetramer has a unique structure with the two vicinal recognition ribbon elements at a short distance, thus inducing strong DNA bend. Further structural resemblance is found with helix-turn-helix regions of unrelated DNA-binding proteins. In contrast to these, however, the bihelical region of CopG has a role in oligomerization instead of DNA recognition. This observation unveils an evolutionary link between ribbon-helix-helix and helix-turn-helix proteins.

An evolutionary link between sporulation and prophage induction in the structure of a repressor:anti-repressor complex

Spore formation is an extreme response of some bacteria to adversity. In Bacillus subtilis the proteins of the sin, sporulation inhibition, region form a component of an elaborate molecular circuitry that regulates the commitment to sporulation. SinR is a tetrameric repressor protein that binds to the promoters of genes essential for entry into sporulation and prevents their transcription. This repression is overcome through the activity of SinI, which disrupts the SinR tetramer through the formation of a SinI-SinR heterodimer. The interactions governing this curious quaternary transition are revealed in the crystal structure of the SinI-SinR complex. The most striking, and unexpected, finding is that the tertiary structure of the DNA-binding domain of SinR is identical with that of the corresponding domains of the repressor proteins, CI and Cro, of bacteriophage 434 that regulate lysis/lysogeny. This structural similarity greatly exceeds that between SinR and any bacterial protein or between the 434 repressor proteins and their homologues in the closely related bacteriophage lambda. The close evolutionary relationship implied by the structures of SinR and the 434 repressors provokes both comparison of their functions and a speculative consideration of the intriguing possibility of an evolutionary link between the two adaptive responses, sporulation and prophage induction.

Engineering an intertwined form of CD2 for stability and assembly

The amino-terminal domain of CD2 has the remarkable ability to fold in two ways: either as a monomer or as an intertwined, metastable dimer. Here we show that it is possible to differentially stabilize either fold by engineering the CD2 sequence, mimicking random mutagenesis events that could occur during molecular evolution. Crystal structures of a hinge-deletion mutant, which is stable as an intertwined dimer, reveal domain rotations that enable the protein to further assemble to a tetramer. These results demonstrate that a variety of folds can be adopted by a single polypeptide sequence, and provide guidance for the design of proteins capable of further assembly.

Structural features of the plasmid pMV158-encoded transcriptional repressor CopG, a protein sharing similarities with both helix-turn-helix and beta-sheet DNA binding proteins

The small transcriptional repressor CopG protein (45 amino acids) encoded by the streptococcal plasmid pMV158 was purified to near homogeneity. Gel filtration chromatography and analytical ultracentrifugation showed that the native protein is a spherical dimer of identical subunits. Circular dichroism measurements of CopG indicated a consensus average content of more than 50% alpha-helix and 10-35% beta-strand and turns, which is compatible with the predicted secondary structure of the protein. CopG exhibited a prolonged intracellular half-life, but deletions in regions other than the C-terminal affected the global structure of the protein, severely reducing the half-lives of the CopG variants. This indicates that CopG has a compact structure, perhaps constituted by a single domain. Molecular modeling of CopG showed a good fitting between the helix-turn-helix motifs of well-known repressor proteins and a bihelical unit of CopG. However, modeling of CopG with ribbon-helix-helix class of DNA binding proteins also exhibited an excellent fit. Eleven out of the 12 replicons belonging to the pMV158 plasmid family could also encode Cop proteins, which share features with both helix-turn-helix and beta-sheet DNA binding proteins.

Refined structure of Cro repressor protein from bacteriophage lambda suggests both flexibility and plasticity

The structure of the Cro repressor protein from phage lambda has been refined to a crystallographic R-value of 19.3% at 2.3 A resolution. The re fined model supports the structure as originally described in 1981 and provides a basis for comparison with the Cro-operator complex described in the accompanying paper. Changes in structure seen in different crystal forms and modifications of Cro suggest that the individual subunits are somewhat plastic in nature. In addition, the dimer of Cro suggests a high degree of flexibility, which may be important in forming the Cro-DNA complex. The structure of the Cro subunit as determined by NMR agrees reasonably well with that in the crystals (root-mean-square discrepancy of about 2 A for all atoms). There are, however, only a limited number of intersubunit distance constraints and, presumably for this reason, the different NMR models for the dimer vary substantially among themselves (discrepancies of 1.3 to 5.5 A). Because of this variation it is not possible to say whether the range of discrepancies between the X-ray and NMR Cro dimers (2.9 to 7.5 A) represent a significant difference between the X-ray and solution structures. It has previously been proposed that substitutions of Tyr26 in Cro increase thermal stability by the "reverse hydrophobic effect", i.e. by exposing 40% more hydrophobic surface to solvent in the folded form than in the unfolded state. The refined structure, however, suggests that Tyr26 is equally solvent exposed in the folded and unfolded states. The most stabilizing substitution is Tyr26-->Asp and in this case it appears that interaction with an alpha-helix dipole is at least partly responsible for the enhanced stability.

Crystal structure of lambda-Cro bound to a consensus operator at 3.0 A resolution

The structure of the Cro protein from bacteriophage lambda in complex with a 19 base-pair DNA duplex that includes the 17 base-pair consensus operator has been determined at 3.0 A resolution. The structure confirms the large changes in the protein and DNA seen previously in a crystallographically distinct low-resolution structure of the complex and, for the first time, reveals the detailed interactions between the side-chains of the protein and the base-pairs of the operator. Relative to the crystal structure of the free protein, the subunits of Cro rotate 53 degrees with respect to each other on binding DNA. At the same time the DNA is bent by 40 degrees through the 19 base-pairs. The intersubunit connection includes a region within the protein core that is structurally reminiscent of the "ball and socket" motif seen in the immunoglobulins and T-cell receptors. The crystal structure of the Cro complex is consistent with virtually all available biochemical and related data. Some of the interactions between Cro and DNA proposed on the basis of model-building are now seen to be correct, but many are different. Tests of the original model by mutagenesis and biochemical analysis corrected some but not all of the errors. Within the limitations of the crystallographic resolution it appears that operator recognition is achieved almost entirely by direct hydrogen-bonding and van der Waals contacts between the protein and the exposed bases within the major groove of the DNA. The discrimination of Cro between the operators OR3 and OR1, which differ in sequence at just three positions, is inferred to result from a combination of small differences, both favorable and unfavorable. A van der Waals contact at one of the positions is of primary importance, while the other two provide smaller, indirect effects. Direct hydrogen bonding is not utilized in this distinction.

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.

A new DNA-binding motif in the Skn-1 binding domain-DNA complex

The DNA-binding domain of Skn-1, a developmental transcription factor that specifies mesoderm in C. elegans, is shown by X-ray crystallography to have a novel fold in which a compact, monomeric, four-helix unit organizes two DNA-contact elements. At the C-terminus, a helix extends from the domain to occupy the major groove of DNA in a manner similar to bZip proteins. Skn-1, however, lacks the leucine zipper found in all bZips. Additional contacts with the DNA are made by a short basic segment at the N-terminus of the domain, reminiscent of the 'homeodomain arm'.

Single-chain lambda Cro repressors confirm high intrinsic dimer-DNA affinity

The overall affinity of the bacteriophage lambda Cro repressor for its operator DNA site is limited by dimer dissociation at submicromolar concentrations. Since Cro dimer-operator complexes form at nanomolar concentrations of Cro subunits where free dimers are rare, these dimers must bind with compensating high affinities. Previous studies of the covalent dimer Cro V55C suggest little change in DNA binding affinity even though the dimeric species is quantitatively populated; this is an apparent contradiction to the expectation of high intrinsic dimer-DNA affinity. In contrast to the disulfide linkage at the center of the dimer interface in Cro V55C, polypeptide linkers that join the two subunits allow single-chain Cro repressors to bind operator DNA with picomolar affinities. A series of five single-chain Cro repressors have been expressed from fused tandem cro genes. Each contains a peptide linker of 8-16 hydrophilic residues that connects the C-terminus of one subunit to the N-terminus of the next. All bind to operator DNA with at least 100-fold higher affinity than Cro V55C. Proteins containing the longest and shortest linkers have been purified and characterized in detail. Both exhibit similar CD spectra to wild-type Cro and enhanced thermal stability. Sedimentation equilibrium experiments show that single-chain Cro repressors do not associate at concentrations up to 30 microM. The rate of dissociation of Cro-DNA complexes is almost unchanged by covalent linkage. Biophysical characterization of Cro variants such as these, where DNA binding is uncoupled from subunit assembly, is necessary for a quantitative understanding of the structural and energetic determinants of DNA recognition in this simple model system.

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.

How Cro and lambda-repressor distinguish between operators: the structural basis underlying a genetic switch

Knowledge of the three-dimensional structures of the lambda-Cro and lambda-repressor proteins in complex with DNA has made it possible to evaluate how these proteins discriminate between different operators in phage lambda. As anticipated in previous studies, the helix-turn-helix units of the respective proteins bind in very different alignments. In Cro the recognition helices are 29 A apart and are tilted by 55 degrees with respect to each other, but bind parallel to the major groove of the DNA. In lambda-repressor [Beamer, L. J. & Pabo, C. O. (1992) J. Mol. Biol. 227, 177-196] the helices are 34 A apart and are essentially parallel to each other, but are inclined to the major grooves. The DNA is much more bent when bound by Cro than in the case with lambda-repressor. The first two amino acids of the recognition helices of the two proteins, Gln-27 and Ser-28 in Cro, and Gln-44 and Ser-45 in lambda-repressor, make very similar interactions with the invariant bps 2 and 4. There are also analogous contacts between the thymine of bp 5 and, respectively, the backbone of Ala-29 of Cro and the backbone of Gly-46 of lambda-repressor. Otherwise, however, unrelated parts of the two proteins are used in sequence-specific recognition. It appears that similar contacts to the invariant or almost invariant bps (especially 2 and 4) are used by both Cro and lambda-repressor to differentiate the operator sites as a group from other sites on the DNA. The discrimination of Cro and lambda-repressor between their different operators is more subtle and seems to be achieved primarily through differences in van der Waals contacts at bp 3', together with weaker, less direct effects at bps 5' and 8', all in the nonconsensus half of the operators. The results provide further support for the idea that there is no simple code for DNA-protein recognition.

Arac/XylS family of transcriptional regulators

The ArC/XylS family of prokaryotic positive transcriptional regulators includes more than 100 proteins and polypeptides derived from open reading frames translated from DNA sequences. Members of this family are widely distributed and have been found in the gamma subgroup of the proteobacteria, low- and high-G + C-content gram-positive bacteria, and cyanobacteria. These proteins are defined by a profile that can be accessed from PROSITE PS01124. Members of the family are about 300 amino acids long and have three main regulatory functions in common: carbon metabolism, stress response, and pathogenesis. Multiple alignments of the proteins of the family define a conserved stretch of 99 amino acids usually located at the C-terminal region of the regulator and connected to a nonconserved region via a linker. The conserved stretch contains all the elements required to bind DNA target sequences and to activate transcription from cognate promoters. Secondary analysis of the conserved region suggests that it contains two potential alpha-helix-turn-alpha-helix DNA binding motifs. The first, and better-fitting motif is supported by biochemical data, whereas existing biochemical data neither support nor refute the proposal that the second region possesses this structure. The phylogenetic relationship suggests that members of the family have recruited the nonconserved domain(s) into a series of existing domains involved in DNA recognition and transcription stimulation and that this recruited domain governs the role that the regulator carries out. For some regulators, it has been demonstrated that the nonconserved region contains the dimerization domain. For the regulators involved in carbon metabolism, the effector binding determinants are also in this region. Most regulators belonging to the AraC/XylS family recognize multiple binding sites in the regulated promoters. One of the motifs usually overlaps or is adjacent to the -35 region of the cognate promoters. Footprinting assays have suggested that these regulators protect a stretch of up to 20 bp in the target promoters, and multiple alignments of binding sites for a number of regulators have shown that the proteins recognize short motifs within the protected region.

The role of DNA bending in Cro protein-DNA interactions

Binding energy of DNA-Cro protein complexes is analyzed in terms of DNA elasticity, using a sequence-dependent anisotropic bendability (SDAB) model of DNA, developed recently [M.M. Gromiha, M.G. Munteanu, A. Gabrielian and S. Pongor, J. Biol. Phys. 22(1996) 227-243.]. The protein is considered to bind aspecifically to DNA that reduces the freedom of movement in the DNA molecule. In cognate DNA, the Cro protein moves on to form specific interactions and bends DNA. A comparison of the experimental data [Y. Takeda, A. Sarai and V.M. Rivera, Proc. Natl. Acad. Sci. U.S.A. 86 (1989) 439-443.] with the calculated DNA stiffness data shows that delta G of the complex formation increases with stiffness of the ligand when the interactions are nonspecific ones, while an opposite trend is observed for specific binding. Both of these trends are in agreement with our approach using the SDAB model. A decomposition of the energy terms suggests that binding energy in the nonspecific case is used maily to compensate the free energy changes due to entropy lost by DNA, while the energy of specific interactions provide enough energy both to bend the DNA molecule and to change the conformation of the Cro protein upon ligand binding.

A folded monomeric intermediate in the formation of lambda Cro dimer-DNA complexes

The folding, dimerization and DNA binding equilibria of the bacteriophage lambda Cro repressor have been characterized. Comparison with four engineered variants shows that a folded monomeric species is substantially populated under conditions used for the formation of dimer-DNA complexes. Although Cro dimers are the only DNA-bound species observed in electrophoretic mobility shift assays, cooperativity in Cro-DNA binding isotherms shows that the predominant free protein species is monomeric at nanomolar concentrations. Micromolar dissociation constants for Cro dimers have been measured in the absence of DNA by sedimentation equilibrium and gel filtration chromatography. Denaturation of Cro dimers in the 10 to 100 micromolar concentration range by guanidine hydrochloride (GdnHCl) is well modeled as a two-state process, with folded dimers and unfolded monomers as the only significantly populated species. However, linear extrapolation of this composite unfolding and dimer dissociation free energy predicts a nanomolar dissociation constant in the absence of denaturant. This extrapolation is clearly inconsistent with the DNA binding and hydrodynamic measurements. Our interpretation of these results is that the monomeric species detected in DNA binding and hydrodynamic experiments is predominantly folded. The stability of the folded monomeric species can be calculated as the difference between the dimerization free energy determined from hydrodynamic measurements and the folding free energy extrapolated from GdnHCl denaturation. The calculated stability of the Cro F58W monomer is greater than that of the wild-type Cro monomer. Thus, residue 58, which makes critical intermolecular contacts across the dimer interface, is also involved in intramolecular stabilization of the monomeric intermediate.

Ribosomal protein S7: a new RNA-binding motif with structural similarities to a DNA architectural factor

BACKGROUND

The ribosome is a ribonucleoprotein complex which performs the crucial function of protein biosynthesis. Its role is to decode mRNAs within the cell and to synthesize the corresponding proteins. Ribosomal protein S7 is located at the head of the small (30S) subunit of the ribosome and faces into the decoding centre. S7 is one of the primary 16S rRNA-binding proteins responsible for initiating the assembly of the head of the 30S subunit. In addition, S7 has been shown to be the major protein component to cross-link with tRNA molecules bound at both the aminoacyl-tRNA (A) and peptidyl-tRNA (P) sites of the ribosome. The ribosomal protein S7 clearly plays an important role in ribosome function. It was hoped that an atomic-resolution structure of this protein would aid our understanding of ribosomal mechanisms.

RESULTS

The structure of ribosomal protein S7 from Bacillus stearothermophilus has been solved at 2.5 A resolution using multiwavelength anomalous diffraction and selenomethionyl-substituted proteins. The molecule consists of a helical hydrophobic core domain and a beta-ribbon arm extending from the hydrophobic core. The helical core domain is composed of a pair of entangled helix-turn-helix motifs; the fold of the core is similar to that of a DNA architectural factor. Highly conserved basic and aromatic residues are clustered on one face of the S7 molecule and create a 16S rRNA contact surface.

CONCLUSIONS

The molecular structure of S7, together with the results of previous cross-linking experiments, suggest how this ribosomal protein binds to the 3' major domain of 16S rRNA and mediates the folding of 16S rRNA to create the ribosome decoding centre.

Conservative Val47 residue of POU homeodomain: role in DNA recognition

Conservative Val47 residue, located in the third recognition helix of the Oct-2 POU domain, was alternately substituted with other 19 amino acids. Affinity and specificity of interaction with oct-site ATGCAAANGA and homeo-specific site ATAANGA were determined for all mutants. The wild type protein (with Val47) has maximal affinity and specificity in POU domain interaction with octamer sequence. However, V47I mutant showed stronger interaction with homeo-specific site. The highest specificity of interaction with homeo-site was recorded for V47S mutant. We conclude that only Val47 provides sequence-specific high-affinity binding of POU proteins with octamer targets other than the homeo-specific site. It is shown also that damages caused by point mutations may be at least partially compensated by participation in the oct-site recognition of both POUh and POUs domains.

Structural trees for protein superfamilies

Structural trees for large protein superfamilies, such as beta proteins with the aligned beta sheet packing, beta proteins with the orthogonal packing of alpha helices, two-layer and three-layer alpha/beta proteins, have been constructed. The structural motifs having unique overall folds and a unique handedness are taken as root structures of the trees. The larger protein structures of each superfamily are obtained by a stepwise addition of alpha helices and/or beta strands to the corresponding root motif, taking into account a restricted set of rules inferred from known principles of the protein structure. Among these rules, prohibition of crossing connections, attention to handedness and compactness, and a requirement for alpha helices to be packed in alpha-helical layers and beta strands in beta layers are the most important. Proteins and domains whose structures can be obtained by stepwise addition of alpha helices and/or beta strands to the same root motif can be grouped into one structural class or a superfamily. Proteins and domains found within branches of a structural tree can be grouped into subclasses or subfamilies. Levels of structural similarity between different proteins can easily be observed by visual inspection. Within one branch, protein structures having a higher position in the tree include the structures located lower. Proteins and domains of different branches have the structure located in the branching point as the common fold.