Recent improvements in prediction of protein structure by global optimization of a potential energy function

Recent improvements of a hierarchical ab initio or de novo approach for predicting both α and β structures of proteins are described. The united-residue energy function used in this procedure includes multibody interactions from a cumulant expansion of the free energy of polypeptide chains, with their relative weights determined by Z-score optimization. The critical initial stage of the hierarchical procedure involves a search of conformational space by the conformational space annealing (CSA) method, followed by optimization of an all-atom model. The procedure was assessed in a recent blind test of protein structure prediction (CASP4). The resulting lowest-energy structures of the target proteins (ranging in size from 70 to 244 residues) agreed with the experimental structures in many respects. The entire experimental structure of a cyclic α-helical protein of 70 residues was predicted to within 4.3 Å α-carbon (Cα) rms deviation (rmsd) whereas, for other α-helical proteins, fragments of roughly 60 residues were predicted to within 6.0 Å Cα rmsd. Whereas β structures can now be predicted with the new procedure, the success rate for α/β- and β-proteins is lower than that for α-proteins at present. For the β portions of α/β structures, the Cα rmsd's are less than 6.0 Å for contiguous fragments of 30–40 residues; for one target, three fragments (of length 10, 23, and 28 residues, respectively) formed a compact part of the tertiary structure with a Cα rmsd less than 6.0 Å. Overall, these results constitute an important step toward the ab initio prediction of protein structure solely from the amino acid sequence.

Top-down morphogenesis of colorectal tumors

One of the fundamental tenets of oncology is that tumors arise from stem cells. In the colon, stem cells are thought to reside at the base of crypts. In the early stages of tumorigenesis, however, dysplastic cells are routinely found at the luminal surface of the crypts whereas the cells at the bases of these same crypts appear morphologically normal. To understand this discrepancy, we evaluated the molecular characteristics of cells isolated from the bases and orifices of the same crypts in small colorectal adenomas. We found that the dysplastic cells at the tops of the crypts often exhibited genetic alterations of adenomatous polyposis coli (APC) and neoplasia-associated patterns of gene expression. In contrast, cells located at the base of these same crypts did not contain such alterations and were not clonally related to the contiguous transformed cells above them. These results imply that development of adenomatous polyps proceeds through a top-down mechanism. Genetically altered cells in the superficial portions of the mucosae spread laterally and downward to form new crypts that first connect to preexisting normal crypts and eventually replace them.

Toward optimized carbohydrate-based anticancer vaccines: Epitope clustering, carrier structure, and adjuvant all influence antibody responses to Lewisy conjugates in mice

The feasibility of using carbohydrate-based vaccines for the immunotherapy of cancer is being actively explored at the present time. Although a number of clinical trials have already been conducted with glycoconjugate vaccines, the optimal design and composition of the vaccines has yet to be determined. Among the candidate antigens being examined is Lewisy (Ley), a blood group-related antigen that is overexpressed on the majority of human carcinomas. Using Ley as a model for specificity, we have examined the role of epitope clustering, carrier structure, and adjuvant on the immunogenicity of Ley conjugates in mice. A glycolipopeptide containing a cluster of three contiguous Ley-serine epitopes and the Pam3Cys immunostimulating moiety was found to be superior to a similar construct containing only one Ley-serine epitope in eliciting antitumor cell antibodies. Because only IgM antibodies were produced by this vaccine, the effect on immunogenicity of coupling the glycopeptide to keyhole limpet hemocyanin was examined; although both IgM and IgG antibodies were formed, the antibodies reacted only with the immunizing structure. Reexamination of the clustered Ley-serine Pam3Cys conjugate with the adjuvant QS-21 resulted in the identification of both IgG and IgM antibodies reacting with tumor cells, thus demonstrating the feasibility of an entirely synthetic carbohydrate-based anticancer vaccine in an animal model.

Grass genomes

For the most part, studies of grass genome structure have been limited to the generation of whole-genome genetic maps or the fine structure and sequence analysis of single genes or gene clusters. We have investigated large contiguous segments of the genomes of maize, sorghum, and rice, primarily focusing on intergenic spaces. Our data indicate that much (>50%) of the maize genome is composed of interspersed repetitive DNAs, primarily nested retrotransposons that insert between genes. These retroelements are less abundant in smaller genome plants, including rice and sorghum. Although 5- to 200-kb blocks of methylated, presumably heterochromatic, retrotransposons flank most maize genes, rice and sorghum genes are often adjacent. Similar genes are commonly found in the same relative chromosomal locations and orientations in each of these three species, although there are numerous exceptions to this collinearity (i.e., rearrangements) that can be detected at the levels of both the recombinational map and cloned DNA. Evolutionarily conserved sequences are largely confined to genes and their regulatory elements. Our results indicate that a knowledge of grass genome structure will be a useful tool for gene discovery and isolation, but the general rules and biological significance of grass genome organization remain to be determined. Moreover, the nature and frequency of exceptions to the general patterns of grass genome structure and collinearity are still largely unknown and will require extensive further investigation.

