Constraint-based, Single-point Approximate Kinetic Energy Functionals

We present a substantial extension of our constraint-based approach for development of orbital-free (OF) kinetic-energy (KE) density functionals intended for the calculation of quantum-mechanical forces in multi-scale molecular dynamics simulations. Suitability for realistic system simulations requires that the OF-KE functional yield accurate forces on the nuclei yet be relatively simple. We therefore require that the functionals be based on DFT constraints, local, dependent upon a small number of parameters fitted to a training set of limited size, and applicable beyond the scope of the training set. Our previous "modified conjoint" generalized-gradient-type functionals were constrained to producing a positive-definite Pauli potential. Though distinctly better than several published GGA-type functionals in that they gave semi-quantitative agreement with Born-Oppenheimer forces from full Kohn-Sham results, those modified conjoint functionals suffer from unphysical singularities at the nuclei. Here we show how to remove such singularities by introducing higher-order density derivatives. We give a simple illustration of such a functional used for the dissociation energy as a function of bond length for selected molecules.Comment: 16 pages, 9 figures, 2 tables, submitted to Phys. Rev.

Selective reduction of layers at low temperature in artificial superlattice thin films

Reduction and oxidation in transition-metal oxides are keys to develop technologies related to energy and the environment. Here we report the selective topochemical reduction observed when artificial superlattices with transition-metal oxides are treated at a temperature below 300 °C with CaH2. [CaFeO2]m/[SrTiO3]n infinite-layer/perovskite artificial superlattice thin films were obtained by low-temperature reduction of [CaFeO2.5]m/[SrTiO3]n brownmillerite/perovskite artificial superlattice thin films. By the reduction only the CaFeO2.5 layers in the artificial superlattices were reduced to the CaFeO2 infinite layers whereas the SrTiO3 layers were unchanged. The observed low-temperature reduction behaviors strongly suggest that the oxygen ion diffusion in the artificial superlattices is confined within the two-dimensional brownmillerite layers. The reduced artificial superlattice could be reoxidized, and thus, the selective reduction and oxidation of the constituent layers in the perovskite-structure framework occur reversibly

DNA Clasping by Mycobacterial HU: The C-Terminal Region of HupB Mediates Increased Specificity of DNA Binding

BACKGROUND: HU a small, basic, histone like protein is a major component of the bacterial nucleoid. E. coli has two subunits of HU coded by hupA and hupB genes whereas Mycobacterium tuberculosis (Mtb) has only one subunit of HU coded by ORF Rv2986c (hupB gene). One noticeable feature regarding Mtb HupB, based on sequence alignment of HU orthologs from different bacteria, was that HupB(Mtb) bears at its C-terminal end, a highly basic extension and this prompted an examination of its role in Mtb HupB function. METHODOLOGY/PRINCIPAL FINDINGS: With this objective two clones of Mtb HupB were generated; one expressing full length HupB protein (HupB(Mtb)) and another which expresses only the N terminal region (first 95 amino acid) of hupB (HupB(MtbN)). Gel retardation assays revealed that HupB(MtbN) is almost like E. coli HU (heat stable nucleoid protein) in terms of its DNA binding, with a binding constant (K(d)) for linear dsDNA greater than 1000 nM, a value comparable to that obtained for the HUalphaalpha and HUalphabeta forms. However CTR (C-terminal Region) of HupB(Mtb) imparts greater specificity in DNA binding. HupB(Mtb) protein binds more strongly to supercoiled plasmid DNA than to linear DNA, also this binding is very stable as it provides DNase I protection even up to 5 minutes. Similar results were obtained when the abilities of both proteins to mediate protection against DNA strand cleavage by hydroxyl radicals generated by the Fenton's reaction, were compared. It was also observed that both the proteins have DNA binding preference for A:T rich DNA which may occur at the regulatory regions of ORFs and the oriC region of Mtb. CONCLUSIONS/SIGNIFICANCE: These data thus point that HupB(Mtb) may participate in chromosome organization in-vivo, it may also play a passive, possibly an architectural role

Concentration-dependent organization of DNA by the dinoflagellate histone-like protein HCc3

