Browse
Submit Data
Databases
API
Help

Dataset Information

25 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Binding free energy decomposition and multiple unbinding paths of buried ligands in a PreQ1 riboswitch.

ABSTRACT: Riboswitches are naturally occurring RNA elements that control bacterial gene expression by binding to specific small molecules. They serve as important models for RNA-small molecule recognition and have also become a novel class of targets for developing antibiotics. Here, we carried out conventional and enhanced-sampling molecular dynamics (MD) simulations, totaling 153.5 μs, to characterize the determinants of binding free energies and unbinding paths for the cognate and synthetic ligands of a PreQ1 riboswitch. Binding free energy analysis showed that two triplets of nucleotides, U6-C15-A29 and G5-G11-C16, contribute the most to the binding of the cognate ligands, by hydrogen bonding and by base stacking, respectively. Mg2+ ions are essential in stabilizing the binding pocket. For the synthetic ligands, the hydrogen-bonding contributions of the U6-C15-A29 triplet are significantly compromised, and the bound state resembles the apo state in several respects, including the disengagement of the C15-A14-A13 and A32-G33 base stacks. The bulkier synthetic ligands lead to significantly loosening of the binding pocket, including extrusion of the C15 nucleobase and a widening of the C15-C30 groove. Enhanced-sampling simulations further revealed that the cognate and synthetic ligands unbind in almost opposite directions. Our work offers new insight for designing riboswitch ligands.

SUBMITTER: Hu G

PROVIDER: S-EPMC8612554 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Similar Datasets

Folding of a transcriptionally acting preQ1 riboswitch.

Project description:7-Aminomethyl-7-deazaguanine (preQ(1)) sensitive mRNA domains belong to the smallest riboswitches known to date. Although recent efforts have revealed the three-dimensional architecture of the ligand-aptamer complex less is known about the molecular details of the ligand-induced response mechanism that modulates gene expression. We present an in vitro investigation on the ligand-induced folding process of the preQ(1) responsive RNA element from Fusobacterium nucleatum using biophysical methods, including fluorescence and NMR spectroscopy of site-specifically labeled riboswitch variants. We provide evidence that the full-length riboswitch domain adopts two different coexisting stem-loop structures in the expression platform. Upon addition of preQ(1), the equilibrium of the competing hairpins is significantly shifted. This system therefore, represents a finely tunable antiterminator/terminator interplay that impacts the in vivo cellular response mechanism. A model is presented how a riboswitch that provides no obvious overlap between aptamer and terminator stem-loop solves this communication problem by involving bistable sequence determinants.

| S-EPMC2890745 | biostudies-literature

The free energy landscape of small molecule unbinding.

Project description:The spontaneous dissociation of six small ligands from the active site of FKBP (the FK506 binding protein) is investigated by explicit water molecular dynamics simulations and network analysis. The ligands have between four (dimethylsulphoxide) and eleven (5-diethylamino-2-pentanone) non-hydrogen atoms, and an affinity for FKBP ranging from 20 to 0.2 mM. The conformations of the FKBP/ligand complex saved along multiple trajectories (50 runs at 310 K for each ligand) are grouped according to a set of intermolecular distances into nodes of a network, and the direct transitions between them are the links. The network analysis reveals that the bound state consists of several subbasins, i.e., binding modes characterized by distinct intermolecular hydrogen bonds and hydrophobic contacts. The dissociation kinetics show a simple (i.e., single-exponential) time dependence because the unbinding barrier is much higher than the barriers between subbasins in the bound state. The unbinding transition state is made up of heterogeneous positions and orientations of the ligand in the FKBP active site, which correspond to multiple pathways of dissociation. For the six small ligands of FKBP, the weaker the binding affinity the closer to the bound state (along the intermolecular distance) are the transition state structures, which is a new manifestation of Hammond behavior. Experimental approaches to the study of fragment binding to proteins have limitations in temporal and spatial resolution. Our network analysis of the unbinding simulations of small inhibitors from an enzyme paints a clear picture of the free energy landscape (both thermodynamics and kinetics) of ligand unbinding.

| S-EPMC3033371 | biostudies-literature

Dynamics Correlation Network for Allosteric Switching of PreQ1 Riboswitch.

Project description:Riboswitches are a class of metabolism control elements mostly found in bacteria. Due to their fundamental importance in bacteria gene regulation, riboswitches have been proposed as antibacterial drug targets. Prequeuosine (preQ1) is the last free precursor in the biosynthetic pathway of queuosine that is crucial for translation efficiency and fidelity. However, the regulation mechanism for the preQ1 riboswitch remains unclear. Here we constructed fluctuation correlation network based on all-atom molecular dynamics simulations to reveal the regulation mechanism. The results suggest that the correlation network in the bound riboswitch is distinctly different from that in the apo riboswitch. The community network indicates that the information freely transfers from the binding site of preQ1 to the expression platform of the P3 helix in the bound riboswitch and the P3 helix is a bottleneck in the apo riboswitch. Thus, a hypothesis of "preQ1-binding induced allosteric switching" is proposed to link riboswitch and translation regulation. The community networks of mutants support this hypothesis. Finally, a possible allosteric pathway of A50-A51-A52-U10-A11-G12-G56 was also identified based on the shortest path algorithm and confirmed by mutations and network perturbation. The novel fluctuation network analysis method can be used as a general strategy in studies of riboswitch structure-function relationship.

