The conformation of the sarcin/ricin loop from 28S ribosomal RNA

RNAMotifProfile: a graph-based approach to build RNA structural motif profiles

RNA structural motifs are the recurrent segments in RNA three-dimensional structures that play a crucial role in the functional diversity of RNAs. Understanding the similarities and variations within these recurrent motif groups is essential for gaining insights into RNA structure and function. While recurrent structural motifs are generally assumed to be composed of the same isosteric base interactions, this consistent pattern is not observed across all examples of these motifs. Existing methods for analyzing and comparing RNA structural motifs may overlook variations in base interactions and associated nucleotides. RNAMotifProfile is a novel profile-to-profile alignment algorithm that generates a comprehensive profile from a group of structural motifs, incorporating all base interactions and associated nucleotides at each position. By structurally aligning input motif instances using a guide-tree-based approach, RNAMotifProfile captures the similarities and variations within recurrent motif groups. Additionally, RNAMotifProfile can function as a motif search tool, enabling the identification of instances of a specific motif family by searching with the corresponding profile. The ability to generate accurate and comprehensive profiles for RNA structural motif families, and to search for these motifs, facilitates a deeper understanding of RNA structure-function relationships and potential applications in RNA engineering and therapeutic design.

A dual ribosomal system in the zebrafish soma and germline

Protein synthesis during vertebrate embryogenesis is driven by ribosomes of two distinct origins: maternal ribosomes synthesized during oogenesis and stored in the egg, and somatic ribosomes, produced by the developing embryo after zygotic genome activation (ZGA). In zebrafish, these two ribosome types are expressed from different genomic loci and also differ in their ribosomal RNA (rRNA) sequence. To characterize this dual ribosome system further, we examined the expression patterns of maternal and somatic rRNAs during embryogenesis and in adult tissues. We found that maternal rRNAs are not only expressed during oogenesis but are continuously produced in the zebrafish germline. Proteomic analyses of maternal and somatic ribosomes unveiled differences in core ribosomal protein composition. Most nucleotide differences between maternal and somatic rRNAs are located in the flexible, structurally not resolved expansion segments. Our in vivo data demonstrated that both maternal and somatic ribosomes can be translationally active in the embryo. Using transgenically tagged maternal or somatic ribosome subunits, we experimentally confirm the presence of hybrid 80S ribosomes composed of 40S and 60S subunits from both origins and demonstrate the preferential in vivo association of maternal ribosomes with germline-specific transcripts. Our study identifies a distinct type of ribosomes in the zebrafish germline and thus presents a foundation for future explorations into possible regulatory mechanisms and functional roles of heterogeneous ribosomes.

1H, 13C, and 15N backbone and methyl group resonance assignments of ricin toxin A subunit

Ricin is a potent plant toxin that targets the eukaryotic ribosome by depurinating an adenine from the sarcin-ricin loop (SRL), a highly conserved stem-loop of the rRNA. As a category-B agent for bioterrorism it is a prime target for therapeutic intervention with antibodies and enzyme blocking inhibitors since no effective therapy exists for ricin. Ricin toxin A subunit (RTA) depurinates the SRL by binding to the P-stalk proteins at a remote site. Stimulation of the N-glycosidase activity of RTA by the P-stalk proteins has been studied extensively by biochemical methods and by X-ray crystallography. The current understanding of RTA's depurination mechanism relies exclusively on X-ray structures of the enzyme in the free state and complexed with transition state analogues. To date we have sparse evidence of conformational dynamics and allosteric regulation of RTA activity that can be exploited in the rational design of inhibitors. Thus, our primary goal here is to apply solution NMR techniques to probe the residue specific structural and dynamic coupling active in RTA as a prerequisite to understand the functional implications of an allosteric network. In this report we present de novo sequence specific amide and sidechain methyl chemical shift assignments of the 267 residue RTA in the free state and in complex with an 11-residue peptide (P11) representing the identical C-terminal sequence of the ribosomal P-stalk proteins. These assignments will facilitate future studies detailing the propagation of binding induced conformational changes in RTA complexed with inhibitors, antibodies, and biologically relevant targets.

Reducing the Immunogenicity of Pulchellin A-Chain, Ribosome-Inactivating Protein Type 2, by Computational Protein Engineering for Potential New Immunotoxins

Pulchellin is a plant biotoxin categorized as a type 2 ribosome-inactivating protein (RIPs) which potentially kills cells at very low concentrations. Biotoxins serve as targeting immunotoxins (IT), consisting of antibodies conjugated to toxins. ITs have two independent protein components, a human antibody and a toxin with a bacterial or plant source; therefore, they pose unique setbacks in immunogenicity. To overcome this issue, the engineering of epitopes is one of the beneficial methods to elicit an immunological response. Here, we predicted the tertiary structure of the pulchellin A-chain (PAC) using five common powerful servers and adopted the best model after refining. Then, predicted structure using four distinct computational approaches identified conformational B-cell epitopes. This approach identified some amino acids as a potential for lowering immunogenicity by point mutation. All mutations were then applied to generate a model of pulchellin containing all mutations (so-called PAM). Mutants’ immunogenicity was assessed and compared to the wild type as well as other mutant characteristics, including stability and compactness, were computationally examined in addition to immunogenicity. The findings revealed a reduction in immunogenicity in all mutants and significantly in N146V and R149A. Furthermore, all mutants demonstrated remarkable stability and validity in Molecular Dynamic (MD) simulations. During docking and simulations, the most homologous toxin to pulchellin, Abrin-A was applied as a control. In addition, the toxin candidate containing all mutations (PAM) disclosed a high level of stability, making it a potential model for experimental deployment. In conclusion, by eliminating B-cell epitopes, our computational approach provides a potential less immunogenic IT based on PAC.

RNA Probes for Visualization of Sarcin/ricin Loop Depurination without Background Fluorescence

Protein synthesis via ribosomes is a fundamental process in all known living organisms. However, it can be completely stalled by removing a single nucleobase (depurination) at the sarcin/ricin loop of the ribosomal RNA. In this work, we describe the preparation and optimization process of a fluorescent probe that can be used to visualize depurination. Starting from a fluorescent thiophene nucleobase analog, various RNA probes that fluoresce exclusively in the presence of a depurinated sarcin/ricin-loop RNA were designed and characterized. The main challenge in this process was to obtain a high fluorescence signal in the hybridized state with an abasic RNA strand, while keeping the background fluorescence low. With our new RNA probes, the fluorescence intensity and lifetime can be used for efficient monitoring of depurinated RNA.

