Genomic Profile of Fatigued Men Receiving Localized Radiation Therapy.

PURPOSE: The purpose of this study was to explore gene expression changes in fatigued men with nonmetastatic prostate cancer receiving localized external beam radiation therapy (EBRT). METHODS: Fatigue was measured in 40 men with prostate cancer (20 receiving EBRT and 20 controls on active surveillance) using the Functional Assessment of Cancer Therapy-Fatigue (FACT-F). EBRT subjects were followed from baseline to midpoint and end point of EBRT, while controls were seen at one time point. EBRT subjects were categorized into high- and low-fatigue groups based on change in FACT-F scores from baseline to EBRT completion. Full genome microarray was performed from peripheral leukocyte RNA to determine gene expression changes related to fatigue phenotypes. Real-time polymerase chain reaction and enzyme-linked immunosorbent assay confirmed the most differentially expressed gene in the microarray experiment. RESULTS: At baseline, mean FACT-F scores were not different between EBRT subjects (44.3 ± 7.16) and controls (46.7 ± 4.32, p = .24). Fatigue scores of EBRT subjects decreased at treatment midpoint (38.6 ± 9.17, p = .01) and completion (37.6 ± 9.9, p = .06), indicating worsening fatigue. Differential expression of 42 genes was observed between fatigue groups when EBRT time points were controlled. Membrane-spanning four domains, subfamily A, member (MS4A1) was the most differentially expressed gene and was associated with fatigue at treatment end point (r = -.46, p = .04). CONCLUSION: Fatigue intensification was associated with MS4A1 downregulation, suggesting that fatigue during EBRT may be related to impairment in B-cell immune response. The 42 differentially expressed fatigue-related genes are associated with glutathione biosynthesis, γ-glutamyl cycle, and antigen presentation pathways.

A Diagnosis of Insomnia Is Associated With Differential Expression of Sleep-Regulating Genes in Military Personnel.

Sleep disturbance is a common and disturbing symptom in military personnel, with many individuals progressing to the development of insomnia, which is characterized by increased arousals, wakefulness after sleep onset, and distorted sleep architecture. The molecular mechanisms underlying insomnia remain elusive, limiting future therapeutic development to address this critical issue. We examined whole gene expression profiles associated with insomnia. We compared subjects with insomnia (n = 25) to controls (n = 13) without insomnia using microarray gene expression profiles obtained from peripheral samples of whole blood obtained from military personnel. Compared to controls, participants with insomnia had differential expression of 44 transcripts from 43 identified genes. Among the identified genes, urotensin 2 was downregulated by more than 6 times in insomnia participants, and the fold-change remained significant after controlling for depression, posttraumatic stress disorder, and medication use. Urotensin 2 is involved in regulation of orexin A and B activity and rapid eye movement during sleep. These findings suggest that differential expression of these sleep-regulating genes contributes to symptoms of insomnia and, specifically, that switching between rapid eye movement and nonrapid eye movement sleep stages underlies insomnia symptoms. Future work to identify therapeutic agents that are able to regulate these pathways may provide novel treatments for insomnia.

Sleep quality, BDNF genotype and gene expression in individuals with chronic abdominal pain.

