Three human RNA polymerases interact with TFIIH via a common RPB6 subunit.

In eukaryotes, three RNA polymerases (RNAPs) play essential roles in the synthesis of various types of RNA: namely, RNAPI for rRNA; RNAPII for mRNA and most snRNAs; and RNAPIII for tRNA and other small RNAs. All three RNAPs possess a short flexible tail derived from their common subunit RPB6. However, the function of this shared N-terminal tail (NTT) is not clear. Here we show that NTT interacts with the PH domain (PH-D) of the p62 subunit of the general transcription/repair factor TFIIH, and present the structures of RPB6 unbound and bound to PH-D by nuclear magnetic resonance (NMR). Using available cryo-EM structures, we modelled the activated elongation complex of RNAPII bound to TFIIH. We also provide evidence that the recruitment of TFIIH to transcription sites through the p62-RPB6 interaction is a common mechanism for transcription-coupled nucleotide excision repair (TC-NER) of RNAPI- and RNAPII-transcribed genes. Moreover, point mutations in the RPB6 NTT cause a significant reduction in transcription of RNAPI-, RNAPII- and RNAPIII-transcribed genes. These and other results show that the p62-RPB6 interaction plays multiple roles in transcription, TC-NER, and cell proliferation, suggesting that TFIIH is engaged in all RNAP systems.

Divide-and-Conquer Linear-Scaling Quantum Chemical Computations.

Fragmentation and embedding schemes are of great importance when applying quantum-chemical calculations to more complex and attractive targets. The divide-and-conquer (DC)-based quantum-chemical model is a fragmentation scheme that can be connected to embedding schemes. This feature article explains several DC-based schemes developed by the authors over the last two decades, which was inspired by the pioneering study of DC self-consistent field (SCF) method by Yang and Lee (J. Chem. Phys. 1995, 103, 5674-5678). First, the theoretical aspects of the DC-based SCF, electron correlation, excited-state, and nuclear orbital methods are described, followed by the two-component relativistic theory, quantum-mechanical molecular dynamics simulation, and the introduction of three programs, including DC-based schemes. Illustrative applications confirmed the accuracy and feasibility of the DC-based schemes.

Species-selective nanoreactor molecular dynamics simulations based on linear-scaling tight-binding quantum chemical calculations.

Here, extensions to quantum chemical nanoreactor molecular dynamics simulations for discovering complex reactive events are presented. The species-selective algorithm, where the nanoreactor effectively works for the selected desired reactants, was introduced to the original scheme. Moreover, for efficient simulations of large model systems with the modified approach, the divide-and-conquer linear-scaling density functional tight-binding method was exploited. Two illustrative applications of the polymerization of propylene and cyclopropane mixtures and the aggregation of sodium chloride from aqueous solutions indicate that species-selective quantum chemical nanoreactor molecular dynamics is a promising method to accelerate the sampling of multicomponent chemical processes proceeding under relatively mild conditions.

Enabling large-scale quantum path integral molecular dynamics simulations through the integration of Dcdftbmd and i-PI codes.

A large-scale quantum chemical calculation program, Dcdftbmd, was integrated with a Python-based advanced atomistic simulation program, i-PI. The implementation of a client-server model enabled hierarchical parallelization with respect to replicas and force evaluations. The established framework demonstrated that quantum path integral molecular dynamics simulations can be executed with high efficiency for systems consisting of a few tens of replicas and containing thousands of atoms. The application of the framework to bulk water systems, with and without an excess proton, demonstrated that nuclear quantum effects are significant for intra- and inter-molecular structural properties, including oxygen-hydrogen bond distance and radial distribution function around the hydrated excess proton.

Structural and dynamical insights into the PH domain of p62 in human TFIIH.

