Search

Summary The FACT (facilitates chromatin transcription) complex, comprising SPT16 and SSRP1, conducts structural alterations during nucleosome unwrapping. Our previous cryoelectron microscopic (cryo-EM) analysis revealed the first intermediate structure of an unwrapped nucleosome with human FACT, in which 112-bp DNA and the phosphorylated intrinsically disordered (pAID) segment of SPT16 jointly wrapped around the histone core instead of 145-bp DNA. Using NMR, here we clarified that the histone H3 N-terminal tails, unobserved in the cryo-EM structure, adopt two different conformations reflecting their asymmetric locations at entry/exit sites: one corresponds to the original nucleosome site buried in two DNA gyres (DNA side), whereas the other, comprising pAID and DNA, is more exposed to the solvent (pAID side). NMR real-time monitoring showed that H3 acetylation is faster on the pAID side than on the DNA side. Our findings highlight that accessible conformations of H3 tails are created by the replacement of nucleosomal DNA with pAID., Graphical Abstract, Highlights • H3 N-tail, restricted to two DNA gyres of nucleosome, is protected from Gcn5 • H3 N-tail is dynamically exposed by replacement of nucleosomal DNA with pAID of FACT • Gcn5 efficiently acetylates accessible H3 N-tail of nucleosome with FACT • FACT acts as a modulator for dynamic behavior of H3 tails in nucleosome, Biochemistry; Molecular Biology; Structural Biology

Background Human ether-à-go-go-related gene potassium channel 1 (hERG) is a voltage-gated potassium channel, the voltage-sensing domain (VSD) of which is targeted by a gating-modifier toxin, APETx1. APETx1 is a 42-residue peptide toxin of sea anemone Anthopleura elegantissima and inhibits hERG by stabilizing the resting state. A previous study that conducted cysteine-scanning analysis of hERG identified two residues in the S3-S4 region of the VSD that play important roles in hERG inhibition by APETx1. However, mutational analysis of APETx1 could not be conducted as only natural resources have been available until now. Therefore, it remains unclear where and how APETx1 interacts with the VSD in the resting state. Results We established a method for preparing recombinant APETx1 and determined the NMR structure of the recombinant APETx1, which is structurally equivalent to the natural product. Electrophysiological analyses using wild type and mutants of APETx1 and hERG revealed that their hydrophobic residues, F15, Y32, F33, and L34, in APETx1, and F508 and I521 in hERG, in addition to a previously reported acidic hERG residue, E518, play key roles in the inhibition of hERG by APETx1. Our hypothetical docking models of the APETx1-VSD complex satisfied the results of mutational analysis. Conclusions The present study identified the key residues of APETx1 and hERG that are involved in hERG inhibition by APETx1. These results would help advance understanding of the inhibitory mechanism of APETx1, which could provide a structural basis for designing novel ligands targeting the VSDs of KV channels.

Dysregulation of repressor-element 1 silencing transcription factor REST/NRSF is related to several neuropathies, including medulloblastoma, glioblastoma, Huntington’s disease, and neuropathic pain. Inhibitors of the interaction between the N-terminal repressor domain of REST/NRSF and the PAH1 domain of its corepressor mSin3 may ameliorate such neuropathies. In-silico screening based on the complex structure of REST/NRSF and mSin3 PAH1 yielded 52 active compounds, including approved neuropathic drugs. We investigated their binding affinity to PAH1 by NMR, and their inhibitory activity toward medulloblastoma cell growth. Interestingly, three antidepressant and antipsychotic medicines, sertraline, chlorprothixene, and chlorpromazine, were found to strongly bind to PAH1. Multivariate analysis based on NMR chemical shift changes in PAH1 residues induced by ligand binding was used to identify compound characteristics associated with cell growth inhibition. Active compounds showed a new chemo-type for inhibitors of the REST/NRSF-mSin3 interaction, raising the possibility of new therapies for neuropathies caused by dysregulation of REST/NRSF.

Bacteriorhodopsin (BR) is a model protein for light-driven proton pumps, where the vectorial active proton transport results in light-energy conversion. To clarify the microscopic mechanism of primary proton transfer from retinal Schiff base (SB) to Asp85 in BR, herein we performed quantum-mechanical metadynamics simulations of the whole BR system (∼3800 atoms). The simulations showed a novel proton transfer mechanism, viz. hydroxide ion mechanism, in which the deprotonation of specific internal water (Wat452) yields the protonation of Asp85 via Thr89, after which the resulting hydroxide ion accepts the remaining proton from retinal SB. Furthermore, systematic investigations adopting four sequential snapshots obtained by the time-resolved serial femtosecond crystallography revealed that proton transfer took 2–5.25 μs on the photocycle. The presence of Wat401, which is the main difference between snapshots at 2 and 5.25 μs, is found to be essential in assisting the primary proton transfer.SIGNIFICANCEBacteriorhodopsin (BR), the benchmark of light-driven proton pumps, has attracted much attention from diverse areas in terms of energy conversion. Despite the significant experimental and theoretical efforts, the microscopic mechanism of the proton transfers in BR is not completely unveiled. In this study, quantum-mechanical molecular dynamics simulations of whole BR system were performed to elucidate the primary proton transfer in the L intermediate state with the latest snapshots obtained from X-ray free electron laser. As a result, it is found that the hydroxide ion originating from the specific internal water, which appears at the active site only in the L state, acts as a carrier for the primary proton transfer, demonstrating the importance of hydroxide ions in proton pumps.