BACKGROUND: Sleep quality and genetics may contribute to the etiology of gastrointestinal (GI) symptoms. Individuals with impaired sleep often have a number of associated symptoms including chronic abdominal pain (CAP). The current study examined the interactions of brain-derived neurotrophic factor (BDNF) genotype with sleep quality in persons with CAP and healthy controls. In addition, associations among sleep quality, BDNF genotype, and gene expression were explored in the participants. METHODS: Data were collected on 59 participants (46% male, 61% White, 26.9 ± 6.6 years; CAP (n=19) and healthy controls (n=40)). Participants with CAP reported poorer sleep quality compared to healthy controls. BDNF genotype, categorized as Val/Val homozygotes versus the Met carriers. RESULTS: Microarray analysis found twenty-four differentially expressed genes by a two-fold magnitude in participants with poor sleep quality compared to good sleep quality, and seven differentially expressed genes comparing CAP to healthy control. Three specific genes in the pain group overlap with sleep quality, including insulin-like growth factor 1 (IGF1), spermatogenesis associated serine-rich 2-like (SPATS2L), and immunoglobulin heavy constant gamma 1 or mu (IGHG1/// IGHM). BDNF was shown to have an interaction effect with GI and sleep symptoms. CONCLUSIONS: Participants with CAP reported poor sleep quality compared to healthy controls. The role of the BDNF Met allele on differential gene expression was not distinct as main factor, but impacted interactions with sleep quality and CAP. Down-regulation of IGF1, SPATS2L, and IGHG1 expression may be related to the etiology of poor sleep quality and CAP. TRIAL REGISTRATION: Clinicaltrial.gov # NCT00824941.

Military personnel with chronic symptoms following blast traumatic brain injury have differential expression of neuronal recovery and epidermal growth factor receptor genes.

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.

Gene expression profiles of fatigued fibromyalgia patients with different categories of pain and catastrophizing: a preliminary report.

BACKGROUND: Fibromyalgia (FM) is a chronic condition characterized by diffused musculoskeletal pain and overwhelming fatigue. PURPOSE: To compare the gene expression profiles of fatigued FM women with different levels of pain and catastrophizing. METHODS: Nine women with FM enrolled in an active Medstar Research Institute protocol were included in the gene expression analyses of peripheral blood RNA via Affymetrix GeneChip Human Genome U133 Plus 2.0 array (Santa Clara, CA). Scores from Brief Pain Inventory, Pain Catastrophizing Scale, and Multidimensional Fatigue Inventory categorized the nine participants into pain (high, n = 3; low, n = 6) and catastrophizing groups (high, n = 5; low, n = 4). DISCUSSION: Differential expression of 107 genes between the high and low pain groups and 139 genes between the high and low catastrophizing groups (over 2.0-fold change, p < .05) were observed. Network analyses showed interferon signaling and interferon regulatory activation factor pathways distinguished between the pain groups whereas dendritic cell maturation delineated between the catastrophizing groups. CONCLUSION: Findings provide preliminary evidence that specific physiological pathways may possibly delineate pain and catastrophizing mechanisms. Further investigation via the use of a larger and more homogenous sample is warranted.

Mechanism of methanol oxidation by quinoprotein methanol dehydrogenase.

At neutral pH, oxidation of CH(3)OH --> CH(2)O by an o-quinone requires general-base catalysis and the reaction is endothermic. The active-site -CO(2)(-) groups of Glu-171 and Asp-297 (Glu-171-CO(2)(-) and Asp-297-CO(2)(-)) have been considered as the required general base catalysts in the bacterial o-quinoprotein methanol dehydrogenase (MDH) reaction. Based on quantum mechanics/molecular mechanics (QM/MM) calculations, the free energy for MeOH reduction of o-PQQ when MeOH is hydrogen bonded to Glu-171-CO(2)(-) and the crystal water (Wat1) is hydrogen bonded to Asp-297-CO(2)(-) is DeltaG++ = 11.7 kcal/mol, which is comparable with the experimental value of 8.5 kcal/mol. The calculated DeltaG++ when MeOH is hydrogen bonded to Asp-297-CO(2)(-) is >50 kcal/mol. The Asp-297-CO(2)(-)...Wat1 complex is very stable. Molecular dynamics (MD) simulations on MDH.PQQ.Wat1 complex in TIP3P water for 5 ns does not result in interchange of Asp-297-CO(2)(-) bound Wat1 for a solvent water. Starting with Wat1 removed and MeOH hydrogen bonded to Asp-297-CO(2)(-), we find that MeOH returns to be hydrogen bonded to Glu-171-CO(2)(-) and Asp-297-CO(2)(-) coordinates to Ca(2+) during 3 ns simulation. The Asp-297-CO(2)(-)...Wat1 of reactant complex does play a crucial role in catalysis. By QM/MM calculation DeltaG++ = 1.1 kcal/mol for Asp-297-CO(2)(-) general-base catalysis of Wat1 hydration of the immediate CH(2)==O product --> CH(2)(OH)(2). By this means, the endothermic oxidation-reduction reaction is pulled such that the overall conversion of MeOH to CH(2)(OH)(2) is exothermic.

An investigation of the ionic and solvation patterns of dsDNG versus dsDNA by use of molecular dynamics simulations.

Molecular dynamics has been employed to analyze the counterion and water atmospheres around the deoxynucleic guanidine (DNG) duplex G(12)-C(12). These features are compared to those of DNA G(12)-C(12). The chloride counterions of cationic DNG demonstrate fewer penetrations and only fleeting residence times in the minor groove, as opposed to the multi-nanosecond visits seen by sodium ions in DNA minor grooves. The 10 ns simulations also show the differences in hydration patterns between the DNG and DNA duplexes.

Comparison of positively charged DNG with DNA duplexes: a computational approach.

Molecular dynamics is used to investigate the structural properties of the cationic DNA analogue deoxynucleic guanidine (DNG), in which a guanidinium group replaces the phosphate moiety of DNA. This study examines the DNG duplex dodecamers d(Ag)(12).d(Tg)(12) and d(Gg)(12).d(Cg)(12), as well as their DNA counterparts. Watson-Crick base-pairing is maintained in the solvated DNG duplex models during the 5ns simulations. The idealized DNG dodecamers assume many parameters characteristic of the corresponding native DNA, assuming B-DNA conformations. Several helical parameters are rather unique to DNG, including buckle, slide, inclination, propeller, and X-displacement. Fewer transitions in backbone torsions occur in the DNG duplexes compared to those of the DNA, which may result from the greater rigidity of the sp(2) hybridized guanidinium group verses the flexible sp(3) phosphate group. The DNG helices have exceptionally shallow major grooves and very deep minor grooves. The major and minor groove widths of DNG are narrower than those of the respective DNA counterparts.

Mechanisms of ammonia activation and ammonium ion inhibition of quinoprotein methanol dehydrogenase: a computational approach.

The mechanism of methanol oxidation by quinoprotein methanol dehydrogenase (MDH.PQQ) in combination with methanol (MDH.PQQ.methanol) involves Glu-171--CO2(-) general base removal of the hydroxyl proton of methanol in concert with hydride equivalent transfer to the >C5=O quinone carbon of pyrroloquinoline quinone (PQQ) and rearrangement to hydroquinone (PQQH2) with release of formaldehyde. Molecular dynamics (MD) studies of the structures of MDH.PQQ.methanol in the presence of activator NH3 and inhibitor NH4(+) have been carried out. In the MD structure of MDH.PQQ.methanol.NH3, the hydrated NH3 resides at a distance of approximately 24 A away from methanol and the ortho-quinone portion of PQQ. As such, influence of NH3 on the oxidation reaction is not probable. We find that NH4(+) competes with the substrate by hydrogen-bonding to Glu-171CO2(-) such that the MDH.PQQ.methanol.NH4(+) complex is not reactive. Ammonia readily forms imines with quinone. Imines are present in solution as neutral (>C5=NH) and protonated (>C5=NH2(+)) species. MD simulations establish that the >C5=NH2(+) derivative of MDH.PQQ(NH2(+).methanol structure is unreactive because of the nonproductive means of methanol binding. The structure obtained by the MD simulations with the neutral >C5=NH imine of MDH.PQQ(NH).methanol structure is similar to the reactive MDH.PQQ.methanol complex. This active site geometry allows for catalysis of hydride equivalent transfer to the >C5=NH of PQQ(NH) by concerted Glu-171CO(2)(-) general-base removal of the H-OCH3 proton and Arg-324H+ general-acid proton transfer to the imine nitrogen. Enzyme-bound

Determination of enzyme mechanisms by molecular dynamics: studies on quinoproteins, methanol dehydrogenase, and soluble glucose dehydrogenase.

Molecular dynamics (MD) simulations have been carried out to study the enzymatic mechanisms of quinoproteins, methanol dehydrogenase (MDH), and soluble glucose dehydrogenase (sGDH). The mechanisms of reduction of the orthoquinone cofactor (PQQ) of MDH and sGDH involve concerted base-catalyzed proton abstraction from the hydroxyl moiety of methanol or from the 1-hydroxyl of glucose, and hydride equivalent transfer from the substrate to the quinone carbonyl carbon C5 of PQQ. The products of methanol and glucose oxidation are formaldehyde and glucolactone, respectively. The immediate product of PQQ reduction, PQQH- [-HC5(O-)-C4(=O)-] and PQQH [-HC5(OH)-C4(=O)-] converts to the hydroquinone PQQH2 [-C5(OH)=C4(OH)-]. The main focus is on MD structures of MDH * PQQ * methanol, MDH * PQQH-, MDH * PQQH, sGDH * PQQ * glucose, sGDH * PQQH- (glucolactone, and sGDH * PQQH. The reaction PQQ-->PQQH- occurs with Glu 171-CO2- and His 144-Im as the base species in MDH and sGDH, respectively. The general-base-catalyzed hydroxyl proton abstraction from substrate concerted with hydride transfer to the C5 of PQQ is assisted by hydrogen-bonding to the C5=O by Wat1 and Arg 324 in MDH and by Wat89 and Arg 228 in sGDH. Asp 297-COOH would act as a proton donor for the reaction PQQH(-)-->PQQH, if formed by transfer of the proton from Glu 171-COOH to Asp 297-CO2- in MDH. For PQQH-->PQQH2, migration of H5 to the C4 oxygen may be assisted by a weak base like water (either by crystal water Wat97 or bulk solvent, hydrogen-bonded to Glu 171-CO2- in MDH and by Wat89 in sGDH).

Mechanism of glucose oxidation by quinoprotein soluble glucose dehydrogenase: insights from molecular dynamics studies.

We have generated 3 ns molecular dynamic (MD) simulations, in aqueous solution, of the bacterial soluble glucose dehydrogenase enzyme.PQQ.glucose complex and intermediates formed in PQQ reduction. In the MD structure of enzyme.PQQ.glucose complex the imidazole of His144 is hydrogen bonded to the hydroxyl hydrogen of H[bond]OC1(H) of glucose. The tightly hydrogen-bonded triad Asp163-His144-glucose (2.70 and 2.91 A) is involved in proton abstraction from glucose concerted with the hydride transfer from the C1[bond]H of glucose to the >C5[double bond]O quinone carbon of PQQ. The reaction is assisted by Arg228 hydrogen bonding to the carbonyl oxygen of >C5[double bond]O. The rearrangement of [bond](H)C5(O-)[bond]C4([double bond]O)[bond] of II to [bond]C5(OH)[double bond]C4(OH)[bond] of PQQH(2) hydroquinone is assisted by general acid protonatation of the >C4[double bond]O oxygen by protonated His144 and hydrogen bonds of Arg228 to the oxyanion O5. The continuous hydrogen bonding of the amide side chain of Asn229 to >C4[double bond]O4 oxygen and that of the O5 oxygen of the cofactor to Wat89 is observed throughout the entire reaction.

Conformations and dynamics of Ets-1 ETS domain-DNA complexes.

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 1500 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 1500 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.

In silico studies of the mechanism of methanol oxidation by quinoprotein methanol dehydrogenase.

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.

Sequence specific DNA binding of Ets-1 transcription factor: molecular dynamics study on the Ets domain--DNA complexes.

Molecular dynamics (MD) simulations for Ets-1 ETS domain-DNA complexes were performed to investigate the mechanism of sequence-specific recognition of the GGAA DNA core by the ETS domain. Employing the crystal structure of the Ets-1 ETS domain-DNA complex as a starting structure we carried out MD simulations of: (i). the complex between Ets-1 ETS domain and a 14 base-pair DNA containing GGAA core sequence (ETS-GGAA); (ii). the complex between the ETS domain and a DNA having single base-pair mutation, GGAG sequence (ETS-GGAG); and (iii). the 14 base-pair DNA alone (GGAA). Comparative analyses of the MD structures of ETS-GGAA and ETS-GGAG reveal that the DNA bending angles and the ETS domain-DNA phosphate interactions are similar in these complexes. These results support that the GGAA core sequence is distinguished from the mutated GGAG sequence by a direct readout mechanism in the Ets-1 ETS domain-DNA complex. Further analyses of the direct contacts in the interface between the helix-3 region of Ets-1 and the major groove of the core DNA sequence clearly show that the highly conserved arginine residues, Arg391 and Arg394, play a critical role in binding to the GGAA core sequence. These arginine residues make bidentate contacts with the nucleobases of GG dinucleotides in GGAA core sequence. In ETS-GGAA, the hydroxyl group of Tyr395 is hydrogen bonded to N7 nitrogen of A(3) (the third adenosine in the GGAA core), while the hydroxyl group makes a contact with N4 nitrogen of C(4') (the complementary nucleotide of the fourth guanosine G(4) in the GGAG sequence) in the ETS-GGAG complex. We have found that this difference in behavior of Tyr395 results in the relatively large motion of helix-3 in the ETS-GGAG complex, causing the collapse of bidentate contacts between Arg391/Arg394 and the GG dinucleotides in the GGAG sequence.

DNA polymorphism: a comparison of force fields for nucleic acids.

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.

Protein engineering of nitrile hydratase activity of papain: molecular dynamics study of a mutant and wild-type enzyme.

The mechanism of hydrolysis of the nitrile (N-acetyl-phenylalanyl-2-amino-propionitrile, I) catalyzed by Gln19Glu mutant of papain has been studied by nanosecond molecular dynamics (MD) simulations. MD simulations of the complex of mutant enzyme with I and of mutant enzyme covalently attached to both neutral (II) and protonated (III) thioimidate intermediates were performed. An MD simulation with the wild-type enzyme.I complex was undertaken as a reference. The ion pair between protonated His159 and thiolate of Cys25 is coplanar, and the hydrogen bonding interaction S(-)(25).HD1-ND1(159) is observed throughout MD simulation of the mutant enzyme.I complex. Such a sustained hydrogen bond is absent in nitrile-bound wild-type papain due to the flexibility of the imidazole ring of His159. The nature of the residue at position 19 plays a critical role in the hydrolysis of the covalent thioimidate intermediate. When position 19 represents Glu, the imidazolium ion of His159-ND1(+).Cys25-S(-) ion pair is distant, on average, from the nitrile nitrogen of substrate I. Near attack conformers (NACs) have been identified in which His159-ImH(+) is positioned to initiate a general acid-catalyzed addition of Cys-S(-) to nitrile. Though Glu19-CO(2)H is distant from nitrile nitrogen in the mutant.I structure, MD simulations of the mutant.II covalent adduct finds Glu19-CO(2)H hydrogen bonded to the thioimide nitrogen of II. This hydrogen bonded species is much less stable than the hydrogen bonded Glu19-CO(2)(-) with mutant-bound protonated thioimidate (III). This observation supports Glu19-CO(2)H general acid catalysis of the formation of mutant.III. This is the commitment step in the Gln19Glu mutant catalysis of nitrile hydrolysis.