In silico footprinting of ligands binding to the minor groove of DNA

The sequence selectivity of small molecules binding to the minor groove of DNA can be predicted by "in silico footprinting". Any potential ligand can be docked in the minor groove and then moved along it using simple simulation techniques. By applying a simple scoring function to the trajectory after energy minimization, the preferred binding site can be identified. We show application to all known noncovalent binding modes, namely 1:1 ligand:DNA binding (including hairpin ligands) and 2:1 side-by-side binding, with various DNA base pair sequences and show excellent agreement with experimental results from X-ray crystallography, NMR, and gel-based footprinting.     

诺氟沙星-DNA复合物的分子动力学模拟

采用分子模建的方法构建了诺氟沙星-DNA复合物的初始结构, 通过2 ns的分子动力学(MD)模拟研究表明: 诺氟沙星能够和双螺旋d[ATATCGATAT]₂形成稳定的复合物, 药物分子可紧密结合在DNA的小沟区域, 并且能够与DNA的鸟嘌呤碱基形成两个稳定的氢键. 在分子水平上提供了诺氟沙星直接与双螺旋DNA相互作用的结构及复合物的动态变化情况.     

Water-mediated binding of agents that target the DNA minor groove

Small molecule complexes with DNA that incorporate linking water molecules are rare, and the DB921-DNA complex has provided a unique and well-defined system for analysis of water-mediated binding in the context of a DNA complex. DB921 has a benzimidazole-biphenyl system with terminal amidines that results in a linear conformation that does not possess the appropriate radius of curvature to match the minor groove shape and represents a new paradigm that does not fit the classical model of minor groove interactions. To better understand the role of the bound water molecule observed in the X-ray crystal structure of the DB921 complex, synthetic modifications have been made in the DB921 structure, and the interactions of the new compounds with DNA AT sites have been evaluated with an array of methods, including DNase I footprinting, biosensor-surface plasmon resonance, isothermal titration microcalorimetry, and circular dichroism. The interaction of a key compound, which has the amidine at the phenyl shifted from the para position in DB921 to the meta position, has also been examined by X-ray crystallography. The detailed structural, thermodynamic, and kinetic results provide valuable new information for incorporation of water molecules in the design of new lead scaffolds for targeting DNA in chemical biology and therapeutic applications.     

DNA小沟结合药物DB818的分子动力学模拟和结合自由能分析

采用分子动力学模拟了DNA小沟结合芳香二脒药物DB818形成的复合物. 通过5 ns的模拟研究表明: DB818药物分子可紧密结合在DNA的AATTC小沟区域, 和双螺旋d[CGCGAATTCGCG]2形成稳定的复合物. 由于噻吩硫原子的弱电负性, 使DB818能够以更大的伸展程度与DNA的小沟结合, 形成更强的结合力. DB818苯并咪唑的氮原子能够与DNA 7位和19位T碱基上的氧原子形成两个稳定的氢键, 同时, DB818末端氨基氮原子分别与DNA 的20位T碱基的氧原子和9位C碱基的氧原子形成两个氢键. 另外, 运用MM_PBSA方法计算了DB293-DNA和DB818-DNA复合物的结合自由能, 计算结合能与实验值能较好的吻合, 通过比较其结合自由能, 从热力学能量角度说明了DB818有较大的熵值与较小的焓值贡献, 从而与DNA小沟结合的结合力比DB293强. 本文在分子水平上提供了DB818直接与双螺旋DNA相互作用的结构及复合物的动态变化情况, 为设计出更高生物活性的DNA小沟结合剂提供一定的理论依据.     

Direct observation of essential DNA dynamics: melting and reformation of the DNA minor groove

The dynamics of bound water and ions present in the minor groove of a dodecamer DNA has been decoupled from that of the long-range twisting/bending of the DNA backbone, using the minor groove binder Hoechst 33258 as a fluorescence reporter in the picosecond-resolved time window. The bound water and ions are essential structural components of the minor groove and are destroyed with the destruction of the minor groove when the dodecamer melts at high temperatures and reforms on subsequent cooling of the melted DNA. The melting and rehybridization of the DNA has been monitored by the changes in secondary structure using circular dichroism (CD) spectroscopy. The change in the relaxation dynamics of the DNA has been studied with picosecond resolution at different temperatures, following the temperature-dependent melting and rehybridization profile of the dodecamer, using time-resolved emission spectra (TRES). At room temperature, the relaxation dynamics of DNA is governed by a 40 ps (30%) and a 12.3 ns (70%) component. The dynamics of bound water and ions present in the minor groove is characterized by the 40 ps component in the relaxation dynamics of the probe bound in the minor groove of the dodecamer DNA. Analyses of the TRES taken at different temperatures show that the contribution of this component decreases and ultimately vanishes with the destruction of the minor groove and reappears again with the reformation of the groove. The dynamical behavior of bound water molecules and ions of a genomic DNA (from salmon testes) at different temperatures is also found to be consistent with that of the dodecamer. The longer component of approximately 10 ns in the DNA dynamics is found to be associated with the long-range bending/twisting of the DNA backbone and the associated counterions. The transition from bound water to free water at the DNA surface, indicative of the change in the hydration number associated with each base pair, has also been ascertained in the case of the genomic DNA at different temperatures by employing densimetric and acoustic techniques.     

