Restoration of Ribozyme Tertiary Contact and Function by Using a Molecular Glue for RNA

Some RNA classes require folding into the proper higher‐order structures to exert their functions. Hammerhead ribozyme (HHR) requires a folding conformation stabilized by tertiary interaction for full activity. A rationally engineered HHR was developed that was inactive, but could be activated by a synthetic RNA‐binding ligand, naphthyridine carbamate tetramer with Z‐stilbene linker (Z‐NCTS). Binding of Z‐NCTS could induce the formation of an active folding structure and thereby restore ribozyme activity, where Z‐NCTS acts as a molecular glue to bring two isolated RNA loops into contact with each other. Next, we designed a Z‐NCTS‐responsive genetic switch using the HHR sequence inserted into the 3′ untranslated region as a cis‐acting element. We demonstrated that the rationally designed ribozyme switch enabled regulation of gene expression by Z‐NCTS and was functional in mammalian cells.     

A Crystal Structure of a Functional RNA Molecule Containing an Artificial Nucleobase Pair

As one of its goals, synthetic biology seeks to increase the number of building blocks in nucleic acids. While efforts towards this goal are well advanced for DNA, they have hardly begun for RNA. Herein, we present a crystal structure for an RNA riboswitch where a stem C:G pair has been replaced by a pair between two components of an artificially expanded genetic‐information system (AEGIS), Z and P , (6‐amino‐5‐nitro‐2(1H)‐pyridone and 2‐amino‐imidazo[1,2‐a]‐1,3,5‐triazin‐4‐(8H)‐one). The structure shows that the Z : P pair does not greatly change the conformation of the RNA molecule nor the details of its interaction with a hypoxanthine ligand. This was confirmed in solution by in‐line probing, which also measured a 3.7 nM affinity of the riboswitch for guanine. These data show that the Z : P pair mimics the natural Watson–Crick geometry in RNA in the first example of a crystal structure of an RNA molecule that contains an orthogonal added nucleobase pair.     

Pseudoknot preorganization of the preQ1 class I riboswitch

To explore folding and ligand recognition of metabolite-responsive RNAs is of major importance to comprehend gene regulation by riboswitches. Here, we demonstrate, using NMR spectroscopy, that the free aptamer of a preQ(1) class I riboswitch preorganizes into a pseudoknot fold in a temperature- and Mg(2+)-dependent manner. The preformed pseudoknot represents a structure that is close to the ligand-bound state and that likely represents the conformation selected by the ligand. Importantly, a defined base pair mutation within the pseudoknot interaction stipulates whether, in the absence of ligand, dimer formation of the aptamer competes with intramolecular pseudoknot formation. This study pinpoints how RNA preorganization is a crucial determinant for the adaptive recognition process of RNA and ligand.     

Binding to an RNA aptamer changes the charge distribution and conformation of malachite green

RNA plays a central role in many biological processes and is therefore an important target for drug development. In recent years an increasing wealth of structural and functional information about RNA-ligand complexes has been obtained using in vitro selected RNAs (aptamers). However, all those studies focused on structure and changes of the nucleic acid and mostly considered the ligand as a rigid target. To develop a detailed picture of ligand structure and dynamics in RNA-small molecule complexes, the malachite green binding aptamer was studied. Isotopically labeled ligand in complex with RNA was analyzed by NMR spectroscopy in solution. The surprisingly asymmetric changes in the (13)C chemical shift of the ligand methyl groups indicate that the dye undergoes changes in its conformation and charge distribution upon binding. The role of the RNA electrostatic field in this interaction was explored using ab initio calculations of the ligand structure and charge distribution. The results indicate that the uneven charge distribution in the RNA binding pocket provides a major contribution to the driving force of the ligand structural changes. The observation that not only the RNA adapts to the ligand, in what is called adaptive binding, but that the ligand itself also undergoes conformational changes ("induced fit") is crucial for the rational design of RNA ligands and for understanding the properties of RNA-ligand complexes.     

Intrinsic contribution of the 2'-hydroxyl to RNA conformational heterogeneity

