A phage display selection of engrailed homeodomain mutants and the importance of residue Q50

Diversifying the structure of zinc finger nucleases for high-precision genome editing

Genome editing for therapeutic applications often requires cleavage within a narrow sequence window. Here, to enable such high-precision targeting with zinc-finger nucleases (ZFNs), we have developed an expanded set of architectures that collectively increase the configurational options available for design by a factor of 64. These new architectures feature the functional attachment of the FokI cleavage domain to the amino terminus of one or both zinc-finger proteins (ZFPs) in the ZFN dimer, as well as the option to skip bases between the target triplets of otherwise adjacent fingers in each zinc-finger array. Using our new architectures, we demonstrate targeting of an arbitrarily chosen 28 bp genomic locus at a density that approaches 1.0 (i.e., efficient ZFNs available for targeting almost every base step). We show that these new architectures may be used for targeting three loci of therapeutic significance with a high degree of precision, efficiency, and specificity.

Reprogramming cells with synthetic proteins

Conversion of one cell type into another cell type by forcibly expressing specific cocktails of transcription factors (TFs) has demonstrated that cell fates are not fixed and that cellular differentiation can be a two-way street with many intersections. These experiments also illustrated the sweeping potential of TFs to “read” genetically hardwired regulatory information even in cells where they are not normally expressed and to access and open up tightly packed chromatin to execute gene expression programs. Cellular reprogramming enables the modeling of diseases in a dish, to test the efficacy and toxicity of drugs in patient-derived cells and ultimately, could enable cell-based therapies to cure degenerative diseases. Yet, producing terminally differentiated cells that fully resemble their in vivo counterparts in sufficient quantities is still an unmet clinical need. While efforts are being made to reprogram cells nongenetically by using drug-like molecules, defined TF cocktails still dominate reprogramming protocols. Therefore, the optimization of TFs by protein engineering has emerged as a strategy to enhance reprogramming to produce functional, stable and safe cells for regenerative biomedicine. Engineering approaches focused on Oct4, MyoD, Sox17, Nanog and Mef2c and range from chimeric TFs with added transactivation domains, designer transcription activator-like effectors to activate endogenous TFs to reprogramming TFs with rationally engineered DNA recognition principles. Possibly, applying the complete toolkit of protein design to cellular reprogramming can help to remove the hurdles that, thus far, impeded the clinical use of cells derived from reprogramming technologies.

Selection of heptapeptides that bind helix 69 of bacterial 23S ribosomal RNA

Helix 69 of Escherichia coli 23S rRNA has important roles in specific steps of translation, such as subunit association, translocation, and ribosome recycling. An M13 phage library was used to identify peptide ligands with affinity for helix 69. One selected sequence, NQVANHQ, was shown through a bead assay to interact with helix 69. Electrospray ionization mass spectroscopy revealed an apparent dissociation constant for the amidated peptide and helix 69 in the low micromolar range. This value is comparable to that of aminoglycoside antibiotics binding to the A site of 16S rRNA or helix 69. Helix 69 variants (human) and unrelated RNAs (helix 31 or A site of 16S rRNA) showed two- to fourfold lower affinity for NQVANHQ-NH(2). These results suggest that the peptide has desirable features for development as a lead compound for novel antimicrobials.

Probing the informational and regulatory plasticity of a transcription factor DNA-binding domain

Transcription factors have two functional constraints on their evolution: (1) their binding sites must have enough information to be distinguishable from all other sequences in the genome, and (2) they must bind these sites with an affinity that appropriately modulates the rate of transcription. Since both are determined by the biophysical properties of the DNA–binding domain, selection on one will ultimately affect the other. We were interested in understanding how plastic the informational and regulatory properties of a transcription factor are and how transcription factors evolve to balance these constraints. To study this, we developed an in vivo selection system in Escherichia coli to identify variants of the helix-turn-helix transcription factor MarA that bind different sets of binding sites with varying degrees of degeneracy. Unlike previous in vitro methods used to identify novel DNA binders and to probe the plasticity of the binding domain, our selections were done within the context of the initiation complex, selecting for both specific binding within the genome and for a physiologically significant strength of interaction to maintain function of the factor. Using MITOMI, quantitative PCR, and a binding site fitness assay, we characterized the binding, function, and fitness of some of these variants. We observed that a large range of binding preferences, information contents, and activities could be accessed with a few mutations, suggesting that transcriptional regulatory networks are highly adaptable and expandable.

