On the structural repertoire of pools of short, random RNA sequences

A detailed knowledge of the mapping between sequence and structure spaces in populations of RNA molecules is essential to better understand their present-day functional properties, to envisage a plausible early evolution of RNA in a prebiotic chemical environment and to improve the design of in vitro evolution experiments, among others. Analysis of natural RNAs, as well as in vitro and computational studies, show that certain RNA structural motifs are much more abundant than others, pointing out a complex relation between sequence and structure. Within this framework, we have investigated computationally the structural properties of a large pool (10 molecules) of single-stranded, 35 nt-long, random RNA sequences. The secondary structures obtained are ranked and classified into structure families. The number of structures in main families is analytically calculated and compared with the numerical results. This permits a quantification of the fraction of structure space covered by a large pool of sequences. We further show that the number of structural motifs and their frequency is highly unbalanced with respect to the nucleotide composition: simple structures such as stem-loops and hairpins arise from sequences depleted in G, while more complex structures require an enrichment of G. In general, we observe a strong correlation between subfamilies-characterized by a fixed number of paired nucleotides-and nucleotide composition. Our results are compared to the structural repertoire obtained in a second pool where isolated base pairs are prohibited. © 2008 Elsevier Ltd. All rights reserved.

Phenotypic effect of mutations in evolving populations of RNA molecules

Background. The secondary structure of folded RNA sequences is a good model to map phenotype onto genotype, as represented by the RNA sequence. Computational studies of the evolution of ensembles of RNA molecules towards target secondary structures yield valuable clues to the mechanisms behind adaptation of complex populations. The relationship between the space of sequences and structures, the organization of RNA ensembles at mutation-selection equilibrium, the time of adaptation as a function of the population parameters, the presence of collective effects in quasispecies, or the optimal mutation rates to promote adaptation all are issues that can be explored within this framework. Results. We investigate the effect of microscopic mutations on the phenotype of RNA molecules during their in silico evolution and adaptation. We calculate the distribution of the effects of mutations on fitness, the relative fractions of beneficial and deleterious mutations and the corresponding selection coefficients for populations evolving under different mutation rates. Three different situations are explored: the mutation-selection equilibrium (optimized population) in three different fitness landscapes, the dynamics during adaptation towards a goal structure (adapting population), and the behavior under periodic population bottlenecks (perturbed population). Conclusions. The ratio between the number of beneficial and deleterious mutations experienced by a population of RNA sequences increases with the value of the mutation rate µ at which evolution proceeds. In contrast, the selective value of mutations remains almost constant, independent of µ, indicating that adaptation occurs through an increase in the amount of beneficial mutations, with little variations in the average effect they have on fitness. Statistical analyses of the distribution of fitness effects reveal that small effects, either beneficial or deleterious, are well described by a Pareto distribution. These results are robust under changes in the fitness landscape, remarkably when, in addition to selecting a target secondary structure, specific subsequences or low-energy folds are required. A population perturbed by bottlenecks behaves similarly to an adapting population, struggling to return to the optimized state. Whether it can survive in the long run or whether it goes extinct depends critically on the length of the time interval between bottlenecks. © 2010 Stich et al; licensee BioMed Central Ltd.

Characterization of the structure of RAMP1 by mutagenesis and molecular modeling

Receptor activity modifying proteins (RAMPs) are a family of single-pass transmembrane proteins that dimerize with G-protein-coupled receptors. They may alter the ligand recognition properties of the receptors (particularly for the calcitonin receptor-like receptor, CLR). Very little structural information is available about RAMPs. Here, an ab initio model has been generated for the extracellular domain of RAMP1. The disulfide bond arrangement (Cys 27-Cys82, Cys40-Cys72, and Cys 57-Cys104) was determined by site-directed mutagenesis. The secondary structure (a-helices from residues 29-51, 60-80, and 87-100) was established from a consensus of predictive routines. Using these constraints, an assemblage of 25,000 structures was constructed and these were ranked using an all-atom statistical potential. The best 1000 conformations were energy minimized. The lowest scoring model was refined by molecular dynamics simulation. To validate our strategy, the same methods were applied to three proteins of known structure; PDB:1HP8, PDB:1V54 chain H (residues 21-85), and PDB:1T0P. When compared to the crystal structures, the models had root mean-square deviations of 3.8 Å, 4.1 Å, and 4.0 Å, respectively. The model of RAMP1 suggested that Phe93, Tyr 100, and Phe101 form a binding interface for CLR, whereas Trp74 and Phe92 may interact with ligands that bind to the CLR/RAMP1 heterodimer. © 2006 by the Biophysical Society.

Forced retroevolution of an RNA bacteriophage