Histone tails, which protrude from nucleosome core particles (NCPs), play crucial roles in the regulation of DNA transcription, replication, and repair. In this study, structural diversity of nucleosomes was investigated in detail by analyzing the observed charge states of nucleosomes reconstituted with various lengths of DNA using positive-mode electrospray ionization mass spectrometry (ESI-MS) and molecular dynamics (MD) simulation. Here, we show that canonical NCPs, having 147 bp DNA closely wrapped around a histone octamer, can be classified into three groups by charge state, with the least-charged group being more populated than the highly charged and intermediate groups. Ions with low charge showed small collision cross sections (CCSs), suggesting that the histone tails are generally compact in the gas phase, whereas the minor populations with higher charges appeared to have more loosened structure. Overlapping dinucleosomes, which contain 14 histone proteins closely packed with 250 bp DNA, showed similar characteristics. In contrast, mononucleosomes reconstituted with a histone octamer and longer DNA (≥250 bp), which have DNA regions uninvolved in the core-structure formation, showed only low-charge ions. This was also true for dinucleosomes with free DNA regions. These results suggest that free DNA regions affect the nucleosome structures. To investigate the possible structures of NCP observed in ESI-MS, computational structural calculations in solution and in vacuo were performed. They suggested that conformers with large CCS values have slightly loosened structure with extended tail regions, which might relate to the biological function of histone tails.

The structural, dynamical, and energetic properties of the excess proton in ice were studied using density-functional tight-binding molecular dynamics simulations. The ice systems investigated herein consisted of low-density hexagonal and cubic crystalline variants (ice Ih and Ic) and high-density structures (ice III and melted ice VI). Analysis of the temperature dependence of radial distribution function and bond order parameters served to characterize the distribution and configuration of hundreds of water molecules in a unit cell. We confirmed that ice Ih and Ic possess higher hexagonal symmetries than ice III and melted ice VI. The estimated Grotthuss shuttling diffusion coefficients in ice were larger than that of liquid water, indicating a slower proton diffusion process in high-density structures than in low-density systems. The energy barriers calculated on the basis of the Arrhenius plot of diffusion coefficients were in reasonable agreement with experimental measurement for ice Ih. Furthermore,...

Predicting pKa values for different types of amine species with high accuracy and efficiency is of critical importance for the design of high performance and economical solvents in carbon capture and storage with aqueous amine solutions. In this study, we demonstrate that density-functional tight-binding-based metadynamics simulations are a promising approach to calculate the free energy difference between the protonated and neutral states of amines in aqueous solution with inexpensive computational cost. The calculated pKa values were in satisfactory agreement with the experimental values, the mean absolute deviation being only 0.09 pKa units for 34 amines commonly used in CO2 scrubbing. Such superior reproducibility and correlation compared to estimations by static quantum mechanical calculations highlight the significant effect of dynamical proton transfer processes in explicit solvent molecules for the improvement of the estimation accuracy.

Divide-and-conquer-type density-functional tight-binding molecular dynamics simulations of the CO2 absorption process in monoethanolamine (MEA) solution have been performed for systems containing thousands of atoms. The formation of carbamate anions has been widely investigated for neutral systems via ab initio molecular dynamics simulations, yet the present study is aimed at identifying the role of hydroxide ions in acid-base equilibrium. The structural and electronic analyses reveal that the hydroxide ion approaches, via Grotthuss-type shuttling, the zwitterionic intermediates and abstracts a proton from the nitrogen atom of MEA. We also estimated the fraction of reacted CO2 and carbamate formed at different initial CO2 concentrations that confirm a high absorbed CO2 concentration decreases the fraction of MEA(C) formed due to the abundance of MEA(Z) in the solution.

