Molecular modeling study of the norfloxacin-DNA complex

Molecular modeling and molecular dynamics were performed to investigate the interaction of norfloxacin with the DNA oligonucleotide 5'-d(ATACGTAT)(2). Eight quinolone-DNA binding structures were built by molecular modeling on the basis of experimental results. A 100ps molecular dynamics calculation was carried out on two groove binding models and six partially intercalating models. The resulting average structures were compared with each other and to free DNA structure as a reference. The favorable binding mode of norfloxacin to a DNA substrate was pursued by structural assess including steric hindrance, presence of hydrogen-bonding, non-bonding energies of the complex and presence of abnormal structural distortion. Although two of the intercalative models showed the highest binding energy and the lowest non-bonding interaction energy, they presented structural features which contrast with experimental results. On the other hand, one groove binding model demonstrated the most acceptable structure when the experimental observation was accounted. In this model, hydrogen bonding of the carbonyl and carboxyl group of the norfloxacin rings with the DNA bases was present, and norfloxacin binds to the amine group of the guanine base which protrudes toward the minor groove of B-DNA.     

Methylene Blue Binding to DNA with Alternating AT Base Sequence: Minor Groove Binding is Favored over Intercalation

Abstract

The results presented in this paper on methylene blue (MB) binding to DNA with AT alternating base sequence complement the data obtained in two former modeling studies of MB binding to GC alternating DNA. In the light of the large amount of experimental data for both systems, this theoretical study is focused on a detailed energetic analysis and comparison in order to understand their different behavior. Since experimental high-resolution structures of the complexes are not available, the analysis is based on energy minimized structural models of the complexes in different binding modes. For both sequences, four different intercalation structures and two models for MB binding in the minor and major groove have been proposed. Solvent electrostatic effects were included in the energetic analysis by using electrostatic continuum theory, and the dependence of MB binding on salt concentration was investigated by solving the non-linear Poisson-Boltzmann equation. We find that the relative stability of the different complexes is similar for the two sequences, in agreement with the interpretation of spectroscopic data. Subtle differences, however, are seen in energy decompositions and can be attributed to the change from symmetric 5′-YpR-3′ intercalation to minor groove binding with increasing salt concentration, which is experimentally observed for the AT sequence at lower salt concentration than for the GC sequence. According to our results, this difference is due to the significantly lower non-electrostatic energy for the minor groove complex with AT alternating DNA, whereas the slightly lower binding energy to this sequence is caused by a higher deformation energy of DNA. The energetic data are in agreement with the conclusions derived from different spectroscopic studies and can also be structurally interpreted on the basis of the modeled complexes. The simple static modeling technique and the neglect of entropy terms and of non-electrostatic solute-solvent interactions, which are assumed to be nearly constant for the compared complexes of MB with DNA, seem to be justified by the results.     

A Study of the Interaction of Cr(III) Complexes and Their Selective Binding with B-DNA: A Molecular Modeling Approach

Abstract

Molecular modeling and energy minimisation calculations have been used to investigate the interaction of chromium(III) complexes in different ligand environments with various sequences of B-DNA. The complexes are [Cr(salen)(H₂O)₂]⁺; salen denotes 1, 2 bis-salicylideneaminoethane, [Cr(salprn)(H₂O)₂]⁺; salprn denotes 1, 3 bis- salicylideneamino-propane, [Cr(phen)₃]³⁺; phen denotes 1, 10 phenanthroline and [Cr(en)₃]³⁺; en denotes eth- ylenediamine. All the chromium(III) complexes are interacted with the minor groove and major groove of d(AT)₁₂, d(CGCGAATTCGCG)₂ and d(GC)₁₂ sequences of DNA. The binding energy and hydrogen bond parameters of DNA-Cr complex adduct in both the groove have been determined using molecular mechanics approach. The binding energy and formation of hydrogen bonds between chromium(III) complex and DNA has shown that all complexes of chromium(III) prefer minor groove interaction as the favourable binding mode.     

