Chrysanthemum chlorotic mottle viroid: Unusual structural properties of a subgroup of self-cleaving viroids with hammerhead ribozymes

The causal agent of chrysanthemum chlorotic mottle (CChM) disease has been identified, cloned, and sequenced. It is a viroid RNA (CChMVd) of 398–399 nucleotides. In vitro transcripts with the complete CChMVd sequence were infectious and induced the typical symptoms of the CChM disease. CChMVd can form hammerhead structures in both polarity strands. Plus and minus monomeric CChMVd RNAs self-cleaved during in vitro transcription and after purification as predicted by these structures, which are stable and most probably act as single hammerhead structures as in peach latent mosaic viroid (PLMVd), but not in avocado sunblotch viroid (ASBVd). Moreover, the plus CChMVd hammerhead structure also appears to be active in vivo, because the 5′ terminus of the linear plus CChMVd RNA isolated from infected tissue is that predicted by the corresponding hammerhead ribozyme. Both hammerhead structures of CChMVd display some peculiarities: the plus self-cleaving domain has an unpaired A after the conserved A9 residue, and the minus one has an unusually long helix II. The most stable secondary structure predicted for CChMVd is a branched conformation that does not fulfill the rod-like or quasi-rod-like model proposed for the in vitro structure of most viroids with the exception of PLMVd, whose proposed secondary structure of lowest free energy also is branched. The unusual conformation of CChMVd and PLMVd is supported by their insolubility in 2 M LiCl, in contrast to ASBVd and a series of representative non-self-cleaving viroids that are soluble under the same high salt conditions. These results support the classification of self-cleaving viroids into two subgroups, one formed by ASBVd and the other one by PLMVd and CChMVd.

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.

Structure of the recombinant full-length hamster prion protein PrP(29–231): The N terminus is highly flexible

The prion diseases seem to be caused by a conformational change of the prion protein (PrP) from the benign cellular form PrPC to the infectious scrapie form PrPSc; thus, detailed information about PrP structure may provide essential insights into the mechanism by which these diseases develop. In this study, the secondary structure of the recombinant Syrian hamster PrP of residues 29–231 [PrP(29–231)] is investigated by multidimensional heteronuclear NMR. Chemical shift index analysis and nuclear Overhauser effect data show that PrP(29–231) contains three helices and possibly one short β-strand. Most striking is the random-coil nature of chemical shifts for residues 30–124 in the full-length PrP. Although the secondary structure elements are similar to those found in mouse PrP fragment PrP(121–231), the secondary structure boundaries of PrP(29–231) are different from those in mouse PrP(121–231) but similar to those found in the structure of Syrian hamster PrP(90–231). Comparison of resonance assignments of PrP(29–231) and PrP(90–231) indicates that there may be transient interactions between the additional residues and the structured core. Backbone dynamics studies done by using the heteronuclear [1H]-15N nuclear Overhauser effect indicate that almost half of PrP(29–231), residues 29–124, is highly flexible. This plastic region could feature in the conversion of PrPC to PrPSc by template-assisted formation of β-structure.

Protein–RNA interactions: a structural analysis

A detailed computational analysis of 32 protein–RNA complexes is presented. A number of physical and chemical properties of the intermolecular interfaces are calculated and compared with those observed in protein–double-stranded DNA and protein–single-stranded DNA complexes. The interface properties of the protein–RNA complexes reveal the diverse nature of the binding sites. van der Waals contacts played a more prevalent role than hydrogen bond contacts, and preferential binding to guanine and uracil was observed. The positively charged residue, arginine, and the single aromatic residues, phenylalanine and tyrosine, all played key roles in the RNA binding sites. A comparison between protein–RNA and protein–DNA complexes showed that whilst base and backbone contacts (both hydrogen bonding and van der Waals) were observed with equal frequency in the protein–RNA complexes, backbone contacts were more dominant in the protein–DNA complexes. Although similar modes of secondary structure interactions have been observed in RNA and DNA binding proteins, the current analysis emphasises the differences that exist between the two types of nucleic acid binding protein at the atomic contact level.

GOBASE: the organelle genome database

