RNAPII driven post-translational modifications of nucleosomal histones.

The current understanding of how specific distributions of histone post-translational modifications (PTMs) are achieved throughout the chromatin remains incomplete. This review focuses on the role of RNA polymerase II (RNAPII) in establishing H2BK120/K123 ubiquitination and H3K4/K36 methylation distribution. The rate of RNAPII transcription is mainly a function of the RNAPII elongation and recruitment rates. Two major mechanisms link RNAPII's transcription rate to the distribution of PTMs. First, the phosphorylation patterns of Ser2P/Ser5P in the C-terminal domain of RNAPII change as a function of time, since the start of elongation, linking them to the elongation rate. Ser2P/Ser5P recruits specific histone PTM enzymes/activators to the nucleosome. Second, multiple rounds of binding and catalysis by the enzymes are required to establish higher methylations (H3K4/36me3). Thus, methylation states are determined by the transcription rate. In summary, the first mechanism determines the location of methylations in the gene, while the second mechanism determines the methylation state.

Torsional stress can regulate the unwrapping of two outer half superhelical turns of nucleosomal DNA.

Torsional stress has a significant impact on the structure and stability of the nucleosome. RNA polymerase imposes torsional stress on the DNA in chromatin and unwraps the DNA from the nucleosome to access the genetic information encoded in the DNA. To understand how the torsional stress affects the stability of the nucleosome, we examined the unwrapping of two half superhelical turns of nucleosomal DNA from either end of the DNA under torsional stress with all-atom molecular dynamics simulations. The free energies for unwrapping the DNA indicate that positive stress that overtwists DNA facilitates a large-scale asymmetric unwrapping of the DNA without a large extension of the DNA. During the unwrapping, one end of the DNA was dissociated from H3 and H2A-H2B, while the other end of the DNA stably remained wrapped. The detailed analysis indicates that this asymmetric dissociation is facilitated by the geometry and bendability of the DNA under positive stress. The geometry stabilized the interaction between the major groove of the twisted DNA and the H3 αN-helix, and the straightened DNA destabilized the interaction with H2A-H2B. Under negative stress, the DNA became more bendable and flexible, which facilitated the binding of the unwrapped DNA to the octamer in a stable state. Consequently, we conclude that the torsional stress has a significant impact on the affinity of the DNA and the octamer through the inherent nature of the DNA and can change the accessibility of regulatory proteins.

Extended ensemble simulations of a SARS-CoV-2 nsp1-5'-UTR complex.

Nonstructural protein 1 (nsp1) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a 180-residue protein that blocks translation of host mRNAs in SARS-CoV-2-infected cells. Although it is known that SARS-CoV-2's own RNA evades nsp1's host translation shutoff, the molecular mechanism underlying the evasion was poorly understood. We performed an extended ensemble molecular dynamics simulation to investigate the mechanism of the viral RNA evasion. Simulation results suggested that the stem loop structure of the SARS-CoV-2 RNA 5'-untranslated region (SL1) binds to both nsp1's N-terminal globular region and intrinsically disordered region. The consistency of the results was assessed by modeling nsp1-40S ribosome structure based on reported nsp1 experiments, including the X-ray crystallographic structure analysis, the cryo-EM electron density map, and cross-linking experiments. The SL1 binding region predicted from the simulation was open to the solvent, yet the ribosome could interact with SL1. Cluster analysis of the binding mode and detailed analysis of the binding poses suggest residues Arg124, Lys47, Arg43, and Asn126 may be involved in the SL1 recognition mechanism, consistent with the existing mutational analysis.

NRF2 DLG Domain Mutations Identified in Japanese Liver Cancer Patients Affect the Transcriptional Activity in HCC Cell Lines.

