Multistep Drug Intercalation: Molecular Dynamics and Free Energy Studies of the Binding of Daunomycin to DNA

Mechanism of Daunomycin Intercalation into DNA from Enhanced Sampling Simulations

Daunomycin is a widely used anticancer drug, yet the mechanism underlying how it binds to DNA remains contested. 469 all-atom trajectories of daunomycin binding to the DNA oligonucleotide d(GCG CAC GTG CGC) were collected using weighted ensemble (WE)-enhanced sampling. Mechanistic insights were revealed through analysis of the ensemble of trajectories. Initially, the binding process involves a ubiquitous hydrogen bond between the DNA backbone and the NH3+ group on daunomycin. During the binding process, most trajectories exhibited similar structural changes to DNA, including DNA base pair rise, bending, and minor groove width changes. Variability within the ensemble of binding trajectories illuminates differences in the orientation of daunomycin as it initially intercalates; around 10% of trajectories needed minimal rearrangement from intercalation to reaching the fully bound configuration, whereas most needed an additional 1-5 ns to rearrange. The results here emphasize the utility of generating an ensemble of trajectories to discern biomolecular binding mechanisms.

Molecular Mechanism of Ligand Binding to the Minor Groove of DNA

Although DNA-ligand binding is pervasive in biology, little is known about molecular-level binding mechanisms. Using all-atom, explicit-solvent molecular dynamics simulations in conjunction with weighted ensemble (WE)-enhanced sampling, an ensemble of 2562 binding trajectories of Hoechst 33258 (H33258) to d(CGC AAA TTT GCG) was generated from which the binding mechanism was extracted. In particular, the electrostatic interaction between the positively charged H33258 and the negatively charged DNA backbone drives the formation of initial H33258-DNA contacts. After this initial contact, a hinge-like intermediate state is formed in which one end of H33258 inserts into the minor groove of DNA. Following hinge state formation is a concerted motion whereby the second end of H33258 swings into the minor groove and the spine of hydration along the minor groove causing dehydration. This study illustrates how WE-enhanced simulations of biomolecular ligation processes can offer novel mechanistic insights by generating ensembles of binding events.

High-affinity binding at quadruplex-duplex junctions: rather the rule than the exception

Quadruplex-duplex (Q-D) junctions constitute unique structural motifs in genomic sequences. Through comprehensive calorimetric as well as high-resolution NMR structural studies, Q-D junctions with a hairpin-type snapback loop coaxially stacked onto an outer G-tetrad were identified to be most effective binding sites for various polycyclic quadruplex ligands. The Q-D interface is readily recognized by intercalation of the ligand aromatic core structure between G-tetrad and the neighboring base pair. Based on the thermodynamic and structural data, guidelines for the design of ligands with enhanced selectivity towards a Q-D interface emerge. Whereas intercalation at Q-D junctions mostly outcompete stacking at the quadruplex free outer tetrad or intercalation between duplex base pairs to varying degrees, ligand side chains considerably contribute to the selectivity for a Q-D target over other binding sites. In contrast to common perceptions, an appended side chain that additionally interacts within the duplex minor groove may confer only poor selectivity. Rather, the Q-D selectivity is suggested to benefit from an extension of the side chain towards the exposed part of the G-tetrad at the junction. The presented results will support the design of selective high-affinity binding ligands for targeting Q-D interfaces in medicinal but also technological applications.

Graphene Quantum Dot Assisted Translocation of Daunomycin through an Ordered Lipid Membrane: A Study by Fluorescence Lifetime Imaging Microscopy and Resonance Energy Transfer

Daunomycin (DN) is a well-known chemotherapy drug frequently used in treating acute myeloid and lymphoblastic leukemia. It needs to be delivered to the therapeutic target by a delivering agent that beats the blood-brain barrier. DN is known to be specifically located at the membrane surface and scantly to the bilayer. Penetration of DN into the membrane bilayer depends on the molecular packing of the lipid. It does not travel promptly to the interior of the cells and needs a carrier to serve the purpose. Here, we have demonstrated, by fluorescence lifetime imaging spectroscopy (FLIM) and resonance energy transfer (RET) phenomenon, that ultrasmall graphene quantum dots (GQDs) can be internalized into the aqueous pool of giant unilamellar vesicles (GUVs) made from dipalmitoylphosphatidylcholine (DPPC) lipids, which, in turn, help in fast translocation of DN through the membrane without any delivery vehicle.

Deciphering the mechanism of action of 5FDQD and the design of new neutral analogues for the FMN riboswitch: a well-tempered metadynamics simulation study

The FMN riboswitch is a novel drug target for the design of new antibiotics, and efforts have been made to design new charged and uncharged ligands. Uncharged ligands have shown advantages of not requiring any transporter for intracellular transport or proteins for their phosphorylation. 5FDQD (5-(3-(4-fluorophenyl)butyl)-7,8-dimethylpyrido(3,4-b)quinoxaline-1,3(2H,5H)-dione) is a recently reported neutral ligand for the FMN riboswitch active against Clostridium difficile infection in mice. However, the crystal structure of the 5FDQD bound FMN riboswitch is not available, and the mechanism of ligand binding and triggering the function of the riboswitch is not well understood. We have examined 5FDQD for its binding affinity with the FMN riboswitch using the well-tempered metadynamics (WT-MtD) simulation technique. The crystal structure of the FMN riboswitch shows that the FMN interacts with the J4/5 region through the phosphate group with G62; however, the uncharged ligands take advantage of π-π stacking interactions with the same residue of the riboswitch observed from the WT-MtD simulation results. The simulation results show that the presence of fluorine on the phenyl ring in 5FDQD is important to enhance the binding affinity of the neutral ligands with the FMN riboswitch. The WT-MtD results showed that the 1,2-difluoro substitution on the phenyl ring in 5FDQD (FMN-difluoro2) and the 1,3 positions in the phenyl ring (FMN-difluoro1) showed weaker binding energy with the FMN riboswitch compared to 5FDQD. The substitution of another fluorine atom at the 5-position of the phenyl ring (FMN-trifluoro) showed a comparable binding affinity (∼-31.4 kcal mol^-1) to 5FDQD. Electron-donating substitution on the phenyl ring such as the amino group also lowered the binding affinity (-28.8 kcal mol^-1) with the riboswitch compared to 5FDQD. The computed results suggest that the position and nature of substitution in the phenyl ring of the uncharged ligands affect the overall binding and such a delicate balance is important to achieve superior binding affinity with the FMN riboswitch.

