Structure and Dynamics of a 197 bp Nucleosome in Complex with Linker Histone H1.

Linker histones associate with nucleosomes to promote the formation of higher-order chromatin structure, but the underlying molecular details are unclear. We investigated the structure of a 197 bp nucleosome bearing symmetric 25 bp linker DNA arms in complex with vertebrate linker histone H1. We determined electron cryo-microscopy (cryo-EM) and crystal structures of unbound and H1-bound nucleosomes and validated these structures by site-directed protein cross-linking and hydroxyl radical footprinting experiments. Histone H1 shifts the conformational landscape of the nucleosome by drawing the two linkers together and reducing their flexibility. The H1 C-terminal domain (CTD) localizes primarily to a single linker, while the H1 globular domain contacts the nucleosome dyad and both linkers, associating more closely with the CTD-distal linker. These findings reveal that H1 imparts a strong degree of asymmetry to the nucleosome, which is likely to influence the assembly and architecture of higher-order structures.

Sequence-specific dynamics of DNA response elements and their flanking sites regulate the recognition by AP-1 transcription factors.

Activator proteins 1 (AP-1) comprise one of the largest families of eukaryotic basic leucine zipper transcription factors. Despite advances in the characterization of AP-1 DNA-binding sites, our ability to predict new binding sites and explain how the proteins achieve different gene expression levels remains limited. Here we address the role of sequence-specific DNA flexibility for stability and specific binding of AP-1 factors, using microsecond-long molecular dynamics simulations. As a model system, we employ yeast AP-1 factor Yap1 binding to three different response elements from two genetic environments. Our data show that Yap1 actively exploits the sequence-specific flexibility of DNA within the response element to form stable protein-DNA complexes. The stability also depends on the four to six flanking nucleotides, adjacent to the response elements. The flanking sequences modulate the conformational adaptability of the response element, making it more shape-efficient to form specific contacts with the protein. Bioinformatics analysis of differential expression of the studied genes supports our conclusions: the stability of Yap1-DNA complexes, modulated by the flanking environment, influences the gene expression levels. Our results provide new insights into mechanisms of protein-DNA recognition and the biological regulation of gene expression levels in eukaryotes.

Michler's hydrol blue elucidates structural differences in prion strains.

Yeast prions provide self-templating protein-based mechanisms of inheritance whose conformational changes lead to the acquisition of diverse new phenotypes. The best studied of these is the prion domain (NM) of Sup35, which forms an amyloid that can adopt several distinct conformations (strains) that confer distinct phenotypes when introduced into cells that do not carry the prion. Classic dyes, such as thioflavin T and Congo red, exhibit large increases in fluorescence when bound to amyloids, but these dyes are not sensitive to local structural differences that distinguish amyloid strains. Here we describe the use of Michler's hydrol blue (MHB) to investigate fibrils formed by the weak and strong prion fibrils of Sup35NM and find that MHB differentiates between these two polymorphs. Quantum mechanical time-dependent density functional theory (TDDFT) calculations indicate that the fluorescence properties of amyloid-bound MHB can be correlated to the change of binding site polarity and that a tyrosine to phenylalanine substitution at a binding site could be detected. Through the use of site-specific mutants, we demonstrate that MHB is a site-specific environmentally sensitive probe that can provide structural details about amyloid fibrils and their polymorphs.

The ancillary N-terminal region of the yeast AP-1 transcription factor Yap8 contributes to its DNA binding specificity.

Activator protein 1 (AP-1) is one of the largest families of basic leucine zipper (bZIP) transcription factors in eukaryotic cells. How AP-1 proteins achieve target DNA binding specificity remains elusive. In Saccharomyces cerevisiae, the AP-1-like protein (Yap) family comprises eight members (Yap1 to Yap8) that display distinct genomic target sites despite high sequence homology of their DNA binding bZIP domains. In contrast to the other members of the Yap family, which preferentially bind to short (7-8 bp) DNA motifs, Yap8 binds to an unusually long DNA motif (13 bp). It has been unclear what determines this unique specificity of Yap8. In this work, we use molecular and biochemical analyses combined with computer-based structural design and molecular dynamics simulations of Yap8-DNA interactions to better understand the structural basis of DNA binding specificity determinants. We identify specific residues in the N-terminal tail preceding the basic region, which define stable association of Yap8 with its target promoter. We propose that the N-terminal tail directly interacts with DNA and stabilizes Yap8 binding to the 13 bp motif. Thus, beside the core basic region, the adjacent N-terminal region contributes to alternative DNA binding selectivity within the AP-1 family.

