Statistical prediction of single-stranded regions in RNA secondary structure and application to predicting effective antisense target sites and beyond

Single-stranded regions in RNA secondary structure are important for RNA–RNA and RNA–protein interactions. We present a probability profile approach for the prediction of these regions based on a statistical algorithm for sampling RNA secondary structures. For the prediction of phylogenetically-determined single-stranded regions in secondary structures of representative RNA sequences, the probability profile offers substantial improvement over the minimum free energy structure. In designing antisense oligonucleotides, a practical problem is how to select a secondary structure for the target mRNA from the optimal structure(s) and many suboptimal structures with similar free energies. By summarizing the information from a statistical sample of probable secondary structures in a single plot, the probability profile not only presents a solution to this dilemma, but also reveals ‘well-determined’ single-stranded regions through the assignment of probabilities as measures of confidence in predictions. In antisense application to the rabbit β-globin mRNA, a significant correlation between hybridization potential predicted by the probability profile and the degree of inhibition of in vitro translation suggests that the probability profile approach is valuable for the identification of effective antisense target sites. Coupling computational design with DNA–RNA array technique provides a rational, efficient framework for antisense oligonucleotide screening. This framework has the potential for high-throughput applications to functional genomics and drug target validation.

Incipient species formation in salamanders of the Ensatina complex

The Ensatina eschscholtzii complex of plethodontid salamanders, a well-known “ring species,” is thought to illustrate stages in the speciation process. Early research, based on morphology and coloration, has been extended by the incorporation of studies of protein variation and mitochondrial DNA sequences. The new data show that the complex includes a number of geographically and genetically distinct components that are at or near the species level. The complex is old and apparently has undergone instances of range contraction, isolation, differentiation, and then expansion and secondary contact. While the hypothesis that speciation is retarded by gene flow around the ring is not supported by molecular data, the general biogeographical hypothesis is supported. There is evidence of a north to south range expansion along two axes, with secondary contact and completion of the ring in southern California. Current research targets regions once thought to show primary intergradation, but which molecular markers reveal to be zones of secondary contact. Here emphasis is on the subspecies E. e. xanthoptica, which is involved in four distinct secondary contacts in central California. There is evidence of renewed genetic interactions upon recontact, with greater genetic differentiation within xanthoptica than between it and some of the interacting populations. The complex presents a full array of intermediate conditions between well-marked species and geographically variable populations. Geographically differentiated segments represent a diversity of depths of time of isolation and admixture, reflecting the complicated geomorphological history of California. Ensatina illustrates the continuing difficulty in making taxonomic assignments in complexes studied during species formation.

Lignin monomer composition is determined by the expression of a cytochrome P450-dependent monooxygenase in Arabidopsis

The phenylpropanoid pathway provides precursors for the biosynthesis of soluble secondary metabolites and lignin in plants. Ferulate-5-hydroxylase (F5H) catalyzes an irreversible hydroxylation step in this pathway that diverts ferulic acid away from guaiacyl lignin biosynthesis and toward sinapic acid and syringyl lignin. This fact led us to postulate that F5H was a potential regulatory step in the determination of lignin monomer composition. To test this hypothesis, we have used Arabidopsis to examine the impact of F5H overexpression. Arabidopsis is a useful model system in which to study lignification because in wild-type plants, guaiacyl and syringyl lignins are deposited in a tissue-specific fashion, while the F5H-deficient fah1 mutant accumulates only guaiacyl lignin. Here we show that ectopic overexpression of F5H in Arabidopsis abolishes tissue-specific lignin monomer accumulation. Surprisingly, overexpression of F5H under the control of the lignification-associated cinnamate-4-hydroxylase promoter, but not the commonly employed cauliflower mosaic virus 35S promoter, generates a lignin that is almost entirely comprised of syringylpropane units. These experiments demonstrate that modification of F5H expression may enable engineering of lignin monomer composition in agronomically important plant species.

Transition-state structure as a unifying basis in protein-folding mechanisms: Contact order, chain topology, stability, and the extended nucleus mechanism