Screening of the Prime bioactive compounds from Aloe vera as potential anti-proliferative agents targeting DNA

BACKGROUND

Aloe vera extract and its bioactive compounds possess anti-proliferative properties against cancer cells. However, no detailed molecular mechanism of action studies has been reported. We have now employed a computational approach to scrutinize the molecular mechanism of lead bioactive compounds from Aloe vera that potentially inhibit DNA synthesis.

METHODS

Initially, the anti-proliferative activity of Aloe vera extract was examined in human breast cancer cells (in vitro/in vivo). Later on, computational screening of bioactive compounds from Aloe vera targeting DNA was performed by molecular docking and molecular dynamics simulation.

RESULTS

In-vitro and in-vivo studies confirm that Aloe vera extract effectively suppresses the growth of breast cancer cells without significant cytotoxicity towards non-cancerous normal immortal cells. Computational screening predicts that growth suppression may be due to the presence of DNA intercalating bioactive compounds (riboflavin, daidzin, aloin, etc.) contained in Aloe vera. MM/PBSA calculation showed that riboflavin has a higher binding affinity at the DNA binding sites compared to standard drug daunorubicin.

CONCLUSIONS

These observations support the hypothesis that riboflavin may be exploited as an anti-proliferative DNA intercalating agent to prevent cancer and is worthy of testing for the management of cancer by performing more extensive pre-clinical and if validated clinical trials.

Long-range DNA-water interactions

DNA functions only in aqueous environments and adopts different conformations depending on the hydration level. The dynamics of hydration water and hydrated DNA leads to rotating and oscillating dipoles that, in turn, give rise to a strong megahertz to terahertz absorption. Investigating the impact of hydration on DNA dynamics and the spectral features of water molecules influenced by DNA, however, is extremely challenging because of the strong absorption of water in the megahertz to terahertz frequency range. In response, we have employed a high-precision megahertz to terahertz dielectric spectrometer, assisted by molecular dynamics simulations, to investigate the dynamics of water molecules within the hydration shells of DNA as well as the collective vibrational motions of hydrated DNA, which are vital to DNA conformation and functionality. Our results reveal that the dynamics of water molecules in a DNA solution is heterogeneous, exhibiting a hierarchy of four distinct relaxation times ranging from ∼8 ps to 1 ns, and the hydration structure of a DNA chain can extend to as far as ∼18 Å from its surface. The low-frequency collective vibrational modes of hydrated DNA have been identified and found to be sensitive to environmental conditions including temperature and hydration level. The results reveal critical information on hydrated DNA dynamics and DNA-water interfaces, which impact the biochemical functions and reactivity of DNA.

Concentration-dependent mode of binding of drug oxatomide with DNA: multi-spectroscopic, voltammetric and metadynamics simulation analysis

The interaction between antihistaminic drug oxatomide (OXT) and calf-thymus DNA (CT-DNA) has been investigated in a physiological buffer (pH 7.4) using UV-Vis, fluorescence, ¹H NMR and circular dichroism spectral techniques coupled with viscosity measurements, KI quenching, voltammetry and in silico molecular modeling studies. OXT binds with CT-DNA in a concentration-dependent manner. At a lower [Drug]/[CT-DNA] molar ratio (0.6-0.1), OXT intercalates into the base pairs of CT-DNA, while at a higher [Drug]/[CT-DNA] molar ratio (13-6), the drug binds in the minor grooves of CT-DNA. The binding constants for the interaction are found to be in the order of 10³-10⁵ M^-1, and the groove binding mode of interaction exhibits a slightly higher binding constant than that of intercalative mode. Thermodynamic analysis of binding constants at three different temperatures suggests that both these modes of binding are mainly driven by hydrophobic interactions (ΔH^o > 0 and ΔS^o > 0). Voltammetric investigations indicate that the electro-reduction of OXT is an adsorption controlled process and shifts in reduction peak potentials reiterate the concentration-dependent mode of binding of the drug with CT-DNA. The free energy landscape obtained at the all-atom level, using metadynamics simulation studies, revealed two major binding forces: partial intercalation and minor groove binding, which corroborate well with the experimental results.Communicated by Ramaswamy H. Sarma.

Fine-tuning the DNA conductance by intercalation of drug molecules

In this work we study the structure-transport property relationships of small ligand intercalated DNA molecules using a multiscale modeling approach where extensive ab initio calculations are performed on numerous MD-simulated configurations of dsDNA and dsDNA intercalated with two different intercalators, ethidium and daunomycin. DNA conductance is found to increase by one order of magnitude upon drug intercalation due to the local unwinding of the DNA base pairs adjacent to the intercalated sites, which leads to modifications of the density of states in the near-Fermi-energy region of the ligand-DNA complex. Our study suggests that the intercalators can be used to enhance or tune the DNA conductance, which opens new possibilities for their potential applications in nanoelectronics.

Potential Unwinding of Double-Stranded DNA upon Binding to a Carbon Nitride Polyaniline (C3N) Nanosheet

Recently, carbon nitride polyaniline (C₃N) had attracted considerable attention from many scientific fields after its successful synthesis. However, thus far, limited efforts were devoted to reveal its potential effect to biomolecules, which correlated intimately with its further utilization. In this study, by using a molecular dynamics (MD) simulation approach, we investigated in detail the interaction between C₃N and a double-stranded DNA (dsDNA) segment to expose the underlying biological effect of C₃N to dsDNA and the corresponding molecular basis. MD simulation results demonstrated that dsDNA presented serious damages upon adsorption onto a C₃N nanosheet with the terminal base pairs denaturized, unwound, and directly packing on the C₃N surface, which implied that C₃N was potentially deleterious to biomolecules. This binding/unwinding process was mainly guided by a combination of van der Waals and π-π stacking interactions together with a continuous lateral migration of dsDNA. Moreover, the nanoscale dewetting also played an important role during the adsorption. These findings revealed the potential bio-effect of the C₃N nanomaterial and its molecular mechanism, which might benefit the future applications of C₃N-based nanostructures.

