Anatomy of specific interactions between lambda repressor and operator DNA

Recognition of specific DNA sequences by proteins is essential for regulation of gene expression. To fully understand the recognition mechanism, it is necessary to understand not only the structure of the specific protein-DNA interactions but also the energetics. We therefore performed a computer analysis in which a phage DNA-binding protein, lambda repressor, was used to examine the changes in binding free energy (DeltaDeltaG) and its energy components caused by single base mutations. We then determined which of the calculated energy components best correlated with the experimental data. The experimental DeltaDeltaG values were well reproduced by the calculations. Component analysis revealed that the electrostatic and hydrogen bond energies were most strongly correlated with the experimental data. Among the 51 single base-substitution mutants examined, positive DeltaDeltaG values, corresponding to weakened binding, were caused by the loss of favorable electrostatic interactions and hydrogen bonds, the introduction of steric collisions and electrostatic repulsion, the loss of favorable interactions with a thymine methyl group, and the increase of unfavorable hydration energy from isolated DNA. This analysis also showed distinct patterns of recognition at A-T and G-C positions, as different combinations of energy components were involved in DeltaDeltaG caused by the two substitution types. We have thus been able to identify the energy components that most strongly correlate with sequence-dependent DeltaDeltaG and determine their contribution to the specificity of DNA sequence recognition by the lambda repressor. Application of this method to other systems should provide additional insight into the molecular mechanism of protein-DNA recognition.     

Structure Based Sequence Dependent Stiffness Scale for Trinucleotides: A Direct Method

A new set of stiffness parameters for all the 32trinucleotide units has been set up directly from thethree dimensional structures of DNA molecules. It wasobserved that GAC/GTC is the stiffest trinucleotideand ACC/GGT is the most flexible one. The averagestiffness values computed for a set of operatorsequences using the new parameters correlate very wellwith the protein-DNA binding specificity and bindingfree energy change of 434 repressor and Cro repressor,respectively. The new structure based stiffness scalecan explain the protein-DNA binding specificity to thelevel of 0.92.     

Two-hybrid system for characterization of protein-protein interactions in E. coli

The yeast two-hybrid system has been used to characterize many protein-protein interactions. A two-hybrid system for E. coli was constructed in which one hybrid protein bound to a specific DNA site recruits another to an adjacent DNA binding site. The first hybrid comprises a test protein, the bait, fused to a chimeric protein containing the 434 repressor DNA binding domain. In the second hybrid, a second test protein, the prey, is fused downstream of a chimeric protein with the DNA binding specificity of the lambda repressor. Reporters were designed to express cat and lacZ under the control of a low-affinity lambda operator. At low expression levels, lambda repressor hybrids weakly repress the reporter genes. A high-affinity operator recognized by 434 repressor was placed nearby, in a position that does not yield repression by 434 repressor alone. If the test proteins interact, the 434 hybrid bound to the 434 operator stabilizes the binding of the lambda repressor hybrid to the lambda operator, causing increased repression of the reporter genes. Reconstruction experiments with the fos and jun leucine zippers detected protein-protein interactions between either homodimeric or heterodimeric leucine zippers.     

High-specific-activity probes and a high-resolution in-gel photo cross-linking assay for protein-DNA complexes

Photochemical cross-linking has been widely employed to identify proteins interacting with specific sites on DNA. Identification of bound proteins usually relies on transfer of a radiolabel from the DNA to the protein by cross-linking. We set out to fine-map a small viral replication preinitiation complex composed of two protein dimers bound to DNA, the bovine papillomavirus E1E2-ori complex. Here we describe a simple method for generating high-specific-activity probes with a phenyl-azide photoactivatible cross-linking group positioned immediately adjacent to a labeled nucleotide. The method is based on the selective destruction of one 5'-phosphorylated strand of a polymerase chain reaction product with lambda exonuclease and reconstitution of the probe with a phosphorothioate-substituted oligonucleotide, an [alpha-(32)P]dNTP, and thermophilic enzymes. We also developed a high-resolution in-gel cross-linking assay to probe defined protein-DNA complexes. With these methods we have obtained structural information for the papillomavirus E1E2-ori preinitiation complex that would otherwise have been hard to obtain. These approaches should be widely applicable to the study of protein-DNA complexes.     

