Dataset Information

Loop-loop interaction in an adenine-sensing riboswitch: a molecular dynamics study.

ABSTRACT: Riboswitches are mRNA-based molecules capable of controlling the expression of genes. They undergo conformational changes upon ligand binding, and as a result, they inhibit or promote the expression of the associated gene. The close connection between structural rearrangement and function makes a detailed knowledge of the molecular interactions an important step to understand the riboswitch mechanism and efficiency. We have performed all-atom molecular dynamics simulations of the adenine-sensing add A-riboswitch to study the breaking of the kissing loop, one key tertiary element in the aptamer structure. We investigated the aptamer domain of the add A-riboswitch in complex with its cognate ligand and in the absence of the ligand. The opening of the hairpins was simulated using umbrella sampling using the distance between two loops as the reaction coordinate. A two-step process was observed in all the simulated systems. First, a general loss of stacking and hydrogen bond interactions is seen. The last interactions that break are the two base pairs G37-C61 and G38-C60, but the break does not affect the energy profile, indicating their pivotal role in the tertiary structure formation but not in the structure stabilization. The junction area is partially organized before the kissing loop formation and residue A24 anchors together the loop helices. Moreover, when the distance between the loops is increased, one of the hairpins showed more flexibility by changing its orientation in the structure, while the other conserved its coaxial arrangement with the rest of the structure.

SUBMITTER: Allner O

PROVIDER: S-EPMC3683926 | biostudies-literature | 2013 Jul

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Loop-loop interaction in an adenine-sensing riboswitch: a molecular dynamics study.

Allnér Olof O Nilsson Lennart L Villa Alessandra A

RNA (New York, N.Y.) 20130528 7

Riboswitches are mRNA-based molecules capable of controlling the expression of genes. They undergo conformational changes upon ligand binding, and as a result, they inhibit or promote the expression of the associated gene. The close connection between structural rearrangement and function makes a detailed knowledge of the molecular interactions an important step to understand the riboswitch mechanism and efficiency. We have performed all-atom molecular dynamics simulations of the adenine-sensing ...[more]

PMID: 23716711

Similar Datasets

Project description:IntroductionRiboswitches are cis-acting regulatory RNA elements prevalently located in the leader sequences of bacterial mRNA. An adenine sensing riboswitch cis-regulates adeninosine deaminase gene (add) in Vibrio vulnificus. The structural mechanism regulating its conformational changes upon ligand binding mostly remains to be elucidated. In this open framework it has been suggested that the ligand stabilizes the interaction of the distal "kissing loop" complex. Using accurate full-atom molecular dynamics with explicit solvent in combination with enhanced sampling techniques and advanced analysis methods it could be possible to provide a more detailed perspective on the formation of these tertiary contacts.MethodsIn this work, we used umbrella sampling simulations to study the thermodynamics of the kissing loop complex in the presence and in the absence of the cognate ligand. We enforced the breaking/formation of the loop-loop interaction restraining the distance between the two loops. We also assessed the convergence of the results by using two alternative initialization protocols. A structural analysis was performed using a novel approach to analyze base contacts.ResultsContacts between the two loops were progressively lost when larger inter-loop distances were enforced. Inter-loop Watson-Crick contacts survived at larger separation when compared with non-canonical pairing and stacking interactions. Intra-loop stacking contacts remained formed upon loop undocking. Our simulations qualitatively indicated that the ligand could stabilize the kissing loop complex. We also compared with previously published simulation studies.Discussion and conclusionsKissing complex stabilization given by the ligand was compatible with available experimental data. However, the dependence of its value on the initialization protocol of the umbrella sampling simulations posed some questions on the quantitative interpretation of the results and called for better converged enhanced sampling simulations.

Project description:Expression of a large fraction of genes in bacteria is controlled by riboswitches, which are found in the untranslated region of mRNA. Structurally riboswitches have a conserved aptamer domain to which a metabolite binds, resulting in a conformational change in the downstream expression platform. Prediction of the functions of riboswitches requires a quantitative description of the folding landscape so that the barriers and time scales for the conformational change in the switching region in the aptamer can be estimated. Using a combination of all atom molecular dynamics (MD) and coarse-grained model simulations we studied the response of adenine (A) binding add and pbuE A-riboswitches to mechanical force. The two riboswitches contain a structurally similar three-way junction formed by three paired helices, P1, P2, and P3, but carry out different functions. Using pulling simulations, with structures generated in MD simulations, we show that after P1 rips the dominant unfolding pathway in the add A-riboswitch is the rupture of P2 followed by unraveling of P3. In the pbuE A-riboswitch, after P1 unfolds P3 ruptures ahead of P2. The order of unfolding of the helices, which is in accord with single molecule pulling experiments, is determined by the relative stabilities of the individual helices. Our results show that the stability of isolated helices determines the order of assembly and response to force in these non-coding regions. We use the simulated free energy profile for the pbuE A-riboswitch to estimate the time scale for allosteric switching, which shows that this riboswitch is under kinetic control lending additional support to the conclusion based on single molecule pulling experiments. A consequence of the stability hypothesis is that a single point mutation (U28C) in the P2 helix of the add A-riboswitch, which increases the stability of P2, would make the folding landscapes of the two riboswitches similar. This prediction can be tested in single molecule pulling experiments.

Dataset Information

Loop-loop interaction in an adenine-sensing riboswitch: a molecular dynamics study.

Publications

Loop-loop interaction in an adenine-sensing riboswitch: a molecular dynamics study.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets