Cells that register logical relationships among proteins

Two-hybrid methods have augmented the classical genetic techniques biologists use to assign function to genes. Here, we describe construction of a two-bait interaction trap that uses yeast cells to register more complex protein relationships than those detected in existing two-hybrid systems. We show that such cells can identify bridge or connecting proteins and peptide aptamers that discriminate between closely related allelic variants. The protein relationships detected by these cells are analogous to classical genetic relationships, but lend themselves to systematic application to the products of entire genomes and combinatorial libraries. We show that, by performing logical operations on the phenotypic outputs of these complex cells and existing two-hybrid cells, we can make inferences about the topology and order of protein interactions. Finally, we show that cells that register such relationships can perform logical operations on protein inputs. Thus these cells will be useful for analysis of gene and allele function, and may also define a path for construction of biological computational devices.

Deciphering the Nuclear Import Pathway for the Cytoskeletal Red Cell Protein 4.1R

The erythroid membrane cytoskeletal protein 4.1 is the prototypical member of a genetically and topologically complex family that is generated by combinatorial alternative splicing pathways and is localized at diverse intracellular sites including the nucleus. To explore the molecular determinants for nuclear localization, we transfected COS-7 cells with epitope-tagged versions of natural red cell protein 4.1 (4.1R) isoforms as well as mutagenized and truncated derivatives. Two distant topological sorting signals were required for efficient nuclear import of the 4.1R80 isoform: a basic peptide, KKKRER, encoded by alternative exon 16 and acting as a weak core nuclear localization signal (4.1R NLS), and an acidic peptide, EED, encoded by alternative exon 5. 4.1R80 isoforms lacking either of these two exons showed decreased nuclear import. Fusion of various 4.1R80 constructs to the cytoplasmic reporter protein pyruvate kinase confirmed a requirement for both motifs for full NLS function. 4.1R80 was efficiently imported in the nuclei of digitonin-permeabilized COS-7 cells in the presence of recombinant Rch1 (human importin α2), importin β, and GTPase Ran. Quantitative analysis of protein–protein interactions using a resonant mirror detection technique showed that 4.1R80 bound to Rch1 in vitro with high affinity (KD = 30 nM). The affinity decreased at least 7- and 20-fold, respectively, if the EED motif in exon 5 or if 4.1R NLS in exon 16 was lacking or mutated, confirming that both motifs were required for efficient importin-mediated nuclear import of 4.1R80.

Functional and Biochemical Analysis of the C2 Domains of Synaptotagmin IV

Synaptotagmins (Syts) are a family of vesicle proteins that have been implicated in both regulated neurosecretion and general membrane trafficking. Calcium-dependent interactions mediated through their C2 domains are proposed to contribute to the mechanism by which Syts trigger calcium-dependent neurotransmitter release. Syt IV is a novel member of the Syt family that is induced by cell depolarization and has a rapid rate of synthesis and a short half-life. Moreover, the C2A domain of Syt IV does not bind calcium. We have examined the biochemical and functional properties of the C2 domains of Syt IV. Consistent with its non–calcium binding properties, the C2A domain of Syt IV binds syntaxin isoforms in a calcium-independent manner. In neuroendocrine pheochromocytoma (PC12) cells, Syt IV colocalizes with Syt I in the tips of the neurites. Microinjection of the C2A domain reveals that calcium-independent interactions mediated through this domain of Syt IV inhibit calcium-mediated neurotransmitter release from PC12 cells. Conversely, the C2B domain of Syt IV contains calcium binding properties, which permit homo-oligomerization as well as hetero-oligomerization with Syt I. Our observation that different combinatorial interactions exist between Syt and syntaxin isoforms, coupled with the calcium stimulated hetero-oligomerization of Syt isoforms, suggests that the secretory machinery contains a vast repertoire of biochemical properties for sensing calcium and regulating neurotransmitter release accordingly.

Molecular computation: RNA solutions to chess problems

We have expanded the field of “DNA computers” to RNA and present a general approach for the solution of satisfiability problems. As an example, we consider a variant of the “Knight problem,” which asks generally what configurations of knights can one place on an n × n chess board such that no knight is attacking any other knight on the board. Using specific ribonuclease digestion to manipulate strands of a 10-bit binary RNA library, we developed a molecular algorithm and applied it to a 3 × 3 chessboard as a 9-bit instance of this problem. Here, the nine spaces on the board correspond to nine “bits” or placeholders in a combinatorial RNA library. We recovered a set of “winning” molecules that describe solutions to this problem.

Multiple genetic modifications of the erythromycin polyketide synthase to produce a library of novel “unnatural” natural products