Combinatorial redesign of the DNA binding specificity of a prokaryotic helix-turn-helix repressor

Redesign of the bacteriophage 434 Cro repressor was accomplished by using an in vivo genetic screening system to identify new variants that specifically bound previously unrecognized DNA sequences. Site-directed, combinatorial mutagenesis of the 434 Cro helix-turn-helix (HTH) motif generated libraries of new variants which were screened for binding to new target sequences. Multiple mutations of 434 Cro that functionally converted wild-type (wt) 434 Cro DNA binding-sequence specificity to that of a lambda bacteriophage-specific repressor were identified. The libraries contained variations within the HTH sequence at only three positions. In vivo and in vitro analysis of several of the identified 434 Cro variants showed that the relatively few changes in the recognition helix of the HTH motif of 434 Cro resulted in specific and tight binding of the target DNA sequences. For the best 434 Cro variant identified, an apparent K(d) for lambda O(R)3 of 1 nM was observed. In competition experiments, this Cro variant was observed to be highly selective. We conclude that functional 434 Cro repressor variants with new DNA binding specificities can be generated from wt 434 Cro by mutating just the recognition helix. Important characteristics of the screening system responsible for the successful identifications are discussed. Application of the techniques presented here may allow the identification of DNA binding protein variants that functionally affect DNA regulatory sequences important in disease and industrial and biotechnological processes.     

The role of the minor groove substituents in indirect readout of DNA sequence by 434 repressor

The sequence of non-contacted bases at the center of the 434 repressor binding site affects the strength of the repressor-DNA complex by influencing the structure and flexibility of DNA (Koudelka, G. B., and Carlson, P. (1992) Nature 355, 89-91). We synthesized 434 repressor binding sites that differ in their central sequence base composition to test the importance of minor groove substituents and/or the number of base pair hydrogen bonds between these base pairs on DNA structure and strength of the repressor-DNA complex. We show here that the number of base pair H-bonds between the central bases apparently has no role in determining the relative affinity of a DNA site for repressor. Instead we find that the affinity of DNA for repressor depends on the absence or presence the N2-NH(2) group on the purine bases at the binding site center. The N2-NH(2) group on bases at the center of the 434 binding site appears to destabilize 434 repressor-DNA complexes by decreasing the intimacy of the specific repressor-DNA contacts, while increasing the reliance on protein contacts to the DNA phosphate backbone. Thus, the presence of an N2-NH(2) group on the purines at the center of a binding site globally alters the precise conformation of the protein-DNA interface.     

Intermolecular and intramolecular readout mechanisms in protein-DNA recognition

Protein-DNA recognition plays an essential role in the regulation of gene expression. Regulatory proteins are known to recognize specific DNA sequences directly through atomic contacts (intermolecular readout) and/or indirectly through the conformational properties of the DNA (intramolecular readout). However, little is known about the respective contributions made by these so-called direct and indirect readout mechanisms. We addressed this question by making use of information extracted from a structural database containing many protein-DNA complexes. We quantified the specificity of intermolecular (direct) readout by statistical analysis of base-amino acid interactions within protein-DNA complexes. The specificity of the intramolecular (indirect) readout due to DNA was quantified by statistical analysis of the sequence-dependent DNA conformation. Systematic comparison of these specificities in a large number of protein-DNA complexes revealed that both intermolecular and intramolecular readouts contribute to the specificity of protein-DNA recognition, and that their relative contributions vary depending upon the protein-DNA complexes. We demonstrated that combination of the intermolecular and intramolecular energies derived from the statistical analyses lead to enhanced specificity, and that the combined energy could explain experimental data on binding affinity changes caused by base mutations. These results provided new insight into the relationship between specificity and structure in the process of protein-DNA recognition, which would lead to prediction of specific protein-DNA binding sites.     

Quantitative parameters for amino acid-base interaction: implications for prediction of protein-DNA binding sites.

