Intercalation of Daunomycin into Stacked DNA Base Pairs. DFT Study of an Anticancer Drug

Staphylococcus epidermidis biofilm in inflammatory breast cancer and its treatment strategies

Bacterial biofilms represent a significant challenge in both clinical and industrial settings because of their robust nature and resistance to antimicrobials. Biofilms are formed by microorganisms that produce an exopolysaccharide matrix, protecting function and supporting for nutrients. Among the various bacterial species capable of forming biofilms, Staphylococcus epidermidis, a commensal organism found on human skin and mucous membranes, has emerged as a prominent opportunistic pathogen, when introduced into the body via medical devices, such as catheters, prosthetic joints, and heart valves. The formation of biofilms by S. epidermidis on these surfaces facilitates colonization and provides protection against host immune responses and antibiotic therapies, leading to persistent and difficult-to-treat infections. The possible involvement of biofilms for breast oncogenesis has recently created the curiosity. This paper therefore delves into S. epidermidis biofilm involvement in breast cancer. S. epidermidis biofilms can create a sustained inflammatory environment through their metabolites and can break DNA in breast tissue, promoting cellular proliferation, angiogenesis, and genetic instability. Preventing biofilm formation primarily involves preventing bacterial proliferation using prophylactic measures and sterilization of medical devices and equipment. In cancer treatment, common modalities include chemotherapy, surgery, immunotherapy, alkylating agents, and various anticancer drugs. Understanding the relationship between anticancer drugs and bacterial biofilms is crucial, especially for those undergoing cancer treatment who may be at increased risk of bacterial infections, for improving patient outcomes. By elucidating these interactions, strategies to prevent or disrupt biofilm formation, thereby reducing the incidence of infections associated with medical devices and implants, can be identified.

UV-Resonance Raman Spectra of Systems in Complex Environments: A Multiscale Modeling Applied to Doxorubicin Intercalated into DNA

UV-Resonance Raman (RR) spectroscopy is a valuable tool to study the binding of drugs to biomolecular receptors. The extraction of information at the molecular level from experimental RR spectra is made much easier and more complete thanks to the use of computational approaches, specifically tuned to deal with the complexity of the supramolecular system. In this paper, we propose a protocol to simulate RR spectra of complex systems at different levels of sophistication, by exploiting a quantum mechanics/molecular mechanics (QM/MM) approach. The approach is challenged to investigate RR spectra of a widely used chemotherapy drug, doxorubicin (DOX) intercalated into a DNA double strand. The computed results show good agreement with experimental data, thus confirming the reliability of the computational protocol.

Absorption Properties of Large Complex Molecular Systems: The DFTB/Fluctuating Charge Approach

We report on the first formulation of a novel polarizable QM/MM approach, where the density functional tight binding (DFTB) is coupled to the fluctuating charge (FQ) force field. The resulting method (DFTB/FQ) is then extended to the linear response within the TD-DFTB framework and challenged to study absorption spectra of large condensed-phase systems.

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.

Intercalation Ability of Novel Monofunctional Platinum Anticancer Drugs: A Key Step in Their Biological Action

Phenanthriplatin (PtPPH) is a monovalent platinum(II)-based complex with a large cytotoxicity against cancer cells. Although the aqua-activated drug has been assumed to be the precursor for DNA damage, it is still under debate whether the way in which that metallodrug attacks to DNA is dominated by a direct binding to a guanine base or rather by an intercalated intermediate product. Aiming to capture the mechanism of action of PtPPH, the present contribution used theoretical tools to systematically assess the sequence of all possible mechanisms on drug activation and reactivity, for example, hydrolysis, intercalation, and covalent damage to DNA. Ab initio quantum mechanical (QM) methods, hybrid QM/QM′ schemes, and independent gradient model approaches are implemented in an unbiased protocol. The performed simulations show that the cascade of reactions is articulated in three well-defined stages: (i) an early and fast intercalation of the complex between the DNA bases, (ii) a subsequent hydrolysis reaction that leads to the aqua-activated form, and (iii) a final formation of the covalent bond between PtPPH and DNA at a guanine site. The permanent damage to DNA is consequently driven by that latter bond to DNA but with a simultaneous π–π intercalation of the phenanthridine into nucleobases. The impact of the DNA sequence and the lateral backbone was also discussed to provide a more complete picture of the forces that anchor the drug into the double helix.

Topoisomerases inhibition and DNA binding mode of daunomycin-oligoarginine conjugate

Cancer is a major health issue adsorbing the attention of a biomedical research. To fight this disease, new drugs are developed, specifically tailored to target biological pathways or peculiar components of the tumour cells. Particularly interesting is the use of intercalating agents as drugs capable to bind DNA and inhibit enzymes involved in DNA metabolism. Anthracyclines are the most commonly used anticancer drugs. In particular, daunomycin is used to cancer treatment by exploiting its ability to intercalate DNA and inhibit the activity of DNA topoisomerases implicated in the replication processes. Unfortunately, clinical application of anthracyclines is limited by their side effects. The conjugation with specific carriers could affect the selectivity and reduce side effect by improving stability and/or cellular uptake properties. We here report the biochemical characterisation of a daunomycin oligopeptide conjugate containing six residues of arginine, by the analysis of its fluorescence properties, DNA interaction and topoisomerases inhibitory effects.

A computational study of intercalation of streptozotocin (STZ) into DNA base pairs

Deoxyribonucleic acid (DNA) drug intercalation is a well-known phenomenon for the treatment of cancer. Streptozotocin (STZ) is a drug agent containing toxic properties that make it good in the pancreatic cancer. The main objective of this study is the intercalation of the anticancer drug into the stacked base pair of DNA sequence with ATGC using a density functional theory (DFT) code named as ADF-Molecule. ADF code implements DFT using the Slater-type orbitals (STO) for computational analysis of atomic and molecular structures. All the calculations were carried out with the GGA and hybrid exchange correlation functional with TZ2P basis sets. It was captivatingly studied that during the intercalation process, the bonds between the DNA base pairs broken. Moreover, during the process of intercalation, the free radicals are considered responsible for disturbance in the base configurations. It was determined that the disturbances that occurred in the base pairs lead to discontinuity in the replication of that particular sequence in the DNA strand.

Role of non-coding RNA networks in leukemia progression, metastasis and drug resistance

Early-stage detection of leukemia is a critical determinant for successful treatment of the disease and can increase the survival rate of leukemia patients. The factors limiting the current screening approaches to leukemia include low sensitivity and specificity, high costs, and a low participation rate. An approach based on novel and innovative biomarkers with high accuracy from peripheral blood offers a comfortable and appealing alternative to patients, potentially leading to a higher participation rate.

Recently, non-coding RNAs due to their involvement in vital oncogenic processes such as differentiation, proliferation, migration, angiogenesis and apoptosis have attracted much attention as potential diagnostic and prognostic biomarkers in leukemia. Emerging lines of evidence have shown that the mutational spectrum and dysregulated expression of non-coding RNA genes are closely associated with the development and progression of various cancers, including leukemia. In this review, we highlight the expression and functional roles of different types of non-coding RNAs in leukemia and discuss their potential clinical applications as diagnostic or prognostic biomarkers and therapeutic targets.

Elucidating the interaction of aminophylline with calf thymus DNA using multispectroscopic and molecular docking approach

Aminophylline (Am) is a methylxanthine compound clinically applied for chronic lung diseases like asthma, bronchitis or emphysema. Chemically, it comprises theophylline and ethylenediamine in a ratio 2:1. For the widening of the therapeutic window of any class of drug or for the designing of the newer therapeutic compound, an insight into the binding mechanics of available drugs with DNA is quite imperative. In view of that, here in this study we have investigated binding mechanics of aminophylline molecule with calf thymus DNA (Ct-DNA) using various spectroscopic techniques as well as molecular docking approach. Spectral analysis employing UV-visible and fluorescence approach confirmed the formation of aminophylline-Ct-DNA complex. The binding constant was calculated as 3.5 × 10⁴ M^-1 with 0.90 as the value of binding site suggestive of minor groove binding mode of aminophylline. The groove binding mode was further confirmed through spectrofluorimetric experiments like competitive displacement assay employing ethidium bromide, hoechst and rhodamine 6 G dyes as well as iodide quenching studies. The circular dichroic spectral evaluation and molecular docking study finally validated the minor groove binding mode of aminophylline with binding energy calculated as -4.5 Kcal/mol.

