HyPaLib: a database of RNAs and RNA structural elements defined by hybrid patterns

The database, called HyPaLib (for Hybrid Pattern Library), contains annotated structural elements characteristic for certain classes of structural and/or functional RNAs. These elements are described in a language specifically designed for this purpose. The language allows convenient specification of hybrid patterns, i.e. motifs consisting of sequence features and structural elements together with sequence similarity and thermodynamic constraints. We are currently developing software tools that allow a user to search sequence databases for any pattern in HyPaLib, thus providing functionality which is similar to PROSITE, but dedicated to the more complex patterns in RNA sequences. HyPaLib is available at http://bibiserv.techfak.uni-bielefeld.de/HyPa/.

Comparative analysis of ribonuclease P RNA using gene sequences from natural microbial populations reveals tertiary structural elements.

PCR amplification of template DNAs extracted from mixed, naturally occurring microbial populations, using oligonucleotide primers complementary to highly conserved sequences, was used to obtain a large collection of diverse RNase P RNA-encoding genes. An alignment of these sequences was used in a comparative analysis of RNase P RNA secondary and tertiary structure. The new sequences confirm the secondary structure model based on sequences from cultivated organisms (with minor alterations in helices P12 and P18), providing additional support for nearly every base pair. Analysis of sequence covariation using the entire RNase P RNA data set reveals elements of tertiary structure in the RNA; the third nucleotides (underlined) of the GNRA tetraloops L14 and L18 are seen to interact with adjacent Watson-Crick base pairs in helix P8, forming A:G/C or G:A/U base triples. These experiments demonstrate one way in which the enormous diversity of natural microbial populations can be used to elucidate molecular structure through comparative analysis.

Structural mimicry of a native protein by a minimized binding domain

The affinity between molecules depends both on the nature and presentation of the contacts. Here, we observe coupling of functional and structural elements when a protein binding domain is evolved to a smaller functional mimic. Previously, a 38-residue form of the 59-residue B-domain of protein A, termed Z38, was selected by phage display. Z38 contains 13 mutations and binds IgG only 10-fold weaker than the native B-domain. We present the solution structure of Z38 and show that it adopts a tertiary structure remarkably similar to that observed for the first two helices of B-domain in the B-domain/Fc complex [Deisenhofer, J. (1981) Biochemistry 20, 2361–2370], although it is significantly less stable. Based on this structure, we have improved on Z38 by designing a 34-residue disulfide-bonded variant (Z34C) that has dramatically enhanced stability and binds IgG with 9-fold higher affinity. The improved stability of Z34C led to NMR spectra with much greater chemical shift dispersion, resulting in a more precisely determined structure. Z34C, like Z38, has a structure virtually identical to the equivalent region from native protein A domains. The well-defined hydrophobic core of Z34C reveals key structural features that have evolved in this small, functional domain. Thus, the stabilized two-helix peptide, about half the size and having one-third of the remaining residues altered, accurately mimics both the structure and function of the native domain.

Molecular recognition of angiogenesis inhibitors fumagillin and ovalicin by methionine aminopeptidase 2

Angiogenesis inhibitors are a novel class of promising therapeutic agents for treating cancer and other human diseases. Fumagillin and ovalicin compose a class of structurally related natural products that potently inhibit angiogenesis by blocking endothelial cell proliferation. A synthetic analog of fumagillin, TNP-470, is currently undergoing clinical trials for treatment of a variety of cancers. A common target for fumagillin and ovalicin recently was identified as the type 2 methionine aminopeptidase (MetAP2). These natural products bind MetAP2 covalently, inhibiting its enzymatic activity. The specificity of this binding is underscored by the lack of inhibition of the closely related type 1 enzyme, MetAP1. The molecular basis of the high affinity and specificity of these inhibitors for MetAP2 has remained undiscovered. To determine the structural elements of these inhibitors and MetAP2 that are involved in this interaction, we synthesized fumagillin analogs in which each of the potentially reactive epoxide groups was removed either individually or in combination. We found that the ring epoxide in fumagillin is involved in the covalent modification of MetAP2, whereas the side chain epoxide group is dispensable. By using a fumagillin analog tagged with fluorescein, His-231 in MetAP2 was identified as the residue that is covalently modified by fumagillin. Site-directed mutagenesis of His-231 demonstrated its importance for the catalytic activity of MetAP2 and confirmed that the same residue is covalently modified by fumagillin. These results, in agreement with a recent structural study, suggest that fumagillin and ovalicin inhibit MetAP2 by irreversible blockage of the active site.

