Reparameterization of RNA χ Torsion Parameters for the AMBER Force Field and Comparison to NMR Spectra for Cytidine and Uridine

A reparameterization of the torsional parameters for the glycosidic dihedral angle, χ, for the AMBER99 force field in RNA nucleosides is used to provide a modified force field, AMBER99χ. Molecular dynamics simulations of cytidine, uridine, adenosine, and guanosine in aqueous solution using the AMBER99 and AMBER99χ force fields are compared with NMR results. For each nucleoside and force field, 10 individual molecular dynamics simulations of 30 ns each were run. For cytidine with AMBER99χ force field, each molecular dynamics simulation time was extended to 120 ns for convergence purposes. Nuclear magnetic resonance (NMR) spectroscopy, including one-dimensional (1D) 1H, steady-state 1D 1H nuclear Overhauser effect (NOE), and transient 1D 1H NOE, was used to determine the sugar puckering and preferred base orientation with respect to the ribose of cytidine and uridine. The AMBER99 force field overestimates the population of syn conformations of the base orientation and of C2′-endo sugar puckering of the pyrimidines, while the AMBER99χ force field’s predictions are more consistent with NMR results. Moreover, the AMBER99 force field prefers high anti conformations with glycosidic dihedral angles around 310° for the base orientation of purines. The AMBER99χ force field prefers anti conformations around 185°, which is more consistent with the quantum mechanical calculations and known 3D structures of folded ribonucleic acids (RNAs). Evidently, the AMBER99χ force field predicts the structural characteristics of ribonucleosides better than the AMBER99 force field and should improve structural and thermodynamic predictions of RNA structures

ATP and its N6-substituted analogues: parameterization, molecular dynamics simulation and conformational analysis

In this work we used a combination of classical molecular dynamics and simulated annealing techniques to shed more light on the conformational flexibility of 12 adenosine triphosphate (ATP) analogues in a water environment. We present simulations in AMBER force field for ATP and 12 published analogues [Shah et al. (1997) Proc Natl Acad Sci USA 94: 3565–3570]. The calculations were carried out using the generalized Born (GB) solvation model in the presence of the cation Mg2+. The ion was placed at a close distance (2 Å) from the charged oxygen atoms of the beta and gamma phosphate groups of the −3 negatively charged ATP analogue molecules. Analysis of the results revealed the distribution of inter-proton distances H8–H1′ and H8–H2′ versus the torsion angle ψ (C4–N9-C1′–O4′) for all conformations of ATP analogues. There are two gaps in the distribution of torsion angle ψ values: the first is between −30 and 30 degrees and is described by cis-conformation; and the second is between 90 and 175 degrees, which mostly covers a region of anti conformation. Our results compare favorably with results obtained in experimental assays [Jiang and Mao (2002) Polyhedron 21:435–438]

Recent advances in the development and evaluation of molecular diagnostics for Ebola virus disease

The 2014-16 outbreak of ebola virus disease (EVD) in West Africa resulted in 11,308 deaths. During the outbreak only 60% of patients were laboratory confirmed and global health authorities have identified the need for accurate and readily deployable molecular diagnostics as an important component of the ideal response to future outbreaks, to quickly identify and isolate patients. Areas covered: Currently PCR-based techniques and rapid diagnostic tests (RDTs) that detect antigens specific to EVD infections dominate the diagnostic landscape, but recent advances in biosensor technologies have led to novel approaches for the development of EVD diagnostics. This review summarises the literature and available performance data of currently available molecular diagnostics for ebolavirus, identifies knowledge gaps and maps out future priorities for research in this field. Expert opinion: While there are now a plethora of diagnostic tests for EVD at various stages of development, there is an acute need for studies to compare their clinical performance, but the sporadic nature of EVD outbreaks makes this extremely challenging, demanding pragmatic new modalities of research funding and ethical/institutional approval, to enable responsive research in outbreak settings. Retrospective head-to-head diagnostic comparisons could also be implemented using biobanked specimens, providing this can be done safely

Characterizing Structural Transitions Using Localized Free Energy Landscape Analysis