TFIIH is a crucial transcription and DNA repair factor consisting of the seven-subunit core. The core subunit p62 contains a pleckstrin homology domain (PH-D), which is essential for locating TFIIH at transcription initiation and DNA damage sites, and two BSD (BTF2-like transcription factors, synapse-associated proteins and DOS2-like proteins) domains. A recent cryo-electron microscopy (cryo-EM) structure of human TFIIH visualized most parts of core, except for the PH-D. Here, by nuclear magnetic resonance spectroscopy we have established the solution structure of human p62 PH-D connected to the BSD1 domain by a highly flexible linker, suggesting the flexibility of PH-D in TFIIH. Based on this dynamic character, the PH-D was modeled in the cryo-EM structure to obtain the whole human TFIIH core structure, which indicates that the PH-D moves around the surface of core with a specific but limited spatial distribution; these dynamic structures were refined by molecular dynamics (MD) simulations. Furthermore, we built models, also refined by MD simulations, of TFIIH in complex with five p62-binding partners, including transcription factors TFIIEα, p53 and DP1, and nucleotide excision repair factors XPC and UVSSA. The models explain why the PH-D is crucially targeted by these factors, which use their intrinsically disordered acidic regions for TFIIH recruitment.

Acetylated histone H4 tail enhances histone H3 tail acetylation by altering their mutual dynamics in the nucleosome.

The structural unit of eukaryotic chromatin is a nucleosome, comprising two histone H2A-H2B heterodimers and one histone (H3-H4)2 tetramer, wrapped around by ∼146 bp of DNA. The N-terminal flexible histone tails stick out from the histone core and have extensive posttranslational modifications, causing epigenetic changes of chromatin. Although crystal and cryogenic electron microscopy structures of nucleosomes are available, the flexible tail structures remain elusive. Using NMR, we have examined the dynamics of histone H3 tails in nucleosomes containing unmodified and tetra-acetylated H4 tails. In unmodified nucleosome, the H3 tail adopts a dynamic equilibrium structure between DNA-contact and reduced-contact states. In acetylated H4 nucleosome, however, the H3 tail equilibrium shifts to a mainly DNA-contact state with a minor reduced-contact state. The acetylated H4 tail is dynamically released from its own DNA-contact state to a reduced-contact state, while the H3 tail DNA-contact state becomes major. Notably, H3 K14 in the acetylated H4 nucleosome is much more accessible to acetyltransferase Gcn5 relative to unmodified nucleosome, possibly due to the formation of a favorable H3 tail conformation for Gcn5. In summary, each histone tail adopts a characteristic dynamic state but regulates one other, probably creating a histone tail network even on a nucleosome.

Characteristic computed tomography features in mesenchymal-epithelial transition exon14 skipping-positive non-small cell lung cancer.

BACKGROUND: Mesenchymal-epithelial transition exon14 (METex14) skipping is one of the therapeutic driver oncogene mutations in non-small cell lung cancer (NSCLC), and can be treated with tepotinib and capmatinib. There is only one report on computed tomography (CT) findings of METex14 skipping-positive NSCLC, which shows that the primary tumor tends to have a large mass in the upper lobe, and extrathoracic metastases are common. This study examined the CT findings of METex14 skipping-positive NSCLC, focusing on the features of the margins and internal structures. METHODS: We consecutively included patients with METex14 skipping-positive NSCLC who were diagnosed between January 2018 and December 2020 at four independent institutions. We retrospectively reviewed the patient demographics and CT findings for tumor margins (invasion into surrounding tissue, lobulation, pleural indentation, spicula, and ground-glass opacity) and internal structures (air bronchograms, cavitation and internal low-density area). RESULTS: Fifteen patients with METex14 skipping-positive NSCLC were identified. Almost half of the patients were men (7/15; 46.7%), and their median age was 75.0 years. More than half were either current or former smokers (9/15; 60.0%). A vast majority of histological subtypes were adenocarcinoma (10/15; 66.7%), followed by pleomorphic carcinoma (3/15; 20.0%) and squamous cell carcinoma (2/15; 13.3%). With regard to CT findings, most primary tumors presented as masses larger than 30 mm (12/15; 80.0%) and were located in the upper lobes (12/15; 80.0%). Invasion into surrounding tissue and presence of internal low-density areas were observed in 60.0% (9/15) and 66.7% (10/15) of the primary tumors, respectively. Additionally, their frequencies increased to 72.7% (8/11) and 90.9% (10/11) in stage III/IV cases, respectively. In lymph node metastasis, internal low-density areas were observed in 8/10 cases (80.0%). Although these two CT features were rarely observed in distant metastases at diagnosis, they became apparent with progression of the metastatic tumor size. CONCLUSIONS: METex14 skipping-positive NSCLC tumors tend to invade surrounding tissue and possess internal low-density areas. These CT findings might be characteristic of METex14 skipping-positive NSCLC.