Stacking Effects on Anthraquinone/DNA Charge-Transfer Electronically Excited States

The design of more efficient photosensitizers is a matter of great importance in the field of cancer treatment by means of photodynamic therapy. One of the main processes involved in the activation of apoptosis in cancer cells is the oxidative stress on DNA once a photosensitizer is excited by light. As a consequence, it is very relevant to investigate in detail the binding modes of the chromophore with DNA, and the nature of the electronically excited states that participate in the induction of DNA damage, for example, charge-transfer states. In this work, we investigate the electronic structure of the anthraquinone photosensitizer intercalated into a double-stranded poly(dG-dC) decamer model of DNA. First, the different geometric configurations are analyzed by means of classical molecular dynamics simulations. Then, the excited states for the most relevant poses of anthraquinone inside the binding pocket are computed by an electrostatic-embedding quantum mechanics/molecular mechanics approach, where anthraquinone and one of the nearby guanine residues are described quantum mechanically to take into account intermolecular charge-transfer states. The excited states are characterized as monomer, exciton, excimer, and charge-transfer states based on the analysis of the transition density matrix, and each of these contributions to the total density of states and absorption spectrum is discussed in terms of the stacking interactions. These results are relevant as they represent the footing for future studies on the reactivity of anthraquinone derivatives with DNA and give insights on possible geometrical configurations that potentially favor the oxidative stress of DNA.

Thermodynamics and Kinetics of Drug-Target Binding by Molecular Simulation

Computational studies play an increasingly important role in chemistry and biophysics, mainly thanks to improvements in hardware and algorithms. In drug discovery and development, computational studies can reduce the costs and risks of bringing a new medicine to market. Computational simulations are mainly used to optimize promising new compounds by estimating their binding affinity to proteins. This is challenging due to the complexity of the simulated system. To assess the present and future value of simulation for drug discovery, we review key applications of advanced methods for sampling complex free-energy landscapes at near nonergodicity conditions and for estimating the rate coefficients of very slow processes of pharmacological interest. We outline the statistical mechanics and computational background behind this research, including methods such as steered molecular dynamics and metadynamics. We review recent applications to pharmacology and drug discovery and discuss possible guidelines for the practitioner. Recent trends in machine learning are also briefly discussed. Thanks to the rapid development of methods for characterizing and quantifying rare events, simulation's role in drug discovery is likely to expand, making it a valuable complement to experimental and clinical approaches.

Nonlinear TeraHertz Transmission by Liquid Water at 1 THz

The solvation properties of liquid water originate from the transient network of hydrogen-bonded molecules. In order to probe the coupling between the different modes of this network, nonlinear terahertz (THz) spectroscopy techniques are required. Ideally, these techniques should use a minimal volume and capitalize on sensitive field-resolved detection. Here we performed open aperture z-scan transmission experiments on static liquid cells, and detect the THz fields with electro-optical techniques. We show that it is possible to quantify the nonlinear response of liquid water at ~1 THz even when large signals originate from the sample holder windows.

Towards Intense THz Spectroscopy on Water: Characterization of Optical Rectification by GaP, OH1, and DSTMS at OPA Wavelengths

Water is the most prominent solvent. The unique properties of water are rooted in the dynamical hydrogen-bonded network. While TeraHertz (THz) radiation can probe directly the collective molecular network, several open issues remain about the interpretation of these highly anharmonic, coupled bands. In order to address this problem, we need intense THz radiation able to drive the liquid into the nonlinear response regime. Firstly, in this study, we summarize the available brilliant THz sources and compare their emission properties. Secondly, we characterize the THz emission by Gallium Phosphide (GaP), 2-{3-(4-hydroxystyryl)-5,5-dimethylcyclohex-2-enylidene}malononitrile (OH1), and 4-N,N-dimethylamino-4'-N'-methyl-stilbazolium 2,4,6-trimethylbenzenesulfonate (DSTMS) crystals pumped by an amplified near-infrared (NIR) laser with tunable wavelength. We found that both OH1 as well as DSTMS could convert NIR laser radiation between 1200 and 2500 nm into THz radiation with high efficiency (> 2 × 10-4), resulting in THz peak fields exceeding 0.1 MV/cm for modest pump excitation (~ mJ/cm2). DSTMS emits the broadest spectrum, covering the entire bandwidth of our detector from ca. 0.5 to ~7 THz, also at a laser wavelength of 2100 nm. Future improvements will require handling the photothermal damage of these delicate organic crystals, and increasing the THz frequency.

Dynamical Recrossing in the Intercalation Process of the Anticancer Agent Proflavine into DNA

Intercalation into DNA is the interaction mode of some anthracycline antibiotics. Recently, the molecular mechanism of this process was explored using the static free energy landscape. Here we explore the dynamical effects in the intercalation of proflavine into DNA by calculating the transmission coefficient κ-providing a measure of the departure from transition state theory for the reaction rate constant-by examination of the recrossing events at the transition state. For that purpose, we first found the accurate transition state of this complex system-as judged by a committor analysis-using a set of all-atom simulations of total length 6.3 ms. In a subsequent calculation of the transmission coefficient κ in another extensive set of simulations the small value κ = 0.1 was found, indicating a significant departure from TST. Comparison of this result with Grote-Hynes and Kramers theories shows that neither theory is able to capture this complex system's recrossing events; the source of this striking failure is discussed, as are related aspects of the mechanism. This study suggests that, for biomolecular processes similar to this, dynamical effects essential for the process are complex in nature and require novel approaches for their elucidation.

Understanding enhanced mechanical stability of DNA in the presence of intercalated anticancer drug: Implications for DNA associated processes

