NMR cross-correlated relaxation rates reveal ion coordination sites in DNA

Structural, Dynamical, and Entropic Differences between SARS-CoV and SARS-CoV-2 s2m Elements Using Molecular Dynamics Simulations

The functional role of the highly conserved stem-loop II motif (s2m) in SARS-CoV and SARS-CoV-2 in the viral lifecycle remains enigmatic and an intense area of research. Structure and dynamics of the s2m are key to establishing a structure-function connection, yet a full set of atomistic resolution coordinates is not available for SARS-CoV-2 s2m. Our work constructs three-dimensional coordinates consistent with NMR solution phase data for SARS-CoV-2 s2m and provides a comparative analysis with its counterpart SARS-CoV s2m. We employed initial coordinates based on PDB ID 1XJR for SARS-CoV s2m and two models for SARS-CoV-2 s2m: one based on 1XJR in which we introduced the mutations present in SARS-CoV-2 s2m and the second based on the available SARS-CoV-2 NMR NOE data supplemented with knowledge-based methods. For each of the three systems, 3.5 μs molecular dynamics simulations were used to sample the structure and dynamics, and principal component analysis (PCA) reduced the ensembles to hierarchal conformational substates for detailed analysis. Dilute solution simulations of SARS-CoV s2m demonstrate that the GNRA-like terminal pentaloop is rigidly defined by base stacking uniquely positioned for possible kissing dimer formation. However, the SARS-CoV-2 s2m simulation did not retain the reported crystallographic SARS-CoV motifs and the terminal loop expands to a highly dynamic "nonaloop." Increased flexibility and structural disorganization are observed for the larger terminal loop, where an entropic penalty is computed to explain the experimentally observed reduction in kissing complex formation. Overall, both SARS-CoV and SARS-CoV-2 s2m elements have a similarly pronounced L-shape due to different motif interactions. Our study establishes the atomistic three-dimensional structure and uncovers dynamic differences that arise from s2m sequence changes, which sets the stage for the interrogation of different mechanistic pathways of suspected biological function.

Towards Profiling of the G-Quadruplex Targeting Drugs in the Living Human Cells Using NMR Spectroscopy

Recently, the 1H-detected in-cell NMR spectroscopy has emerged as a unique tool allowing the characterization of interactions between nucleic acid-based targets and drug-like molecules in living human cells. Here, we assess the application potential of 1H and 19F-detected in-cell NMR spectroscopy to profile drugs/ligands targeting DNA G-quadruplexes, arguably the most studied class of anti-cancer drugs targeting nucleic acids. We show that the extension of the original in-cell NMR approach is not straightforward. The severe signal broadening and overlap of 1H in-cell NMR spectra of polymorphic G-quadruplexes and their complexes complicate their quantitative interpretation. Nevertheless, the 1H in-cell NMR can be used to identify drugs that, despite strong interaction in vitro, lose their ability to bind G-quadruplexes in the native environment. The in-cell NMR approach is adjusted to a recently developed 3,5-bis(trifluoromethyl)phenyl probe to monitor the intracellular interaction with ligands using 19F-detected in-cell NMR. The probe allows dissecting polymorphic mixture in terms of number and relative populations of individual G-quadruplex species, including ligand-bound and unbound forms in vitro and in cellulo. Despite the probe's discussed limitations, the 19F-detected in-cell NMR appears to be a promising strategy to profile G-quadruplex-ligand interactions in the complex environment of living cells.

HERMES - A Software Tool for the Prediction and Analysis of Magnetic-Field-Induced Residual Dipolar Couplings in Nucleic Acids

Field-Induced Residual Dipolar Couplings (fiRDC) are a valuable source of long-range information on structure of nucleic acids (NA) in solution. A web application (HERMES) was developed for structure-based prediction and analysis of the (fiRDCs) in NA. fiRDC prediction is based on input 3D model structure(s) of NA and a built-in library of nucleobase-specific magnetic susceptibility tensors and reference geometries. HERMES allows three basic applications: (i) the prediction of fiRDCs for a given structural model of NAs, (ii) the validation of experimental or modeled NA structures using experimentally derived fiRDCs, and (iii) assessment of the oligomeric state of the NA fragment and/or the identification of a molecular NA model that is consistent with experimentally derived fiRDC data. Additionally, the program's built-in routine for rigid body modeling allows the evaluation of relative orientation of domains within NA that is in agreement with experimental fiRDCs.

