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.

The area code hypothesis revisited: Olfactory receptors and other related transmembrane receptors may function as the last digits in a cell surface code for assembling embryos

Recent evidence emerging from several laboratories, integrated with new data obtained by searching the genome databases, suggests that the area code hypothesis provides a good heuristic model for explaining the remarkable specificity of cell migration and tissue assembly that occurs throughout embryogenesis. The area code hypothesis proposes that cells assemble organisms, including their brains and nervous systems, with the aid of a molecular-addressing code that functions much like the country, area, regional, and local portions of the telephone dialing system. The complexity of the information required to code cells for the construction of entire organisms is so enormous that we assume that the code must make combinatorial use of members of large multigene families. Such a system would reuse the same receptors as molecular digits in various regions of the embryo, thus greatly reducing the total number of genes required. We present the hypothesis that members of the very large families of olfactory receptors and vomeronasal receptors fulfill the criteria proposed for area code molecules and could serve as the last digits in such a code. We discuss our evidence indicating that receptors of these families are expressed in many parts of developing embryos and suggest that they play a key functional role in cell recognition and targeting not only in the olfactory system but also throughout the brain and numerous other organs as they are assembled.

In vitro generated antibodies specific for telomeric guanine-quadruplex DNA react with Stylonychia lemnae macronuclei

Most eukaryotic telomeres contain a repeating motif with stretches of guanine residues that form a 3′-terminal overhang extending beyond the telomeric duplex region. The telomeric repeat of hypotrichous ciliates, d(T4G4), forms a 16-nucleotide 3′-overhang. Such sequences can adopt parallel-stranded as well as antiparallel-stranded quadruplex conformations in vitro. Although it has been proposed that guanine-quadruplex conformations may have important cellular roles including telomere function, recombination, and transcription, evidence for the existence of this DNA structure in vivo has been elusive to date. We have generated high-affinity single-chain antibody fragment (scFv) probes for the guanine-quadruplex formed by the Stylonychia telomeric repeat, by ribosome display from the Human Combinatorial Antibody Library. Of the scFvs selected, one (Sty3) had an affinity of Kd = 125 pM for the parallel-stranded guanine-quadruplex and could discriminate with at least 1,000-fold specificity between parallel or antiparallel quadruplex conformations formed by the same sequence motif. A second scFv (Sty49) bound both the parallel and antiparallel quadruplex with similar (Kd = 3–5 nM) affinity. Indirect immunofluorescence studies show that Sty49 reacts specifically with the macronucleus but not the micronucleus of Stylonychia lemnae. The replication band, the region where replication and telomere elongation take place, was also not stained, suggesting that the guanine-quadruplex is resolved during replication. Our results provide experimental evidence that the telomeres of Stylonychia macronuclei adopt in vivo a guanine-quadruplex structure, indicating that this structure may have an important role for telomere functioning.

Tbx19, a tissue-selective regulator of POMC gene expression

Pituitary cell types arise in a temporally and spatially specific fashion, in response to combinatorial actions of transcription factors induced by transient signaling gradients. The critical transcriptional determinants of the two pituitary cell types that express the pro-opiomelanocortin (POMC) gene, the anterior lobe corticotropes, producing adrenocorticotropin, and the intermediate lobe melanotropes, producing melanocyte-stimulating hormone (MSHα), have remained unknown. Here, we report that a member of the T-box gene family, Tbx19, which is expressed only in the rostral ventral diencephalon and pituitary gland, commencing on e11.5, marks pituitary cells that will subsequently express the POMC gene and is capable of altering progression of ventral cell types and inducing adrenocorticotropin in rostral tip cells. It is suggested that Tbx19, depending on the presence of synergizing transcription factors, can activate POMC gene expression and repress the α glycoprotein subunit and thyroid-stimulating hormone β promoters.

NOMAD: a versatile strategy for in vitro DNA manipulation applied to promoter analysis and vector design.

Molecular analysis of complex modular structures, such as promoter regions or multi-domain proteins, often requires the creation of families of experimental DNA constructs having altered composition, order, or spacing of individual modules. Generally, creation of every individual construct of such a family uses a specific combination of restriction sites. However, convenient sites are not always available and the alternatives, such as chemical resynthesis of the experimental constructs or engineering of different restriction sites onto the ends of DNA fragments, are costly and time consuming. A general cloning strategy (nucleic acid ordered assembly with directionality, NOMAD; WWW resource locator http:@Lmb1.bios.uic.edu/NOMAD/NOMAD.htm l) is proposed that overcomes these limitations. Use of NOMAD ensures that the production of experimental constructs is no longer the rate-limiting step in applications that require combinatorial rearrangement of DNA fragments. NOMAD manipulates DNA fragments in the form of "modules" having a standardized cohesive end structure. Specially designed "assembly vectors" allow for sequential and directional insertion of any number of modules in an arbitrary predetermined order, using the ability of type IIS restriction enzymes to cut DNA outside of their recognition sequences. Studies of regulatory regions in DNA, such as promoters, replication origins, and RNA processing signals, construction of chimeric proteins, and creation of new cloning vehicles, are among the applications that will benefit from using NOMAD.

