Reparameterization of RNA chi 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, chi, for the AMBER99 force field in RNA nucleosides is used to provide a modified force field, AMBER99chi. Molecular dynamics simulations of cytidine, uridine, adenosine, and guanosine in aqueous solution using the AMBER99 and AMBER99chi 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 AMBER99chi force field, each molecular dynamics simulation time was extended to 120 ns for convergence purposes. Nuclear magnetic resonance (NMR) spectroscopy, including one-dimensional (1D) (1)H, steady-state 1D (1)H nuclear Overhauser effect (NOE), and transient 1D (1)H 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 AMBER99chi 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 degrees for the base orientation of purines. The AMBER99chi force field prefers anti conformations around 185 degrees , which is more consistent with the quantum mechanical calculations and known 3D structures of folded ribonucleic acids (RNAs). Evidently, the AMBER99chi force field predicts the structural characteristics of ribonucleosides better than the AMBER99 force field and should improve structural and thermodynamic predictions of RNA structures.

TRANSCRIPTION. Allosteric transcriptional regulation via changes in the overall topology of the core promoter

Many transcriptional activators act at a distance from core promoter elements and work by recruiting RNA polymerase through protein-protein interactions. We show here how the prokaryotic regulatory protein CueR both represses and activates transcription by differentially modulating local DNA structure within the promoter. Structural studies reveal that the repressor state slightly bends the promoter DNA, precluding optimal RNA polymerase-promoter recognition. Upon binding a metal ion in the allosteric site, CueR switches into an activator conformation. It maintains all protein-DNA contacts but introduces torsional stresses that kink and undertwist the promoter, stabilizing an A-form DNA-like conformation. These factors switch on and off transcription by exerting dynamic control of DNA stereochemistry, reshaping the core promoter and making it a better or worse substrate for polymerase.

Stacking in RNA: NMR of Four Tetramers Benchmark Molecular Dynamics

Molecular dynamics (MD) simulations for RNA tetramers r(AAAA), r(CAAU), r(GACC), and r(UUUU) are benchmarked against 1H-1H NOESY distances and 3J scalar couplings to test effects of RNA torsion parametrizations. Four different starting structures were used for r(AAAA), r(CAAU), and r(GACC), while five starting structures were used for r(UUUU). On the basis of X-ray structures, criteria are reported for quantifying stacking. The force fields, AMBER ff99, parmbsc0, parm99χ_Yil, ff10, and parmTor, all predict experimentally unobserved stacks and intercalations, e.g., base 1 stacked between bases 3 and 4, and incorrect χ, ϵ, and sugar pucker populations. The intercalated structures are particularly stable, often lasting several microseconds. Parmbsc0, parm99χ_Yil, and ff10 give similar agreement with NMR, but the best agreement is only 46%. Experimentally unobserved intercalations typically are associated with reduced solvent accessible surface area along with amino and hydroxyl hydrogen bonds to phosphate nonbridging oxygens. Results from an extensive set of MD simulations suggest that recent force field parametrizations improve predictions, but further improvements are necessary to provide reasonable agreement with NMR. In particular, intramolecular stacking and hydrogen bonding interactions may not be well balanced with the TIP3P water model. NMR data and the scoring method presented here provide rigorous benchmarks for future changes in force fields and MD methods.

Benchmarking AMBER force fields for RNA: comparisons to NMR spectra for single-stranded r(GACC) are improved by revised χ torsions

Accurately modeling unpaired regions of RNA is important for predicting structure, dynamics, and thermodynamics of folded RNA. Comparisons between NMR data and molecular dynamics simulations provide a test of force fields used for modeling. Here, NMR spectroscopy, including NOESY, (1)H-(31)P HETCOR, DQF-COSY, and TOCSY, was used to determine conformational preferences for single-stranded GACC RNA. The spectra are consistent with a conformational ensemble containing major and minor A-form-like structures. In a series of 50 ns molecular dynamics (MD) simulations with the AMBER99 force field in explicit solvent, initial A-form-like structures rapidly evolve to disordered conformations. A set of 50 ns simulations with revised χ torsions (AMBER99χ force field) gives two primary conformations, consistent with the NMR spectra. A single 1.9 μs MD simulation with the AMBER99χ force field showed that the major and minor conformations are retained for almost 68% of the time in the first 700 ns, with multiple transformations from A-form to non-A-form conformations. For the rest of the simulation, random-coil structures and a stable non-A-form conformation inconsistent with NMR spectra were seen. Evidently, the AMBER99χ force field improves structural predictions for single-stranded GACC RNA compared to the AMBER99 force field, but further force field improvements are needed.

