On the existence and implications of an inverse folding code in proteins.

The existence of a code relating the set of possible sequences at a given position in a protein backbone to the local structure at that location is investigated. It is shown that only 73% of 4-C alpha structure fragments in a sample of 114 protein structures exhibit a preference for a particular set of sequences. The remaining structures can accommodate essentially any sequence. The structures that encode specific sequence distributions include the classical "secondary" structures, with the notable exception of planar (beta) bends. It is suggested that this has implications as to the mechanism of folding in proteins with extensive sheet/barrel structure. The possible role of structures that do not encode specific sequences as mutation hot spots is noted.

Atom density in protein structures

The residue environment in protein structures is studied with respect to the density of carbon (C), oxygen (O), and nitrogen (N) atoms within a certain distance (say 5 Å) of each residue. Two types of environments are evaluated: one based on side-chain atom contacts (abbreviated S-S) and the other based on all atom (side-chain + backbone) contacts (abbreviated A-A). Different atom counts are observed about nine-residue structural categories defined by three solvent accessibility levels and three secondary structure states. Among the structural categories, the S-S atom count ratios generally vary more than the A-A atom count ratios because of the fact that the backbone (O) and (N) atoms contribute equal counts. Secondary structure affects the (C) density for the A-A contacts whereas secondary structure has little influence on the (C) density for the S-S contacts. For S-S contacts, a greater density of (O) over (N) atom neighbors stands out in the environment of most amino acid types. By contrast, for A-A contacts, independent of the solvent accessibility levels, the ratio (O)/(N) is ≈1 in helical states, consistent with the geometry of α-helical residues whose side-chains tilt oppositely to the amino to carboxy α-helical axis. The highest ratio of neighbor (O)/(N) is achieved under solvent exposed conditions. This (O) vs. (N) prevalence is advantageous at the protein surface that generally exhibits an acid excess that helps to enhance protein solubility in the cell and to avoid nonspecific interactions with phosphate groups of DNA, RNA, and other plasma constituents.

A trans-acting RNA as a control switch in Escherichia coli: DsrA modulates function by forming alternative structures

DsrA is an 87-nucleotide regulatory RNA of Escherichia coli that acts in trans by RNA–RNA interactions with two different mRNAs, hns and rpoS. DsrA has opposite effects on these transcriptional regulators. H-NS levels decrease, whereas RpoS (σs) levels increase. Here we show that DsrA enhances hns mRNA turnover yet stabilizes rpoS mRNA, either directly or via effects on translation. Computational and RNA footprinting approaches led to a refined structure for DsrA, and a model in which DsrA interacts with the hns mRNA start and stop codon regions to form a coaxial stack. Analogous bipartite interactions exist in eukaryotes, albeit with different regulatory consequences. In contrast, DsrA base pairs in discrete fashion with the rpoS RNA translational operator. Thus, different structural configurations for DsrA lead to opposite regulatory consequences for target RNAs.

Association of RNase mitochondrial RNA processing enzyme with ribonuclease P in higher ordered structures in the nucleolus: a possible coordinate role in ribosome biogenesis.

RNase mitochondrial RNA processing enzyme (MRP) is a nucleolar ribonucleoprotein particle that participates in 5.8S ribosomal RNA maturation in eukaryotes. This enzyme shares a polypeptide and an RNA structural motif with ribonuclease P (RNase P), a nuclear endoribonuclease originally described in the nucleus that processes RNA transcripts to generate their mature 5' termini. Both enzymes are also located in mitochondria. This report further characterizes the relationship between RNase MRP and RNase P. Antisense affinity selection with biotinylated 2'-O-methyl oligoribonucleotides and glycerol gradient fractionation experiments demonstrated that small subpopulations of RNase MRP and RNase P associate with each other in vivo in macromolecular complex, possibly 60-80S preribosomes. This latter notion was supported by fluorescence in situ hybridization experiments with antisense oligonucleotides that localized that RNA components of RNase MRP and RNase P to the nucleolus and to discrete cytoplasmic structures. These findings suggest that small subpopulations of RNase MRP and RNase P are physically associated, and that both may function in ribosomal RNA maturation or ribosome assembly.

Direct evidence for secondary loss of mitochondria in Entamoeba histolytica.

