Short but Weak: The Z‐DNA Lone‐Pair⋅⋅⋅π Conundrum Challenges Standard Carbon Van der Waals Radii

Light-Up Fluorescence and Circularly Polarized Luminescence in Achiral Interlocked Framework via Adaptive Lone Pair-π Interaction Confinement

Interactions between lone pairs and aromatic π systems are significant across biology and self-assembled materials. Herein, employing an achiral confinement metal-organic framework (MOF) encapsulates guest molecules, it is successfully realized that lone pair (lp)-π interaction induces fluorescence "turn-on" and circularly polarized luminescence for the first time. The MOFs synthesized based on naphthalenediimide show nearly non-emissive, which can be light-up by introducing acetone or ester guests containing lone pairs-π interaction. Furthermore, the introduction of a series of lp-rich chiral esters induces supramolecular chirality as well as circularly polarized luminescence in achiral MOFs, while also observing chiral adaptability. This work first demonstrates the luminescence and chiral induction via lone pair electrons-π interactions, presenting a fresh paradigm for the advancement of chiroptical materials.

RNA dynamics from experimental and computational approaches

Conformational dynamics is crucial for the biological function of RNA molecules and for their potential as therapeutic targets. This meeting report outlines key "take-home" messages that emerged from the presentations and discussions during the CECAM workshop "RNA dynamics from experimental and computational approaches" in Paris, June 26-28, 2023.

Lone-Pair-π Bond Strength Unveiled by a Combined Experimental and Computational Study

A series of naphthalene-diimide (NDI) and perylene-diimide (PDI) connected bis-N-heterocyclic carbene complexes of iridium(III) have been prepared and fully characterized. The analysis of their NMR spectroscopic features, together with their molecular structures show that these species display lone-pair-π interactions between the chloride ligands of the Ir(III) complex and the heterocycles of the NDI/PDI moieties. The detection of this type of interaction in solution is due to the formation of two atropisomers, which are formed as a result of the restricted rotation about the Ir-Ccarbene bond imposed by the (Cl)lp⋅⋅⋅π interaction. Variable-temperature 1H NMR analysis allowed the determination of the strength of this non-covalent interaction, which lies between ΔH=6.6 and 10 kcal/mol. The computational studies performed fully support the experimental findings.

Zα Domain of ADAR1 Binds to an A-Form-like Nucleic Acid Duplex with Low Micromolar Affinity

The left-handed Z-conformation of nucleic acids can be adopted by both DNA and RNA when bound by Zα domains found within a variety of viral and innate immune response proteins. While Z-form adoption is preferred by certain sequences, such as the commonly studied (CpG)n repeats, Zα has been reported to bind to a wide range of sequence contexts. Studying how Zα interacts with B-/A-form helices prior to their conversion to the Z-conformation is challenging as binding coincides with Z-form adoption. Here, we studied the binding of Zα fromHomo sapiens ADAR1 to a locked "A-type" version of the (CpG)3 construct (LNA (CpG)3) where the sugar pucker is locked into the C3'-endo/C2'-exo conformation, which prevents the duplex from adopting the alternating C2'/C3'-endo sugar puckers found in the Z-conformation. Using NMR and other biophysical techniques, we find that ZαADAR1 binds to the LNA (CpG)3 using a similar interface as for Z-form binding, with a dissociation constant (KD) of ∼4 μM. In contrast to Z-DNA/Z-RNA, where two ZαADAR1 bind to every 6 bp stretch, our data suggests that ZαADAR1 binds to multiple LNA molecules, indicating a completely different binding mode. Because ZαADAR1 binds relatively tightly to a non-Z-form model, its binding to B/A-form helices may need to be considered when experiments are carried out which attempt to identify the Z-form targets of Zα domains. The use of LNA constructs may be beneficial in experiments where negative controls for Z-form adoption are needed.

Z-Form Adoption of Nucleic Acid is a Multi-Step Process Which Proceeds through a Melted Intermediate