Identification of specific Rp-phosphate oxygens in the tRNA anticodon loop required for ribosomal P-site binding

tRNA binding to the ribosomal P site is dependent not only on correct codon–anticodon interaction but also involves identification of structural elements of tRNA by the ribosome. By using a phosphorothioate substitution–interference approach, we identified specific nonbridging Rp-phosphate oxygens in the anticodon loop of tRNAPhe from Escherichia coli which are required for P-site binding. Stereo-specific involvement of phosphate oxygens at these positions was confirmed by using synthetic anticodon arm analogues at which single Rp- or Sp-phosphorothioates were incorporated. Identical interference results with yeast tRNAPhe and E. coli tRNAfMet indicate a common backbone conformation or common recognition elements in the anticodon loop of tRNAs. N-ethyl-N-nitrosourea modification–interference experiments with natural tRNAs point to the importance of the same phosphates in the loop. Guided by the crystal structure of tRNAPhe, we propose that specific Rp-phosphate oxygens are required for anticodon loop (“U-turn”) stabilization or are involved in interactions with the ribosome on correct tRNA–mRNA complex formation.

Crystal structure of the ribosomal RNA domain essential for binding elongation factors

The structure of a 29-nucleotide RNA containing the sarcin/ricin loop (SRL) of rat 28 S rRNA has been determined at 2.1 Å resolution. Recognition of the SRL by elongation factors and by the ribotoxins, sarcin and ricin, requires a nearly universal dodecamer sequence that folds into a G-bulged cross-strand A stack and a GAGA tetraloop. The juxtaposition of these two motifs forms a distorted hairpin structure that allows direct recognition of bases in both grooves as well as recognition of nonhelical backbone geometry and two 5′-unstacked purines. Comparisons with other RNA crystal structures establish the cross-strand A stack and the GNRA tetraloop as defined and modular RNA structural elements. The conserved region at the top is connected to the base of the domain by a region presumed to be flexible because of the sparsity of stabilizing contacts. Although the conformation of the SRL RNA previously determined by NMR spectroscopy is similar to the structure determined by x-ray crystallography, significant differences are observed in the “flexible” region and to a lesser extent in the G-bulged cross-strand A stack.

Intron positions correlate with module boundaries in ancient proteins

We analyze the three-dimensional structure of proteins by a computer program that finds regions of sequence that contain module boundaries, defining a module as a segment of polypeptide chain bounded in space by a specific given distance. The program defines a set of “linker regions” that have the property that if an intron were to be placed into each linker region, the protein would be dissected into a set of modules all less than the specified diameter. We test a set of 32 proteins, all of ancient origin, and a corresponding set of 570 intron positions, to ask if there is a statistically significant excess of intron positions within the linker regions. For 28-Å modules, a standard size used historically, we find such an excess, with P < 0.003. This correlation is neither due to a compositional or sequence bias in the linker regions nor to a surface bias in intron positions. Furthermore, a subset of 20 introns, which can be putatively identified as old, lies even more explicitly within the linker regions, with P < 0.0003. Thus, there is a strong correlation between intron positions and three-dimensional structural elements of ancient proteins as expected by the introns-early approach. We then study a range of module diameters and show that, as the diameter varies, significant peaks of correlation appear for module diameters centered at 21.7, 27.6, and 32.9 Å. These preferred module diameters roughly correspond to predicted exon sizes of 15, 22, and 30 residues. Thus, there are significant correlations between introns, modules, and a quantized pattern of the lengths of polypeptide chains, which is the prediction of the “Exon Theory of Genes.”

