A detailed binding free energy study of 2:1 ligand-DNA complex formation by experiment and simulation

Structure-based drug design of DNA minor groove binders and evaluation of their antibacterial and anticancer properties

Antimicrobial and chemotherapy resistance are escalating medical problem of paramount importance. Yet, research for novel antimicrobial and anticancer agents remains lagging behind. With their reported medical applications, DNA minor groove binders (MGBs) are worthy of exploration. In this study, the approach of structure-based drug design was implemented to generate 11 MGB compounds including a novel class of bioactive alkyne-linked MGBs. The NCI screening protocol was utilized to evaluate the antitumor activity of the target MGBs. Furthermore, a variety of bactericidal, cytopathogenicity, MIC90, and cytotoxicity assays were carried out using these MGBs against 6 medically relevant bacteria: Salmonella enterica, Escherichia coli, Serratia marcescens, Bacillus cereus, Streptococcus pneumoniae and Streptococcus pyogenes. Moreover, molecular docking, molecular dynamic simulations, DNA melting, and isothermal titration calorimetry (ITC) analyses were utilized to explore the binding mode and interactions between the most potent MGBs and the DNA duplex d(CGACTAGTCG)2. NCI results showed that alkyne-linked MGBs (26 & 28) displayed the most significant growth inhibition among the NCI-60 panel. In addition, compounds MGB3, MGB4, MGB28, and MGB32 showed significant bactericidal effects, inhibited B. cereus and S. enterica-mediated cytopathogenicity, and exhibited low cytotoxicity. MGB28 and MGB32 demonstrated significant inhibition of S. pyogenes, whereas MGB28 notably inhibited S. marcescens and all four minor groove binders significantly inhibited B. cereus. The ability of these compounds to bind with DNA and distort its groove dimensions provides the molecular basis for the allosteric perturbation of proteins-DNA interactions by MGBs. This study shed light on the mechanism of action of MGBs and revealed the important structural features for their antitumor and antibacterial activities, which are important to guide future development of MGB derivatives as novel antibacterial and anticancer agents.

Interaction of an anticancer benzopyrane derivative with DNA: Biophysical, biochemical, and molecular modeling studies

BACKGROUND

SIMR1281 is a potent anticancer lead candidate with multi- target activity against several proteins; however, its mechanism of action at the molecular level is not fully understood. Revealing the mechanism and the origin of multitarget activity is important for the rational identification and optimization of multitarget drugs.

METHODS

We have used a variety of biophysical (circular dichroism, isothermal titration calorimetry, viscosity, and UV DNA melting), biochemical (topoisomerase I & II assays) and computational (molecular docking and MD simulations) methods to study the interaction of SIMR1281 with duplex DNA structures.

RESULTS

The biophysical results revealed that SIMR1281 binds to dsDNA via an intercalation-binding mode with an average binding constant of 3.1 × 10⁶ M^-1. This binding mode was confirmed by the topoisomerases' inhibition assays and molecular modeling simulations, which showed the intercalation of the benzopyrane moiety between DNA base pairs, while the remaining moieties (thiazole and phenyl rings) sit in the minor groove and interact with the flanking base pairs adjacent to the intercalation site.

CONCLUSIONS

The DNA binding characteristics of SIMR1281, which can disrupt/inhibit DNA function as confirmed by the topoisomerases' inhibition assays, indicate that the observed multi-target activity might originate from ligand intervention at nucleic acids level rather than due to direct interactions with multiple biological targets at the protein level.

GENERAL SIGNIFICANCE

The findings of this study could be helpful to guide future optimization of benzopyrane-based ligands for therapeutic purposes.

Efficient DNA-ligand interaction framework using fuzzy C-means clustering based glowworm swarm optimization (FCMGSO) method

