Deflating the RNA Mg2+ bubble. Stereochemistry to the rescue!

Proper evaluation of the ionic structure of biomolecular systems through X ray and cryo-EM techniques remains challenging but is essential for advancing our understanding of the underlying structure/activity/solvent relationships. However, numerous studies overestimate the number of Mg2+ in deposited structures due to assignment errors finding their origin in improper consideration of stereochemical rules. Herein, to tackle such issues, we re-evaluate the PDBid 6QNR and 6SJ6 models of the ribosome ionic structure. We establish that stereochemical principles need to be carefully pondered when evaluating ion binding features, even when K+ anomalous signals are available as it is the case for the 6QNR PDB entry. For ribosomes, assignment errors can result in misleading conceptions of their solvent structure. For instance, present stereochemical analysis result in a significant decrease of the number of assigned Mg2+ in 6QNR, suggesting that K+ and not Mg2+ is the prevalent ion in the ribosome 1st solvation shell. We stress that the use of proper stereochemical guidelines in combination or not with other identification techniques, such as those pertaining to the detection of transition metals, of some anions and of K+ anomalous signals, is critical for deflating the current Mg2+ bubble witnessed in many ribosome and other RNA structures. We also stress that for the identification of lighter ions such as Mg2+, Na+, …, for which no anomalous signals can be detected, stereochemistry coupled with high resolution structures (<2.4 Å) remain the best currently available option.

Nucleobase carbonyl groups are poor Mg2+ inner-sphere binders but excellent monovalent ion binders-a critical PDB survey.

Precise knowledge of Mg2+ inner-sphere binding site properties is vital for understanding the structure and function of nucleic acid systems. Unfortunately, the PDB, which represents the main source of Mg2+ binding sites, contains a substantial number of assignment issues that blur our understanding of the functions of these ions. Here, following a previous study devoted to Mg2+ binding to nucleobase nitrogens, we surveyed nucleic acid X-ray structures from the PDB with resolutions ≤2.9 Å to classify the Mg2+ inner-sphere binding patterns to nucleotide carbonyl, ribose hydroxyl, cyclic ether, and phosphodiester oxygen atoms. From this classification, we derived a set of "prior-knowledge" nucleobase Mg2+ binding sites. We report that crystallographic examples of trustworthy nucleobase Mg2+ binding sites are fewer than expected since many of those are associated with misidentified Na+ or K+ We also emphasize that binding of Na+ and K+ to nucleic acids is much more frequent than anticipated. Overall, we provide evidence derived from X-ray structures that nucleobases are poor inner-sphere binders for Mg2+ but good binders for monovalent ions. Based on strict stereochemical criteria, we propose an extended set of guidelines designed to help in the assignment and validation of ions directly contacting nucleobase and ribose atoms. These guidelines should help in the interpretation of X-ray and cryo-EM solvent density maps. When borderline Mg2+ stereochemistry is observed, alternative placement of Na+, K+, or Ca2+ must be considered. We also critically examine the use of lanthanides (Yb3+, Tb3+) as Mg2+ substitutes in crystallography experiments.

Identification of receptors for UNCG and GNRA Z-turns and their occurrence in rRNA.

In contrast to GNRA tetraloop receptors that are common in RNA, receptors for the more thermostable UNCG loops have remained elusive for almost three decades. An analysis of all RNA structures with resolution ≤3.0 Å from the PDB allowed us to identify three previously unnoticed receptors for UNCG and GNRA tetraloops that adopt a common UNCG fold, named 'Z-turn' in agreement with our previously published nomenclature. These receptors recognize the solvent accessible second Z-turn nucleotide in different but specific ways. Two receptors participating in a complex network of tertiary interactions are associated with the rRNA UUCG and GAAA Z-turns capping helices H62 and H35a in rRNA large subunits. Structural comparison of fully assembled ribosomes and comparative sequence analysis of >6500 rRNA sequences helped us recognize that these motifs are almost universally conserved in rRNA, where they may contribute to organize the large subunit around the subdomain-IV four-way junction. The third UCCG receptor was identified in a rRNA/protein construct crystallized at acidic pH. These three non-redundant Z-turn receptors are relevant for our understanding of the assembly of rRNA and other long-non-coding RNAs, as well as for the design of novel folding motifs for synthetic biology.

Short but Weak: The Z-DNA Lone-Pair⋠⋠⋠π Conundrum Challenges Standard Carbon Van der Waals Radii.

Current interest in lone-pair⋅⋅⋅π (lp⋅⋅⋅π) interactions is gaining momentum in biochemistry and (supramolecular) chemistry. However, the physicochemical origin of the exceptionally short (ca. 2.8 Å) oxygen-to-nucleobase plane distances observed in prototypical Z-DNA CpG steps remains unclear. High-level quantum mechanical calculations, including SAPT2+3 interaction energy decompositions, demonstrate that lp⋅⋅⋅π contacts do not result from n→π* orbital overlaps but from weak dispersion and electrostatic interactions combined with stereochemical effects imposed by the locally strained structural context. They also suggest that the carbon van der Waals (vdW) radii, originally derived for sp3 carbons, should not be used for smaller sp2 carbons attached to electron-withdrawing groups. Using a more adapted carbon vdW radius results in these lp⋅⋅⋅π contacts being no longer of the sub-vdW type. These findings challenge the whole lp⋅⋅⋅π concept that refers to elusive orbital interactions that fail to explain short interatomic contact distances.

