Conservation of mRNA secondary structures may filter out mutations in Escherichia coli evolution

Recent reports indicate that mutations in viral genomes tend to preserve RNA secondary structure, and those mutations that disrupt secondary structural elements may reduce gene expression levels, thereby serving as a functional knockout. In this article, we explore the conservation of secondary structures of mRNA coding regions, a previously unknown factor in bacterial evolution, by comparing the structural consequences of mutations in essential and nonessential Escherichia coli genes accumulated over 40 000 generations in the course of the ‘long-term evolution experiment’. We monitored the extent to which mutations influence minimum free energy (MFE) values, assuming that a substantial change in MFE is indicative of structural perturbation. Our principal finding is that purifying selection tends to eliminate those mutations in essential genes that lead to greater changes of MFE values and, therefore, may be more disruptive for the corresponding mRNA secondary structures. This effect implies that synonymous mutations disrupting mRNA secondary structures may directly affect the fitness of the organism. These results demonstrate that the need to maintain intact mRNA structures imposes additional evolutionary constraints on bacterial genomes, which go beyond preservation of structure and function of the encoded proteins.     

Context-Dependent Remodeling of Structure in Two Large Protein Fragments

Protein folding involves the formation of secondary structural elements from the primary sequence and their association with tertiary assemblies. The relation of this primary sequence to a specific folded protein structure remains a central question in structural biology. An increasing body of evidence suggests that variations in homologous sequence ranging from point mutations to substantial insertions or deletions can yield stable proteins with markedly different folds. Here we report the structural characterization of domain IV (D4) and ΔD4 (polypeptides with 222 and 160 amino acids, respectively) that differ by virtue of an N-terminal deletion of 62 amino acids (28% of the overall D4 sequence). The high-resolution crystal structures of the monomeric D4 and the dimeric ΔD4 reveal substantially different folds despite an overall conservation of secondary structure. These structures show that the formation of tertiary structures, even in extended polypeptide sequences, can be highly context dependent, and they serve as a model for structural plasticity in protein isoforms.     

Relation between sequence and structure of HIV-1 protease inhibitor complexes: a model system for the analysis of protein flexibility.

The flexibility of different regions of HIV-1 protease was examined by using a database consisting of 73 X-ray structures that differ in terms of sequence, ligands or both. The root-mean-square differences of the backbone for the set of structures were shown to have the same variation with residue number as those obtained from molecular dynamics simulations, normal mode analyses and X-ray B-factors. This supports the idea that observed structural changes provide a measure of the inherent flexibility of the protein, although specific interactions between the protease and the ligand play a secondary role. The results suggest that the potential energy surface of the HIV-1 protease is characterized by many local minima with small energetic differences, some of which are sampled by the different X-ray structures of the HIV-1 protease complexes. Interdomain correlated motions were calculated from the structural fluctuations and the results were also in agreement with molecular dynamics simulations and normal mode analyses. Implications of the results for the drug-resistance engendered by mutations are discussed briefly.     

The human immunodeficiency virus type 1 encapsidation site is a multipartite RNA element composed of functional hairpin structures.

We analyzed the leader region of human immunodeficiency virus type 1 (HIV-1) RNA to decipher the nature of the cis-acting E/psi element required for encapsidation of viral RNA into virus particles. Our data indicate that, for RNA encapsidation, there are at least two functional subregions in the leader region. One subregion is located at a position immediately proximal to the major splice donor, and the second is located between the splice donor and the beginning of the gag gene. This suggests that at least two discrete cis-acting elements are recognition signals for encapsidation. To determine whether specific putative RNA secondary structures serve as the signal(s) for encapsidation, we constructed primary base substitution mutations that would be expected to destabilize these potential structures and second-site compensatory mutations that would restore secondary structure. Analysis of these mutants allowed the identification of two discrete hairpins that facilitate RNA encapsidation in vivo. Thus, the HIV-1 E/psi region is a multipartite element composed of specific and functional RNA secondary structures. Compensation of the primary mutations by the second-site mutations could not be attained in trans. This indicates that interstrand base pairing between these two stem regions within the hairpins does not appear to be the basis for HIV-1 RNA dimer formation. Comparison of the hypothetical RNA secondary structures from 10 replication-competent HIV-1 strains suggests that a subset of the hydrogen-bonded base pairs within the stems of the hairpins is likely to be required for function in cis.     

Sequence swapping does not result in conformation swapping for the beta4/beta5 and beta8/beta9 beta-hairpin turns in human acidic fibroblast growth factor