Importance of model size in quantum mechanical studies of DNA intercalation

The convergence of DFT-computed interaction energies with increasing binding site model size was assessed. The data show that while accurate intercalator interaction energies can be derived from binding site models featuring only the flanking nucleotides for uncharged intercalators that bind parallel to the DNA base pairs, errors remain significant even when including distant nucleotides for intercalators that are charged, exhibit groove-binding tails that engage in noncovalent interactions with distant nucleotides, or that bind perpendicular to the DNA base pairs. Consequently, binding site models that include at least three adjacent nucleotides are required to consistently predict converged binding energies. The computationally inexpensive HF-3c method is shown to provide reliable interaction energies and can be routinely applied to such large models.

Monocyclic Quinone Structure-Activity Patterns: Synthesis of Catalytic Inhibitors of Topoisomerase II with Potent Antiproliferative Activity

The monocyclic 1,4-benzoquinone, HU-331, the direct oxidation product of cannabidiol, inhibits the catalytic activity of topoisomerase II but without inducing DNA strand breaks or generating free radicals, and unlike many fused-ring quinones exhibits minimal cardiotoxicity. Thus, monocyclic quinones have potential as anticancer agents, and investigation of the structural origins of their biological activity is warranted. New syntheses of cannabidiol and (±)-HU-331 are here reported. Integrated synthetic protocols afforded a wide range of polysubstituted resorcinol derivatives; many of the corresponding novel 2-hydroxy-1,4-benzoquinone derivatives are potent inhibitors of the catalytic activity of topoisomerase II, some more so than HU-331, whose monoterpene unit replaced by a 3-cycloalkyl unit conferred increased antiproliferative properties in cell lines with IC₅₀ values extending below 1 mM, and greater stability in solution than HU-331. The principal pharmacophore of quinones related to HU-331 was identified. Selected monocyclic quinones show potential for the development of new anticancer agents.

Graphene quantum dots mediated electron transfer in DNA base pairs

Graphene quantum dots (GQDs) were connected to [Ru(bpy)3]2+ to sense DNA-mediated charge transfer. Interaction between abasic site double stranded DNA (Abasic-DNA) and [Ru(bpy)3-GQD]2+ was investigated by absorption spectroscopy, gel electrophoresis, circular dichroism, and melting temperature measurements. The results indicate that [Ru(bpy)3-GQD]2+ could be intercalated into double stranded DNA. Using [Ru(bpy)3-GQD]2+ as a signal molecule, the charge transfer performance of DNA-intercalated [Ru(bpy)3-GQD]2+ was determined using electrochemical and electrochemiluminescence measurements. Various DNA types were immobilized on Au electrodes via Au-S bonds. Electrochemiluminescence and electrochemical measurements indicate that [Ru(bpy)3-GQD]2+ could enhance DNA-mediated charge transfer when intercalated into an abasic site of double stranded DNA. And comparing with [Ru(bpy)3]2+, it can be concluded that GQDs intercalate into the DNA duplex by acting as a base analog, thus enhancing DNA charge transfer. These findings suggest that the DNA-GQD structure could aid the development of molecular devices and electric drivers, and broaden the application of DNA charge transfer.

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.

DNA physical interaction mediated b-lymphoma treatment offered by tetra benzimidazole-substituted zinc (ii) phthalocyanine derivative

Role of heterocyclic compounds with nitrogen substitution in therapeutic frontiers is well established. The efforts made in this study are directed to dissect the biological significance of benzimidazole-substituted zinc phthalocyanine derivative. Its capacity to act as an anticancer agent against the 2 B-lymphoma cell lines (low-grade and high-grade malignancy) was found out by recording florescence using Alamar blue dye. Further cytotoxic effect at the DNA level was analyzed by performing agarose gel electrophoresis. Molecular docking studies made mechanistic details crystal clear by showing potential dual binding modes employed for interaction with DNA that include minor groove binding and intercalation between bases. This advocates this derivative as potential anticancer agent and deserves further rounds of mechanistic study for its final journey to serve as a marketed drug.

microRNAs and Acute Myeloid Leukemia Chemoresistance: A Mechanistic Overview