Archezoan protists are though to represent lineages that diverged from other eukaryotes before acquisition of the mitochondrion and other organelles. The parasite Entamoeba histolytica was originally included in this group. Ribosomal RNA based phylogenies, however, place E. histolytica on a comparatively recent branch of the eukaryotic tree, implying that its ancestors had these structures. In this study, direct evidence for secondary loss of mitochondrial function was obtained by isolating two E. histolytica genes encoding proteins that in other eukaryotes are localized in the mitochondrion: the enzyme pyridine nucleotide transhydrogenase and the chaperonin cpn60. Phylogenetic analysis of the E. histolytica homolog of cpn60 confirmed that it is specifically related to the mitochondrial lineage. The data suggest that a mitochondrial relic may persist in this organism. Similar studies are needed in archezoan protists to ascertain which, if any, eukaryotic lineages primitively lack mitochondria.

A secondary-structure model for the self-cleaving region of Neurospora VS RNA.

Neurospora VS RNA performs an RNA-mediated self-cleavage reaction whose products contain 2',3'-cyclic phosphate and 5'-hydroxyl termini. This reaction is similar to those of hammerhead, hairpin, and hepatitis delta virus ribozymes; however, VS RNA is not similar in sequence to these other self-cleaving motifs. Here we propose a model for the secondary structure of the self-cleaving region of VS RNA, supported by site-directed mutagenesis and chemical modification structure probing data. The secondary structure of VS RNA is distinct from those of the other naturally occurring RNA self-cleaving domains. In addition to a unique secondary structure, several Mg-dependent interactions occur during the folding of VS RNA into its active tertiary conformation.

Sequence-specific polypeptoids: A diverse family of heteropolymers with stable secondary structure

We have synthesized and characterized a family of structured oligo-N-substituted-glycines (peptoids) up to 36 residues in length by using an efficient solid-phase protocol to incorporate chemically diverse side chains in a sequence-specific fashion. We investigated polypeptoids containing side chains with a chiral center adjacent to the main chain nitrogen. Some of these sequences have stable secondary structure, despite the achirality of the polymer backbone and its lack of hydrogen bond donors. In both aqueous and organic solvents, peptoid oligomers as short as five residues give rise to CD spectra that strongly resemble those of peptide α-helices. Differential scanning calorimetry and CD measurements show that polypeptoid secondary structure is highly stable and that unfolding is reversible and cooperative. Thermodynamic parameters obtained for unfolding are similar to those obtained for the α-helix to coil transitions of peptides. This class of biomimetic polymers may enable the design of self-assembling macromolecules with novel structures and functions.

Deleterious mutations destabilize ribosomal RNA in endosymbiotic bacteria

In populations that are small and asexual, mutations with slight negative effects on fitness will drift to fixation more often than in large or sexual populations in which they will be eliminated by selection. If such mutations occur in substantial numbers, the combined effects of long-term asexuality and small population size may result in substantial accumulation of mildly deleterious substitutions. Prokaryotic endosymbionts of animals that are transmitted maternally for very long periods are effectively asexual and experience smaller effective population size than their free-living relatives. The contrast between such endosymbionts and related free-living bacteria allows us to test whether a population structure imposing frequent bottlenecks and asexuality does lead to an accumulation of slightly deleterious substitutions. Here we show that several independently derived insect endosymbionts, each with a long history of maternal transmission, have accumulated destabilizing base substitutions in the highly conserved 16S rRNA. Stabilities of Domain I of this subunit are 15–25% lower in endosymbionts than in closely related free-living bacteria. By mapping destabilizing substitutions onto a reconstructed phylogeny, we show that decreased ribosomal stability has evolved separately in each endosymbiont lineage. Our phylogenetic approach allows us to demonstrate statistical significance for this pattern: becoming endosymbiotic predictably results in decreased stability of rRNA secondary structure.

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.

RNA folding kinetics regulates translation of phage MS2 maturation gene