Sequence-dependent response of DNA to torsional stress: a potential biological regulation mechanism.

Torsional restraints on DNA change in time and space during the life of the cell and are an integral part of processes such as gene expression, DNA repair and packaging. The mechanical behavior of DNA under torsional stress has been studied on a mesoscopic scale, but little is known concerning its response at the level of individual base pairs and the effects of base pair composition. To answer this question, we have developed a geometrical restraint that can accurately control the total twist of a DNA segment during all-atom molecular dynamics simulations. By applying this restraint to four different DNA oligomers, we are able to show that DNA responds to both under- and overtwisting in a very heterogeneous manner. Certain base pair steps, in specific sequence environments, are able to absorb most of the torsional stress, leaving other steps close to their relaxed conformation. This heterogeneity also affects the local torsional modulus of DNA. These findings suggest that modifying torsional stress on DNA could act as a modulator for protein binding via the heterogeneous changes in local DNA structure.

A stretched conformation of DNA with a biological role?

We have discovered a well-defined extended conformation of double-stranded DNA, which we call Σ-DNA, using laser-tweezers force-spectroscopy experiments. At a transition force corresponding to free energy change ΔG = 1·57 ± 0·12 kcal (mol base pair)-1 60 or 122 base-pair long synthetic GC-rich sequences, when pulled by the 3'-3' strands, undergo a sharp transition to the 1·52 ± 0·04 times longer Σ-DNA. Intriguingly, the same degree of extension is also found in DNA complexes with recombinase proteins, such as bacterial RecA and eukaryotic Rad51. Despite vital importance to all biological organisms for survival, genome maintenance and evolution, the recombination reaction is not yet understood at atomic level. We here propose that the structural distortion represented by Σ-DNA, which is thus physically inherent to the nucleic acid, is related to how recombination proteins mediate recognition of sequence homology and execute strand exchange. Our hypothesis is that a homogeneously stretched DNA undergoes a 'disproportionation' into an inhomogeneous Σ-form consisting of triplets of locally B-like perpendicularly stacked bases. This structure may ensure improved fidelity of base-pair recognition and promote rejection in case of mismatch during homologous recombination reaction. Because a triplet is the length of a gene codon, we speculate that the structural physics of nucleic acids may have biased the evolution of recombinase proteins to exploit triplet base stacks and also the genetic code.

Orientation of aromatic residues in amyloid cores: structural insights into prion fiber diversity.

Structural conversion of one given protein sequence into different amyloid states, resulting in distinct phenotypes, is one of the most intriguing phenomena of protein biology. Despite great efforts the structural origin of prion diversity remains elusive, mainly because amyloids are insoluble yet noncrystalline and therefore not easily amenable to traditional structural-biology methods. We investigate two different phenotypic prion strains, weak and strong, of yeast translation termination factor Sup35 with respect to angular orientation of tyrosines using polarized light spectroscopy. By applying a combination of alignment methods the degree of fiber orientation can be assessed, which allows a relatively accurate determination of the aromatic ring angles. Surprisingly, the strains show identical average orientations of the tyrosines, which are evenly spread through the amyloid core. Small variations between the two strains are related to the local environment of a fraction of tyrosines outside the core, potentially reflecting differences in fibril packing.

ATP Hydrolysis in the RecA-DNA Filament Promotes Structural Changes at the Protein-DNA Interface.