| S-EPMC4971525 | biostudies-literature

Comparison of a preQ1 riboswitch aptamer in metabolite-bound and free states with implications for gene regulation.

Project description:Riboswitches are RNA regulatory elements that govern gene expression by recognition of small molecule ligands via a high affinity aptamer domain. Molecular recognition can lead to active or attenuated gene expression states by controlling accessibility to mRNA signals necessary for transcription or translation. Key areas of inquiry focus on how an aptamer attains specificity for its effector, the extent to which the aptamer folds prior to encountering its ligand, and how ligand binding alters expression signal accessibility. Here we present crystal structures of the preQ(1) riboswitch from Thermoanaerobacter tengcongensis in the preQ(1)-bound and free states. Although the mode of preQ(1) recognition is similar to that observed for preQ(0), surface plasmon resonance revealed an apparent K(D) of 2.1 ± 0.3 nm for preQ(1) but a value of 35.1 ± 6.1 nm for preQ(0). This difference can be accounted for by interactions between the preQ(1) methylamine and base G5 of the aptamer. To explore conformational states in the absence of metabolite, the free-state aptamer structure was determined. A14 from the ceiling of the ligand pocket shifts into the preQ(1)-binding site, resulting in "closed" access to the metabolite while simultaneously increasing exposure of the ribosome-binding site. Solution scattering data suggest that the free-state aptamer is compact, but the "closed" free-state crystal structure is inadequate to describe the solution scattering data. These observations are distinct from transcriptional preQ(1) riboswitches of the same class that exhibit strictly ligand-dependent folding. Implications for gene regulation are discussed.

| S-EPMC3137038 | biostudies-literature

Cooperative and directional folding of the preQ1 riboswitch aptamer domain.

Project description:Riboswitches are cis-acting RNA fragments that regulate gene expression by sensing cellular levels of the associated small metabolites. In bacteria, the class I preQ(1) riboswitch allows the fine-tuning of queuosine biosynthesis in response to the intracellular concentration of the queuosine anabolic intermediate preQ(1). When binding preQ(1), the aptamer domain undergoes a significant degree of secondary and tertiary structural rearrangement and folds into an H-type pseudoknot. Conformational "switching" of the riboswitch aptamer domain upon recognizing its cognate metabolite plays a key role in the regulatory mechanism of the preQ(1) riboswitch. We investigate the folding mechanism of the preQ(1) riboswitch aptamer domain using all-atom Go?-model simulations. The folding pathway of such a single domain is found to be cooperative and sequentially coordinated, as the folding proceeds in the 5' ? 3' direction. This kinetically efficient folding mechanism suggests a fast ligand-binding response in competition with RNA elongation.

| S-EPMC3109358 | biostudies-literature

Structural, functional, and taxonomic diversity of three preQ1 riboswitch classes.

Project description:Previously, two riboswitch classes have been identified that sense and respond to the hypermodified nucleobase called prequeuosine1 (preQ1). The enormous expansion of available genomic DNA sequence data creates new opportunities to identify additional representatives of the known riboswitch classes and to discover novel classes. We conducted bioinformatics searches on microbial genomic DNA data sets to discover numerous additional examples belonging to the two previously known riboswitch classes for preQ1 (classes preQ1-I and preQ1-II), including some structural variants that further restrict ligand specificity. Additionally, we discovered a third preQ1-binding riboswitch class (preQ1-III) that is structurally distinct from previously known classes. These findings demonstrate that numerous organisms monitor the concentrations of this modified nucleobase by exploiting one or more riboswitch classes for this widespread compound.

| S-EPMC4145258 | biostudies-literature

Molecular mechanism of preQ1 riboswitch action: a molecular dynamics study.

Project description:Riboswitches often occur in the 5'-untranslated regions of bacterial mRNA where they regulate gene expression. The preQ(1) riboswitch controls the biosynthesis of a hypermodified nucleoside queuosine in response to binding the queuosine metabolic intermediate. Structures of the ligand-bound and ligand-free states of the preQ(1) riboswitch from Thermoanaerobacter tengcongensis were determined recently by X-ray crystallography. We used multiple, microsecond-long molecular dynamics simulations (29 ?s in total) to characterize the structural dynamics of preQ(1) riboswitches in both states. We observed different stabilities of the stem in the bound and free states, resulting in different accessibilities of the ribosome-binding site. These differences are related to different stacking interactions between nucleotides of the stem and the associated loop, which itself adopts different conformations in the bound and free states. We suggest that the loop not only serves to bind preQ(1) but also transmits information about ligand binding from the ligand-binding pocket to the stem, which has implications for mRNA accessibility to the ribosome. We explain functional results obscured by a high salt crystallization medium and help to refine regions of disordered electron density, which demonstrates the predictive power of our approach. Besides investigating the functional dynamics of the riboswitch, we have also utilized this unique small folded RNA system for analysis of performance of the RNA force field on the ?s time scale. The latest AMBER parmbsc0?(OL3) RNA force field is capable of providing stable trajectories of the folded molecule on the ?s time scale. On the other hand, force fields that are not properly balanced lead to significant structural perturbations on the sub-?s time scale, which could easily lead to inappropriate interpretation of the simulation data.