The gene for the maturation protein of the single-stranded RNA coliphage MS2 is preceded by an untranslated leader of 130 nt, which folds into a cloverleaf, i.e., three stem–loop structures enclosed by a long distance interaction (LDI). This LDI prevents translation because its 3′ moiety contains the Shine–Dalgarno sequence of the maturation gene. Previously, several observations suggested that folding of the cloverleaf is kinetically delayed, providing a time window for ribosomes to access the RNA. Here we present direct evidence for this model. In vitro experiments show that ribosome binding to the maturation gene is faster than refolding of the denatured cloverleaf. This folding delay appears related to special properties of the leader sequence. We have replaced the three stem–loop structures by a single five nt loop. This change does not affect the equilibrium structure of the LDI. Nevertheless, in this construct, the folding delay has virtually disappeared, suggesting that now the RNA folds faster than ribosomes can bind. Perturbation of the cloverleaf by an insertion makes the maturation start permanently accessible. A pseudorevertant that evolved from an infectious clone carrying the insertion had overcome this defect. It showed a wild-type folding delay before closing down the maturation gene. This experiment reveals the biological significance of retarded cloverleaf formation.

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.

Subtle atomic group discrimination in the RNA minor groove

As a problem in molecular recognition and for drug discovery, great interest has developed around the possibility that RNA structures could be discriminated by peptides and other small molecules. Although small peptides have been shown to have the capacity to discriminate specific bulges and loops in RNA molecules, discrimination of double helical regions by a peptide binder has not been reported. Indeed, the most accessible part of an RNA helix is the minor groove, and fundamental stereochemical considerations have suggested that discrimination of at least some base pairs would be difficult in the minor groove. Here we report the design and isolation of a peptide binder that manifests the most subtle kind of discrimination of base pair differences in the RNA minor groove. Functional discrimination of a single atomic group is demonstrated as well as the difference between two different angular orientations of the same group. This report of RNA helix discrimination by a peptide binder suggests a richer potential for RNA minor groove recognition than previously thought.

Molding a peptide into an RNA site by in vivo peptide evolution

Short peptides corresponding to the arginine-rich domains of several RNA-binding proteins are able to bind to their specific RNA sites with high affinities and specificities. In the case of the HIV-1 Rev-Rev response element (RRE) complex, the peptide forms a single α-helix that binds deeply in a widened, distorted RNA major groove and makes a substantial set of base-specific and backbone contacts. Using a reporter system based on antitermination by the bacteriophage λ N protein, it has been possible to identify novel arginine-rich peptides from combinatorial libraries that recognize the RRE with affinities and specificities similar to Rev but that appear to bind in nonhelical conformations. Here we have used codon-based mutagenesis to evolve one of these peptides, RSG-1, into an even tighter binder. After two rounds of evolution, RSG-1.2 bound the RRE with 7-fold higher affinity and 15-fold higher specificity than the wild-type Rev peptide, and in vitro competition experiments show that RSG-1.2 completely displaces the intact Rev protein from the RRE at low peptide concentrations. By fusing RRE-binding peptides to the activation domain of HIV-1 Tat, we show that the peptides can deliver Tat to the RRE site and activate transcription in mammalian cells, and more importantly, that the fusion proteins can inhibit the activity of Rev in chloramphenicol acetyltransferase reporter assays. The evolved peptides contain proline and glutamic acid mutations near the middle of their sequences and, despite the presence of a proline, show partial α-helix formation in the absence of RNA. These directed evolution experiments illustrate how readily complex peptide structures can be evolved within the context of an RNA framework, perhaps reflecting how early protein structures evolved in an “RNA world.”

Crystal structure of the RNA-dependent RNA polymerase of hepatitis C virus

We report the crystal structure of the RNA-dependent RNA polymerase of hepatitis C virus, a major human pathogen, to 2.8-Å resolution. This enzyme is a key target for developing specific antiviral therapy. The structure of the catalytic domain contains 531 residues folded in the characteristic fingers, palm, and thumb subdomains. The fingers subdomain contains a region, the “fingertips,” that shares the same fold with reverse transcriptases. Superposition to the available structures of the latter shows that residues from the palm and fingertips are structurally equivalent. In addition, it shows that the hepatitis C virus polymerase was crystallized in a closed fingers conformation, similar to HIV-1 reverse transcriptase in ternary complex with DNA and dTTP [Huang H., Chopra, R., Verdine, G. L. & Harrison, S. C. (1998) Science 282, 1669–1675]. This superposition reveals the majority of the amino acid residues of the hepatitis C virus enzyme that are likely to be implicated in binding to the replicating RNA molecule and to the incoming NTP. It also suggests a rearrangement of the thumb domain as well as a possible concerted movement of thumb and fingertips during translocation of the RNA template-primer in successive polymerization rounds.