Clinical significance of BIM deletion polymorphism in chemoradiotherapy for non-small cell lung cancer.

The standard treatment for locally advanced non-small cell lung cancer (NSCLC) is chemoradiotherapy (CRT) followed by anti-programmed cell death-ligand 1 (anti-PD-L1) treatment. BIM deletion polymorphism induces the suppression of apoptosis resulting from epidermal growth factor (EGFR)-tyrosine kinase inhibitors in EGFR-mutated NSCLC patients. We aimed to examine the effects of BIM polymorphism on CRT and anti-PD-L1/PD-1 treatment in NSCLC patients. In this retrospective study of 1312 patients with unresectable NSCLC treated at Higashi-Hiroshima Medical Center and Hiroshima University Hospital between April 1994 and October 2019, we enrolled those who underwent CRT or chemotherapy using carboplatin + paclitaxel or cisplatin + vinorelbine, or anti-PD-L1/PD-1 treatment. Of 1312 patients, 88, 80, and 74 underwent CRT, chemotherapy, and anti-PD-L1/PD-1 treatment, respectively, and 17.0%, 15.2% and 17.6% of these patients showed BIM polymorphism. Among patients receiving CRT, the progression-free survival was significantly shorter in those with BIM deletion than in those without. In the multivariate analyses, BIM polymorphism was an independent factor of poor anti-tumor effects. These results were not observed in the chemotherapy and anti-PD-L1/PD-1 treatment groups. In in vitro experiments, BIM expression suppression using small interfering RNA in NSCLC cell lines showed a significantly suppressed anti-tumor effect and apoptosis after irradiation but not chemotherapy. In conclusion, we showed that BIM polymorphism was a poor-predictive factor for anti-tumor effects in NSCLC patients who underwent CRT, specifically radiotherapy. In the implementation of CRT in patients with BIM polymorphism, we should consider subsequent treatment, keeping in mind that CRT may be insufficient.

Density-Functional Tight-Binding Parameters for Bulk Zirconium: A Case Study for Repulsive Potentials.

Density-functional tight-binding (DFTB) parameters are presented for the simulation of the bulk phases of zirconium. Electronic parameters were obtained using a band structure fitting strategy, while two-center repulsive potentials were created by particle swarm optimization. As objective functions for the repulsive potential fitting, we employed the Birch-Murnaghan equations of state for hexagonal close-packed (HCP), body-centered cubic (BCC) and ω phases of Zr from density-functional theory (DFT). When fractional atomic coordinates are not allowed to change in the generation of the equation-of-state curves, long-range repulsive DFTB potentials are able to almost perfectly reproduce equilibrium structures, relative DFT energies of the bulk phases, and bulk moduli. However, the same potentials lead to artifacts in the DFTB potential energy surfaces when atom positions in the unit cell are allowed to fully relax during the change of unit cell parameters. Conventional short-range repulsive DFTB potentials, while inferior in their ability to reproduce DFT bulk energetics, are able to correctly reproduce the qualitative shape of the DFT potential energy surfaces, including the location of global minima, and can therefore be considered more transferable.

Hierarchical parallelization of divide-and-conquer density functional tight-binding molecular dynamics and metadynamics simulations.

Massively parallel divide-and-conquer density functional tight-binding (DC-DFTB) molecular dynamics and metadynamics simulations are efficient approaches for describing various chemical reactions and dynamic processes of large complex systems via quantum mechanics. In this study, DC-DFTB simulations were combined with multi-replica techniques. Specifically, multiple walkers metadynamics, replica exchange molecular dynamics, and parallel tempering metadynamics methods were implemented hierarchically into the in-house Dcdftbmd program. Test simulations in an aqueous phase of the internal rotation of formamide and conformational changes of dialanine showed that the newly developed extensions increase the sampling efficiency and the exploration capabilities in DC-DFTB configuration space.

Spin-flip approach within time-dependent density functional tight-binding method: Theory and applications.

