Chromatin methylation activity of Dnmt3a and Dnmt3a/3L is guided by interaction of the ADD domain with the histone H3 tail

Using peptide arrays and binding to native histone proteins, we show that the ADD domain of Dnmt3a specifically interacts with the H3 histone 1–19 tail. Binding is disrupted by di- and trimethylation of K4, phosphorylation of T3, S10 or T11 and acetylation of K4. We did not observe binding to the H4 1–19 tail. The ADD domain of Dnmt3b shows the same binding specificity, suggesting that the distinct biological functions of both enzymes are not related to their ADD domains. To establish a functional role of the ADD domain binding to unmodified H3 tails, we analyzed the DNA methylation of in vitro reconstituted chromatin with Dnmt3a2, the Dnmt3a2/Dnmt3L complex, and the catalytic domain of Dnmt3a. All Dnmt3a complexes preferentially methylated linker DNA regions. Chromatin substrates with unmodified H3 tail or with H3K9me3 modification were methylated more efficiently by full-length Dnmt3a and full-length Dnmt3a/3L complexes than chromatin trimethylated at H3K4. In contrast, the catalytic domain of Dnmt3a was not affected by the H3K4me3 modification. These results demonstrate that the binding of the ADD domain to H3 tails unmethylated at K4 leads to the preferential methylation of DNA bound to chromatin with this modification state. Our in vitro results recapitulate DNA methylation patterns observed in genome-wide DNA methylation studies.     

Characterization of suramin binding sites on the human group IIA secreted phospholipase A2 by site-directed mutagenesis and molecular dynamics simulation

Suramin is a polysulphonated naphthylurea with inhibitory activity against the human secreted group IIA phospholipase A(2) (hsPLA2GIIA), and we have investigated suramin binding to recombinant hsPLA2GIIA using site-directed mutagenesis and molecular dynamics (MD) simulations. The changes in suramin binding affinity of 13 cationic residue mutants of the hsPLA2GIIA was strongly correlated with alterations in the inhibition of membrane damaging activity of the protein. Suramin binding to hsPLA2GIIA was also studied by MD simulations, which demonstrated that altered intermolecular potential energy of the suramin/mutant complexes was a reliable indicator of affinity change. Although residues in the C-terminal region play a major role in the stabilization of the hsPLA2GIIA/suramin complex, attractive and repulsive hydrophobic and electrostatic interactions with residues throughout the protein together with the adoption of a bent suramin conformation, all contribute to the stability of the complex. Analysis of the hsPLA2GIIA/suramin interactions allows the prediction of the properties of suramin analogues with improved binding and higher affinities which may be candidates for novel phospholipase A(2) inhibitors.     

Sliding of Proteins Non-specifically Bound to DNA: Brownian Dynamics Studies with Coarse-Grained Protein and DNA Models

DNA binding proteins efficiently search for their cognitive sites on long genomic DNA by combining 3D diffusion and 1D diffusion (sliding) along the DNA. Recent experimental results and theoretical analyses revealed that the proteins show a rotation-coupled sliding along DNA helical pitch. Here, we performed Brownian dynamics simulations using newly developed coarse-grained protein and DNA models for evaluating how hydrodynamic interactions between the protein and DNA molecules, binding affinity of the protein to DNA, and DNA fluctuations affect the one dimensional diffusion of the protein on the DNA. Our results indicate that intermolecular hydrodynamic interactions reduce 1D diffusivity by 30%. On the other hand, structural fluctuations of DNA give rise to steric collisions between the CG-proteins and DNA, resulting in faster 1D sliding of the protein. Proteins with low binding affinities consistent with experimental estimates of non-specific DNA binding show hopping along the CG-DNA. This hopping significantly increases sliding speed. These simulation studies provide additional insights into the mechanism of how DNA binding proteins find their target sites on the genome.     

Antigen-specific suppressor T lymphocytes in man make use of the same set of surface molecules as do cytolytic T lymphocytes: roles of Leu-2/T8, Leu-4/T3, Leu-5/T11, LFA-1 molecules

When cultured with autologous antigen-primed Leu-3+ lymphoblasts, Leu-2+ cells differentiate into suppressor T cells (Ts) that specifically inhibit the responses of fresh autologous Leu-3+ cells to the priming antigen. We have shown previously that the Leu-4/T3 (CD-3) molecular complex and HLA-A,B molecules on the surface of Leu-3+ inducer blasts are recognized by Leu-2+ Ts during their differentiation. This study examines the role of various cell surface molecules expressed by Leu-2+ Ts during the inductive and effector phases of suppression. Leu-2+ cells were treated in the absence of complement with a variety of monoclonal antibodies recognizing distinct human lymphoid antigens either before or after their activation with alloantigen-primed Leu-3+ blasts. Antibodies to Leu-2/T8 (CD-8) and lymphocyte function-associated antigen-1 (LFA-1) (CDw-18) molecules inhibited not only the generation but also the effector function of Leu-2+ Ts. Although antibodies to Leu-4/T3 (CD-3) and Leu-5/T11 (CD-2) molecules caused profound inhibition of the activation of Ts, these antibodies failed to inhibit the effector function of Ts. On the contrary, anti-Leu-4 antibody consistently augmented the suppressor effect of Ts. Antibodies directed against Leu-1/T1 (CD-5), Leu-3/T4 (CD-4), LFA-3, and class I (HLA-A,B,C) and class II (HLA-DR,DQ) major histocompatibility complex molecules had no effect on either the generation or the effector function of Ts. These results suggest the involvement of Leu-2/T8 (CD-8), Leu-4/T3 (CD-3), Leu-5/T11 (CD-2), and LFA-1 (CDw-18) molecules on the surfaces of Leu-2+ cells in the activation and effector functions of Ts.     

Ubiquitin-mediated degradation of histone H3 does not require the substrate-binding ubiquitin protein ligase, E3, or attachment of polyubiquitin chains

Radioiodinated histone H3 was incubated with ubiquitin, the ubiquitin-activating enzyme E1, and one of three ubiquitin carrier proteins, reticulocyte E2(20K) or E2(32K) or the yeast RAD6 product. Although the resulting ubiquitin-histone conjugates were synthesized in the absence of the substrate-binding protein E3, they were nevertheless degraded by purified rabbit reticulocyte 26 S protease. In contrast, unmodified histone H3 remained intact upon challenge with the 26 S ubiquitin/ATP-dependent enzyme. Conjugates produced by the RAD6 protein were better proteolytic substrates than those formed by reticulocyte E2 unless ubiquitin molecules with altered lysines were used for conjugate synthesis. Substitution of methylated ubiquitin or ubiquitin molecules in which lysine 48 was converted to arginine by site-directed mutation produced histone conjugates that were degraded at slow but measurable rates. Since methylated ubiquitin molecules are incapable of forming branched polyubiquitin chains, these results demonstrate that neither ubiquitin "trees" nor the substrate binding factor E3 is absolutely required for ubiquitin-dependent degradation of histone H3 in vitro.     

Mechanistic insights into enhancement or inhibition of phase separation by different polyubiquitin chains

Ubiquitin‐binding shuttle UBQLN2 mediates crosstalk between proteasomal degradation and autophagy, likely via interactions with K48‐ and K63‐linked polyubiquitin chains, respectively. UBQLN2 comprises self‐associating regions that drive its homotypic liquid–liquid phase separation (LLPS). Specific interactions between one of these regions and ubiquitin inhibit UBQLN2 LLPS. Here, we show that, unlike ubiquitin, the effects of multivalent polyubiquitin chains on UBQLN2 LLPS are highly dependent on chain types. Specifically, K11‐Ub4 and K48‐Ub4 chains generally inhibit UBQLN2 LLPS, whereas K63‐Ub4, M1‐Ub4 chains, and a designed tetrameric ubiquitin construct significantly enhance LLPS. We demonstrate that these opposing effects stem from differences in chain conformations but not in affinities between chains and UBQLN2. Chains with extended conformations and increased accessibility to the ubiquitin‐binding surface promote UBQLN2 LLPS by enabling a switch between homotypic to partially heterotypic LLPS that is driven by both UBQLN2 self‐interactions and interactions between multiple UBQLN2 units with each polyubiquitin chain. Our study provides mechanistic insights into how the structural and conformational properties of polyubiquitin chains contribute to heterotypic LLPS with ubiquitin‐binding shuttles and adaptors.     

Water and ion binding around RNA and DNA (C,G) oligomers

The dynamics, hydration, and ion-binding features of two duplexes, the A(r(CG)(12)) and the B(d(CG)(12)), in a neutralizing aqueous environment with 0.25 M added KCl have been investigated by molecular dynamics (MD) simulations. The regular repeats of the same C=G base-pair motif have been exploited as a statistical alternative to long MD simulations in order to extend the sampling of the conformational space. The trajectories demonstrate the larger flexibility of DNA compared to RNA helices. This flexibility results in less well defined hydration patterns around the DNA than around the RNA backbone atoms. Yet, 22 hydration sites are clearly characterized around both nucleic acid structures. With additional results from MD simulations, the following hydration scale for C=G pairs can be deduced: A-DNA相似文献   

An In Silico Analysis of the Binding Modes and Binding Affinities of Small Molecule Modulators of PDZ-Peptide Interactions

Inhibitors of PDZ-peptide interactions have important implications in a variety of biological processes including treatment of cancer and Parkinson’s disease. Even though experimental studies have reported characterization of peptidomimetic inhibitors of PDZ-peptide interactions, the binding modes for most of them have not been characterized by structural studies. In this study we have attempted to understand the structural basis of the small molecule-PDZ interactions by in silico analysis of the binding modes and binding affinities of a set of 38 small molecules with known K_i or K_d values for PDZ2 and PDZ3 domains of PSD-95 protein. These two PDZ domains show differential selectivity for these compounds despite having a high degree of sequence similarity and almost identical peptide binding pockets. Optimum binding modes for these ligands for PDZ2 and PDZ3 domains were identified by using a novel combination of semi-flexible docking and explicit solvent molecular dynamics (MD) simulations. Analysis of the binding modes revealed most of the peptidomimectic ligands which had high K_i or K_d moved away from the peptide binding pocket, while ligands with high binding affinities remained in the peptide binding pocket. The differential specificities of the PDZ2 and PDZ3 domains primarily arise from differences in the conformation of the loop connecting βB and βC strands, because this loop interacts with the N-terminal chemical moieties of the ligands. We have also computed the MM/PBSA binding free energy values for these 38 compounds with both the PDZ domains from multiple 5 ns MD trajectories on each complex i.e. a total of 228 MD trajectories of 5 ns length each. Interestingly, computational binding free energies show good agreement with experimental binding free energies with a correlation coefficient of approximately 0.6. Thus our study demonstrates that combined use of docking and MD simulations can help in identification of potent inhibitors of PDZ-peptide complexes.     

Two-site binding of beta-cyclodextrin to the Alzheimer Abeta(1-40) peptide measured with combined PFG-NMR diffusion and induced chemical shifts

The interactions of Alzheimer's amyloid beta-peptide with cyclodextrins were studied by (1)H NMR: the translational diffusion coefficient of the peptide and chemical shift changes were studied by the presence of variable concentrations of cyclodextrins. For the full-length peptide, Abeta(1-40), the combined results of translational diffusion and chemical shift changes are consistent with a model where aromatic side chains interact with beta-cyclodextrin with dissociation constants in the millimolar range. The diffusion data were consistent with two beta-cyclodextrin molecules bound per peptide. The binding occurs at two sites, at F(19) and/or F(20) and at Y(10), with dissociation constants K(d)(F) = 4.7 mM and K(d)(Y) = 6.6 mM, respectively, in 10 mM sodium phosphate, pH 7.4 and 298 K. Shorter Alzheimer peptide fragments were studied to measure specific affinities for different binding sites. The N-terminal fragment Abeta(1-9) with a putative binding site at F(4) does not show measurable affinity for beta-cyclodextrin. The fragment Abeta(12-28) has similar apparent affinity (K(d) = 3.8 mM) to beta-cyclodextrin as the full-length peptide Abeta(1-40). Here, the diffusion data suggests a one-to-one stoichiometry, and the binding site is F(19) and/or F(20). Both diffusion results and chemical shift changes give the same affinity. A variant Abeta(12-28)G(19)G(20) without phenylalanines does not bind to beta-cyclodextrin. Other potential ligands, alpha-cyclodextrin, gamma-cyclodextrin, nicotine, and nornicotine do not bind to the Abeta(12-28) fragment. This study shows that combined (1)H NMR diffusion and chemical shift changes may be used to quantitatively determine affinities and stoichiometries of weak interactions, using unlabeled ligands and hosts of comparable sizes.     

Comparative molecular dynamics simulations of HIV-1 integrase and the T66I/M154I mutant: binding modes and drug resistance to a diketo acid inhibitor

HIV-1 IN is an essential enzyme for viral replication and an interesting target for the design of new pharmaceuticals for use in multidrug therapy of AIDS. L-731,988 is one of the most active molecules of the class of beta-diketo acids. Individual and combined mutations of HIV-1 IN at residues T66, S153, and M154 confer important degrees of resistance to one or more inhibitors belonging to this class. In an effort to understand the molecular mechanism of the resistance of T66I/M154I IN to the inhibitor L-731,988 and its specific binding modes, we have carried out docking studies, explicit solvent MD simulations, and binding free energy calculations. The inhibitor was docked against different protein conformations chosen from prior MD trajectories, resulting in 2 major orientations within the active site. MD simulations have been carried out for the T66I/M154I DM IN, DM IN in complex with L-731,988 in 2 different orientations, and 1QS4 IN in complex with L-731,988. The results of these simulations show a similar dynamical behavior between T66I/M154I IN alone and in complex with L-731,988, while significant differences are observed in the mobility of the IN catalytic loop (residues 138-149). Water molecules bridging the inhibitor to residues from the active site have been identified, and residue Gln62 has been found to play an important role in the interactions between the inhibitor and the protein. This work provides information about the binding modes of L-731,988, as well as insight into the mechanism of inhibitor-resistance in HIV-1 integrase.     

A comparative study of trypsin specificity based on QM/MM molecular dynamics simulation and QM/MM GBSA calculation

Hydrogen bonding and polar interactions play a key role in identification of protein-inhibitor binding specificity. Quantum mechanics/molecular mechanics molecular dynamics (QM/MM MD) simulations combined with DFT and semi-empirical Hamiltonian (AM1d, RM1, PM3, and PM6) methods were performed to study the hydrogen bonding and polar interactions of two inhibitors BEN and BEN1 with trypsin. The results show that the accuracy of treating the hydrogen bonding and polar interactions using QM/MM MD simulation of PM6 can reach the one obtained by the DFT QM/MM MD simulation. Quantum mechanics/molecular mechanics generalized Born surface area (QM/MM-GBSA) method was applied to calculate binding affinities of inhibitors to trypsin and the results suggest that the accuracy of binding affinity prediction can be significantly affected by the accurate treatment of the hydrogen bonding and polar interactions. In addition, the calculated results also reveal the binding specificity of trypsin: (1) the amidinium groups of two inhibitors generate favorable salt bridge interaction with Asp189 and form hydrogen bonding interactions with Ser190 and Gly214, (2) the phenyl of inhibitors can produce favorable van der Waals interactions with the residues His58, Cys191, Gln192, Trp211, Gly212, and Cys215. This systematic and comparative study can provide guidance for the choice of QM/MM MD methods and the designs of new potent inhibitors targeting trypsin.     

U1A protein-stem loop 2 RNA recognition: prediction of structural differences from protein mutations

Molecular dynamics (MD) simulations were carried out to compare the free and bound structures of wild type U1A protein with several Phe56 mutant U1A proteins that bind the target stem loop 2 (SL2) RNA with a range of affinities. The simulations indicate the free U1A protein is more flexible than the U1A-RNA complex for both wild type and Phe56 mutant systems. A complete analysis of the hydrogen-bonding (HB) and non-bonded (VDW) interactions over the course of the MD simulations suggested that changes in the interactions in the free U1A protein caused by the Phe56Ala and Phe56Leu mutations may stabilize the helical character in loop 3, and contribute to the weak binding of these proteins to SL2 RNA. Compared with wild type, changes in HB and VDW interactions in Phe56 mutants of the free U1A protein are global, and include differences in β-sheet, loop 1 and loop 3 interactions. In the U1A-RNA complex, the Phe56Ala mutation leads to a series of differences in interactions that resonate through the complex, while the Phe56Leu and Phe56Trp mutations cause local differences around the site of mutation. The long-range networks of interactions identified in the simulations suggest that direct interactions and dynamic processes in both the free and bound forms contribute to complex stability.     

Cationic 5,10,15,20-tetrakis(N-methylpyridinium-4-yl)porphyrin fully intercalates at 5'-CG-3' steps of duplex DNA in solution

The interaction of 5,10,15, 20-tetrakis(N-methylpyridinium-4-yl)porphyrin (T4MPyP(4+)) with the oligonucleotide DNA duplex [d(GCACGTGC)](2) was studied by two-dimensional (1)H NMR spectroscopy, optical absorbance, circular dichroism, and molecular dynamics simulation employing particle mesh Ewald methods. T4MPyP(4+) is one of the largest aromatic molecules for which intercalative binding to DNA has been proposed, although this has been called into question by recent X-ray crystallographic work [Lipscomb et al. (1996) Biochemistry 35, 2818-2823]. T4MPyP(4+) binding to [d(GCACGTGC)](2) produced a single set of (mostly) upfield-shifted DNA resonances in slow exchange with the resonances of the free DNA. Intra- and intermolecular NOEs observed in the complex showed that the porphyrin intercalates at the central 5'-CG-3' step of the DNA duplex without disrupting the flanking base pairs. Absorption and circular dichroism spectra of the complex also support intercalative binding. Molecular dynamics simulations (using explicit solvent and PME methods), carried out for fully and partially intercalated complexes, yielded stable trajectories and plausible structures, but only the symmetrical, fully intercalated model agreed with NOESY data. Stable hydrogen bonding was observed during 600 ps of MD simulation for the base pairs flanking the binding site.     

Effects of phospholipid unsaturation on the membrane/water interface: a molecular simulation study.

Molecular dynamics (MD) simulations of fully hydrated bilayers in the liquid-crystalline state made of 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) or 1-palmitoyl-2-elaidoyl-phosphatidylcholine (PEPC) were carried out to investigate the effect of the incorporation of a double bond in the phosphatidylcholine (PC) beta-chain (cis or trans) on the membrane/water interface. The bilayers reached thermal equilibrium after 3 and 1 ns of MD simulations, respectively, and productive runs were carried out for 3 ns for each bilayer. As reference systems, the 1,2-dimyristoyl-phosphatidylcholine (DMPC) bilayer (M. Pasenkiewicz-Gierula, Y. Takaoka, H. Miyagawa, K. Kitamura, and A. Kusumi, 1999, Biophys. J. 76:1228-1240) and DMPC-cholesterol (Chol) bilayer containing 22 mol % Chol (M. Pasenkiewicz-Gierula, T. Róg, K. Kitamura, A. and Kusumi, 2000, Biophys. J. 78:1376-1389) were used. The study shows that at the interface of POPC, PEPC, and DMPC-Chol bilayers, average numbers of PC-water and PC-PC interactions are similar and, respectively, greater and smaller than in the DMPC bilayer. The average area/PC in mono-unsaturated bilayers is approximately 4 A(2) larger than in the DMPC bilayer; nevertheless, a strong correlation was found between a single molecular area (SMA) of a PC and the number of interactions this PC makes; i.e., PCs (either saturated or unsaturated) with the same SMA form similar numbers of intermolecular links. The numbers and corresponding SMAs are distributed about averages pertinent to each bilayer. No significant difference between cis and trans bonds was found.     

Mouse Dnmt3a preferentially methylates linker DNA and is inhibited by histone H1

In mammals, DNA methylation is crucial for embryonic development and germ cell differentiation. The DNA methylation patterns are created by de novo-type DNA methyltransferases (Dnmts) 3a and 3b. Dnmt3a is crucial for global methylation, including that of imprinted genes in germ cells. In eukaryotic nuclei, genomic DNA is packaged into multinucleosomes with linker histone H1, which binds to core nucleosomes, simultaneously making contacts in the linker DNA that separates adjacent nucleosomes. In the present study, we prepared oligonucleosomes from HeLa nuclei with or without linker histone H1 and used them as a substrate for Dnmt3a. Removal of histone H1 enhanced the DNA methylation activity. Furthermore, Dnmt3a preferentially methylated the linker between the two nucleosome core regions of reconstituted dinucleosomes, and the binding of histone H1 inhibited the DNA methylation activity of Dnmt3a towards the linker DNA. Since an identical amount of histone H1 did not inhibit the activity towards naked DNA, the inhibitory effect of histone H1 was not on the Dnmt3a catalytic activity but on its preferential location in the linker DNA of the dinucleosomes. The central globular domain and C-terminal tail of the histone H1 molecule were indispensable for inhibition of the DNA methylation activity of Dnmt3a. We propose that the binding and release of histone H1 from the linker portion of chromatin may regulate the local DNA methylation of the genome by Dnmt3a, which is expressed ubiquitously in somatic cells in vivo.     

Molecular Analysis of the Interaction between Staphylococcal Virulence Factor Sbi-IV and Complement C3d

Staphylococcus aureus expresses numerous virulence factors that aid in immune evasion. The four-domain staphylococcal immunoglobulin binding (Sbi) protein interacts with complement component 3 (C3) and its thioester domain (C3d)-containing fragments. Recent structural data suggested two possible modes of binding of Sbi domain IV (Sbi-IV) to C3d, but the physiological binding mode remains unclear. We used a computational approach to provide insight into the C3d-Sbi-IV interaction. Molecular dynamics (MD) simulations showed that the first binding mode (PDB code 2WY8) is more robust than the second (PDB code 2WY7), with more persistent polar and nonpolar interactions, as well as conserved interfacial solvent accessible surface area. Brownian dynamics and steered MD simulations revealed that the first binding mode has faster association kinetics and maintains more stable intermolecular interactions compared to the second binding mode. In light of available experimental and structural data, our data confirm that the first binding mode represents Sbi-IV interaction with C3d (and C3) in a physiological context. Although the second binding mode is inherently less stable, we suggest a possible physiological role. Both binding sites may serve as a template for structure-based design of novel complement therapeutics.