Targeted DNA oxidation and trajectory of radical DNA using DFT based QM/MM dynamics.

Molecular insight into electronic rearrangements and structural trajectories arising from oxidative damages to DNA backbone is of crucial importance in understanding the effect of ionizing radiation, developing DNA biosensors and designing effective DNA cleaving molecules. Employing a Density Functional Theory based multi-scale Quantum-Mechanical-Molecular-Mechanical (QM/MM) simulation and a suitable partitioning of the Hamiltonian on solvated nucleotide, and single-, and double-stranded DNA, we mimic hydrogen transfer reactions from the backbone by OH radicals and report structural trajectories arising from on-the-fly electronic charge- and spin-density redistribution in these three different structural topologies of DNA. Trajectories reveal that H4' abstraction can disrupt the deoxyribose moiety through the formation of C4'=O4' ketone and a π-bond with base at C1'-N9 in a nucleotide versus only partial ketone formation in single- and double-stranded DNA, where the orientation of the base is topologically restrained. However, H5' abstraction can lead DNA cleavage at 5' end through the formation of C5'=O5' ketone and breakage of P-O5' bond. Results demonstrate that structural damages from oxidative reactions are restrained by base stacking and base-pair hydrogen bonding. The methodology can be suitably used to study targeted DNA and RNA damages from radicals and radiomimetic drugs to design DNA cleaving molecules for chemotherapy.

Deciphering the Complex Mineralogy of River Sand Deposits through Clustering and Quantification of Hyperspectral X-Ray Maps.

Alluvial mineral sands rank among the most complex subjects for mineral characterization due to the diverse range of minerals present in the sediments, which may collectively contain a daunting number of elements (>20) in major or minor concentrations (>1 wt%). To comprehensively characterize the phase abundance and chemistry of these complex mineral specimens, a method was developed using hyperspectral x-ray and cathodoluminescence mapping in an electron probe microanalyser (EPMA), coupled with automated cluster analysis and quantitative analysis of clustered x-ray spectra. This method proved successful in identifying and quantifying over 40 phases from mineral sand specimens, including unexpected phases with low modal abundance (<0.1%). The standard-based quantification method measured compositions in agreement with expected stoichiometry, with elemental detection limits in the range of <10­1,000 ppm, depending on phase abundance, and proved reliable even for challenging mineral species, such as the multi-rare earth element (REE) bearing mineral xenotime [(Y,REE)PO4] for which 24 elements were analyzed, including 12 overlapped REEs. The mineral identification procedure was also capable of characterizing mineral groups that exhibit significant compositional variability due to the substitution of multiple elements, such as garnets (Mg, Ca, Fe, Mn, Cr), pyroxenes (Mg, Ca, Fe), and amphiboles (Na, Mg, Ca, Fe, Al).

Salmonella FraE, an Asparaginase Homolog, Contributes to Fructose-Asparagine but Not Asparagine Utilization.