Canonical duplex RNA assumes only the A-form conformation at the secondary structure level while, in contrast, a wide range of noncanonical, tertiary conformations of RNA occur. Here, we show how the 2'-hydroxyl controls RNA conformational properties. Quantum mechanical calculations reveal that the orientation of the 2'-hydroxyl significantly alters the intrinsic flexibility of the phosphodiester backbone, favoring the A-form in duplex RNA when it is in the base orientation and facilitating sampling of a wide range of noncanonical, tertiary structures when it is in the O3' orientation. Influencing the orientation of the 2'-hydroxyl are interactions with the environment, as evidenced by crystallographic survey data, indicating the 2'-hydroxyl to sample more of the O3' orientation in noncanonical RNA structures. These results indicate that the 2'-hydroxyl acts as a "switch", both limiting the conformation of RNA to the A-form at the secondary structure level and allowing RNA to sample a wide range of noncanonical tertiary conformations.     

Formation of a Ligand‐Assisted Complex of Two RNA Hairpin Loops

The hairpin structure is one of the most common secondary structures in RNA and holds a central position in the stream of RNA folding from a non‐structured RNA to structurally complex and functional ribonucleoproteins. Since the RNA secondary structure is strongly correlated to the function and can be modulated by the binding of small molecules, we have investigated the modulation of RNA folding by a ligand‐assisted formation of loop–loop complexes of two RNA hairpin loops. With a ligand (NCT6), designed based on the ligand binding to the G–G mismatches in double‐stranded DNA, we successfully demonstrated the formation of both inter‐ and intra‐molecular NCT6‐assisted complex of two RNA hairpin loops. NCT6 selectively bound to the two hairpin loops containing (CGG)₃ in the loop region. Native polyacrylamide gel electrophoresis analysis of two doubly‐labeled RNA hairpin loops clearly showed the formation of intermolecular NCT6‐assisted loop–loop complex. Förster resonance energy‐transfer studies of RNA constructs containing two hairpin loops, in which each hairpin was labeled with Alexa488 and Cy3 fluorophores, showed the conformational change of the RNA constructs upon binding of NCT6. These experimental data showed that NCT6 simultaneously bound to two hairpin RNAs at the loop region, and can induce the conformational change of the RNA molecule. These data strongly support that NCT6 functions as molecular glue for two hairpin RNAs.     

Structure of the Ribosomal RNA Decoding Site Containing a Selenium‐Modified Responsive Fluorescent Ribonucleoside Probe

Comprehensive understanding of the structure–function relationship of RNA both in real time and at atomic level will have a profound impact in advancing our understanding of RNA functions in biology. Here, we describe the first example of a multifunctional nucleoside probe, containing a conformation‐sensitive fluorophore and an anomalous X‐ray diffraction label (5‐selenophene uracil), which enables the correlation of RNA conformation and recognition under equilibrium and in 3D. The probe incorporated into the bacterial ribosomal RNA decoding site, fluorescently reports antibiotic binding and provides diffraction information in determining the structure without distorting native RNA fold. Further, by comparing solution binding data and crystal structure, we gained insight on how the probe senses ligand‐induced conformational change in RNA. Taken together, our nucleoside probe represents a new class of biophysical tool that would complement available tools for functional RNA investigations.     

DNA vibrational coupling revealed with two-dimensional infrared spectroscopy: insight into why vibrational spectroscopy is sensitive to DNA structure

Two-dimensional infrared (2D IR) spectroscopy was used to study the carbonyl vibrational modes of guanine and cytosine bases in A- and B-form DNA. Located between 1600 and 1700 cm(-1), these modes are often used to monitor DNA secondary structure with traditional infrared spectroscopies such as FTIR, but traditional spectroscopies lack the necessary observables to unravel the coupling mechanisms that make these modes sensitive to secondary structure. By using 2D IR spectroscopy and electronic structure calculations on d(G(5)C(5)) and d(GC)(8) model nucleic acids, we find that hydrogen-bonded guanine/cytosine base pairs are primarily electrostatically coupled and that the coupling between these modes can be modeled with a transition dipole density approach. In comparison, electrostatics is insufficient to model stacked bases because of cooperative charge-sharing effects, but the coupling can be accurately calculated using a finite difference method. We find that the coupling is very strong for both hydrogen-bonded and stacked base geometries, creating vibrational modes that extend both across the base pairs and along the lengths of the helices. Our results provide a physical basis for understanding how strong coupling gives rise to the empirically established relationship between infrared spectroscopy and DNA/RNA secondary structure.     

A molecular mechanics/grid method for evaluation of ligand–receptor interactions