The beta-turn is the most common type of nonrepetitive structure in globular proteins, comprising ~25% of all residues; however, a detailed understanding of effects of specific residues upon beta-turn stability and conformation is lacking. Human acidic fibroblast growth factor (FGF-1) is a member of the beta-trefoil superfold and contains a total of five beta-hairpin structures (antiparallel beta-sheets connected by a reverse turn). beta-Turns related by the characteristic threefold structural symmetry of this superfold exhibit different primary structures, and in some cases, different secondary structures. As such, they represent a useful system with which to study the role that turn sequences play in determining structure, stability, and folding of the protein. Two turns related by the threefold structural symmetry, the beta4/beta5 and beta8/beta9 turns, were subjected to both sequence-swapping and poly-glycine substitution mutations, and the effects upon stability, folding, and structure were investigated. In the wild-type protein these turns are of identical length, but exhibit different conformations. These conformations were observed to be retained during sequence-swapping and glycine substitution mutagenesis. The results indicate that the beta-turn structure at these positions is not determined by the turn sequence. Structural analysis suggests that residues flanking the turn are a primary structural determinant of the conformation within the turn.     

Binding of Amphipathic Cell Penetrating Peptide p28 to Wild Type and Mutated p53 as studied by Raman,Atomic Force and Surface Plasmon Resonance spectroscopies

BackgroundMutations within the DNA binding domain (DBD) of the tumor suppressor p53 are found in > 50% of human cancers and may significantly modify p53 secondary structure impairing its function. p28, an amphipathic cell-penetrating peptide, binds to the DBD through hydrophobic interaction and induces a posttranslational increase in wildtype and mutant p53 restoring functionality. We use mutation analyses to explore which elements of secondary structure may be critical to p28 binding.MethodsMolecular modeling, Raman spectroscopy, Atomic Force Spectroscopy (AFS) and Surface Plasmon Resonance (SPR) were used to identify which secondary structure of site-directed and naturally occurring mutant DBDs are potentially altered by discrete changes in hydrophobicity and the molecular interaction with p28.ResultsWe show that specific point mutations that alter hydrophobicity within non-mutable and mutable regions of the p53 DBD alter specific secondary structures. The affinity of p28 was positively correlated with the β-sheet content of a mutant DBD, and reduced by an increase in unstructured or random coil that resulted from a loss in hydrophobicity and redistribution of surface charge.ConclusionsThese results help refine our knowledge of how mutations within p53-DBD alter secondary structure and provide insight on how potential structural alterations in p28 or similar molecules improve their ability to restore p53 function.General significanceRaman spectroscopy, AFS, SPR and computational modeling are useful approaches to characterize how mutations within the p53DBD potentially affect secondary structure and identify those structural elements prone to influence the binding affinity of agents designed to increase the functionality of p53.     

Secondary structure of the 3'-noncoding region of flavivirus genomes: comparative analysis of base pairing probabilities.

The prediction of the complete matrix of base pairing probabilities was applied to the 3'' noncoding region (NCR) of flavivirus genomes. This approach identifies not only well-defined secondary structure elements, but also regions of high structural flexibility. Flaviviruses, many of which are important human pathogens, have a common genomic organization, but exhibit a significant degree of RNA sequence diversity in the functionally important 3''-NCR. We demonstrate the presence of secondary structures shared by all flaviviruses, as well as structural features that are characteristic for groups of viruses within the genus reflecting the established classification scheme. The significance of most of the predicted structures is corroborated by compensatory mutations. The availability of infectious clones for several flaviviruses will allow the assessment of these structural elements in processes of the viral life cycle, such as replication and assembly.     

Functional roles of loops 3 and 4 in the cyclic nucleotide binding domain of cyclic AMP receptor protein from Escherichia coli

Cyclic AMP is a ubiquitous secondary message that regulates a large variety of functions. The protein structural motif that binds cAMP is highly conserved with the exception of loops 3 and 4, whose structure and length are variable. The cAMP receptor protein of Escherichia coli, CRP, was employed as a model system to elucidate the functional roles of these loops. Based on the sequence differences between CRP and cyclic nucleotide gated channel, three mutants of CRP were constructed: deletion (residues 54-56 in loop 3 were deleted), insertion (loop 4 was lengthened by 5 residues between Glu-78 and Gly-79) and double mutants. The effects of these mutations on the structure and function of CRP were monitored. Results show that the deletion and insertion mutations do not significantly change the secondary structure of CRP, although the tertiary and quaternary structures are perturbed. The functional data indicate that loop 3 modulates the binding affinities of cAMP and DNA. Although the lengthened loop 4 may have some fine-tuning functions, the specific function of the original loop 4 of CRP remains uncertain. The function consequences of mutation in loop 3 of CRP are similar to that of site A and site B in the regulatory subunits of cyclic AMP-dependent protein kinases. Thus, the roles played by loop 3 in CRP may represent a more common mechanism employed by cyclic nucleotide binding domain in modulating ligand binding affinity and intramolecular communication.     

