JGromacs: A Java Package for Analyzing Protein Simulations

UNLABELLED: In this paper, we introduce JGromacs, a Java API (Application Programming Interface) that facilitates the development of cross-platform data analysis applications for Molecular Dynamics (MD) simulations. The API supports parsing and writing file formats applied by GROMACS (GROningen MAchine for Chemical Simulations), one of the most widely used MD simulation packages. JGromacs builds on the strengths of object-oriented programming in Java by providing a multilevel object-oriented representation of simulation data to integrate and interconvert sequence, structure, and dynamics information. The easy-to-learn, easy-to-use, and easy-to-extend framework is intended to simplify and accelerate the implementation and development of complex data analysis algorithms. Furthermore, a basic analysis toolkit is included in the package. The programmer is also provided with simple tools (e.g., XML-based configuration) to create applications with a user interface resembling the command-line interface of GROMACS applications. AVAILABILITY: JGromacs and detailed documentation is freely available from http://sbcb.bioch.ox.ac.uk/jgromacs under a GPLv3 license

Subcellular localization and formation of huntingtin aggregates correlates with symptom onset and progression in a Huntington's disease model

Huntington's disease is caused by the expansion of a CAG repeat within exon 1 of the HTT gene, which is unstable, leading to further expansion, the extent of which is brain region and peripheral tissue specific. The identification of DNA repair genes as genetic modifiers of Huntington's disease, that were known to abrogate somatic instability in Huntington's disease mouse models, demonstrated that somatic CAG expansion is central to disease pathogenesis, and that the CAG repeat threshold for pathogenesis in specific brain cells might not be known. We have previously shown that the HTT gene is incompletely spliced generating a small transcript that encodes the highly pathogenic exon 1 HTT protein. The longer the CAG repeat, the more of this toxic fragment is generated, providing a pathogenic consequence for somatic expansion. Here, we have used the R6/2 mouse model to investigate the molecular and behavioural consequences of expressing exon 1 HTT with 90 CAGs, a mutation that causes juvenile Huntington's disease, compared to R6/2 mice carrying ∼200 CAGs, a repeat expansion of a size rarely found in Huntington's disease patient's blood, but which has been detected in post-mortem brains as a consequence of somatic CAG repeat expansion. We show that nuclear aggregation occurred earlier in R6/2(CAG)(90) mice and that this correlated with the onset of transcriptional dysregulation. Whereas in R6/2(CAG)(200) mice, cytoplasmic aggregates accumulated rapidly and closely tracked with the progression of behavioural phenotypes and with end-stage disease. We find that aggregate species formed in the R6/2(CAG)(90) brains have different properties to those in the R6/2(CAG)(200) mice. Within the nucleus, they retain a diffuse punctate appearance throughout the course of the disease, can be partially solubilized by detergents and have a greater seeding potential in young mice. In contrast, aggregates from R6/2(CAG)(200) brains polymerize into larger structures that appear as inclusion bodies. These data emphasize that a subcellular analysis, using multiple complementary approaches, must be undertaken in order to draw any conclusions about the relationship between HTT aggregation and the onset and progression of disease phenotypes

Proteomic Profile of Reversible Protein Oxidation Using PROP, Purification of Reversibly Oxidized Proteins

Signal transduction pathways that are modulated by thiol oxidation events are beginning to be uncovered, but these discoveries are limited by the availability of relatively few analytical methods to examine protein oxidation compared to other signaling events such as protein phosphorylation. We report here the coupling of PROP, a method to purify reversibly oxidized proteins, with the proteomic identification of the purified mixture using mass spectrometry. A gene ontology (GO), KEGG enrichment and Wikipathways analysis of the identified proteins indicated a significant enrichment in proteins associated with both translation and mRNA splicing. This methodology also enabled the identification of some of the specific cysteine residue targets within identified proteins that are reversibly oxidized by hydrogen peroxide treatment of intact cells. From these identifications, we determined a potential consensus sequence motif associated with oxidized cysteine residues. Furthermore, because we identified proteins and specific sites of oxidation from both abundant proteins and from far less abundant signaling proteins (e.g. hepatoma derived growth factor, prostaglandin E synthase 3), the results suggest that the PROP procedure was efficient. Thus, this PROP-proteomics methodology offers a sensitive means to identify biologically relevant redox signaling events that occur within intact cells

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