A spin-flip time-dependent density functional tight-binding (SF-TDDFTB) method is developed that describes target states as spin-flipping excitation from a high-spin reference state obtained by the spin-restricted open shell treatment. Furthermore, the SF-TDDFTB formulation is extended to long-range correction (LC), denoted as SF-TDLCDFTB. The LC technique corrects the overdelocalization of electron density in systems such as charge-transfer systems, which is typically found in conventional DFTB calculations as well as density functional theory calculations using pure functionals. The numerical assessment of the SF-TDDFTB method shows smooth potential curves for the bond dissociation of hydrogen fluoride and the double-bond rotation of ethylene and the double-cone shape of H3 as the simplest degenerate systems. In addition, numerical assessments of SF-TDDFTB and SF-TDLCDFTB for 39 S0 /S1 minimum energy conical intersection (MECI) structures are performed. The SF-TDDFTB and SF-TDLCDFTB methods drastically reduce the computational cost with accuracy for MECI structures compared with SF-TDDFT.

Intrinsic nucleic acid-binding activity of Chp1 chromodomain is required for heterochromatic gene silencing.

Centromeric heterochromatin assembly in fission yeast requires the RNAi pathway. Chp1, a chromodomain (CD) protein, forms the Ago1-containing RNA-induced transcriptional silencing (RITS) complex and recruits siRNA-bound RITS to methylated histone H3 lysine 9 (H3K9me) via its CD. Here, we show that the CD of Chp1 (Chp1-CD) possesses unique nucleic acid-binding activities that are essential for heterochromatic gene silencing. Detailed electrophoretic-mobility shift analyses demonstrated that Chp1 binds to RNA via the CD in addition to its central RNA-recognition motif. Interestingly, robust RNA- and DNA-binding activity of Chp1-CD was strongly enhanced when it was bound to H3K9me, which was revealed to involve a positively charged domain within the Chp1-CD by structural analyses. These results demonstrate a role for the CD that provides a link between RNA, DNA, and methylated histone tails to ensure heterochromatic gene silencing.

Dcdftbmd: Divide-and-Conquer Density Functional Tight-Binding Program for Huge-System Quantum Mechanical Molecular Dynamics Simulations.

Dcdftbmd is a Fortran 90/95 program that enables efficient quantum mechanical molecular dynamics (MD) simulations using divide-and-conquer density functional tight-binding (DC-DFTB) method. Based on the remarkable performance of previous massively parallel DC-DFTB energy and gradient calculations for huge systems, the code has been specialized to MD simulations. Recent implementations and modifications including DFTB extensions, improved computational speed in the DC-DFTB computational steps, algorithms for efficient initial guess charge prediction, and free energy calculations via metadynamics technique have enhanced the capability to obtain atomistic insights in novel applications to nanomaterials and biomolecules. The energy, structure, and other molecular properties are also accessible through the single-point calculation, geometry optimization, and vibrational frequency analysis. The available functionalities are outlined together with efficiency tests and simulation examples. © 2019 Wiley Periodicals, Inc.

GPU-Accelerated Large-Scale Excited-State Simulation Based on Divide-and-Conquer Time-Dependent Density-Functional Tight-Binding.

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.

Development of Divide-and-Conquer Density-Functional Tight-Binding Method for Theoretical Research on Li-Ion Battery.

The density-functional tight-binding (DFTB) method is one of the useful quantum chemical methods, which provides a good balance between accuracy and computational efficiency. In this account, we reviewed the basis of the DFTB method, the linear-scaling divide-and-conquer (DC) technique, as well as the parameterization process. We also provide some refinement, modifications, and extension of the existing parameters that can be applicable for lithium-ion battery systems. The diffusion constants of common electrolyte molecules and LiTFSA salt in solution have been estimated using DC-DFTB molecular dynamics simulation with our new parameters. The resulting diffusion constants have good agreement to the experimental diffusion constants.

Common TFIIH recruitment mechanism in global genome and transcription-coupled repair subpathways.

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.

Sodium- and Potassium-Hydrate Melts Containing Asymmetric Imide Anions for High-Voltage Aqueous Batteries.

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.

Structural Diversity of Nucleosomes Characterized by Native Mass Spectrometry.

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.

Parallel implementation of efficient charge-charge interaction evaluation scheme in periodic divide-and-conquer density-functional tight-binding calculations.

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.

Density-Functional Tight-Binding Molecular Dynamics Simulations of Excess Proton Diffusion in Ice Ih, Ice Ic, Ice III, and Melted Ice VI Phases.

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, the energy barriers for high-density structures were several times larger than those of low-density systems. The simulation results were likely related to the suppression of proton transfer in disordered water configurations, in particular, ice with low hexagonal symmetry.