The operator hairpin ahead of the replicase gene in RNA bacteriophage MS2 contains overlapping signals for binding the coat protein and ribosomes. Coat protein binding inhibits further translation of the gene and forms the first step in capsid formation. The hairpin sequence was partially randomized to assess the importance of this structure element for the bacteriophage and to monitor alternative solutions that would evolve on the passaging of mutant phages. The evolutionary reconstruction of the operator failed in the majority of mutants. Instead, a poor imitation developed containing only some of the recognition signals for the coat protein. Three mutants were of particular interest in that they contained double nonsense codons in the lysis reading frame that runs through the operator hairpin. The simultaneous reversion of two stop codons into sense codons has a very low probability of occurring. Therefore the phage solved the problem by deleting the nonsense signals and, in fact, the complete operator, except for the initiation codon of the replicase gene. Several revertants were isolated with activities ranging from 1% to 20% of wild type. The operator, long thought to be a critical regulator, now appears to be a dispensable element. In addition, the results indicate how RNA viruses can be forced to step back to an attenuated form.

Long-range translational coupling in single-stranded RNA bacteriophages: an evolutionary analysis

In coliphage MS2 RNA a long-distance interaction (LDI) between an internal segment of the upstream coat gene and the start region of the replicase gene prevents initiation of replicase synthesis in the absence of coat gene translation. Elongating ribosomes break up the repressor LDI and thus activate the hidden initiation site. Expression studies on partial MS2 cDNA clones identified base pairing between 1427-1433 and 1738-1744, the so-called Min Jou (MJ) interaction, as the molecular basis for the long-range coupling mechanism. Here, we examine the biological significance of this interaction for the control of replicase gene translation. The LDI was disrupted by mutations in the 3'-side and the evolutionary adaptation was monitored upon phage passaging. Two categories of pseudorevertants emerged. The first type had restored the MJ interaction but not necessarily the native sequence. The pseudorevertants of the second type acquired a compensatory substitution some 80 nt downstream of the MJ interaction that stabilizes an adjacent LDI. In one examined case we confirmed that the second site mutations had restored coat-replicase translational coupling. Our results show the importance of translational control for fitness of the phage. They also reveal that the structure that buries the replicase start extends to structure elements bordering the MJ interaction.

Structure, dynamics, and energetics of siRNA-cationic vector complexation:a molecular dynamics study

The design and synthesis of safe and efficient nonviral vectors for gene delivery has attracted significant attention in recent years. Previous experiments have revealed that the charge density of a polycation (the carrier) plays a crucial role in complexation and the release of the gene from the complex in the cytosol. In this work, we adopt an atomistic molecular dynamics simulation approach to study the complexation of short strand duplex RNA with six cationic carrier systems of varying charge and surface topology. The simulations reveal detailed molecular-level pictures of the structures and dynamics of the RNA-polycation complexes. Estimates for the binding free energy indicate that electrostatic contributions are dominant followed by van der Waals interactions. The binding free energy between the 8(+)polymers and the RNA is found to be larger than that of the 4(+)polymers, in general agreement with previously published data. Because reliable binding free energies provide an effective index of the ability of the polycationic carrier to bind the nucleic acid and also carry implications for the process of gene release within the cytosol, these novel simulations have the potential to provide us with a much better understanding of key mechanistic aspects of gene-polycation complexation and thereby advance the rational design of nonviral gene delivery systems.

Flexibility of short-strand RNA in aqueous solution as revealed by molecular dynamics simulation:are A′-RNA and A-RNA distinct conformational structures?

We use molecular dynamics simulations to compare the conformational structure and dynamics of a 21-base pair RNA sequence initially constructed according to the canonical A-RNA and A'-RNA forms in the presence of counterions and explicit water. Our study aims to add a dynamical perspective to the solid-state structural information that has been derived from X-ray data for these two characteristic forms of RNA. Analysis of the three main structural descriptors commonly used to differentiate between the two forms of RNA namely major groove width, inclination and the number of base pairs in a helical twist over a 30 ns simulation period reveals a flexible structure in aqueous solution with fluctuations in the values of these structural parameters encompassing the range between the two crystal forms and more. This provides evidence to suggest that the identification of distinct A-RNA and A'-RNA structures, while relevant in the crystalline form, may not be generally relevant in the context of RNA in the aqueous phase. The apparent structural flexibility observed in our simulations is likely to bear ramifications for the interactions of RNA with biological molecules (e.g. proteins) and non-biological molecules (e.g. non-viral gene delivery vectors). © CSIRO 2009.

Structure and dynamics of multiple cationic vectors-siRNA complexation by all-atomic molecular dynamics simulations

Understanding the molecular mechanism of gene condensation is a key component to rationalizing gene delivery phenomena, including functional properties such as the stability of the gene-vector complex and the intracellular release of the gene. In this work, we adopt an atomistic molecular dynamics simulation approach to study the complexation of short strand duplex RNA with four cationic carrier systems of varying charge and surface topology at different charge ratios. At lower charge ratios, polymers bind quite effectively to siRNA, while at high charge ratios, the complexes are saturated and there are free polymers that are unable to associate with RNA. We also observed reduced fluctuations in RNA structures when complexed with multiple polymers in solution as compared to both free siRNA in water and the single polymer complexes. These novel simulations provide a much better understanding of key mechanistic aspects of gene-polycation complexation and thereby advance progress toward rational design of nonviral gene delivery systems.

Coalescence of secondary dispersions