We present a computational method for prediction of the conformation of a ligand when bound to a macromolecular receptor. The method is intended for use in systems in which the approximate location of the binding site is known and no large-scale rearrangements of the receptor are expected upon formation of the complex. The ligand is initially placed in the vicinity of the binding site and the atomic motions of the ligand and binding site are explicitly simulated, with solvent represented by an implicit solvation model and using a grid representation for the bulk of the receptor protein. These two approximations make the method computationally efficient and yet maintain accuracy close to that of an all-atom calculation. For the benzamidine/trypsin system, we ran 100 independent simulations, in many of which the ligand settled into the low-energy conformation observed in the crystal structure of the complex. The energy of these conformations was lower than and well-separated from that of others sampled. Extensions of this method are also discussed. © 1995 by John Wiley & Sons, Inc.     

RNA architecture dictates the conformations of a bound peptide.

BACKGROUND: The biological function of several viral and bacteriophage proteins, and their arginine-rich subdomains, involves RNA-mediated interactions. It has been shown recently that bound peptides adopt either beta-hairpin or alpha-helical conformations in viral and phage peptide-RNA complexes. We have compared the structures of the arginine-rich peptide domain of HIV-1 Rev bound to two RNA aptamers to determine whether RNA architecture can dictate the conformations of a bound peptide. RESULTS: The core-binding segment of the HIV-1 Rev peptide class II RNA aptamer complex spans the two-base bulge and hairpin loop of the bound RNA and the carboxy-terminal segment of the bound peptide. The bound peptide is anchored in place by backbone and sidechain intermolecular hydrogen bonding and van der Waals stacking interactions. One of the bulge bases participates in U*(A*U) base triple formation, whereas the other is looped out and flaps over the bound peptide in the complex. The seven-residue hairpin loop is closed by a sheared G*A mismatch pair with several pyrimidines looped out of the hairpin fold. CONCLUSIONS: Our structural studies establish that RNA architecture dictates whether the same HIV-1 Rev peptide folds into an extended or alpha-helical conformation on complex formation. Arginine-rich peptides can therefore adapt distinct secondary folds to complement the tertiary folds of their RNA targets. This contrasts with protein-RNA complexes in which elements of RNA secondary structure adapt to fit within the tertiary folds of their protein targets.     

Fingerprinting Noncanonical and Tertiary RNA Structures by Differential SHAPE Reactivity

Many RNA structures are composed of simple secondary structure elements linked by a few critical tertiary interactions. SHAPE chemistry has made interrogation of RNA dynamics at single-nucleotide resolution straightforward. However, de novo identification of nucleotides involved in tertiary interactions remains a challenge. Here we show that nucleotides that form noncanonical or tertiary contacts can be detected by comparing information obtained using two SHAPE reagents, N-methylisatoic anhydride (NMIA) and 1-methyl-6-nitroisatoic anhydride (1M6). Nucleotides that react preferentially with NMIA exhibit slow local nucleotide dynamics and usually adopt the less common C2'-endo ribose conformation. Experiments and first-principles calculations show that 1M6 reacts preferentially with nucleotides in which one face of the nucleobase allows an unhindered stacking interaction with the reagent. Differential SHAPE reactivities were used to detect noncanonical and tertiary interactions in four RNAs with diverse structures and to identify preformed noncanonical interactions in partially folded RNAs. Differential SHAPE reactivity analysis will enable experimentally concise, large-scale identification of tertiary structure elements and ligand binding sites in complex RNAs and in diverse biological environments.     

A quantitative reconstruction of the amide I contour in the IR spectra of globular proteins: from structure to spectrum

The Amide I contours of six globular proteins of varied secondary structure content along with a peptide model for collagen and pulmonary surfactant protein C have been simulated very closely by using a modified GF matrix method. The starting point for the method uses the three-dimensional structure as obtained from the Protein Data Bank. Elements of the interactions between peptide groups (e.g., transition dipole coupling) are very sensitive to tertiary structure, thus the current formalism demonstrates that the Amide I contour may be useful for a more detailed probe of 3-D conformation that goes beyond the traditional use of this band to probe the percentages of particular elements of secondary structure. For example, postulated changes to a known structure can be tested by comparing the new simulated band to the experimental band. A number of refinements to the transition dipole interaction calculation have been made. Most of the important interactions between the C=O oscillators that define the Amide I mode appear to have been identified, including through space transition dipole coupling, through valence bond and through hydrogen bond coupling. The eigenvector matrix produced by the method permits the contribution of each peptide group to the spectrum to be precisely determined. Analysis of the results shows that the often-used structure-frequency correlations are at best approximate and at worst misleading. The subbands from helices, sheets, turns, and loops are much broader and more overlapped than has been commonly assumed. Furthermore, the traditional alpha-helical marker band may be substantially distorted in short segments. Difference spectra based on isotope editing, a technique thought capable of revealing the spectral contributions of individual peptide groups, are shown to be prone to misinterpretation.     