Up until the early 2000s, a functional role for microRNAs (miRNAs) was yet to be elucidated. With the advent of increasingly high-throughput and precise RNA-sequencing techniques within the last two decades, it has become well established that miRNAs can regulate almost all cellular processes through their ability to post-transcriptionally regulate a majority of protein-coding genes and countless other non-coding genes. In cancer, miRNAs have been demonstrated to play critical roles by modifying or controlling all major hallmarks including cell division, self-renewal, invasion, and DNA damage among others. Before the introduction of anthracyclines and cytarabine in the 1960s, acute myeloid leukemia (AML) was considered a fatal disease. In decades since, prognosis has improved substantially; however, long-term survival with AML remains poor. Resistance to chemotherapy, whether it is present at diagnosis or induced during treatment is a major therapeutic challenge in the treatment of this disease. Certain mechanisms such as DNA damage response and drug targeting, cell cycling, cell death, and drug trafficking pathways have been shown to be further dysregulated in treatment resistant cancers. miRNAs playing key roles in the emergence of these drug resistance phenotypes have recently emerged and replacement or inhibition of these miRNAs may be a viable treatment option. Herein, we describe the roles miRNAs can play in drug resistant AML and we describe miRNA-transcript interactions found within other cancer states which may be present within drug resistant AML. We describe the mechanisms of action of these miRNAs and how they can contribute to a poor overall survival and outcome as well. With the precision of miRNA mimic- or antagomir-based therapies, miRNAs provide an avenue for exquisite targeting in the therapy of drug resistant cancers.

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.

DrrC protein of Streptomyces peucetius removes daunorubicin from intercalated dnrI promoter

DrrC is a DNA-binding protein of Streptomyces peucetius that provides self-resistance against daunorubicin, the antibiotic produced by the organism. DrrC was expressed in E.coli and purified by using N-terminal MBP-tag which retained DNA-binding property in spite of the tag. Mobility shift assay confirmed the interaction of 313bp DNA that has the dnrI promoter, daunorubicin and MBP-DrrC in the presence of ATP. Biotinylated and immobilized 313bp DNA was intercalated with daunorubicin to observe the release of the drug when MBP-DrrC is allowed to act on the DNA. The release of daunorubicin was recorded by absorption and fluorescence spectroscopy. The experiments proved that daunorubicin was released from DNA in the presence of MBP-DrrC. Fluorescence emission of daunorubicin had a maximum peak at 591nm. However, emission spectrum of released daunorubicin showed hypochromism with a maximum peak at 584nm that is possibly because it is in complex with MBP-DrrC. We propose that DrrC naturally binds at intercalated sites to eject daunorubicin; in the process both drug and protein are dislodged from DNA. Like UvrA, DrrC possibly scans the DNA for intercalated daunorubicin. When it encounters daunorubicin, DrrC dislodges it, thereby allowing DNA replication and transcription to go on unhindered. Thus a novel self resistance mechanism by DNA repair is mediated by DrrC.

Direct Determination of Metal Complexes' Interaction with DNA by Atomic Telemetry and Multiscale Molecular Dynamics

The lack of molecular mechanistic understanding of the interaction between metal complexes and biomolecules hampers their potential medical use. Herein we present a robust procedure combining resonant X-ray emission spectroscopy and multiscale molecular dynamics simulations, which allows for straightforward elucidation of the precise interaction mechanism at the atomic level. The report unveils an unforeseen hydrolysis process and DNA binding of [Pt{N(p-HC₆F₄)CH₂}₂py₂] (Pt103), which showed potential cytotoxic activity in the past. Pt103 preferentially coordinates to adjacent adenine sites, instead of guanine sites as in cisplatin, because of its hydrogen bond ability. Comparison with previous research on cisplatin suggests that selective binding to guanine or adenine may be achieved by controlling the acidity of the compound.

Effects of oxygenation on the intercalation of 1,10-phenanthroline-5,6/4,7-dione between DNA base pairs: a computational study

The effects of oxygen in positions 4,7 and 5,6 of phenanthroline have been studied computationally when this ligand intercalates between DNA base pairs. Our results indicate that solvation energy could be the driving force of the process and thus, it can be also related with the cytotoxicity of the drug.

Saturation Binding of Nicotine to Synthetic Neuromelanin Demonstrated by Fluorescence Spectroscopy

