Redesigning secondary structure to invert coenzyme specificity in isopropylmalate dehydrogenase.

Rational engineering of enzymes involves introducing key amino acids guided by a knowledge of protein structure to effect a desirable change in function. To date, all successful attempts to change specificity have been limited to substituting individual amino acids within a protein fold. However, the infant field of protein engineering will only reach maturity when changes in function can be generated by rationally engineering secondary structures. Guided by x-ray crystal structures and molecular modeling, site-directed mutagenesis has been used to systematically invert the coenzyme specificity of Thermus thermophilus isopropylmalate dehydrogenase from a 100-fold preference for NAD to a 1000-fold preference for NADP. The engineered mutant, which is twice as active as wild type, contains four amino acid substitutions and an alpha-helix and loop that replaces the original beta-turn. These results demonstrate that rational engineering of secondary structures to produce enzymes with novel properties is feasible.

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.

Determining RNA solution structure by segmental isotopic labeling and NMR: application to Caenorhabditis elegans spliced leader RNA 1.

Recent developments in multidimensional heteronuclear NMR spectroscopy and large-scale synthesis of uniformly 13C- and 15N-labeled oligonucleotides have greatly improved the prospects for determination of the solution structure of RNA. However, there are circumstances in which it may be advantageous to label only a segment of the entire RNA chain. For example, in a larger RNA molecule the structural question of interest may reside in a localized domain. Labeling only the corresponding nucleotides simplifies the spectrum and resonance assignments because one can filter proton spectra for coupling to 13C and 15N. Another example is in resolving alternative secondary structure models that are indistinguishable in imino proton connectivities. Here we report a general method for enzymatic synthesis of quantities of segmentally labeled RNA molecules required for NMR spectroscopy. We use the method to distinguish definitively two competing secondary structure models for the 5' half of Caenorhabditis elegans spliced leader RNA by comparison of the two-dimensional [15N] 1H heteronuclear multiple quantum correlation spectrum of the uniformly labeled sample with that of a segmentally labeled sample. The method requires relatively small samples; solutions in the 200-300 microM concentration range, with a total of 30 nmol or approximately 40 micrograms of RNA in approximately 150 microliters, give strong NMR signals in a short accumulation time. The method can be adapted to label an internal segment of a larger RNA chain for study of localized structural problems. This definitive approach provides an alternative to the more common enzymatic and chemical footprinting methods for determination of RNA secondary structure.

Periodicity of polar and nonpolar amino acids is the major determinant of secondary structure in self-assembling oligomeric peptides.

The tendency of a polypeptide chain to form alpha-helical or beta-strand secondary structure depends upon local and nonlocal effects. Local effects reflect the intrinsic propensities of the amino acid residues for particular secondary structures, while nonlocal effects reflect the positioning of the individual residues in the context of the entire amino acid sequence. In particular, the periodicity of polar and nonpolar residues specifies whether a given sequence is consistent with amphiphilic alpha-helices or beta-strands. The importance of intrinsic propensities was compared to that of polar/nonpolar periodicity by a direct competition. Synthetic peptides were designed using residues with intrinsic propensities that favored one or the other type of secondary structure. The polar/nonpolar periodicities of the peptides were designed either to be consistent with the secondary structure favored by the intrinsic propensities of the component residues or in other cases to oppose these intrinsic propensities. Characterization of the synthetic peptides demonstrated that in all cases the observed secondary structure correlates with the periodicity of the peptide sequence--even when this secondary structure differs from that predicted from the intrinsic propensities of the component amino acids. The observed secondary structures are concentration dependent, indicating that oligomerization of the amphiphilic peptides is responsible for the observed secondary structures. Thus, for self-assembling oligomeric peptides, the polar/nonpolar periodicity can overwhelm the intrinsic propensities of the amino acid residues and serves as the major determinant of peptide secondary structure.

Operational RNA code for amino acids: species-specific aminoacylation of minihelices switched by a single nucleotide.