The liquid crystalline chromosomes of dinoflagellates are the alternative to the nucleosome-based organization of chromosomes in the eukaryotes. These nucleosome-less chromosomes have to devise novel ways to maintain active parts of the genome. The dinoflagellate histone-like protein HCc3 has significant sequence identity with the bacterial DNA-binding protein HU. HCc3 also has a secondary structure resembling HU in silico. We have examined HCc3 in its recombinant form. Experiments on DNA-cellulose revealed its DNA-binding activity is on the C-terminal domain. The N-terminal domain is responsible for intermolecular oligomerization as demonstrated by cross-linking studies. However, HCc3 could not complement Escherichia coli HU-deficient mutants, suggesting functional differences. In ligation assays, HCc3-induced DNA concatenation but not ring closure as the DNA-bending HU does. The basic HCc3 was an efficient DNA condensing agent, but it did not behave like an ordinary polycationic compound. HCc3 also induced specific structures with DNA in a concentration-dependent manner, as demonstrated by atomic force microscopy (AFM). At moderate concentration of HCc3, DNA bridging and bundling were observed; at high concentrations, the complexes were even more condensed. These results are consistent with a biophysical role for HCc3 in maintaining extended DNA loops at the periphery of liquid crystalline chromosomes

Concentration-dependent organization of DNA by the dinoflagellate histone-like protein HCc3

The liquid crystalline chromosomes of dinoflagellates are the alternative to the nucleosome-based organization of chromosomes in the eukaryotes. These nucleosome-less chromosomes have to devise novel ways to maintain active parts of the genome. The dinoflagellate histone-like protein HCc3 has significant sequence identity with the bacterial DNA-binding protein HU. HCc3 also has a secondary structure resembling HU in silico. We have examined HCc3 in its recombinant form. Experiments on DNA-cellulose revealed its DNA-binding activity is on the C-terminal domain. The N-terminal domain is responsible for intermolecular oligomerization as demonstrated by cross-linking studies. However, HCc3 could not complement Escherichia coli HU-deficient mutants, suggesting functional differences. In ligation assays, HCc3-induced DNA concatenation but not ring closure as the DNA-bending HU does. The basic HCc3 was an efficient DNA condensing agent, but it did not behave like an ordinary polycationic compound. HCc3 also induced specific structures with DNA in a concentration-dependent manner, as demonstrated by atomic force microscopy (AFM). At moderate concentration of HCc3, DNA bridging and bundling were observed; at high concentrations, the complexes were even more condensed. These results are consistent with a biophysical role for HCc3 in maintaining extended DNA loops at the periphery of liquid crystalline chromosomes

Interplay of Protein and DNA Structure Revealed in Simulations of the lac Operon

The E. coli Lac repressor is the classic textbook example of a protein that attaches to widely spaced sites along a genome and forces the intervening DNA into a loop. The short loops implicated in the regulation of the lac operon suggest the involvement of factors other than DNA and repressor in gene control. The molecular simulations presented here examine two likely structural contributions to the in-vivo looping of bacterial DNA: the distortions of the double helix introduced upon association of the highly abundant, nonspecific nucleoid protein HU and the large-scale deformations of the repressor detected in low-resolution experiments. The computations take account of the three-dimensional arrangements of nucleotides and amino acids found in crystal structures of DNA with the two proteins, the natural rest state and deformational properties of protein-free DNA, and the constraints on looping imposed by the conformation of the repressor and the orientation of bound DNA. The predicted looping propensities capture the complex, chain-length-dependent variation in repression efficacy extracted from gene expression studies and in vitro experiments and reveal unexpected chain-length-dependent variations in the uptake of HU, the deformation of repressor, and the folding of DNA. Both the opening of repressor and the presence of HU, at levels approximating those found in vivo, enhance the probability of loop formation. HU affects the global organization of the repressor and the opening of repressor influences the levels of HU binding to DNA. The length of the loop determines whether the DNA adopts antiparallel or parallel orientations on the repressor, whether the repressor is opened or closed, and how many HU molecules bind to the loop. The collective behavior of proteins and DNA is greater than the sum of the parts and hints of ways in which multiple proteins may coordinate the packaging and processing of genetic information. © 2013 Czapla et al

Transcriptional correlates of the pathological phenotype in a Huntington’s disease mouse model