The main role of transcription factors is to modulate the expression levels of functionally related genes in response to environmental and cellular cues. For this process to be precise, the transcription factor needs to locate and bind specific DNA sequences in the genome and needs to bind these sites with a strength that appropriately adjusts the amount of gene expressed. Both specific protein–DNA interactions and transcription factor activity are intimately coupled, because they are both dependent upon the biochemical properties of the DNA–binding domain. Here we experimentally probe how variable these properties are using a novel in vivo selection assay. We observed that the specific binding preferences for the transcription factor MarA and its transcriptional activity can be altered over a large range with a few mutations and that selection on one function will impact the other. This work helps us to better understand the mechanism of transcriptional regulation and its evolution, and may prove useful for the engineering of transcription factors and regulatory networks.

The Homeodomain Resource: a comprehensive collection of sequence, structure, interaction, genomic and functional information on the homeodomain protein family

The Homeodomain Resource is a curated collection of sequence, structure, interaction, genomic and functional information on the homeodomain family. The current version builds upon previous versions by the addition of new, complete sets of homeodomain sequences from fully sequenced genomes, the expansion of existing curated homeodomain information and the improvement of data accessibility through better search tools and more complete data integration. This release contains 1534 full-length homeodomain-containing sequences, 93 experimentally derived homeodomain structures, 101 homeodomain protein-protein interactions, 107 homeodomain DNA-binding sites and 206 homeodomain proteins implicated in human genetic disorders.Database URL: The Homeodomain Resource is freely available and can be accessed at http://research.nhgri.nih.gov/homeodomain/

A one-dimensional reaction coordinate for identification of transition states from explicit solvent P(fold)-like calculations

A properly identified transition state ensemble (TSE) in a molecular dynamics (MD) simulation can reveal a tremendous amount about how a protein folds and offer a point of comparison to experimentally derived Phi(F) values, which reflect the degree of structure in these transient states. In one such method of TSE identification, dubbed P(fold), MD simulations of individual protein structures taken from an unfolding trajectory are used to directly assess an input structure's probability of folding before unfolding, and P(fold) is, by definition, 0.5 for the TSE. Other, less computationally intensive methods, such as multidimensional scaling (MDS) of the pairwise root mean-squared deviation (RMSD) matrix of the conformations sampled in a thermal unfolding trajectory, have also been used to identify the TSE. Identification of the TSE is made from the original MD simulation without the need to run further simulations. Here we present a P(fold)-like study and describe methods for identification of the TSE through the derivation of a high fidelity, bounded, one-dimensional reaction coordinate for protein folding. These methods are applied to the engrailed homeodomain. The TSE identified by this approach is essentially identical to the TSE identified previously by MDS of the pairwise RMSD matrix. However, the cost of performing P(fold), or even our reduced P(fold)-like calculations, is at least 36,000 times greater than the MDS method.

Profiling the DNA-binding specificities of engineered Cys2His2 zinc finger domains using a rapid cell-based method

The C2H2 zinc finger is the most commonly utilized framework for engineering DNA-binding domains with novel specificities. Many different selection strategies have been developed to identify individual fingers that possess a particular DNA-binding specificity from a randomized library. In these experiments, each finger is selected in the context of a constant finger framework that ensures the identification of clones with a desired specificity by properly positioning the randomized finger on the DNA template. Following a successful selection, multiple zinc-finger clones are typically recovered that share similarities in the sequences of their DNA-recognition helices. In principle, each of the clones isolated from a selection is a candidate for assembly into a larger multi-finger protein, but to date a high-throughput method for identifying the most specific candidates for incorporation into a final multi-finger protein has not been available. Here we describe the development of a specificity profiling system that facilitates rapid and inexpensive characterization of engineered zinc-finger modules. Moreover, we demonstrate that specificity data collected using this system can be employed to rationally design zinc fingers with improved DNA-binding specificities.