Structure of neurolysin reveals a deep channel that limits substrate access

The zinc metallopeptidase neurolysin is shown by x-ray crystallography to have large structural elements erected over the active site region that allow substrate access only through a deep narrow channel. This architecture accounts for specialization of this neuropeptidase to small bioactive peptide substrates without bulky secondary and tertiary structures. In addition, modeling studies indicate that the length of a substrate N-terminal to the site of hydrolysis is restricted to approximately 10 residues by the limited size of the active site cavity. Some structural elements of neurolysin, including a five-stranded β-sheet and the two active site helices, are conserved with other metallopeptidases. The connecting loop regions of these elements, however, are much extended in neurolysin, and they, together with other open coil elements, line the active site cavity. These potentially flexible elements may account for the ability of the enzyme to cleave a variety of sequences.

Similarities between the DNA replication initiators of Gram-negative bacteria plasmids (RepA) and eukaryotes (Orc4p)/archaea (Cdc6p)

The proteins responsible for the initiation of DNA replication are thought to be essentially unrelated in bacteria and archaea/eukaryotes. Here we show that RepA, the initiator from the Pseudomonas plasmid pPS10, and the C-terminal domain of ScOrc4p, a subunit of Saccharomyces cerevisiae (Sc) origin recognition complex (ORC), share sequence similarities. Based on biochemical and spectroscopic evidence, these similarities include common structural elements, such as a winged-helix domain and a leucine-zipper dimerization motif. We have also found that ScOrc4p, as previously described for RepA-type initiators, interacts with chaperones of the Hsp70 family both in vitro and in vivo, most probably to regulate the assembly of active ORC. In evolutionary terms, our results are compatible with the recruitment of the same protein module for initiation of DNA replication by the ancestors of present-day Gram-negative bacteria plasmids, archaea, and eukaryotes.

Structure of melanoma inhibitory activity protein, a member of a recently identified family of secreted proteins

Melanoma inhibitory activity (MIA) is a 12-kDa protein that is secreted from both chondrocytes and malignant melanoma cells. MIA has been reported to have effects on cell growth and adhesion, and it may play a role in melanoma metastasis and cartilage development. We report the 1.4-Å crystal structure of human MIA, which consists of an Src homology 3 (SH3)-like domain with N- and C-terminal extensions of about 20 aa each. The N- and C-terminal extensions add additional structural elements to the SH3 domain, forming a previously undescribed fold. MIA is a representative of a recently identified family of proteins and is the first structure of a secreted protein with an SH3 subdomain. The structure also suggests a likely protein interaction site and suggests that, unlike conventional SH3 domains, MIA does not recognize polyproline helices.

Electrostatic steering and ionic tethering in enzyme–ligand binding: Insights from simulations

To bind at an enzyme’s active site, a ligand must diffuse or be transported to the enzyme’s surface, and, if the binding site is buried, the ligand must diffuse through the protein to reach it. Although the driving force for ligand binding is often ascribed to the hydrophobic effect, electrostatic interactions also influence the binding process of both charged and nonpolar ligands. First, electrostatic steering of charged substrates into enzyme active sites is discussed. This is of particular relevance for diffusion-influenced enzymes. By comparing the results of Brownian dynamics simulations and electrostatic potential similarity analysis for triose-phosphate isomerases, superoxide dismutases, and β-lactamases from different species, we identify the conserved features responsible for the electrostatic substrate-steering fields. The conserved potentials are localized at the active sites and are the primary determinants of the bimolecular association rates. Then we focus on a more subtle effect, which we will refer to as “ionic tethering.” We explore, by means of molecular and Brownian dynamics simulations and electrostatic continuum calculations, how salt links can act as tethers between structural elements of an enzyme that undergo conformational change upon substrate binding, and thereby regulate or modulate substrate binding. This is illustrated for the lipase and cytochrome P450 enzymes. Ionic tethering can provide a control mechanism for substrate binding that is sensitive to the electrostatic properties of the enzyme’s surroundings even when the substrate is nonpolar.