An Optimized Ni-catalyzed Chan-Lam Type Coupling: Enantioretentive Access to Chiral N-Aryl Sulfinamides

C−N bond formation takes on a critical significance in reactions of organic synthesis, material production and pharmaceutical manufactory. Chan-Lam has proposed a useful methodology to furnish secondary arylamides under mild conditions. However, when chiral sulfinamides serve as the coupling precursors, the Cu-catalyzed coupling reaction is found with low efficacy. Complex side-products are generated under classic conditions. Moreover, it led to the racemization of the coupling product. In this study, an optimized Ni-catalyzed Chan-Lam type coupling conditions were proposed, which resulted in clean conversion from chiral sulfinamides and arylboronic acids to offer N-aryl sulfinamides efficiently and enantioretentively. The trans-N¹,N²-dimethylcyclohexane-1,2-diamine was proven as the most efficient ligand. Under the optimized conditions, a series of chiral N-aryl sulfinamides was prepared with high chemical yield without racemization. Furthermore, a plausible and novel mechanism was proposed. Interestingly, the method could efficiently furnish a wide variety of C−X bonds by coupling arylboronic acids with different nucleophiles.     

NMR structure of the peptidyl transferase RNA inhibitor antibiotic amicetin

In this communication, we report the solution state NMR structure determination of the peptidyl transferase RNA inhibitor antibiotic amicetin. We have successfully characterised the NMR spectrum of amicetin using a range of homo- and heteronuclear NMR techniques. Using experimental ROE-based distance and 1H--1H scalar coupling derived dihedral angle geometrical constraints as input into the three-dimensional structure determination protocol, we have generated an energy-minimised average structure of the antibiotic. Amicetin adopts a stable well-folded conformation in solution, mediated by a network of hydrogen bonds caused by proton donor and acceptor groups at either end of the molecule. The NMR structure of amicetin shows that the cytosine moiety occupies the critical turn position within the fold, which may be structurally significant for interaction with peptidyl transferase ribosomal RNA. The structure is distinctly different from the published X-ray crystal structure of amicetin in which it adopts a linear, extended chain-like conformation with a number of intermolecular hydrogen bonds. In addition to structure, we have probed the dynamics of amicetin in solution and have observed retarded exchange of the amide proton involved in folding. We have also characterised the ionisation properties of amicetin by carrying out NMR pH titration and measuring the pKa of the primary and tertiary amino groups, 7.27 and 7.52, respectively, which are in agreement with the reported values in literature. Solving the NMR structure of amicetin provides a valuable opportunity to determine the structure of its complex with RNA in solution state.     

GalaxyDock3: Protein–ligand docking that considers the full ligand conformational flexibility

Predicting conformational changes of both the protein and the ligand is a major challenge when a protein–ligand complex structure is predicted from the unbound protein and ligand structures. Herein, we introduce a new protein–ligand docking program called GalaxyDock3 that considers the full ligand conformational flexibility by explicitly sampling the ligand ring conformation and allowing the relaxation of the full ligand degrees of freedom, including bond angles and lengths. This method is based on the previous version (GalaxyDock2) which performs the global optimization of a designed score function. Ligand ring conformation is sampled from a ring conformation library constructed from structure databases. The GalaxyDock3 score function was trained with an additional bonded energy term for the ligand on a large set of complex structures. The performance of GalaxyDock3 was improved compared to GalaxyDock2 when predicted ligand conformation was used as the input for docking, especially when the input ligand conformation differs significantly from the crystal conformation. GalaxyDock3 also compared favorably with other available docking programs on two benchmark tests that contained diverse ligand rings. The program is freely available at http://galaxy.seoklab.org/softwares/galaxydock.html . © 2019 Wiley Periodicals, Inc.     

Mutual induced-fit controlled signal amplification of hydrogen-bonded capsule formation