Immune-type receptor genes in zebrafish share genetic and functional properties with genes encoded by the mammalian leukocyte receptor cluster

An extensive, highly diversified multigene family of novel immune-type receptor (nitr) genes has been defined in Danio rerio (zebrafish). The genes are predicted to encode type I transmembrane glycoproteins consisting of extracellular variable (V) and V-like C2 (V/C2) domains, a transmembrane region and a cytoplasmic tail. All of the genes examined encode immunoreceptor tyrosine-based inhibition motifs in the cytoplasmic tail. Radiation hybrid panel mapping and analysis of a deletion mutant line (b240) indicate that a minimum of ≈40 nitr genes are contiguous in the genome and span ≈0.6 Mb near the top of zebrafish linkage group 7. One flanking region of the nitr gene complex shares conserved synteny with a region of mouse chromosome 7, which shares conserved synteny with human 19q13.3-q13.4 that encodes the leukocyte receptor cluster. Antibody-induced crosslinking of Nitrs that have been introduced into a human natural killer cell line inhibits the phosphorylation of mitogen-activated protein kinase that is triggered by natural killer-sensitive tumor target cells. Nitrs likely represent intermediates in the evolution of the leukocyte receptor cluster.

A statistical mechanical model for β-hairpin kinetics

Understanding the mechanism of protein secondary structure formation is an essential part of the protein-folding puzzle. Here, we describe a simple statistical mechanical model for the formation of a β-hairpin, the minimal structural element of the antiparallel β-pleated sheet. The model accurately describes the thermodynamic and kinetic behavior of a 16-residue, β-hairpin-forming peptide, successfully explaining its two-state behavior and apparent negative activation energy for folding. The model classifies structures according to their backbone conformation, defined by 15 pairs of dihedral angles, and is further simplified by considering only the 120 structures with contiguous stretches of native pairs of backbone dihedral angles. This single sequence approximation is tested by comparison with a more complete model that includes the 215 possible conformations and 15 × 215 possible kinetic transitions. Finally, we use the model to predict the equilibrium unfolding curves and kinetics for several variants of the β-hairpin peptide.

The highly recombinogenic bz locus lies in an unusually gene-rich region of the maize genome

The bronze (bz) locus exhibits the highest rate of recombination of any gene in higher plants. To investigate the possible basis of this high rate of recombination, we have analyzed the physical organization of the region around the bz locus. Two adjacent bacterial artificial chromosome clones, comprising a 240-kb contig centered around the Bz-McC allele, were isolated, and 60 kb of contiguous DNA spanning the two bacterial artificial chromosome clones was sequenced. We find that the bz locus lies in an unusually gene-rich region of the maize genome. Ten genes, at least eight of which are shown to be transcribed, are contained in a 32-kb stretch of DNA that is uninterrupted by retrotransposons. We have isolated nearly full length cDNAs corresponding to the five proximal genes in the cluster. The average intertranscript distance between them is just 1 kb, revealing a surprisingly compact packaging of adjacent genes in this part of the genome. At least 11 small insertions, including several previously described miniature inverted repeat transposable elements, were detected in the introns and 3′ untranslated regions of genes and between genes. The gene-rich region is flanked at the proximal and distal ends by retrotransposon blocks. Thus, the maize genome appears to have scattered regions of high gene density similar to those found in other plants. The unusually high rate of intragenic recombination seen in bz may be related to the very high gene density of the region.

Molecular cytogenetic delineation of a novel critical genomic region in chromosome bands 11q22.3-923.1 in lymphoproliferative disorders.