Assessment of DNA and ligand interaction is a great challenge to the medical researchers and drug industries since the accurate mapping of DNA and ligand plays an important role in associating drugs for suitable diseases. The primary objective of this research work is to develop an efficient model for predicting the best DNA and Ligand mapping. In this research work, 500 instances of DNA and drugs used for cancer and non-cancer diseases from the National Centre for Biotechnology Information (NCBI) were considered for analysis. Binding energy is one of the important measures to predict and finalize the best DNA and ligand interaction. Existing methods used for the docking process such as Simulated Annealing (SA), Lamarckian Genetic Algorithm (LGA), Genetic Clustering (GC), Fuzzy C-means clustering (FCM), and Genetic Clustering with Multi swarm Optimization (GCMSO) were applied for all 500 instances. These algorithms failed to produce better binding energy due to a lack of optimization in the existing approaches. Optimization methods play a major role in predicting accurate DNA ligand docking. Hence, this research proposes an efficient architecture using Fuzzy C-Means Clustering with Glowworm Swarm (FCMGSO) optimization method for accurate analysis of the DNA-ligand docking process. Results are proving that the proposed FCMGSO algorithm shows less binding energy than other existing methods in all instances of samples considered from the NCBI dataset.Communicated by Ramaswamy H. Sarma.

Multitargeted anti-infective drugs: resilience to resistance in the antimicrobial resistance era

The standard drug discovery paradigm of single molecule – single biological target – single biological effect is perhaps particularly unsuitable for anti-infective drug discovery. This is due to the rapid evolution of resistance likely to be observed with single target drugs. Multitargeted anti-infective drugs are likely to be superior due to their lower susceptibility to target-related resistance mechanisms. Strathclyde minor groove binders are a class of compounds which have been developed by adopting the multitargeted anti-infective drugs paradigm, and their effectiveness against a wide range of pathogenic organisms is discussed. The renaming of this class to Strathclyde nucleic acid binders is also presented due to their likely targets including both DNA and RNA.

RLDOCK method for predicting RNA-small molecule binding modes

RNA molecules play critical roles in cellular functions at the level of gene expression and regulation. The intricate 3D structures and the functional roles of RNAs make RNA molecules ideal targets for therapeutic drugs. The rational design of RNA-targeted drug requires accurate modeling of RNA-ligand interactions. Recently a new computational tool, RLDOCK, was developed to predict ligand binding sites and binding poses. Using an iterative multiscale sampling and search algorithm and a energy-based evaluation of ligand poses, the method enables efficient and accurate predictions for RNA-ligand interactions. Here we present a detailed illustration of the computational procedure for the practical implementation of the RLDOCK method. Using Flavin mononucleotide (FMN) docking to F. nucleatum FMN riboswitch as an example, we illustrate the computational protocol for RLDOCK-based prediction of RNA- ligand interactions. The RLDOCK software is freely accessible at http://https://github.com/Vfold-RNA/RLDOCK.

RLDOCK: A New Method for Predicting RNA-Ligand Interactions

The ability to accurately predict the binding site, binding pose, and binding affinity for ligand-RNA binding is important for RNA-targeted drug design. Here, we describe a new computational method, RLDOCK, for predicting the binding site and binding pose for ligand-RNA binding. By developing an energy-based scoring function, we sample exhaustively all of the possible binding sites with flexible ligand conformations for a ligand-RNA pair based on the geometric and energetic scores. The model distinguishes from other approaches in three notable features. First, the model enables exhaustive scanning of all of the possible binding sites, including multiple alternative or coexisting binding sites, for a given ligand-RNA pair. Second, the model is based on a new energy-based scoring function developed here. Third, the model employs a novel multistep screening algorithm to improve computational efficiency. Specifically, first, for each binding site, we use a gird-based energy map to rank the binding sites according to the minimum Lennard-Jones potential energy for the different ligand poses. Second, for a given selected binding site, we predict the ligand pose using a two-step algorithm. In the first step, we quickly identify the probable ligand poses using a coarse-grained simplified energy function. In the second step, for each of the probable ligand poses, we predict the ligand poses using a refined energy function. Tests of the RLDOCK for a set of 230 RNA-ligand-bound structures indicate that RLDOCK can successfully predict ligand poses for 27.8, 58.3, and 69.6% of all of the test cases with the root-mean-square deviation within 1.0, 2.0, and 3.0 Å, respectively, for the top three predicted docking poses. The computational method presented here may enable the development of a new, more comprehensive framework for the prediction of ligand-RNA binding with an ensemble of RNA conformations and the metal-ion effects.