Salmonella enterica can utilize fructose-asparagine (F-Asn) as a source of carbon and nitrogen. This capability has been attributed to five genes in the fra locus. Previously, we determined that mutations in fraB (deglycase), fraD (kinase), or fraA (transporter) eliminated the ability of Salmonella to grow on F-Asn, while a mutation in fraE allowed partial growth. We hypothesized that FraE, a putative periplasmic fructose-asparaginase, converts F-Asn to NH4 + and fructose-aspartate (F-Asp). FraA could then transport F-Asp into the cytoplasm for subsequent catabolism. Here, we report that growth of the fraE mutant on F-Asn is caused by a partially redundant activity provided by AnsB, a periplasmic asparaginase. Indeed, a fraE ansB double mutant is unable to grow on F-Asn. Moreover, biochemical assays using periplasmic extracts of mutants that express only FraE or AnsB confirmed that each of these enzymes converts F-Asn to F-Asp and NH4 + However, FraE does not contribute to growth on asparagine. We tested and confirmed the hypothesis that a fraE ansB mutant can grow on F-Asp, while mutants lacking fraA, fraD, or fraB cannot. This finding provides strong evidence that FraA transports F-Asp but not F-Asn from the periplasm to the cytoplasm. Previously, we determined that F-Asn is toxic to a fraB mutant due to the accumulation of the FraB substrate, 6-phosphofructose-aspartate (6-P-F-Asp). Here, we found that, as expected, a fraB mutant is also inhibited by F-Asp. Collectively, these findings contribute to a better understanding of F-Asn utilization by Salmonella IMPORTANCE Salmonella is able to utilize fructose-asparagine (F-Asn) as a nutrient. We recently reported that the disruption of a deglycase enzyme in the F-Asn utilization pathway inhibits the growth of Salmonella in mice and recognized this pathway as a novel and specific drug target. Here, we characterize the first step in the pathway wherein FraE hydrolyzes F-Asn to release NH4 + and F-Asp in the periplasm of the cell. A fraE mutant continues to grow slowly on F-Asn due to asparaginase activity encoded by ansB.

Characterization of a Salmonella sugar kinase essential for the utilization of fructose-asparagine.

Salmonella can utilize fructose-asparagine (F-Asn), a naturally occurring Amadori product, as its sole carbon and nitrogen source. Conversion of F-Asn to the common intermediates glucose-6-phosphate, aspartate, and ammonia was predicted to involve the sequential action of an asparaginase, a kinase, and a deglycase. Mutants lacking the deglycase are highly attenuated in mouse models of intestinal inflammation owing to the toxic build-up of the deglycase substrate. The limited distribution of this metabolic pathway in the animal gut microbiome raises the prospects for antibacterial discovery. We report the biochemical characterization of the kinase that was expected to transform fructose-aspartate to 6-phosphofructose-aspartate during F-Asn utilization. In addition to confirming its anticipated function, we determined through studies of fructose-aspartate analogues that this kinase exhibits a substrate-specificity with greater tolerance to changes to the amino acid (including the d-isomer of aspartate) than to the sugar.

Structural insights into Resveratrol's antagonist and partial agonist actions on estrogen receptor alpha.

BACKGROUND: Resveratrol, a naturally occurring stilbene, has been categorized as a phytoestrogen due to its ability to compete with natural estrogens for binding to estrogen receptor alpha (ERα) and modulate the biological responses exerted by the receptor. Biological effects of resveratrol (RES) on estrogen receptor alpha (ERα) remain highly controversial, since both estrogenic and anti-estrogenic properties were observed. RESULTS: Here, we provide insight into the structural basis of the agonist/antagonist effects of RES on ERα ligand binding domain (LBD). Using atomistic simulation, we found that RES bound ERα monomer in antagonist conformation, where Helix 12 moves away from the ligand pocket and orients into the co-activator binding groove of LBD, is more stable than RES bound ERα in agonist conformation, where Helix 12 lays over the ligand binding pocket. Upon dimerization, the agonistic conformation of RES-ERα dimer becomes more stable compared to the corresponding monomer but still remains less stable compared to the corresponding dimer in antagonist conformation. Interestingly, while the binding pocket and the binding contacts of RES to ERα are similar to those of pure agonist diethylstilbestrol (DES), the binding energy is much less and the hydrogen bonding contacts also differ providing clues for the partial agonistic character of RES on ERα. CONCLUSIONS: Our Molecular Dynamics simulation of RES-ERα structures with agonist and antagonist orientations of Helix 12 suggests RES action is more similar to Selective Estrogen Receptor Modulator (SERM) opening up the importance of cellular environment and active roles of co-regulator proteins in a given system. Our study reveals that potential co-activators must compete with the Helix 12 and displace it away from the activator binding groove to enhance the agonistic activity.

Simulation of multiphase systems utilizing independent force fields to control intraphase and interphase behavior.