The left-handed Z-conformation of nucleic acids can be adopted by both DNA and RNA when bound by Zα domains found within a variety of innate immune response proteins. Zα domains stabilize this higher-energy conformation by making specific interactions with the unique geometry of Z-DNA/Z-RNA. However, the mechanism by which a right-handed helix contorts to become left-handed in the presence of proteins, including the intermediate steps involved, is poorly understood. Through a combination of nuclear magnetic resonance (NMR) and other biophysical measurements, we have determined that in the absence of Zα, under low salt conditions at room temperature, d(CpG) and r(CpG) constructs show no observable evidence of transient Z-conformations greater than 0.5% on either the intermediate or slow NMR time scales. At higher temperatures, we observed a transient unfolded intermediate. The ease of melting a nucleic acid duplex correlates with Z-form adoption rates in the presence of Zα. The largest contributing factor to the activation energies of Z-form adoption as calculated by Arrhenius plots is the ease of flipping the sugar pucker, as required for Z-DNA and Z-RNA. Together, these data validate the previously proposed "zipper model" for Z-form adoption in the presence of Zα. Overall, Z-conformations are more likely to be adopted by double-stranded DNA and RNA regions flanked by less stable regions and by RNAs experiencing torsional/mechanical stress.

Z-RNA biology: a central role in the innate immune response?

Z-RNA is a higher-energy, left-handed conformation of RNA, whose function has remained elusive. A growing body of work alludes to regulatory roles for Z-RNA in the immune response. Here, we review how Z-RNA features present in cellular RNAs-especially containing retroelements-could be recognized by a family of winged helix proteins, with an impact on host defense. We also discuss how mutations to specific Z-contacting amino acids disrupt their ability to stabilize Z-RNA, resulting in functional losses. We end by highlighting knowledge gaps in the field, which, if addressed, would significantly advance this active area of research.

Structure and Formation of Z-DNA and Z-RNA

Despite structural differences between the right-handed conformations of A-RNA and B-DNA, both nucleic acids adopt very similar, left-handed Z-conformations. In contrast to their structural similarities and sequence preferences, RNA and DNA exhibit differences in their ability to adopt the Z-conformation regarding their hydration shells, the chemical modifications that promote the Z-conformation, and the structure of junctions connecting them to right-handed segments. In this review, we highlight the structural and chemical properties of both Z-DNA and Z-RNA and delve into the potential factors that contribute to both their similarities and differences. While Z-DNA has been extensively studied, there is a gap of knowledge when it comes to Z-RNA. Where such information is lacking, we try and extend the principles of Z-DNA stability and formation to Z-RNA, considering the inherent differences of the nucleic acids.

Multiscale Modeling of Phosphate···π Contacts in RNA U-Turns Exposes Differences between Quantum-Chemical and AMBER Force Field Descriptions

Phosphate···π, also called anion···π, contacts occur between nucleobases and anionic phosphate oxygens (OP2) in r(GNRA) and r(UNNN) U-turn motifs (N = A,G,C,U; R = A,G). These contacts were investigated using state-of-the-art quantum-chemical methods (QM) to characterize their physicochemical properties and to serve as a reference to evaluate AMBER force field (AFF) performance. We found that phosphate···π interaction energies calculated with the AFF for dimethyl phosphate···nucleobase model systems are less stabilizing in comparison with double-hybrid DFT and that minimum contact distances are larger for all nucleobases. These distance stretches are also observed in large-scale AFF vs QM/MM computations and classical molecular dynamics (MD) simulations on several r(gcGNRAgc) tetraloop hairpins when compared to experimental data extracted from X-ray/cryo-EM structures (res. ≤ 2.5 Å) using the WebFR3D bioinformatic tool. MD simulations further revealed shifted OP2/nucleobase positions. We propose that discrepancies between the QM and AFF result from a combination of missing polarization in the AFF combined with too large AFF Lennard-Jones (LJ) radii of nucleobase carbon atoms in addition to an exaggerated short-range repulsion of the r^-12 LJ repulsive term. We compared these results with earlier data gathered on lone pair···π contacts in CpG Z-steps occurring in r(UNCG) tetraloops. In both instances, charge transfer calculations do not support any significant n → π* donation effects. We also investigated thiophosphate···π contacts that showed reduced stabilizing interaction energies when compared to phosphate···π contacts. Thus, we challenge suggestions that the experimentally observed enhanced thermodynamic stability of phosphorothioated r(GNRA) tetraloops can be explained by larger London dispersion.

Occurrence and stability of anion-π interactions between phosphate and nucleobases in functional RNA molecules

