Accurate ab initio prediction of NMR chemical shifts of nucleic acids and nucleic acids/protein complexes

NMR chemical shift predictions based on empirical methods are nowadays indispensable tools during resonance assignment and 3D structure calculation of proteins. However, owing to the very limited statistical data basis, such methods are still in their infancy in the field of nucleic acids, especially when non-canonical structures and nucleic acid complexes are considered. Here, we present an ab initio approach for predicting proton chemical shifts of arbitrary nucleic acid structures based on state-of-the-art fragment-based quantum chemical calculations. We tested our prediction method on a diverse set of nucleic acid structures including double-stranded DNA, hairpins, DNA/protein complexes and chemically-modified DNA. Overall, our quantum chemical calculations yield highly/very accurate predictions with mean absolute deviations of 0.3–0.6 ppm and correlation coefficients (r²) usually above 0.9. This will allow for identifying misassignments and validating 3D structures. Furthermore, our calculations reveal that chemical shifts of protons involved in hydrogen bonding are predicted significantly less accurately. This is in part caused by insufficient inclusion of solvation effects. However, it also points toward shortcomings of current force fields used for structure determination of nucleic acids. Our quantum chemical calculations could therefore provide input for force field optimization.     

Protein Structure Validation and Refinement Using Amide Proton Chemical Shifts Derived from Quantum Mechanics

We present the ProCS method for the rapid and accurate prediction of protein backbone amide proton chemical shifts - sensitive probes of the geometry of key hydrogen bonds that determine protein structure. ProCS is parameterized against quantum mechanical (QM) calculations and reproduces high level QM results obtained for a small protein with an RMSD of 0.25 ppm (r = 0.94). ProCS is interfaced with the PHAISTOS protein simulation program and is used to infer statistical protein ensembles that reflect experimentally measured amide proton chemical shift values. Such chemical shift-based structural refinements, starting from high-resolution X-ray structures of Protein G, ubiquitin, and SMN Tudor Domain, result in average chemical shifts, hydrogen bond geometries, and trans-hydrogen bond (^h3 J_NC'') spin-spin coupling constants that are in excellent agreement with experiment. We show that the structural sensitivity of the QM-based amide proton chemical shift predictions is needed to obtain this agreement. The ProCS method thus offers a powerful new tool for refining the structures of hydrogen bonding networks to high accuracy with many potential applications such as protein flexibility in ligand binding.     

The proto-Nucleic Acid Builder: a software tool for constructing nucleic acid analogs