Neuromelanin (NM) has long been considered as an aging pigment, perhaps an unavoidable and undesirable byproduct of dopaminergic neural transmission. However, NM is carefully packaged into double membrane-bound structures within cells of the substantia nigra and other neural tissues, suggesting a beneficial function to maintaining these stores. It is well established that NM is able to concentrate toxic xenobiotics within pigmented cells due to its unique chemical environment. In doing so, such agents may confer susceptibility to Parkinson's disease (PD) as illustrated by model PD-inducing neurotoxins such as methyl-phenyl-pyridinium ion. It is possible that high-affinity binding interactions toward NM may contribute to the adverse effects of PD-inducing toxins, as well as neuroprotective agents. Here we aim to develop a generalized assay capable of elucidating the binding constants of chemical agents to synthetic and natural neuromelanins. Toward this end, a model neuromelanin synthesized from dopamine and cysteine was prepared according to published procedure. Using a UV/Visible spectroscopic assay, we show that dopamine, 6-hydroxy dopamine, and nicotine bind to the synthetic neuromelanin, while caffeine did not. More importantly, nicotine was further found to induce a fluorescence signal in the presence of NM which was used to establish a binding constant estimated at 0.65 mM. Dopamine appears to enhance this signal, also in a saturable manner, with an estimated Kd of 0.05 mM in our isolated chemical system. In summary, the micro-scale fluorescence assay described herein will allow us to overcome many of the problems inherent in the study of chemical interaction with NM through traditional spectroscopic means. Using a single standardized signal, it should now be possible to rank a number of PD-related toxins based on NM-binding affinity and shed further light on this important problem.

Modeling techniques and fluorescence imaging investigation of the interactions of an anthraquinone derivative with HSA and ctDNA

A new anthraquinone derivative (AORha) was synthesized. Its interactions with human serum albumin (HSA) and calf thymus DNA (ctDNA) were investigated by fluorescence spectroscopy, UV-visible absorption spectroscopy and molecular modeling. Cell viability assay and cell imaging experiment were performed using cervical cancer cells (HepG2 cells). The fluorescence results revealed that the quenching mechanism was static quenching. At different temperatures (290, 300, 310 K), the binding constants (K) and the number of binding sites (n) were determined, respectively. The positive ΔH and ΔS values showed that the binding of AORha with HSA was hydrophobic force, which was identical with the molecular docking result. Studying the fluorescence spectra, UV spectra and molecular modeling also verified that the binding mode of AORha and ctDNA might be intercalative. When HepG2 cells were treated with AORha, the fluorescence became brighter and turned green, which could be used for bioimaging.

Interaction of doxorubicin with polynucleotides. A spectroscopic study

The interaction of doxorubicin (DX) with model polynucleotides poly(dG-dC)·poly(dG-dC) (polyGC), poly(dA-dT)·poly(dA-dT) (polyAT), and calf thymus DNA has been studied by several spectroscopic techniques in phosphate buffer aqueous solutions. UV-vis, circular dichroism, and fluorescence spectroscopic data confirm that intercalation is the prevailing mode of interaction, and also reveal that the interaction with AT-rich regions leads to the transfer of excitation energy to DX not previously documented in the literature. Moreover, the DX affinity for AT sites has been found to be on the same order of magnitude as that reported for GC sites.

B-DNA Structure and Stability as Function of Nucleic Acid Composition: Dispersion-Corrected DFT Study of Dinucleoside Monophosphate Single and Double Strands

We have computationally investigated the structure and stability of all 16 combinations of two out of the four natural DNA bases A, T, G and C in a di-2′-deoxyribonucleoside-monophosphate model DNA strand as well as in 10 double-strand model complexes thereof, using dispersion-corrected density functional theory (DFT-D). Optimized geometries with B-DNA conformation were obtained through the inclusion of implicit water solvent and, in the DNA models, of sodium counterions, to neutralize the negative charge of the phosphate groups. The results obtained allowed us to compare the relative stability of isomeric single and double strands. Moreover, the energy of the Watson–Crick pairing of complementary single strands to form double-helical structures was calculated. The latter furnished the following increasing stability trend of the double-helix formation energy: d(TpA)₂ <d(CpA)₂ <d(ApT)₂ <d(ApA)₂ <d(GpT)₂ <d(GpA)₂ <d(ApG)₂ <d(CpG)₂ <d(GpG)₂ <d(GpC)₂, where the energy differences between the last four dimers, d(ApG)₂, d(CpG)₂, d(GpG)₂ and d(GpC)₂, is within 4.0 kcal mol⁻¹, and the energy between the most and the least stable isomers is 13.4 kcal mol⁻¹. This trend shows that the formation energy essentially increases with the number of hydrogen bonds per base pair, that is two between A and T and three between G and C. Superimposed on this main trend are more subtle effects that depend on the order in which bases occur within a strand from the 5’- to the 3’-end.