Most of the anticancer drugs bind to double-stranded DNA (dsDNA) by intercalative-binding mode. Although experimental studies have become available recently, a molecular-level understanding of the interactions between the drug and dsDNA that lead to the stability of the intercalated drug is lacking. Of particular interest are the modifications of the mechanical properties of dsDNA observed in experiments. The latter could affect many biological functions, such as DNA transcription and replication. Here, we probe, via all-atom molecular dynamics (MD) simulations, the change in the mechanical properties of intercalated drug-DNA complexes for two intercalators, daunomycin and ethidium. We find that, upon drug intercalation, the stretch modulus of DNA increases significantly, whereas its persistence length and bending modulus decrease. Steered MD simulations reveal that it requires higher forces to stretch the intercalated dsDNA complexes than the normal dsDNA. Adopting various pulling protocols to study force-induced DNA melting, we find that the dissociation of dsDNA becomes difficult in the presence of intercalators. The results obtained here provide a plausible mechanism of function of the anticancer drugs, i.e., via altering the mechanical properties of DNA. We also discuss long-time consequences of using these drugs, which require further in vivo investigations.

Origin of Slow Solvation Dynamics in DNA: DAPI in Minor Groove of Dickerson-Drew DNA

The measurement and understanding of collective solvation dynamics in DNA have vital biological implications, as protein and ligand binding to DNA can be directly controlled by complex electrostatic interactions of anionic DNA and surrounding dipolar water, and ions. Time-resolved fluorescence Stokes shift (TRFSS) experiments revealed anomalously slow solvation dynamics in DNA much beyond 100 ps that follow either power-law or slow multiexponential decay over several nanoseconds. The origin of such dispersed dynamics remains difficult to understand. Here we compare results of TRFSS experiments to molecular dynamics (MD) simulations of well-known 4',6-diamidino-2-phenylindole (DAPI)/Dickerson-Drew DNA complex over five decades of time from 100 fs to 10 ns to understand the origin of such dispersed dynamics. We show that the solvation time-correlation function (TCF) calculated from 200 ns simulation trajectory (total 800 ns) captures most features of slow dynamics as measured in TRFSS experiments. Decomposition of TCF into individual components unravels that slow dynamics originating from dynamically coupled DNA-water motion, although contribution from coupled water-Na⁺ motion is non-negligible. The analysis of residence time of water molecules around the probe (DAPI) reveals broad distribution from ∼6 ps to ∼3.5 ns: Several (49 nos.) water molecules show residences time greater than 500 ps, of which at least 14 water molecules show residence times of more than 1 ns in the first solvation shell of DAPI. Most of these slow water molecules are found to occupy two hydration sites in the minor groove near DAPI binding site. The residence time of Na⁺, however, is found to vary within ∼17-120 ps. Remarkably, we find that freezing the DNA fluctuations in simulation eliminates slower dynamics beyond ∼100 ps, where water and Na⁺ dynamics become faster, although strong anticorrelation exists between them. These results indicate that primary origin of slow dynamics lies within the slow fluctuations of DNA parts that couple with nearby slow water and ions to control the dispersed collective solvation dynamics in DNA minor groove.

5-Nitro-Thiophene-Thiosemicarbazone Derivatives Present Antitumor Activity Mediated by Apoptosis and DNA Intercalation

BACKGROUND

Considering the need for the development of new antitumor drugs, associated with the great antitumor potential of thiophene and thiosemicarbazonic derivatives, in this work we promote molecular hybridization approach to synthesize new compounds with increased anticancer activity.

OBJECTIVE

Investigate the antitumor activity and their likely mechanisms of action of a series of N-substituted 2-(5-nitro-thiophene)-thiosemicarbazone derivatives.

METHODS

Methods were performed in vitro (cytotoxicity, cell cycle progression, morphological analysis, mitochondrial membrane potential evaluation and topoisomerase assay), spectroscopic (DNA interaction studies), and in silico studies (docking and molecular modelling).

RESULTS

Most of the compounds presented significant inhibitory activity; the NCIH-292 cell line was the most resistant, and the HL-60 cell line was the most sensitive. The most promising compound was LNN-05 with IC50 values ranging from 0.5 to 1.9 µg.mL-1. The in vitro studies revealed that LNN-05 was able to depolarize (dose-dependently) the mitochondrial membrane, induceG1 phase cell cycle arrest noticeably, promote morphological cell changes associated with apoptosis in chronic human myelocytic leukaemia (K-562) cells, and presented no topoisomerase II inhibition. Spectroscopic UV-vis and molecular fluorescence studies showed that LNN compounds interact with ctDNA forming supramolecular complexes. Intercalation between nitrogenous bases was revealed through KI quenching and competitive ethidium bromide assays. Docking and Molecular Dynamics suggested that 5-nitro-thiophene-thiosemicarbazone compounds interact against the larger DNA groove, and corroborating the spectroscopic results, may assume an intercalating interaction mode.

CONCLUSION

Our findings highlight 5-nitro-thiophene-thiosemicarbazone derivatives, especially LNN-05, as a promising new class of compounds for further studies to provide new anticancer therapies.

Surface Topography Effects of Globular Biomolecules on Hydration Water

Hydration water serves as a microscopic manifestation of structural stability and functions of biomolecules. To develop bio-nanomaterials in applications, it is important to study how the surface topography and heterogeneity of biomolecules result in their diversity of the hydration dynamics and energetics. We here performed molecular dynamics simulations combined with the steered molecular dynamics and umbrella sampling to investigate the dynamics and escape process associated with the free energy change of water molecules close to a globular biomolecule, i.e., hemoglobin (Hb) and G-quadruplex DNA (GDNA). The residence time, power of long-time tail, and dipole relaxation time were found to display drastic changes within the averaged hydration shell of 3.0-5.0 Å. Compared with bulk water, in the inner hydration shell, the water dipole moment displays a slower relaxation process and is more oriented toward GDNA than toward Hb, forming a hedgehog-like structure when it surrounds GDNA. In particular, a spine water structure is observed in the GDNA narrow groove. The water isotope effect not only prolongs the dynamic time scales of libration motion in the inner hydration shell and the dipole relaxation processes in the bulk but also strengthens the DNA spine water structure. The potential of the mean force profile reflects the integrity of the hydration shell structure and enables us to obtain detailed insights into the structures formed by water, such as the caged H-bond network and the edge bridge structures; it also reveals that local hydration shell free energy (LHSFE) depends on H-bond rupture processes and ranges from 0.2 to 4.2 kcal/mol. Our results demonstrate that the surface topography of a biomolecule influences the integrity of the hydration shell structure and LHSFE. Our studies are able to identify various further applications in the areas of microfluid devices and nano-dewetting on bioinspired surfaces.