We present a systematic structural and energetic characterization of phosphate(OP)-nucleobase anion…π stacking interactions in RNAs. We observed OP-nucleobase stacking contacts in a variety of structural motifs other than regular helices and spanning broadly diverse sequence distances. Apart from the stacking between a phosphate and a guanine or a uracil two-residue upstream in specific U-turns, such interactions in RNA have been scarcely characterized to date. Our QM calculations showed an energy minimum at a distance between the OP atom and the nucleobase plane centroid slightly below 3 Å for all the nucleobases. By sliding the OP atom over the nucleobase plane we localized the optimal mutual positioning of the stacked moieties, corresponding to an energy minimum below -6 kcal•mol-1, for all the nucleobases, consistently with the projections of the OP atoms over the different π-rings we observed in experimental occurrences. We also found that the strength of the interaction clearly correlates with its electrostatic component, pointing to it as the most relevant contribution. Finally, as OP-uracil and OP-guanine interactions represent together 86% of the instances we detected, we also proved their stability under dynamic conditions in model systems simulated by state-of-the art DFT-MD calculations.

Lone Pair…π Contacts and Structure Signatures of r(UNCG) Tetraloops, Z-Turns, and Z-Steps: A WebFR3D Survey

Z-DNA and Z-RNA have long appeared as oddities to nucleic acid scientists. However, their Z-step constituents are recurrently observed in all types of nucleic acid systems including ribosomes. Z-steps are NpN steps that are isostructural to Z-DNA CpG steps. Among their structural features, Z-steps are characterized by the presence of a lone pair…π contact that involves the stacking of the ribose O4′ atom of the first nucleotide with the 3′-face of the second nucleotide. Recently, it has been documented that the CpG step of the ubiquitous r(UNCG) tetraloops is a Z-step. Accordingly, such r(UNCG) conformations were called Z-turns. It has also been recognized that an r(GAAA) tetraloop in appropriate conditions can shapeshift to an unusual Z-turn conformation embedding an ApA Z-step. In this report, we explore the multiplicity of RNA motifs based on Z-steps by using the WebFR3D tool to which we added functionalities to be able to retrieve motifs containing lone pair…π contacts. Many examples that underscore the diversity and universality of these motifs are provided as well as tutorial guidance on using WebFR3D. In addition, this study provides an extensive survey of crystallographic, cryo-EM, NMR, and molecular dynamics studies on r(UNCG) tetraloops with a critical view on how to conduct database searches and exploit their results.

Toward Convergence in Folding Simulations of RNA Tetraloops: Comparison of Enhanced Sampling Techniques and Effects of Force Field Modifications

Atomistic molecular dynamics simulations represent an established technique for investigation of RNA structural dynamics. Despite continuous development, contemporary RNA simulations still suffer from suboptimal accuracy of empirical potentials (force fields, ffs) and sampling limitations. Development of efficient enhanced sampling techniques is important for two reasons. First, they allow us to overcome the sampling limitations, and second, they can be used to quantify ff imbalances provided they reach a sufficient convergence. Here, we study two RNA tetraloops (TLs), namely the GAGA and UUCG motifs. We perform extensive folding simulations and calculate folding free energies (ΔG_fold^°) with the aim to compare different enhanced sampling techniques and to test several modifications of the nonbonded terms extending the AMBER OL3 RNA ff. We demonstrate that replica-exchange solute tempering (REST2) simulations with 12-16 replicas do not show any sign of convergence even when extended to a timescale of 120 μs per replica. However, the combination of REST2 with well-tempered metadynamics (ST-MetaD) achieves good convergence on a timescale of 5-10 μs per replica, improving the sampling efficiency by at least 2 orders of magnitude. Effects of ff modifications on ΔG_fold^° energies were initially explored by the reweighting approach and then validated by new simulations. We tested several manually prepared variants of the gHBfix potential which improve stability of the native state of both TLs by ∼2 kcal/mol. This is sufficient to conveniently stabilize the folded GAGA TL while the UUCG TL still remains under-stabilized. Appropriate adjustment of van der Waals parameters for C-H···O5' base-phosphate interaction may further stabilize the native states of both TLs by ∼0.6 kcal/mol.

Short-Range Imbalances in the AMBER Lennard-Jones Potential for (Deoxy)Ribose···Nucleobase Lone-Pair···π Contacts in Nucleic Acids