Aberrations of the long arm of chromosome 11 are among the most common chromosome abnormalities in lymphoproliferative disorders (LPD). Translocations involving BCL1 at 11q13 are strongly associated with mantle cell lymphoma. other nonrandom aberrations, especially deletions and, less frequently, translocations, involving bands 11q21-923 have been identified by chromosome banding analysis. To date, the critical genomic segment and candidate genes involved in these deletions have not been identified. In the present study, we have analyzed tumors from 43 patients with LPD (B-cell chronic lymphocytic leukemia, n = 40; mantle cell lymphoma, n = 3) showing aberrations of bands 11q21-923 by fluorescence in situ hybridization. As probes we used Alu-PCR products from 17 yeast artificial chromosome clones spanning chromosome bands 11q14.3-923.3, including a panel of yeast artificial chromosome clones recognizing a contiguous genomic DNA fragment of approximately 9-10 Mb in bands 11q22.3-923.3. In the 41 tumors exhibiting deletions, we identified a commonly deleted segment in band 11q22.3-923.1; this region is approximately 2-3 Mb in size and contains the genes coding for ATM (ataxia telangiectasia mutated), RDX (radixin), and FDX1 (ferredoxin 1). Furthermore, two translocation break-points were localized to a 1.8-Mb genomic fragment contained within the commonly deleted segment. Thus, we have identified a single critical region of 2-3 Mb in size in which 11q14-923 aberrations in LPD cluster. This provides the basis for the identification of the gene(s) at 11q22.3-923.1 that are involved in the pathogenesis of LPD.

The [(G/C)3NN]n motif: a common DNA repeat that excludes nucleosomes.

Nucleosomes, the basic structural elements of chromosomes, consist of 146 bp of DNA coiled around an octamer of histone proteins, and their presence can strongly influence gene expression. Considerations of the anisotropic flexibility of nucleotide triplets containing 3 cytosines or guanines suggested that a [5'(G/C)3 NN3']n motif might resist wrapping around a histone octamer. To test this, DNAs were constructed containing a 5'-CCGNN-3' pentanucleotide repeat with the Ns varied. Using in vitro nucleosome reconstitution and electron microscopy, a plasmid with 48 contiguous CCGNN repeats strongly excluded nucleosomes in the repeat region. Competitive reconstitution gel retardation experiments using DNA fragments containing 12, 24, or 48 CCGNN repeats showed that the propensity to exclude nucleosomes increased with the length of the repeat. Analysis showed that a 268-bp DNA containing a (CCGNN)48 block is 4.9 +/- 0.6-fold less efficient in nucleosome assembly than a similar length pUC19 fragment and approximately 78-fold less efficient than a similar length (CTG)n sequence, based on results from previous studies. Computer searches against the GenBank database for matches with a [(G/C)3NN]48 sequence revealed numerous examples that frequently were present in the control regions of "TATA-less" genes, including the human ETS-2 and human dihydrofolate reductase genes. In both cases the (G/C)3NN repeat, present in the promoter region, co-maps with loci previously shown to be nuclease hypersensitive sites.

Operant conditioning of H-reflex changes synaptic terminals on primate motoneurons.

Operant conditioning of the primate triceps surae H-reflex, the electrical analog of the spinal stretch reflex, creates a memory trace that includes changes in the spinal cord. To define the morphological correlates of this plasticity, we analyzed the synaptic terminal coverage of triceps surae motoneurons from animals in which the triceps surae H-reflex in one leg had been increased (HRup mode) or decreased (HRdown mode) by conditioning and compared them to each other and to motoneurons from unconditioned animals. Motoneurons were labeled by intramuscular injection of cholera toxin-horseradish peroxidase. A total of 5055 terminals on the cell bodies and proximal dendrites of 114 motoneurons from 14 animals were studied by electron microscopy. Significant differences were found between HRup and HRdown animals and between HRup and naive (i.e., unconditioned) animals. F terminals (i.e., putative inhibitory terminals) were smaller and their active zone coverage on the cell body was lower on motoneurons from the conditioned side of HRup animals than on motoneurons from the conditioned side of HRdown animals. C terminals (i.e., terminals associated with postsynaptic cisterns and rough endoplasmic reticulum) were smaller and the number of C terminals in each C complex (i.e., a group of contiguous C terminals) was larger on motoneurons from the conditioned side of HRup animals than on motoneurons either from the conditioned side of HRdown animals or from naive animals. Because the treatment of HRup and HRdown animals differed only in the reward contingency, the results imply that the two contingencies had different effects on motoneuron synaptic terminals. In combination with other recent data, they show that H-reflex conditioning produces a complex pattern of spinal cord plasticity that includes changes in motoneuron physiological properties as well as in synaptic terminals. Further delineation of this pattern should reveal the contribution of the structural changes described here to the learned change in behavior.

DNA analysis and diagnostics on oligonucleotide microchips.

We present a further development in the technology of sequencing by hybridization to oligonucleotide microchips (SHOM) and its application to diagnostics for genetic diseases. A robot has been constructed to manufacture sequencing "microchips." The microchip is an array of oligonucleotides immobilized into gel elements fixed on a glass plate. Hybridization of the microchip with fluorescently labeled DNA was monitored in real time simultaneously for all microchip elements with a two-wavelength fluorescent microscope equipped with a charge-coupled device camera. SHOM has been used to detect beta-thalassemia mutations in patients by hybridizing PCR-amplified DNA with the microchips. A contiguous stacking hybridization technique has been applied for the detection of mutations; it can simplify medical diagnostics and enhance its reliability. The use of multicolor monitoring of contiguous stacking hybridization is suggested for large-scale diagnostics and gene polymorphism studies. Other applications of the SHOM technology are discussed.