Conformational stability of ageritin, a metal binding ribotoxin-like protein of fungal origin

Ageritin is a ribotoxin-like protein of biotechnological interest, belonging to a family of ribonucleases from edible mushrooms. Its enzymatic activity is explicated through the hydrolysis of a single phosphodiester bond, located in the sarcin/ricin loop of ribosomes. Unlike other ribotoxins, ageritin activity requires divalent cations (Zn²⁺). Here we investigated the conformational stability of ageritin in the pH range 4.0-7.4, using calorimetric and spectroscopic techniques. We observed a high protein thermal stability at all pHs with a denaturation temperature of 78 °C. At pH 5.0 we calculated a value of 36 kJ mol^-1 for the unfolding Gibbs energy at 25 °C. We also analysed the thermodynamic and catalytic behaviour of S-pyridylethylated form, obtained by alkylating the single Cys18 residue, which is predicted to bind Zn²⁺. We show that this form possesses the same activity and structure of ageritin, but lower stability. In fact, the corresponding values of 52 °C and 14 kJ mol^-1 were found. Conservation of activity is consistent with the location of alkylation site on the opposite site of the catalytic site cleft. Inasmuch as Cys18 is part of a structurally stabilizing zinc-binding site, disrupted by cysteine alkylation, our results point to an important role of metal ions in ageritin stability.

Cryo neutron crystallography demonstrates influence of RNA 2'-OH orientation on conformation, sugar pucker and water structure

The ribose 2′-hydroxyl is the key chemical difference between RNA and DNA and primary source of their divergent structural and functional characteristics. Macromolecular X-ray diffraction experiments typically do not reveal the positions of hydrogen atoms. Thus, standard crystallography cannot determine 2′-OH orientation (H2′-C2′-O2′-HO2′ torsion angle) and its potential roles in sculpting the RNA backbone and the expansive fold space. Here, we report the first neutron crystal structure of an RNA, the Escherichia coli rRNA Sarcin-Ricin Loop (SRL). 2′-OD orientations were established for all 27 residues and revealed O-D bonds pointing toward backbone (O3′, 13 observations), nucleobase (11) or sugar (3). Most riboses in the SRL stem region show a 2′-OD backbone-orientation. GAGA-tetraloop riboses display a 2′-OD base-orientation. An atypical C2′-endo sugar pucker is strictly correlated with a 2′-OD sugar-orientation. Neutrons reveal the strong preference of the 2′-OH to donate in H-bonds and that 2′-OH orientation affects both backbone geometry and ribose pucker. We discuss 2′-OH and water molecule orientations in the SRL neutron structure and compare with results from a solution phase 10 μs MD simulation. We demonstrate that joint cryo-neutron/X-ray crystallography offers an all-in-one approach to determine the complete structural properties of RNA, i.e. geometry, conformation, protonation state and hydration structure.

Genetic and Functional Analyses of Archaeal ATP-Dependent RNA Ligase in C/D Box sRNA Circularization and Ribosomal RNA Processing

RNA ligases play important roles in repairing and circularizing RNAs post-transcriptionally. In this study, we generated an allelic knockout of ATP-dependent RNA ligase (Rnl) in the hyperthermophilic archaeon Thermococcus kodakarensis to identify its biological targets. A comparative analysis of circular RNA reveals that the Rnl-knockout strain represses circularization of C/D box sRNAs without affecting the circularization of tRNA and rRNA processing intermediates. Recombinant archaeal Rnl could circularize C/D box sRNAs with a mutation in the conserved C/D box sequence element but not when the terminal stem structures were disrupted, suggesting that proximity of the two ends could be critical for intramolecular ligation. Furthermore, T. kodakarensis accumulates aberrant RNA fragments derived from ribosomal RNA in the absence of Rnl. These results suggest that Rnl is responsible for C/D box sRNA circularization and may also play a role in ribosomal RNA processing.

RNAMotifContrast: a method to discover and visualize RNA structural motif subfamilies

Understanding the 3D structural properties of RNAs will play a critical role in identifying their functional characteristics and designing new RNAs for RNA-based therapeutics and nanotechnology. While several existing computational methods can help in the analysis of RNA properties by recognizing structural motifs, they do not provide the means to compare and contrast those motifs extensively. We have developed a new method, RNAMotifContrast, which focuses on analyzing the similarities and variations of RNA structural motif characteristics. In this method, a graph is formed to represent the similarities among motifs, and a new traversal algorithm is applied to generate visualizations of their structural properties. Analyzing the structural features among motifs, we have recognized and generalized the concept of motif subfamilies. To asses its effectiveness, we have applied RNAMotifContrast on a dataset of known RNA structural motif families. From the results, we observed that the derived subfamilies possess unique structural variations while holding standard features of the families. Overall, the visualization approach of this method presents a new perspective to observe the relation among motifs more closely, and the discovered subfamilies provide opportunities to achieve valuable insights into RNA’s diverse roles.

Toxicity and membrane perturbation properties of the ribotoxin-like protein Ageritin

Ageritin is the prototype of a new ribotoxin-like protein family, which has been recently identified also in basidiomycetes. The protein exhibits specific RNase activity through the cleavage of a single phosphodiester bond located at sarcin/ricin loop of the large rRNA, thus inhibiting protein biosynthesis at early stages. Conversely to other ribotoxins, its activity requires the presence of divalent cations. In the present study, we report the activity of Ageritin on both prokaryotic and eukaryotic cells showing that the protein has a prominent effect on cancer cells viability and no effects on eukaryotic and bacterial cells. In order to rationalize these findings, the ability of the protein to interact with various liposomes mimicking normal, cancer and bacterial cell membranes was explored. The collected results indicate that Ageritin can interact with DPPC/DPPS/Chol vesicles, used as a model of cancer cell membranes, and with DPPC/DPPG vesicles, used as a model of bacterial cell membranes, suggesting a selective interaction with anionic lipids. However, a different perturbation of the two model membranes, mediated by cholesterol redistribution, was observed and this might be at the basis of Ageritin selective toxicity towards cancer cells.

