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.     

Single exchanges of amino acids in the basic region change the specificity of N-Myc.

We exchanged specific amino acids in the basic region of the murine N-Myc protein and tested the mutant proteins for their DNA binding specificity. The amino acids we exchanged were chosen in analogy to residues of the homologous basic regions of bHLH and bZIP proteins. Mutant N-Myc peptides were expressed in Escherichia coli and specific DNA binding was monitored by gel shift experiments. For this we used palindromic target sequences with systematic base pair exchanges. Several mutants with altered DNA binding specificity were identified. Amino acid exchanges of residues -14 or -10 of the basic region lead to specificity changes (we define leucine 402 of N-Myc as +1; comparable to GCN4 see (1)). The palindromic N-Myc recognition sequence 5'CACGTG is no longer recognized by the mutant proteins, but DNA fragments with symmetrical exchanges of the target sequence are. Exchanges at position -15 broaden the binding specificity. These data were used to build a computer based model of the putative interactions of the N-Myc basic DNA binding region with its target sequence.     

Tailored amino acid diversity for the evolution of antibody affinity

Antibodies are a unique class of proteins with the ability to adapt their binding sites for high affinity and high specificity to a multitude of antigens. Many analyses have been performed on antibody sequences and structures to elucidate which amino acids have a predominant role in antibody interactions with antigens. These studies have generally not distinguished between amino acids selected for broad antigen specificity in the primary immune response and those selected for high affinity in the secondary immune response. By studying a large data set of affinity matured antibodies derived from in vitro directed evolution experiments, we were able to specifically highlight a subset of amino acids associated with affinity improvements. In a comparison of affinity maturations using either tailored or full amino acid diversification, the tailored approach was found to be at least as effective at improving affinity while requiring fewer mutagenesis libraries than the traditional method. The resulting sequence data also highlight the potential for further reducing amino acid diversity for high affinity binding interactions.     

Tailored amino acid diversity for the evolution of antibody affinity

Antibodies are a unique class of proteins with the ability to adapt their binding sites for high affinity and high specificity to a multitude of antigens. Many analyses have been performed on antibody sequences and structures to elucidate which amino acids have a predominant role in antibody interactions with antigens. These studies have generally not distinguished between amino acids selected for broad antigen specificity in the primary immune response and those selected for high affinity in the secondary immune response. By studying a large data set of affinity matured antibodies derived from in vitro directed evolution experiments, we were able to specifically highlight a subset of amino acids associated with affinity improvements. In a comparison of affinity maturations using either tailored or full amino acid diversification, the tailored approach was found to be at least as effective at improving affinity while requiring fewer mutagenesis libraries than the traditional method. The resulting sequence data also highlight the potential for further reducing amino acid diversity for high affinity binding interactions.     

A simple physical model for the prediction and design of protein-DNA interactions

Protein-DNA interactions are crucial for many biological processes. Attempts to model these interactions have generally taken the form of amino acid-base recognition codes or purely sequence-based profile methods, which depend on the availability of extensive sequence and structural information for specific structural families, neglect side-chain conformational variability, and lack generality beyond the structural family used to train the model. Here, we take advantage of recent advances in rotamer-based protein design and the large number of structurally characterized protein-DNA complexes to develop and parameterize a simple physical model for protein-DNA interactions. The model shows considerable promise for redesigning amino acids at protein-DNA interfaces, as design calculations recover the amino acid residue identities and conformations at these interfaces with accuracies comparable to sequence recovery in globular proteins. The model shows promise also for predicting DNA-binding specificity for fixed protein sequences: native DNA sequences are selected correctly from pools of competing DNA substrates; however, incorporation of backbone movement will likely be required to improve performance in homology modeling applications. Interestingly, optimization of zinc finger protein amino acid sequences for high-affinity binding to specific DNA sequences results in proteins with little or no predicted specificity, suggesting that naturally occurring DNA-binding proteins are optimized for specificity rather than affinity. When combined with algorithms that optimize specificity directly, the simple computational model developed here should be useful for the engineering of proteins with novel DNA-binding specificities.     

Sequence-specific DNA-protein interaction: the lac represser

A model is proposed for lac repressor-lac operator binding which accounts for the tetrameric subunit structure of the lac repressor and for factors involved in the strength, specificity and regulation of repressor-operator interaction. The model employs a π-helix in the amino terminal 25 residues of the lac repressor whereby three tyrosine residues of each subunit intercalate between base pairs of the lac operator. For the outer palindromic sequences of the operator, base specificity is provided by amino acids adjacent to the carboxyl sides of the tyrosine residues of two of the subunits. The inner palindromic sequences which bind the other two subunits form stems of hairpin loops in the operator. Base specificity for these two subunits is provided by amino acids adjacent to the amino sides of the tyrosine residues. In addition to 12 intercalated tyrosine residues, the model provides for a total of at least eight electrostatic interactions and ten sequence-specific hydrogen bonds.     

Stacking interactions in PUF-RNA complexes