Enzymatic degradation of extracellular DNA exposed to chlorpyrifos and chlorpyrifos-methyl in an aqueous system

The persistence of extracellular DNA (eDNA) is crucial for ensuring species diversity and ecological function in aquatic systems. However, scarce information exists about the impact of pesticides on eDNA, although they often co-exist in the aquatic environment. Using a variety of spectroscopic analyses, eDNA degradation and the associated alterations in DNA secondary structure was investigated by exposing DNase I to tested DNA in the presence of chlorpyrifos, a commonly used organophosphate pesticide. Molecular dynamics simulation was used to explore the weak interactions between the tested DNA and chlorpyrifos. The results indicated that chlorpyrifos significantly enhanced DNA degradation without affecting the enzyme activity of DNase I in an aqueous system. Spectroscopic experiments confirmed that chlorpyrifos and the analog chlorpyrifos-methyl could bind with DNA to cause the bases noncovalent stacking interaction. Molecular simulations further demonstrated that pesticide binding with DNA molecules caused widening of the DNA grooves and destruction of the hydrated layer, which enhanced DNA degradation. The findings presented herein provide novel insight into the genotoxicity and ecotoxicity of chlorpyrifos and chlorpyrifos-methyl, as well as their impacts on DNA persistence in aquatic environments.

Accurate Estimation of the Standard Binding Free Energy of Netropsin with DNA

DNA is the target of chemical compounds (drugs, pollutants, photosensitizers, etc.), which bind through non-covalent interactions. Depending on their structure and their chemical properties, DNA binders can associate to the minor or to the major groove of double-stranded DNA. They can also intercalate between two adjacent base pairs, or even replace one or two base pairs within the DNA double helix. The subsequent biological effects are strongly dependent on the architecture of the binding motif. Discriminating between the different binding patterns is of paramount importance to predict and rationalize the effect of a given compound on DNA. The structural characterization of DNA complexes remains, however, cumbersome at the experimental level. In this contribution, we employed all-atom molecular dynamics simulations to determine the standard binding free energy of DNA with netropsin, a well-characterized antiviral and antimicrobial drug, which associates to the minor groove of double-stranded DNA. To overcome the sampling limitations of classical molecular dynamics simulations, which cannot capture the large change in configurational entropy that accompanies binding, we resort to a series of potentials of mean force calculations involving a set of geometrical restraints acting on collective variables.

Resolving the Ligand-Binding Specificity in c-MYC G-Quadruplex DNA: Absolute Binding Free Energy Calculations and SPR Experiment

We report the absolute binding free energy calculation and surface plasmon resonance (SPR) experiment for ligand binding with the c-MYC G-quadruplex DNA. The unimolecular parallel DNA G-quadruplex formed in nuclease hypersensitivity element III₁ of the c-MYC gene promoter regulates the c-MYC transcription and is recognized as an emerging drug target for cancer therapy. Quindoline derivatives have been shown to stabilize the G-quadruplex and inhibit the c-MYC expression in cancer cells. NMR revealed two binding sites located at the 5' and 3' termini of the G-quadruplex. Questions about which site is more favored and the basis for the ligand-induced binding site formation remain unresolved. Here, we employ two absolute binding free energy methods, the double decoupling and the potential of mean force methods, to dissect the ligand-binding specificity in the c-MYC G-quadruplex. The calculated absolute binding free energies are in general agreement with the SPR result and suggest that quindoline has a slight preference for the 5' site. The flanking residues around the two sites undergo significant reorganization as the ligand unbinds, which provides evidence for ligand-induced binding pocket formation. The results help interpret experimental data and inform rational design of small molecules targeting the c-MYC G-quadruplex.