DNA小沟结合二脒:水分子的识别作用

采用分子动力学模拟了DB921-DNA复合物, 通过7 ns的模拟研究表明: DB921一端的氨基氮原子与一个水分子形成氢键, 同时, 水分子又与DNA的5位A碱基的氮原子形成一个氢键. 水分子在DB921与DNA小沟结合中起了桥连的作用, 使得直线型的芳香二脒化合物DB921通过水桥与DNA小沟结合, 水分子诱导DB921分子与DNA的小沟域构型相适应, 与DNA小沟域的AATTC碱基有较强的结合作用. 在分子水平上提供了DB921与双螺旋DNA相互作用的结构及复合物的动态变化情况, 指出水分子在DNA小沟结合二脒化合物中的识别作用, 为设计出更高生物活性的DNA小沟结合剂提供一定的理论依据.     

Exploring DNA groove water dynamics through hydrogen bond lifetime and orientational relaxation

Dynamics of water molecules in the grooves of DNA are of great interest both for practical (functionality of DNA) and fundamental (as examples of confined systems) interest. Here the authors employ atomistic molecular dynamics simulations to understand varying water dynamics at the minor and the major grooves of a 38 base-pair long DNA duplex in water. In order to understand and quantify the diversity in the nature of hydrogen bond due to many hydrogen bond donors and acceptors present in the four bases, they have undertaken study of hydrogen bond lifetime (HBLT) correlation functions of all the specific hydrogen bonds between the base atoms and water molecules. They find that the HBLT correlation functions are in general multiexponential, with the average lifetime depending significantly on the specificity and may thus be biologically relevant. The average hydrogen bond lifetime is longer in the minor groove than that in the major groove by almost a factor of 2. Analysis further shows that water hydrogen bonds with phosphate oxygen have substantially shorter lifetimes than those with the groove atoms. They also compute two different orientational time correlation functions (OTCFs) of the water molecules present at the major and the minor grooves and attempt to correlate OTCF with HBLT correlation function. The OTCFs in the minor groove exhibit three time scales, with the time constant of the slowest component one to two orders of magnitude longer than what is observed for bulk water. A slow component is also present for the major groove water but with shorter time constant. Interestingly, correlation between reformations allowed HBLT correlation function [C(HB)(t)] and the OTCF markedly deviates from each other in the grooves, indicating enhanced rigidity of water molecules in the grooves.     

Role of aromacity, ring strain, and stereochemistry of selected disulfides and their congeners in the oxidative linkage of DNA strands at the major groove. A computational study

Three aromatic disulfides and their two selected congeners were studied in order to see if they could be initiators or catalysts in a process that aims to tie together the two strands of a DNA double helix. Thermodynamic functions were calculated from first principle molecular computations on deoxyribose models and the disulfide congener. The gas reaction of one of the disulfides turned out to be the closest to the thermoneutral process. Solvation study caused a relatively minor energetic modification. The disulfides were subjected to docking studies, and the stabilities of the complexes obtained were further analyzed by molecular dynamics simulations. Only one of the two atropic disulfide enantiomers was the matching partner of the chiral DNA double helix, and therefore this remained in the major groove of the DNA. The other enantiomer dissociated from the major groove as it was incapable of forming a stable complex with the chiral double helix. A mechanistic scheme for the reaction has been investigated.     

Elucidating the intermolecular interactions within a desolvated protein-ligand complex. An experimental and computational study

The first detailed study of the intermolecular hydrogens bonds (H-bonds) within a desolvated, noncovalent protein-ligand complex is reported. Using both experimental and computational methods, the intermolecular H-bonds stabilizing protonated and deprotonated ions of a complex composed of a single chain fragment (scFv) of a monoclonal antibody and its native trisaccharide ligand, alphaGal[alphaAbe] alphaMan (1), are characterized. Using the blackbody infrared radiative dissociation-functional group replacement (BIRD/FGR) technique, three H-bond donor-acceptor pairs within the gaseous (scFv + 1)n+ ions are identified and quantified. Additional sites of interaction on the protein and ligand, for which the binding partner could not be elucidated, are also identified. Comparison of the gas-phase interaction maps with the crystal structure suggests that at least two of the specific H-bonds are conserved upon transfer of the complex from solution to the gas phase by electrospray ionization. However, new (nonspecific) interactions can also form in the gas phase. Notably, the nature and strength of the intermolecular interactions can vary significantly with charge state, and striking differences in the structures of the (scFv + 1)n+ and (scFv + 1)n- ions are evident. Intermolecular H-bonds are also identified from molecular dynamics (MD) simulations performed at the +8 and -8 charge states. Agreement is found for a majority of intermolecular interactions predicted for the (scFv + 1)8+ ion by the MD simulation and BIRD/FGR method; the agreement is less favorable in the case of the (scFv + 1)8- ion. However, both the computational and experimental results point to structural differences between the +8 and -8 ions. The computational results also provide insights into the structural changes that accompany the loss of interfacial waters from the complex.     