Inspection of the amino acid-base interactions in protein-DNA complexes is essential to the understanding of specific recognition of DNA target sites by regulatory proteins. The accumulation of information on protein-DNA co-crystals challenges the derivation of quantitative parameters for amino acid-base interaction based on these data. Here we use the coordinates of 53 solved protein-DNA complexes to extract all non-homologous pairs of amino acid-base that are in close contact, including hydrogen bonds and hydrophobic interactions. By comparing the frequency distribution of the different pairs to a theoretical distribution and calculating the log odds, a quantitative measure that expresses the likelihood of interaction for each pair of amino acid-base could be extracted. A score that reflects the compatibility between a protein and its DNA target can be calculated by summing up the individual measures of the pairs of amino acid-base involved in the complex, assuming additivity in their contributions to binding. This score enables ranking of different DNA binding sites given a protein binding site and vice versa and can be used in molecular design protocols. We demonstrate its validity by comparing the predictions using this score with experimental binding results of sequence variants of zif268 zinc fingers and their DNA binding sites.     

ParaDock: a flexible non-specific DNA--rigid protein docking algorithm

Accurate prediction of protein-DNA complexes could provide an important stepping stone towards a thorough comprehension of vital intracellular processes. Few attempts were made to tackle this issue, focusing on binding patch prediction, protein function classification and distance constraints-based docking. We introduce ParaDock: a novel ab initio protein-DNA docking algorithm. ParaDock combines short DNA fragments, which have been rigidly docked to the protein based on geometric complementarity, to create bent planar DNA molecules of arbitrary sequence. Our algorithm was tested on the bound and unbound targets of a protein-DNA benchmark comprised of 47 complexes. With neither addressing protein flexibility, nor applying any refinement procedure, CAPRI acceptable solutions were obtained among the 10 top ranked hypotheses in 83% of the bound complexes, and 70% of the unbound. Without requiring prior knowledge of DNA length and sequence, and within <2?h per target on a standard 2.0?GHz single processor CPU, ParaDock offers a fast ab initio docking solution.     

Expression, purification, and functional characterization of the carboxyl-terminal domain fragment of bacteriophage 434 repressor.

The repressor protein of bacteriophage 434 binds to DNA as a dimer of identical subunits. Its strong dimerization is mediated by the carboxyl-terminal domain. Cooperative interactions between the C-terminal domains of two repressor dimers bound at adjacent sites can stabilize protein-DNA complexes formed with low-affinity binding sites. We have constructed a plasmid, pCT1, which directs the overproduction of the carboxyl-terminal domain of 434 repressor. The protein encoded by this plasmid is called CT-1. Cells transformed with pCT1 are unable to be lysogenized by wild-type 434 phage, whereas control cells are lysogenized at an efficiency of 1 to 5%. The CT-1-mediated interference with lysogen formation presumably results from formation of heteromeric complexes between the phage-encoded repressor and the plasmid-encoded carboxyl-terminal domain fragment. These heteromers are unable to bind DNA and thereby inhibit the repressor's activity in promoting lysogen formation. Two lines of evidence support this conclusion. First, DNase I footprinting experiments show that at a 2:1 ratio of CT-1 to intact 434 repressor, purified CT-1 protein prevents the formation of complexes between 434 repressor and its OR1 binding site. Second, cross-linking experiments reveal that only a specific heterodimeric complex forms between CT-1 and intact 434 repressor. This latter observation indicates that CT-1 interferes with 434 repressor-operator complex formation by preventing dimerization and not by altering the conformation of the DNA-bound repressor dimer. Our other evidence is also consistent with this suggestion. We have used deletion analysis in an attempt to define the region which mediates the 434 repressor-CT-1 interaction. CT-1 proteins which have more than the last 14 amino acids removed are unable to interfere with 434 repressor action in vivo.     

Structure-based prediction of DNA target sites by regulatory proteins