To address the mechanistic roles of ATP hydrolysis in RecA-promoted strand exchange reaction in homologous recombination, quantum mechanical calculations are performed on key parts of the RecA-DNA complex. We find that ATP hydrolysis may induce changes at the protein-DNA interface, resulting in the rearrangement of the hydrogen bond network connecting the ATP and the DNA binding sites.

Enhanced cellular uptake of antisecretory peptide AF-16 through proteoglycan binding.

Peptide AF-16, which includes the active site of Antisecretory Factor protein, has antisecretory and anti-inflammatory properties, making it a potent drug candidate for treatment of secretory and inflammatory diseases such as diarrhea, inflammatory bowel diseases, and intracranial hypertension. Despite remarkable physiological effects and great pharmaceutical need for drug discovery, very little is yet understood about AF-16 mechanism of action. In order to address interaction mechanisms, we investigated the binding of AF-16 to sulfated glycosaminoglycan, heparin, with focus on the effect of pH and ionic strength, and studied the influence of cell-surface proteoglycans on cellular uptake efficiency. Confocal laser scanning microscopy and flow cytometry experiments on wild type and proteoglycan-deficient Chinese hamster ovary cells reveal an endocytotic nature of AF-16 cellular uptake that is, however, less efficient for the cells lacking cell-surface proteoglycans. Isothermal titration calorimetry provides quantitative thermodynamic data and evidence for that the peptide affinity to heparin increases at lower pH and ionic strength. Experimental data, supported by theoretical modeling, of peptide-glycosaminoglycan interaction indicate that it has a large electrostatic contribution, which will be enhanced in diseases accompanied by decreased pH and ionic strength. These observations show that cell-surface proteoglycans are of general and crucial importance for the antisecretory and anti-inflammatory activities of AF-16.

Ca2+ improves organization of single-stranded DNA bases in human Rad51 filament, explaining stimulatory effect on gene recombination.

Human RAD51 protein (HsRad51) catalyses the DNA strand exchange reaction for homologous recombination. To clarify the molecular mechanism of the reaction in vitro being more effective in the presence of Ca(2+) than of Mg(2+), we have investigated the effect of these ions on the structure of HsRad51 filament complexes with single- and double-stranded DNA, the reaction intermediates. Flow linear dichroism spectroscopy shows that the two ionic conditions induce significantly different structures in the HsRad51/single-stranded DNA complex, while the HsRad51/double-stranded DNA complex does not demonstrate this ionic dependence. In the HsRad51/single-stranded DNA filament, the primary intermediate of the strand exchange reaction, ATP/Ca(2+) induces an ordered conformation of DNA, with preferentially perpendicular orientation of nucleobases relative to the filament axis, while the presence of ATP/Mg(2+), ADP/Mg(2+) or ADP/Ca(2+) does not. A high strand exchange activity is observed for the filament formed with ATP/Ca(2+), whereas the other filaments exhibit lower activity. Molecular modelling suggests that the structural variation is caused by the divalent cation interfering with the L2 loop close to the DNA-binding site. It is proposed that the larger Ca(2+) stabilizes the loop conformation and thereby the protein-DNA interaction. A tight binding of DNA, with bases perpendicularly oriented, could facilitate strand exchange.

Structural dynamics of therapeutic nucleic acids with phosphorothioate backbone modifications.

Antisense oligonucleotides (ASOs) offer ground-breaking possibilities for selective pharmacological intervention for any gene product-related disease. Therapeutic ASOs contain extensive chemical modifications that improve stability to enzymatic cleavage and modulate binding affinity relative to natural RNA/DNA. Molecular dynamics (MD) simulation can provide valuable insights into how such modifications affect ASO conformational sampling and target binding. However, force field parameters for chemically modified nucleic acids (NAs) are still underdeveloped. To bridge this gap, we developed parameters to allow simulations of ASOs with the widely applied phosphorothioate (PS) backbone modification, and validated these in extensive all-atom MD simulations of relevant PS-modified NA systems representing B-DNA, RNA, and DNA/RNA hybrid duplex structures. Compared to the corresponding natural NAs, single PS substitutions had marginal effects on the ordered DNA/RNA duplex, whereas substantial effects of phosphorothioation were observed in single-stranded RNA and B-DNA, corroborated by the experimentally derived structure data. We find that PS-modified NAs shift between high and low twist states, which could affect target recognition and protein interactions for phosphorothioated oligonucleotides. Furthermore, conformational sampling was markedly altered in the PS-modified ssRNA system compared to that of the natural oligonucleotide, indicating sequence-dependent effects on conformational preference that may in turn influence duplex formation.