Experimental Probing and Molecular Dynamics Simulation of the Molecular Recognition of DNA Duplexes by the Flavonoid Luteolin

Luteolin (C₁₅H₁₀O₆) is an important flavonoid found in many fruits, plants, medicinal herbs, and vegetables exhibiting many pharmacological properties. The anticancer, antitumor, antioxidant, and anti-inflammatory activities of luteolin have been reported. The pharmacological action of small molecules is dependent upon its interaction with biomacromolecules. The interactions of small molecules with DNA play a major role in the transcription and translation process. In this work, we explored the energetic profile of DNA-luteolin interaction by isothermal titration calorimetry (ITC). The effect of temperature and salt concentration on DNA binding was examined by UV-Vis method. The mode of interaction was further probed by UV melting temperature analysis and differential scanning calorimetry. An atomic level insight on the recognition of luteolin with DNA was achieved by employing molecular dynamics (MD) simulation on luteolin in complex with AT- and GC-rich DNA sequences. AMBER force field proves to be appropriate in providing an understanding on the binding mode and specificity of luteolin with duplex DNA. MD results suggest a minor groove binding of luteolin with DNA and the binding free energy obtained is in agreement with the experimental results.

Predicting Molecular Crowding Effects in Ion-RNA Interactions

We develop a new statistical mechanical model to predict the molecular crowding effects in ion-RNA interactions. By considering discrete distributions of the crowders, the model can treat the main crowder-induced effects, such as the competition with ions for RNA binding, changes of electrostatic interaction due to crowder-induced changes in the dielectric environment, and changes in the nonpolar hydration state of the crowder-RNA system. To enhance the computational efficiency, we sample the crowder distribution using a hybrid approach: For crowders in the close vicinity of RNA surface, we sample their discrete distributions; for crowders in the bulk solvent away from the RNA surface, we use a continuous mean-field distribution for the crowders. Moreover, using the tightly bound ion (TBI) model, we account for ion fluctuation and correlation effects in the calculation for ion-RNA interactions. Applications of the model to a variety of simple RNA structures such as RNA helices show a crowder-induced increase in free energy and decrease in ion binding. Such crowding effects tend to contribute to the destabilization of RNA structure. Further analysis indicates that these effects are associated with the crowder-ion competition in RNA binding and the effective decrease in the dielectric constant. This simple ion effect model may serve as a useful framework for modeling more realistic crowders with larger, more complex RNA structures.

Thermodynamics of DNA: sensitizer recognition. Characterizing binding motifs with all-atom simulations

We report the investigation of the thermochemical properties of benzophenone interacting with B-DNA studied by all-atom molecular dynamic simulations.

Evaluation of the Gibbs Free Energy Changes and Melting Temperatures of DNA/DNA Duplexes Using Hybridization Enthalpy Calculated by Molecular Dynamics Simulation

A molecular dynamics simulation approach was applied for the prediction of the thermal stability of oligonucleotide duplexes. It was shown that the enthalpy of the DNA/DNA complex formation could be calculated using this approach. We have studied the influence of various simulation parameters on the secondary structure and the hybridization enthalpy value of Dickerson-Drew dodecamer. The optimal simulation parameters for the most reliable prediction of the enthalpy values were determined. The thermodynamic parameters (enthalpy and entropy changes) of a duplex formation were obtained experimentally for 305 oligonucleotides of various lengths and GC-content. The resulting database was studied with molecular dynamics (MD) simulation using the optimized simulation parameters. Gibbs free energy changes and the melting temperatures were evaluated using the experimental correlation between enthalpy and entropy changes of the duplex formation and the enthalpy values calculated by the MD simulation. The average errors in the predictions of enthalpy, the Gibbs free energy change, and the melting temperature of oligonucleotide complexes were 11%, 10%, and 4.4 °C, respectively. We have shown that the molecular dynamics simulation gives a possibility to calculate the thermal stability of native DNA/DNA complexes a priori with an unexpectedly high accuracy.