The structures of complex polyketide natural products, such as erythromycin, are programmed by multifunctional polyketide synthases (PKSs) that contain modular arrangements of functional domains. The colinearity between the activities of modular PKS domains and structure of the polyketide product portends the generation of novel organic compounds—“unnatural” natural products—by genetic manipulation. We have engineered the erythromycin polyketide synthase genes to effect combinatorial alterations of catalytic activities in the biosynthetic pathway, generating a library of >50 macrolides that would be impractical to produce by chemical methods. The library includes examples of analogs with one, two, and three altered carbon centers of the polyketide products. The manipulation of multiple biosynthetic steps in a PKS is an important milestone toward the goal of producing large libraries of unnatural natural products for biological and pharmaceutical applications.

Small antibody-like proteins with prescribed ligand specificities derived from the lipocalin fold

We demonstrate that the ligand pocket of a lipocalin from Pieris brassicae, the bilin-binding protein (BBP), can be reshaped by combinatorial protein design such that it recognizes fluorescein, an established immunological hapten. For this purpose 16 residues at the center of the binding site, which is formed by four loops on top of an eight-stranded β-barrel, were subjected to random mutagenesis. Fluorescein-binding BBP variants were then selected from the mutant library by bacterial phage display. Three variants were identified that complex fluorescein with high affinity, exhibiting dissociation constants as low as 35.2 nM. Notably, one of these variants effects almost complete quenching of the ligand fluorescence, similarly as an anti-fluorescein antibody. Detailed ligand-binding studies and site-directed mutagenesis experiments indicated (i) that the molecular recognition of fluorescein is specific and (ii) that charged residues at the center of the pocket are responsible for tight complex formation. Sequence comparison of the BBP variants directed against fluorescein with the wild-type protein and with further variants that were selected against several other ligands revealed that all of the randomized amino acid positions are variable. Hence, a lipocalin can be used for generating molecular pockets with a diversity of shapes. We term this class of engineered proteins “anticalins.” Their one-domain scaffold makes them a promising alternative to antibodies to create a stable receptor protein for a ligand of choice.

Parent–progeny recognition as a function of MHC odortype identity

The several linked polymorphic genes of the MHC, which has been proposed as a prime determinant of sensed genetic individuality within species, is known to operate in mice by olfactory recognition in aspects of reproductive behavior that concern mate selection, thereby favoring outbreeding and heterozygosity, and also concern the maintenance of pregnancy. A single base-change can alter an individual MHC odortype, and the potential range of combinatorial MHC-determined odortypes is clearly vast. Following our findings that newborn mice already express their MHC odortype (which is detectable at 9 days of gestational age), we sought to determine whether MHC is involved in behavioral aspects of early development, such as rearing. In the studies presented herein, we report the ability and proclivity of mothers to recognize and preferentially retrieve syngeneic (genetically identical) pups from other pups differing only for MHC. Reciprocally, we report the ability of pups to recognize their familial environment, regardless of whether they had been nursed by their biological mothers or by foster mothers. Early learning experiences of the MHC environment are apparently a key element in survival, assuring maternal protection and promoting outbreeding.

Stepwise in vitro affinity maturation of Vitaxin, an αvβ3-specific humanized mAb

A protein engineering strategy based on efficient and focused mutagenesis implemented by codon-based mutagenesis was developed. Vitaxin, a humanized version of the antiangiogenic antibody LM609 directed against a conformational epitope of the αvβ3 integrin complex, was used as a model system. Specifically, focused mutagenesis was used in a stepwise fashion to rapidly improve the affinity of the antigen binding fragment by greater than 90-fold. In the complete absence of structural information about the Vitaxin-αvβ3 interaction, phage-expressed antibody libraries for all six Ig heavy and light chain complementarity-determining regions were expressed and screened by a quantitative assay to identify variants with improved binding to αvβ3. The Vitaxin variants in these libraries each contained a single mutation, and all 20 amino acids were introduced at each complementarity-determining region residue, resulting in the expression of 2,336 unique clones. Multiple clones displaying 2- to 13-fold improved affinity were identified. Subsequent expression and screening of a library of 256 combinatorial variants of the optimal mutations identified from the primary libraries resulted in the identification of multiple clones displaying greater than 50-fold enhanced affinity. These variants inhibited ligand binding to receptor more potently as demonstrated by inhibition of cell adhesion and ligand competition assays. Because of the limited mutagenesis and combinatorial approach, Vitaxin variants with enhanced affinity were identified rapidly and required the synthesis of only 2,592 unique variants. The use of such small focused libraries obviates the need for phage affinity selection approaches typically used, permitting the use of functional assays and the engineering of proteins expressed in mammalian cell culture.