Predicted changes in pre-mRNA secondary structure vary in their association with exon skipping for mutations in exons 2, 4, and 8 of the Hprt gene and exon 51 of the fibrillin gene

Exon skipping that accompanies exonic mutation might be caused by an effect of the mutation on pre-mRNA secondary structure. Previous attempts to associate predicted secondary structure of pre-mRNA with exon skipping have been hindered by either a small number of available mutations, sub-optimal structures, or weak effects on exon skipping. This report identifies more extensive sets of mutations from the human and hamster Hprt gene whose association with exon skipping is clear. Optimal secondary structures of the wild-type and mutant pre-mRNA surrounding each exon were predicted by energy minimization and were compared by energy dot plots. A significant association was found between the occurrence of exon skipping and the disruption of a stem containing the acceptor site consensus sequences of exon 8 of the human Hprt gene. However, no change in secondary structure was associated with skipping of exon 4 of the hamster Hprt gene. Using updated energy parameters we found a different structure than that previously reported for exon 2 of the hamster Hprt gene. In contrast to the previously reported structure, no significant association was found between predicted structural changes and skipping of exon 2. For all three Hprt exons studied, there was a significantly greater number of deoxythymidine substitutions among mutations accompanied by exon skipping than among mutations without exon skipping. For exon 8, deoxythymidine substitution was also associated with structural changes in the stem containing the acceptor site consensus sequences. For exon 51 of the human fibrillin gene, structural differences from wild type were predicted for all four mutations accompanied by exon skipping that were not were predicted for a single mutation without exon skipping. Our results suggest that both primary and secondary pre-mRNA structure contribute to definition of Hprt exons, which may involve exonic splicing enhancers.     

An RNA tertiary structure in the 3' untranslated region of enteroviruses is necessary for efficient replication.

RNA tertiary structures, such as pseudoknots, are known to be biologically significant in a number of virus systems. The 3' untranslated regions of the RNA genomes of all members of the Enterovirus genus of Picornaviridae exhibit a potential, pseudoknot-like, tertiary structure interaction of an unusual type. This is formed by base pairing between loop regions of two secondary structure domains. It is distinct from a potential, conventional pseudoknot, studied previously in poliovirus, which is less conserved phylogenetically. We have analyzed the tertiary structure feature in one enterovirus, coxsackievirus A9, using specific mutagenesis. A double mutant in which the potential interaction was destroyed was nonviable, and viability was restored by introducing compensating mutations, predicted to allow the interaction to reform. Phenotypic pseudorevertants of virus mutants, having mutations designed to disrupt the interaction, were all found to have acquired nucleotide changes which restored the potential interaction. Analysis of one mutant containing a single-base mutation indicated a greatly increased temperature sensitivity due to a step early in replication. The results show that, in addition to secondary structures, tertiary RNA structural interactions can play an important role in the biology of picornaviruses.     

Search for characteristic structural features of mammalian mitochondrial tRNAs

A number of mitochondrial (mt) tRNAs have strong structural deviations from the classical tRNA cloverleaf secondary structure and from the conventional L-shaped tertiary structure. As a consequence, there is a general trend to consider all mitochondrial tRNAs as "bizarre" tRNAs. Here, a large sequence comparison of the 22 tRNA genes within 31 fully sequenced mammalian mt genomes has been performed to define the structural characteristics of this specific group of tRNAs. Vertical alignments define the degree of conservation/variability of primary sequences and secondary structures and search for potential tertiary interactions within each of the 22 families. Further horizontal alignments ascertain that, with the exception of serine-specific tRNAs, mammalian mt tRNAs do fold into cloverleaf structures with mostly classical features. However, deviations exist and concern large variations in size of the D- and T-loops. The predominant absence of the conserved nucleotides G18G19 and T54T55C56, respectively in these loops, suggests that classical tertiary interactions between both domains do not take place. Classification of the tRNA sequences according to their genomic origin (G-rich or G-poor DNA strand) highlight specific features such as richness/poorness in mismatches or G-T pairs in stems and extremely low G-content or C-content in the D- and T-loops. The resulting 22 "typical" mammalian mitochondrial sequences built up a phylogenetic basis for experimental structural and functional investigations. Moreover, they are expected to help in the evaluation of the possible impacts of those point mutations detected in human mitochondrial tRNA genes and correlated with pathologies.     

Correlation between the structural stability and aggregation propensity of proteins