The genetic code is based on aminoacylation reactions where specific amino acids are attached to tRNAs bearing anticodon trinucleotides. However, the anticodon-independent specific aminoacylation of RNA minihelix substrates by bacterial and yeast tRNA synthetases suggested an operational RNA code for amino acids whereby specific RNA sequences/structures in tRNA acceptor stems correspond to specific amino acids. Because of the possible significance of the operational RNA code for the development of the genetic code, we investigated aminoacylation of synthetic RNA minihelices with a human enzyme to understand the sequences needed for that aminoacylation compared with those needed for a microbial system. We show here that the species-specific aminoacylation of glycine tRNAs is recapitulated by a species-specific aminoacylation of minihelices. Although the mammalian and Escherichia coli minihelices differ at 6 of 12 base pairs, two of the three nucleotides essential for aminoacylation by the E. coli enzyme are conserved in the mammalian minihelix. The two conserved nucleotides were shown to be also important for aminoacylation of the mammalian minihelix by the human enzyme. A simple interchange of the differing nucleotide enabled the human enzyme to now charge the bacterial substrate and not the mammalian minihelix. Conversely, this interchange made the bacterial enzyme specific for the mammalian substrate. Thus, the positional locations (if not the actual nucleotides) for the operational RNA code for glycine appear conserved from bacteria to mammals.

Structure-based conformational preferences of amino acids

Proteins can be very tolerant to amino acid substitution, even within their core. Understanding the factors responsible for this behavior is of critical importance for protein engineering and design. Mutations in proteins have been quantified in terms of the changes in stability they induce. For example, guest residues in specific secondary structures have been used as probes of conformational preferences of amino acids, yielding propensity scales. Predicting these amino acid propensities would be a good test of any new potential energy functions used to mimic protein stability. We have recently developed a protein design procedure that optimizes whole sequences for a given target conformation based on the knowledge of the template backbone and on a semiempirical potential energy function. This energy function is purely physical, including steric interactions based on a Lennard-Jones potential, electrostatics based on a Coulomb potential, and hydrophobicity in the form of an environment free energy based on accessible surface area and interatomic contact areas. Sequences designed by this procedure for 10 different proteins were analyzed to extract conformational preferences for amino acids. The resulting structure-based propensity scales show significant agreements with experimental propensity scale values, both for α-helices and β-sheets. These results indicate that amino acid conformational preferences are a natural consequence of the potential energy we use. This confirms the accuracy of our potential and indicates that such preferences should not be added as a design criterion.

Intron self-complementarity enforces exon inclusion in a yeast pre-mRNA

Skipping of internal exons during removal of introns from pre-mRNA must be avoided for proper expression of most eukaryotic genes. Despite significant understanding of the mechanics of intron removal, mechanisms that ensure inclusion of internal exons in multi-intron pre-mRNAs remain mysterious. Using a natural two-intron yeast gene, we have identified distinct RNA–RNA complementarities within each intron that prevent exon skipping and ensure inclusion of internal exons. We show that these complementarities are positioned to act as intron identity elements, bringing together only the appropriate 5′ splice sites and branchpoints. Destroying either intron self-complementarity allows exon skipping to occur, and restoring the complementarity using compensatory mutations rescues exon inclusion, indicating that the elements act through formation of RNA secondary structure. Introducing new pairing potential between regions near the 5′ splice site of intron 1 and the branchpoint of intron 2 dramatically enhances exon skipping. Similar elements identified in single intron yeast genes contribute to splicing efficiency. Our results illustrate how intron secondary structure serves to coordinate splice site pairing and enforce exon inclusion. We suggest that similar elements in vertebrate genes could assist in the splicing of very large introns and in the evolution of alternative splicing.

Molecular Cloning and Characterization of a Radial Spoke Head Protein of Sea Urchin Sperm Axonemes: Involvement of the Protein in the Regulation of Sperm Motility