Probing the relationship between anti-Pneumocystis carinii activity and DNA binding of bisamidines by molecular dynamics simulations

The anti-Pneumocystis carinii activity of 13 synthetic pentamidine analogs was analyzed. The experimental differences in melting points of DNA dodecamer 5'-(CGCGAATTCGCG)₂-3' complexes (ΔT_m), and in the biological activity measured using ATP bioluminescence assay (IC₅₀) together with the theoretical free energy of DNA-ligand binding estimated by the proposed computational protocol, showed that the experimental activity of the tested pentamidines appeared to be due to the binding to the DNA minor groove with extended AT sequences. The effect of heteroatoms in the aliphatic linker, and the sulfonamide or methoxy substituents on the compound inducing changes in the interactions with the DNA minor groove was examined and was correlated with biological activity. In computational analysis, the explicit solvent approximation with the discrete water molecules was taken into account, and the role of water molecules in the DNA-ligand complexes was defined.

Role of DNA conformation & energetic insights in Msx-1-DNA recognition as revealed by molecular dynamics studies on specific and nonspecific complexes

In most of homeodomain-DNA complexes, glutamine or lysine is present at 50th position and interacts with 5th and 6th nucleotide of core recognition region. Molecular dynamics simulations of Msx-1-DNA complex (Q50-TG) and its variant complexes, that is specific (Q50K-CC), nonspecific (Q50-CC) having mutation in DNA and (Q50K-TG) in protein, have been carried out. Analysis of protein-DNA interactions and structure of DNA in specific and nonspecific complexes show that amino acid residues use sequence-dependent shape of DNA to interact. The binding free energies of all four complexes were analysed to define role of amino acid residue at 50th position in terms of binding strength considering the variation in DNA on stability of protein-DNA complexes. The order of stability of protein-DNA complexes shows that specific complexes are more stable than nonspecific ones. Decomposition analysis shows that N-terminal amino acid residues have been found to contribute maximally in binding free energy of protein-DNA complexes. Among specific protein-DNA complexes, K50 contributes more as compared to Q50 towards binding free energy in respective complexes. The sequence dependence of local conformation of DNA enables Q50/Q50K to make hydrogen bond with nucleotide(s) of DNA. The changes in amino acid sequence of protein are accommodated and stabilized around TAAT core region of DNA having variation in nucleotides.

Recognition of the DNA minor groove by thiazotropsin analogues

Solution-phase self-association characteristics and DNA molecular-recognition properties are reported for three close analogues of minor-groove-binding ligands from the thiazotropsin class of lexitropsin molecules; they incorporate isopropyl thiazole as a lipophilic building block. Thiazotropsin B (AcImPy(iPr) ThDp) shows similar self-assembly characteristics to thiazotropsin A (FoPyPy(iPr) ThDp), although it is engineered, by incorporation of imidazole in place of N-methyl pyrrole, to swap its DNA recognition target from 5'-ACTAGT-3' to 5'-ACGCGT-3'. Replacement of the formamide head group in thiazotropsin A by nicotinamide in AIK-18/51 results in a measureable difference in solution-phase self-assembly character and substantially enhanced DNA association characteristics. The structures and associated thermodynamic parameters of self-assembled ligand aggregates and their complexes with their respective DNA targets are considered in the context of cluster targeting of DNA by minor-groove complexes.

Atomistic account of structural and dynamical changes induced by small binders in the double helix of a short DNA