Early stage intercalation of doxorubicin to DNA fragments observed in molecular dynamics binding simulations

The intercalation mode between doxorubicin (an anticancer drug) and two 6-base-pair DNA model fragments (d(CGATCG)₂ and d(CGTACG)₂) has been well studied by X-ray crystallography and NMR experimental methods. Yet, the detailed intercalation pathway at molecular level remains elusive. In this study, we conducted molecular dynamics binding simulations of these two systems using AMBER DNA (parmbsc0) and drug (GAFF) force fields starting from the unbound state. We observed outside binding (minor groove binding or end-binding) in all six independent binding simulations (three for each DNA fragment), followed by the complete intercalation of a drug molecule in two simulations (one for each DNA fragment). First, our data directly supported that the minor groove binding is the dominant pre-intercalation step. Second, we observed that the opening and flipping of a local base pair (A3-T10 for d(CGATCG)₂ and C1-G12 for d(CGTACG)₂) in the two intercalation trajectories. This locally cooperative flipping-intercalation mechanism was different from the previously proposed rise-insertion mechanism by which the distance between two neighboring intact base pairs increases to create a space for the drug insertion. Third, our simulations provided the first set of data to support the applicability of the AMBER DNA and drug force fields in drug-DNA atomistic binding simulations. Implications on the kinetics pathway and drug action are also discussed.

Homology modeling of Cry10Aa toxin from B. thuringiensis israelensis and B. thuringiensis subsp. LDC-9

A three dimensional model was developed for Cry10Aa protein sequence of B. thuringiensis LDC-9 and B. thuringiensis israelensis that has not been solved empirically by X-ray crystallography or NMR. Homology modeling was employed for the structure prediction using Cry2Aa as template protein, a high-resolution X-ray crystallography structure. The model predicted for the B. thuringiensis LDC-9 Cry10Aa protein reveals a partial N-terminal domain only due to its partial sequence of 104 amino acids. B. thuringiensis israelensis Cry10Aa model contains three domains such as domain I, a bundle of eight alpha helices with the central relatively hydrophobic helix surrounded by amphipathic helices while domain II and III contain mostly beta-sheets. Significant structural differences within domain II in this model among all Cry protein structures indicates that it is involved in recognition and binding to cell surfaces. Comparison of B. thuringiensis israelensis predicted structure with available experimentally determined Cry structures reveals identical folds. The distribution of electrostatic potential on the surface of the molecules in the model is non-uniform and identifies one side of the alpha-helical domain as negatively charged indicating orientation of toxic molecules toward the cell membrane during the initial binding with a cell surface receptor. The collective knowledge of Cry toxin structures will lead to a more critical understanding of the structural basis for receptor binding and pore formation, as well as allowing the scope of diversity to be better appreciated. This model will serve as a starting point for the design of mutagenesis experiments aimed to improve the toxicity and to provide a new tool for the elucidation of the mechanism of action of these mosquitocidal proteins.

Structural model of the pre-pore ring-like structure of Panton-Valentine leukocidin: providing dimensionality to biophysical and mutational data

Panton-Valentine leukocidin (PVL) is a bipartite toxin that plays an important role in the pathogenesis of methicillin-resistant Staphylococcus aureus. Recent clinical data suggest a correlation between PVL and severe cases of S. aureus pneumonia. A clear understanding of the structure and function of PVL is critical to the development of novel, effective treatments. Here, we report an all-atom model of the macromolecular structure of Panton-Valentine leukocidin in its octameric, pre-pore conformation that confirms and extends our understanding of the toxin's mechanism of action.