A study has been made of the coalescence of secondary dispersions in beds of woven meshes. The variables investigated were superficial velocity, bed depth, mesh geometry and fibre material; the effects of presoaking the bed in the dispersed phase before operation were also considered. Equipment was design~d to generate a 0.1% phase ratio toluene in water dispersion whose mean drop size was determined using a Coulter Counter. The coalesced drops were sized by photography and a novel holographic technique was developed to evaluate the mean diameter of the effluent secondary drops. Previous models describing single phase flow in porous media are reviewed and it was found that the experimental data obtained in this study is best represented by Keller's equation which is based on a physical model similar to the internal structure of the meshes. Statistical analysis of two phase data produced a correlation, for each mesh tested, relating the pressure drop to superficial velocity and bed depth. The flow parameter evaluated from the single phase model is incorporated into a theoretical comparison of drop capture mechanisms which indicated that direct and indirect interception are predominant. The resulting equation for drop capture efficiericy is used to predict the initial, local drop capture rate in a coalescer. A mathematical description of the saturation profiles was formulated and verified by average saturation data. Based 6n the Blake-Kozeny equation, an expression is derived analytically to predict the two phase pressure drop using the parameters which characterise the saturation profiles. By specifying the local saturation at the inlet face for a given velocity, good agreement between experimental pressure drop data and the model predictions was obtained.

AIDS: structure, replication and spread of the AIDS virus

AIDS (Acquired Immune Deficiency Syndrome)was first, described as a new disease of humans in 1981. The origins of the disease are controversial. AIDS is caused by a retrovirus, a type of virus which rarely attacks human cells. The first virus of this type recorded in humans is reponsible for a type of leukaemia and was identified in 1978. AIDS is thus the third type of human retrovirus to be discovered and hence, is referred to as T-lymphotrophic virus III (HTLV-III). For viruses to replicate, they have to invade a host cell which in this case is a T4-lymphocyte, a type of white blood cell that regulates the immune system. The problems of the disease result directly from the death of these cells. As a consequence, the immune system is compromised leading to a number of opportunistic secondary infections and other disorders.

The dawn of the RNA World:toward functional complexity through ligation of random RNA oligomers

A main unsolved problem in the RNA World scenario for the origin of life is how a template-dependent RNA polymerase ribozyme emerged from short RNA oligomers obtained by random polymerization on mineral surfaces. A number of computational studies have shown that the structural repertoire yielded by that process is dominated by topologically simple structures, notably hairpin-like ones. A fraction of these could display RNA ligase activity and catalyze the assembly of larger, eventually functional RNA molecules retaining their previous modular structure: molecular complexity increases but template replication is absent. This allows us to build up a stepwise model of ligation- based, modular evolution that could pave the way to the emergence of a ribozyme with RNA replicase activity, step at which information-driven Darwinian evolution would be triggered. Copyright © 2009 RNA Society.

Responsive hydrophobically associating polymers:a review of structure and properties

Responsive hydrophobically associating polymers can in many ways be considered to be analogous to proteins in their ability to form compact molecules with a defined secondary structure, and hence, functionality. These molecules are characterized by the presence of alternating charged and hydrophobic groups. The balance between charge repulsion and hydrophobic interactions is sensitive to environmental pH and therefore changes in pH produce controllable conformational changes. The change from a charged extended chain to a collapsed uncharged coil structure is sometimes referred to as hypercoiling behaviour and enables the polymer to act as a simple switch between an 'on' and 'off' state. The purpose of this review is to illustrate the structure and behaviour of polymers that exhibit hypercoiling behaviour and to highlight their potential pharmaceutical applications, which in terms of drug delivery is likely to be related to their surface behaviour and solubilizing activity. © 2001 Elsevier Science B.V. All rights reserved.

Modifying the hierarchical porosity of SBA-15 via mild-detemplation followed by secondary treatments

Fenton-chemistry-based detemplation combined with secondary treatments offers options to tune the hierarchical porosity of SBA-15. This approach has been studied on a series of SBA-15 mesophases and has been compared to the conventional calcination. The as-synthesized and detemplated materials were studied with regard to their template content (TGA, CHN), structure (SAXS, TEM), surface hydroxylation (Blin-Carterets approach), and texture (high-resolution argon physisorption). Fenton detemplation achieves 99% of template removal, leading to highly hydroxylated materials. The structure is better preserved when a secondary treatment is applied after the Fenton oxidation, due to the intense capillary forces during drying in water. Two successful approaches are presented: drying in a low-surface-tension solvent (such as n-BuOH) and a hydrothermal stabilization to further condense the structure and make it structurally more robust. Both approaches give rise to remarkably low structural shrinkage, lower than calcination and the direct water-dried Fenton. Interestingly, the derived textural features are remarkably different. The n-BuOH exchange route gives rise to highly hierarchical structures with enhanced interconnecting pores and the highest surface areas. The hydrothermal stabilization produces large-pore SBA-15 structures with high pore volume, intermediate interconnectivity, and minimal micropores. Therefore, the hierarchical texture can be fine-tuned in these two fashions while the template is removed under mild conditions.