Fixed-charge empirical force fields have been developed and widely used over the past three decades for all-atom molecular simulations. Most simulation programs providing these methods enable only one set of force field parameters to be used for the entire system. Whereas this is generally suitable for single-phase systems, the molecular environment at the interface between two phases may be sufficiently different from the individual phases to require a different set of parameters to be used to accurately represent the system. Recently published simulations of peptide adsorption to material surfaces using the CHARMM force field have clearly demonstrated this issue by revealing that calculated values of adsorption free energy substantially differ from experimental results. Whereas nonbonded parameters could be adjusted to correct this problem, this cannot be done without also altering the conformational behavior of the peptide in solution, for which CHARMM has been carefully tuned. We have developed a dual-force-field approach (Dual-FF) to address this problem and implemented it in the CHARMM simulation package. This Dual-FF method provides the capability to use two separate sets of nonbonded force field parameters within the same simulation: one set to represent intraphase interactions and a separate set to represent interphase interactions. Using this approach, we show that interfacial parameters can be adjusted to correct errors in peptide adsorption free energy without altering peptide conformational behavior in solution. This program thus provides the capability to enable both intraphase and interphase molecular behavior to be accurately and efficiently modeled in the same simulation.

Development of a high-throughput assay for rapid screening of butanologenic strains.

We report a Thermotoga hypogea (Th) alcohol dehydrogenase (ADH)-dependent spectrophotometric assay for quantifying the amount of butanol in growth media, an advance that will facilitate rapid high-throughput screening of hypo- and hyper-butanol-producing strains of solventogenic Clostridium species. While a colorimetric nitroblue tetrazolium chloride-based assay for quantitating butanol in acetone-butanol-ethanol (ABE) fermentation broth has been described previously, we determined that Saccharomyces cerevisiae (Sc) ADH used in this earlier study exhibits approximately 13-fold lower catalytic efficiency towards butanol than ethanol. Any Sc ADH-dependent assay for primary quantitation of butanol in an ethanol-butanol mixture is therefore subject to "ethanol interference". To circumvent this limitation and better facilitate identification of hyper-butanol-producing Clostridia, we searched the literature for native ADHs that preferentially utilize butanol over ethanol and identified Th ADH as a candidate. Indeed, recombinant Th ADH exhibited a 6-fold higher catalytic efficiency with butanol than ethanol, as measured using the reduction of NADP+ to NADPH that accompanies alcohol oxidation. Moreover, the assay sensitivity was not affected by the presence of acetone, acetic acid or butyric acid (typical ABE fermentation products). We broadened the utility of our assay by adapting it to a high-throughput microtiter plate-based format, and piloted it successfully in an ongoing metabolic engineering initiative.

Electronic structure, ionization potential, and electron affinity of the enzyme cofactor (6R)-5,6,7,8-tetrahydrobiopterin in the gas phase, solution, and protein environments.

(6R)-5,6,7,8-Tetrahydrobiopterin (BH(4)) is a key cofactor involved in the electron transfer to the P(450) heme of nitric oxide synthase. We calculated the electronic structure of the neutral, cationic, and anionic forms of BH(4) in the gas phase, in solution (both dielectric and explicit water), and in the protein environment using density functional theory (B3LYP/6-31+G(d,p)). Subsequently, we derived the ionization potential (IP) and electron affinity (EA) of the cofactor in these chemical environments. We found that the electronic structure of BH(4) is susceptible to the presence of an external electric field and that conformational changes in the structure of BH(4) alone do not affect its electronic structure significantly. In the gas phase, water, and protein environments neutral BH(4) is the most stable species, while in the dielectric environment the anion becomes the most stable species. The IP of BH(4) in the protein environment is about half of that in the gas phase, and its EA is about 5 times smaller than that in the gas phase. Our results indicate that changes in the external electric field created by moving charged amino acid residues around BH(4) may lead to configurations that have the BH(4) ion as stable as or more stable than the neutral form, thus facilitating the electron transfer.