Novel druggable space in human KRAS G13D discovered using structural bioinformatics and a P-loop targeting monoclonal antibody.

KRAS belongs to a family of small GTPases that act as binary switches upstream of several signalling cascades, controlling proliferation and survival of cells. Mutations in KRAS drive oncogenesis, especially in pancreatic, lung, and colorectal cancers (CRC). Although historic attempts at targeting mutant KRAS with small molecule inhibitors have proven challenging, there are recent successes with the G12C, and G12D mutations. However, clinically important RAS mutations such as G12V, G13D, Q61L, and A146T, remain elusive drug targets, and insights to their structural landscape is of critical importance to develop novel, and effective therapeutic concepts. We present a fully open, P-loop exposing conformer of KRAS G13D by X-ray crystallography at 1.4-2.4 Å resolution in Mg2+-free phosphate and malonate buffers. The G13D conformer has the switch-I region displaced in an upright position leaving the catalytic core fully exposed. To prove that this state is druggable, we developed a P-loop-targeting monoclonal antibody (mAb). The mAb displayed high-affinity binding to G13D and was shown using high resolution fluorescence microscopy to be spontaneously taken up by G13D-mutated HCT 116 cells (human CRC derived) by macropinocytosis. The mAb inhibited KRAS signalling in phosphoproteomic and genomic studies. Taken together, the data propose novel druggable space of G13D that is reachable in the cellular context. It is our hope that these findings will stimulate attempts to drug this fully open state G13D conformer using mAbs or other modalities.

Initial DNA interactions of the binuclear threading intercalator Λ,Λ-[µ-bidppz(bipy)4Ru2]4+: an NMR study with [d(CGCGAATTCGCG)]2.

Binuclear polypyridine ruthenium compounds have been shown to slowly intercalate into DNA, following a fast initial binding on the DNA surface. For these compounds, intercalation requires threading of a bulky substituent, containing one Ru(II), through the DNA base-pair stack, and the accompanying DNA duplex distortions are much more severe than with intercalation of mononuclear compounds. Structural understanding of the process of intercalation may greatly gain from a characterisation of the initial interactions between binuclear Ru(II) compounds and DNA. We report a structural NMR study on the binuclear Ru(II) intercalator Λ,Λ-B (Λ,Λ-[µ-bidppz(bipy)4Ru2](4+); bidppz = 11,11'-bis(dipyrido[3,2-a:2',3'-c]phenazinyl, bipy = 2,2'-bipyridine) mixed with the palindromic DNA [d(CGCGAATTCGCG)]2. Threading of Λ,Λ-B depends on the presence and length of AT stretches in the DNA. Therefore, the latter was selected to promote initial binding, but due to the short stretch of AT base pairs, final intercalation is prevented. Structural calculations provide a model for the interaction: Λ,Λ-B is trapped in a well-defined surface-bound state consisting of an eccentric minor-groove binding. Most of the interaction enthalpy originates from electrostatic and van der Waals contacts, whereas intermolecular hydrogen bonds may help to define a unique position of Λ,Λ-B. Molecular dynamics simulations show that this minor-groove binding mode is stable on a nanosecond scale. To the best of our knowledge, this is the first structural study by NMR spectroscopy on a binuclear Ru compound bound to DNA. In the calculated structure, one of the positively charged Ru(2+) moieties is near the central AATT region; this is favourable in view of potential intercalation as observed by optical methods for DNA with longer AT stretches. Circular dichroism (CD) spectroscopy suggests that a similar binding geometry is formed in mixtures of Λ,Λ-B with natural calf thymus DNA. The present minor-groove binding mode is proposed to represent the initial surface interactions of binuclear Ru(II) compounds prior to intercalation into AT-rich DNA.