Geographically, East Asia had the highest liver cancer burden in 2017. Besides this, liver cancer-related deaths were high in Japan, accounting for 3.90% of total deaths. The development of liver cancer is influenced by several factors, and genetic alteration is one of the critical factors among them. Therefore, the detailed mechanism driving the oncogenic transformation of liver cells needs to be elucidated. Recently, many researchers have focused on investigating the liver cancer genome and identified somatic mutations (MTs) of several transcription factors. In this line, next-generation sequencing of the cancer genome identified that oxidative stress-related transcription factor NRF2 (NFE2L2) is mutated in different cancers, including hepatocellular carcinoma (HCC). Here, we demonstrated that NRF2 DLG motif mutations (NRF2 D29A and L30F), found in Japanese liver cancer patients, upregulate the transcriptional activity of NRF2 in HCC cell lines. Moreover, the transcriptional activity of NRF2 mutations is not suppressed by KEAP1, presumably because NRF2 MTs disturb proper NRF2-KEAP1 binding and block KEAP1-mediated degradation of NRF2. Additionally, we showed that both MTs upregulate the transcriptional activity of NRF2 on the MMP9 promoter in Hepa1-6 and Huh7 cells, suggesting that MT derived gain-of-function of NRF2 may be important for liver tumor progression. We also found that ectopic overexpression of oncogenic BRAF WT and V600E increases the transcriptional activity of NRF2 WT on both the 3xARE reporter and MMP9 promoter. Interestingly, NRF2 D29A and L30F MTs with oncogenic BRAF V600E MT synergistically upregulate the transcription activity of NRF2 on the 3xARE reporter and MMP9 promoter in Hepa1-6 and Huh7 cells. In summary, our findings suggest that MTs in NRF2 have pathogenic effects, and that NRF2 MTs together with oncogenic BRAF V600E MT synergistically cause more aberrant transcriptional activity. The high activity of NRF2 MTs in HCC with BRAF MT warrants further exploration of the potential diagnostic, prognostic, and therapeutic utility of this pathway in HCC.

Structural Studies of Overlapping Dinucleosomes in Solution.

An overlapping dinucleosome (OLDN) is a structure composed of one hexasome and one octasome and appears to be formed through nucleosome collision promoted by nucleosome remodeling factor(s). In this study, the solution structure of the OLDN was investigated through the integration of small-angle x-ray and neutron scattering (SAXS and SANS, respectively), computer modeling, and molecular dynamics simulations. Starting from the crystal structure, we generated a conformational ensemble based on normal mode analysis and searched for the conformations that reproduced the SAXS and SANS scattering curves well. We found that inclusion of histone tails, which are not observed in the crystal structure, greatly improved model quality. The obtained structural models suggest that OLDNs adopt a variety of conformations stabilized by histone tails situated at the interface between the hexasome and octasome, simultaneously binding to both the hexasomal and octasomal DNA. In addition, our models define a possible direction for the conformational changes or dynamics, which may provide important information that furthers our understanding of the role of chromatin dynamics in gene regulation.

Correction: Free energy profiles for unwrapping the outer superhelical turn of nucleosomal DNA.

[This corrects the article DOI: 10.1371/journal.pcbi.1006024.].

MNase, as a probe to study the sequence-dependent site exposures in the +1 nucleosomes of yeast.

The first nucleosomes in the downstream of transcription starting sites are called +1 nucleosomes, which are expected to be readily unwrapped for DNA transcription. To investigate DNA accessibility in +1 nucleosomes, MNase-seq experiments were carried out with 20 reconstituted +1 nucleosomes of budding yeast. Although MNase has been known for its sequence preference in DNA digestions, we confirmed that this sequence preference is overwhelmed by DNA accessibility by identifying the sequence-driven and accessibility-driven cleavages. Specifically, we find that sequences favoured by MNase at the end regions such as TA dinucleotide are prohibited from cleavage at the internal sites in the early stage of digestion. Nevertheless, sequences less favoured by MNase at the end regions such as AA/TT dinucleotide are predominantly cleaved at the internal sites in the early stage of digestion. Since AA/TT is known as a rigid dinucleotide step resistant to DNA bending, these internal cleavages reflect the local site exposures induced by DNA mechanics. As the DNA entry site of +1 nucleosomes in yeast is found AA/TT-rich, this sequence element may play a role in gene activation by reducing DNA-histone affinities along the direction of DNA transcription.