Molecular insight into the differential anti-androgenic activity of resveratrol and its natural analogs: in silico approach to understand biological actions.

The androgen receptor (AR) is a therapeutic target for the treatment of prostate cancer. Androgen receptor reactivation during the androgen-independent stage of prostate cancer is mediated by numerous mechanisms including expression of AR mutants and splice variants that become non-responsive to conventional anti-androgenic agents. Resveratrol and its natural analogs exhibit varying degrees of anti-androgenic effects on tumor growth suppression in prostate cancer. However, the structural basis for the observed differential activity remains unknown. Here, anti-androgenic activities of resveratrol and its natural analogs, namely, pterostilbene, piceatannol and trimethoxy-resveratrol were studied in LNCaP cells expressing T877A mutant AR and atomistic simulations were employed to establish the structure activity relationship. Interestingly, essential hydrogen bonding contacts and the binding energies of resveratrol analogs with AR ligand binding domain (LBD), emerge as key differentiating factors for varying anti-androgenic action. Among all the analogs, pterostilbene exhibited strongest anti-androgenic activity and its binding energy and hydrogen bonding interactions pattern closely resembled pure anti-androgen, flutamide. Principal component analysis of our simulation studies revealed that androgenic compounds bind more strongly to AR LBD compared to anti-androgenic compounds and provide conformational stabilization of the receptor in essential subspace. The present study provides critical insight into the structure-activity relationship of the anti-androgenic action of resveratrol analogs, which can be translated further to design novel highly potent anti-androgenic stilbenes.

Cleavage of Model Substrates by Arabidopsis thaliana PRORP1 Reveals New Insights into Its Substrate Requirements.