A low-computational-cost algorithm and its parallel implementation for periodic divide-and-conquer density-functional tight-binding (DC-DFTB) calculations are presented. The developed algorithm enables rapid computation of the interaction between atomic partial charges, which is the bottleneck for applications to large systems, by means of multipole- and interpolation-based approaches for long- and short-range contributions. The numerical errors of energy and forces with respect to the conventional Ewald-based technique can be under the control of the multipole expansion order, level of unit cell replication, and interpolation grid size. The parallel performance of four different evaluation schemes combining previous approaches and the proposed one are assessed using test calculations of a cubic water box on the K computer. The largest benchmark system consisted of 3,295,500 atoms. DC-DFTB energy and forces for this system were obtained in only a few minutes when the proposed algorithm was activated and parallelized over 16,000 nodes in the K computer. The high performance using a single node workstation was also confirmed. In addition to liquid water systems, the feasibility of the present method was examined by testing solid systems such as diamond form of carbon, face-centered cubic form of copper, and rock salt form of sodium chloride. © 2017 Wiley Periodicals, Inc.

The neuron-restrictive silencing factor NRSF/REST binds to neuron-restrictive silencing elements in neuronal genes and recruits corepressors such as mSin3 to inhibit epigenetically neuronal gene expression. Because dysregulation of NRSF/REST is related to neuropathic pain, here, we have designed compounds to target neuropathic pain based on the mSin3-binding helix structure of NRSF/REST and examined their ability to bind to mSin3 by NMR. One compound, mS-11, binds strongly to mSin3 with a binding mode similar to that of NRSF/REST. In a mouse model of neuropathic pain, mS-11 was found to ameliorate abnormal pain behavior and to reverse lost peripheral morphine analgesia. Furthermore, even in the less well epigenetically defined case of fibromyalgia, mS-11 ameliorated symptoms in a mouse model, suggesting that fibromyalgia is related to the dysfunction of NRSF/REST. Taken together, these findings show that the chemically optimized mimetic mS-11 can inhibit mSin3-NRSF/REST binding and successfully reverse lost peripheral and central morphine analgesia in mouse models of pain.

Nucleotide excision repair is initiated by two different damage recognition subpathways, global genome repair (GGR) and transcription-coupled repair (TCR). In GGR, XPC detects DNA lesions and recruits TFIIH via interaction with the pleckstrin homology (PH) domain of TFIIH subunit p62. In TCR, an elongating form of RNA Polymerase II detects a lesion on the transcribed strand and recruits TFIIH by an unknown mechanism. Here, we found that the TCR initiation factor UVSSA forms a stable complex with the PH domain of p62 via a short acidic string in the central region of UVSSA, and determined the complex structure by NMR. The acidic string of UVSSA binds strongly to the basic groove of the PH domain by inserting Phe408 and Val411 into two pockets, highly resembling the interaction mechanism of XPC with p62. Mutational binding analysis validated the structure and identified residues crucial for binding. TCR activity was markedly diminished in UVSSA-deficient cells expressing UVSSA mutated at Phe408 or Val411. Thus, a common TFIIH recruitment mechanism is shared by UVSSA in TCR and XPC in GGR.

Nucleosomes in contact In eukaryotic cells, genomic DNA must be compacted to fit inside the nucleus. A key player in DNA packaging is the nucleosome, which comprises a segment of DNA wrapped around an octamer of histone proteins. During replication and transcription, nucleosomes must reposition themselves on the DNA. In this process, nucleosomes can collide to form a dinucleosome. Kato et al. report a high-resolution crystal structure of a dinucleosome. One of the octamers has lost a histone dimer so that the dinucleosome comprises an octamer and a hexamer. The structure may represent an intermediate during chromatin remodeling. Science , this issue p. 205

The physisorption of an acetone molecule on hexagonal graphene nanoflakes with increasing size has been investigated using a variety of quantum chemical methods capable of describing weak intermolecular interactions: coupled-cluster theory (CCSD and CCSD(T)), second-order Moller–Plesset perturbation theory with and without spin-component scaling (SCS-MP2 or standard MP2), long-range corrected density functional theory combined with a van der Waals functional (LC-BOP+ALL), meta-generalized gradient approximation functionals (M06-2X and M05-2X), and the dispersion augmented self-consistent-charge density functional tight-binding (SCC-DFTB-D) method. Our benchmark results for model systems as large as dicircumcoronene C96H24 have confirmed the suitability of the SCS-MP2 method for this specific system and the satisfactory performance of the computationally much more economical semiempirical SCC-DFTB-D method. The latter delivers a qualitatively accurate description of physisorption for flakes containing more...

The diffusion of the hydroxide ion in bulk water was examined by linear-scaling divide-and-conquer density-functional tight-binding molecular dynamics (DC-DFTB-MD) simulations using three different-sized unit cells that contained 522, 1050, and 4999 water molecules as well as one hydroxide ion. The repulsive potential for the oxygen-oxygen pair was improved by iterative Boltzmann inversion, which adjusted the radial distribution function of DFTB-MD simulations to that of the reference density functional theory-MD one. The calculated diffusion coefficients and the Arrhenius diffusion barrier were in good agreement with experimental results. The results of the hydroxide ion coordination number distribution and potential of mean force analyses supported a dynamical hypercoordination diffusion mechanism.