Stability prediction of canonical and non-canonical structures of nucleic acids in various molecular environments and cells

This review provides the biophysicochemical background and recent advances in stability prediction of canonical and non-canonical structures of nucleic acids in various molecular environments and cells.

Advances in the cellular structural biology of nucleic acids

Conventional biophysical and chemical biology approaches for delineating relationships between the structure and biological function of nucleic acids (NAs) abstract NAs from their native biological context. However, cumulative experimental observations have revealed that the structure, dynamics and interactions of NAs might be strongly influenced by a broad spectrum of specific and nonspecific physical-chemical environmental factors. This consideration has recently sparked interest in the development of novel tools for structural characterization of NAs in the native cellular context. Here, we review the individual methods currently being employed for structural characterization of NA structure in a native cellular environment with a focus on recent advances and developments in the emerging fields of in-cell NMR and electron paramagnetic resonance spectroscopy and in-cell single-molecule FRET of NAs.

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.

Solution NMR views of dynamical ordering of biomacromolecules

BACKGROUND

To understand the mechanisms related to the 'dynamical ordering' of macromolecules and biological systems, it is crucial to monitor, in detail, molecular interactions and their dynamics across multiple timescales. Solution nuclear magnetic resonance (NMR) spectroscopy is an ideal tool that can investigate biophysical events at the atomic level, in near-physiological buffer solutions, or even inside cells.

SCOPE OF REVIEW

In the past several decades, progress in solution NMR has significantly contributed to the elucidation of three-dimensional structures, the understanding of conformational motions, and the underlying thermodynamic and kinetic properties of biomacromolecules. This review discusses recent methodological development of NMR, their applications and some of the remaining challenges.

MAJOR CONCLUSIONS

Although a major drawback of NMR is its difficulty in studying the dynamical ordering of larger biomolecular systems, current technologies have achieved considerable success in the structural analysis of substantially large proteins and biomolecular complexes over 1MDa and have characterised a wide range of timescales across which biomolecular motion exists. While NMR is well suited to obtain local structure information in detail, it contributes valuable and unique information within hybrid approaches that combine complementary methodologies, including solution scattering and microscopic techniques.

GENERAL SIGNIFICANCE

For living systems, the dynamic assembly and disassembly of macromolecular complexes is of utmost importance for cellular homeostasis and, if dysregulated, implied in human disease. It is thus instructive for the advancement of the study of the dynamical ordering to discuss the potential possibilities of solution NMR spectroscopy and its applications. This article is part of a Special Issue entitled "Biophysical Exploration of Dynamical Ordering of Biomolecular Systems" edited by Dr. Koichi Kato.

Impact of nucleic acid self-alignment in a strong magnetic field on the interpretation of indirect spin-spin interactions

Heteronuclear and homonuclear direct (D) and indirect (J) spin-spin interactions are important sources of structural information about nucleic acids (NAs). The Hamiltonians for the D and J interactions have the same functional form; thus, the experimentally measured apparent spin-spin coupling constant corresponds to a sum of J and D. In biomolecular NMR studies, it is commonly presumed that the dipolar contributions to Js are effectively canceled due to random molecular tumbling. However, in strong magnetic fields, such as those employed for NMR analysis, the tumbling of NA fragments is anisotropic because the inherent magnetic susceptibility of NAs causes an interaction with the external magnetic field. This motional anisotropy is responsible for non-zero D contributions to Js. Here, we calculated the field-induced D contributions to 33 structurally relevant scalar coupling constants as a function of magnetic field strength, temperature and NA fragment size. We identified two classes of Js, namely (1)JCH and (3)JHH couplings, whose quantitative interpretation is notably biased by NA motional anisotropy. For these couplings, the magnetic field-induced dipolar contributions were found to exceed the typical experimental error in J-coupling determinations by a factor of two or more and to produce considerable over- or under-estimations of the J coupling-related torsion angles, especially at magnetic field strengths >12 T and for NA fragments longer than 12 bp. We show that if the non-zero D contributions to J are not properly accounted for, they might cause structural artifacts/bias in NA studies that use solution NMR spectroscopy.

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.

Unique C. elegans telomeric overhang structures reveal the evolutionarily conserved properties of telomeric DNA