Unraveling principles of lead discovery: from unfrustrated energy landscapes to novel molecular anchors.

The search for novel leads is a critical step in the drug discovery process. Computational approaches to identify new lead molecules have focused on discovering complete ligands by evaluating the binding affinity of a large number of candidates, a task of considerable complexity. A new computational method is introduced in this work based on the premise that the primary molecular recognition event in the protein binding site may be accomplished by small core fragments that serve as molecular anchors, providing a structurally stable platform that can be subsequently tailored into complete ligands. To fulfill its role, we show that an effective molecular anchor must meet both the thermodynamic requirement of relative energetic stability of a single binding mode and its consistent kinetic accessibility, which may be measured by the structural consensus of multiple docking simulations. From a large number of candidates, this technique is able to identify known core fragments responsible for primary recognition by the FK506 binding protein (FKBP-12), along with a diverse repertoire of novel molecular cores. By contrast, absolute energetic criteria for selecting molecular anchors are found to be promiscuous. A relationship between a minimum frustration principle of binding energy landscapes and receptor-specific molecular anchors in their role as "recognition nuclei" is established, thereby unraveling a mechanism of lead discovery and providing a practical route to receptor-biased computational combinatorial chemistry.

Gene recognition via spliced sequence alignment.

Gene recognition is one of the most important problems in computational molecular biology. Previous attempts to solve this problem were based on statistics, and applications of combinatorial methods for gene recognition were almost unexplored. Recent advances in large-scale cDNA sequencing open a way toward a new approach to gene recognition that uses previously sequenced genes as a clue for recognition of newly sequenced genes. This paper describes a spliced alignment algorithm and software tool that explores all possible exon assemblies in polynomial time and finds the multiexon structure with the best fit to a related protein. Unlike other existing methods, the algorithm successfully recognizes genes even in the case of short exons or exons with unusual codon usage; we also report correct assemblies for genes with more than 10 exons. On a test sample of human genes with known mammalian relatives, the average correlation between the predicted and actual proteins was 99%. The algorithm correctly reconstructed 87% of genes and the rare discrepancies between the predicted and real exon-intron structures were caused either by short (less than 5 amino acids) initial/terminal exons or by alternative splicing. Moreover, the algorithm predicts human genes reasonably well when the homologous protein is nonvertebrate or even prokaryotic. The surprisingly good performance of the method was confirmed by extensive simulations: in particular, with target proteins at 160 accepted point mutations (PAM) (25% similarity), the correlation between the predicted and actual genes was still as high as 95%.

A strategy of exon shuffling for making large peptide repertoires displayed on filamentous bacteriophage.

It has been suggested that recombination and shuffling between exons has been a key feature in the evolution of proteins. We propose that this strategy could also be used for the artificial evolution of proteins in bacteria. As a first step, we illustrate the use of a self-splicing group I intron with inserted lox-Cre recombination site to assemble a very large combinatorial repertoire (> 10(11) members) of peptides from two different exons. Each exon comprised a repertoire of 10 random amino acids residues; after splicing, the repertoires were joined together through a central five-residue spacer to give a combinatorial repertoire of 25-residue peptides. The repertoire was displayed on filamentous bacteriophage by fusion to the pIII phage coat protein and selected by binding to several proteins, including beta-glucuronidase. One of the peptides selected against beta-glucuronidase was chemically synthesized and shown to inhibit the enzymatic activity (inhibition constant: 17 nM); by further exon shuffling, an improved inhibitor was isolated (inhibition constant: 7 nM). Not only does this approach provide the means for making very large peptide repertoires, but we anticipate that by introducing constraints in the sequences of the peptides and of the linker, it may be possible to evolve small folded peptides and proteins.

Mutational analysis of NM23-H2/NDP kinase identifies the structural domains critical to recognition of a c-myc regulatory element.