The Hairpin Form of r(G4C2)exp in c9ALS/FTD Is Repeat-Associated Non-ATG Translated and a Target for Bioactive Small Molecules

The most common genetic cause of amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD) is an expanded G₄C₂ repeat [(G₄C₂)^exp] in C9ORF72. ALS/FTD-associated toxicity has been traced to the RNA transcribed from the repeat expansion [r(G₄C₂)^exp], which sequesters RNA-binding proteins (RBPs) and undergoes repeat-associated non-ATG (RAN) translation to generate toxic dipeptide repeats. Using in vitro and cell-based assays, we identified a small molecule (4) that selectively bound r(G₄C₂)^exp, prevented sequestration of an RBP, and inhibited RAN translation. Indeed, biophysical characterization showed that 4 selectively bound the hairpin form of r(G₄C₂)^exp, and nuclear magnetic resonance spectroscopy studies and molecular dynamics simulations defined this molecular recognition event. Cellular imaging revealed that 4 localized to r(G₄C₂)^exp cytoplasmic foci, the putative sites of RAN translation. Collectively, these studies highlight that the hairpin structure of r(G₄C₂)^exp is a therapeutically relevant target and small molecules that bind it can ameliorate c9ALS/FTD-associated toxicity.

Testing the nearest neighbor model for canonical RNA base pairs: revision of GU parameters

Thermodynamic parameters for GU pairs are important for predicting the secondary structures of RNA and for finding genomic sequences that code for structured RNA. Optical melting curves were measured for 29 RNA duplexes with GU pairs to improve nearest neighbor parameters for predicting stabilities of helixes. The updated model eliminates a prior penalty assumed for terminal GU pairs. Six additional duplexes with the 5'GG/3'UU motif were added to the single representation in the previous database. This revises the ΔG°(37) for the 5'GG/3'UU motif from an unfavorable 0.5 kcal/mol to a favorable -0.2 kcal/mol. Similarly, the ΔG°(37) for the 5'UG/3'GU motif changes from 0.3 to -0.6 kcal/mol. The correlation coefficients between predicted and experimental ΔG°(37), ΔH°, and ΔS° for the expanded database are 0.95, 0.89, and 0.87, respectively. The results should improve predictions of RNA secondary structure.

Resistance to extended-spectrum cephalosporins, caused by PER-1 beta-lactamase, in Salmonella typhimurium from Istanbul, Turkey

Two Salmonella typhimurium isolates were studied, one as a representative from a series of neonatal meningitis cases treated at an Istanbul teaching hospital, the other from a gastro-enteritis case seen at a different Istanbul hospital. Both isolates were resistant to extended-spectrum cephalosporins, as well as penicillins, aminoglycosides and chloramphenicol. Cephalosporin resistance depended on production of PER-1 beta-lactamase, which is an extended-spectrum class A enzyme that is only distantly related to TEM and SHV enzymes, and which was previously known only from Pseudomonas aeruginosa isolates. The PER-1 gene was carried by an 81-MDa plasmid, which also determined resistance to aminoglycosides and chloramphenicol. Although it was not self-transmissible to Escherichia coli, this element did transfer if mobilised with plasmid pUZ8. The two S. typhimurium isolates gave indistinguishable DNA restriction patterns and, in addition to their 81-MDa plasmid, also contained 52- and 2.8-MDa plasmids, the last of these encoded TEM-1 enzyme. The two isolates were identical in serotype, antibiogram and plasmid-profile but nevertheless differed in phage type, and, therefore, represented distinct strains. The emergence of cefotaxime and ceftriaxone resistance in salmonellae is disturbing, since these agents are preferred therapy for neonatal meningitis caused by members of the genus.