The heterodimeric transcription factor E2F1-DP1 plays crucial roles in coordinating gene expression during G1/S cell cycle progression. For transcriptional activation, the transactivation domain (TAD) of E2F1 is known to interact with the TATA-binding protein of TFIID and the p62 subunit of TFIIH. It is generally believed that DP1 facilitates E2F1 binding to target DNA and does not possess a TAD. Here, we show that an acidic region of DP1, whose function has remained elusive, binds to the plekstrin homology (PH) domain of p62 with higher affinity than that of E2F1 and contributes to transcriptional activation. The structure of the complex revealed that DP1 forms a twisted U-shaped, string-like conformation and binds to the surface of the PH domain by anchoring Phe403 into a pocket in the PH domain. The transcriptional activity of E2F1-DP1 was reduced when Phe403 of DP1 was mutated. These findings indicate that the acidic region of DP1 acts as a TAD by contacting TFIIH.

Aqueous Na- or K-ion batteries could virtually eliminate the safety and cost concerns raised from Li-ion batteries, but their widespread applications have generally suffered from narrow electrochemical potential window (ca. 1.23 V) of aqueous electrolytes that leads to low energy density. Herein, by exploring optimized eutectic systems of Na and K salts with asymmetric imide anions, we discovered, for the first time, room-temperature hydrate melts for Na and K systems, which are the second and third alkali metal hydrate melts reported since the first discovery of Li hydrate melt by our group in 2016. The newly discovered Na- and K- hydrate melts could significantly extend the potential window up to 2.7 and 2.5 V (at Pt electrode), respectively, owing to the merit that almost all water molecules participate in the Na+ or K+ hydration shells. As a proof-of-concept, a prototype Na3 V2 (PO4 )2 F3 |NaTi2 (PO4 )3 aqueous Na-ion full-cell with the Na-hydrate-melt electrolyte delivers an average discharge voltage of 1.75 V, that is among the highest value ever reported for all aqueous Na-ion batteries.

The present study implemented the divide-and-conquer time-dependent density-functional tight-binding (DC-TDDFTB) code on a graphical processing unit (GPU). The DC method, which is a linear-scaling scheme, divides a total system into several fragments. By separately solving local equations in individual fragments, the DC method could reduce slow central processing unit (CPU)-GPU memory access, as well as computational cost, and avoid shortfalls of GPU memory. Numerical applications confirmed that the present code on GPU significantly accelerated the TDDFTB calculations, while maintaining accuracy. Furthermore, the DC-TDDFTB simulation of 2-acetylindan-1,3-dione displays excited-state intramolecular proton transfer and provides reasonable absorption and fluorescence energies with the corresponding experimental values. © 2019 Wiley Periodicals, Inc.

In this study, the divide-and-conquer (DC) method was extended to time-dependent density functional tight-binding (TDDFTB) theory to enable excited-state calculations of large systems and is denoted by DC-TDDFTB. The efficient diagonalization algorithms of TDDFTB and DC-TDDFTB methods were implemented into our in-house program. Test calculations of polyethylene aldehyde and p-coumaric acid, a pigment in photoactive yellow protein, in water demonstrate the high accuracy and efficiency of the developed DC-TDDFTB method. Furthermore, the (TD)DFTB metadynamics simulations of acridinium in the ground and excited states give reasonable p Ka values compared with the corresponding experimental values.

General transcription factor II E (TFIIE) contains an acid-rich region (residues 378–393) in its α -subunit, comprising 13 acidic and two hydrophobic (Phe387 and Val390) residues. Upon binding to the p62 subunit of TFIIH, the acidic region adopts an extended string-like structure on the basic groove of the pleckstrin homology domain (PHD) of p62, and inserts Phe387 and Val390 into two shallow pockets in the groove. Here, we have examined the dynamics of this interaction by NMR and molecular dynamics (MD) simulations. Although alanine substitution of Phe387 and/or Val390 greatly reduced binding to PHD, the binding mode of the mutants was similar to that of the wild-type, as judged by the chemical-shift changes of the PHD. NMR relaxation dispersion profiles of the interaction exhibited large amplitudes for residues in the C-terminal half-string in the acidic region (Phe387, Glu388, Val390, Ala391, and Asp392), indicating a two-site binding mode: one corresponding to the final complex structure, and one to an off-pathway minor complex. To probe the off-pathway complex structure, an atomically detailed free-energy landscape of the binding mode was computed by all-atom multicanonical MD. The most thermodynamically stable cluster corresponded to the final complex structure. One of the next stable clusters was the off-pathway structure cluster, showing the reversed orientation of the C-terminal half-string on the PHD groove, as compared with the final structure. MD calculations elucidated that the C-terminal half-acidic-string forms encounter complexes mainly around the positive groove region with nearly two different orientations of the string, parallel and antiparallel to the final structure. Interestingly, the most encountered complexes exhibit a parallel-like orientation, suggesting that the string has a tendency to bind around the groove in the proper orientation with the aid of Phe387 and/or Val390 to proceed smoothly to the final complex structure.