Unbinding forces of single antibody-antigen complexes correlate with their thermal dissociation rates

Point mutants of three unrelated antifluorescein antibodies were constructed to obtain nine different single-chain Fv fragments, whose on-rates, off-rates, and equilibrium binding affinities were determined in solution. Additionally, activation energies for unbinding were estimated from the temperature dependence of the off-rate in solution. Loading rate-dependent unbinding forces were determined for single molecules by atomic force microscopy, which extrapolated at zero force to a value close to the off-rate measured in solution, without any indication for multiple transition states. The measured unbinding forces of all nine mutants correlated well with the off-rate in solution, but not with the temperature dependence of the reaction, indicating that the same transition state must be crossed in spontaneous and forced unbinding and that the unbinding path under load cannot be too different from the one at zero force. The distance of the transition state from the ground state along the unbinding pathway is directly proportional to the barrier height, regardless of the details of the binding site, which most likely reflects the elasticity of the protein in the unbinding process. Atomic force microscopy thus can be a valuable tool for the characterization of solution properties of protein-ligand systems at the single molecule level, predicting relative off-rates, potentially of great value for combinatorial chemistry and biology.

Rapid generation of incremental truncation libraries for protein engineering using α-phosphothioate nucleotides

Incremental truncation for the creation of hybrid enzymes (ITCHY) is a novel tool for the generation of combinatorial libraries of hybrid proteins independent of DNA sequence homology. We herein report a fundamentally different methodology for creating incremental truncation libraries using nucleotide triphosphate analogs. Central to the method is the polymerase catalyzed, low frequency, random incorporation of α-phosphothioate dNTPs into the region of DNA targeted for truncation. The resulting phosphothioate internucleotide linkages are resistant to 3′→5′ exonuclease hydrolysis, rendering the target DNA resistant to degradation in a subsequent exonuclease III treatment. From an experimental perspective the protocol reported here to create incremental truncation libraries is simpler and less time consuming than previous approaches by combining the two gene fragments in a single vector and eliminating additional purification steps. As proof of principle, an incremental truncation library of fusions between the N-terminal fragment of Escherichia coli glycinamide ribonucleotide formyltransferase (PurN) and the C-terminal fragment of human glycinamide ribonucleotide formyltransferase (hGART) was prepared and successfully tested for functional hybrids in an auxotrophic E.coli host strain. Multiple active hybrid enzymes were identified, including ones fused in regions of low sequence homology.

CKAAPs DB: a conserved key amino acid positions database

The Conserved Key Amino Acid Positions DataBase (CKAAPs DB) provides access to an analysis of structurally similar proteins with dissimilar sequences where key residues within a common fold are identified. The derivation and significance of CKAAPs starting from pairwise structure alignments is described fully in Reddy et al. [Reddy,B.V.B., Li,W.W., Shindyalov,I.N. and Bourne,P.E. (2000) Proteins, in press]. The CKAAPs identified from this theoretical analysis are provided to experimentalists and theoreticians for potential use in protein engineering and modeling. It has been suggested that CKAAPs may be crucial features for protein folding, structural stability and function. Over 170 substructures, as defined by the Combinatorial Extension (CE) database, which are found in approximately 3000 representative polypeptide chains have been analyzed and are available in the CKAAPs DB. CKAAPs DB also provides CKAAPs of the representative set of proteins derived from the CE and FSSP databases. Thus the database contains over 5000 representative poly­peptide chains, covering all known structures in the PDB. A web interface to a relational database permits fast retrieval of structure-sequence alignments, CKAAPs and associated statistics. Users may query by PDB ID, protein name, function and Enzyme Classification number. Users may also submit protein alignments of their own to obtain CKAAPs. An interface to display CKAAPs on each structure from a web browser is also being implemented. CKAAPs DB is maintained by the San Diego Supercomputer Center and accessible at the URL http://ckaaps.sdsc.edu.

Identification and dissection of an enhancer controlling epithelial gene expression in skin