I attempt to reconcile apparently conflicting factors and mechanisms that have been proposed to determine the rate constant for two-state folding of small proteins, on the basis of general features of the structures of transition states. Φ-Value analysis implies a transition state for folding that resembles an expanded and distorted native structure, which is built around an extended nucleus. The nucleus is composed predominantly of elements of partly or well-formed native secondary structure that are stabilized by local and long-range tertiary interactions. These long-range interactions give rise to connecting loops, frequently containing the native loops that are poorly structured. I derive an equation that relates differences in the contact order of a protein to changes in the length of linking loops, which, in turn, is directly related to the unfavorable free energy of the loops in the transition state. Kinetic data on loop extension mutants of CI2 and α-spectrin SH3 domain fit the equation qualitatively. The rate of folding depends primarily on the interactions that directly stabilize the nucleus, especially those in native-like secondary structure and those resulting from the entropy loss from the connecting loops, which vary with contact order. This partitioning of energy accounts for the success of some algorithms that predict folding rates, because they use these principles either explicitly or implicitly. The extended nucleus model thus unifies the observations of rate depending on both stability and topology.

Mutation R120G in αB-crystallin, which is linked to a desmin-related myopathy, results in an irregular structure and defective chaperone-like function

αB-crystallin, a member of the small heat shock protein family, possesses chaperone-like function. Recently, it has been shown that a missense mutation in αB-crystallin, R120G, is genetically linked to a desmin-related myopathy as well as to cataracts [Vicart, P., Caron, A., Guicheney, P., Li, A., Prevost, M.-C., Faure, A., Chateau, D., Chapon, F., Tome, F., Dupret, J.-M., et al. (1998) Nat. Genet. 20, 92–95]. By using α-lactalbumin, alcohol dehydrogenase, and insulin as target proteins, in vitro assays indicated that R120G αB-crystallin had reduced or completely lost chaperone-like function. The addition of R120G αB-crystallin to unfolding α-lactalbumin enhanced the kinetics and extent of its aggregation. R120G αB-crystallin became entangled with unfolding α-lactalbumin and was a major portion of the resulting insoluble pellet. Similarly, incubation of R120G αB-crystallin with alcohol dehydrogenase and insulin also resulted in the presence of R120G αB-crystallin in the insoluble pellets. Far and near UV CD indicate that R120G αB-crystallin has decreased β-sheet secondary structure and an altered aromatic residue environment compared with wild-type αB-crystallin. The apparent molecular mass of R120G αB-crystallin, as determined by gel filtration chromatography, is 1.4 MDa, which is more than twice the molecular mass of wild-type αB-crystallin (650 kDa). Images obtained from cryoelectron microscopy indicate that R120G αB-crystallin possesses an irregular quaternary structure with an absence of a clear central cavity. The results of this study show, through biochemical analysis, that an altered structure and defective chaperone-like function of αB-crystallin are associated with a point mutation that leads to a desmin-related myopathy and cataracts.

Identification of eIF2Bγ and eIF2γ as cofactors of hepatitis C virus internal ribosome entry site-mediated translation using a functional genomics approach

The 5′-untranslated region of hepatitis C virus (HCV) is highly conserved, folds into a complex secondary structure, and functions as an internal ribosome entry site (IRES) to initiate translation of HCV proteins. We have developed a selection system based on a randomized hairpin ribozyme gene library to identify cellular factors involved in HCV IRES function. A retroviral vector ribozyme library with randomized target recognition sequences was introduced into HeLa cells, stably expressing a bicistronic construct encoding the hygromycin B phosphotransferase gene and the herpes simplex virus thymidine kinase gene (HSV-tk). Translation of the HSV-tk gene was mediated by the HCV IRES. Cells expressing ribozymes that inhibit HCV IRES-mediated translation of HSV-tk were selected via their resistance to both ganciclovir and hygromycin B. Two ribozymes reproducibly conferred the ganciclovir-resistant phenotype and were shown to inhibit IRES-mediated translation of HCV core protein but did not inhibit cap-dependent protein translation or cell growth. The functional targets of these ribozymes were identified as the gamma subunits of human eukaryotic initiation factors 2B (eIF2Bγ) and 2 (eIF2γ), respectively. The involvement of eIF2Bγ and eIF2γ in HCV IRES-mediated translation was further validated by ribozymes directed against additional sites within the mRNAs of these genes. In addition to leading to the identification of cellular IRES cofactors, ribozymes obtained from this cellular selection system could be directly used to specifically inhibit HCV viral translation, thereby facilitating the development of new antiviral strategies for HCV infection.