The lone-pair···π (lp···π) (deoxy)ribose···nucleobase stacking is a recurring interaction in Z-DNA and RNAs that is characterized by sub-van der Waals lp···π contacts (<3.0 Å). It is a part of the structural signature of CpG Z-step motifs in Z-DNA and r(UNCG) tetraloops that are known to behave poorly in molecular dynamics (MD) simulations. Although the exact origin of the MD simulation issues remains unclear, a significant part of the problem might be due to an imbalanced description of nonbonded interactions, including the characteristic lp···π stacking. To gain insights into the links between lp···π stacking and MD, we present an in-depth comparison between accurate large-basis-set double-hybrid Kohn-Sham density functional theory calculations DSD-BLYP-D3/ma-def2-QZVPP (DHDF-D3) and data obtained with the nonbonded potential of the AMBER force field (AFF) for NpN Z-steps (N = G, A, C, and U). Among other differences, we found that the AFF overestimates the DHDF-D3 lp···π distances by ∼0.1-0.2 Å, while the deviation between the DHDF-D3 and AFF descriptions sharply increases in the short-range region of the interaction. Based on atom-in-molecule polarizabilities and symmetry-adapted perturbation theory analysis, we inferred that the DHDF-D3 versus AFF differences partly originate in identical nucleobase carbon atom Lennard-Jones (LJ) parameters despite the presence/absence of connected electron-withdrawing groups that lead to different effective volumes or vdW radii. Thus, to precisely model the very short CpG lp···π contact distances, we recommend revision of the nucleobase atom LJ parameters. Additionally, we suggest that the large discrepancy between DHDF-D3 and AFF short-range repulsive part of the interaction energy potential may significantly contribute to the poor performances of MD simulations of nucleic acid systems containing Z-steps. Understanding where, and if possible why, the point-charge-type effective potentials reach their limits is vital for developing next-generation FFs and for addressing specific issues in contemporary MD simulations.

Solution NMR backbone assignments of the N-terminal Zα-linker-Zβ segment from Homo sapiens ADAR1p150

Adenosine-to-inosine (A-to-I) editing of a subset of RNAs in a eukaryotic cell is required in order to avoid triggering the innate immune system. Editing is carried out by ADAR1, which exists as short (p110) and long (p150) isoforms. ADAR1p150 is mostly cytoplasmic, possesses a Z-RNA binding domain (Zα), and is only expressed during the innate immune response. A structurally homologous domain to Zα, the Zβ domain, is separated by a long linker from Zα on the N-terminus of ADAR1 but its function remains unknown. Zβ does not bind to RNA in isolation, yet the binding kinetics of the segment encompassing Zα, Zβ and the 95-residue linker between the two domains (Zα-Zβ) are markedly different compared to Zα alone. Here we present the solution NMR backbone assignment of Zα-Zβ from H. Sapiens ADAR1. The predicted secondary structure of Zα-Zβ based on chemical shifts is in agreement with previously determined structures of Zα and Zβ in isolation, and indicates that the linker is intrinsically disordered. Comparison of the chemical shifts between the individual Zα and Zβ domains to the full Zα-Zβ construct suggests that Zβ may interact with the linker, the function of which is currently unknown.

Switching Aromatic Character by Complexation: π to π* Change Seen in Molecular Rotation Spectra

A dominating F···π*_aromatic interaction is found to govern the benzaldehyde···tetrafluoromethane complex, as revealed by this rotational spectroscopic investigation. Secondary F···π*_-C=O- and C σ*_CF4···π_aromatic interactions also contribute to the stability of the observed isomer. Narrow splittings have been observed in the rotational spectrum originating from a 3-fold internal rotation of CF₄ above the aromatic moiety, and a corresponding V₃ barrier was determined to be 1.572(14) kJ mol^-1. This is the first rotational spectroscopic evidence in the literature implying that the aromatic π* antibonding orbital can be activated not only by electron-withdrawing substituents but also by complexation partners containing atoms with high electronegativity, like CF₄. The results emphasize the partner molecules' role to modulate the π electron structure and show a change in the orbital character (π or π*) when participating in the formation of noncovalent interactions.

Recognition of non-CpG repeats in Alu and ribosomal RNAs by the Z-RNA binding domain of ADAR1 induces A-Z junctions