The effect of dimerization on the interaction of ibuprofen drug with calf thymus DNA: Molecularmodeling and spectroscopic investigation

The interaction between the dimer structure of ibuprofen drug (D-IB) and calf thymus DNA under simulative physiological conditions was investigated with the use of Hoechst 33258 and methylene blue dye as spectral probes by the methods of UV-visible absorption, fluorescence spectroscopy, circular dichroism spectroscopy and molecular modeling study.Using the Job's plot, a single class of binding sites for theD-IB on DNA was put in evidence. The Stern–Volmer analysis of fluorescence quenching data shows the presence of both the static and dynamic quenching mechanisms. The binding constants, K_b were calculated at different temperatures, and the thermodynamic parameters ?G^°, ?H^° and ?S^° were given. The experimental results showed that D-IB molecules could bind with DNA via groove binding mode as evidenced by: I. DNA binding constant from spectrophotometric studies of the interaction of D-IB with DNA is comparable to groove binding drugs. II. Competitive fluorimetric studies with Hoechst 33258 have shown that D-IB exhibits the ability of this complex to displace with DNA-bounded Hoechst, indicating that it binds to DNA in strong competition with Hoechst for the groove binding. III. There is no significantly change in the absorption of the MB-DNA system upon adding the D-IB, indicates that MB molecules are not released from the DNA helix after addition of the D-IB and are indicative of a non-intercalative mode of binding. IV. Small changes in DNA viscosity in the presence of D-IB, indicating weak link to DNA, which is consistent with DNA groove binding. As well as, induced CD spectral changes, and the docking results revealed that groove mechanism is followed by D-IB to bind with DNA.     

DNA-binding study of anticancer drug cytarabine by spectroscopic and molecular docking techniques

The interaction of anticancer drug cytarabine with calf thymus DNA (CT-DNA) was investigated in vitro under simulated physiological conditions by multispectroscopic techniques and molecular modeling study. The fluorescence spectroscopy and UV absorption spectroscopy indicated drug interacted with CT-DNA in a groove-binding mode, while the binding constant of UV-vis and the number of binding sites were 4.0 ± 0.2 × 10⁴ L mol^?1 and 1.39, respectively. The fluorimetric studies showed that the reaction between the drugs with CT-DNA is exothermic. Circular dichroism spectroscopy was employed to measure the conformational change of DNA in the presence of cytarabine. Furthermore, the drug induces detectable changes in its viscosity for DNA interaction. The molecular modeling results illustrated that cytarabine strongly binds to groove of DNA by relative binding energy of docked structure ?20.61 KJ mol^?1. This combination of multiple spectroscopic techniques and molecular modeling methods can be widely used in the investigation on the interaction of small molecular pollutants and drugs with biomacromolecules for clarifying the molecular mechanism of toxicity or side effect in vivo.     

Binding free energy of selected anticancer compounds to DNA--theoretical calculations

Many studies have elucidated structures and thermodynamics of complexes formed by different ligands with DNA. However, in most cases structural and free energy binding studies were not correlated with each other because of the problem of identifying which experimental free energy of binding corresponds to which experimental DNA-ligand structure. In the present work, Poisson-Boltzmann and solvent-accessible surface area methods were used to predict unknown modes of interaction between DNA and three different ligands: mitoxantrone and two pyrimidoacridine derivatives. In parallel, experimental measurements of binding free energy for the studied complexes were performed to compare experimental and calculated values. Our studies showed that the calculated values of free energy are only close to experimental data for some models of interaction between ligands and DNA. Based on this correlation, the most likely models of DNA-ligand complexes were postulated: (i) mitoxantrone and one derivative of pyrimidoacridine, both with two charged side chains, intercalate from the minor groove of DNA and bind with both chains in this groove; (ii) pyrimidoacridine, with only one side chain, very likely does not intercalate into DNA at all. Additionally, the non-electrostatic and electrostatic parts of the calculated binding free energy for the DNA-ligands studied are discussed.     

DNA binding studies of antibiotic drug cephalexin using spectroscopic and molecular docking techniques

The purpose of this study was to explore an accurate characterization of the binding interaction of antibiotic drug cephalexin with calf thymus DNA (CT-DNA) as a relevant biological target by using UV absorption, fluorescence spectroscopy and circular dichroism (CD) in vitro under simulated physiological conditions (pH = 7.4) and also through a molecular modeling study. The results showed that the drug interacts with the DNA helix via a minor groove binding mode. The thermodynamic parameters were calculated and showed that the reaction between the drug and CT-DNA was exothermic. In addition, the drug enforced traceable changes in the viscosity of DNA. The molecular modeling results indicated that cephalexin forcefully binds to the minor groove of DNA with a relative binding energy of ?21.02?kJ mol^?1. The obtained theoretical results were in good agreement with those obtained from experimental studies.     

Binding properties of palmatine to DNA: spectroscopic and molecular modeling investigations

Palmatine, an isoquinoline alkaloid, is an important medicinal herbal extract with diverse pharmacological and biological properties. In this work, spectroscopic and molecular modeling approaches were employed to reveal the interaction between palmatine and DNA isolated from herring sperm. The absorption spectra and iodide quenching results indicated that groove binding was the main binding mode of palmatine to DNA. Fluorescence studies indicated that the binding constant (K) of palmatine and DNA was ~ 10⁴ L·mol^?1. The associated thermodynamic parameters, ΔG, ΔH, and ΔS, indicated that hydrogen bonds and van der Waals forces played major roles in the interaction. The effects of chemical denaturant, thermal denaturation and pH on the interaction were investigated and provided further support for the groove binding mode. In addition to experimental approaches, molecular modeling was conducted to verify binding pattern of palmatine–DNA. Copyright © 2015 John Wiley & Sons, Ltd.     

Molecular docking and dynamics simulations on the interaction of cationic porphyrin–anthraquinone hybrids with DNA G-quadruplexes

A series of cationic porphyrin–anthraquinone hybrids bearing either pyridine, imidazole, or pyrazole rings at the meso-positions have been investigated for their interaction with DNA G-quadruplexes by employing molecular docking and molecular dynamics simulations. Three types of DNA G-quadruplexes were utilized, which comprise parallel, antiparallel, and mixed hybrid topologies. The porphyrin hybrids have a preference to bind with parallel and mixed hybrid structures compared to the antiparallel structure. This preference arises from the end stacking of porphyrin moiety following G-stem and loop binding of anthraquinone tail, which is not found in the antiparallel due to the presence of diagonal and lateral loops that crowd the G-quartet. The binding to the antiparallel, instead, occurred with poorer affinity through both the loop and wide groove. All sites of porphyrin binding were confirmed by 6 ns molecular dynamics simulation, as well as by the negative value of the total binding free energies that were calculated using the MMPBSA method. Free energy analysis shows that the favorable contribution came from the electrostatic term, which supposedly originated from the interaction of either cationic pyridinium, pyrazole, or imidazole groups and the anionic phosphate backbone, and also from the van der Waals energy, which primarily contributed through end stacking interaction.     

Influence of divalent magnesium ion on DNA: molecular dynamics simulation studies

A large amount of experimental evidence is available on the effect of magnesium ions on the structure and stability of DNA double helix. Less is known, however, on how these ions affect the stability and dynamics of the molecule. The static time average pictures from X-ray structures or the quantum chemical energy minimized structures lack understanding of the dynamic DNA–ion interaction. The present work addresses these questions by molecular dynamics simulation studies on two DNA duplexes and their interaction with magnesium ions. Results show typical B-DNA character with occasional excursions to deviated states. We detected expected stability of the duplexes in terms of backbone conformations and base pair parameter by the CHARMM-27 force field. Ion environment analysis shows that Mg²⁺ retains the coordination sphere throughout the simulation with a preference for major groove over minor. An extensive analysis of the influence of the Mg²⁺ ion shows no evidence of the popular predictions of groove width narrowing by dipositive metal ion. The major groove atoms show higher occupancy and residence time compared to minor groove for magnesium, where no such distinction is found for the charge neutralizing Na⁺ ions. The determining factor of Mg²⁺ ion’s choice in DNA binding site evolves as the steric hindrance faced by the bulky hexahydrated cation where wider major groove gets the preference. We have shown that in case of binding of Mg²⁺ to DNA non electrostatic contributions play a major role.

An animated Interactive 3D Complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:5     

Binding site identification of anticancer drug gefitinib to HSA and DNA in the presence of five different probes

This study was carried out to evaluate the binding interaction of gefitinib (GEF) with human serum albumin (HSA) and calf thymus DNA (ct-DNA) using fluorescence, UV–Visible, zeta potential measurements and molecular docking methods in order to understand its pharmacokinetic mechanism. By increasing the temperature, a steady decrease in Stern–Volmer quenching constants was observed for HSA binding properties; this indicates a static type of fluorescence quenching. Negative values were calculated for Gibbs free energy (ΔG), enthalpy (ΔH), and entropy (ΔS) changes, indicating that the reaction is spontaneous and enthalpy-driven. Probe competitive experimental results showed that GEF contains the same binding site as warfarin and are consistent with modeling results. The zeta potential of the HSA increased with increasing GEF, which represents the presence of electrostatic interactions in the system. DNA binding properties were investigated in the presence of three probes. The experimental results showed that by increasing GEF to DNA-AO (acridine-orange) and DNA-MB (methylene-blue) system, the fluorescence intensity and absorbance spectra had no considerable change. Furthermore, with the addition of GEF to DNA, the zeta potential decreased gradually, indicating that the hydrophobic interaction between the GEF and the bases of DNA is the major factor. Thus, GEF can bind to DNA via a groove binding mode. It was also found that GEF entered into the minor groove in the A–T rich region of DNA fragment and bind via van der-Waals forces and three H-bond with double strands of DNA. This is in good agreement with experimental results.     

Energetic contributions of amino acid residues and its cross-talk to delineate ligand-binding mechanism

Receptor-based QSAR approaches can enumerate the energetic contributions of amino acid residues toward ligand binding only when experimental binding affinity is associated. The structural data of protein-ligand complexes are witnessing a tremendous growth in the Protein Data Bank deposited with a few entries on binding affinity. We present here a new approach to compute the E nergetic CONT ributions of A mino acid residues and its possible C ross-T alk (ECONTACT) to study ligand binding using per-residue energy decomposition, molecular dynamics simulations and rescoring method without the need for experimental binding affinity. This approach recognizes potential cross-talks among amino acid residues imparting a nonadditive effect to the binding affinity with evidence of correlative motions in the dynamics simulations. The protein-ligand interaction energies deduced from multiple structures are decomposed into per-residue energy terms, which are employed as variables to principal component analysis and generated cross-terms. Out of 16 cross-talks derived from eight datasets of protein-ligand systems, the ECONTACT approach is able to associate 10 potential cross-talks with site-directed mutagenesis, free energy, and dynamics simulations data strongly. We modeled these key determinants of ligand binding using joint probability density function (jPDF) to identify cross-talks in protein structures. The top two cross-talks identified by ECONTACT approach corroborated with the experimental findings. Furthermore, virtual screening exercise using ECONTACT models better discriminated known inhibitors from decoy molecules. This approach proposes the jPDF metric to estimate the probability of observing cross-talks in any protein-ligand complex. The source code and related resources to perform ECONTACT modeling is available freely at https://www.gujaratuniversity.ac.in/econtact /.     

Structure-specific binding of [Co(phen)2(HPIP)] to a DNA duplex containing sheared G:A mismatch base pairs

The binding of a Co(III) complex to the decanucleotide d(CCGAATGAGG)₂ containing two pairs of G:A mismatches was studied by 2D-NMR, UV absorption, and molecular modeling. NMR investigations indicate that racemic [Co(phen)₂(HPIP)]Cl₃ [HPIP = 2-(2-hydroxyphenyl) imidazo [4,5-f][1,10] phenanthroline] binds the decanucleotide by intercalation: the HPIP ligand selectively inserts between the stacked bases from the minor groove at the terminal regions and from the major groove at the sheared region. Further, molecular modeling revealed that the recognition shows strong enantioselectivity: the Λ-isomer preferentially intercalates into the T₆G₇:A₅A₄ region from the DNA major groove, while Δ-isomer favors the terminal C₁C₂:G₁₀G₉ region and intercalates from the minor groove. Detailed energy analysis suggests that the steric interaction, especially the electrostatic effect, is the primary determinants of the recognition event. Melting experiments indicate that binding stabilizes the DNA duplex and increases the melting temperature by 9.5 °C. The intrinsic binding constant of the complex to the mismatched duplex was determined to be 3.5 × 105 M⁻¹.     

Novel Hoogsteen-like bases for configurational recognition of the T-A base pair by DNA triplex formation

Effective sequence-specific recognition of duplex DNA is possible by triplex formation with natural oligonucleotides via Hoogsteen H-bonding. However, triplex formation is in practice limited to pyrimidine oligonucleotides binding duplex A-T or G-C base-pair DNA sequences specifically at homopurine sites in the major groove as T·A-T and C⁺·G-C triplets. Here we report the successful modeling of novel unnatural nucleosides that recognize the T-A DNA base pair by Hoogsteen interaction. Since the DNA triplex can be considered to assume an A-type or B-type conformation, these novel Hoogsteen nucleotides are tested within model A-type and B-type conformation triplex structures. A triplet consisting of the T-A base pair and one of the novel Hoogsteen nucleotides replaces the central T·A-T triplet in the triplex using the same deoxyribose-phosphodiester and base-deoxyribose dihedral angle configuration. The entire triplex is energy minimized and the presence of any structural or energetic perturbations due to the central triplet is assessed with respect to the unmodified energy-minimized (T·A-T)₁₁ proposed starting structures. Incorporation of these novel triplets into both A-type and B-type natural triplex structures provokes minimal change in the configuration of the central and adjacent triplets. The plan is to produce a series of Hoogsteen-like bases that preferentially bind the T-A major groove in either an A-type or B-type conformation. Selective recognition of the T-A major groove with respect to the G-C major groove, which presents similar keto and amine placement, is also assessed with configurational preference. Evaluation of the triplex solution structure by using these unnatural bases as binding conformational probes is a prerequisite to the further design of triplet forming bases. © 1996 John Wiley & Sons, Inc.     

Alteration of the groove width of DNA induced by the multimodal hydrogen bonding of denaturants with DNA bases in its grooves affects their stability

BackgroundDenaturants, namely, urea and guanidinium chloride (GdmCl) affect the stability as well as structure of DNA. Critical assessment of the role of hydrogen bonding of these denaturants with the different regions of DNA is essential in terms of its stability and structural aspect. However, the understanding of the mechanistic aspects of structural change of DNA induced by the denaturants is not yet well understood.MethodsIn this study, various spectroscopic along with molecular dynamics (MD) simulation techniques were employed to understand the role of hydrogen bonding of these denaturants with DNA bases in their stability and structural change.Results and conclusionIt has been found that both, GdmCl and urea intrude into groove region of DNA by striping surrounding water. The hydrogen bonding pattern of Gdm⁺ and urea with DNA bases in its groove region is multimodal and distinctly different from each other. The interaction of GdmCl with DNA is stabilized by electrostatic interaction whereas electrostatic and Lennard-Jones interactions both contribute for urea. Gdm⁺ forms direct hydrogen bond with the bases in the minor groove of DNA whereas direct and water assisted hydrogen bond takes place with urea. The hydrogen bond formed between Gdm⁺ with bases in the groove region of DNA is stronger than urea due to strong electrostatic interaction along with less self-aggregation of Gdm⁺ than urea. The distinct hydrogen bonding capability of Gdm⁺ and urea with DNA bases in its groove region affects its width differently. The interaction of Gdm⁺ decreases the width of the minor and major groove which probably increases the strength of hydrogen bond between the Watson-Crick base pairs of DNA leading to its stability. In contrast, the interaction of urea does not affect much to the width of the grooves except the marginal increase in the minor groove width which probably decreases the strength of hydrogen bond between Watson Crick base pairs leading to the destabilization of DNA.General significanceOur study clearly depicts the role of hydrogen bonding between DNA bases and denaturants in their stability and structural change which can be used further for designing of the guanidinium based drug molecules.     

Molecular modeling and biophysical analysis of the c-MYC NHE-III<Subscript>1</Subscript> silencer element

G-Quadruplex and i-Motif-forming sequences in the promoter regions of several oncogenes show promise as targets for the regulation of oncogenes. In this study, molecular models were created for the c-MYC NHE-III₁ (nuclease hypersensitivity element III₁) from two 39-base complementary sequences. The NHE modeled here consists of single folded conformers of the polypurine intramolecular G-Quadruplex and the polypyrimidine intramolecular i-Motif structures, flanked by short duplex DNA sequences. The G-Quadruplex was based on published NMR structural data for the c-MYC 1:2:1 loop isomer. The i-Motif structure is theoretical (with five cytosine–cytosine pairs), where the central intercalated cytosine core interactions are based on NMR structural data obtained for a tetramolecular [d(A₂C₄)₄] model i-Motif. The loop structures are in silico predictions of the c-MYC i-motif loops. The porphyrin meso-tetra(N-methyl-4-pyridyl)porphine (TMPyP4), as well as the ortho and meta analogs TMPyP2 and TMPyP3, were docked to six different locations in the complete c-MYC NHE. Comparisons are made for drug binding to the NHE and the isolated G-Quadruplex and i-Motif structures. NHE models both with and without bound cationic porphyrin were simulated for 100 ps using molecular dynamics techniques, and the non-bonded interaction energies between the DNA and porphyrins calculated for all of the docking interactions. Figure Molecular models of the average structure of the final 20 ps of the molecular dynamics simulation of the c-MYC NHE-III₁ (nuclease hypersensitivity element III₁) “silencer” element. The G-Quadruplex structure is at the top-center, and the i-Motif is at the bottom-center of each picture. a “Rotation #1” of the G-Quadruplex, with the T15 loop at the top and rear and the G19/A20 loop at the top and front of the picture. b “Rotation #2” of the G-Quadruplex, with the T15 loop at the top and front of the image, and the G19/A20 loop at the front and adjacent to the G-Quadruplex/i-Motif interface     

Role of Mg2+ in Chromomycin A3 – DNA Interaction: A Molecular Modeling Study

Chromomycin A₃ (CHR) is an antitumor antibiotic that inhibits macromolecular biosynthesis by reversibly binding to double stranded DNA via the minor groove, with GC-base specificity. At and above physiological pH when CHR is anionic, interaction of CHR with DNA requires the presence of divalent metal ions like Mg²⁺. However, at acidic pHthe molecule is neutral and it binds DNA even in absence of Mg²⁺. Molecular dynamics simulation studies at 300K of neutral CHR and 1:1 CHR:Mg²⁺ complexes formed at pH 5.2 and 8.0 show that hydrophobicity of CHR:Mg²⁺ complex formed with the neutral drug is greater than that of the two other species. Interactions of CHR with DNA in presence and absence of Mg²⁺ have been studied by simulated annealing to understand the role of Mg²⁺ in the DNA binding potential of CHR. This shows that the antibiotic has the structural potential to bind to DNA even in the absence of metal ion. Evaluation of the direct interaction energy between the ligand and DNA does not explain the observed GC-base specificity of the antibiotic. When energy contributions from structural alteration of the interacting ligand and DNA as a sequel to complex formation are taken into account, atrue picture of the theoretical binding propensity emerges. This implies that DNA and/or the ligand undergo significant structural alterations during the process of association, particularly in presence of Mg²⁺. Accessible surface area calculations give idea about the entropy contribution to the binding free energy which is found to be different depending upon the presence and absence of Mg²⁺.     

On the importance of van der Waals interaction in the groove binding of DNA with ligands: restrained molecular dynamics study

Comparison of interaction energy between an oligonucleotide and a DNA-binding ligand in the minor and major groove modes was made by use of restrained molecular dynamics. Distortion in DNA was found for the major groove mode whereas less significant changes for both ligand and DNA were detected for the minor groove binding after molecular dynamics simulation. The conformation of the ligand obtained from the major groove mode resembles that computed with the ligand soaked in water. The van der Waals contact energy was found to be as significant as electrostatic energy and more important for difference in binding energy between these two binding modes. The importance of van der Waals force in groove binding was supported by computations on the complex formed by the repressor peptide fragment from the bacteriophage 434 and its operator oligonucleotide.     

Novel Hoogsteen-like Bases for Recognition of the C-G Base Pair by DNA Triplex Formation

Effective sequence-specific recognition of duplex DNA is possible by triplex formation with natural oligonucleotides via Hoogsteen H-bonding. However, triplex formation is in practice limited to pyrimidine oligonucleotides that bind duplex A-T or G-C base pair DNA sequences specifically at homopurine sites in the major groove as T·A-T and C⁺ ·G-C triplets. Here we report the successful modelling of novel unnatural nucleosides that recognize the C-G DNA base pair by Hoogsteen-like major groove interaction. These novel Hoogsteen nucleotides are examined within model A-type and B-type conformation triplex structures since the DNA triplex can be considered to incorporate A-type and/or B-type configurational properties. Using the same deoxyribose-phosphodiester and base-deoxyribose dihedral angle configuration, a triplet comprised of a C-G base pair and the novel Hoogsteen nucleotide, Y2, replaces the central T·A-T triplet in the triplex. The presence of any structural or energetic perturbations due to the central triplet in the energy-minimized triplex is assessed with respect to the unmodified energy minimized (T·A-T)₁₁ starting structures. Incorporation of this novel triplet into both A-type and B-type natural triplex structures provokes minimal change in the configuration of the central and adjacent triplets.     

Towards an Understanding of DNA Recognition by the Methyl-CpG Binding Domain 1

Abstract

CpG methylation determines a variety of biological functions of DNA. The methylation signal is interpreted by proteins containing a methyl-CpG binding domain (MBDs). Based on the NMR structure of MBD1 complexed with methylated DNA we analysed the recognition mode by means of molecular dynamics simulations.

As the protein is monomeric and recognizes a symmetrically methylated CpG step, the recognition mode is an asymmetric one. We find that the two methyl groups do not contribute equally to the binding energy. One methyl group is associated with the major part of the binding energy and the other one nearly does not contribute at all. The contribution of the two cytosine methyl groups to binding energy is calculated to be ?3.6 kcal/mol. This implies a contribution of greater than two orders of magnitude to the binding constant. The conserved amino acid Asp32 is known to be essential for DNA binding by MBD1, but so far no direct contact with DNA has been observed. We detected a direct DNA base contact to Asp32. This could be the main reason for the importance of this amino acid. MBD contacts DNA exclusively in the major groove, the minor groove is reserved for histone contacts. We found a deformation of the minor groove shape due to complexation by MBD1, which indicates an information transfer between the major and the minor groove.