Mutual induced-fit is well known process in biological system and playing crucial role in enzymatic catalysis, peptide-bond synthesis, protein synthesis, etc…, proceeds by mutual interactions of two or more than two components. However, in an artificial system, to mimic such process is challenging, especially by purely organic components. In this report, we have shown importance of mutual induced-fit in signal amplification of hydrogen-bonded capsule formation. Two different types of highly flexible ligands (one N-bridged and another triazine-bridged) are used. N-bridged tripodal ligand acting as clip over the adduct obtained from triazine-bridged ligand, bringing the resultant adduct in the cone-shape conformation. The resulting cone-shape conformation undergo solvent polarity dependent hydrogen-bonded capsule formation as a sole product. In absence of N-bridged tripodal ligand, only 50% capsule formation is observed through ¹H NMR at 100 mM concentration. However, in the presence of N-bridged tripodal ligand, organic components undergo mutual interactions to opt the cone-shape conformation and 100% capsule is formed (independent of concentration). The complete process is characterized by IR-spectra, ¹H and ¹³C NMR spectra, concentration dependent ¹H NMR titration spectra, ¹H-¹H COSY, ¹H-¹H NOESY, DOSY NMR, high resolution ESI mass spectra, and also by energy-minimized structure.     

The MPSim-Dock hierarchical docking algorithm: application to the eight trypsin inhibitor cocrystals

To help improve the accuracy of protein-ligand docking as a useful tool for drug discovery, we developed MPSim-Dock, which ensures a comprehensive sampling of diverse families of ligand conformations in the binding region followed by an enrichment of the good energy scoring families so that the energy scores of the sampled conformations can be reliably used to select the best conformation of the ligand. This combines elements of DOCK4.0 with molecular dynamics (MD) methods available in the software, MPSim. We test here the efficacy of MPSim-Dock to predict the 64 protein-ligand combinations formed by starting with eight trypsin cocrystals, and crossdocking the other seven ligands to each protein conformation. We consider this as a model for how well the method would work for one given target protein structure. Using as a criterion that the structures within 2 kcal/mol of the top scoring include a conformation within a coordinate root mean square (CRMS) of 1 A of the crystal structure, we find that 100% of the 64 cases are predicted correctly. This indicates that MPSim-Dock can be used reliably to identify strongly binding ligands, making it useful for virtual ligand screening.     

Terminal‐Selective Functionalization of Alkyl Chains by Regioconvergent Cross‐Coupling

Hydrocarbons are still the most important precursors of functionalized organic molecules, which has stirred interest in the discovery of new C?H bond functionalization methods. We describe herein a new step‐economical approach that enables C?C bonds to be constructed at the terminal position of linear alkanes. First, we show that secondary alkyl bromides can undergo in situ conversion into alkyl zinc bromides and regioconvergent Negishi coupling with aryl or alkenyl triflates. The use of a suitable phosphine ligand favoring Pd migration enabled the selective formation of the linear cross‐coupling product. Subsequently, mixtures of secondary alkyl bromides were prepared from linear alkanes by standard bromination, and regioconvergent cross‐coupling then provided access to the corresponding linear arylation product in only two steps.     

MD simulations of homomorphous PNA, DNA, and RNA single strands: characterization and comparison of conformations and dynamics

MD simulations of homomorphous single-stranded PNA, DNA, and RNA with the same base sequence have been performed in aqueous solvent. For each strand two separate simulations were performed starting from a (i) helical conformation and (ii) random coiled state. Comparisons of the simulations with the single-stranded helices (case i) show that the differences in the covalent nature of the backbones cause significant differences in the structural and dynamical properties of the strands. It is found that the PNA strand maintains its nice base-stacked initial helical structure throughout the 1.5-ns MD simulation at 300 K, while DNA/RNA show relatively larger fluctuations in the structures with a few local unstacking events during -ns MD simulation each. It seems that the weak physical coupling between the bases and the backbone in PNA causes a loss of correlation between the dynamics of the bases and the backbone compared to the DNA/RNA and helps maintain the base-stacked helical conformation. The global flexibility of a single-stranded PNA helix was also found to be lowest, while RNA appears to be the most flexible single-stranded helix. The sugar pucker of several nucleotides in single-stranded DNA and RNA was found to adopt both C2'-endo and C3'-endo conformations for significant times. This effect is more pronounced for single strands in completely coiled states. The simulations with single-stranded coils as the initial structure also indicate that a PNA can adopt a more compact globular structure, while DNA/RNA of the same size adopts a more extended coil structure. This allows even a short PNA in the coiled state to form a significantly stable nonsequentially base-stacked globular structure in solution. Due to the hydrophobic nature of the PNA backbone, it interacts with surrounding water rather weakly compared to DNA/RNA.