GOBASE (http://megasun.bch.umontreal.ca/gobase/) is a network-accessible biological database, which is unique in bringing together diverse biological data on organelles with taxonomically broad coverage, and in furnishing data that have been exhaustively verified and completed by experts. So far, we have focused on mitochondrial data: GOBASE contains all published nucleotide and protein sequences encoded by mitochondrial genomes, selected RNA secondary structures of mitochondria-encoded molecules, genetic maps of completely sequenced genomes, taxonomic information for all species whose sequences are present in the database and organismal descriptions of key protistan eukaryotes. All of these data have been integrated and organized in a formal database structure to allow sophisticated biological queries using terms that are inherent in biological concepts. Most importantly, data have been validated, completed, corrected and standardized, a prerequisite of meaningful analysis. In addition, where critical data are lacking, such as genetic maps and RNA secondary structures, they are generated by the GOBASE team and collaborators, and added to the database. The database is implemented in a relational database management system, but features an object-oriented view of the biological data through a Web/Genera-generated World Wide Web interface. Finally, we have developed software for database curation (i.e. data updates, validation and correction), which will be described in some detail in this paper.

Formation of a GNRA tetraloop in P5abc can disrupt an interdomain interaction in the Tetrahymena group I ribozyme

The secondary structure of a truncated P5abc subdomain (tP5abc, a 56-nucleotide RNA) of the Tetrahymena thermophila group I intron ribozyme changes when its tertiary structure forms. We have now used heteronuclear NMR spectroscopy to determine its conformation in solution. The tP5abc RNA that contains only secondary structure is extended compared with the tertiary folded form; both forms coexist in slow chemical exchange (the interconversion rate constant is slower than 1 s−1) in the presence of magnesium. Kinetic experiments have shown that tertiary folding of the P5abc subdomain is one of the earliest folding transitions in the group I intron ribozyme, and that it leads to a metastable misfolded intermediate. Previous mutagenesis studies suggest that formation of the extended P5abc structure described here destabilize a misfolded intermediate. This study shows that the P5abc RNA subdomain containing a GNRA tetraloop in P5c (in contrast to the five-nucleotide loop P5c in the tertiary folded ribozyme) can disrupt the base-paired interdomain (P14) interaction between P5c and P2.

Targeting nucleic acid secondary structures by antisense oligonucleotides designed through in vitro selection.

Using an in vitro selection approach, we have isolated oligonucleotides that can bind to a DNA hairpin structure. Complex formation of these oligonucleotides with the target hairpin involves some type of triple-stranded structure with noncanonical interaction, as indicated by bandshift assays and footprinting studies. The selected oligomers can block restriction endonuclease cleavage of the target hairpin in a sequence-specific manner. We demonstrate that in vitro selection can extend the antisense approach to functional targeting of secondary structure motifs. This could provide a basis for interfering with regulatory processes mediated by a variety of nucleic acid structures.

NMR determination of the major solution conformation of a peptoid pentamer with chiral side chains

Polymers of N-substituted glycines (“peptoids”) containing chiral centers at the α position of their side chains can form stable structures in solution. We studied a prototypical peptoid, consisting of five para-substituted (S)-N-(1-phenylethyl)glycine residues, by NMR spectroscopy. Multiple configurational isomers were observed, but because of extensive signal overlap, only the major isomer containing all cis-amide bonds was examined in detail. The NMR data for this molecule, in conjunction with previous CD spectroscopic results, indicate that the major species in methanol is a right-handed helix with cis-amide bonds. The periodicity of the helix is three residues per turn, with a pitch of ≈6 Å. This conformation is similar to that anticipated by computational studies of a chiral peptoid octamer. The helical repeat orients the amide bond chromophores in a manner consistent with the intensity of the CD signal exhibited by this molecule. Many other chiral polypeptoids have similar CD spectra, suggesting that a whole family of peptoids containing chiral side chains is capable of adopting this secondary structure motif. Taken together, our experimental and theoretical studies of the structural properties of chiral peptoids lay the groundwork for the rational design of more complex polypeptoid molecules, with a variety of applications, ranging from nanostructures to nonviral gene delivery systems.

Isolation of a myoglobin molten globule by selective cobalt(III)-induced unfolding

Reaction of the Schiff-base complex [Co(acetylacetonate-ethylenediimine)(NH3)2]+ with metmyoglobin at pH 6.5 yields a partially folded protein containing six Co(III) complexes. Although half of its α-helical secondary structure is retained, absorption and CD spectra indicate that the tertiary structure in both B-F and AGH domains is disrupted in the partially folded protein. In analogy to proton-induced unfolding, it is likely that the loss of tertiary structure is triggered by metal-ion binding to histidines. Cobalt(III)-induced unfolding of myoglobin is unique in its selectivity (other proteins are unaffected) and in allowing the isolation of the partially folded macromolecule (the protein does not refold or aggregate upon removal of free denaturant).

Measures of residue density in protein structures

A hierarchy of residue density assessments and packing properties in protein structures are contrasted, including a regular density, a variety of charge densities, a hydrophobic density, a polar density, and an aromatic density. These densities are investigated by alternative distance measures and also at the interface of multiunit structures. Amino acids are divided into nine structural categories according to three secondary structure states and three solvent accessibility levels. To take account of amino acid abundance differences across protein structures, we normalize the observed density by the expected density defining a density index. Solvent accessibility levels exert the predominant influence in determinations of the regular residue density. Explicitly, the regular density values vary approximately linearly with respect to solvent accessibility levels, the linearity parameters depending on the amino acid. The charge index reveals pronounced inequalities between lysine and arginine in their interactions with acidic residues. The aromatic density calculations in all structural categories parallel the regular density calculations, indicating that the aromatic residues are distributed as a random sample of all residues. Moreover, aromatic residues are found to be over-represented in the neighborhood of all amino acids. This result might be attributed to nucleation sites and protein stability being substantially associated with aromatic residues.

Isolation and characterization of an amino acid-selective channel protein present in the chloroplastic outer envelope membrane

The reconstituted pea chloroplastic outer envelope protein of 16 kDa (OEP16) forms a slightly cation-selective, high-conductance channel with a conductance of Λ = 1,2 nS (in 1 M KCl). The open probability of OEP16 channel is highest at 0 mV (Popen = 0.8), decreasing exponentially with higher potentials. Transport studies using reconstituted recombinant OEP16 protein show that the OEP16 channel is selective for amino acids but excludes triosephosphates or uncharged sugars. Crosslinking indicates that OEP16 forms a homodimer in the membrane. According to its primary sequence and predicted secondary structure, OEP16 shows neither sequence nor structural homologies to classical porins. The results indicate that the intermembrane space between the two envelope membranes might not be as freely accessible as previously thought.

Exploring the folding free energy surface of a three-helix bundle protein

The multidimensional free energy surface for a small fast folding helical protein is explored based on first-principle calculations. The model represents the 46-residue segment from fragment B of staphylococcal protein A. The relationship between collapse and tertiary structure formation, and the order of collapse and secondary structure formation, are investigated. We find that the initial collapse process gives rise to a transition state with about 30% of the native tertiary structure and 50–70% of the native helix content. We also observe two distinct distributions of native helix in this collapsed state (Rg ≈ 12 Å), one with about 20% of the native helical hydrogen bonds, the other with near 70%. The former corresponds to a local minimum. The barrier from this metastable state to the native state is about 2 kBT. In the latter case, folding is essentially a downhill process involving topological assembly. In addition, the order of formation of secondary structure among the three helices is examined. We observe cooperative formation of the secondary structure in helix I and helix II. Secondary structure in helix III starts to form following the formation of certain secondary structure in both helix I and helix II. Comparisons of our results with those from theory and experiment are made.

Molecular and functional characterization of a novel low-affinity cation transporter (LCT1) in higher plants

The transport of cations across membranes in higher plants plays an essential role in many physiological processes including mineral nutrition, cell expansion, and the transduction of environmental signals. In higher plants the coordinated expression of transport mechanisms is essential for specialized cellular processes and for adaptation to variable environmental conditions. To understand the molecular basis of cation transport in plant roots, a Triticum aestivum cDNA library was used to complement a yeast mutant deficient in potassium (K+) uptake. Two genes were cloned that complemented the mutant: HKT1 and a novel cDNA described in this report encoding a cation transporter, LCT1 (low-affinity cation transporter). Analysis of the secondary structure of LCT1 suggests that the protein contains 8–10 transmembrane helices and a hydrophilic amino terminus containing sequences enriched in Pro, Ser, Thr, and Glu (PEST). The transporter activity was assayed using radioactive isotopes in yeast cells expressing the cDNA. LCT1 mediated low-affinity uptake of the cations Rb+ and Na+, and possibly allowed Ca2+ but not Zn2+ uptake. LCT1 is expressed in low abundance in wheat roots and leaves. The precise functional role of this cation transporter is not known, although the competitive inhibition of cation uptake by Ca2+ has parallels to whole plant and molecular studies that have shown the important role of Ca2+ in reducing Na+ uptake and ameliorating Na+ toxicity. The structure of this higher plant ion transport protein is unique and contains PEST sequences.

A structural census of the current population of protein sequences

We examine the occurrence of the ≈300 known protein folds in different groups of organisms. To do this, we characterize a large fraction of the currently known protein sequences (≈140,000) in structural terms, by matching them to known structures via sequence comparison (or by secondary-structure class prediction for those without structural homologues). Overall, we find that an appreciable fraction of the known folds are present in each of the major groups of organisms (e.g., bacteria and eukaryotes share 156 of 275 folds), and most of the common folds are associated with many families of nonhomologous sequences (i.e., >10 sequence families for each common fold). However, different groups of organisms have characteristically distinct distributions of folds. So, for instance, some of the most common folds in vertebrates, such as globins or zinc fingers, are rare or absent in bacteria. Many of these differences in fold usage are biologically reasonable, such as the folds of metabolic enzymes being common in bacteria and those associated with extracellular transport and communication being common in animals. They also have important implications for database-based methods for fold recognition, suggesting that an unknown sequence from a plant is more likely to have a certain fold (e.g., a TIM barrel) than an unknown sequence from an animal.