Revision of AMBER Torsional Parameters for RNA Improves Free Energy Predictions for Tetramer Duplexes with GC and iGiC Base Pairs

All-atom force fields are important for predicting thermodynamic, structural, and dynamic properties of RNA. In this paper, results are reported for thermodynamic integration calculations of free energy differences of duplex formation when CG pairs in the RNA duplexes r(CCGG)(2), r(GGCC)(2), r(GCGC)(2), and r(CGCG)(2) are replaced by isocytidine-isoguanosine (iCiG) pairs. Agreement with experiment was improved when ε/ζ, α/γ, β, and χ torsional parameters in the AMBER99 force field were revised on the basis of quantum mechanical calculations. The revised force field, AMBER99TOR, brings free energy difference predictions to within 1.3, 1.4, 2.3, and 2.6 kcal/mol at 300 K, respectively, compared to experimental results for the thermodynamic cycles of CCGG → iCiCiGiG, GGCC → iGiGiCiC, GCGC → iGiCiGiC, and CGCG → iCiGiCiG. In contrast, unmodified AMBER99 predictions for GGCC → iGiGiCiC and GCGC → iGiCiGiC differ from experiment by 11.7 and 12.6 kcal/mol, respectively. In order to test the dynamic stability of the above duplexes with AMBER99TOR, four individual 50 ns molecular dynamics (MD) simulations in explicit solvent were run. All except r(CCGG)(2) retained A-form conformation for ≥82% of the time. This is consistent with NMR spectra of r(iGiGiCiC)(2), which reveal an A-form conformation. In MD simulations, r(CCGG)(2) retained A-form conformation 52% of the time, suggesting that its terminal base pairs may fray. The results indicate that revised backbone parameters improve predictions of RNA properties and that comparisons to measured sequence dependent thermodynamics provide useful benchmarks for testing force fields and computational methods.

Structure of the myotonic dystrophy type 2 RNA and designed small molecules that reduce toxicity

Myotonic dystrophy type 2 (DM2) is an incurable neuromuscular disorder caused by a r(CCUG) expansion (r(CCUG)(exp)) that folds into an extended hairpin with periodically repeating 2×2 nucleotide internal loops (5'CCUG/3'GUCC). We designed multivalent compounds that improve DM2-associated defects using information about RNA-small molecule interactions. We also report the first crystal structure of r(CCUG) repeats refined to 2.35 Å. Structural analysis of the three 5'CCUG/3'GUCC repeat internal loops (L) reveals that the CU pairs in L1 are each stabilized by one hydrogen bond and a water-mediated hydrogen bond, while CU pairs in L2 and L3 are stabilized by two hydrogen bonds. Molecular dynamics (MD) simulations reveal that the CU pairs are dynamic and stabilized by Na(+) and water molecules. MD simulations of the binding of the small molecule to r(CCUG) repeats reveal that the lowest free energy binding mode occurs via the major groove, in which one C residue is unstacked and the cross-strand nucleotides are displaced. Moreover, we modeled the binding of our dimeric compound to two 5'CCUG/3'GUCC motifs, which shows that the scaffold on which the RNA-binding modules are displayed provides an optimal distance to span two adjacent loops.

Methods to enable the design of bioactive small molecules targeting RNA

RNA is an immensely important target for small molecule therapeutics or chemical probes of function. However, methods that identify, annotate, and optimize RNA-small molecule interactions that could enable the design of compounds that modulate RNA function are in their infancies. This review describes recent approaches that have been developed to understand and optimize RNA motif-small molecule interactions, including structure-activity relationships through sequencing (StARTS), quantitative structure-activity relationships (QSAR), chemical similarity searching, structure-based design and docking, and molecular dynamics (MD) simulations. Case studies described include the design of small molecules targeting RNA expansions, the bacterial A-site, viral RNAs, and telomerase RNA. These approaches can be combined to afford a synergistic method to exploit the myriad of RNA targets in the transcriptome.