Protein aggregation, being an outcome of improper protein folding, is largely dependent on the folding kinetics of a protein. Previous studies have reported a positive correlation between the stability of the secondary structural elements of a protein and their rate of folding/unfolding. In this in silico study, the secondary and tertiary structures of proteins a) that form inclusion bodies on overexpression in Escherichia coli, b) that form amyloid fibrils and c) that are soluble on overexpression in E. coli are analyzed for certain features that are known to be associated with structural stability. The study revealed that the soluble proteins seem to have a higher rate of folding (based on contact order) and a lower percentage of exposed hydrophobic residues as compared to the inclusion body forming or amyloidogenic proteins. The soluble proteins also seem to have a more favored helix and strand composition (based on the known secondary structural propensities of amino acids). The secondary structure analyses also reveal that the evolutionary pressure is directed against protein aggregation. This understanding of the positive correlation between structural stability and solubility, along with the other parameters known to influence aggregation, could be exploited in the design of mutations aimed at reducing the aggregation propensity of the proteins.     

同义密码子携带多少蛋白质二级结构信息

应用信息论方法考察了大肠杆菌人两种生物的同义密码子用语和蛋白质二级结构的关联情况。研究结果表明：大肠杆菌和人的基因组中都存在着一些同义密码子明显携带有蛋白质二级结构信息,尽管这些信息量都很小;同义密码子与蛋白质二级结构的关联是种属特异性。     

An algorithm for searching RNA motifs in genomic sequences

RNA molecules, which are found in all living cells, fold into characteristic structures that account for their diverse functional activities. Many of these RNA structures consist of a collection of fundamental RNA motifs. The various combinations of RNA basic components form different RNA classes and define their unique structural and functional properties. The availability of many genome sequences makes it possible to search computationally for functional RNAs. Biological experiments indicate that functional RNAs have characteristic RNA structural motifs represented by specific combinations of base pairings and conserved nucleotides in the loop regions. The searching for those well-ordered RNA structures and their homologues in genomic sequences is very helpful for the understanding of RNA-based gene regulation. In this paper, we consider the following problem: given an RNA sequence with a known secondary structure, efficiently determine candidate segments in genomic sequences that can potentially form RNA secondary structures similar to the given RNA secondary structure. Our new bottom-up approach searches all potential stem-loops similar to ones of the given RNA secondary structure first, and then based on located stem-loops, detects potential homologous structural RNAs in genomic sequences.     

RNA structures with pseudo-knots: Graph-theoretical, combinatorial, and statistical properties

The secondary structures of nucleic acids form a particularly important class of contact structures. Many important RNA molecules, however, contain pseudo-knots, a structural feature that is excluded explicitly from the conventional definition of secondary structures. We propose here a generalization of secondary structures incorporating ‘non-nested’ pseudo-knots, which we call bi-secondary structures, and discuss measures for the complexity of more general contact structures based on their graph-theoretical properties. Bi-secondary structures are planar trivalent graphs that are characterized by special embedding properties. We derive exact upper bounds on their number (as a function of the chain length n) implying that there are fewer different structures than sequences. Computational results show that the number of bi-secondary structures grows approximately like 2.35ⁿ. Numerical studies based on kinetic folding and a simple extension of the standard energy model show that the global features of the sequence-structure map of RNA do not change when pseudo-knots are introduced into the secondary structure picture. We find a large fraction of neutral mutations and, in particular, networks of sequences that fold into the same shape. These neutral networks percolate through the entire sequence space.     

Temperature-controlled structural alterations of an RNA thermometer

Thermoresponsive structures in the 5'-untranslated region of mRNA are known to control translation of heat shock and virulence genes. Expression of many rhizobial heat shock genes is regulated by a conserved sequence element called ROSE for repression of heat shock gene expression. This cis-acting, untranslated mRNA is thought to prevent ribosome access at low temperature through an extended secondary structure, which partially melts when the temperature rises. We show here by a series of in vivo and in vitro approaches that ROSE is a sensitive thermometer responding in the physiologically relevant temperature range between 30 and 40 degrees C. Point mutations predicted to disrupt base pairing enhanced expression at 30 degrees C. Compensatory mutations restored repression, emphasizing the importance of secondary structures in the sensory RNA. Only moderate inducibility of a 5'-truncated ROSE variant suggests that interactions between individual stem loops coordinate temperature sensing. In the presence of a complementary oligonucleotide, the functionally important stem loop of ROSE was rendered susceptible to RNase H treatment at heat shock temperatures. Since major structural rearrangements were not observed during UV and CD spectroscopy, subtle structural changes involving the Shine-Dalgarno sequence are proposed to mediate translational control. Temperature perception by the sensory RNA is an ordered process that most likely occurs without the aid of accessory factors.