| S-EPMC3505677 | biostudies-literature

Effect of mutations on binding of ligands to guanine riboswitch probed by free energy perturbation and molecular dynamics simulations.

Project description:Riboswitches can regulate gene expression by direct and specific interactions with ligands and have recently attracted interest as potential drug targets for antibacterial. In this work, molecular dynamics (MD) simulations, free energy perturbation (FEP) and molecular mechanics generalized Born surface area (MM-GBSA) methods were integrated to probe the effect of mutations on the binding of ligands to guanine riboswitch (GR). The results not only show that binding free energies predicted by FEP and MM-GBSA obtain an excellent correlation, but also indicate that mutations involved in the current study can strengthen the binding affinity of ligands GR. Residue-based free energy decomposition was applied to compute ligand-nucleotide interactions and the results suggest that mutations highly affect interactions of ligands with key nucleotides U22, U51 and C74. Dynamics analyses based on MD trajectories indicate that mutations not only regulate the structural flexibility but also change the internal motion modes of GR, especially for the structures J12, J23 and J31, which implies that the aptamer domain activity of GR is extremely plastic and thus readily tunable by nucleotide mutations. This study is expected to provide useful molecular basis and dynamics information for the understanding of the function of GR and possibility as potential drug targets for antibacterial.

| S-EPMC6649850 | biostudies-literature

Structural determinants for ligand capture by a class II preQ1 riboswitch.

Project description:Prequeuosine (preQ1) riboswitches are RNA regulatory elements located in the 5' UTR of genes involved in the biosynthesis and transport of preQ1, a precursor of the modified base queuosine universally found in four tRNAs. The preQ1 class II (preQ1-II) riboswitch regulates preQ1 biosynthesis at the translational level. We present the solution NMR structure and conformational dynamics of the 59 nucleotide Streptococcus pneumoniae preQ1-II riboswitch bound to preQ1. Unlike in the preQ1 class I (preQ1-I) riboswitch, divalent cations are required for high-affinity binding. The solution structure is an unusual H-type pseudoknot featuring a P4 hairpin embedded in loop 3, which forms a three-way junction with the other two stems. (13)C relaxation and residual dipolar coupling experiments revealed interhelical flexibility of P4. We found that the P4 helix and flanking adenine residues play crucial and unexpected roles in controlling pseudoknot formation and, in turn, sequestering the Shine-Dalgarno sequence. Aided by divalent cations, P4 is poised to act as a "screw cap" on preQ1 recognition to block ligand exit and stabilize the binding pocket. Comparison of preQ1-I and preQ1-II riboswitch structures reveals that whereas both form H-type pseudoknots and recognize preQ1 using one A, C, or U nucleotide from each of three loops, these nucleotides interact with preQ1 differently, with preQ1 inserting into different grooves. Our studies show that the preQ1-II riboswitch uses an unusual mechanism to harness exquisite control over queuosine metabolism.

| S-EPMC3926045 | biostudies-literature

Molecular mechanism for preQ1-II riboswitch function revealed by molecular dynamics.

Project description:Riboswitches are RNA molecules that regulate gene expression using conformational change, affected by binding of small molecule ligands. A crystal structure of a ligand-bound class II preQ1 riboswitch has been determined in a previous structural study. To gain insight into the dynamics of this riboswitch in solution, eight total molecular dynamic simulations, four with and four without ligand, were performed using the Amber force field. In the presence of ligand, all four of the simulations demonstrated rearranged base pairs at the 3' end, consistent with expected base-pairing from comparative sequence analysis in a prior bioinformatic analysis; this suggests the pairing in this region was altered by crystallization. Additionally, in the absence of ligand, three of the simulations demonstrated similar changes in base-pairing at the ligand binding site. Significantly, although most of the riboswitch architecture remained intact in the respective trajectories, the P3 stem was destabilized in the ligand-free simulations in a way that exposed the Shine-Dalgarno sequence. This work illustrates how destabilization of two major groove base triples can influence a nearby H-type pseudoknot and provides a mechanism for control of gene expression by a fold that is frequently found in bacterial riboswitches.

| S-EPMC4604430 | biostudies-literature

OmicsDI is part of the ELIXIR infrastructure

OmicsDI is an Elixir interoperability service. Learn more ›

Tweets

OmicsDI Databases

PRIDE
PeptideAtlas
MassIVE
JPOST Repository
Physiome Model Repository

EGA
EVA
ENA
LINCS
PAXDB
Cell Collective

MetaboLights
Metabolomics Workbench
MetabolomeExpress
GNPS
BioModels
FAIRDOMHub

ArrayExpress
dbGaP
ExpressionAtlas
GEO
NODE

Information

Databases
Help
API
Contact us
Code on GitHub
Terms of use
Submit Data