Stacking interactions between amino acids and bases are common in RNA-protein interactions. Many proteins that regulate mRNAs interact with single-stranded RNA elements in the 3' UTR (3'-untranslated region) of their targets. PUF proteins are exemplary. Here we focus on complexes formed between a Caenorhabditis elegans PUF protein, FBF, and its cognate RNAs. Stacking interactions are particularly prominent and involve every RNA base in the recognition element. To assess the contribution of stacking interactions to formation of the RNA-protein complex, we combine in vivo selection experiments with site-directed mutagenesis, biochemistry, and structural analysis. Our results reveal that the identities of stacking amino acids in FBF affect both the affinity and specificity of the RNA-protein interaction. Substitutions in amino acid side chains can restrict or broaden RNA specificity. We conclude that the identities of stacking residues are important in achieving the natural specificities of PUF proteins. Similarly, in PUF proteins engineered to bind new RNA sequences, the identity of stacking residues may contribute to "target" versus "off-target" interactions, and thus be an important consideration in the design of proteins with new specificities.     

A single amino acid can determine the DNA binding specificity of homeodomain proteins

Many Drosophila developmental genes contain a DNA binding domain encoded by the homeobox. This homeodomain contains a region distantly homologous to the helix-turn-helix motif present in several prokaryotic DNA binding proteins. We investigated the nature of homeodomain-DNA interactions by making a series of mutations in the helix-turn-helix motif of the Drosophila homeodomain protein Paired (Prd). This protein does not recognize sequences bound by the homeodomain proteins Fushi tarazu (Ftz) or Bicoid (Bcd). We show that changing a single amino acid at the C-terminus of the recognition helix is both necessary and sufficient to confer the DNA binding specificity of either Ftz or Bcd on Prd. This simple rule indicates that the amino acids that determine the specificity of homeodomains are different from those mediating protein-DNA contacts in prokaryotic proteins. We further show that Prd contains two DNA binding activities. The Prd homeodomain is responsible for one of them while the other is not dependent on the recognition helix.     

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.     

Protein-DNA interactions: amino acid conservation and the effects of mutations on binding specificity

We investigate the conservation of amino acid residue sequences in 21 DNA-binding protein families and study the effects that mutations have on DNA-sequence recognition. The observations are best understood by assigning each protein family to one of three classes: (i) non-specific, where binding is independent of DNA sequence; (ii) highly specific, where binding is specific and all members of the family target the same DNA sequence; and (iii) multi-specific, where binding is also specific, but individual family members target different DNA sequences. Overall, protein residues in contact with the DNA are better conserved than the rest of the protein surface, but there is a complex underlying trend of conservation for individual residue positions. Amino acid residues that interact with the DNA backbone are well conserved across all protein families and provide a core of stabilising contacts for homologous protein-DNA complexes. In contrast, amino acid residues that interact with DNA bases have variable levels of conservation depending on the family classification. In non-specific families, base-contacting residues are well conserved and interactions are always found in the minor groove where there is little discrimination between base types. In highly specific families, base-contacting residues are highly conserved and allow member proteins to recognise the same target sequence. In multi-specific families, base-contacting residues undergo frequent mutations and enable different proteins to recognise distinct target sequences. Finally, we report that interactions with bases in the target sequence often follow (though not always) a universal code of amino acid-base recognition and the effects of amino acid mutations can be most easily understood for these interactions.     

Repertoire selection of variant single-chain Cro: toward directed DNA-binding specificity of helix-turn-helix proteins

A single-chain derivative of the lambda Cro repressor (scCro) has been randomly mutated in amino acid residues critical for specific DNA recognition to create libraries of protein variants. Utilizing phage display-afforded affinity selection, scCro variants have been isolated for binding to synthetic DNA ligands. Isolated scCro variants were analyzed functionally, both in fusion with phage particles and after expression of the corresponding free proteins. The binding properties with regard to specificity and affinity in binding to different DNA ligands were investigated by inhibition studies and determination of equilibrium dissociation constants for formed complexes. Variant proteins with altered DNA-sequence specificity were identified, which favored binding of targeted synthetic DNA sequences over a consensus operator sequence, bound with high affinity by wild-type Cro. The specificities were relatively modest (2-3-fold, as calculated from K(D) values), which can be attributed to the inherent properties in the design of the selection system; one half-site of the synthetic DNA sequences maintains the consensus operator sequence, and one "subunit" of the variant single-chain Cro dimers was conserved as wild-type sequence. The anticipated interaction between the wild-type subunit and the consensus DNA half-site of target DNA ligands is, hence, expected to contribute to the overlap in sequence discrimination. The binding affinity for the synthetic DNA sequences, however, was improved 10-30-fold in selected variant proteins as compared to "wild-type" scCro.     

Rational design of DNA sequence-specific zinc fingers

We developed a rational scheme for designing DNA binding proteins. The scheme was applied for a zinc finger protein and the designed sequences were experimentally characterized with high DNA sequence specificity. Starting with the backbone of a known finger structure, we initially calculated amino acid sequences compatible with the expected structure and the secondary structures of the designed fingers were then experimentally confirmed. The DNA-binding function was added to the designed finger by reconsidering a section of the amino acid sequence and computationally selecting amino acids to have the lowest protein-DNA interaction energy for the target DNA sequences. Among the designed proteins, one had a gap between the lowest and second lowest protein-DNA interaction energies that was sufficient to give DNA sequence-specificity.