Free energy profiles for unwrapping the outer superhelical turn of nucleosomal DNA.

The eukaryotic genome is packaged into a nucleus in the form of chromatin. The fundamental structural unit of chromatin is a protein-DNA complex, the nucleosome, where 146 or 147 base pairs of DNA wrap 1.75 times around a histone core. To function in cellular processes, however, nucleosomal DNA must be unwrapped. Although this unwrapping has been experimentally investigated, details of the process at an atomic level are not yet well understood. Here, we used molecular dynamics simulation with an enhanced sampling method to calculate the free energy profiles for unwrapping the outer superhelical turn of nucleosomal DNA. A free energy change of about 11.5 kcal/mol for the unwrapping agrees well with values obtained in single molecule experiments. This simulation revealed a variety of conformational states, indicating there are many potential paths to outer superhelicdal turn unwrapping, but the dominant path is likely asymmetric. At one end of the DNA, the first five bps unwrap, after which a second five bps unwrap at the same end with no increase in free energy. The unwrapping then starts at the other end of the DNA, where 10 bps are unwrapped. During further unwrapping of 15 bps, the unwrapping advances at one of the ends, after which the other end of the DNA unwraps to complete the unwrapping of the outer superhelical turn. These results provide insight into the construction, disruption, and repositioning of nucleosomes, which are continuously ongoing during cellular processes.

H4 Tails Potentially Produce the Diversity in the Orientation of Two Nucleosomes.

Histone tails play an important role in internucleosomal interaction and chromatin compaction. To understand how the H4 tails are involved in the internucleosomal interaction, an adaptively biased molecular dynamics simulation of 63 models of two stacked nucleosomes, each with the H4 tails in different locations, was carried out. This simulation generated a variety of orientations of the separated nucleosomes depending on the formation of the H4 tail bridge between the H4 tails and the DNA of the neighboring nucleosomes. For the models that showed distinctive orientations of the two nucleosomes, the free energies of the separation of the nucleosomes were further investigated using umbrella sampling simulations. The attractive force between the nucleosomes was estimated from the free energies; the force when two H4 tail bridges formed varied from 36 to 63 pN, depending on the formation of the H4 tail-bridge and the interfacial interaction, whereas the force reduced to 15-18 pN after either one of the H4 tail bridges had broken, regardless of the conformation of the H4 tail. Additional simulations of the nucleosomes show that when the H4 tail was truncated, the force between the nucleosomes became repulsive (from-3 to -7 pN). We concluded that the H4 tails potentially produce the diversity in the orientation of the two nucleosomes, which would contribute to the polymorphism of the chromatin structure.

H3 Histone Tail Conformation within the Nucleosome and the Impact of K14 Acetylation Studied Using Enhanced Sampling Simulation.

Acetylation of lysine residues in histone tails is associated with gene transcription. Because histone tails are structurally flexible and intrinsically disordered, it is difficult to experimentally determine the tail conformations and the impact of acetylation. In this work, we performed simulations to sample H3 tail conformations with and without acetylation. The results show that irrespective of the presence or absence of the acetylation, the H3 tail remains in contact with the DNA and assumes an α-helix structure in some regions. Acetylation slightly weakened the interaction between the tail and DNA and enhanced α-helix formation, resulting in a more compact tail conformation. We inferred that this compaction induces unwrapping and exposure of the linker DNA, enabling DNA-binding proteins (e.g., transcription factors) to bind to their target sequences. In addition, our simulation also showed that acetylated lysine was more often exposed to the solvent, which is consistent with the fact that acetylation functions as a post-translational modification recognition site marker.