Decoding the dual recognition mechanism of the glucocorticoid receptor for DNA and RNA: sequence versus shape.

Transcription factors (TFs) regulate eukaryotic transcription through selective DNA-binding, can also specifically interact with RNA, which may present another layer of transcriptional control. The mechanisms of the TFs-DNA recognition are often well-characterised, while the details of TFs-RNA complexation are less understood. Here we investigate the dual recognition mechanism of the glucocorticoid receptor (GR), which interacts with similar affinities with consensus DNA and diverse RNA hairpin motifs but discriminates against uniform dsRNA. Using atomic molecular dynamics simulations, we demonstrate that the GR binding to nucleic acids requires a wide and shallow groove pocket. The protein effectively moulds its binding site within DNA major groove, which enables base-specific interactions. Contrary, the GR binding has little effect on the grooves geometry of RNA systems, most notably in uniform dsRNA. Instead, a hairpin motif in RNA yields a wide and shallow major groove pocket, allowing the protein to anchor itself through nonspecific electrostatic contacts with RNA backbone. Addition of a bulge increases RNA hairpin flexibility, which leads to a greater number of GR-RNA contacts and, thus, higher affinity. Thus, the combination of structural motifs defines the GR-RNA selective binding: a recognition mechanism, which may be shared by other zinc finger TFs.

Structure of human Rad51 protein filament from molecular modeling and site-specific linear dichroism spectroscopy.

To get mechanistic insight into the DNA strand-exchange reaction of homologous recombination, we solved a filament structure of a human Rad51 protein, combining molecular modeling with experimental data. We build our structure on reported structures for central and N-terminal parts of pure (uncomplexed) Rad51 protein by aid of linear dichroism spectroscopy, providing angular orientations of substituted tyrosine residues of Rad51-dsDNA filaments in solution. The structure, validated by comparison with an electron microscopy density map and results from mutation analysis, is proposed to represent an active solution structure of the nucleo-protein complex. An inhomogeneously stretched double-stranded DNA fitted into the filament emphasizes the strategic positioning of 2 putative DNA-binding loops in a way that allows us speculate about their possibly distinct roles in nucleo-protein filament assembly and DNA strand-exchange reaction. The model suggests that the extension of a single-stranded DNA molecule upon binding of Rad51 is ensured by intercalation of Tyr-232 of the L1 loop, which might act as a docking tool, aligning protein monomers along the DNA strand upon filament assembly. Arg-235, also sitting on L1, is in the right position to make electrostatic contact with the phosphate backbone of the other DNA strand. The L2 loop position and its more ordered compact conformation makes us propose that this loop has another role, as a binding site for an incoming double-stranded DNA. Our filament structure and spectroscopic approach open the possibility of analyzing details along the multistep path of the strand-exchange reaction.

Homologous basic helix-loop-helix transcription factors induce distinct deformations of torsionally-stressed DNA: a potential transcription regulation mechanism.

Changing torsional restraints on DNA is essential for the regulation of transcription. Torsional stress, introduced by RNA polymerase, can propagate along chromatin facilitating topological transitions and modulating the specific binding of transcription factors (TFs) to DNA. Despite the importance, the mechanistic details on how torsional stress impacts the TFs-DNA complexation remain scarce. Herein, we address the impact of torsional stress on DNA complexation with homologous human basic helix-loop-helix (BHLH) hetero- and homodimers: MycMax, MadMax and MaxMax. The three TF dimers exhibit specificity towards the same DNA consensus sequence, the E-box response element, while regulating different transcriptional pathways. Using microseconds-long atomistic molecular dynamics simulations together with the torsional restraint that controls DNA total helical twist, we gradually over- and underwind naked and complexed DNA to a maximum of ± 5°/bp step. We observe that the binding of the BHLH dimers results in a similar increase in DNA torsional rigidity. However, under torsional stress the BHLH dimers induce distinct DNA deformations, characterised by changes in DNA grooves geometry and a significant asymmetric DNA bending. Supported by bioinformatics analyses, our data suggest that torsional stress may contribute to the execution of differential transcriptional programs of the homologous TFs by modulating their collaborative interactions.