Objectives Mustard gas (MG) has been the most widely used chemical warfare agent in the past century. However, few but conflicting data exist on the effects of MG exposure on long-term mortality. We investigated MG-related mortality in retired workers at a poisonous gas factory. Methods We assessed mortality rates among 2392 male and 1226 female workers, whose vital status could be determined through 31 December 2009, at a poisonous gas factory operating from 1929 to 1945 in Okuno-jima, Japan. The analysis employed standardised mortality ratios (SMRs) calculated using national and prefectural references and a Cox proportional hazard regression model. Applying the Kaplan-Meier method, we compared cumulative death rates in the study cohort stratified by an ‘Okuno-jima MG Index’ which represented the product of HRs derived for job category and length of service. Results Among male workers, we found significant excesses in mortality from upper respiratory tract cancer (SMR 3.06), liver cancer (1.67), lung cancer (2.01) and chronic bronchitis/emphysema (4.84) compared with the national population, as well as stomach cancer (1.20) versus the Hiroshima Prefecture population. When stratified into 3 subgroups by the Okuno-jima MG Index, those with a higher Okuno-jima MG Index had significantly higher cumulative rates of death from respiratory cancer and chronic bronchitis/emphysema. Conclusions MG exposure significantly increases the long-term risk of death from respiratory cancer and chronic bronchitis/emphysema. The Okuno-jima MG Index may be a useful indicator for estimating cumulative MG exposure.

Nucleosome assembly protein 1 (NAP1) binds both the (H3-H4)2 tetramer and two H2A-H2B dimers, mediating their sequential deposition on DNA. NAP1 contains a C-terminal acidic domain (CTAD) and a core domain that promotes dimer formation. Here, we have investigated the roles of the core domain and CTAD of human NAP1 in binding to H2A-H2B and H3-H4 by isothermal calorimetry and native mass spectrometry and compared them with the roles of yeast NAP1. We show that the hNAP1 and yNAP1 dimers bind H2A-H2B by two different modes: a strong endothermic interaction and a weak exothermic interaction. A mutant hNAP1, but not yNAP1, dimer lacking CTAD loses the exothermic interaction and shows greatly reduced H2A-H2B binding activity. The isolated CTAD of hNAP1 binds H2A-H2B only exothermically with relatively stronger binding as compared with the exothermic interaction observed for the full-length hNAP1 dimer. Thus, the two CTADs in the hNAP1 dimer seem to provide binding assistance for the strong endothermic interaction of the core domain with H2A-H2B. By contrast, in the relatively weaker binding of hNAP1 to H3-H4 as compared with yNAP1, CTAD of hNAP1 has no significant role. To our knowledge, this is the first distinct role identified for the hNAP1 CTAD.

SummaryIn global genome repair (GGR), XPC detects damaged nucleotides and recruits TFIIH complex. The small acidic region of XPC binds to the pleckstrin homology (PH) domain of TFIIH subunit p62; however, the recognition mechanism remains elusive. Here, we use nuclear magnetic resonance to present the tertiary structure of XPC bound to the PH domain. The XPC acidic region forms a long string stabilized by insertion of Trp133 and Val136 into two separate hollows of the PH domain, coupled with extensive electrostatic contacts. Analysis of several XPC mutants revealed that particularly Trp133 is essential for binding to the PH domain. In cell lines stably expressing mutant XPC, alanine substitution at Trp133 or Trp133/Val136 compromised UV resistance, recruitment of TFIIH to DNA damage, and removal of UV-induced photoproducts from genomic DNA. These findings show how TFIIH complex is recruited by XPC to damaged DNA, advancing our understanding of the early stage of GGR.

Summary: RAD52 mediates homologous recombination by annealing cDNA strands. However, the detailed mechanism of DNA annealing promoted by RAD52 has remained elusive. Here we report two crystal structures of human RAD52 single-stranded DNA (ssDNA) complexes that probably represent key reaction intermediates of RAD52-mediated DNA annealing. The first structure revealed a “wrapped” conformation of ssDNA around the homo-oligomeric RAD52 ring, in which the edges of the bases involved in base pairing are exposed to the solvent. The ssDNA conformation is close to B-form and appears capable of engaging in Watson-Crick base pairing with the cDNA strand. The second structure revealed a “trapped” conformation of ssDNA between two RAD52 rings. This conformation is stabilized by a different RAD52 DNA binding site, which promotes the accumulation of multiple RAD52 rings on ssDNA and the aggregation of ssDNA. These structures provide a structural framework for understanding the mechanism of RAD52-mediated DNA annealing. : Biomolecules; Molecular Genetics; Structural Biology Subject Areas: Biomolecules, Molecular Genetics, Structural Biology