Sequence specificity of single-stranded DNA-binding proteins: a novel DNA microarray approach

We have developed a novel DNA microarray-based approach for identification of the sequence-specificity of single-stranded nucleic-acid-binding proteins (SNABPs). For verification, we have shown that the major cold shock protein (CspB) from Bacillus subtilis binds with high affinity to pyrimidine-rich sequences, with a binding preference for the consensus sequence, 5′-GTCTTTG/T-3′. The sequence was modelled onto the known structure of CspB and a cytosine-binding pocket was identified, which explains the strong preference for a cytosine base at position 3. This microarray method offers a rapid high-throughput approach for determining the specificity and strength of ss DNA–protein interactions. Further screening of this newly emerging family of transcription factors will help provide an insight into their cellular function.

Structure and properties of a re-engineered homeodomain protein-DNA interface

The homeodomain (HD)-DNA interface has been conserved over 500 million years of evolution. Despite this conservation, we have successfully re-engineered the engrailed HD to specifically recognize an unnatural nucleotide using a phage display selection. Here we report the synthesis of novel nucleosides and the selection of mutant HDs that bind these nucleotides using phage display. The high-resolution crystal structure of one mutant in complex with modified and unmodified DNA demonstrates that, even with the substantial perturbation to the interface, this selected mutant retains a canonical HD structure. Dissection of the contributions due to each of the selected mutations reveals that the majority of the modification-specific binding is accomplished by a single mutation (I47G) but that the remaining mutations retune the stability of the HD. These results afford a detailed look at a re-engineered protein-DNA interaction and provide insight into the opportunities for re-engineering highly conserved interfaces.

Dissecting the Engrailed homeodomain-DNA interaction by phage-displayed shotgun scanning

Phage-displayed alanine shotgun scanning was used to dissect contributions by engrailed homedomain (En-HD) residues 17 through 46, which indirectly influence recognition of DNA. The relative contributions of such indirect contacts, quantified by shotgun scanning, highlight previously unexplored En-HD residues. Two motifs dominate En-HD function in this region. First, two surface-exposed aromatic residues (F20 and Y25) bracket the hydrophobic core. Second, two sets of turn-forming residues are highlighted, including carboxamide-requiring residues E22/N23 and a leucine/isoleucine splint. The En-HD hydrophobic core exhibits a surprising degree of malleability, as demonstrated by homolog shotgun scanning. Most selectants from in vitro shotgun scanning mirror the consensus human homeodomain sequence. Thus, natural evolution and in vitro selection use similar selection criteria: affinity, specificity, and stability. However, homolog shotgun scanning identified mutations capable of improving the affinity and specificity of En-HD.

Homeodomain revisited: a lesson from disease-causing mutations

The homeodomain is a highly conserved DNA-binding motif that is found in numerous transcription factors throughout a large variety of species from yeast to humans. These gene-specific transcription factors play critical roles in development and adult homeostasis, and therefore, any germline mutations associated with these proteins can lead to a number of congenital abnormalities. Although much has been revealed concerning the molecular architecture and the mechanism of homeodomain-DNA interactions, the study of disease-causing mutations can further provide us with instructive information as to the role of particular residues in a conserved mode of action. In this paper, I have compiled the homeodomain missense mutations found in various human diseases and re-examined the functional role of the mutational "hot spot" residues in light of the structures obtained from crystallography. These findings should be useful in understanding the essential components of the homeodomain and in attempts to design agonist or antagonists to modulate their activity and to reverse the effects caused by the mutations.