Monoclonal antibodies raised against axonemal proteins of sea urchin spermatozoa have been used to study regulatory mechanisms involved in flagellar motility. Here, we report that one of these antibodies, monoclonal antibody D-316, has an unusual perturbating effect on the motility of sea urchin sperm models; it does not affect the beat frequency, the amplitude of beating or the percentage of motile sperm models, but instead promotes a marked transformation of the flagellar beating pattern which changes from a two-dimensional to a three-dimensional type of movement. On immunoblots of axonemal proteins separated by SDS-PAGE, D-316 recognized a single polypeptide of 90 kDa. This protein was purified following its extraction by exposure of axonemes to a brief heat treatment at 40°C. The protein copurified and coimmunoprecipitated with proteins of 43 and 34 kDa, suggesting that it exists as a complex in its native form. Using D-316 as a probe, a full-length cDNA clone encoding the 90-kDa protein was obtained from a sea urchin cDNA library. The sequence predicts a highly acidic (pI = 4.0) protein of 552 amino acids with a mass of 62,720 Da (p63). Comparison with protein sequences in databases indicated that the protein is related to radial spoke proteins 4 and 6 (RSP4 and RSP6) of Chlamydomonas reinhardtii, which share 37% and 25% similarity, respectively, with p63. However, the sea urchin protein possesses structural features distinct from RSP4 and RSP6, such as the presence of three major acidic stretches which contains 25, 17, and 12 aspartate and glutamate residues of 34-, 22-, and 14-amino acid long stretches, respectively, that are predicted to form α-helical coiled-coil secondary structures. These results suggest a major role for p63 in the maintenance of a planar form of sperm flagellar beating and provide new tools to study the function of radial spoke heads in more evolved species.

Mode matches and their locations in the hydrophobic free energy sequences of peptide ligands and their receptor eigenfunctions

Patterns in sequences of amino acid hydrophobic free energies predict secondary structures in proteins. In protein folding, matches in hydrophobic free energy statistical wavelengths appear to contribute to selective aggregation of secondary structures in “hydrophobic zippers.” In a similar setting, the use of Fourier analysis to characterize the dominant statistical wavelengths of peptide ligands’ and receptor proteins’ hydrophobic modes to predict such matches has been limited by the aliasing and end effects of short peptide lengths, as well as the broad-band, mode multiplicity of many of their frequency (power) spectra. In addition, the sequence locations of the matching modes are lost in this transformation. We make new use of three techniques to address these difficulties: (i) eigenfunction construction from the linear decomposition of the lagged covariance matrices of the ligands and receptors as hydrophobic free energy sequences; (ii) maximum entropy, complex poles power spectra, which select the dominant modes of the hydrophobic free energy sequences or their eigenfunctions; and (iii) discrete, best bases, trigonometric wavelet transformations, which confirm the dominant spectral frequencies of the eigenfunctions and locate them as (absolute valued) moduli in the peptide or receptor sequence. The leading eigenfunction of the covariance matrix of a transmembrane receptor sequence locates the same transmembrane segments seen in n-block-averaged hydropathy plots while leaving the remaining hydrophobic modes unsmoothed and available for further analyses as secondary eigenfunctions. In these receptor eigenfunctions, we find a set of statistical wavelength matches between peptide ligands and their G-protein and tyrosine kinase coupled receptors, ranging across examples from 13.10 amino acids in acid fibroblast growth factor to 2.18 residues in corticotropin releasing factor. We find that the wavelet-located receptor modes in the extracellular loops are compatible with studies of receptor chimeric exchanges and point mutations. A nonbinding corticotropin-releasing factor receptor mutant is shown to have lost the signatory mode common to the normal receptor and its ligand. Hydrophobic free energy eigenfunctions and their transformations offer new quantitative physical homologies in database searches for peptide-receptor matches.

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.

Trans-spliced leader addition to mRNAs in a cnidarian

A search of databases with the sequence from the 5′ untranslated region of a Hydra cDNA clone encoding a receptor protein-tyrosine kinase revealed that a number of Hydra cDNAs contain one of two different sequences at their 5′ ends. This finding suggested the possibility that mRNAs in Hydra receive leader sequences by trans-splicing. This hypothesis was confirmed by the finding that the leader sequences are transcribed as parts of small RNAs encoded by genes located in the 5S rRNA clusters of Hydra. The two spliced leader (SL) RNAs (SL-A and -B) contain splice donor dinucleotides at the predicted positions, and genes that receive SLs contain splice acceptor dinucleotides at the predicted positions. Both of the SL RNAs are bound by antibody against trimethylguanosine, suggesting that they contain a trimethylguanosine cap. The predicted secondary structures of the Hydra SL RNAs show significant differences from the structures predicted for the SLs of other organisms. Messenger RNAs have been identified that can receive either SL-A or -B, although the impact of the two different SLs on the function of the mRNA is unknown. The presence and features of SL addition in the phylum Cnidaria raise interesting questions regarding the evolution of this process.

Rescue of arrested replication forks by homologous recombination