Adenosine-to-inosine (A-to-I) editing of eukaryotic cellular RNAs is essential for protection against auto-immune disorders. Editing is carried out by ADAR1, whose innate immune response-specific cytoplasmic isoform possesses a Z-DNA binding domain (Zα) of unknown function. Zα also binds to CpG repeats in RNA, which are a hallmark of Z-RNA formation. Unexpectedly, Zα has been predicted — and in some cases even shown — to bind to specific regions within mRNA and rRNA devoid of such repeats. Here, we use NMR, circular dichroism, and other biophysical approaches to demonstrate and characterize the binding of Zα to mRNA and rRNA fragments. Our results reveal a broad range of RNA sequences that bind to Zα and adopt Z-RNA conformations. Binding is accompanied by destabilization of neighboring A-form regions which is similar in character to what has been observed for B-Z-DNA junctions. The binding of Zα to non-CpG sequences is specific, cooperative and occurs with an affinity in the low micromolar range. This work allows us to propose a model for how Zα could influence the RNA binding specificity of ADAR1.

ADAR1 is an interferon-induced enzyme that catalyzes editing of adenine to inosine across the transcriptome as part of the immune response. Here the authors establish how ADAR1 recognizes non-CpG RNA sequences to facilitate the formation of A-Z junctions.

When are two hydrogen bonds better than one? Accurate first-principles models explain the balance of hydrogen bond donors and acceptors found in proteins

Hydrogen bonds (HBs) play an essential role in the structure and catalytic action of enzymes, but a complete understanding of HBs in proteins challenges the resolution of modern structural (i.e., X-ray diffraction) techniques and mandates computationally demanding electronic structure methods from correlated wavefunction theory for predictive accuracy. Numerous amino acid sidechains contain functional groups (e.g., hydroxyls in Ser/Thr or Tyr and amides in Asn/Gln) that can act as either HB acceptors or donors (HBA/HBD) and even form simultaneous, ambifunctional HB interactions. To understand the relative energetic benefit of each interaction, we characterize the potential energy surfaces of representative model systems with accurate coupled cluster theory calculations. To reveal the relationship of these energetics to the balance of these interactions in proteins, we curate a set of 4000 HBs, of which >500 are ambifunctional HBs, in high-resolution protein structures. We show that our model systems accurately predict the favored HB structural properties. Differences are apparent in HBA/HBD preference for aromatic Tyr versus aliphatic Ser/Thr hydroxyls because Tyr forms significantly stronger O–H⋯O HBs than N–H⋯O HBs in contrast to comparable strengths of the two for Ser/Thr. Despite this residue-specific distinction, all models of residue pairs indicate an energetic benefit for simultaneous HBA and HBD interactions in an ambifunctional HB. Although the stabilization is less than the additive maximum due both to geometric constraints and many-body electronic effects, a wide range of ambifunctional HB geometries are more favorable than any single HB interaction.

Correlated wavefunction theory predicts and high-resolution crystal structure analysis confirms the important, stabilizing effect of simultaneous hydrogen bond donor and acceptor interactions in proteins.

UUCG RNA Tetraloop as a Formidable Force-Field Challenge for MD Simulations

Explicit solvent atomistic molecular dynamics (MD) simulations represent an established technique to study structural dynamics of RNA molecules and an important complement for diverse experimental methods. However, performance of molecular mechanical (MM) force fields (ff's) remains far from satisfactory even after decades of development, as apparent from a problematic structural description of some important RNA motifs. Actually, some of the smallest RNA molecules belong to the most challenging systems for MD simulations and, among them, the UUCG tetraloop is saliently difficult. We report a detailed analysis of UUCG MD simulations, depicting the sequence of events leading to the loss of the UUCG native state during MD simulations. The total amount of MD simulation data analyzed in this work is close to 1.3 ms. We identify molecular interactions, backbone conformations, and substates that are involved in the process. Then, we unravel specific ff deficiencies using diverse quantum mechanical/molecular mechanical (QM/MM) and QM calculations. Comparison between the MM and QM methods shows discrepancies in the description of the 5'-flanking phosphate moiety and both signature sugar-base interactions. Our work indicates that poor behavior of the UUCG tetraloop in simulations is a complex issue that cannot be attributed to one dominant and straightforwardly correctable factor. Instead, there is a concerted effect of multiple ff inaccuracies that are coupled and amplifying each other. We attempted to improve the simulation behavior by some carefully tailored interventions, but the results were still far from satisfactory, underlying the difficulties in development of accurate nucleic acid ff's.