Keratins 14 and 5 are the structural hallmarks of the basal keratinocytes of the epidermis and outer root sheath (ORS) of the hair follicle. Their genes are controlled in a tissue-specific manner and thus serve as useful tools to elucidate the regulatory mechanisms involved in keratinocyte-specific transcription. Previously we identified several keratinocyte-specific DNase I hypersensitive sites (HSs) in the 5′ regulatory sequences of the K14 gene and showed that a 700-bp regulatory domain encompassing HSs II and III can confer epidermal and ORS-specific gene expression in transgenic mice in vivo. Although HS II harbored much of the transactivation activity in vitro, it was not sufficient to restrict expression to keratinocytes in vivo. We now explore the HS III regulatory element. Surprisingly, this element on its own confers gene expression to the keratinocytes of the inner root sheath (IRS) of the hair follicle, whereas a 275-bp DNA fragment containing both HSs II and III shifts the expression from the IRS to the basal keratinocytes and ORS in vivo. Electrophoretic mobility-shift assays and mutational studies of HSs III reveal a role for CACCC-box binding proteins, Sp1 family members, and other factors adding to the list of previously described factors that are involved in keratinocyte-specific gene expression. These studies highlight a cooperative interaction of the two HSs domains and strengthen the importance of combinatorial play of transcription factors that govern keratinocyte-specific gene regulation.

Reverse engineering the (β/α)8 barrel fold

The (β/α)8 barrel is the most commonly occurring fold among protein catalysts. To lay a groundwork for engineering novel barrel proteins, we investigated the amino acid sequence restrictions at 182 structural positions of the prototypical (β/α)8 barrel enzyme triosephosphate isomerase. Using combinatorial mutagenesis and functional selection, we find that turn sequences, α-helix capping and stop motifs, and residues that pack the interface between β-strands and α-helices are highly mutable. Conversely, any mutation of residues in the central core of the β-barrel, β-strand stop motifs, and a single buried salt bridge between amino acids R189 and D227 substantially reduces catalytic activity. Four positions are effectively immutable: conservative single substitutions at these four positions prevent the mutant protein from complementing a triosephosphate isomerase knockout in Escherichia coli. At 142 of the 182 positions, mutation to at least one amino acid of a seven-letter amino acid alphabet produces a triosephosphate isomerase with wild-type activity. Consequently, it seems likely that (β/α)8 barrel structures can be encoded with a subset of the 20 amino acids. Such simplification would greatly decrease the computational burden of (β/α)8 barrel design.

An N-methyl-d-aspartate receptor channel blocker with neuroprotective activity

Excitotoxicity, resulting from sustained activation of glutamate receptors of the N-methyl-d-aspartate (NMDA) subtype, is considered to play a causative role in the etiology of ischemic stroke and several neurodegenerative diseases. The NMDA receptor is therefore a target for the development of neuroprotective agents. Here, we identify an N-benzylated triamine (denoted as NBTA) as a highly selective and potent NMDA-receptor channel blocker selected by screening a reduced dipeptidomimetic synthetic combinatorial library. NBTA blocks recombinant NMDA receptors expressed in Xenopus laevis oocytes with a mean IC50 of 80 nM; in contrast, it does not block GluR1, a glutamate receptor of the non-NMDA subtype. The blocking activity of NBTA on NMDA receptors exhibits the characteristics of an open-channel blocker: (i) no competition with agonists, (ii) voltage dependence, and (iii) use dependence. Significantly, NBTA protects rodent hippocampal neurons from NMDA receptor, but not kainate receptor-mediated excitotoxic cell death, in agreement with its selective action on the corresponding recombinant receptors. Mutagenesis data indicate that the N site, a key asparagine on the M2 transmembrane segment of the NR1 subunit, is the main determinant of the blocker action. The results highlight the potential of this compound as a neuroprotectant.

A selective peroxisome proliferator-activated receptor δ agonist promotes reverse cholesterol transport

The peroxisome proliferator-activated receptors (PPARs) are dietary lipid sensors that regulate fatty acid and carbohydrate metabolism. The hypolipidemic effects of the fibrate drugs and the antidiabetic effects of the glitazone drugs in humans are due to activation of the α (NR1C1) and γ (NR1C3) subtypes, respectively. By contrast, the therapeutic potential of the δ (NR1C2) subtype is unknown, due in part to the lack of selective ligands. We have used combinatorial chemistry and structure-based drug design to develop a potent and subtype-selective PPARδ agonist, GW501516. In macrophages, fibroblasts, and intestinal cells, GW501516 increases expression of the reverse cholesterol transporter ATP-binding cassette A1 and induces apolipoprotein A1-specific cholesterol efflux. When dosed to insulin-resistant middle-aged obese rhesus monkeys, GW501516 causes a dramatic dose-dependent rise in serum high density lipoprotein cholesterol while lowering the levels of small-dense low density lipoprotein, fasting triglycerides, and fasting insulin. Our results suggest that PPARδ agonists may be effective drugs to increase reverse cholesterol transport and decrease cardiovascular disease associated with the metabolic syndrome X.