DNA-binding mechanism of spiropyran photoswitches: the role of electrostatics

The binding mechanism of the protonated open form of three spiropyran derivatives into a 12-mer (poly-dAT)₂ has been unveiled by means of computational methods.

The binding mechanism of the protonated open form of three spiropyran derivatives into a 12-mer (poly-dAT)₂ has been unveiled by means of computational methods. It is found that the electrostatic term in the probe:DNA binding energy, modulates the binding mode, providing new guidelines for the design of spiropyran photoswitches with specific binding modes to DNA.

Molecular mechanism and binding free energy of doxorubicin intercalation in DNA

The intercalation process of binding doxorubicin (DOX) in DNA is studied by extensive MD simulations. Many molecular factors that control the binding affinity of DOX to DNA to form a stable complex are inspected and quantified by employing continuum solvation models for estimating the binding free energy. The modified MM-PB(GB)SA methodology provides a complete energetic profile of ΔGele, ΔGvDW, ΔGpolar, ΔGnon-polar, TΔStotal, ΔGdeform, ΔGcon, and ΔGion. To identify the sequence specificity of DOX, two different DNA sequences, d(CGATCG) or DNA1 and d(CGTACG) or DNA2, with one molecule (1 : 1 complex) or two molecule (2 : 1 complex) configurations of DOX were selected in this study. Our results show that the DNA deformation energy (ΔGdeform), the energy cost from translational and rotational entropic contributions (TΔStran+rot), the total electrostatic interactions (ΔGpolar-PB/GB + ΔGele) of incorporation, the intramolecular electrostatic interactions (ΔGele) and electrostatic polar solvation interactions (ΔGpolar-PB/GB) are all unfavorable to the binding of DOX to DNA. However, they are overcome by at least five favorable interactions: the van der Waals interactions (ΔGvDW), the non-polar solvation interaction (ΔGnon-polar), the vibrational entropic contribution (TΔSvib), and the standard concentration dependent free energies of DOX (ΔGcon) and the ionic solution (ΔGion). Specifically, the van der Waals interaction appears to be the major driving force to form a stable DOX-DNA complex. We also predict that DOX has stronger binding to DNA1 than DNA2. The DNA deformation penalty and entropy cost in the 2 : 1 complex are less than those in the 1 : 1 complex, thus they indicate that the 2 : 1 complex is more stable than the 1 : 1 complex. We have calculated the total binding free energy (BFE) (ΔGt-sim) using both MM-PBSA and MM-GBSA methods, which suggests a more stable DOX-DNA complex at lower ionic concentration. The calculated BFE from the modified MM-GBSA method for DOX-DNA1 and DOX-DNA2 in the 1 : 1 complex is -9.1 and -5.1 kcal mol-1 respectively. The same quantities from the modified MM-PBSA method are -12.74 and -8.35 kcal mol-1 respectively. The value of the total BFE ΔGt-sim in the 1 : 1 complex is in reasonable agreement with the experimental value of -7.7 ± 0.3 kcal mol-1.

Tuning the Chemoselectivity of Silyl Protected Rhamnals by Temperature and Brønsted Acidity: Kinetically Controlled 1,2-Addition vs Thermodynamically Controlled Ferrier Rearrangement

An acidity- and temperature-dependent chemoselective glycosylation of silyl-protected rhamnals with alcohols has been revealed. The reaction undergoes a 1,2-addition pathway with (±)-CSA as the catalyst at rt, affording kinetically controlled 2-deoxyl rhamnosides. In contrast, only thermodynamically controlled 2,3-unsaturated rhamnosides are formed via Ferrier rearrangement when elevating reaction temperature to 85 °C or using CF₃SO₃H instead. This tunable glycosylation allows facile and practical access to both 2-deoxyl and 2,3-unsaturated rhamnosides with excellent yields and high α-stereoselectivity.

Overcoming antibiotic resistance: Is siderophore Trojan horse conjugation an answer to evolving resistance in microbial pathogens?

Comparative study of siderophore biosynthesis pathway in pathogens provides potential targets for antibiotics and host drug delivery as a part of computationally feasible microbial therapy. Iron acquisition using siderophore models is an essential and well established model in all microorganisms and microbial infections a known to cause great havoc to both plant and animal. Rapid development of antibiotic resistance in bacterial as well as fungal pathogens has drawn us at a verge where one has to get rid of the traditional way of obstructing pathogen using single or multiple antibiotic/chemical inhibitors or drugs. 'Trojan horse' strategy is an answer to this imperative call where antibiotic are by far sneaked into the pathogenic cell via the siderophore receptors at cell and outer membrane. This antibiotic once gets inside, generates a 'black hole' scenario within the opportunistic pathogens via iron scarcity. For pathogens whose siderophore are not compatible to smuggle drug due to their complex conformation and stiff valence bonds, there is another approach. By means of the siderophore biosynthesis pathways, potential targets for inhibition of these siderophores in pathogenic bacteria could be achieved and thus control pathogenic virulence. Method to design artificial exogenous siderophores for pathogens that would compete and succeed the battle of intake is also covered with this review. These manipulated siderophore would enter pathogenic cell like any other siderophore but will not disperse iron due to which iron inadequacy and hence pathogens control be accomplished. The aim of this review is to offer strategies to overcome the microbial infections/pathogens using siderophore.

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.