Nucleic acids are flexible molecules and their dynamical properties play a key role in molecular recognition events. Small binders interacting with DNA fragments induce both structural and dynamical changes in the double helix. We study the dynamics of a DNA dodecamer and of its complexes with Hoechst 33258, which is a minor groove binder, and with the ethidium cation, which is an intercalator, by molecular dynamics simulation. The thermodynamics of DNA-drug interaction is evaluated in connection with the structure and the dynamics of the resulting complexes. We identify and characterize the relevant changes in the configurational distribution of the DNA helix and relate them to the corresponding entropic contributions to the binding free energy. The binder Hoechst locks the breathing motion of the minor groove inducing a reduction of the configurational entropy of the helix, which amounts to 20 kcal mol(-1). In contrast, intercalations with the ethidium cation enhance the flexibility of the double helix. We show that the balance between the energy required to deform the helix for the intercalation and the gain in configurational entropy is the origin of cooperativity in the binding of a second ethidium and of anti-cooperativity in the binding of a third one. The results of our study provide an understanding of the relation between structure, dynamics and energetics in the interaction between DNA fragments and small binders, highlighting the role of dynamical changes and consequent variation of the configurational entropy of the DNA double helix for both types of binders.

Complexation of DNA with ruthenium organometallic compounds: the high complexation ratio limit

Organometallic compounds possess two modes of interaction with DNA: intercalation at low complexation ratios and electrostatic adsorption at high ratios.

Groove binding mechanism of ionic liquids: a key factor in long-term stability of DNA in hydrated ionic liquids?

Nucleic acid sample storage is of paramount importance in biotechnology and forensic sciences. Very recently, hydrated ionic liquids (ILs) have been identified as ideal media for long-term DNA storage. Hence, understanding the binding characteristics and molecular mechanism of interactions of ILs with DNA is of both practical and fundamental interest. Here, we employ molecular dynamics simulations and spectroscopic experiments to unravel the key factors that stabilize DNA in hydrated ILs. Both simulation and experimental results show that DNA maintains the native B-conformation in ILs. Simulation results further suggest that, apart from the electrostatic association of IL cations with the DNA backbone, groove binding of IL cations through hydrophobic and polar interactions contributes significantly to DNA stability. Circular dichroism spectral measurements and fluorescent dye displacement assay confirm the intrusion of IL molecules into the DNA minor groove. Very interestingly, the IL ions were seen to disrupt the water cage around DNA, including the spine of hydration in the minor groove. This partial dehydration by ILs likely prevents the hydrolytic reactions that denature DNA and helps stabilize DNA for the long term. The detailed understanding of IL-DNA interactions provided here could guide the future development of novel ILs, specific for nucleic acid solutes.

Design, synthesis and antibacterial activity of minor groove binders: the role of non-cationic tail groups

The design and synthesis of a new class of minor groove binder (MGBs) in which, the cationic tail group has been replaced by a neutral, polar variant including cyanoguanidine, nitroalkene, and trifluoroacetamide groups. Antibacterial activity (against Gram positive bacteria) was found for both the nitroalkene and trifluoroacetamide groups. For the case of the nitroalkene tail group, strong binding of a minor groove binder containing this tail group was demonstrated by both DNA footprinting and melting temperature measurements, showing a correlation between DNA binding and antibacterial activity. The compounds have also been evaluated for binding to the hERG ion channel to determine whether non-cationic but polar substituents might have an advantage compared with conventional cationic tail groups in avoiding hERG binding. In this series of compounds, it was found that whilst non-cationic compounds generally had lower affinity to the hERG ion channel, all of the compounds studied bound weakly to the hERG ion channel, probably associated with the hydrophobic head groups.

Unprocessed viral DNA could be the primary target of the HIV-1 integrase inhibitor raltegravir