The structural, dynamical, and energetic properties of the excess proton in ice were studied using density-functional tight-binding molecular dynamics simulations. The ice systems investigated herein consisted of low-density hexagonal and cubic crystalline variants (ice I

Predicting pK

It is well known that various modifications of histone tails play important roles in the regulation of transcription initiation. In this study, some lysine (Lys) and arginine (Arg) residues were acetylated and deiminated, respectively, in the histone H2A/H2B dimer, and charge-neutralization effects on the dimer structure were studied by native mass spectrometry. Given that both acetylation and deimination neutralize the positive charges of basic amino acid residues, it had been expected that these modifications would correspondingly affect the gas-phase behavior of the histone H2A/H2B dimer. Contrary to this expectation, it was found that Arg deimination led to greater difficulty of dissociation of the dimer by gas-phase collision, whereas acetylation of Lys residues did not cause such a drastic change in the dimer stability. In contrast, ion mobility-mass spectrometry (IM-MS) experiments showed that arrival times in the mobility cell both of acetylated and of deiminated dimer ions changed little from those of the unmodified dimer ions, indicating that the sizes of the dimer ions did not change by modification. Charge neutralization of Arg, basicity of which is higher than Lys, might have triggered some alteration of the dimer structure that cannot be found in IM-MS but can be detected by collision in the gas phase.

Telomeres are DNA-protein complexes located at the ends of linear eukaryotic chromosomes, and are essential for chromosome stability and maintenance. In most organisms, telomeres consist of tandemly repeated sequences of guanine-clusters. In higher eukaryotes, most of the telomeric repeat regions are tightly packaged into nucleosomes, even though telomeric repeats act as nucleosome-disfavoring sequences. Although telomeres were considered to be condensed heterochromatin structures, recent studies revealed that the chromatin structures in telomeres are actually dynamic. The dynamic properties of telomeric chromatin are considered to be important for the structural changes between the euchromatic and heterochromatic states during the cell cycle and in cellular differentiation. We propose that the nucleosome-disfavoring property of telomeric repeats is a crucial determinant for the lability of telomeric nucleosomes, and provides a platform for chromatin dynamics in telomeres. Furthermore, we discuss the influences of telomeric components on the nucleosome organization and chromatin dynamics in telomeres.

The histone H2A/H2B dimer is a component of nucleosome core particles (NCPs). The structure of the dimer at the atomic level has not yet been revealed. A possible reason for this is that the dimer has three intrinsically disordered tail regions: the N- and C-termini of H2A and the N-terminus of H2B. To investigate the role of the tail regions of the H2A/H2B dimer structure, we characterized behaviors of the H2A/H2B mutant dimers, in which these functionally important disordered regions were depleted, using mass spectrometry (MS). After verifying that the acetylation of Lys residues in the tail regions had little effect on the gas-phase conformations of the wild-type dimer, we prepared two histone H2A/H2B dimer mutants: an H2A/H2B dimer depleted of both N-termini (dN-H2A/dN-H2B) and a dimer with the N- and C-termini of H2A and the N-terminus of H2B depleted (dNC-H2A/dN-H2B). We analyzed these mutants using ion mobility-mass spectrometry (IM-MS) and hydrogen/deuterium exchange mass spectrometry (HDX-MS). With IM-MS, reduced structural diversity was observed for each of the tail-truncated H2A/H2B mutants. In addition, global HDX-MS proved that the dimer mutant dNC-H2A/dN-H2B was susceptible to deuteration, suggesting that its structure in solution was somewhat loosened. A partial relaxation of the mutant's structure was demonstrated also by IM-MS. In this study, we characterized the relationship between the tail lengths and the conformations of the H2A/H2B dimer in solution and gas phases, and demonstrated, using mass spectrometry, that disordered tail regions play an important role in stabilizing the conformation of the core region of the dimer in both phases.

Heterochromatin protein 1 (HP1) is an evolutionarily conserved chromosomal protein that binds to lysine 9-methylated histone H3 (H3K9me), a hallmark of heterochromatin. Although HP1 phosphorylation has been described in several organisms, the biological implications of this modification remain largely elusive. Here we show that HP1's phosphorylation has a critical effect on its nucleosome binding properties. By in vitro phosphorylation assays and conventional chromatography, we demonstrated that casein kinase II (CK2) is the kinase primarily responsible for phosphorylating the N-terminus of human HP1α. Pull-down assays using in vitro-reconstituted nucleosomes showed that unmodified HP1α bound H3K9-methylated and H3K9-unmethylated nucleosomes with comparable affinity, whereas CK2-phosphorylated HP1α showed a high specificity for H3K9me3-modified nucleosomes. Electrophoretic mobility shift assays showed that CK2-mediated phosphorylation diminished HP1α's intrinsic DNA binding, which contributed to its H3K9me-independent nucleosome binding. CK2-mediated phosphorylation had a similar effect on the nucleosome-binding specificity of fly HP1a and S. pombe Swi6. These results suggested that HP1 phosphorylation has an evolutionarily conserved role in HP1's recognition of H3K9me-marked nucleosomes.

We report the preliminary results of the development of density-functional tight-binding (DFTB) repulsive potentials for the Zr – O element pair. The repulsive potentials were created using a computer code based on the particle swarm optimization. The potentials were tested on a set of systems for high temperature phases of bulk ZrO2, namely cubic and tetragonal. The potential sets were primarily developed for simulation of zirconia phase transitions at elevated temperatures.

Back electron transfer (BET) is one of the important processes that govern the decay of generated ion pairs in intermolecular photoinduced electron transfer reactions. Unfortunately, a detailed mechanism of BET reactions remains largely unknown in spite of their importance for the development of molecular photovoltaic structures. Here, we examine the BET reaction of pyrene (Py) and 1,4-dicyanobenzene (DCB) in acetonitrile (ACN) by using time-resolved near- and mid-IR spectroscopy. The Py dimer radical cation (Py2(·+)) and DCB radical anion (DCB(·-)) generated after photoexcitation of Py show asynchronous decay kinetics. To account for this observation, we propose a reaction mechanism that involves electron transfer from DCB(·-) to the solvent and charge recombination between the resulting ACN dimer anion and Py2(·+). The unique role of ACN as a charge mediator revealed herein could have implications for strategies that retard charge recombination in dye-sensitized solar cells.

The SN1-type hydrolysis reaction of cellobiose in ionic liquids (ILs) was theoretically investigated. First principles and ab initio quantum chemical methods were used in conjunction with the ‘reference interaction site model self-consistent field with spatial electron density distribution’ (RISM-SCF-SEDD) method. Reaction mechanism pathways are discussed and compared to calculations in gas phase and in aqueous solution. Analysis of solvation effects indicates strong interaction between hydrogen atoms of glucose hydroxyl groups and the anions in ILs, contributing to large stabilization of the reaction product. The calculated activation energy in ILs (24.5 kcal/mol) agrees quantitatively with the experimental value (26.5 kcal/mol).

Telomeric DNAs consist of tandem repeats of G-clusters such as TTAGGG and TG1-3, which are the human and yeast repeat sequences, respectively. In the yeast Saccharomyces cerevisiae, the telomeric repeats are non-nucleosomal, whereas in humans, they are organized in tightly packaged nucleosomes. However, previous in vitro studies revealed that the binding affinities of human and yeast telomeric repeat sequences to histone octamers in vitro were similar, which is apparently inconsistent with the differences in the human and yeast telomeric chromatin structures. To further investigate the relationship between telomeric sequences and chromatin structure, we examined the effect of telomeric repeats on the formation of positioned nucleosomes in vivo by indirect end-label mapping, primer extension mapping and nucleosome repeat analyses, using a defined minichromosome in yeast cells. We found that the human and yeast telomeric repeat sequences both disfavour nucleosome assembly and alter nucleosome positioning in the yeast minichromosome. We further demonstrated that the G-clusters in the telomeric repeats are required for the nucleosome-disfavouring properties. Thus, our results suggest that this inherent structural feature of the telomeric repeat sequences is involved in the functional dynamics of the telomeric chromatin structure.

The minimum structural unit of chromatin is the nucleosome core particle (NCP), consisting of 146 bp of DNA wrapped around a histone octamer, which itself contains two H2A/H2B dimers and one (H3/H4)2 tetramer. These multimers possess functionally important tail regions that are intrinsically disordered. In order to elucidate the mechanisms behind NCP assembly and disassembly processes, which are highly related to gene expression, structural characterization of the H2A/H2B dimer and (H3/H4)2 tetramer will be of importance. In the present study, human histone multimers with disordered tail regions were characterized by electrospray ionization (ESI) ion mobility-mass spectrometry (IM-MS) and molecular dynamics (MD) simulation. Experimentally obtained arrival times of these histone multimer ions showed rather wide distributions, implying that multiple conformers exist for each histone multimer in the gas phase. To examine their structures, MD simulations of the histone multimers were performed first in solution and then in vacuo at four temperatures, resulting in a variety of histone multimer structures. Theoretical collision cross-section (CCS) values calculated for the simulated structures revealed that structural models with smaller CCS values had more compact tail regions than those with larger CCS values. This implied that variation of the CCS values of the histone multimers were primarily due to the random behaviors of the tail regions in the gas phase. The combination of IM-MS and MD simulation enabled clear and comprehensive characterization of the gas-phase structures of histone multimers containing disordered tails.

Carbon nanotubes have long been described as rolled-up graphene sheets. It is only fairly recently observed that longitudinal cleavage of carbon nanotubes, using chemical, catalytical and electrical approaches, unzips them into thin graphene strips of various widths, the so-called graphene nanoribbons. In contrast, rolling up these flimsy ribbons into tubes in a real experiment has not been possible. Theoretical studies conducted by Kit et al. recently demonstrated the tube formation through twisting of graphene nanoribbon, an idea very different from the rolling-up postulation. Here we report the first experimental evidence of a thermally induced self-intertwining of graphene nanoribbons for the preferential synthesis of (7, 2) and (8, 1) tubes within parent-tube templates. Through the tailoring of ribbon’s width and edge, the present finding adds a radically new aspect to the understanding of carbon nanotube formation, shedding much light on not only the future chirality tuning, but also contemporary nanomaterials engineering., Carbon nanotubes can be considered as rolled-up small sheets of graphene. Here Lim and colleagues demonstrate this process, by fabricating carbon nanotubes through a thermally induced process of self-intertwining of graphene nanoribbons.

During chromatin-regulated processes, the histone H2A-H2B heterodimer functions dynamically in and out of the nucleosome. Although detailed crystal structures of nucleosomes have been established, that of the isolated full-length H2A-H2B heterodimer has remained elusive. Here, we have determined the solution structure of human H2A-H2B by NMR coupled with CS-Rosetta. H2A and H2B each contain a histone fold, comprising four α-helices and two β-strands (α1–β1–α2–β2–α3–αC), together with the long disordered N- and C-terminal H2A tails and the long N-terminal H2B tail. The N-terminal αN helix, C-terminal β3 strand and 310 helix of H2A observed in the H2A-H2B nucleosome structure are disordered in isolated H2A-H2B. In addition, the H2A α1 and H2B αC helices are not well fixed in the heterodimer and the H2A and H2B tails are not completely random coils. Comparison of hydrogen-deuterium exchange, fast hydrogen exchange and {1H}-15N hetero-nuclear NOE data with the CS-Rosetta structure indicates that there is some conformation in the H2A 310 helical and H2B Lys11 regions, while the repression domain of H2B (residues 27–34) exhibits an extended string-like structure. This first structure of the isolated H2A-H2B heterodimer provides insight into its dynamic functions in chromatin.

The linear-scaling divide-and-conquer (DC) quantum chemical methodology is applied to the density-functional tight-binding (DFTB) theory to develop a massively parallel program that achieves on-the-fly molecular reaction dynamics simulations of huge systems from scratch. The functions to perform large scale geometry optimization and molecular dynamics with DC-DFTB potential energy surface are implemented to the program called DC-DFTB-K. A novel interpolation-based algorithm is developed for parallelizing the determination of the Fermi level in the DC method. The performance of the DC-DFTB-K program is assessed using a laboratory computer and the K computer. Numerical tests show the high efficiency of the DC-DFTB-K program, a single-point energy gradient calculation of a one-million-atom system is completed within 60 s using 7290 nodes of the K computer. © 2016 Wiley Periodicals, Inc.

The chromodomain of HP1α binds directly to lysine 9-methylated histone H3 (H3K9me). This interaction is enhanced by phosphorylation of serine residues in the N-terminal tail of HP1α by unknown mechanism. Here we show that phosphorylation modulates flexibility of HP1α’s N-terminal tail, which strengthens the interaction with H3. NMR analysis of HP1α’s chromodomain with N-terminal tail reveals that phosphorylation does not change the overall tertiary structure, but apparently reduces the tail dynamics. Small angle X-ray scattering confirms that phosphorylation contributes to extending HP1α’s N-terminal tail. Systematic analysis using deletion mutants and replica exchange molecular dynamics simulations indicate that the phosphorylated serines and following acidic segment behave like an extended string and dynamically bind to H3 basic residues; without phosphorylation, the most N-terminal basic segment of HP1α inhibits interaction of the acidic segment with H3. Thus, the dynamic string-like behavior of HP1α’s N-terminal tail underlies the enhancement in H3 binding due to phosphorylation.

We present an efficient numerical integration scheme (TWOCENT) to be used in the context of automatized parameterization of the density-functional tight-binding (DFTB) method. The accuracy of the integration process is assessed and its range of applicability is discussed. The functionality of the developed code is tested by reproducing the electronic portion of the existing mio parameter sets and by reproducing a series of reference DFT band structures of elemental solids.