Encapsulating Active Pharmaceutical Ingredients in Self-Assembling Adamantanes with Short DNA Zippers

Formulating pharmaceutically active ingredients for drug delivery is a challenge. There is a need for new drug delivery systems that take up therapeutic molecules and release them into biological systems. We propose a novel mode of encapsulation that involves matrices formed through co-assembly of drugs with adamantane hybrids that feature four CG dimers as sticky ends. Such adamantanes are accessible via inexpensive solution-phase syntheses, and the resulting materials show attractive properties for controlled release. This is demonstrated for two different hybrids and a series of drugs, including anticancer drugs, antibiotics, and cyclosporin. Up to 20 molar equivalents of active pharmaceutical ingredients (APIs) are encapsulated in hybrid materials. Encapsulation is demonstrated for DNA-binding and several non-DNA binding compounds. Nanoparticles were detected that range in size from 114-835 nm average diameter, and ζ potentials were found to be between -29 and +28 mV. Release of doxorubicin into serum at near-constant rates for 10 days was shown, demonstrating the potential for slow release. The encapsulation and release in self-assembling matrices of dinucleotide-bearing adamantanes appears to be broadly applicable and may thus lead to new drug delivery systems for APIs.

Interstrand cross-linking implies contrasting structural consequences for DNA: insights from molecular dynamics

Oxidatively-generated interstrand cross-links rank among the most deleterious DNA lesions. They originate from abasic sites, whose aldehyde group can form a covalent adduct after condensation with the exocyclic amino group of purines, sometimes with remarkably high yields. We use explicit solvent molecular dynamics simulations to unravel the structures and mechanical properties of two DNA sequences containing an interstrand cross-link. Our simulations palliate the absence of experimental structural and stiffness information for such DNA lesions and provide an unprecedented insight into the DNA embedding of lesions that represent a major challenge for DNA replication, transcription and gene regulation by preventing strand separation. Our results based on quantum chemical calculations also suggest that the embedding of the ICL within the duplex can tune the reaction profile, and hence can be responsible for the high difference in yields of formation.

Binding of anticancer drug daunomycin to a TGGGGT G-quadruplex DNA probed by all-atom molecular dynamics simulations: additional pure groove binding mode and implications on designing more selective G-quadruplex ligands

DNA G-quadruplex structures are emerging cancer-specific targets for chemotherapeutics. Ligands that bind to and stabilize DNA G-quadruplexes have the potential to be anti-cancer drugs. Lack of binding selectivity to DNA G-quadruplex over DNA duplex remains a major challenge when attempting to develop G-quadruplex ligands into successful anti-cancer drugs. Thorough understanding of the binding nature of existing non-selective ligands that bind to both DNA quadruplex and DNA duplex will help to address this challenge. Daunomycin and doxorubicin, two commonly used anticancer drugs, are examples of non-selective DNA ligands. In this study, we extended our early all-atom binding simulation studies between doxorubicin and a DNA duplex (d(CGATCG)₂) to probe the binding between daunomycin and a parallel DNA quadruplex (d(TGGGGT)₄) and DNA duplex. In addition to the end stacking mode, which mimics the mode in the crystal structure, a pure groove binding mode was observed in our free binding simulations. The dynamic and energetic properties of these two binding modes are thoroughly examined, and a detailed comparison is made between DNA quadruplex binding modes and DNA duplex binding modes. Implications on the design of more selective DNA quadruplex ligands are also discussed. Graphical abstract Top stacking and groov binding modes from the MD simulations.

Metadynamics supports molecular dynamics simulation-based binding affinities of eucalyptol and beta-cyclodextrin inclusion complexes

The development of various molecular dynamics methods enables the detailed investigation of association processes, like host–guest complexes, including their dynamics and, additionally, the release of the guest compound.

Immersive virtual reality in computational chemistry: Applications to the analysis of QM and MM data

The role of Virtual Reality (VR) tools in molecular sciences is analyzed in this contribution through the presentation of the Caffeine software to the quantum chemistry community. Caffeine, developed at Scuola Normale Superiore, is specifically tailored for molecular representation and data visualization with VR systems, such as VR theaters and helmets. Usefulness and advantages that can be gained by exploiting VR are here reported, considering few examples specifically selected to illustrate different level of theory and molecular representation.

Fifty years of DNA "breathing": Reflections on old and new approaches

The coding sequences for genes, and much other regulatory information involved in genome expression, are located 'inside' the DNA duplex. Thus the "macromolecular machines" that read-out this information from the base sequence of the DNA must somehow access the DNA "interior." Double-stranded (ds) DNA is a highly structured and cooperatively stabilized system at physiological temperatures, but is also only marginally stable and undergoes a cooperative "melting phase transition" at temperatures not far above physiological. Furthermore, due to its length and heterogeneous sequence, with AT-rich segments being less stable than GC-rich segments, the DNA genome 'melts' in a multistate fashion. Therefore the DNA genome must also manifest thermally driven structural ("breathing") fluctuations at physiological temperatures that should reflect the heterogeneity of the dsDNA stability near the melting temperature. Thus many of the breathing fluctuations of dsDNA are likely also to be sequence dependent, and could well contain information that should be "readable" and useable by regulatory proteins and protein complexes in site-specific binding reactions involving dsDNA "opening." Our laboratory has been involved in studying the breathing fluctuations of duplex DNA for about 50 years. In this "Reflections" article we present a relatively chronological overview of these studies, starting with the use of simple chemical probes (such as hydrogen exchange, formaldehyde, and simple DNA "melting" proteins) to examine the local stability of the dsDNA structure, and culminating in sophisticated spectroscopic approaches that can be used to monitor the breathing-dependent interactions of regulatory complexes with their duplex DNA targets in "real time."

Structure and dynamics of proflavine association around DNA