Structural changes in molecules are frequently observed during biological processes like replication, transcription and translation. These structural changes can usually be traced to specific distortions in the backbones of the macromolecules involved. Quantitative energetic characterization of such distortions can greatly advance the atomic-level understanding of the dynamic character of these biological processes.Molecular dynamics simulations combined with a variation of the Weighted Histogram Analysis Method for potential of mean force determination are applied to characterize localized structural changes for the test case of cytosine (underlined) base flipping in a GTCAGCGCATGG DNA duplex. Free energy landscapes for backbone torsion and sugar pucker degrees of freedom in the DNA are used to understand their behavior in response to the base flipping perturbation. By simplifying the base flipping structural change into a two-state model, a free energy difference of upto 14 kcal/mol can be attributed to the flipped state relative to the stacked Watson-Crick base paired state. This two-state classification allows precise evaluation of the effect of base flipping on local backbone degrees of freedom.The calculated free energy landscapes of individual backbone and sugar degrees of freedom expectedly show the greatest change in the vicinity of the flipping base itself, but specific delocalized effects can be discerned upto four nucleotide positions away in both 5' and 3' directions. Free energy landscape analysis thus provides a quantitative method to pinpoint the determinants of structural change on the atomic scale and also delineate the extent of propagation of the perturbation along the molecule. In addition to nucleic acids, this methodology is anticipated to be useful for studying conformational changes in all macromolecules, including carbohydrates, lipids, and proteins

Aberrant Cyclization Affords a C-6 Modified Cyclic Adenosine 5′-Diphosphoribose Analogue with Biological Activity in Jurkat T Cells

*S Supporting Information ABSTRACT: Two nicotinamide adenine dinucleotide (NAD +) analogues modified at the 6 position of the purine ring were synthesized, and their substrate properties toward Aplysia californica ADP-ribosyl cyclase were investigated. 6-N-Methyl NAD + (6-N-methyl nicotinamide adenosine 5′-dinucleotide 10) hydrolyzes to give the linear 6-N-methyl ADPR (adenosine 5′-diphosphoribose, 11), whereas 6-thio NHD + (nicotinamide 6-mercaptopurine 5′-dinucleotide, 17) generates a cyclic dinucleotide. Surprisingly, NMR correlation spectra confirm this compound to be the N1 cyclic product 6-thio N1-cIDPR (6-thio cyclic inosine 5′-diphosphoribose, 3), although the corresponding 6-oxo analogue is well-known to cyclize at N7. In Jurkat T cells, unlike the parent cyclic inosine 5′-diphosphoribose N1-cIDPR 2, 6-thio N1-cIDPR antagonizes both cADPR- and N1cIDPR-induced Ca 2+ release but possesses weak agonist activity at higher concentration. 3 is thus identified as the first C-6 modified cADPR (cyclic adenosine 5′-diphosphoribose) analogue antagonist; it represents the first example of a fluorescent N1cyclized cADPR analogue and is a new pharmacological tool for intervention in the cADPR pathway of cellular signaling

Bulge-out structures in the single-stranded trimer AUA and in the duplex (CUGGUGCGG).(CCGCCCAG). A model-building and NMR study.

Model-building studies were carried out on the trimer AUA. Bulge-out structures which allow incorporation into a continuous RNA helix were generated and energy-minimized. All geometrical features obtained by previous NMR studies on purine-pyrimidine-purine sequences are accounted for in these models. One of the models was used to fit into a double helical fragment. Only minor changes were necessary to construct a central bulge-out in an otherwise intact duplex. NMR and model-building studies were performed on the duplex (CUGGUGCGG).(CCGCCCAG) which contains an unpaired uridine residue. NOE data, chemical-shift profiles and imino-proton resonances provided evidence that the extra U is bulged out of the duplex. The relatively small dispersion in 31P chemical shifts (approximately equal to 0.7 ppm) indicate the absence of t/g or g/t combinations for the phosphodiester angles zeta/alpha. An energy-minimized model of the duplex, which fits the present collection of data, is presented