The helical structures of DNA and RNA were originally revealed by experimental data. Likewise, the development of programs for modeling these natural polymers was guided by known structures. These nucleic acid polymers represent only two members of a potentially vast class of polymers with similar structural features, but that differ from DNA and RNA in the backbone or nucleobases. Xeno nucleic acids (XNAs) incorporate alternative backbones that affect the conformational, chemical, and thermodynamic properties of XNAs. Given the vast chemical space of possible XNAs, computational modeling of alternative nucleic acids can accelerate the search for plausible nucleic acid analogs and guide their rational design. Additionally, a tool for the modeling of nucleic acids could help reveal what nucleic acid polymers may have existed before RNA in the early evolution of life. To aid the development of novel XNA polymers and the search for possible pre-RNA candidates, this article presents the proto-Nucleic Acid Builder (https://github.com/GT-NucleicAcids/pnab), an open-source program for modeling nucleic acid analogs with alternative backbones and nucleobases. The torsion-driven conformation search procedure implemented here predicts structures with good accuracy compared to experimental structures, and correctly demonstrates the correlation between the helical structure and the backbone conformation in DNA and RNA.     

The crystal structure of an ‘All Locked’ nucleic acid duplex

‘Locked nucleic acids’ (LNAs) are known to introduce enhanced bio- and thermostability into natural nucleic acids rendering them powerful tools for diagnostic and therapeutic applications. We present the 1.9 Å X-ray structure of an ‘all LNA’ duplex containing exclusively modified β-d-2′-O-4′C-methylene ribofuranose nucleotides. The helix illustrates a new type of nucleic acid geometry that contributes to the understanding of the enhanced thermostability of LNA duplexes. A notable decrease of several local and overall helical parameters like twist, roll and propeller twist influence the structure of the LNA helix and result in a widening of the major groove, a decrease in helical winding and an enlarged helical pitch. A detailed structural comparison to the previously solved RNA crystal structure with the corresponding base pair sequence underlines the differences in conformation. The surrounding water network of the RNA and the LNA helix shows a similar hydration pattern.     

Proton and phosphorus NMR studies of d-CpG(pCpG)n duplexes in solution. Helix-coil transition and complex formation with actinomycin-D.

The Watson–Crick imino and amino exchangeable protons, the nonexchangeable base and sugar protons, and the backbone phosphates for d-CpG(pCpG)_n, n = 1 and 2, have been monitored by high-resolution nmr spectroscopy in aqueous solution over the temperature range 0°–90°C. The temperature dependence of the chemical shifts of the tetramer and hexamer resonances is consistent with the formation of stable duplexes at low temperature in solution. Comparison of the spectral characteristics of the tetranucleotide with those of the hexanucleotide with temperature permits the differentiation and assignment of the cytosine proton resonances on base pairs located at the end of the helix from those in an interior position. There is fraying at the terminal base pairs in the tetranucleotide and hexanucleotide duplexes. The Watson–Crick ring imino protons exchange at a faster rate than the Watson–Crick side-chain amino protons, with exchange occurring by transient opening of the double helix. The structure of the d-CpG(pCpG)_n double helices has been probed by proton relaxation time measurements, sugar proton coupling constants, and the proton chemical shift changes associated with the helix–coil transition. The experimental data support a structural model in solution, which incorporates an anti conformation about the glycosyl bonds, C(3) exo sugar ring pucker, and base overlap geometries similar to the B-DNA helix. Rotational correlation times of 1.7 and 0.9 × 10^?9 sec have been computed for the hexanucleotide and tetranucleotide duplexes in 0.1 M salt, D₂O, pH 6.25 at 27°C. The well-resolved ³¹P resonances for the internucleotide phosphates of the tetramer and hexamer sequences at superconducting fields shift upfield by 0.2–0.5 ppm on helix formation. These shifts reflect a conformational change about the ω,ω′ phosphodiester bonds from gauche-gauche in the duplex structure to a distribution of gauche-trans states in the coil structure. Significant differences are observed in the transition width and midpoint of the chemical shift versus temperature profiles plotted in differentiated form for the various base and sugar proton and internucleotide phosphorous resonances monitoring the d-CpG(pCpG)_n helix–coil transition. The twofold symmetry of the d-CpGpCpG duplex is removed on complex formation with the antibiotic actinomycin-D. Two phosphorous resonances are shifted downfield by ～2.6 ppm and ～1.6 ppm on formation of the 1:2 Act-D:d-CpGpCpG complex in solution. Model studies on binding of the antibiotic to dinucleotides of varying sequence indicate that intercalation of the actinomycin-D occurs at the GpC site in the d-CpGpCpG duplex and that the magnitude of the downfield shifts reflects strain at the O-P-O backbone angles and hydrogen bonding between the phenoxazone and the phosphate oxygens. Actinomycin-D is known to bind to nucleic acids that exhibit a B-DNA conformation; this suggests that the d-CpG(pCpG)_n duplexes exhibit a B-DNA conformation in solution.     

2DCSi: identification of protein secondary structure and redox state using 2D cluster analysis of NMR chemical shifts

Chemical shifts of amino acids in proteins are the most sensitive and easily obtainable NMR parameters that reflect the primary, secondary, and tertiary structures of the protein. In recent years, chemical shifts have been used to identify secondary structure in peptides and proteins, and it has been confirmed that ¹H^α, ¹³C^α, ¹³C^β, and ¹³C′ NMR chemical shifts for all 20 amino acids are sensitive to their secondary structure. Currently, most of the methods are purely based on one-dimensional statistical analyses of various chemical shifts for each residue to identify protein secondary structure. However, it is possible to achieve an increased accuracy from the two-dimensional analyses of these chemical shifts. The 2DCSi approach performs two-dimension cluster analyses of ¹H^α, ¹H^N, ¹³C^α, ¹³C^β, ¹³C′, and ¹⁵N^H chemical shifts to identify protein secondary structure and the redox state of cysteine residue. For the analysis of paired chemical shifts of 6 data sets, each of the 20 amino acids has its own 15 two-dimension cluster scattering diagrams. Accordingly, the probabilities for identifying helix and extended structure were calculated by using our scoring matrix. Compared with existing the chemical shift-based methods, it appears to improve the prediction accuracy of secondary structure identification, particularly in the extended structure. In addition, the probability of the given residue to be helix or extended structure is displayed, allows the users to make decisions by themselves. Electronic Supplementary Material The online version of this article (doi:) contains supplementary material, which is available to authorized users. Grant sponsor: National Science Council of ROC; Grant numbers: NSC-94-2323-B006- 001, NSC-93-2212-E-006.     

Practical use of chemical shift databases for protein solid-state NMR: 2D chemical shift maps and amino-acid assignment with secondary-structure information

We introduce a Python-based program that utilizes the large database of ¹³C and ¹⁵N chemical shifts in the Biological Magnetic Resonance Bank to rapidly predict the amino acid type and secondary structure from correlated chemical shifts. The program, called PACSYlite Unified Query (PLUQ), is designed to help assign peaks obtained from 2D ¹³C–¹³C, ¹⁵N–¹³C, or 3D ¹⁵N–¹³C–¹³C magic-angle-spinning correlation spectra. We show secondary-structure specific 2D ¹³C–¹³C correlation maps of all twenty amino acids, constructed from a chemical shift database of 262,209 residues. The maps reveal interesting conformation-dependent chemical shift distributions and facilitate searching of correlation peaks during amino-acid type assignment. Based on these correlations, PLUQ outputs the most likely amino acid types and the associated secondary structures from inputs of experimental chemical shifts. We test the assignment accuracy using four high-quality protein structures. Based on only the Cα and Cβ chemical shifts, the highest-ranked PLUQ assignments were 40–60 % correct in both the amino-acid type and the secondary structure. For three input chemical shifts (CO–Cα–Cβ or N–Cα–Cβ), the first-ranked assignments were correct for 60 % of the residues, while within the top three predictions, the correct assignments were found for 80 % of the residues. PLUQ and the chemical shift maps are expected to be useful at the first stage of sequential assignment, for combination with automated sequential assignment programs, and for highly disordered proteins for which secondary structure analysis is the main goal of structure determination.     

Ned Seeman and the prediction of amino acid-basepair motifs mediating protein-nucleic acid recognition

Fifty years ago, the first atomic-resolution structure of a nucleic acid double helix, the mini-duplex (ApU)₂, revealed details of basepair geometry, stacking, sugar conformation, and backbone torsion angles, thereby superseding earlier models based on x-ray fiber diffraction, including the original DNA double helix proposed by Watson and Crick. Just 3 years later, in 1976, Ned Seeman, John Rosenberg, and Alex Rich leapt from their structures of mini-duplexes and H-bonding motifs between bases in small-molecule structures and transfer RNA to predicting how proteins could sequence specifically recognize double helix nucleic acids. They proposed interactions between amino acid side chains and nucleobases mediated by two hydrogen bonds in the major or minor grooves. One of these, the arginine-guanine pair, emerged as the most favored amino acid-base interaction in experimental structures of protein-nucleic acid complexes determined since 1986. In this brief review we revisit the pioneering work by Seeman et al. and discuss the importance of the arginine-guanine pairing motif.     

High-throughput single-nucleotide structural mapping by capillary automated footprinting analysis

The use of capillary electrophoresis with fluorescently labeled nucleic acids revolutionized DNA sequencing, effectively fueling the genomic revolution. We present an application of this technology for the high-throughput structural analysis of nucleic acids by chemical and enzymatic mapping ('footprinting'). We achieve the throughput and data quality necessary for genomic-scale structural analysis by combining fluorophore labeling of nucleic acids with novel quantitation algorithms. We implemented these algorithms in the CAFA (capillary automated footprinting analysis) open-source software that is downloadable gratis from https://simtk.org/home/cafa. The accuracy, throughput and reproducibility of CAFA analysis are demonstrated using hydroxyl radical footprinting of RNA. The versatility of CAFA is illustrated by dimethyl sulfate mapping of RNA secondary structure and DNase I mapping of a protein binding to a specific sequence of DNA. Our experimental and computational approach facilitates the acquisition of high-throughput chemical probing data for solution structural analysis of nucleic acids.     

The Potential Utility of Predicted One Bond Carbon-Proton Coupling Constants in the Structure Elucidation of Small Organic Molecules by NMR Spectroscopy

NMR spectroscopy is the most popular technique used for structure elucidation of small organic molecules in solution, but incorrect structures are regularly reported. One-bond proton-carbon J-couplings provide additional information about chemical structure because they are determined by different features of molecular structure than are proton and carbon chemical shifts. However, these couplings are not routinely used to validate proposed structures because few software tools exist to predict them. This study assesses the accuracy of Density Functional Theory for predicting them using 396 published experimental observations from a diverse range of small organic molecules. With the B3LYP functional and the TZVP basis set, Density Functional Theory calculations using the open-source software package NWChem can predict one-bond CH J-couplings with good accuracy for most classes of small organic molecule. The root-mean-square deviation after correction is 1.5 Hz for most sp³ CH pairs and 1.9 Hz for sp² pairs; larger errors are observed for sp³ pairs with multiple electronegative substituents and for sp pairs. These results suggest that prediction of one-bond CH J-couplings by Density Functional Theory is sufficiently accurate for structure validation. This will be of particular use in strained ring systems and heterocycles which have characteristic couplings and which pose challenges for structure elucidation.     

The relationship between amide proton chemical shifts and secondary structure in proteins

Summary The parameters for HN chemical shift calculations of proteins have been determined using data from high-resolution crystal structures of 15 proteins. Employing these chemical shift calculations for HN protons, the observed secondary structure chemical shift trends of HN protons, i.e., upfield shifts on helix formation and downfield shifts on -sheet formation, are discussed. Our calculations suggest that the main reason for the difference in NH chemical shifts in helices and sheets is not an effect from the directly hydrogen-bonded carbonyl, which gives rise to downfield shifts in both cases, but arises from an additional upfield shift predicted in helices and originating in residues i-2 and i-3. The calculations also explain the well-known relationship between amide proton shifts and hydrogen-bond lengths. In addition, the HN chemical shifts of the distorted amphipathic helices of the GCN4 leucine zipper are calculated and used to characterise the solution structure of the helices. By comparing the calculated and experimental shifts, it is shown that in general the agreement is good between residues 15 and 28. The most interesting observation is that in the N-terminal half of the zipper, although both calculated and experimental shifts show clear periodicity, they are no longer in phase. This suggests that for the N-terminal half, in the true average solution structure the period of the helix coil is longer by roughly one residue compared to the NMR structures.     

Identification of helix capping and β-turn motifs from NMR chemical shifts

We present an empirical method for identification of distinct structural motifs in proteins on the basis of experimentally determined backbone and ¹³C^β chemical shifts. Elements identified include the N-terminal and C-terminal helix capping motifs and five types of β-turns: I, II, I′, II′ and VIII. Using a database of proteins of known structure, the NMR chemical shifts, together with the PDB-extracted amino acid preference of the helix capping and β-turn motifs are used as input data for training an artificial neural network algorithm, which outputs the statistical probability of finding each motif at any given position in the protein. The trained neural networks, contained in the MICS (motif identification from chemical shifts) program, also provide a confidence level for each of their predictions, and values ranging from ca 0.7–0.9 for the Matthews correlation coefficient of its predictions far exceed those attainable by sequence analysis. MICS is anticipated to be useful both in the conventional NMR structure determination process and for enhancing on-going efforts to determine protein structures solely on the basis of chemical shift information, where it can aid in identifying protein database fragments suitable for use in building such structures.     

Structure of single-stranded nucleic acids in the presence of ribosomal protein S1.

We have developed a new method for mounting nucleic acids and nucleic acidprotein complexes for high-resolution electron microscopy, and have used it to characterize the interaction between ribosomal protein S1 and single-stranded nucleic acids. We find that SI unwinds most, but not all of the secondary structure present in MS2 RNA and øX174 viral DNA. The binding of S1 to DNA and RNA is not highly co-operative, and has a stoichiometry of one S1 per 10 to 15 nucleotides. We have not observed any tendency for S1 nucleic acid complexes to form aggregates in either 0·01 m-Na⁺ or 0·1 m-Na⁺. An analogous protein isolated from the 30 S ribosomal subunit of Caulobacter crescentus is indistinguishable from Escherichia coli S1 in these studies. The mono-N-ethylmaleimide derivative of E. coli S1 will bind to both MS2 RNA and øX174 viral DNA with a stoichiometry of one N-ethylmaleimide-S1 per 10 to 15 nucleotides, but will not unwind the secondary structure of either of them.     

Structural basis of nucleic acid binding by Nicotiana tabacum glycine-rich RNA-binding protein: implications for its RNA chaperone function

Glycine-rich RNA-binding proteins (GR-RBPs) are involved in cold shock response of plants as RNA chaperones facilitating mRNA transport, splicing and translation. GR-RBPs are bipartite proteins containing a RNA recognition motif (RRM) followed by a glycine-rich region. Here, we studied the structural basis of nucleic acid binding of full-length Nicotiana tabacum GR-RBP1. NMR studies of NtGR-RBP1 show that the glycine-rich domain, while intrinsically disordered, is responsible for mediating self-association by transient interactions with its RRM domain (NtRRM). Both NtGR-RBP1 and NtRRM bind specifically and with low micromolar affinity to RNA and single-stranded DNA. The solution structure of NtRRM shows that it is a canonical RRM domain. A HADDOCK model of the NtRRM–RNA complex, based on NMR chemical shift and NOE data, shows that nucleic acid binding results from a combination of stacking and electrostatic interactions with conserved RRM residues. Finally, DNA melting experiments demonstrate that NtGR-RBP1 is more efficient in melting CTG containing nucleic acids than isolated NtRRM. Together, our study supports the model that self-association of GR-RBPs by the glycine-rich region results in cooperative unfolding of non-native substrate structures, thereby enhancing its chaperone function.     

Thermal difference spectra: a specific signature for nucleic acid structures

We show that nucleic acid structures may be conveniently and inexpensively characterized by their UV thermal difference spectra. A thermal difference spectrum (TDS) is obtained for a nucleic acid by simply recording the ultraviolet absorbance spectra of the unfolded and folded states at temperatures above and below its melting temperature (T_m). The difference between these two spectra is the TDS. The TDS has a specific shape that is unique for each type of nucleic acid structure, a conclusion that is based on a comparison of >900 spectra from 200 different sequences. The shape of the TDS reflects the subtleties of base stacking interactions that occur uniquely within each type of nucleic acid structure. TDS provides a simple, inexpensive and rapid method to obtain structural insight into nucleic acid structures, which is applicable to both DNA and RNA from short oligomers to polynucleotides. TDS complements circular dichroism as a tool for the structural characterization of nucleic acids in solution.     

Assignment of resonances in the Escherichia coli 5 S RNA fragment proton NMR spectrum using uniform nitrogen-15 enrichment

The downfield proton NMR spectrum of aqueous uniformly nitrogen-15 enriched 5 S RNA fragment is presented. Selective nitrogen-15 decoupling difference proton spectroscopy revealed nitrogen-15 chemical shifts of fragment imino nitrogens. Nitrogen chemical shifts of nucleic acid guanine and uracil imino nitrogens have separate small ranges. Nitrogen-15 and proton chemical shift correlation by the heteronuclear decoupling permitted the identification of the base type of some previously unassigned imino proton resonances in the 5 S RNA fragment spectrum. Corresponding resonances in the natural isotopic abundance 5 S RNA fragment spectrum are assigned to base types by comparison with the enriched sample spectrum.     

De novo protein structure generation from incomplete chemical shift assignments

NMR chemical shifts provide important local structural information for proteins. Consistent structure generation from NMR chemical shift data has recently become feasible for proteins with sizes of up to 130 residues, and such structures are of a quality comparable to those obtained with the standard NMR protocol. This study investigates the influence of the completeness of chemical shift assignments on structures generated from chemical shifts. The Chemical-Shift-Rosetta (CS-Rosetta) protocol was used for de novo protein structure generation with various degrees of completeness of the chemical shift assignment, simulated by omission of entries in the experimental chemical shift data previously used for the initial demonstration of the CS-Rosetta approach. In addition, a new CS-Rosetta protocol is described that improves robustness of the method for proteins with missing or erroneous NMR chemical shift input data. This strategy, which uses traditional Rosetta for pre-filtering of the fragment selection process, is demonstrated for two paramagnetic proteins and also for two proteins with solid-state NMR chemical shift assignments. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Heteroassociation of antibiotic norfloxacin with aromatic vitamins in aqueous solution

Heteroassociation of antibacterial antibiotic norfloxacin with aromatic vitamins nicotinamide and flavin mononucleotide in aqueous solution was studied by ¹H NMR spectroscopy (500 MHz). Equilibrium constants, induced proton chemical shifts, and thermodynamic parameters (ΔH, ΔS) for the reactions of heteroassociation of the molecules were determined on the basis of the concentration and temperature dependences of proton chemical shifts for interacting aromatic molecules. The analysis of the results obtained indicates the formation of heterocomplexes between vitamin molecules and norfloxacin owing to stacking interactions between aromatic chromophores and additional intermolecular hydrogen bonding in norfloxacin-nicotinamide. The most probable spatial structures of 1:1 norfloxacin-flavin mononucleotide and norfloxacin-nicotinamide heterocomplexes were determined by molecular modeling methods using X-PLOR software on the basis of analysis of induced proton chemical shifts.