Regulatory proteins play a critical role in controlling complex spatial and temporal patterns of gene expression in higher organism, by recognizing multiple DNA sequences and regulating multiple target genes. Increasing amounts of structural data on the protein-DNA complex provides clues for the mechanism of target recognition by regulatory proteins. The analyses of the propensities of base-amino acid interactions observed in those structural data show that there is no one-to-one correspondence in the interaction, but clear preferences exist. On the other hand, the analysis of spatial distribution of amino acids around bases shows that even those amino acids with strong base preference such as Arg with G are distributed in a wide space around bases. Thus, amino acids with many different geometries can form a similar type of interaction with bases. The redundancy and structural flexibility in the interaction suggest that there are no simple rules in the sequence recognition, and its prediction is not straightforward. However, the spatial distributions of amino acids around bases indicate a possibility that the structural data can be used to derive empirical interaction potentials between amino acids and bases. Such information extracted from structural databases has been successfully used to predict amino acid sequences that fold into particular protein structures. We surmised that the structures of protein-DNA complexes could be used to predict DNA target sites for regulatory proteins, because determining DNA sequences that bind to a particular protein structure should be similar to finding amino acid sequences that fold into a particular structure. Here we demonstrate that the structural data can be used to predict DNA target sequences for regulatory proteins. Pairwise potentials that determine the interaction between bases and amino acids were empirically derived from the structural data. These potentials were then used to examine the compatibility between DNA sequences and the protein-DNA complex structure in a combinatorial "threading" procedure. We applied this strategy to the structures of protein-DNA complexes to predict DNA binding sites recognized by regulatory proteins. To test the applicability of this method in target-site prediction, we examined the effects of cognate and noncognate binding, cooperative binding, and DNA deformation on the binding specificity, and predicted binding sites in real promoters and compared with experimental data. These results show that target binding sites for several regulatory proteins are successfully predicted, and our data suggest that this method can serve as a powerful tool for predicting multiple target sites and target genes for regulatory proteins.     

Nucleotide sequence of the cro-cII-oop region of bacteriophage 434 DNA.

The nucleotide sequence of a 869 bp segment of phage 434 DNA including the regulatory genes cro and cII is presented and compared with the corresponding part of the phage lambda DNA sequence. The 434 cro protein as deduced from the DNA sequence is a highly basic protein of 71 amino acid residues with a calculated molecular weight of 8089. While the cro gene sequences of phage 434 and lambda DNA are very different, the nuleotide sequences to the right of the lambda imm434 boundary show differences only at 11 out of 512 positions. Nucleotide substitutions in the cII gene occur with one exception in the third positions of the respective codons and only one out of several DNA regulatory signals located in this region of the phage genomes is affected by these nucleotide substitutions.     

Novel method for identifying sequence-specific DNA-binding proteins.

We developed a general method for the enrichment and identification of sequence-specific DNA-binding proteins. A well-characterized protein-DNA interaction is used to isolate from crude cellular extracts or fractions thereof proteins which bind to specific DNA sequences; the method is based solely on this binding property of the proteins. The DNA sequence of interest, cloned adjacent to the lac operator DNA segment is incubated with a lac repressor-beta-galactosidase fusion protein which retains full operator and inducer binding properties. The DNA fragment bound to the lac repressor-beta-galactosidase fusion protein is precipitated by the addition of affinity-purified anti-beta-galactosidase immobilized on beads. This forms an affinity matrix for any proteins which might interact specifically with the DNA sequence cloned adjacent to the lac operator. When incubated with cellular extracts in the presence of excess competitor DNA, any protein(s) which specifically binds to the cloned DNA sequence of interest can be cleanly precipitated. When isopropyl-beta-D-thiogalactopyranoside is added, the lac repressor releases the bound DNA, and thus the protein-DNA complex consisting of the specific restriction fragment and any specific binding protein(s) is released, permitting the identification of the protein by standard biochemical techniques. We demonstrate the utility of this method with the lambda repressor, another well-characterized DNA-binding protein, as a model. In addition, with crude preparations of the yeast mitochondrial RNA polymerase, we identified a 70,000-molecular-weight peptide which binds specifically to the promoter region of the yeast mitochondrial 14S rRNA gene.