Revisiting GNRA and UNCG folds: U-turns versus Z-turns in RNA hairpin loops.

When thinking about RNA three-dimensional structures, coming across GNRA and UNCG tetraloops is perceived as a boon since their folds have been extensively described. Nevertheless, analyzing loop conformations within RNA and RNP structures led us to uncover several instances of GNRA and UNCG loops that do not fold as expected. We noticed that when a GNRA does not assume its "natural" fold, it adopts the one we typically associate with a UNCG sequence. The same folding interconversion may occur for loops with UNCG sequences, for instance within tRNA anticodon loops. Hence, we show that some structured tetranucleotide sequences starting with G or U can adopt either of these folds. The underlying structural basis that defines these two fold types is the mutually exclusive stacking of a backbone oxygen on either the first (in GNRA) or the last nucleobase (in UNCG), generating an oxygen-π contact. We thereby propose to refrain from using sequences to distinguish between loop conformations. Instead, we suggest using descriptors such as U-turn (for "GNRA-type" folds) and a newly described Z-turn (for "UNCG-type" folds). Because tetraloops adopt for the largest part only two (inter)convertible turns, we are better able to interpret from a structural perspective loop interchangeability occurring in ribosomes and viral RNA. In this respect, we propose a general view on the inclination for a given sequence to adopt (or not) a specific fold. We also suggest how long-noncoding RNAs may adopt discrete but transient structures, which are therefore hard to predict.

Mg2+ ions: do they bind to nucleobase nitrogens?

Given the many roles proposed for Mg2+ in nucleic acids, it is essential to accurately determine their binding modes. Here, we surveyed the PDB to classify Mg2+ inner-sphere binding patterns to nucleobase imine N1/N3/N7 atoms. Among those, purine N7 atoms are considered to be the best nucleobase binding sites for divalent metals. Further, Mg2+ coordination to N7 has been implied in several ribozyme catalytic mechanisms. We report that Mg2+ assigned near imine nitrogens derive mostly from poor interpretations of electron density patterns and are most often misidentified Na+, K+, NH4+ ions, water molecules or spurious density peaks. Consequently, apart from few documented exceptions, Mg2+ ions do not bind to N7 atoms. Without much of a surprise, Mn2+, Zn2+ and Cd2+, which have a higher affinity for nitrogens, may contact N7 atoms when present in crystallization buffers. In this respect, we describe for the first time a potential Zn2+ ribosomal binding site involving two purine N7 atoms. Further, we provide a set of guidelines to help in the assignment of Mg2+ in crystallographic, cryo-EM, NMR and model building practices and discuss implications of our findings related to ion substitution experiments.

'Z-DNA like' fragments in RNA: a recurring structural motif with implications for folding, RNA/protein recognition and immune response.

Since the work of Alexander Rich, who solved the first Z-DNA crystal structure, we have known that d(CpG) steps can adopt a particular structure that leads to forming left-handed helices. However, it is still largely unrecognized that other sequences can adopt 'left-handed' conformations in DNA and RNA, in double as well as single stranded contexts. These 'Z-like' steps involve the coexistence of several rare structural features: a C2'-endo puckering, a syn nucleotide and a lone pair-π stacking between a ribose O4' atom and a nucleobase. This particular arrangement induces a conformational stress in the RNA backbone, which limits the occurrence of Z-like steps to ≈0.1% of all dinucleotide steps in the PDB. Here, we report over 600 instances of Z-like steps, which are located within r(UNCG) tetraloops but also in small and large RNAs including riboswitches, ribozymes and ribosomes. Given their complexity, Z-like steps are probably associated with slow folding kinetics and once formed could lock a fold through the formation of unique long-range contacts. Proteins involved in immunologic response also specifically recognize/induce these peculiar folds. Thus, characterizing the conformational features of these motifs could be a key to understanding the immune response at a structural level.

Structure- and machine learning-guided engineering demonstrate that a non-canonical disulfide in an anti-PD-1 rabbit antibody does not impede antibody developability.