Electrospray ionization mass spectrometry analysis of changes in phospholipids in RBL-2H3 mastocytoma cells during degranulation

Biological membranes contain an extraordinary diversity of lipids. Phospholipids function as major structural elements of cellular membranes, and analysis of changes in the highly heterogeneous mixtures of lipids found in eukaryotic cells is central to understanding the complex functions in which lipids participate. Phospholipase-catalyzed hydrolysis of phospholipids often follows cell surface receptor activation. Recently, we demonstrated that granule fusion is initiated by addition of exogenous, nonmammalian phospholipases to permeabilized mast cells. To pursue this finding, we use positive and negative mode Fourier-transform ion cyclotron resonance mass spectrometry (FTICR-MS) to measure changes in the glycerophospholipid composition of total lipid extracts of intact and permeabilized RBL-2H3 (mucosal mast cell line) cells. The low energy of the electrospray ionization results in efficient production of molecular ions of phospholipids uncomplicated by further fragmentation, and changes were observed that eluded conventional detection methods. From these analyses we have spectrally resolved more than 130 glycerophospholipids and determined changes initiated by introduction of exogenous phospholipase C, phospholipase D, or phospholipase A2. These exogenous phospholipases have a preference for phosphatidylcholine with long polyunsaturated alkyl chains as substrates and, when added to permeabilized mast cells, produce multiple species of mono- and polyunsaturated diacylglycerols, phosphatidic acids, and lysophosphatidylcholines, respectively. The patterns of changes of these lipids provide an extraordinarily rich source of data for evaluating the effects of specific lipid species generated during cellular processes, such as exocytosis.

Mechanical behavior in living cells consistent with the tensegrity model

Alternative models of cell mechanics depict the living cell as a simple mechanical continuum, porous filament gel, tensed cortical membrane, or tensegrity network that maintains a stabilizing prestress through incorporation of discrete structural elements that bear compression. Real-time microscopic analysis of cells containing GFP-labeled microtubules and associated mitochondria revealed that living cells behave like discrete structures composed of an interconnected network of actin microfilaments and microtubules when mechanical stresses are applied to cell surface integrin receptors. Quantitation of cell tractional forces and cellular prestress by using traction force microscopy confirmed that microtubules bear compression and are responsible for a significant portion of the cytoskeletal prestress that determines cell shape stability under conditions in which myosin light chain phosphorylation and intracellular calcium remained unchanged. Quantitative measurements of both static and dynamic mechanical behaviors in cells also were consistent with specific a priori predictions of the tensegrity model. These findings suggest that tensegrity represents a unified model of cell mechanics that may help to explain how mechanical behaviors emerge through collective interactions among different cytoskeletal filaments and extracellular adhesions in living cells.

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.

Global properties of the mapping between local amino acid sequence and local structure in proteins.

Local protein structure prediction efforts have consistently failed to exceed approximately 70% accuracy. We characterize the degeneracy of the mapping from local sequence to local structure responsible for this failure by investigating the extent to which similar sequence segments found in different proteins adopt similar three-dimensional structures. Sequence segments 3-15 residues in length from 154 different protein families are partitioned into neighborhoods containing segments with similar sequences using cluster analysis. The consistency of the sequence-to-structure mapping is assessed by comparing the local structures adopted by sequence segments in the same neighborhood in proteins of known structure. In the 154 families, 45% and 28% of the positions occur in neighborhoods in which one and two local structures predominate, respectively. The sequence patterns that characterize the neighborhoods in the first class probably include virtually all of the short sequence motifs in proteins that consistently occur in a particular local structure. These patterns, many of which occur in transitions between secondary structural elements, are an interesting combination of previously studied and novel motifs. The identification of sequence patterns that consistently occur in one or a small number of local structures in proteins should contribute to the prediction of protein structure from sequence.