Integration of HIV DNA into host chromosome requires a 3'-processing (3'-P) and a strand transfer (ST) reactions catalyzed by virus integrase (IN). Raltegravir (RAL), commonly used in AIDS therapy, belongs to the family of IN ST inhibitors (INSTIs) acting on IN-viral DNA complexes (intasomes). However, studies show that RAL fails to bind IN alone, but nothing has been reported on the behaviour of RAL toward free viral DNA. Here, we assessed whether free viral DNA could be a primary target for RAL, assuming that the DNA molecule is a receptor for a huge number of pharmacological agents. Optical spectroscopy, molecular dynamics and free energy calculations, showed that RAL is a tight binder of both processed and unprocessed LTR (long terminal repeat) ends. Complex formation involved mainly van der Waals forces and was enthalpy driven. Dissociation constants (Kds) revealed that RAL affinity for unbound LTRs was stronger than for bound LTRs. Moreover, Kd value for binding of RAL to LTRs and IC50 value (half concentration for inhibition) were in same range, suggesting that RAL binding to DNA and ST inhibition are correlated events. Accommodation of RAL into terminal base-pairs of unprocessed LTR is facilitated by an extensive end fraying that lowers the RAL binding energy barrier. The RAL binding entails a weak damping of fraying and correlatively of 3'-P inhibition. Noteworthy, present calculated RAL structures bound to free viral DNA resemble those found in RAL-intasome crystals, especially concerning the contacts between the fluorobenzyl group and the conserved 5'C(4)pA(3)3' step. We propose that RAL inhibits IN, in binding first unprocessed DNA. Similarly to anticancer drug poisons acting on topoisomerases, its interaction with DNA does not alter the cut, but blocks the subsequent joining reaction. We also speculate that INSTIs having viral DNA rather IN as main target could induce less resistance.

From multiply active natural product to candidate drug? Antibacterial (and other) minor groove binders for DNA

Natural products that bind to DNA in the minor groove are valuable templates for drug design. Examples include distamycin, netropsin, duocarmycin and anthramycin. Anticancer and anti-infective drugs feature strongly amongst their derivatives. The structures and activities of chemotypes with various therapeutic actions are discussed in the context of the broader field of therapeutically active minor groove binders. The evolution of a series of exceptionally potent and nontoxic antibacterial compounds is discussed using the general design principle of introducing additional hydrophobicity into the distamycin template to increase the strength of binding to DNA. As well as potent antibacterial compounds, antifungal and antiparasitic compounds with exceptional cellular activity against trypanosomes have been identified. Possible mechanisms of action including gene regulation and topoisomerase inhibition are discussed with the need in mind to understand selective toxicity in the series to support future drug discovery.

Molecular insight into conformational transition of amyloid β-peptide 42 inhibited by (-)-epigallocatechin-3-gallate probed by molecular simulations

Considerable experimental evidence indicates that (-)-epigallocatechin-3-gallate (EGCG) inhibits the fibrillogenesis of Aβ(42) and alleviates its associated cytotoxicity. However, the molecular mechanism of the inhibition effect of EGCG on the conformational transition of Aβ(42) remains unclear due to the limitations of current experimental techniques. In this work, molecular dynamics simulations and molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) analysis were coupled to better understand the issue. It was found that the direct interactions between EGCG and the peptide are the origin of its inhibition effects. Specifically, EGCG molecules expel water from the surface of the Aβ(42), cluster with each other, and interact directly with the peptide. The results of free energy decomposition calculated by MM-PBSA indicate that the nonpolar term contributes more than 71% to the binding free energy of the EGCG-Aβ(42) complex, while polar interactions (i.e., hydrogen bonding) play a minor role. It was identified that there are 12 important residues of Aβ(42) that strongly interact with EGCG (Phe4, Arg5, Phe19, Phe20, Glu22, Lys28, Gly29, Leu34-Gly37, and Ile41), while nonpolar interactions are mainly provided by the side chains of some hydrophobic residues (Phe, Met and Ile) and the main chains of some nonhydrophobic residues (Lys28 and Gly29). On the contrary, polar interactions are mainly formed by the main chain of Aβ(42), of which the main chains of Gly29 and Gly37 contribute greatly. The work has thus elucidated the molecular mechanism of the inhibition effect of EGCG on the conformational transition of Aβ(42), and the findings are considered critical for exploring more effective agents for the inhibition of Aβ(42) fibrillogenesis.

Amide isosteres in structure-activity studies of antibacterial minor groove binders