Spotting the difference in molecular dynamics simulations of biomolecules.

Comparing two trajectories from molecular simulations conducted under different conditions is not a trivial task. In this study, we apply a method called Linear Discriminant Analysis with ITERative procedure (LDA-ITER) to compare two molecular simulation results by finding the appropriate projection vectors. Because LDA-ITER attempts to determine a projection such that the projections of the two trajectories do not overlap, the comparison does not suffer from a strong anisotropy, which is an issue in protein dynamics. LDA-ITER is applied to two test cases: the T4 lysozyme protein simulation with or without a point mutation and the allosteric protein PDZ2 domain of hPTP1E with or without a ligand. The projection determined by the method agrees with the experimental data and previous simulations. The proposed procedure, which complements existing methods, is a versatile analytical method that is specialized to find the "difference" between two trajectories.

Adaptive lambda square dynamics simulation: an efficient conformational sampling method for biomolecules.

A novel, efficient sampling method for biomolecules is proposed. The partial multicanonical molecular dynamics (McMD) was recently developed as a method that improved generalized ensemble (GE) methods to focus sampling only on a part of a system (GEPS); however, it was not tested well. We found that partial McMD did not work well for polylysine decapeptide and gave significantly worse sampling efficiency than a conventional GE. Herein, we elucidate the fundamental reason for this and propose a novel GEPS, adaptive lambda square dynamics (ALSD), which can resolve the problem faced when using partial McMD. We demonstrate that ALSD greatly increases the sampling efficiency over a conventional GE. We believe that ALSD is an effective method and is applicable to the conformational sampling of larger and more complicated biomolecule systems.

A new method for evaluating the specificity of indirect readout in protein-DNA recognition.

Proteins recognize a specific DNA sequence not only through direct contact (direct readout) with base pairs but also through sequence-dependent conformation and/or flexibility of DNA (indirect readout). However, it is difficult to assess the contribution of indirect readout to the sequence specificity. What is needed is a straightforward method for quantifying its contributions to specificity. Using Bayesian statistics, we derived the probability of a particular sequence for a given DNA structure from the trajectories of molecular dynamics (MD) simulations of DNAs containing all possible tetramer sequences. Then, we quantified the specificity of indirect readout based on the information entropy associated with the probability. We tested this method with known structures of protein-DNA complexes. This method enabled us to correctly predict those regions where experiments suggested the involvement of indirect readout. The results also indicated new regions where the indirect readout mechanism makes major contributions to the recognition. The present method can be used to estimate the contribution of indirect readout without approximations to the distributions in the conformational ensembles of DNA, and would serve as a powerful tool to study the mechanism of protein-DNA recognition.

Benchmarking of force fields to characterize the intrinsically disordered R2-FUS-LC region.

Intrinsically Disordered Proteins (IDPs) play crucial roles in numerous diseases like Alzheimer's and ALS by forming irreversible amyloid fibrils. The effectiveness of force fields (FFs) developed for globular proteins and their modified versions for IDPs varies depending on the specific protein. This study assesses 13 FFs, including AMBER and CHARMM, by simulating the R2 region of the FUS-LC domain (R2-FUS-LC region), an IDP implicated in ALS. Due to the flexibility of the region, we show that utilizing multiple measures, which evaluate the local and global conformations, and combining them together into a final score are important for a comprehensive evaluation of force fields. The results suggest c36m2021s3p with mTIP3p water model is the most balanced FF, capable of generating various conformations compatible with known ones. In addition, the mTIP3P water model is computationally more efficient than those of top-ranked AMBER FFs with four-site water models. The evaluation also reveals that AMBER FFs tend to generate more compact conformations compared to CHARMM FFs but also more non-native contacts. The top-ranking AMBER and CHARMM FFs can reproduce intra-peptide contacts but underperform for inter-peptide contacts, indicating there is room for improvement.

Structural and Dynamic Changes of Nucleosome upon GATA3 Binding.