DNA synthesis is an accurate and very processive phenomenon; nevertheless, replication fork progression on chromosomes can be impeded by DNA lesions, DNA secondary structures, or DNA-bound proteins. Elements interfering with the progression of replication forks have been reported to induce rearrangements and/or render homologous recombination essential for viability, in all organisms from bacteria to human. Arrested replication forks may be the target of nucleases, thereby providing a substrate for double-strand break repair enzyme. For example in bacteria, direct fork breakage was proposed to occur at replication forks blocked by a bona fide replication terminator sequence, a specific site that arrests bacterial chromosome replication. Alternatively, an arrested replication fork may be transformed into a recombination substrate by reversal of the forked structures. In reversed forks, the last duplicated portions of the template strands reanneal, allowing the newly synthesized strands to pair. In bacteria, this reaction was proposed to occur in replication mutants, in which fork arrest is caused by a defect in a replication protein, and in UV irradiated cells. Recent studies suggest that it may also occur in eukaryote organisms. We will review here observations that link replication hindrance with DNA rearrangements and the possible underlying molecular processes.

Intrinsic compressibility and volume compression in solvated proteins by molecular dynamics simulation at high pressure.

Constant pressure and temperature molecular dynamics techniques have been employed to investigate the changes in structure and volumes of two globular proteins, superoxide dismutase and lysozyme, under pressure. Compression (the relative changes in the proteins' volumes), computed with the Voronoi technique, is closely related with the so-called protein intrinsic compressibility, estimated by sound velocity measurements. In particular, compression computed with Voronoi volumes predicts, in agreement with experimental estimates, a negative bound water contribution to the apparent protein compression. While the use of van der Waals and molecular volumes underestimates the intrinsic compressibilities of proteins, Voronoi volumes produce results closer to experimental estimates. Remarkably, for two globular proteins of very different secondary structures, we compute identical (within statistical error) protein intrinsic compressions, as predicted by recent experimental studies. Changes in the protein interatomic distances under compression are also investigated. It is found that, on average, short distances compress less than longer ones. This nonuniform contraction underlines the peculiar nature of the structural changes due to pressure in contrast with temperature effects, which instead produce spatially uniform changes in proteins. The structural effects observed in the simulations at high pressure can explain protein compressibility measurements carried out by fluorimetric and hole burning techniques. Finally, the calculation of the proteins static structure factor shows significant shifts in the peaks at short wavenumber as pressure changes. These effects might provide an alternative way to obtain information concerning compressibilities of selected protein regions.

Smoothness within ruggedness: the role of neutrality in adaptation.

RNA secondary structure folding algorithms predict the existence of connected networks of RNA sequences with identical structure. On such networks, evolving populations split into subpopulations, which diffuse independently in sequence space. This demands a distinction between two mutation thresholds: one at which genotypic information is lost and one at which phenotypic information is lost. In between, diffusion enables the search of vast areas in genotype space while still preserving the dominant phenotype. By this dynamic the success of phenotypic adaptation becomes much less sensitive to the initial conditions in genotype space.

A nuclear gene of eubacterial origin in Euglena gracilis reflects cryptic endosymbioses during protist evolution.

Genes for glycolytic and Calvin-cycle glyceraldehyde-3-phosphate dehydrogenase (GAPDH) of higher eukaryotes derive from ancient gene duplications which occurred in eubacterial genomes; both were transferred to the nucleus during the course of endosymbiosis. We have cloned cDNAs encoding chloroplast and cytosolic GAPDH from the early-branching photosynthetic protist Euglena gracilis and have determined the structure of its nuclear gene for cytosolic GAPDH. The gene contains four introns which possess unusual secondary structures, do not obey the GT-AG rule, and are flanked by 2- to 3-bp direct repeats. A gene phylogeny for these sequences in the context of eubacterial homologues indicates that euglenozoa, like higher eukaryotes, have obtained their GAPDH genes from eubacteria via endosymbiotic (organelle-to-nucleus) gene transfer. The data further suggest that the early-branching protists Giardia lamblia and Entamoeba histolytica--which lack mitochondria--and portions of the trypanosome lineage have acquired GAPDH genes from eubacterial donors which did not ultimately give rise to contemporary membrane-bound organelles. Evidence that "cryptic" (possibly ephemeral) endosymbioses during evolution may have entailed successful gene transfer is preserved in protist nuclear gene sequences.