Huntington disease (HD) is a fatal neurodegenerative disorder without a cure that is caused by an aberrant expansion of CAG repeats in exon 1 of the huntingtin (HTT) gene. Although a negative correlation between the number of CAG repeats and the age of disease onset is established, additional factors may contribute to the high heterogeneity of the complex manifestation of symptoms among patients. This variability is also observed in mouse models, even under controlled genetic and environmental conditions. To better understand this phenomenon, we analysed the R6/1 strain in search of potential correlates between pathological motor/cognitive phenotypical traits and transcriptional alterations. HD-related genes (e.g., Penk, Plk5, Itpka), despite being downregulated across the examined brain areas (the prefrontal cortex, striatum, hippocampus and cerebellum), exhibited tissue-specific correlations with particular phenotypical traits that were attributable to the contribution of the brain region to that trait (e.g., striatum and rotarod performance, cerebellum and feet clasping). Focusing on the striatum, we determined that the transcriptional dysregulation associated with HD was partially exacerbated in mice that showed poor overall phenotypical scores, especially in genes with relevant roles in striatal functioning (e.g., Pde10a, Drd1, Drd2, Ppp1r1b). However, we also observed transcripts associated with relatively better outcomes, such as Nfya (CCAAT-binding transcription factor NF-Y subunit A) plus others related to neuronal development, apoptosis and differentiation. In this study, we demonstrated that altered brain transcription can be related to the manifestation of HD-like symptoms in mouse models and that this can be extrapolated to the highly heterogeneous population of HD patients

Analysis of In-Vivo LacR-Mediated Gene Repression Based on the Mechanics of DNA Looping

Interactions of E. coli lac repressor (LacR) with a pair of operator sites on the same DNA molecule can lead to the formation of looped nucleoprotein complexes both in vitro and in vivo. As a major paradigm for loop-mediated gene regulation, parameters such as operator affinity and spacing, repressor concentration, and DNA bending induced by specific or non-specific DNA-binding proteins (e.g., HU), have been examined extensively. However, a complete and rigorous model that integrates all of these aspects in a systematic and quantitative treatment of experimental data has not been available. Applying our recent statistical-mechanical theory for DNA looping, we calculated repression as a function of operator spacing (58–156 bp) from first principles and obtained excellent agreement with independent sets of in-vivo data. The results suggest that a linear extended, as opposed to a closed v-shaped, LacR conformation is the dominant form of the tetramer in vivo. Moreover, loop-mediated repression in wild-type E. coli strains is facilitated by decreased DNA rigidity and high levels of flexibility in the LacR tetramer. In contrast, repression data for strains lacking HU gave a near-normal value of the DNA persistence length. These findings underscore the importance of both protein conformation and elasticity in the formation of small DNA loops widely observed in vivo, and demonstrate the utility of quantitatively analyzing gene regulation based on the mechanics of nucleoprotein complexes

Metabolic labeling of RNA uncovers principles of RNA production and degradation dynamics in mammalian cells

available in PMC 2011 November 01.Cellular RNA levels are determined by the interplay of RNA production, processing and degradation. However, because most studies of RNA regulation do not distinguish the separate contributions of these processes, little is known about how they are temporally integrated. Here we combine metabolic labeling of RNA at high temporal resolution with advanced RNA quantification and computational modeling to estimate RNA transcription and degradation rates during the response of mouse dendritic cells to lipopolysaccharide. We find that changes in transcription rates determine the majority of temporal changes in RNA levels, but that changes in degradation rates are important for shaping sharp 'peaked' responses. We used sequencing of the newly transcribed RNA population to estimate temporally constant RNA processing and degradation rates genome wide. Degradation rates vary significantly between genes and contribute to the observed differences in the dynamic response. Certain transcripts, including those encoding cytokines and transcription factors, mature faster. Our study provides a quantitative approach to study the integrative process of RNA regulation.Human Frontier Science Program (Strasbourg, France)Howard Hughes Medical InstituteBurroughs Wellcome Fund (Career Award at the Scientific Interface

RhoGTPase Regulators Orchestrate Distinct Stages of Synaptic Development

Small RhoGTPases regulate changes in post-synaptic spine morphology and density that support learning and memory. They are also major targets of synaptic disorders, including Autism. Here we sought to determine whether upstream RhoGTPase regulators, including GEFs, GAPs, and GDIs, sculpt specific stages of synaptic development. The majority of examined molecules uniquely regulate either early spine precursor formation or later matura- tion. Specifically, an activator of actin polymerization, the Rac1 GEF β-PIX, drives spine pre- cursor formation, whereas both FRABIN, a Cdc42 GEF, and OLIGOPHRENIN-1, a RhoA GAP, regulate spine precursor elongation. However, in later development, a novel Rac1 GAP, ARHGAP23, and RhoGDIs inactivate actomyosin dynamics to stabilize mature synap- ses. Our observations demonstrate that specific combinations of RhoGTPase regulatory pro- teins temporally balance RhoGTPase activity during post-synaptic spine development