The improvements of the force fields and the more accurate treatment of long-range interactions are providing more reliable molecular dynamics simulations of nucleic acids. The abilities of certain nucleic acid force fields to represent the structural and conformational properties of nucleic acids in solution are compared. The force fields are AMBER 4.1, BMS, CHARMM22, and CHARMM27; the comparison of the latter two is the primary focus of this paper. The performance of each force field is evaluated first on its ability to reproduce the B-DNA decamer d(CGATTAATCG)(2) in solution with simulations in which the long-range electrostatics were treated by the particle mesh Ewald method; the crystal structure determined by Quintana et al. (1992) is used as the starting point for all simulations. A detailed analysis of the structural and solvation properties shows how well the different force fields can reproduce sequence-specific features. The results are compared with data from experimental and previous theoretical studies.

The mechanism of bacterial methanol dehydrogenase involves hydride equivalent transfer from substrate to the ortho-quinone PQQ to provide a C5-reduced intermediate that subsequently rearranges to the hydroquinone PQQH(2). We have studied the PQQ reduction by molecular dynamic (MD) simulations in aqueous solution. Among the five simulated structures, either Asp297 or Glu171 or both are ionized. Reasonable structures are obtained only when both carboxyl groups are ionized. This is not unexpected since the kinetic pH optimum is 9.0. In the structure of the enzyme.PQQ.HOCH(3) complex, the hydrogen bonded Glu171-CO(2)(-).H-OCH(3) is in a position to act as a general base catalyst for hydride equivalent transfer to C5 of PQQ. We thus suggest that Glu171 plays the role of general base catalyst in PQQ reduction rather than Asp297 as previously suggested. The reduction is assisted by Arg324, which hydrogen bonds to the ortho-quinone moiety of PQQ. The rearrangement of the C5-reduced intermediate to provide hydroquinone PQQH(2) is also assisted by proton abstraction by Glu171-CO(2)(-) and the continuous hydrogen bonding of Arg324 throughout the entire reaction. These features as well as the mapping of the channel for substrate and water into the active site entrance are the observations of major importance.

Objective: Approximately one-quarter of military personnel who deployed to combat stations sustained one or more blast-related, closed-head injuries. Blast injuries result from the detonation of an explosive device. The mechanisms associated with blast exposure that give rise to traumatic brain injury (TBI), and place military personnel at high risk for chronic symptoms of post-concussive disorder (PCD), post-traumatic stress disorder (PTSD), and depression are not elucidated.Methods: To investigate the mechanisms of persistent blast-related symptoms, we examined expression profiles of transcripts across the genome to determine the role of gene activity in chronic symptoms following blast-TBI. Active duty military personnel with (1) a medical record of a blast-TBI that occurred during deployment (n = 19) were compared to control participants without TBI (n = 17). Controls were matched to cases on demographic factors including age, gender, and race, and also in diagnoses of sleep disturbance, and symptoms of PTSD and depression. Due to the high number of PCD symptoms in the TBI+ group, we did not match on this variable. Using expression profiles of transcripts in microarray platform in peripheral samples of whole blood, significantly differentially expressed gene lists were generated. Statistical threshold is based on criteria of 1.5 magnitude fold-change (up or down) and p-values with multiple test correction (false discovery rate <0.05).Results: There were 34 transcripts in 29 genes that were differentially regulated in blast-TBI participants compared to controls. Up-regulated genes included epithelial cell transforming sequence and zinc finger proteins, which are necessary for astrocyte differentiation following injury. Tensin-1, which has been implicated in neuronal recovery in pre-clinical TBI models, was down-regulated in blast-TBI participants. Protein ubiquitination genes, such as epidermal growth factor receptor, were also down-regulated and identified as the central regulators in the gene network determined by interaction pathway analysis.Conclusion: In this study, we identified a gene-expression pathway of delayed neuronal recovery in military personnel a blast-TBI and chronic symptoms. Future work is needed to determine if therapeutic agents that regulate these pathways may provide novel treatments for chronic blast-TBI-related symptoms.

Molecular dynamics studies have been performed for 3.5 ns on the ETS domain of Ets-1 transcription factor bound to the 14-bp DNA, d(AGTGCCGGAAATGT), comprising the core sequence of high-affinity (GGAA), ETS-GGAA. In like manner, molecular dynamics simulations have been carried out for 3.9 ns on the mutant low-affinity core sequence, GGAG (ETS-GGAG). Analyses of the DNA backbone of ETS-GGAG show conformational interconversions from BI to BII substates. Also, crank shaft motions are noticed at the mutated nucleotide base pair step after 1,500 ps of dynamics. The corresponding nucleotide of ETS-GGAA is characteristic of a BI conformation and no crank shaft motions are observed. The single mutation of ETS-GGAA to ETS-GGAG also results in variations of helical parameters and solvent-accessible surface area around the major and minor grooves of the DNA. The presence of water contacts during the entire simulation proximal to the fourth base pair step of core DNA sequence is a characteristic feature of ETS-GGAA. Such waters are more mobile in ETS-GGAG at 100 ps and distant after 1,500 ps. Anticorrelated motions between certain amino acids of Ets-1 protein are predominant in ETS-GGAA but less so or absent in the mutant. These motions are reflected in the flexibility of amino acid residues of the protein backbone. We consider that these conformational features and water contacts are involved in stabilizing the hydrogen bond interactions between helix-3 residues of Ets-1 and DNA during the transcription process. The Ets protein family of transcription factors includes species interacting with various genes that code for transcriptional activators and inhibitors involved in cell proliferation and differentiation (1, 2). The regulation of the initiation of gene transcription arises from the combined activity of different transcriptional regulators (2, 3). Ets family members found in species from invertebrates to humans share a conserved sequence of 85 amino acids, named the ETS domain. The ETS domain folds into a winged helix-turn-helix motif and binds to a consensus DNA sequence centered on the core GGAA motif, named the Ets-binding site. The sequences flanking this core motif (in the major groove) are variable and characterize the specificity of binding of the Ets transcription factor. Ets proteins have also been implicated in several types of cancer and other human diseases (4). Detailed conformational preferences that influence the sequence specificity of Ets proteins are essential for the design of anticancer drugs.The high-affinity DNA contains the GGAA core sequence, ETS-GGAA (Fig. 1). The low-affinity DNA is the single-base-pair mutant, ETS-GGAG. Recently we reported molecular dynamics (MD) studies (5) dealing with the binding of the ETS domain of Ets-1 protein to the high-and low-affinity 14-bp DNA structures. We have observed that the most conserved residues Arg-391, Arg-394, along with Tyr-395 of Ets-1, jointly contribute to recognize the GGAA or GGAG core DNA sequences (5). The differential hydrogen bond interact...