Retrovirus-mediated gene transfer of MLL-ELL transforms primary myeloid progenitors and causes acute myeloid leukemias in mice

The MLL-ELL fusion gene results from the translocation t(11;19)(q23;p13.1) that is associated with de novo and therapy-related acute myeloid leukemia. To study its transforming properties, we retrovirally transduced primary murine hematopoietic progenitors and assessed their growth properties both in vitro and in vivo. MLL-ELL increased the proliferation of myeloid colony-forming cells in methylcellulose cultures upon serial replating, whereas overexpression of ELL alone had no effect. We reconstituted lethally irradiated congenic mice with bone marrow progenitors transduced with MLL-ELL or the control MIE vector encoding the enhanced green fluorescent protein. When the peripheral blood of the mice was analyzed 11–13 weeks postreconstitution, we found that the engraftment of the MLL-ELL-transduced cells was superior to that of the MIE controls. At this time point, the contribution of the donor cells was normally distributed among the myeloid and nonmyeloid compartments. Although all of the MIE animals (n = 10) remained healthy for more than a year, all of the MLL-ELL mice (n = 20) succumbed to monoclonal or pauciclonal acute myeloid leukemias within 100–200 days. The leukemic cells were readily transplantable to secondary recipients and could be established as immortalized cell lines in liquid cultures. These studies demonstrate the enhancing effect of MLL-ELL on the proliferative potential of myeloid progenitors as well as its causal role in the genesis of acute myeloid leukemias.

Etoposide induces heritable chromosomal aberrations and aneuploidy during male meiosis in the mouse

Etoposide, a topoisomerase II inhibitor widely used in cancer therapy, is suspected of inducing secondary tumors and affecting the genetic constitution of germ cells. A better understanding of the potential heritable risk of etoposide is needed to provide sound genetic counseling to cancer patients treated with this drug in their reproductive years. We used a mouse model to investigate the effects of clinical doses of etoposide on the induction of chromosomal abnormalities in spermatocytes and their transmission to zygotes by using a combination of chromosome painting and 4′,6-diamidino-2-phenylindole staining. High frequencies of chromosomal aberrations were detected in spermatocytes within 64 h after treatment when over 30% of the metaphases analyzed had structural aberrations (P < 0.01). Significant increases in the percentages of zygotic metaphases with structural aberrations were found only for matings that sampled treated pachytene (28-fold, P < 0.0001) and preleptotene spermatocytes (13-fold, P < 0.001). Etoposide induced mostly acentric fragments and deletions, types of aberrations expected to result in embryonic lethality, because they represent loss of genetic material. Chromosomal exchanges were rare. Etoposide treatment of pachytene cells induced aneuploidy in both spermatocytes (18-fold, P < 0.01) and zygotes (8-fold, P < 0.05). We know of no other report of an agent for which paternal exposure leads to an increased incidence of aneuploidy in the offspring. Thus, we found that therapeutic doses of etoposide affect primarily meiotic germ cells, producing unstable structural aberrations and aneuploidy, effects that are transmitted to the progeny. This finding suggests that individuals who undergo chemotherapy with etoposide may be at a higher risk for abnormal reproductive outcomes especially within the 2 months after chemotherapy.

Alloantigen-induced unresponsiveness in cord blood T lymphocytes is associated with defective activation of Ras