Two broad classes of RNase P trim the 5' leader of precursor tRNAs (pre-tRNAs): ribonucleoprotein (RNP)- and proteinaceous (PRORP)-variants. These two RNase P types, which use different scaffolds for catalysis, reflect independent evolutionary paths. While the catalytic RNA-based RNP form is present in all three domains of life, the PRORP family is restricted to eukaryotes. To obtain insights on substrate recognition by PRORPs, we examined the 5' processing ability of recombinant Arabidopsis thaliana PRORP1 (AtPRORP1) using a panel of pre-tRNASer variants and model hairpin-loop derivatives (pATSer type) that consist of the acceptor-T-stem stack and the T-/D-loop. Our data indicate the importance of the identity of N-1 (the residue immediately 5' to the cleavage site) and the N-1:N+73 base pair for cleavage rate and site selection of pre-tRNASer and pATSer. The nucleobase preferences that we observed mirror the frequency of occurrence in the complete suite of organellar pre-tRNAs in eight algae/plants that we analyzed. The importance of the T-/D-loop in pre-tRNASer for tight binding to AtPRORP1 is indicated by the 200-fold weaker binding of pATSer compared to pre-tRNASer, while the essentiality of the T-loop for cleavage is reflected by the near-complete loss of activity when a GAAA-tetraloop replaced the T-loop in pATSer. Substituting the 2'-OH at N-1 with 2'-H also resulted in no detectable cleavage, hinting at the possible role of this 2'-OH in coordinating Mg2+ ions critical for catalysis. Collectively, our results indicate similarities but also key differences in substrate recognition by the bacterial RNase P RNP and AtPRORP1: while both forms exploit the acceptor-T-stem stack and the elbow region in the pre-tRNA, the RNP form appears to require more recognition determinants for cleavage-site selection.

Characterization of diadzein-hemoglobin binding using optical spectroscopy and molecular dynamics simulations.

The present study establishes the effectiveness of natural drug delivery mechanisms and investigates the interactions between drug and its natural carrier. The binding between the isoflavone diadzein (DZN) and the natural carrier hemoglobin (HbA) was studied using optical spectroscopy and molecular dynamics simulations. The inherent fluorescence emission characteristics of DZN along with that of tryptophan (Trp) residues of the protein HbA were exploited to elucidate the binding location and other relevant parameters of the drug inside its delivery vehicle HbA. Stern-Volmer studies at different temperatures indicate that static along with collisional quenching mechanisms are responsible for the quenching of protein fluorescence by the drug. Molecular dynamics and docking studies supported the hydrophobic interactions between ligand and protein, as was observed from spectroscopy. DZN binds between the subunits of HbA, ∼15 Å away from the closest heme group of chain α1, emphasizing the fact that the drug does not interfere with oxygen binding site of HbA.

Numerical simulation and graphical analysis of in vitro benign tumor growth: application of single-particle state bosonic matter equation with length scaling.

We describe the application of a non-linear single-particle state bosonic condensate equation to simulate multicellular tumor growth by treating it as a coupling of two classical wave equations with real components. With one component representing the amplitude of the cells in their volume growth phase and the other representing the amplitude of the cells in their proliferation or mitosis phase, the two components of the coupled equation feed each other during the time evolution and are coupled together through diffusion and other linear and non-linear terms. The features of quiescent and necrotic cells, which result from poor nutrient diffusion into a tumor, have been found to correspond quite well to experimental data when they are modeled as depending on higher cell density. Classical hallmarks of benign tumor growth, such as the initial rapid growth, followed by a dramatic collapse in the proliferating cell count and a strong re-growth thereafter appear quite encouragingly in the theoretical results. A tool for graphical analysis of the tumor simulation results has been developed to provide morphological information about tumors at various growth stages. The model and the graphical analysis can be extended further to create an effective tool to predict/monitor tumor growth.

Effects of acclimation temperature and cadmium exposure on cellular energy budgets in the marine mollusk Crassostrea virginica: linking cellular and mitochondrial responses.

In order to understand the role of metabolic regulation in environmental stress tolerance, a comprehensive analysis of demand-side effects (i.e. changes in energy demands for basal maintenance) and supply-side effects (i.e. metabolic capacity to provide ATP to cover the energy demand) of environmental stressors is required. We have studied the effects of temperature (12, 20 and 28 degrees C) and exposure to a trace metal, cadmium (50 microg l(-1)), on the cellular energy budget of a model marine poikilotherm, Crassostrea virginica (eastern oysters), using oxygen demand for ATP turnover, protein synthesis, mitochondrial proton leak and non-mitochondrial respiration in isolated gill and hepatopancreas cells as demand-side endpoints and mitochondrial oxidation capacity, abundance and fractional volume as supply-side endpoints. Cadmium exposure and high acclimation temperatures resulted in a strong increase of oxygen demand in gill and hepatopancreas cells of oysters. Cd-induced increases in cellular energy demand were significant at 12 and 20 degrees C but not at 28 degrees C, possibly indicating a metabolic capacity limitation at the highest temperature. Elevated cellular demand in cells from Cd-exposed oysters was associated with a 2-6-fold increase in protein synthesis and, at cold acclimation temperatures, with a 1.5-fold elevated mitochondrial proton leak. Cellular aerobic capacity, as indicated by mitochondrial oxidation capacity, abundance and volume, did not increase in parallel to compensate for the elevated energy demand. Mitochondrial oxidation capacity was reduced in 28 degrees C-acclimated oysters, and mitochondrial abundance decreased in Cd-exposed oysters, with a stronger decrease (by 20-24%) in warm-acclimated oysters compared with cold-acclimated ones (by 8-13%). These data provide a mechanistic basis for synergism between temperature and cadmium stress on metabolism of marine poikilotherms. Exposure to combined temperature and cadmium stress may result in a strong energy deficiency due to the elevated energy demand on one hand and a reduced mitochondrial capacity to cover this demand on the other hand, which may have important implications for surviving seasonally and/or globally elevated temperatures in polluted estuaries.