Rapid appearance and asymmetric distribution of glucose transporter SGTP4 at the apical surface of intramammalian-stage Schistosoma mansoni.

Adult Schistosoma mansoni blood flukes reside in the mesenteric veins of their vertebrate hosts, where they absorb immense quantities of glucose through their tegument by facilitated diffusion. Previously, we obtained S. mansoni cDNAs encoding facilitated-diffusion schistosome glucose transporter proteins 1 and 4 (SGTP1 and SGTP4) and localized SGTP1 to the basal membranes of the tegument and the underlying muscle. In this study, we characterize the expression and localization of SGTP4 during the schistosome life cycle. Antibodies specific to SGTP4 appear to stain only the double-bilayer, apical membranes of the adult parasite tegument, revealing an asymmetric distribution relative to the basal transporter SGTP1. On living worms, SGTP4 is available to surface biotinylation, suggesting that it is exposed at the hose-parasite interface. SGTP4 is detected shortly after the transformation of free-living, infectious cercariae into schistosomula and coincides with the appearance of the double membrane. Within 15 min after transformation, anti-SGTP4 staining produces a bright, patchy distribution at the surface of schistosomula, which becomes contiguous over the entire surface of the schistosomula by 24 hr after transformation. SGTP4 is not detected in earlier developmental stages (eggs, sporocysts, and cercariae) that do not possess the specialized double membrane. Thus, SGTP4 appears to be expressed only in the mammalian stages of the parasite's life cycle and specifically localized within the host-interactive, apical membranes of the tegument.

Foldons, protein structural modules, and exons.

Foldons, which are kinetically competent, quasi-independently folding units of a protein, may be defined using energy landscape analysis. Foldons can be identified by maxima in a scan of the ratio of a contiguous segment's energetic stability gap to the energy variance of that segment's molten globule states, reflecting the requirement of minimal frustration. The predicted foldons are compared with the exons and structural modules for 16 of the 30 proteins studied. Statistical analysis indicates a strong correlation between the energetically determined foldons and Go's geometrically defined structural modules, but there are marked sequence-dependent effects. There is only a weak correlation of foldons to exons. For gammaII-crystallin, myoglobin, barnase, alpha-lactalbumin, and cytochrome c the foldons and some noncontiguous clusters of foldons compare well with intermediates observed in experiment.

Targeted gene correction of episomal DNA in mammalian cells mediated by a chimeric RNA.DNA oligonucleotide.

An experimental strategy to facilitate correction of single-base mutations of episomal targets in mammalian cells has been developed. The method utilizes a chimeric oligonucleotide composed of a contiguous stretch of RNA and DNA residues in a duplex conformation with double hairpin caps on the ends. The RNA/DNA sequence is designed to align with the sequence of the mutant locus and to contain the desired nucleotide change. Activity of the chimeric molecule in targeted correction was tested in a model system in which the aim was to correct a point mutation in the gene encoding the human liver/bone/kidney alkaline phosphatase. When the chimeric molecule was introduced into cells containing the mutant gene on an extrachromosomal plasmid, correction of the point mutation was accomplished with a frequency approaching 30%. These results extend the usefulness of the oligonucleotide-based gene targeting approaches by increasing specific targeting frequency. This strategy should enable the design of antiviral agents.

A two-subunit type I DNA topoisomerase (reverse gyrase) from an extreme hyperthermophile.

A recently described reverse gyrase from the hyperthermophilic methanogen Methanopyrus kandleri is the only known example of a heterodimeric type I topoisomerase. The enzyme is made up of a 42-kDa subunit which covalently interacts with DNA (RgyA) and a 138-kDa subunit which binds ATP (RgyB). We have now cloned and sequenced the genes for both subunits of this enzyme. Surprisingly, the universally conserved type I topoisomerase domain [Lima, C. D., Wang, J. C. & Mondragon, A. (1994) Nature (London) 367, 138-146] which has been found as a contiguous polypeptide in the prokaryotes and eukaryotes is shared between the protomers. The subdomain with the active-site tyrosine is entirely within RgyA, whereas the subdomain implicated in noncovalent binding of the cleaved DNA strand is contained entirely in RgyB. The appearance of this unique structure in a highly conserved enzyme family supports the hypothesis that the methanogens branched from other prokaryotes and eukaryotes very early in evolution.