Human umbilical cord blood T lymphocytes (CBTL) respond to primary allostimulation but they do not proliferate upon rechallenge with alloantigen. Using PKH-26-labeled cells created a proliferative block that was observed only in CBTL that have divided during primary stimulation (PKH-26dim) but not in unstimulated (PKH-26bright) CBTL. CBTL’s secondary unresponsiveness resembles anergy and can be overcome by treatment with phorbol myristate acetate (PMA) and ionomycin or by high doses (50–100 units/ml) of interleukin 2. Addition of interleukin 2 to the primary cultures does not prevent the induction of secondary unresponsiveness. Defective Ras activation is detected in PKH-26dim CBTL during secondary response to alloantigen or after antibody-mediated T cell receptor stimulation whereas Ras is activated and proliferation is induced in CBTL during primary alloantigenic stimulation. Upon stimulation with PMA plus ionomycin, PMA plus alloantigen, but not alloantigen plus ionomycin, Ras is activated in PKH-26dim CBTL, and the block in proliferation is overcome. Correction of PKH-26dim CBTL’s proliferative defect correlates with PMA-induced Ras activation, suggesting a defect in the signaling pathway leading to Ras. Ras-independent signals, necessary but not sufficient to induce PKH-26dim CBTL proliferation, are provided by alloantigen exposure, as evident by the ability of PMA plus alloantigen but not PMA alone to overcome the proliferative block. Functional signal transduction through CD28 in PKH-26dim CBTL is supported by detectable CD28-mediated PI-3 kinase activation after PKH-26dim CBTL’s exposure to alloantigen or CD28 cross-linking. These results suggest that defective activation of Ras plays a key role in PKH-26dim CBTL’s secondary unresponsiveness and point to a defect along the T cell receptor rather than the CD28 signaling pathway.

Overexpression of a glutamate receptor (GluR2) ligand binding domain in Escherichia coli: Application of a novel protein folding screen

Expression of the S1S2 ligand binding domain [Kuusinen, A., Arvola, M. & Keinänen, K. (1995) EMBO J. 14, 6327–6332] of the rat α-amino-3-hydroxy-5-methylisoxazole-4-propionic acid-selective glutamate receptor GluR2 in Escherichia coli under control of a T7 promoter leads to production of >100 mg/liter of histidine-tagged S1S2 protein (HS1S2) in the form of inclusion bodies. Using a novel fractional factorial folding screen and a rational, step-by-step approach, multiple conditions were determined for the folding of the HS1S2 α-amino-3-hydroxy-5-methylisoxazole-4-propionic acid binding domain. Characterization of the HS1S2 ligand binding domain showed that it is water-soluble, monomeric, has significant secondary structure, and is sensitive to trypsinolysis at sites close to the beginning of the putative transmembrane regions. Application of a fractional factorial folding screen to other proteins may provide a useful means to evaluate E. coli as an economical and convenient expression host.

tmRDB (tmRNA database)

The tmRNA database (tmRDB) is maintained at the University of Texas Health Science Center at Tyler, Texas, and accessible on the World Wide Web at the URL http://psyche.uthct.edu/dbs/tmRDB/tmRDB.html. Mirror sites are located at Auburn University, Auburn, Alabama (http://www.ag.auburn.edu/mirror/tmRDB/) and the Institute of Biological Sciences, Aarhus, Denmark (http://www.bioinf.au.dk/tmRDB/). The tmRDB provides information and citation links about tmRNA, a molecule that combines functions of tRNA and mRNA in trans-translation. tmRNA is likely to be present in all bacteria and has been found in algae chloroplasts, the cyanelle of Cyanophora paradoxa and the mitochondrion of the flagellate Reclinomonas americana. This release adds 26 new sequences and corresponding predicted tmRNA-encoded tag peptides for a total of 86 tmRNAs, ordered alphabetically and phylogenetically. Secondary structures and three-dimensional models in PDB format for representative molecules are being made available. tmRNA alignments prove individual base pairs and are generated manually assisted by computational tools. The alignments with their corresponding structural annotation can be obtained in various formats, including a new column format designed to improve and simplify computational usability of the data.

Systematic Reverse Genetics of Transfer-DNA-Tagged Lines of Arabidopsis1 : Isolation of Mutations in the Cytochrome P450 Gene Superfamily