Proflavine is a small molecule that intercalates into DNA and, thereby, acts as an anticancer agent. Intercalation of proflavine is shown to be a two-step process in which the first step is believed to be the formation of a pre-intercalative outside bound state. Experimental studies so far have been unable to capture the nature of the outside bound state. However, the sub-millisecond timescale observed in fluorescence kinetic experiments is often attributed to the binding of proflavine outside of DNA. Here, we have performed molecular dynamics simulations with multiple proflavine molecules to study the structure and dynamics of the formation of the outside bound state of DNA at different ion concentrations. We observed that the timescale of the outside bound state formation is, at least, five orders of magnitude faster (in nanoseconds) than the experimentally reported timescale (sub-milliseconds) attributed to binding outside DNA. Moreover, we also observed the stacked arrangement of proflavine all around DNA, which is different from the experimentally predicted stacking arrangement perpendicular to the helical axis of DNA in the close vicinity of the phosphate groups. This study, therefore, provides insight into the molecular structure and dynamics of the pre-intercalative outside bound state and will help in understanding the overall intercalation mechanism.

Nemorubicin and doxorubicin bind the G-quadruplex sequences of the human telomeres and of the c-MYC promoter element Pu22

BACKGROUND

Intra-molecular G-quadruplex structures are present in the guanine rich regions of human telomeres and were found to be prevalent in gene promoters. More recently, the targeting of c-MYC transcriptional control has been suggested, because the over expression of the c-MYC oncogene is one of the most common aberration found in a wide range of human tumors.

METHODS

The interaction of nemorubicin and doxorubicin with DNA G-quadruplex structures has been studied by NMR, ESI-MS and molecular modelling, in order to obtain further information about the complex and the multiple mechanisms of action of these drugs.

RESULTS AND CONCLUSIONS

Nemorubicin intercalates between A3 and G4 of d(TTAGGGT)4 and form cap-complex at the G6pT7 site. The presence of the adenine in this sequence is important for the stabilization of the complex, as was shown by the interaction with d(TTGGGTT)4 and d(TTTGGGT)4, which form only a 1:1 complex. The interaction of doxorubicin with d(TTAGGGT)4 is similar, but the complex appears less stable. Nemorubicin also binds with high efficiency the c-MYC G-quadruplex sequence Pu22, to form a very well defined complex. Two nemorubicin molecules bind to the 3'-end and to the 5'-end, forming an additional plane of stacking over each external G-tetrad. The wild type c-MYCPu22 sequence forms with nemorubicin the same complex.

GENERAL SIGNIFICANCE

Nemorubicin and doxorubicin, not only intercalate into the duplex DNA, but also result in significant ligands for G-quadruplex DNA segments, stabilizing their structure; this may in part explain the multiple mechanisms of action of their antitumor activity.

Daunorubicin-Loaded DNA Origami Nanostructures Circumvent Drug-Resistance Mechanisms in a Leukemia Model

Many cancers show primary or acquired drug resistance due to the overexpression of efflux pumps. A novel mechanism to circumvent this is to integrate drugs, such as anthracycline antibiotics, with nanoparticle delivery vehicles that can bypass intrinsic tumor drug-resistance mechanisms. DNA nanoparticles serve as an efficient binding platform for intercalating drugs (e.g., anthracyclines doxorubicin and daunorubicin, which are widely used to treat acute leukemias) and enable precise structure design and chemical modifications, for example, for incorporating targeting capabilities. Here, DNA nanostructures are utilized to circumvent daunorubicin drug resistance at clinically relevant doses in a leukemia cell line model. The fabrication of a rod-like DNA origami drug carrier is reported that can be controllably loaded with daunorubicin. It is further directly verified that nanostructure-mediated daunorubicin delivery leads to increased drug entry and retention in cells relative to free daunorubicin at equal concentrations, which yields significantly enhanced drug efficacy. Our results indicate that DNA origami nanostructures can circumvent efflux-pump-mediated drug resistance in leukemia cells at clinically relevant drug concentrations and provide a robust DNA nanostructure design that could be implemented in a wide range of cellular applications due to its remarkably fast self-assembly (≈5 min) and excellent stability in cell culture conditions.

Searching target sites on DNA by proteins: Role of DNA dynamics under confinement

DNA-binding proteins (DBPs) rapidly search and specifically bind to their target sites on genomic DNA in order to trigger many cellular regulatory processes. It has been suggested that the facilitation of search dynamics is achieved by combining 3D diffusion with one-dimensional sliding and hopping dynamics of interacting proteins. Although, recent studies have advanced the knowledge of molecular determinants that affect one-dimensional search efficiency, the role of DNA molecule is poorly understood. In this study, by using coarse-grained simulations, we propose that dynamics of DNA molecule and its degree of confinement due to cellular crowding concertedly regulate its groove geometry and modulate the inter-communication with DBPs. Under weak confinement, DNA dynamics promotes many short, rotation-decoupled sliding events interspersed by hopping dynamics. While this results in faster 1D diffusion, associated probability of missing targets by jumping over them increases. In contrast, strong confinement favours rotation-coupled sliding to locate targets but lacks structural flexibility to achieve desired specificity. By testing under physiological crowding, our study provides a plausible mechanism on how DNA molecule may help in maintaining an optimal balance between fast hopping and rotation-coupled sliding dynamics, to locate target sites rapidly and form specific complexes precisely.

Integrating multi-scale data on homologous recombination into a new recognition mechanism based on simulations of the RecA-ssDNA/dsDNA structure

RecA protein is the prototypical recombinase. Members of the recombinase family can accurately repair double strand breaks in DNA. They also provide crucial links between pairs of sister chromatids in eukaryotic meiosis. A very broad outline of how these proteins align homologous sequences and promote DNA strand exchange has long been known, as are the crystal structures of the RecA-DNA pre- and postsynaptic complexes; however, little is known about the homology searching conformations and the details of how DNA in bacterial genomes is rapidly searched until homologous alignment is achieved. By integrating a physical model of recognition to new modeling work based on docking exploration and molecular dynamics simulation, we present a detailed structure/function model of homology recognition that reconciles extremely quick searching with the efficient and stringent formation of stable strand exchange products and which is consistent with a vast body of previously unexplained experimental results.

Molecular Origin of DNA Kinking by Transcription Factors