Pioneer factors, which can directly bind to nucleosomes, have been considered to change chromatin conformations. However, the binding impact on the nucleosome is little known. Here, we show how the pioneer factor GATA3 binds to nucleosomal DNA and affects the conformation and dynamics of nucleosomes by using a combination of SAXS, molecular modeling, and molecular dynamics simulations. Our structural models, consistent with the SAXS data, indicate that only one of the two DNA binding domains, N- and C-fingers, of GATA3 binds to an end of the DNA in solution. Our MD simulations further showed that the other unbound end of the DNA increases the fluctuation and enhances the DNA dissociation from the histone core when the N-finger binds to a DNA end, a site near the entry or exit of the nucleosome. However, this was not true for the binding of the C-finger that binds to a location about 15 base pairs distant from the DNA end. In this case, DNA dissociation occurred on the bound end. Taken together, we suggest that the N-finger and C-finger bindings of GATA3 commonly enhance DNA dissociation at one of the two DNA ends (the bound end for the C-finger binding and the unbound end for the N-finger binding), leading to triggering a conformational change in the chromatin.

Designer Adaptor Proteins for Functional Conversion of Peptides to Small-Molecule Ligands toward In-Cell Catalytic Protein Modification.

Peptides are privileged ligands for diverse biomacromolecules, including proteins; however, their utility is often limited due to low membrane permeability and in-cell instability. Here, we report peptide ligand-inserted eDHFR (PLIED) fusion protein as a universal adaptor for targeting proteins of interest (POI) with cell-permeable and stable synthetic functional small molecules (SFSM). PLIED binds to POI through the peptide moiety, properly orienting its eDHFR moiety, which then recruits trimethoprim (TMP)-conjugated SFSM to POI. Using a lysine-acylating BAHA catalyst as SFSM, we demonstrate that POI (MDM2 and chromatin histone) are post-translationally and synthetically acetylated at specific lysine residues. The residue-selectivity is predictable in an atomic resolution from molecular dynamics simulations of the POI/PLIED/TMP-BAHA (MTX was used as a TMP model) ternary complex. This designer adaptor approach universally enables functional conversion of impermeable peptide ligands to permeable small-molecule ligands, thus expanding the in-cell toolbox of chemical biology.

Molecular structure of isolated MvspI, a variable surface protein of the fish pathogen Mycoplasma mobile.

Mycoplasma mobile is a parasitic bacterium that causes necrosis in the gills of freshwater fishes. This study examines the molecular structure of its variable surface protein, MvspI, whose open reading frame encodes 2,002 amino acids. MvspI was isolated from mycoplasma cells by a biochemical procedure to 92% homogeneity. Gel filtration and analytical ultracentrifugation suggested that this protein is a cylinder-shaped monomer with axes of 66 and 2.7 nm. Rotary shadowing transmission electron microscopy of MvspI showed that the molecule is composed of two rods 30 and 45 nm long; the latter rod occasionally features a bulge. Immuno-electron microscopy and epitope mapping showed that the bulge end of the molecular image corresponds to the C terminus of the amino acid sequence. Partial digestion by various proteases suggested that the N-terminal part, comprised of 697 amino acids, is flexible. Analysis of the predicted amino acid sequence showed that the molecule features a lipoprotein and 16 repeats of about 90 residues; 15 positions exist between residues 88 and 1479, and the other position is between residues 1725 and 1807. The amino acid sequence of MvspI was mapped onto a molecular image obtained by electron microscopy. The present study is the first to elucidate the molecular shape of a variable surface protein of mycoplasma.

Phase retrieval from single biomolecule diffraction pattern.