Antibacterial minor groove binders related to the natural product, distamycin, are development candidates for novel antibiotics. Alkenes have been found to be effective substitutes for the isosteric amide links in some positions and alkyl groups larger than methyl have been found to increase binding to DNA in both selectivity and affinity. However the impact of other isosteres such as diazenes and the position of an alkyl group with respect to DNA binding and antibacterial activity are not known. The effects of some systematic variations in the structure of polyamide minor groove binders are investigated. Isosteres of the amide link (alkenes and diazenes) are compared: it is shown that all three are competent for binding to DNA but that alkene links give the tightest binding and highest antibacterial activity; no significant antibacterial activity was found for compounds with a diazene link. Within a series of alkene linked compounds, the effect of branched N-alkyl substituents on binding to DNA and antibacterial activity is investigated: it was found that C3 and C4 branched chains are acceptable at the central pyrrole residue but that at the pyrrole ring adjacent to the basic tail group, a C4 branched chain was too large both for DNA binding and for antibacterial activity. The active branched alkyl chain compounds were found to be especially active against Mycobacterium aurum, a bacterium related to the causative agent of tuberculosis.

Polyamide-scorpion cyclam lexitropsins selectively bind AT-rich DNA independently of the nature of the coordinated metal

Cyclam was attached to 1-, 2- and 3-pyrrole lexitropsins for the first time through a synthetically facile copper-catalyzed "click" reaction. The corresponding copper and zinc complexes were synthesized and characterized. The ligand and its complexes bound AT-rich DNA selectively over GC-rich DNA, and the thermodynamic profile of the binding was evaluated by isothermal titration calorimetry. The metal, encapsulated in a scorpion azamacrocyclic complex, did not affect the binding, which was dominated by the organic tail.

Ranking ligand affinity for the DNA minor groove by experiment and simulation

The structural and thermodynamic basis for the strength and selectivity of the interactions of minor groove binders (MGBs) with DNA is not fully understood. In 2003, we reported the first example of a thiazole-containing MGB that bound in a phase-shifted pattern that spanned six base pairs rather than the usual four (for tricyclic distamycin-like compounds). Since then, using DNA footprinting, NMR spectroscopy, isothermal titration calorimetry, and molecular dynamics, we have established that the flanking bases around the central four being read by the ligand have subtle effects on recognition. We have investigated the effect of these flanking sequences on binding and the reasons for the differences and established a computational method to rank ligand affinity against varying DNA sequences.

Revisiting the association of cationic groove-binding drugs to DNA using a Poisson-Boltzmann approach

Proper modeling of nonspecific salt-mediated electrostatic interactions is essential to understanding the binding of charged ligands to nucleic acids. Because the linear Poisson-Boltzmann equation (PBE) and the more approximate generalized Born approach are applied routinely to nucleic acids and their interactions with charged ligands, the reliability of these methods is examined vis-à-vis an efficient nonlinear PBE method. For moderate salt concentrations, the negative derivative, SK(pred), of the electrostatic binding free energy, DeltaG(el), with respect to the logarithm of the 1:1 salt concentration, [M(+)], for 33 cationic minor groove drugs binding to AT-rich DNA sequences is shown to be consistently negative and virtually constant over the salt range considered (0.1-0.4 M NaCl). The magnitude of SK(pred) is approximately equal to the charge on the drug, as predicted by counterion condensation theory (CCT) and observed in thermodynamic binding studies. The linear PBE is shown to overestimate the magnitude of SK(pred), whereas the nonlinear PBE closely matches the experimental results. The PBE predictions of SK(pred) were not correlated with DeltaG(el) in the presence of a dielectric discontinuity, as would be expected from the CCT. Because this correlation does not hold, parameterizing the PBE predictions of DeltaG(el) against the reported experimental data is not possible. Moreover, the common practice of extracting the electrostatic and nonelectrostatic contributions to the binding of charged ligands to biopolyelectrolytes based on the simple relation between experimental SK values and the electrostatic binding free energy that is based on CCT is called into question by the results presented here. Although the rigid-docking nonlinear PB calculations provide reliable predictions of SK(pred), at least for the charged ligand-nucleic acid complexes studied here, accurate estimates of DeltaG(el) will require further development in theoretical and experimental approaches.