Rabbits produce robust antibody responses and have unique features in their antibody repertoire that make them an attractive alternative to rodents for in vivo discovery. However, the frequent occurrence of a non-canonical disulfide bond between complementarity-determining region (CDR) H1 (C35a) and CDRH2 (C50) is often seen as a liability for therapeutic antibody development, despite limited reports of its effect on antibody binding, function, and stability. Here, we describe the discovery and humanization of a human-mouse cross-reactive anti-programmed cell death 1 (PD-1) monoclonal rabbit antibody, termed h1340.CC, which possesses this non-canonical disulfide bond. Initial removal of the non-canonical disulfide resulted in a loss of PD-1 affinity and cross-reactivity, which led us to explore protein engineering approaches to recover these. First, guided by the sequence of a related clone and the crystal structure of h1340.CC in complex with PD-1, we generated variant h1340.SA.LV with a potency and cross-reactivity similar to h1340.CC, but only partially recovered affinity. Side-by-side developability assessment of both h1340.CC and h1340.SA.LV indicate that they possess similar, favorable properties. Next, and prompted by recent developments in machine learning (ML)-guided protein engineering, we used an unbiased ML- and structure-guided approach to rapidly and efficiently generate a different variant with recovered affinity. Our case study thus indicates that, while the non-canonical inter-CDR disulfide bond found in rabbit antibodies does not necessarily constitute an obstacle to therapeutic antibody development, combining structure- and ML-guided approaches can provide a fast and efficient way to improve antibody properties and remove potential liabilities.

Pytheas: a software package for the automated analysis of RNA sequences and modifications via tandem mass spectrometry.

Mass spectrometry is an important method for analysis of modified nucleosides ubiquitously present in cellular RNAs, in particular for ribosomal and transfer RNAs that play crucial roles in mRNA translation and decoding. Furthermore, modifications have effect on the lifetimes of nucleic acids in plasma and cells and are consequently incorporated into RNA therapeutics. To provide an analytical tool for sequence characterization of modified RNAs, we developed Pytheas, an open-source software package for automated analysis of tandem MS data for RNA. The main features of Pytheas are flexible handling of isotope labeling and RNA modifications, with false discovery rate statistical validation based on sequence decoys. We demonstrate bottom-up mass spectrometry characterization of diverse RNA sequences, with broad applications in the biology of stable RNAs, and quality control of RNA therapeutics and mRNA vaccines.

Anions in Nucleic Acid Crystallography.

Nucleic acid crystallization buffers contain a large variety of chemicals fitting specific needs. Among them, anions are often solely considered for pH-regulating purposes and as cationic co-salts while their ability to directly bind to nucleic acid structures is rarely taken into account. Here we review current knowledge related to the use of anions in crystallization buffers along with data on their biological prevalence. Chloride ions are frequently identified in crystal structures but display low cytosolic concentrations. Hence, they are thought to be distant from nucleic acid structures in the cell. Sulfate ions are also frequently identified in crystal structures but their localization in the cell remains elusive. Nevertheless, the characterization of the binding properties of these ions is essential for better interpreting the solvent structure in crystals and consequently, avoiding mislabeling of electron densities. Furthermore, understanding the binding properties of these anions should help to get clues related to their potential effects in crowded cellular environments.

Sodium and Potassium Interactions with Nucleic Acids.

Metal ions are essential cofactors for the structure and functions of nucleic acids. Yet, the early discovery in the 70s of the crucial role of Mg(2+) in stabilizing tRNA structures has occulted for a long time the importance of monovalent cations. Renewed interest in these ions was brought in the late 90s by the discovery of specific potassium metal ions in the core of a group I intron. Their importance in nucleic acid folding and catalytic activity is now well established. However, detection of K(+) and Na(+) ions is notoriously problematic and the question about their specificity is recurrent. Here we review the different methods that can be used to detect K(+) and Na(+) ions in nucleic acid structures such as X-ray crystallography, nuclear magnetic resonance or molecular dynamics simulations. We also discuss specific versus non-specific binding to different structures through various examples.

A comprehensive classification and nomenclature of carboxyl-carboxyl(ate) supramolecular motifs and related catemers: implications for biomolecular systems.

Carboxyl and carboxylate groups form important supramolecular motifs (synthons). Besides carboxyl cyclic dimers, carboxyl and carboxylate groups can associate through a single hydrogen bond. Carboxylic groups can further form polymeric-like catemer chains within crystals. To date, no exhaustive classification of these motifs has been established. In this work, 17 association types were identified (13 carboxyl-carboxyl and 4 carboxyl-carboxylate motifs) by taking into account the syn and anti carboxyl conformers, as well as the syn and anti lone pairs of the O atoms. From these data, a simple rule was derived stating that only eight distinct catemer motifs involving repetitive combinations of syn and anti carboxyl groups can be formed. Examples extracted from the Cambridge Structural Database (CSD) for all identified dimers and catemers are presented, as well as statistical data related to their occurrence and conformational preferences. The inter-carboxyl(ate) and carboxyl(ate)-water hydrogen-bond properties are described, stressing the occurrence of very short (strong) hydrogen bonds. The precise characterization and classification of these supramolecular motifs should be of interest in crystal engineering, pharmaceutical and also biomolecular sciences, where similar motifs occur in the form of pairs of Asp/Glu amino acids or motifs involving ligands bearing carboxyl(ate) groups. Hence, we present data emphasizing how the analysis of hydrogen-containing small molecules of high resolution can help understand structural aspects of larger and more complex biomolecular systems of lower resolution.