Modeling and Predicting RNA Three-Dimensional Structures

Modeling the three-dimensional structure of RNAs is a milestone toward better understanding and prediction of nucleic acids molecular functions. Physics-based approaches and molecular dynamics simulations are not tractable on large molecules with all-atom models. To address this issue, coarse-grained models of RNA three-dimensional structures have been developed. In this chapter, we describe a graphical modeling based on the Leontis-Westhof extended base pair classification. This representation of RNA structures enables us to identify highly conserved structural motifs with complex nucleotide interactions in structure databases. We show how to take advantage of this knowledge to quickly predict three-dimensional structures of large RNA molecules and present the RNA-MoIP web server (http://rnamoip.cs.mcgill.ca) that streamlines the computational and visualization processes. Finally, we show recent advances in the prediction of local 3D motifs from sequence data with the BayesPairing software and discuss its impact toward complete 3D structure prediction.

Effects of Noncanonical Base Pairing on RNA Folding: Structural Context and Spatial Arrangements of G·A Pairs

Noncanonical base pairs play important roles in assembling the three-dimensional structures critical to the diverse functions of RNA. These associations contribute to the looped segments that intersperse the canonical double-helical elements within folded, globular RNA molecules. They stitch together various structural elements, serve as recognition elements for other molecules, and act as sites of intrinsic stiffness or deformability. This work takes advantage of new software (DSSR) designed to streamline the analysis and annotation of RNA three-dimensional structures. The multiscale structural information gathered for individual molecules, combined with the growing number of unique, well-resolved RNA structures, makes it possible to examine the collective features deeply and to uncover previously unrecognized patterns of chain organization. Here we focus on a subset of noncanonical base pairs involving guanine and adenine and the links between their modes of association, secondary structural context, and contributions to tertiary folding. The rigorous descriptions of base-pair geometry that we employ facilitate characterization of recurrent geometric motifs and the structural settings in which these arrangements occur. Moreover, the numerical parameters hint at the natural motions of the interacting bases and the pathways likely to connect different spatial forms. We draw attention to higher-order multiplexes involving two or more G·A pairs and the roles these associations appear to play in bridging different secondary structural units. The collective data reveal pairing propensities in base organization, secondary structural context, and deformability and serve as a starting point for further multiscale investigations and/or simulations of RNA folding.

Enzymatic Transition States and Drug Design

Transition state theory teaches that chemically stable mimics of enzymatic transition states will bind tightly to their cognate enzymes. Kinetic isotope effects combined with computational quantum chemistry provides enzymatic transition state information with sufficient fidelity to design transition state analogues. Examples are selected from various stages of drug development to demonstrate the application of transition state theory, inhibitor design, physicochemical characterization of transition state analogues, and their progress in drug development.

A ruthenium-platinum metal complex that binds to sarcin ricin loop RNA and lowers mRNA expression

IT127 is a dinuclear transition metal complex that contains a Pt(ii) and a Ru(iii) metal center. We have shown that IT127 is significantly more effective in binding the 29-base sarcin ricin loop (SRL) RNA in comparison to Cisplatin, a hallmark anticancer agent. Binding site analysis shows that IT127 prefers purine bases and the GAGA tetraloop region of SRL RNA. Our results with a dihydrofolate reductase (DHFR) model system reveal that IT127 binding to mRNA reduces translation of DHFR enzyme and that the Ru(iii) and Pt(ii) centers in IT127 appear to work in a synergistic manner.

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

Structural insights into species-specific features of the ribosome from the human pathogen Mycobacterium tuberculosis

Ribosomes from Mycobacterium tuberculosis (Mtb) possess species-specific ribosomal RNA (rRNA) expansion segments and ribosomal proteins (rProtein). Here, we present the near-atomic structures of the Mtb 50S ribosomal subunit and the complete Mtb 70S ribosome, solved by cryo-electron microscopy. Upon joining of the large and small ribosomal subunits, a 100-nt long expansion segment of the Mtb 23S rRNA, named H54a or the ‘handle’, switches interactions from with rRNA helix H68 and rProtein uL2 to with rProtein bS6, forming a new intersubunit bridge ‘B9’. In Mtb 70S, bridge B9 is mostly maintained, leading to correlated motions among the handle, the L1 stalk and the small subunit in the rotated and non-rotated states. Two new protein densities were discovered near the decoding center and the peptidyl transferase center, respectively. These results provide a structural basis for studying translation in Mtb as well as developing new tuberculosis drugs.

NMR experiments for the rapid identification of P=O···H-X type hydrogen bonds in nucleic acids

Hydrogen bonds involving the backbone phosphate groups occur with high frequency in functional RNA molecules. They are often found in well-characterized tertiary structural motifs presenting powerful probes for the rapid identification of these motifs for structure elucidation purposes. We have shown recently that stable hydrogen bonds to the phosphate backbone can in principle be detected by relatively simple NMR-experiments, providing the identity of both the donor hydrogen and the acceptor phosphorous within the same experiment (Duchardt-Ferner et al., Angew Chem Int Ed Engl 50:7927-7930, 2011). However, for imino and hydroxyl hydrogen bond donor groups rapidly exchanging with the solvent as well as amino groups broadened by conformational exchange experimental sensitivity is severely hampered by extensive line broadening. Here, we present improved methods for the rapid identification of hydrogen bonds to phosphate groups in nucleic acids by NMR. The introduction of the SOFAST technique into ¹H,³¹P-correlation experiments as well as a BEST-HNP experiment exploiting ^3hJ_N,P rather than ^2hJ_H,P coupling constants enables the rapid and sensitive identification of these hydrogen bonds in RNA. The experiments are applicable for larger RNAs (up to ~ 100-nt), for donor groups influenced by conformational exchange processes such as amino groups and for hydrogen bonds with rather labile hydrogens such as 2'-OH groups as well as for moderate sample concentrations. Interestingly, the size of the through-hydrogen bond scalar coupling constants depends not only on the type of the donor group but also on the structural context. The largest coupling constants were measured for hydrogen bonds involving the imino groups of protonated cytosine nucleotides as donors.

Expression of distinct maternal and somatic 5.8S, 18S, and 28S rRNA types during zebrafish development

There is mounting evidence that the ribosome is not a static translation machinery, but a cell-specific, adaptive system. Ribosomal variations have mostly been studied at the protein level, even though the essential transcriptional functions are primarily performed by rRNAs. At the RNA level, oocyte-specific 5S rRNAs are long known for Xenopus. Recently, we described for zebrafish a similar system in which the sole maternal-type 5S rRNA present in eggs is replaced completely during embryonic development by a somatic-type. Here, we report the discovery of an analogous system for the 45S rDNA elements: 5.8S, 18S, and 28S. The maternal-type 5.8S, 18S, and 28S rRNA sequences differ substantially from those of the somatic-type, plus the maternal-type rRNAs are also replaced by the somatic-type rRNAs during embryogenesis. We discuss the structural and functional implications of the observed sequence differences with respect to the translational functions of the 5.8S, 18S, and 28S rRNA elements. Finally, in silico evidence suggests that expansion segments (ES) in 18S rRNA, previously implicated in ribosome-mRNA interaction, may have a preference for interacting with specific mRNA genes. Taken together, our findings indicate that two distinct types of ribosomes exist in zebrafish during development, each likely conducting the translation machinery in a unique way.

Nuclear Magnetic Resonance Structure of an 8 × 8 Nucleotide RNA Internal Loop Flanked on Each Side by Three Watson-Crick Pairs and Comparison to Three-Dimensional Predictions

The prediction of RNA three-dimensional structure from sequence alone has been a long-standing goal. High-resolution, experimentally determined structures of simple noncanonical pairings and motifs are critical to the development of prediction programs. Here, we present the nuclear magnetic resonance structure of the (5'CCAGAAACGGAUGGA)₂ duplex, which contains an 8 × 8 nucleotide internal loop flanked by three Watson-Crick pairs on each side. The loop is comprised of a central 5'AC/3'CA nearest neighbor flanked by two 3RRs motifs, a known stable motif consisting of three consecutive sheared GA pairs. Hydrogen bonding patterns between base pairs in the loop, the all-atom root-mean-square deviation for the loop, and the deformation index were used to compare the structure to automated predictions by MC-sym, RNA FARFAR, and RNAComposer.

Synthesis of RNA 5'-Azides from 2'-O-Pivaloyloxymethyl-Protected RNAs and Their Reactivity in Azide-Alkyne Cycloaddition Reactions

Commercially available 2'-O-pivaloyloxymethyl (PivOM) phosphoramidites were employed in an SPS protocol for RNA 5' azides. The utility of the N₃-RNAs in CuAAC and SPAAC was demonstrated by RNA 5' labeling, chemical ligation including fragment joining and cyclization, and bioconjugation. As a result, several new RNA conjugates that may be valuable tools for studies on biological events such as innate immune response (cyclic dinucleotides), post-transcriptional gene regulation (circular RNAs), or mRNA turnover (m⁷G capped RNAs) were obtained.

RNA Structure: Advances and Assessment of 3D Structure Prediction

Biological functions of RNA molecules are dependent upon sustained specific three-dimensional (3D) structures of RNA, with or without the help of proteins. Understanding of RNA structure is frequently based on 2D structures, which describe only the Watson-Crick (WC) base pairs. Here, we hierarchically review the structural elements of RNA and how they contribute to RNA 3D structure. We focus our analysis on the non-WC base pairs and on RNA modules. Several computer programs have now been designed to predict RNA modules. We describe the RNA-Puzzles initiative, which is a community-wide, blind assessment of RNA 3D structure prediction programs to determine the capabilities and bottlenecks of current predictions. The assessment metrics used in RNA-Puzzles are briefly described. The detection of RNA 3D modules from sequence data and their automatic implementation belong to the current challenges in RNA 3D structure prediction.

RNA damage in biological conflicts and the diversity of responding RNA repair systems

RNA is targeted in biological conflicts by enzymatic toxins or effectors. A vast diversity of systems which repair or ‘heal’ this damage has only recently become apparent. Here, we summarize the known effectors, their modes of action, and RNA targets before surveying the diverse systems which counter this damage from a comparative genomics viewpoint. RNA-repair systems show a modular organization with extensive shuffling and displacement of the constituent domains; however, a general ‘syntax’ is strongly maintained whereby systems typically contain: a RNA ligase (either ATP-grasp or RtcB superfamilies), nucleotidyltransferases, enzymes modifying RNA-termini for ligation (phosphatases and kinases) or protection (methylases), and scaffold or cofactor proteins. We highlight poorly-understood or previously-uncharacterized repair systems and components, e.g. potential scaffolding cofactors (Rot/TROVE and SPFH/Band-7 modules) with their respective cognate non-coding RNAs (YRNAs and a novel tRNA-like molecule) and a novel nucleotidyltransferase associating with diverse ligases. These systems have been extensively disseminated by lateral transfer between distant prokaryotic and microbial eukaryotic lineages consistent with intense inter-organismal conflict. Components have also often been ‘institutionalized’ for non-conflict roles, e.g. in RNA-splicing and in RNAi systems (e.g. in kinetoplastids) which combine a distinct family of RNA-acting prim-pol domains with DICER-like proteins.

Modeling and predicting RNA three-dimensional structures

Modeling the three-dimensional structure of RNAs is a milestone toward better understanding and prediction of nucleic acids molecular functions. Physics-based approaches and molecular dynamics simulations are not tractable on large molecules with all-atom models. To address this issue, coarse-grained models of RNA three-dimensional structures have been developed. In this chapter, we describe a graphical modeling based on the Leontis-Westhof extended base-pair classification. This representation of RNA structures enables us to identify highly conserved structural motifs with complex nucleotide interactions in structure databases. Further, we show how to take advantage of this knowledge to quickly and simply predict three-dimensional structures of large RNA molecules.

An introduction to recurrent nucleotide interactions in RNA

RNA secondary structure diagrams familiar to molecular biologists summarize at a glance the folding of RNA chains to form Watson–Crick paired double helices. However, they can be misleading: First of all, they imply that the nucleotides in loops and linker segments, which can amount to 35% to 50% of a structured RNA, do not significantly interact with other nucleotides. Secondly, they give the impression that RNA molecules are loosely organized in three-dimensional (3D) space. In fact, structured RNAs are compactly folded as a result of numerous long-range, sequence-specific interactions, many of which involve loop or linker nucleotides. Here, we provide an introduction for students and researchers of RNA on the types, prevalence, and sequence variations of inter-nucleotide interactions that structure and stabilize RNA 3D motifs and architectures, using Escherichia coli (E. coli) 16S ribosomal RNA as a concrete example. The picture that emerges is that almost all nucleotides in structured RNA molecules, including those in nominally single-stranded loop or linker regions, form specific interactions that stabilize functional structures or mediate interactions with other molecules. The small number of noninteracting, ‘looped-out’ nucleotides make it possible for the RNA chain to form sharp turns. Base-pairing is the most specific interaction in RNA as it involves edge-to-edge hydrogen bonding (H-bonding) of the bases. Non-Watson–Crick base pairs are a significant fraction (30% or more) of base pairs in structured RNAs.

STAR3D: a stack-based RNA 3D structural alignment tool

The various roles of versatile non-coding RNAs typically require the attainment of complex high-order structures. Therefore, comparing the 3D structures of RNA molecules can yield in-depth understanding of their functional conservation and evolutionary history. Recently, many powerful tools have been developed to align RNA 3D structures. Although some methods rely on both backbone conformations and base pairing interactions, none of them consider the entire hierarchical formation of the RNA secondary structure. One of the major issues is that directly applying the algorithms of matching 2D structures to the 3D coordinates is particularly time-consuming. In this article, we propose a novel RNA 3D structural alignment tool, STAR3D, to take into full account the 2D relations between stacks without the complicated comparison of secondary structures. First, the 3D conserved stacks in the inputs are identified and then combined into a tree-like consensus. Afterward, the loop regions are compared one-to-one in accordance with their relative positions in the consensus tree. The experimental results show that the prediction of STAR3D is more accurate for both non-homologous and homologous RNAs than other state-of-the-art tools with shorter running time.

Deuterated nucleotides as chemical probes of RNA structure: a detailed protocol for the enzymatic synthesis of a complete set of nucleotides specifically deuterated at ribose carbons

Abstract We describe here a detailed protocol for the synthesis of ribonucleotides specifically deuterated at each ribose carbon atom. We synthesized 20 specifically deuterated ribonucleotides: ATP, CTP, GTP, and UTP, each of which contained one of five deuterated riboses (either 1′-D, 2″-D, 3′-D, 4′-D, or 5′,5″-D2). Our synthetic approach is inspired by the pioneering work of Tolbert and Williamson, who developed a method for the convenient one-pot enzymatic synthesis of nucleotides (Tolbert, T. J. and Williamson, J. R. (1996) J. Am. Chem. Soc. 118, 7929–7940). Our protocol consists of a comprehensive list of required chemical and enzymatic reagents and equipment, detailed procedures for enzymatic assays and nucleotide synthesis, and chromatographic procedures for purification of deuterated nucleotides. As an example of the utility of specifically deuterated nucleotides, we used them to synthesize specifically deuterated sarcin/ricin loop (SRL) RNA and measured the deuterium kinetic isotope effect on hydroxyl radical cleavage of the SRL.

Chemical probing of RNA with the hydroxyl radical at single-atom resolution

While hydroxyl radical cleavage is widely used to map RNA tertiary structure, lack of mechanistic understanding of strand break formation limits the degree of structural insight that can be obtained from this experiment. Here, we determine how individual ribose hydrogens of sarcin/ricin loop RNA participate in strand cleavage. We find that substituting deuterium for hydrogen at a ribose 5'-carbon produces a kinetic isotope effect on cleavage; the major cleavage product is an RNA strand terminated by a 5'-aldehyde. We conclude that hydroxyl radical abstracts a 5'-hydrogen atom, leading to RNA strand cleavage. We used this approach to obtain structural information for a GUA base triple, a common tertiary structural feature of RNA. Cleavage at U exhibits a large 5' deuterium kinetic isotope effect, a potential signature of a base triple. Others had noted a ribose-phosphate hydrogen bond involving the G 2'-OH and the U phosphate of the GUA triple, and suggested that this hydrogen bond contributes to backbone rigidity. Substituting deoxyguanosine for G, to eliminate this hydrogen bond, results in a substantial decrease in cleavage at G and U of the triple. We conclude that this hydrogen bond is a linchpin of backbone structure around the triple.

Role of S-turn2 in the structure, dynamics, and function of mitochondrial ribosomal A-site. A bioinformatics and molecular dynamics simulation study

The mRNA decoding site (A-site) in the small ribosomal subunit controls fidelity of the translation process. Here, using molecular dynamics simulations and bioinformatic analyses, we investigated the structural dynamics of the human mitochondrial A-site (native and A1490G mutant) and compared it with the dynamics of the bacterial A-site. We detected and characterized a specific RNA backbone configuration, S-turn2, which occurs in the human mitochondrial but not in the bacterial A-site. Mitochondrial and bacterial A-sites show different propensities to form S-turn2 that may be caused by different base-pairing patterns of the flanking nucleotides. Also, the S-turn2 structural stability observed in the simulations supports higher accuracy and lower speed of mRNA decoding in mitochondria in comparison with bacteria. In the mitochondrial A-site, we observed collective movement of stacked nucleotides A1408·C1409·C1410, which may explain the known differences in aminoglycoside antibiotic binding affinities toward the studied A-site variants.

NMR localization of divalent cations at the active site of the Neurospora VS ribozyme provides insights into RNA-metal-ion interactions

Metal cations represent key elements of RNA structure and function. In the Neurospora VS ribozyme, metal cations play diverse roles; they are important for substrate recognition, formation of the active site, and shifting the pKa's of two key nucleobases that contribute to the general acid-base mechanism. Recently, we determined the NMR structure of the A730 loop of the VS ribozyme active site (SLVI) that contributes the general acid (A756) in the enzymatic mechanism of the cleavage reaction. Our studies showed that magnesium (Mg(2+)) ions are essential to stabilize the formation of the S-turn motif within the A730 loop that exposes the A756 nucleobase for catalysis. In this article, we extend these NMR investigations by precisely mapping the Mg(2+)-ion binding sites using manganese-induced paramagnetic relaxation enhancement and cadmium-induced chemical-shift perturbation of phosphorothioate RNAs. These experiments identify five Mg(2+)-ion binding sites within SLVI. Four Mg(2+) ions in SLVI are associated with known RNA structural motifs, including the G-U wobble pair and the GNRA tetraloop, and our studies reveal novel insights about Mg(2+) ion binding to these RNA motifs. Interestingly, one Mg(2+) ion is specifically associated with the S-turn motif, confirming its structural role in the folding of the A730 loop. This Mg(2+) ion is likely important for formation of the active site and may play an indirect role in catalysis.

Novel class of potential therapeutics that target ricin retrograde translocation

Ricin toxin, an A-B toxin from Ricinus communis, induces cell death through the inhibition of protein synthesis. The toxin binds to the cell surface via its B chain (RTB) followed by its retrograde trafficking through intracellular compartments to the ER where the A chain (RTA) is transported across the membrane and into the cytosol. Ricin A chain is transported across the ER membrane utilizing cellular proteins involved in the disposal of aberrant ER proteins by a process referred to as retrograde translocation. Given the current lack of therapeutics against ricin intoxication, we developed a high-content screen using an enzymatically attenuated RTA chimera engineered with a carboxy-terminal enhanced green fluorescent protein (RTA(E177Q)egfp) to identify compounds that target RTA retrograde translocation. Stabilizing RTA(E177Q)egfp through the inclusion of proteasome inhibitor produced fluorescent peri-nuclear granules. Quantitative analysis of the fluorescent granules provided the basis to discover compounds from a small chemical library (2080 compounds) with known bioactive properties. Strikingly, the screen found compounds that stabilized RTA molecules within the cell and several compounds limited the ability of wild type RTA to suppress protein synthesis. Collectively, a robust high-content screen was developed to discover novel compounds that stabilize intracellular ricin and limit ricin intoxication.

Prediction of RNA 1H and 13C chemical shifts: a structure based approach

The use of NMR-derived chemical shifts in protein structure determination and prediction has received much attention, and, as such, many methods have been developed to predict protein chemical shifts from three-dimensional (3D) coordinates. In contrast, little attention has been paid to predicting chemical shifts from RNA coordinates. Using the random forest machine learning approach, we developed RAMSEY, which is capable of predicting both (1)H and protonated (13)C chemical shifts from RNA coordinates. In this report, we introduce RAMSEY, assess its accuracy, and demonstrate the sensitivity of RAMSEY-predicted chemical shifts to RNA 3D structure.

RNA structural motif recognition based on least-squares distance

RNA structural motifs are recurrent structural elements occurring in RNA molecules. RNA structural motif recognition aims to find RNA substructures that are similar to a query motif, and it is important for RNA structure analysis and RNA function prediction. In view of this, we propose a new method known as RNA Structural Motif Recognition based on Least-Squares distance (LS-RSMR) to effectively recognize RNA structural motifs. A test set consisting of five types of RNA structural motifs occurring in Escherichia coli ribosomal RNA is compiled by us. Experiments are conducted for recognizing these five types of motifs. The experimental results fully reveal the superiority of the proposed LS-RSMR compared with four other state-of-the-art methods.

An innate twist between Crick's wobble and Watson-Crick base pairs

Non-Watson-Crick pairs like the G·U wobble are frequent in RNA duplexes. Their geometric dissimilarity (nonisostericity) with the Watson-Crick base pairs and among themselves imparts structural variations decisive for biological functions. Through a novel circular representation of base pairs, a simple and general metric scheme for quantification of base-pair nonisostericity, in terms of residual twist and radial difference that can also envisage its mechanistic effect, is proposed. The scheme is exemplified by G·U and U·G wobble pairs, and their predicable local effects on helical twist angle are validated by MD simulations. New insights into a possible rationale for contextual occurrence of G·U and other non-WC pairs, as well as the influence of a G·U pair on other non-Watson-Crick pair neighborhood and RNA-protein interactions are obtained from analysis of crystal structure data. A few instances of RNA-protein interactions along the major groove are documented in addition to the well-recognized interaction of the G·U pair along the minor groove. The nonisostericity-mediated influence of wobble pairs for facilitating helical packing through long-range interactions in ribosomal RNAs is also reviewed.

Structure and function of the T-loop structural motif in noncoding RNAs

The T-loop is a frequently occurring five-nucleotide motif found in the structure of noncoding RNAs where it is commonly assumed to play an important role in stabilizing the tertiary RNA structure by facilitating long-range interactions between different regions of the molecule. T-loops were first identified in tRNA(Phe) and a formal consensus sequence for this motif was formulated and later revised based on analyses of the crystal structures of prokaryotic ribosomal RNAs and RNase P and the corresponding primary sequence of their orthologues. In the past decade, several new structures of large RNA molecules have been added to the RCSB Protein Data Bank, including the eukaryotic ribosome, a self-splicing group II intron, numerous synthetases in complex with their cognate transfer RNAs (tRNAs), transfer-messenger RNA (tmRNA) in complex with SmpB, several riboswitches, and a complex of bacterial RNase P bound to its tRNA substrate. In this review, the search for T-loops is extended to these new RNA molecules based on the previously established structure-based criteria. The review highlights and discusses the function and additional roles of T-loops in four broad categories of RNA molecules, namely tRNAs, ribosomal RNAs (rRNAs), P RNAs, and RNA genetic elements. Additionally, the potential application for T-loops as interaction modules is also discussed.

Functional role of the sarcin-ricin loop of the 23S rRNA in the elongation cycle of protein synthesis

The sarcin-ricin loop (SRL) is one of the longest conserved sequences in the 23S ribosomal RNA. The SRL has been accepted as crucial for the activity of the ribosome because it is targeted by cytotoxins such as α-sarcin and ricin that completely abolish translation. Nevertheless, the precise functional role of the SRL in translation is not known. Recent biochemical and structural studies indicate that the SRL is critical for triggering GTP hydrolysis on elongation factor Tu (EF-Tu) and elongation factor G (EF-G). To determine the functional role of the SRL in the elongation stage of protein synthesis, we analyzed mutations in the SRL that are known to abolish protein synthesis and are lethal to cells. Here, we show that the SRL is not critical for GTP hydrolysis on EF-Tu and EF-G. The SRL also is not essential for peptide bond formation. Our results, instead, suggest that the SRL is crucial for anchoring EF-G on the ribosome during mRNA-tRNA translocation.

Identification of a minimal region of the HIV-1 5'-leader required for RNA dimerization, NC binding, and packaging

Assembly of human immunodeficiency virus type 1 (HIV-1) particles is initiated in the cytoplasm by the formation of a ribonucleoprotein complex comprising the dimeric RNA genome and a small number of viral Gag polyproteins. Genomes are recognized by the nucleocapsid (NC) domains of Gag, which interact with packaging elements believed to be located primarily within the 5'-leader (5'-L) of the viral RNA. Recent studies revealed that the native 5'-L exists as an equilibrium of two conformers, one in which dimer-promoting residues and NC binding sites are sequestered and packaging is attenuated, and one in which these sites are exposed and packaging is promoted. To identify the elements within the dimeric 5'-L that are important for packaging, we generated HIV-1 5'-L RNAs containing mutations and deletions designed to eliminate substructures without perturbing the overall structure of the leader and examined effects of the mutations on RNA dimerization, NC binding, and packaging. Our findings identify a 159-residue RNA packaging signal that possesses dimerization and NC binding properties similar to those of the intact 5'-L and contains elements required for efficient RNA packaging.

Neutrons, magnets, and photons: a career in structural biology

The purpose of Reflections articles, it seems, is to give elderly scientists a chance to write about the "good old days," when everyone walked to school in the snow. They enjoy this activity so much that your editor, Martha Fedor, must have known that I would accept her invitation to write such an article, no matter how much I demurred at first. As everyone knows, flattery will get you everywhere. It may comfort the apprehensive reader to learn that there is not going to be much walking to school in the snow in this story. On the contrary, rather than thinking how hard I had it during my scientific career, I find it inconceivable that anyone could have had a smoother ride. At the time I began my career, science was an expanding enterprise in the United States that welcomed the young. Only in such an opportunity-rich environment would someone like me have stood a chance. The contrast between that world and the dog-eat-dog world young scientists confront today is stark.

Clustering RNA structural motifs in ribosomal RNAs using secondary structural alignment

RNA structural motifs are the building blocks of the complex RNA architecture. Identification of non-coding RNA structural motifs is a critical step towards understanding of their structures and functionalities. In this article, we present a clustering approach for de novo RNA structural motif identification. We applied our approach on a data set containing 5S, 16S and 23S rRNAs and rediscovered many known motifs including GNRA tetraloop, kink-turn, C-loop, sarcin–ricin, reverse kink-turn, hook-turn, E-loop and tandem-sheared motifs, with higher accuracy than the state-of-the-art clustering method. We also identified a number of potential novel instances of GNRA tetraloop, kink-turn, sarcin–ricin and tandem-sheared motifs. More importantly, several novel structural motif families have been revealed by our clustering analysis. We identified a highly asymmetric bulge loop motif that resembles the rope sling. We also found an internal loop motif that can significantly increase the twist of the helix. Finally, we discovered a subfamily of hexaloop motif, which has significantly different geometry comparing to the currently known hexaloop motif. Our discoveries presented in this article have largely increased current knowledge of RNA structural motifs.

Adapting yeast as model to study ricin toxin a uptake and trafficking

The plant A/B toxin ricin represents a heterodimeric glycoprotein belonging to the family of ribosome inactivating proteins, RIPs. Its toxicity towards eukaryotic cells results from the depurination of 28S rRNA due to the N-glycosidic activity of ricin toxin A chain, RTA. Since the extention of RTA by a mammalian-specific endoplasmic reticulum (ER) retention signal (KDEL) significantly increases RTA in vivo toxicity against mammalian cells, we here analyzed the phenotypic effect of RTA carrying the yeast-specific ER retention motif HDEL. Interestingly, such a toxin (RTA^HDEL) showed a similar cytotoxic effect on yeast as a corresponding RTA^KDEL variant on HeLa cells. Furthermore, we established a powerful yeast bioassay for RTA in vivo uptake and trafficking which is based on the measurement of dissolved oxygen in toxin-treated spheroplast cultures of S. cerevisiae. We show that yeast spheroplasts are highly sensitive against external applied RTA and further demonstrate that its toxicity is greatly enhanced by replacing the C-terminal KDEL motif by HDEL. Based on the RTA resistant phenotype seen in yeast knock-out mutants defective in early steps of endocytosis (∆end3) and/or in RTA depurination activity on 28S rRNA (∆rpl12B) we feel that the yeast-based bioassay described in this study is a powerful tool to dissect intracellular A/B toxin transport from the plasma membrane through the endosomal compartment to the ER.

Small-molecule inhibitor leads of ribosome-inactivating proteins developed using the doorstop approach

Ribosome-inactivating proteins (RIPs) are toxic because they bind to 28S rRNA and depurinate a specific adenine residue from the α-sarcin/ricin loop (SRL), thereby inhibiting protein synthesis. Shiga-like toxins (Stx1 and Stx2), produced by Escherichia coli, are RIPs that cause outbreaks of foodborne diseases with significant morbidity and mortality. Ricin, produced by the castor bean plant, is another RIP lethal to mammals. Currently, no US Food and Drug Administration-approved vaccines nor therapeutics exist to protect against ricin, Shiga-like toxins, or other RIPs. Development of effective small-molecule RIP inhibitors as therapeutics is challenging because strong electrostatic interactions at the RIP•SRL interface make drug-like molecules ineffective in competing with the rRNA for binding to RIPs. Herein, we report small molecules that show up to 20% cell protection against ricin or Stx2 at a drug concentration of 300 nM. These molecules were discovered using the doorstop approach, a new approach to protein•polynucleotide inhibitors that identifies small molecules as doorstops to prevent an active-site residue of an RIP (e.g., Tyr80 of ricin or Tyr77 of Stx2) from adopting an active conformation thereby blocking the function of the protein rather than contenders in the competition for binding to the RIP. This work offers promising leads for developing RIP therapeutics. The results suggest that the doorstop approach might also be applicable in the development of other protein•polynucleotide inhibitors as antiviral agents such as inhibitors of the Z-DNA binding proteins in poxviruses. This work also calls for careful chemical and biological characterization of drug leads obtained from chemical screens to avoid the identification of irrelevant chemical structures and to avoid the interference caused by direct interactions between the chemicals being screened and the luciferase reporter used in screening assays.

The ribotoxin restrictocin recognizes its RNA substrate by selective engagement of active site residues

Restrictocin and related fungal endoribonucleases from the α-sarcin family site-specifically cleave the sarcin/ricin loop (SRL) on the ribosome to inhibit translation and ultimately trigger cell death. Previous studies showed that the SRL folds into a bulged-G motif and tetraloop, with restrictocin achieving a specificity of ∼1000-fold by recognizing both motifs only after the initial binding step. Here, we identify contacts within the protein−RNA interface and determine the extent to which each one contributes to enzyme specificity by examining the effect of protein mutations on the cleavage of the SRL substrate compared to a variety of other RNA substrates. As with other biomolecular interfaces, only a subset of contacts contributes to specificity. One contact of this subset is critical, with the H49A mutation resulting in quantitative loss of specificity. Maximum catalytic activity occurs when both motifs of the SRL are present, with the major contribution involving the bulged-G motif recognized by three lysine residues located adjacent to the active site: K110, K111, and K113. Our findings support a kinetic proofreading mechanism in which the active site residues H49 and, to a lesser extent, Y47 make greater catalytic contributions to SRL cleavage than to suboptimal substrates. This systematic and quantitative analysis begins to elucidate the principles governing RNA recognition by a site-specific endonuclease and may thus serve as a mechanistic model for investigating other RNA modifying enzymes.

NMR structure of the A730 loop of the Neurospora VS ribozyme: insights into the formation of the active site

The Neurospora VS ribozyme is a small nucleolytic ribozyme with unique primary, secondary and global tertiary structures, which displays mechanistic similarities to the hairpin ribozyme. Here, we determined the high-resolution NMR structure of a stem–loop VI fragment containing the A730 internal loop, which forms part of the active site. In the presence of magnesium ions, the A730 loop adopts a structure that is consistent with existing biochemical data and most likely reflects its conformation in the VS ribozyme prior to docking with the cleavage site internal loop. Interestingly, the A730 loop adopts an S-turn motif that is also present in loop B within the hairpin ribozyme active site. The S-turn appears necessary to expose the Watson–Crick edge of a catalytically important residue (A756) so that it can fulfill its role in catalysis. The A730 loop and the cleavage site loop of the VS ribozyme display structural similarities to internal loops found in the active site of the hairpin ribozyme. These similarities provided a rationale to build a model of the VS ribozyme active site based on the crystal structure of the hairpin ribozyme.

Structure of the K-turn U4 RNA: a combined NMR and SANS study

K-turn motifs are universal RNA structural elements providing a binding platform for proteins in several cellular contexts. Their characteristic is a sharp kink in the phosphate backbone that puts the two helical stems of the protein-bound RNA at an angle of 60°. However, to date no high-resolution structure of a naked K-turn motif is available. Here, we present the first structural investigation at atomic resolution of an unbound K-turn RNA (the spliceosomal U4-Kt RNA) by a combination of NMR and small-angle neutron scattering data. With this study, we wish to address the question whether the K-turn structural motif assumes the sharply kinked conformation in the absence of protein binders and divalent cations. Previous studies have addressed this question by fluorescence resonance energy transfer, biochemical assays and molecular dynamics simulations, suggesting that the K-turn RNAs exist in equilibrium between a kinked conformation, which is competent for protein binding, and a more extended conformation, with the population distribution depending on the concentration of divalent cations. Our data shows that the U4-Kt RNA predominantly assumes the more extended conformation in the absence of proteins and divalent cations. The internal loop region is well structured but adopts a different conformation from the one observed in complex with proteins. Our data suggests that the K-turn consensus sequence does not per se code for the kinked conformation; instead the sharp backbone kink requires to be stabilized by protein binders.

RNAMotifScan: automatic identification of RNA structural motifs using secondary structural alignment

Recent studies have shown that RNA structural motifs play essential roles in RNA folding and interaction with other molecules. Computational identification and analysis of RNA structural motifs remains a challenging task. Existing motif identification methods based on 3D structure may not properly compare motifs with high structural variations. Other structural motif identification methods consider only nested canonical base-pairing structures and cannot be used to identify complex RNA structural motifs that often consist of various non-canonical base pairs due to uncommon hydrogen bond interactions. In this article, we present a novel RNA structural alignment method for RNA structural motif identification, RNAMotifScan, which takes into consideration the isosteric (both canonical and non-canonical) base pairs and multi-pairings in RNA structural motifs. The utility and accuracy of RNAMotifScan is demonstrated by searching for kink-turn, C-loop, sarcin-ricin, reverse kink-turn and E-loop motifs against a 23S rRNA (PDBid: 1S72), which is well characterized for the occurrences of these motifs. Finally, we search these motifs against the RNA structures in the entire Protein Data Bank and the abundances of them are estimated. RNAMotifScan is freely available at our supplementary website (http://genome.ucf.edu/RNAMotifScan).

Syntheses, spectroscopic study and X-ray crystallography of some new phosphoramidates and lanthanide(III) complexes of N-(4-nitrobenzoyl)-N′,N′′-bis(morpholino)phosphoric triamide

New phosphoramidates with the formula RC(O)N(H)P(O)X 2, R = 2-NO2—C6H4, 3-NO2—C6H4 and 4-NO2—C6H4, X = N(CH2CH3) (1)–(3), NC4H8 (4)–(6), and NC4H8O (7)–(9) were synthesized and characterized by 1H, 13C, 31P NMR and IR spectroscopy, and elemental analysis. The reaction of (9) with hydrated lanthanide(III) nitrate leads to ten- or nine-coordinated complexes, (10)–(13). The crystal structure has been determined for (3), (5), (9), (10) and (13). In contrast to all of the previously reported similar phosphoramidate compounds, the —C(O)—N(H)—P(O) skeleton in the free ligand (9) shows a cisoid conformation, with the C=O and P=O double bonds adopting a nearly syn conformation. Quantum chemical calculations were applied for clarifying this exceptional conformational behavior. The monodentate neutral ligand (9) is coordinated to the metal ions via the phosphoryl O atom, adopting the usual anti conformation between the C=O and P=O groups.