In this paper, we propose the SPR (sparse phase retrieval) method, which is a new phase retrieval method for coherent x-ray diffraction imaging (CXDI). Conventional phase retrieval methods effectively solve the problem for high signal-to-noise ratio measurements, but would not be sufficient for single biomolecular imaging which is expected to be realized with femto-second x-ray free electron laser pulses. The SPR method is based on the Bayesian statistics. It does not need to set the object boundary constraint that is required by the commonly used hybrid input-output (HIO) method, instead a prior distribution is defined with an exponential distribution and used for the estimation. Simulation results demonstrate that the proposed method reconstructs the electron density under a noisy condition even some central pixels are masked.

Free Energy Landscape of H2A-H2B Displacement From Nucleosome.

Nucleosome reconstitution plays an important role in many cellular functions. As an initial step, H2A-H2B dimer displacement, which is accompanied by disruption of many of the interactions within the nucleosome, should occur. To understand how H2A-H2B dimer displacement occurs, an adaptively biased molecular dynamics (ABMD) simulation was carried out to generate a variety of displacements of the H2A-H2B dimer from the fully wrapped to partially unwrapped nucleosome structures. With regards to these structures, the free energy landscape of the dimer displacement was investigated using umbrella sampling simulations. We found that the main contributors to the free energy were the docking domain of H2A and the C-terminal of H4. There were various paths for the dimer displacement which were dependent on the extent of nucleosomal DNA wrapping, suggesting that modulation of the intra-nucleosomal interaction by external factors such as histone chaperones could control the path for the H2A-H2B dimer displacement. Key residues which contributed to the free energy have also been reported to be involved in the mutations and posttranslational modifications (PTMs) which are important for assembling and/or reassembling the nucleosome at the molecular level and are found in cancer cells at the phenotypic level. Our results give insight into how the H2A-H2B dimer displacement proceeds along various paths according to different interactions within the nucleosome.

HNF1A POU Domain Mutations Found in Japanese Liver Cancer Patients Cause Downregulation of HNF4A Promoter Activity with Possible Disruption in Transcription Networks.

Hepatocyte nuclear factor 1A (HNF1A) is the master regulator of liver homeostasis and organogenesis and regulates many aspects of hepatocyte functions. It acts as a tumor suppressor in the liver, evidenced by the increased proliferation in HNF1A knockout (KO) hepatocytes. Hence, we postulated that any loss-of-function variation in the gene structure or composition (mutation) could trigger dysfunction, including disrupted transcriptional networks in liver cells. From the International Cancer Genome Consortium (ICGC) database of cancer genomes, we identified several HNF1A mutations located in the functional Pit-Oct-Unc (POU) domain. In our biochemical analysis, we found that the HNF1A POU-domain mutations Y122C, R229Q and V259F suppressed HNF4A promoter activity and disrupted the binding of HNF1A to its target HNF4A promoter without any effect on the nuclear localization. Our results suggest that the decreased transcriptional activity of HNF1A mutants is due to impaired DNA binding. Through structural simulation analysis, we found that a V259F mutation was likely to affect DNA interaction by inducing large conformational changes in the N-terminal region of HNF1A. The results suggest that POU-domain mutations of HNF1A downregulate HNF4A gene expression. Therefore, to mimic the HNF1A mutation phenotype in transcription networks, we performed siRNA-mediated knockdown (KD) of HNF4A. Through RNA-Seq data analysis for the HNF4A KD, we found 748 differentially expressed genes (DEGs), of which 311 genes were downregulated (e.g., HNF1A, ApoB and SOAT2) and 437 genes were upregulated. Kyoto Encyclopedia of Genes and Genomes (KEGG) mapping revealed that the DEGs were involved in several signaling pathways (e.g., lipid and cholesterol metabolic pathways). Protein-protein network analysis suggested that the downregulated genes were related to lipid and cholesterol metabolism pathways, which are implicated in hepatocellular carcinoma (HCC) development. Our study demonstrates that mutations of HNF1A in the POU domain result in the downregulation of HNF1A target genes, including HNF4A, and this may trigger HCC development through the disruption of HNF4A-HNF1A transcriptional networks.