We have developed an efficient reverse-genetics protocol that uses expedient pooling and hybridization strategies to identify individual transfer-DNA insertion lines from a collection of 6000 independently transformed lines in as few as 36 polymerase chain reactions. We have used this protocol to systematically isolate Arabidopsis lines containing insertional mutations in individual cytochrome P450 genes. In higher plants P450 genes encode enzymes that perform an exceptionally wide range of functions, including the biosynthesis of primary metabolites necessary for normal growth and development, the biosynthesis of secondary products, and the catabolism of xenobiotics. Despite their importance, progress in assigning enzymatic function to individual P450 gene products has been slow. Here we report the isolation of the first 12 such lines, including one (CYP83B1-1) that displays a runt phenotype (small plants with hooked leaves), and three insertions in abundantly expressed genes. The DNAs used in this study are publicly available and can be used to systematically isolate mutants in Arabidopsis.

Developmental Regulation of Intercellular Protein Trafficking through Plasmodesmata in Tobacco Leaf Epidermis1

Plasmodesmata mediate direct cell-to-cell communication in plants. One of their significant features is that primary plasmodesmata formed at the time of cytokinesis often undergo structural modifications, by the de novo addition of cytoplasmic strands across cell walls, to become complex secondary plasmodesmata during plant development. Whether such modifications allow plasmodesmata to gain special transport functions has been an outstanding issue in plant biology. Here we present data showing that the cucumber mosaic virus 3a movement protein (MP):green fluorescent protein (GFP) fusion was not targeted to primary plasmodesmata in the epidermis of young or mature leaves in transgenic tobacco (Nicotiana tabacum) plants constitutively expressing the 3a:GFP fusion gene. Furthermore, the cucumber mosaic virus 3a MP:GFP fusion protein produced in planta by biolistic bombardment of the 3a:GFP fusion gene did not traffic between cells interconnected by primary plasmodesmata in the epidermis of a young leaf. In contrast, the 3a MP:GFP was targeted to complex secondary plasmodesmata and trafficked from cell to cell when a leaf reached a certain developmental stage. These data provide the first experimental evidence, to our knowledge, that primary and complex secondary plasmodesmata have different protein-trafficking functions and suggest that complex secondary plasmodesmata may be formed to traffic specific macromolecules that are important for certain stages of leaf development.

Indole-3-Acetic Acid Controls Cambial Growth in Scots Pine by Positional Signaling1

The vascular cambium produces secondary xylem and phloem in plants and is responsible for wood formation in forest trees. In this study we used a microscale mass-spectrometry technique coupled with cryosectioning to visualize the radial concentration gradient of endogenous indole-3-acetic acid (IAA) across the cambial meristem and the differentiating derivatives in Scots pine (Pinus sylvestris L.) trees that had different rates of cambial growth. This approach allowed us to investigate the relationship between growth rate and the concentration of endogenous IAA in the dividing cells. We also tested the hypothesis that IAA is a positional signal in xylem development (C. Uggla, T. Moritz, G. Sandberg, B. Sundberg [1996] Proc Natl Acad Sci USA 93: 9282–9286). This idea postulates that the width of the radial concentration gradient of IAA regulates the radial number of dividing cells in the cambial meristem, which is an important component for determining cambial growth rate. The relationship between IAA concentration in the dividing cells and growth rate was poor, although the highest IAA concentration was observed in the fastest-growing cambia. The radial width of the IAA concentration gradient showed a strong correlation with cambial growth rate. The results indicate that IAA gives positional information in plants.

Prediction of protein folding class using global description of amino acid sequence.

We present a method for predicting protein folding class based on global protein chain description and a voting process. Selection of the best descriptors was achieved by a computer-simulated neural network trained on a data base consisting of 83 folding classes. Protein-chain descriptors include overall composition, transition, and distribution of amino acid attributes, such as relative hydrophobicity, predicted secondary structure, and predicted solvent exposure. Cross-validation testing was performed on 15 of the largest classes. The test shows that proteins were assigned to the correct class (correct positive prediction) with an average accuracy of 71.7%, whereas the inverse prediction of proteins as not belonging to a particular class (correct negative prediction) was 90-95% accurate. When tested on 254 structures used in this study, the top two predictions contained the correct class in 91% of the cases.