Abnormal methylation in the NDUFA13 gene promoter of breast cancer cells breaks the cooperative DNA recognition by transcription factors.

Selective DNA binding by transcription factors (TFs) is crucial for the correct regulation of DNA transcription. In healthy cells, promoters of active genes are hypomethylated. A single CpG methylation within a TF response element (RE) may change the binding preferences of the protein, thus causing the dysregulation of transcription programs. Here, we investigate a molecular mechanism driving the downregulation of the NDUFA13 gene, due to hypermethylation, which is associated with multiple cancers. Using bioinformatic analyses of breast cancer cell line MCF7, we identify a hypermethylated region containing the binding sites of two TFs dimers, CEBPB and E2F1-DP1, located 130 b.p. from the gene transcription start site. All-atom extended MD simulations of wild type and methylated DNA alone and in complex with either one or both TFs dimers provide mechanistic insights into the cooperative asymmetric binding order of the two dimers; the CEBPB binding should occur first to facilitate the E2F1-DP1-DNA association. The CpG methylation within the E2F1-DP1 RE and the linker decrease the cooperativity effects and renders the E2F1-DP1 binding site less recognizable by the TF dimer. Taken together, the identified CpG methylation site may contribute to the downregulation of the NDUFA13 gene.

Rational antibody design for undruggable targets using kinetically controlled biomolecular probes.

Several important drug targets, e.g., ion channels and G protein-coupled receptors, are extremely difficult to approach with current antibody technologies. To address these targets classes, we explored kinetically controlled proteases as structural dynamics-sensitive druggability probes in native-state and disease-relevant proteins. By using low-Reynolds number flows, such that a single or a few protease incisions are made, we could identify antibody binding sites (epitopes) that were translated into short-sequence antigens for antibody production. We obtained molecular-level information of the epitope-paratope region and could produce high-affinity antibodies with programmed pharmacological function against difficult-to-drug targets. We demonstrate the first stimulus-selective monoclonal antibodies targeting the transient receptor potential vanilloid 1 (TRPV1) channel, a clinically validated pain target widely considered undruggable with antibodies, and apoptosis-inducing antibodies selectively mediating cytotoxicity in KRAS-mutated cells. It is our hope that this platform will widen the scope of antibody therapeutics for the benefit of patients.

Specifically bound BZIP transcription factors modulate DNA supercoiling transitions.

Torsional stress on DNA, introduced by molecular motors, constitutes an important regulatory mechanism of transcriptional control. Torsional stress can modulate specific binding of transcription factors to DNA and introduce local conformational changes that facilitate the opening of promoters and nucleosome remodelling. Using all-atom microsecond scale molecular dynamics simulations together with a torsional restraint that controls the total twist of a DNA fragment, we address the impact of torsional stress on DNA complexation with a human BZIP transcription factor, MafB. We gradually over- and underwind DNA alone and in complex with MafB by 0.5° per dinucleotide step, starting from the relaxed state to a maximum of 5° per dinucleotide step, monitoring the evolution of the protein-DNA contacts at different degrees of torsional strain. Our computations show that MafB changes the DNA sequence-specific response to torsional stress. The dinucleotide steps that are susceptible to absorbing most of the torsional stress become more torsionally rigid, as they are involved in protein-DNA contacts. Also, the protein undergoes substantial conformational changes to follow the stress-induced DNA deformation, but mostly maintains the specific contacts with DNA. This results in a significant asymmetric increase of free energy of DNA twisting transitions, relative to free DNA, where overtwisting is more energetically unfavourable. Our data suggest that specifically bound BZIP factors could act as torsional stress insulators, modulating the propagation of torsional stress along the chromatin fibre, which might promote cooperative binding of collaborative DNA-binding factors.