DNA与二脒DB293复合物的分子动力学模拟

采用分子动力学模拟了DNA小沟与芳香二脒化合物DB293结合形成的复合物,通过5ns的模拟研究表明,DB293分子可紧密结合在DNA的AATT小沟区域,和双螺旋d[CGCGAATTCGCG]2形成稳定的复合物。DB293苯并咪唑的氮原子N2能够与DNA胸腺嘧啶碱基T7的O2原子和T19的O2原子形成两个较强的氢键,同时,其末端氨基的N3原子和T20的O2原子形成一个较弱的氢键。本文在分子水平上提供了DB293直接与双螺旋DNA相互作用的结构及复合物的动态变化情况,为设计出更高活性的芳香二脒类DNA小沟结合剂提供一定的理论依据。     

Competitive Na(+) and Rb(+) binding in the minor groove of DNA

Sequence-dependent coordination of alkali ions to the nucleotide bases in the minor groove of AT-tract B-DNA has recently been inferred from X-ray crystallography, solution NMR and computer simulations. Here, we present new (23)Na and (87)Rb magnetic relaxation dispersion (MRD) data that demonstrate competitive and long-lived binding of Na(+) and Rb(+) ions in the minor groove of the B-DNA duplex [d(CGCGAATTCGCG)](2). The Na(+)/Rb(+) selectivity of the minor groove is found to be weak, consistent with local structural flexibility. The ion occupancies derived from the MRD data are substantially higher than previously reported, suggesting that groove-bound ions significantly influence the energetics and structural polymorphism of DNA in vivo. For example, in the presence of 0.20 M Na(+) and 0.56 M Rb(+) at 4 degrees C, the ApT site in the minor groove is occupied by a Rb(+) ion, a Na(+) ion, or a water molecule 40, 10, and 50% of the time, respectively. In the absence of Rb(+), the Na(+) occupancy increases to 50%. At 4 degrees C, the mean residence time of groove-bound ions is 0.2 +/- 0.1 micros for Rb(+) and 10 ns to 100 micros for Na(+). A shorter correlation time of 2 ns is attributed to counterions bridging cross-strand phosphate groups.     

Molecular modeling on the recognition of DNA sequence and conformational repair of sheared DNA by novel chiral metal complex D, L-[Co(phen)_2hpip]~(3 )

Recent studies revealed that DNA, once considered as a very stable macromolecular, is rather unstable. Familiar factors, like heavy metal, microbe, high fre-quency electromagnetic radiation and so on, could easily damage the structure of DNA in different …     

Influence of the dynamic positions of cations on the structure of the DNA minor groove: sequence-dependent effects

Different models for minor groove structures predict that the conformation is essentially fixed by sequence and has an influence on local ion distribution or alternatively that temporal positions of ions around the minor groove can affect the structure if they neutralize cross-strand phosphate charges. Our previous studies show that the minor groove in an AATT dodecamer responds to local sodium ion positions and is narrow when ions neutralize cross-strand phosphate-phosphate charges [J. Am. Chem. Soc. 2000, 122, 10513-10520]. Previous results from a number of laboratories have shown that G-tracts often have a wider minor groove than A-tracts, but they do not indicate whether this is due to reduced flexibility or differences in ion interactions. We have undertaken a molecular dynamics study of a d(TATAGGCCTATA) duplex to answer this question. The results show that the G-tract has the same amplitude of minor groove fluctuations as the A-tract sequence but that it has fewer ion interactions that neutralize cross-strand phosphate charges. These results demonstrate that differences in time-average groove width between A- and G-tracts are due to differences in ion interactions at the minor groove. When ions neutralize the cross-strand phosphates, the minor groove is narrow. When there are no neutralizing ion interactions, the minor groove is wide. The population of structures with no ion interactions is larger with the GGCC than with the AATT duplex, and GGCC has a wider time-average minor groove in agreement with experiment.     

Energetics of ligand binding to the DNA minor groove