A dynamic structural model of expanded RNA CAG repeats: a refined X-ray structure and computational investigations using molecular dynamics and umbrella sampling simulations

One class of functionally important RNA is repeating transcripts that cause disease through various mechanisms. For example, expanded CAG repeats can cause Huntington's and other disease through translation of toxic proteins. Herein, a crystal structure of r[5'UUGGGC(CAG)3GUCC]2, a model of CAG expanded transcripts, refined to 1.65 Å resolution is disclosed that shows both anti-anti and syn-anti orientations for 1 × 1 nucleotide AA internal loops. Molecular dynamics (MD) simulations using AMBER force field in explicit solvent were run for over 500 ns on the model systems r(5'GCGCAGCGC)2 (MS1) and r(5'CCGCAGCGG)2 (MS2). In these MD simulations, both anti-anti and syn-anti AA base pairs appear to be stable. While anti-anti AA base pairs were dynamic and sampled multiple anti-anti conformations, no syn-anti ↔ anti-anti transformations were observed. Umbrella sampling simulations were run on MS2, and a 2D free energy surface was created to extract transformation pathways. In addition, an explicit solvent MD simulation over 800 ns was run on r[5'GGGC(CAG)3GUCC]2, which closely represents the refined crystal structure. One of the terminal AA base pairs (syn-anti conformation), transformed to anti-anti conformation. The pathway followed in this transformation was the one predicted by umbrella sampling simulations. Further analysis showed a binding pocket near AA base pairs in syn-anti conformations. Computational results combined with the refined crystal structure show that global minimum conformation of 1 × 1 nucleotide AA internal loops in r(CAG) repeats is anti-anti but can adopt syn-anti depending on the environment. These results are important to understand RNA dynamic-function relationships and to develop small molecules that target RNA dynamic ensembles.

Interplay of LNA and 2'-O-methyl RNA in the structure and thermodynamics of RNA hybrid systems: a molecular dynamics study using the revised AMBER force field and comparison with experimental results

When used in nucleic acid duplexes, locked nucleic acid (LNA) and 2'-O-methyl RNA residues enhance the duplex stabilities, and this makes it possible to create much better RNA aptamers to target specific molecules in cells. Thus, LNA and 2'-O-methyl RNA residues are finding increasingly widespread use in RNA-based therapeutics. Herein, we utilize molecular dynamics (MD) simulations and UV melting experiments to investigate the structural and thermodynamic properties of 13 nucleic acid duplexes, including full DNA, RNA, LNA, and 2'-O-methyl RNA duplexes as well as hybrid systems such as LNA:RNA, 2'-O-methyl RNA:RNA, LNA/2'-O-methyl RNA:RNA, and RNA/2'-O-methyl RNA:RNA duplexes. The MD simulations are based on a version of the Amber force field revised specifically for RNA and LNA residues. Our results indicate that LNA and 2'-O-methyl RNA residues have two different hybridization mechanisms when included in hybrid duplexes with RNA wherein the former underwinds while the latter overwinds the duplexes. These computational predictions are supported by X-ray structures of LNA and 2'-O-methyl RNA duplexes that were recently presented by different groups, and there is also good agreement with the measured thermal stabilities of the duplexes. We find out that the "underwinding" phenomenon seen in LNA and LNA:RNA hybrid duplexes happens due to expansion of the major groove widths (Mgw) of the duplexes that is associated with decrease in the slide and twist values in base-pair steps. In contrast, 2'-O-methyl RNA residues in RNA duplexes slightly overwind the duplexes while the backbone is forced to stay in C3'-endo. Moreover, base-pair stacking in the LNA and LNA:RNA hybrid systems is gradually reduced with the inclusion of LNA residues in the duplexes while no such effect is seen in the 2'-O-methyl RNA systems. Our results show how competition between base stacking and structural rigidity in these RNA hybrid systems influences structures and stabilities. Even though both LNA and 2'-O-methyl RNA residues have C3'-endo sugar puckering, structurally LNA residues have a frozen sugar backbone which provides entropic enhancement of stabilities while the 2'-O-methyl RNA residues are more flexible and maintain base stacking that is almost untouched compared to RNA. Thus, enhancement of the structural stabilities of RNA duplexes by 2'-O-methyl RNA modifications is smaller than for the corresponding LNA modifications. Indeed, our experimental measurements show that on average each 2'-O-methyl RNA and LNA substitution in a RNA duplex enhances duplex stability by 0.2 and 1.4 kcal/mol, respectively. Our computational binding free energy predictions are qualitatively in line with these results. The only exception is for the full 2'-O-methyl RNA duplex, which is overstabilized, implying that further force field revisions are needed. Collectively, the results presented in this paper explain the atomistic details of the structural and thermodynamic roles of LNA and 2'-O-methyl RNA residues in RNA hybrid duplexes, shedding light on the mechanism behind targeting endogenous micro RNA (miRNA) in order to regulate mRNA activity and inhibit gene expression in the cell.

Characterization of multiple-antibiotic-resistant Salmonella typhimurium stains: molecular epidemiology of PER-1-producing isolates and evidence for nosocomial plasmid exchange by a clone

We characterized epidemiologic and genetic features of nosocomially originated multiple-antibiotic-resistant Salmonella typhimurium isolates from two hospitals. A total of 32 multiply resistant strains, isolated during a 28-month period, were studied. Four resistance phenotypes were distinguished on the basis of the results of disc diffusion tests. Group 1 was resistant to chloramphenicol, gentamicin, tobramycin, amikacin, and the newer cephalosporins because of the production of an extended-spectrum beta-lactamase (PER-1). Group 2 exhibited the same pattern plus resistance to sulfamethoxazole-trimethoprim (Sxt). Except for Sxt resistance, dominant phenotypes of both groups were transferred on an identical plasmid, pSTI1 (81 MDa). Group 3 was resistant to ampicillin, chloramphenicol, gentamicin, tobramycin, and Sxt. This pattern was also transferred on an 81-MDa plasmid (pSTI2) which differed from pSTI1 on the basis of EcoRI and HindIII restriction fragments. Group 4 was resistant to ampicillin, chloramphenicol, and tetracycline, and a 74-MDa nonconjugative plasmid was detected. Restriction fragment length polymorphism of RNA-encoding DNA and arbitrarily primed PCR tests revealed that bacteria from groups 1, 2, and 3 were clonally related. Epidemiologic data also supported the clonal-dissemination hypothesis. We concluded that S. typhimurium isolates acquire and exchange multiple-resistance plasmids in hospital microflora.

An RNA Molecular Switch: Intrinsic Flexibility of 23S rRNA Helices 40 and 68 5'-UAA/5'-GAN Internal Loops Studied by Molecular Dynamics Methods

Functional RNA molecules such as ribosomal RNAs (rRNAs) frequently contain highly conserved internal loops with a 5'-UAA/5'-GAN (UAA/GAN) consensus sequence. The UAA/GAN internal loops adopt a distinctive structure inconsistent with secondary structure predictions. The structure has a narrow major groove and forms a trans Hoogsteen/Sugar edge (tHS) A/G base pair followed by an unpaired stacked adenine, a trans Watson-Crick/Hoogsteen (tWH) U/A base pair, and finally a bulged nucleotide (N). The structure is further stabilized by a three-adenine stack and base-phosphate interaction. In the ribosome, the UAA/GAN internal loops are involved in extensive tertiary contacts, mainly as donors of A-minor interactions. Further, this sequence can adopt an alternative 2D/3D pattern stabilized by a four-adenine stack involved in a smaller number of tertiary interactions. The solution structure of an isolated UAA/GAA internal loop shows substantially rearranged base pairing with three consecutive non-Watson-Crick base pairs. Its A/U base pair adopts an incomplete cis Watson-Crick/Sugar edge (cWS) A/U conformation instead of the expected Watson-Crick arrangement. We performed 3.1 μs of explicit solvent molecular dynamics (MD) simulations of the X-ray and NMR UAA/GAN structures, supplemented by molecular mechanics, Poisson-Boltzmann, and surface area free energy calculations; locally enhanced sampling (LES) runs; targeted MD (TMD); and nudged elastic band (NEB) analysis. We compared parm99 and parmbsc0 force fields and net-neutralizing Na(+) versus excess salt KCl ion environments. Both force fields provide a similar description of the simulated structures, with the parmbsc0 leading to modest narrowing of the major groove. The excess salt simulations also cause a similar effect. While the NMR structure is entirely stable in simulations, the simulated X-ray structure shows considerable widening of the major groove, a loss of base-phosphate interaction, and other instabilities. The alternative X-ray geometry even undergoes a conformational transition toward the solution 2D structure. Free energy calculations confirm that the X-ray arrangement is less stable than the solution structure. LES, TMD, and NEB provide a rather consistent pathway for interconversion between the X-ray and NMR structures. In simulations, the incomplete cWS A/U base pair of the NMR structure is water-mediated and alternates with the canonical A-U base pair, which is not indicated by the NMR data. Completion of the full cWS A/U base pair is prevented by the overall internal loop arrangement. In summary, the simulations confirm that the UAA/GAN internal loop is a molecular switch RNA module that adopts its functional geometry upon specific tertiary contexts.

Contributions of stacking, preorganization, and hydrogen bonding to the thermodynamic stability of duplexes between RNA and 2'-O-methyl RNA with locked nucleic acids

Locked nucleic acids (LNA) considerably enhance the thermodynamic stability of DNA and RNA duplexes. We report the thermodynamic stabilities of LNA-2'-O-methyl RNA/RNA duplexes designed to provide insight into the contributions of stacking and hydrogen bonding interactions to the enhanced stability. The results show that hydrogen bonding of LNA nucleotides is similar to that of 2'-O-methyl RNA nucleotides, whereas the 3'-stacking interactions are on average approximately 0.7 kcal/mol more favorable at 37 degrees C than for 2'-O-methyl or RNA nucleotides. Moreover, NMR spectra suggest helical preorganization of the single-stranded tetramer, C(L)A(M)A(L)U(M), probably due to restriction of some torsion angles. Thus, enhanced stacking interactions and helical preorganization of single-stranded oligonucleotides contribute to the extraordinary stabilization of duplexes by LNA nucleotides.

RNA challenges for computational chemists

Some experimental results for the thermodynamics of RNA folding cannot be explained by simple pairwise hydrogen-bonding models. Such effects include the stabilities of isoguanosine-isocytidine (iG-iC) base pairs and of various 2 x 2 nucleotide internal loops. Presumably, these results can be explained by base stacking effects, which can be partitioned into Coulombic and overlap effects. We review experimental measurements that provide benchmarks for testing the approximations and theories used for modeling nucleic acids. Quantitative agreement between experiment and theory will indicate understanding of the interactions determining RNA stability and structure.

Surveillance of pneumococcal colonization and invasive pneumococcal disease reveals shift in prevalent carriage serotypes in Massachusetts' children to relatively low invasiveness

BACKGROUND

Following the introduction of pneumococcal conjugate vaccines (PCV), overall nasopharyngeal colonization rates have not changed significantly, however a dramatic and sustained decline in invasive pneumococcal disease (IPD) in children was observed in every setting where the PCVs were implemented. We aimed to describe the differences in invasive disease potential of serotypes that are common colonizers in pre- and post-vaccine eras in order to provide further insight in our understanding of dynamic epidemiology of pneumococcal diseases.

METHODS

Using data from surveillance of nasopharyngeal carriage and enhanced surveillance for IPD, a serotype specific "invasive capacity (IC)" was computed by dividing the incidence of IPD due to serotype x by the carriage prevalence of that same serotype in children <7years of age in Massachusetts. We have evaluated the serotype specific invasive capacity in two periods; pre-PCV13 (2001/02, 2003/04, 2006/07, 2008/09) and post-PCV13 (2010/11 and 2013/14), and by age groups; <24monthsvs. ≥24months.

RESULTS

An approximate 50-fold variation in the point estimate was observed between the serotypes having the highest (7F, 38, 19A, 3, 33F) and the lowest (6C, 35B, 21, 11A, 23B and 23A) computed serotype specific invasive disease potential. In the post-PCV13 era (6C, 35B, 11A, 23B and 23A), 5 of the 7 most common serotypes colonizing the nasopharynx were serotypes with the lowest invasive capacity. Serotype specific invasive capacity trended down in older children for majority of the serotypes, and serotypes 3, 10A and 19A had significantly lower invasive capacity in children older than 24months of age compared to younger children.

CONCLUSION

Invasive capacity differs among serotypes and likely by age. Point estimates of IC for most of the common serotypes colonizing children in Massachusetts in post-PCV13 era were low and likely explain the continued reduction in IPD from the pre-PCV era in the absence of specific protection against these serotypes.

Effects of Restrained Sampling Space and Nonplanar Amino Groups on Free-Energy Predictions for RNA with Imino and Sheared Tandem GA Base Pairs Flanked by GC, CG, iGiC or iCiG Base Pairs

Guanine-adenine (GA) base pairs play important roles in determining the structure, dynamics, and stability of RNA. In RNA internal loops, GA base pairs often occur in tandem arrangements and their structure is context and sequence dependent. Calculations reported here test the thermodynamic integration (TI) approach with the amber99 force field by comparing computational predictions of free energy differences with the free energy differences expected on the basis of NMR determined structures of the RNA motifs (5′-GCGGACGC-3′)₂, (5′-GCiGGAiCGC-3′)₂, (5′-GGCGAGCC-3′)₂, and (5′-GGiCGAiGCC-3′)₂. Here, iG and iC denote isoguanosine and isocytidine, which have amino and carbonyl groups transposed relative to guanosine and cytidine. The NMR structures show that the GA base pairs adopt either imino (cis Watson−Crick/Watson−Crick A-G) or sheared (trans Hoogsteen/Sugar edge A-G) conformations depending on the identity and orientation of the adjacent base pair. A new mixing function for the TI method is developed that allows alchemical transitions in which atoms can disappear in both the initial and final states. Unrestrained calculations gave ΔG° values 2−4 kcal/mol different from expectations based on NMR data. Restraining the structures with hydrogen bond restraints did not improve the predictions. Agreement with NMR data was improved by 0.7 to 1.5 kcal/mol, however, when structures were restrained with weak positional restraints to sample around the experimentally determined NMR structures. The amber99 force field was modified to partially include pyramidalization effects of the unpaired amino group of guanosine in imino GA base pairs. This provided little or no improvement in comparisons with experiment. The marginal improvement is observed when the structure has potential cross-strand out-of-plane hydrogen bonding with the G amino group. The calculations using positional restraints and a nonplanar amino group reproduce the signs of ΔG° from the experimental results and are, thus, capable of providing useful qualitative insights complementing the NMR experiments. Decomposition of the terms in the calculations reveals that the dominant terms are from electrostatic and interstrand interactions other than hydrogen bonds in the base pairs. The results suggest that a better description of the backbone is key to reproducing the experimental free energy results with computational free energy predictions.

Ribonuclease recruitment using a small molecule reduced c9ALS/FTD r(G4C2) repeat expansion in vitro and in vivo ALS models

[Figure: see text].

Optimization of an AMBER force field for the artificial nucleic acid, LNA, and benchmarking with NMR of L(CAAU)

Locked Nucleic Acids (LNAs) are RNA analogues with an O2'-C4' methylene bridge which locks the sugar into a C3'-endo conformation. This enhances hybridization to DNA and RNA, making LNAs useful in microarrays and potential therapeutics. Here, the LNA, L(CAAU), provides a simplified benchmark for testing the ability of molecular dynamics (MD) to approximate nucleic acid properties. LNA χ torsions and partial charges were parametrized to create AMBER parm99_LNA. The revisions were tested by comparing MD predictions with AMBER parm99 and parm99_LNA against a 200 ms NOESY NMR spectrum of L(CAAU). NMR indicates an A-Form equilibrium ensemble. In 3000 ns simulations starting with an A-form structure, parm99_LNA and parm99 provide 66% and 35% agreement, respectively, with NMR NOE volumes and (3)J-couplings. In simulations of L(CAAU) starting with all χ torsions in a syn conformation, only parm99_LNA is able to repair the structure. This implies methods for parametrizing force fields for nucleic acid mimics can reasonably approximate key interactions and that parm99_LNA will improve reliability of MD studies for systems with LNA. A method for approximating χ population distribution on the basis of base to sugar NOEs is also introduced.

Growth performance and quality properties of meat from broiler chickens reared under different monochromatic light sources

1. Quality properties of breast and drumstick muscles of female broiler chickens reared under different light sources were evaluated using a total of 480 female chicks (Ross 308). 2. After hatch, the chicks were randomly divided into three lighting treatment groups: control (daylight; mini incandescent lamps), G-B lighting group (green light monochromatic (MC) lamps for first 3 weeks switching to blue MC lamps for remaining 3 weeks), G-GB mixed lighting group (Green MC light for first 3 weeks, switching to Green MC light + Blue MC light mixture for remaining 3 weeks). 3. Feed consumption, body weight and total muscle weight values of the muscles from G-B and G-GB mix lighting groups were significantly higher than those of incandescent (control) lighting groups. 4. The breast and drumstick muscles from control lighting groups had a lower pH and water-holding capacity, but higher cooking loss values than those from G-B and G-GB mix lighting groups. The muscles from G-GB mix lighting groups exhibited the softer structure than those from other lighting groups. 5. We suggest that G-B or G-GB mix lighting during the rearing period of female broilers would improve body and muscle growth and then meat quality properties.

A Small Molecule that Binds an RNA Repeat Expansion Stimulates Its Decay via the Exosome Complex

Many diseases are caused by toxic RNA repeats. Herein, we designed a lead small molecule that binds the structure of the r(CUG) repeat expansion [r(CUG)^exp] that causes myotonic dystrophy type 1 (DM1) and Fuchs endothelial corneal dystrophy (FECD) and rescues disease biology in patient-derived cells and in vivo. Interestingly, the compound's downstream effects are different in the two diseases, owing to the location of the repeat expansion. In DM1, r(CUG)^exp is harbored in the 3' untranslated region, and the compound has no effect on the mRNA's abundance. In FECD, however, r(CUG)^exp is located in an intron, and the small molecule facilitates excision of the intron, which is then degraded by the RNA exosome complex. Thus, structure-specific, RNA-targeting small molecules can act disease specifically to affect biology, either by disabling the gain-of-function mechanism (DM1) or by stimulating quality control pathways to rid a disease-affected cell of a toxic RNA (FECD).

Improving Computational Predictions of Single-Stranded RNA Tetramers with Revised α/γ Torsional Parameters for the Amber Force Field

With current advancements in RNA based therapeutics, it is becoming crucial to utilize theoretical and computational methods to describe properly the physical properties of RNA molecules. NMR and X-ray crystallography are two powerful techniques for investigating structural properties. However, if the RNA molecules are complex or dynamic, these methods might not be adequate. For computational approaches, the quality of the force field will determine accuracy of our predictions. In this contribution, we revise the α/γ torsional parameters of RNA for amber force field using a model system representing an RNA dimer backbone. Combined with revised χ torsional parameters, previously shown to improve computational predictions, we benchmarked the revised force field on five single-stranded RNA (ssRNA) tetramers, three RNA dodecamer duplexes, and an RNA hairpin. A total of 60 μs of molecular dynamics (MD) simulations were run. We also employ the discrete path sampling (DPS) approach to compare the predictions for the revised amber force field with those for amber10. Our results indicate that the unphysical states observed with amber10 in ssRNA MD simulations are suppressed for the revised amber force field. In line with NMR experimental observations, incorporation of the revised α/γ and χ torsional parameters leads to A-form-like conformational states as the most favorable ssRNA tetramer conformations. Furthermore, the revised force field maintains the A-form geometry in regular RNA duplexes. Our revised amber force field for RNA should therefore improve structural and thermodynamic predictions for challenging RNA systems.