Binding of transcription factor (TF) proteins with DNA may cause severe kinks in the latter. Here, we investigate the molecular origin of the DNA kinks observed in the TF-DNA complexes using small molecule intercalation pathway, crystallographic analysis, and free energy calculations involving four different transcription factor (TF) protein-DNA complexes. We find that although protein binding may bend the DNA, bending alone is not sufficient to kink the DNA. We show that partial, not complete, intercalation is required to form the kink at a particular place in the DNA. It turns out that while amino acid alone can induce the desired kink through partial intercalation, protein provides thermodynamic stabilization of the kinked state in TF-DNA complexes.

Selectivity of pyoverdine recognition by the FpvA receptor of Pseudomonas aeruginosa from molecular dynamics simulations

The Gram-negative bacterium Pseudomonas aeruginosa, a ubiquitous human opportunistic pathogen, has developed resistances to multiple antibiotics. It uses its primary native siderophore, pyoverdine, to scavenge the iron essential to its growth in the outside medium and transport it back into its cytoplasm. The FpvA receptor on the bacterial outer membrane recognizes and internalizes pyoverdine bearing its iron payload, but can also bind pyoverdines from other Pseudomonads or synthetic analogues. Pyoverdine derivatives could therefore be used as vectors to deliver antibiotics into the bacterium. In this study, we use molecular dynamics and free energy calculations to characterize the mechanisms and thermodynamics of the recognition of the native pyoverdines of P. aeruginosa and P. fluorescens by FpvA. Based on these results, we delineate the features that pyoverdines with high affinity for FpvA should possess. In particular, we show that (i) the dynamics and interaction of the unbound pyoverdines with water should be optimized with equal care as the interface contacts in the complex with FpvA; (ii) the C-terminal extremity of the pyoverdine chain, which appears to play no role in the bound complex, is involved in the intermediate stages of recognition; and (iii) the length and cyclicity of the pyoverdine chain can be used to fine-tune the kinetics of the recognition mechanism.

Understanding DNA under oxidative stress and sensitization: the role of molecular modeling

DNA is constantly exposed to damaging threats coming from oxidative stress, i.e., from the presence of free radicals and reactive oxygen species. Sensitization from exogenous and endogenous compounds that strongly enhance the frequency of light-induced lesions also plays an important role. The experimental determination of DNA lesions, though a difficult subject, is somehow well established and allows to elucidate even extremely rare DNA lesions. In parallel, molecular modeling has become fundamental to clearly understand the fine mechanisms related to DNA defects induction. Indeed, it offers an unprecedented possibility to get access to an atomistic or even electronic resolution. Ab initio molecular dynamics may also describe the time-evolution of the molecular system and its reactivity. Yet the modeling of DNA (photo-)reactions does necessitate elaborate multi-scale methodologies to tackle a damage induction reactivity that takes place in a complex environment. The double-stranded DNA environment is first characterized by a very high flexibility, but also a strongly inhomogeneous electrostatic embedding. Additionally, one aims at capturing more subtle effects, such as the sequence selectivity which is of critical important for DNA damage. The structure and dynamics of the DNA/sensitizers complexes, as well as the photo-induced electron- and energy-transfer phenomena taking place upon sensitization, should be carefully modeled. Finally the factors inducing different repair ratios for different lesions should also be rationalized. In this review we will critically analyze the different computational strategies used to model DNA lesions. A clear picture of the complex interplay between reactivity and structural factors will be sketched. The use of proper multi-scale modeling leads to the in-depth comprehension of DNA lesions mechanisms and also to the rational design of new chemo-therapeutic agents.

Intercalation processes of copper complexes in DNA

The family of anticancer complexes that include the transition metal copper known as Casiopeínas® shows promising results. Two of these complexes are currently in clinical trials. The interaction of these compounds with DNA has been observed experimentally and several hypotheses regarding the mechanism of action have been developed, and these include the generation of reactive oxygen species, phosphate hydrolysis and/or base-pair intercalation. To advance in the understanding on how these ligands interact with DNA, we present a molecular dynamics study of 21 Casiopeínas with a DNA dodecamer using 10 μs of simulation time for each compound. All the complexes were manually inserted into the minor groove as the starting point of the simulations. The binding energy of each complex and the observed representative type of interaction between the ligand and the DNA is reported. With this extended sampling time, we found that four of the compounds spontaneously flipped open a base pair and moved inside the resulting cavity and four compounds formed stacking interactions with the terminal base pairs. The complexes that formed the intercalation pocket led to more stable interactions.

Validation of a computational docking methodology to identify the non-covalent binding site of ligands to DNA

Despite the biomedical consequences of carcinogen-DNA interactions and the potential of DNA as a drug target in medicinal chemistry, only a small number of studies have validated or used docking methods for the prediction of the physical binding of small molecules to DNA. Knowledge of the DNA-physically-bound ligand geometry can lead to the elucidation of the molecular-level mechanism of drugs as well as predicting the subsequent chemical interactions that lead to DNA damage from carcinogens. We sought to validate AutoDock 4.2, a docking method that includes a physics-based free energy function and a Lamarckian Genetic Algorithm, for the prediction of ligand geometries upon physical binding to DNA. We performed simulations by systematically changing the length of the search process for a comprehensive set of 32 ligand-DNA molecular systems with different physico-chemical properties, and we used a free-energy-based convergence criterion to terminate our simulations. For 11 out of 28 molecular systems for which convergence was achieved, the lowest binding free energy geometries were within 2 Å of the experimentally determined geometry. Considering all predicted sites with free energy changes within 20% of the lowest binding free energy site, we found a site within 2 Å of the experimentally determined geometry for 24 out of the 28 systems. However, the predicted hydrogen bonding interactions were different for most molecular systems compared to the same interactions in the experimentally determined geometry. We discuss reasons for the successes and failures, implications, and the importance of ensuring an adequate search in docking calculations. Overall, we concluded that AutoDock 4.2 can be used to predict the non-covalent binding geometry of a small molecule to DNA with some limitations.