NM23-H2, a presumed regulator of tumor metastasis in humans, is a hexameric protein with both enzymatic (NDP kinase) and regulatory (transcriptional activation) activity. While the structure and catalytic mechanisms have been well characterized, the mode of DNA binding is not known. We examined this latter function in a site-directed mutational study and identified residues and domains essential for the recognition of a c-myc regulatory sequence. Three amino acids, Arg-34, Asn-69, and Lys-135, were found among 30 possibilities to be critical for DNA binding. Two of these, Asn-69 and Lys-135, are not conserved between NM23 variants differing in DNA-binding potential, suggesting that DNA recognition resides partly in nonconserved amino acids. All three DNA-binding defective mutant proteins are active enzymatically and appear to be stable hexamers, suggesting that they perform at the level of DNA recognition and that separate functional domains exist for enzyme catalysis and DNA binding. In the context of the known crystal structure of NM23-H2, the DNA-binding residues are located within distinct structural motifs in the monomer, which are exposed to the surface near the 2-fold axis of adjacent subunits in the hexamer. These findings are explained by a model in which NM23-H2 binds DNA with a combinatorial surface consisting of the "outer" face of the dimer. Chemical crosslinking data support a dimeric DNA-binding mode by NM23-H2.

Calcium-dependent oligonucleotide antagonists specific for L-selectin.

The selectins are calcium-dependent C-type lectins that recognize complex anionic carbohydrate ligands, initiating many cell-cell interactions in the vascular system. Selectin blockade shows therapeutic promise in a variety of inflammatory and postischemic pathologies. However, the available oligosaccharide ligand mimetics have low affinities and show cross-reaction among the three selectins, precluding efficient and specific blockade. The SELEX (systematic evolution of ligands by exponential enrichment) process uses combinatorial chemistry and in vitro selection to yield high affinity oligonucleotides with unexpected binding specificities. Nuclease-stabilized randomized oligonucleotides subjected to SELEX against recombinant L-selectin yielded calcium-dependent antagonists with approximately 10(5) higher affinity than the conventional oligosaccharide ligand sialyl LewisX. Most of the isolated ligands shared a common consensus sequence. Unlike sialyl LewisX, these antagonists show little binding to E- or P-selectin. Moreover, they show calcium-dependent binding to native L-selectin on peripheral blood lymphocytes and block L-selectin-dependent interactions with the natural ligands on high endothelial venules.

Application of capillary electrophoresis-electrospray ionization mass spectometry in the determination of molecular diversity.

By means of capillary electrophoresis coupled online to electrospray ionization MS, a library of theoretically 171 disubstituted xanthene derivatives was analyzed. The method allowed the purity and makeup of the library to be determined: 160 of the expected compounds were found to be present, and 12 side-products were also detected in the mixture. Due to the ability of capillary electrophoresis to separate analytes on the basis of charge, most of the xanthene derivatives could be resolved by simple capillary electrophoresis-MS procedures even though 124 of the 171 theoretical compounds were isobaric with at least one other molecule in the mixture. Any remaining unresolved peaks were resolved by MS/MS experiments. The method shows promise for the analysis of small combinatorial libraries with fewer than 1000 components.

Dimerization specificity of Arabidopsis MADS domain homeotic proteins APETALA1, APETALA3, PISTILLATA, and AGAMOUS.

The MADS domain homeotic proteins APETALA1 (AP1), APETALA3 (AP3), PISTILLATA (PI), and AGAMOUS (AG) act in a combinatorial manner to specify the identity of Arabidopsis floral organs. The molecular basis for this combinatorial mode of action was investigated. Immunoprecipitation experiments indicate that all four proteins are capable of interacting with each other. However, these proteins exhibit "partner-specificity" for the formation of DNA-binding dimers; only AP1 homodimers, AG homodimers, and AP3/PI heterodimers are capable of binding to CArG-box sequences. Both the AP3/PI heterodimer and the AP1 or AG homodimers are formed when the three corresponding proteins are present together. The use of chimeric proteins formed by domain swapping indicates that the L region (which follows the MADS box) constitutes a key molecular determinant for the selective formation of DNA-binding dimers. The implications of these results for the ABC genetic model of flower development are discussed.

Hyperplane arrangements, interval orders, and trees.

A hyperplane arrangement is a finite set of hyperplanes in a real affine space. An especially important arrangement is the braid arrangement, which is the set of all hyperplanes xi - xj = 1, 1 combinatorial properties of certain deformations of the braid arrangement are surveyed. In particular, there are unexpected connections with the theory of interval orders and with the enumeration of trees. For instance, the number of labeled interval orders that can be obtained from n intervals I1,..., In of generic lengths is counted. There is also discussed an arrangement due to N. Linial whose number of regions is the number of alternating (or intransitive) trees, as defined by Gelfand, Graev, and Postnikov [Gelfand, I. M., Graev, M. I., and Postnikov, A. (1995), preprint]. Finally, a refinement is given, related to counting labeled trees by number of inversions, of a result of Shi [Shi, J.-Y. (1986), Lecture Notes in Mathematics, no. 1179, Springer-Verlag] that a certain deformation of the braid arrangement has (n + 1)n-1 regions.