In the present work the decomposition of the total Gibbs free energy of ligand-DNA binding onto various physical terms was accomplished for the group of nine DNA minor groove binders (MGB ligands) differing in both structure and charge state. The decomposition protocol includes the analysis of the most complete set of physical factors known to contribute to the complexation process, viz. the net change in the number of degrees of freedom (translational, rotational, vibrations of the chemical bonds and vibrations of the ligand as a whole within the binding site), the conformational entropy, van der Waals, electrostatic and hydrophobic interactions, the polyelectrolyte contribution and the net effect of changes in the number of hydrogen bonds. All of these processes are further decomposed into the interaction with the solvent and the interaction of the ligand with DNA. The principal outcome of the decomposition is the possibility of performing a comparative analysis of the energetic contribution of various physical terms and provide an answer to the question concerning what physical factors stabilize or destabilize the complexes of MGB ligands with DNA.     

Protein–ligand interfaces are polarized: discovery of a strong trend for intermolecular hydrogen bonds to favor donors on the protein side with implications for predicting and designing ligand complexes

Understanding how proteins encode ligand specificity is fascinating and similar in importance to deciphering the genetic code. For protein–ligand recognition, the combination of an almost infinite variety of interfacial shapes and patterns of chemical groups makes the problem especially challenging. Here we analyze data across non-homologous proteins in complex with small biological ligands to address observations made in our inhibitor discovery projects: that proteins favor donating H-bonds to ligands and avoid using groups with both H-bond donor and acceptor capacity. The resulting clear and significant chemical group matching preferences elucidate the code for protein-native ligand binding, similar to the dominant patterns found in nucleic acid base-pairing. On average, 90% of the keto and carboxylate oxygens occurring in the biological ligands formed direct H-bonds to the protein. A two-fold preference was found for protein atoms to act as H-bond donors and ligand atoms to act as acceptors, and 76% of all intermolecular H-bonds involved an amine donor. Together, the tight chemical and geometric constraints associated with satisfying donor groups generate a hydrogen-bonding lock that can be matched only by ligands bearing the right acceptor-rich key. Measuring an index of H-bond preference based on the observed chemical trends proved sufficient to predict other protein–ligand complexes and can be used to guide molecular design. The resulting Hbind and Protein Recognition Index software packages are being made available for rigorously defining intermolecular H-bonds and measuring the extent to which H-bonding patterns in a given complex match the preference key.     

Molecular modeling on the recognition of DNA sequence and conformational repair of sheared DNA by novel chiral metal complex D,L-[Co(phen)2hpip]3+

A study on the recognition of DNA sequence and conformational repair of sheared DNA by Novel Chiral Metal complex D,L-[Co(phen)₂hpip]³⁺ (phen=1,10 phenanthroline, hpip=2-[2-hydroxyphenyl] imidazole [4,5-f][1,10] phenanthroline) is carried out with molecular simulations. The results reveal that two isomers of the complex could both recognize the normal DNA in the minor groove orientation, while recognize the sheared DNA in the major groove orientation and both isomers could convert the conformation of mismatched bases from sheared form to parallel form. Further analysis shows that the steric details of complex’s intercalation to base stack determine the results of recognition, which is induced by the steric collision among ancillary ligand phen, bases and DNA backbone, and by the steric crowding occurring in the process of structural expansion of bases and DNA backbone. Detailed analysis reveals that the conformational repair of mismatched bases relates not only to the steric interactions, but also the π-π stack among normal bases, mismatched bases and hpip ligand.     

Circular dichroism spectroscopic studies reveal pH dependent binding of curcumin in the minor groove of natural and synthetic nucleic acids

For the first time, an interaction between the non-toxic, cancer chemopreventive agent curcumin and both natural and synthetic DNA duplexes has been demonstrated by using circular dichroism (CD) and absorption spectroscopy techniques. Upon addition of curcumin to calf thymus DNA, poly(dG-dC).poly(dG-dC) and poly(dA-dT).poly(dA-dT) solutions, an intense positive induced CD band centered around 460-470 nm was observed depending on the actual pH and Na+ ion concentration of the medium; no CD signal was obtained, however, with single stranded poly(dC). Interaction of curcumin with calf thymus DNA was observed already at pH 6.5 in contrast with poly(dG-dC).poly(dG-dC) which induces no extrinsic Cotton effect above a pH value of 5. The protonated, Hoogsteen base-paired structure of poly(dG-dC).poly(dG-dC) is necessary for curcumin binding while the alternating AT-rich polymer formed complexes with curcumin only at certain Na+ concentrations. Evaluation of the spectral data and molecular modeling calculations suggested that curcumin, this dietary polyphenolic compound binds in the minor groove of the double helix. The mechanism of the induced CD activity, the effects of the pH and Na+ ions on the ligand binding and conformation of the double helix are discussed in detail. As well as being an essentially new phenolic minor groove binder agent curcumin is also a promising molecular probe to study biologically important, pH and cation induced conformational polymorphisms of nucleic acids.