There are two basic mechanisms that are associated with the maintenance of the telomere length, which endows cancer cells with unlimited proliferative potential. One mechanism, referred to as alternative lengthening of telomeres (ALT), accounts for approximately 10–15% of all human cancers. Tumours engaged in the ALT pathway are characterised by the presence of the single stranded 5′-C-rich telomeric overhang (C-overhang). This recently identified hallmark of ALT cancers distinguishes them from healthy tissues and renders the C-overhang as a clear target for anticancer therapy. We analysed structures of the 5′-C-rich and 3′-G-rich telomeric overhangs from human and Caenorhabditis elegans, the recently established multicellular in vivo model of ALT tumours. We show that the telomeric DNA from C. elegans and humans forms fundamentally different secondary structures. The unique structural characteristics of C. elegans telomeric DNA that are distinct not only from those of humans but also from those of other multicellular eukaryotes allowed us to identify evolutionarily conserved properties of telomeric DNA. Differences in structural organisation of the telomeric DNA between the C. elegans and human impose limitations on the use of the C. elegans as an ALT tumour model.

Investigation of quadruplex structure under physiological conditions using in-cell NMR

In this chapter we describe the application of in-cell NMR spectroscopy to the investigation of G-quadruplex structures inside living Xenopus laevis oocytes and in X. laevis egg extract. First, in-cell NMR spectroscopy of nucleic acids (NA) is introduced and applications and limitations of the approach are discussed. In the following text the application of in-cell NMR spectroscopy to investigation of G-quadruplexes are reviewed. Special emphasis is given to the discussion of the influence of the intracellular environmental factors such as low molecular weight compounds, molecular crowding, and hydration on structural behavior of G-quadruplexes. Finally, future perspectives of in-cell NMR spectroscopy for quantitative characterization of G-quadruplexes and NA are discussed.

Probe Position-Dependent Counterion Dynamics in DNA: Comparison of Time-Resolved Stokes Shift of Groove-Bound to Base-Stacked Probes in the Presence of Different Monovalent Counterions

Time-resolved fluorescence Stokes shifts (TRFSS) of 4',6-diamidino-2-phenylindole (DAPI) inside the minor groove of DNA are measured in the presence of three different monovalent counterions: sodium (Na(+)), rubidium (Rb(+)), and tetrabutylammonium (TBA(+)). Fluorescence up-conversion and time-correlated single photon counting are combined to obtain the time-resolved emission spectra (TRES) of DAPI in DNA from 100 fs to 10 ns. Time-resolved Stokes shift data suggest that groove-bound DAPI can not sense the counterion dynamics because the ions are displaced by DAPI far from the probe-site. However, when these results are compared to the earlier base-stacked coumarin data, the same ions are found to affect the nanosecond dynamics significantly. This suggests that the ions come close to the probe-site, such that they can affect the dynamics when measured by base-stacked coumarin. These results support previous molecular dynamics (MD) simulation data of groove-bound and base-stacked probes inside DNA.

Magnesium Ion-Water Coordination and Exchange in Biomolecular Simulations

Magnesium ions have an important role in the structure and folding mechanism of ribonucleic acid systems. To properly simulate these biophysical processes, the applied molecular models should reproduce, among other things, the kinetic properties of the ions in water solution. Here, we have studied the kinetics of the binding of magnesium ions with water molecules and nucleic acid systems using molecular dynamics simulation in detail. We have validated the parameters used in biomolecular force fields, such as AMBER and CHARMM, for Mg(2+) ions and also for the biologically relevant ions Na(+), K(+), and Ca(2+) together with three different water models (TIP3P, SPC/E, and TIP5P). The results show that Mg(2+) ions have a slower exchange rate than Na(+), K(+), and Ca(2+) in agreement with the experimental trend, but the simulated value underestimates the experimentally observed Mg(2+)-water exchange rate by several orders of magnitude, irrespective of the force field and water model. A new set of parameters for Mg(2+) was developed to reproduce the experimental kinetic data. This set also leads to better reproduction of structural data than existing models. We have applied the new parameter set to Mg(2+) binding with a monophosphate model system and with the purine riboswitch, add A-riboswitch. In line with the Mg(2+)-water results, the newly developed parameters show a better description of the structure and kinetics of the Mg(2+)-phosphate binding than all other models. The characterization of the ion binding to the riboswitch system shows that the new parameter set does not affect the global structure of the ribonucleic acid system or the number of ions involved in direct or indirect binding. A slight decrease in the number of water-bridged contacts between A-riboswitch and the Mg(2+) ion is observed. The results support the ability of the newly developed parameters to improve the kinetic description of the Mg(2+) and phosphate ions and their applicability in nucleic acid simulation.