Theoretical studies on interaction of anticancer drugs (dacarbazine,procarbazine and triethylenemelamine) with normal (AT and GC) and mismatch (GG,CC, AA and TT) base pairs

This study is an attempt to gain a better understanding of the physicochemical interaction between novel anticancer drugs and DNA bases. We have employed quantum chemical tools to explore the interaction of a few anticancer drugs [namely procarbazine (PR), dacarbazine (DC) and triethylenemelamine (TR)] with isolated normal (GC and AT) and mismatch (AA, CC, GG and TT) base pairs. The molecular geometries, electronic structural stability, vibrational energies, chemical reactivity and other electronic properties were studied using MP2/6-311+G^**, B3LYP/6-311+G^** and M05-2X/6-311+G^** methods. The optimised geometries of the usual and mismatch base pairs are almost planar whereas the geometries of drug-interacting complexes deviate from planarity. The presence of steric hindrance and π-bond overlaps between C–C bonds in the complexes has distorted the planarity of the four- and five-member rings in the base pairs. Among the three drugs chosen, DC and PR bond well with normal and mismatch base pairs with large interaction energy. The electron density (ED) difference maps of the most stable GG–DC, GG–PR and GG–TR drug-interacting complexes show the information about sharing of ED and gain or loss of ED within the interacting molecules. The stabilisation energy of the charge transfer interaction between the relevant donor–acceptor orbital of GG–DC and GC–DC complexes has been found to be around 16 kcal/mol and GG–PR and GC–PR complexes has been found to be around 12 kcal/mol. But, for the GG–TR and GC–TR complexes, the stabilisation energy is found to be less than 6 kcal/mol. Moreover, the topological analysis of hydrogen bond network of DC and PR drug-interacting complexes have high electron and Laplacian density with structural stability at the bond critical points (BCPs), while compared TR drug-interacting complexes by atoms in molecules and natural bond orbital analysis. Finally, we may conclude that the drugs DC and PR are highly efficient drugs to target normal and mismatch base pair for control and inhibition of DNA replication.     

Studies on tautomeric forms of Guanine-Cytosine base pairs of nucleic acids and their interactions with water molecules

The relative stabilities of Guanine-Cytosine (G-C) DNA bare base pairs, its tautomeric forms and microhydrated base pairs are theoretically investigated with a focus on the keto-enol tautomerism as well as on the cis-trans isomerism using ab initio and density functional theory methods. The stabilities of the G-C bare base pairs, its tautomeric forms and microhydrated base pairs were affected by various factors including keto-enol tautomerization, cis-trans enol isomerization, and steric hindrance between the base pair and water molecules. The Atoms in Molecules theory (AIM) is employed to investigate H-bonding patterns both in bare and microhydrated base pairs. From the above topological results, an excellent linear correlation is shown between electron density [rho(r)], and its Laplacian [V2rho(r)] at the bond critical points. NBO analysis has been carried out to study the charge transfer between proton acceptor to the antibonding orbital of the X-H bond both in bare and microhydrated base pairs.     

Intercalation of daunomycin into stacked DNA base pairs. DFT study of an anticancer drug

We have computationally studied the intercalation of the antitumor drug daunomycin into six stacks of Watson-Crick DNA base pairs (i.e., AT-AT, AT-TA, GC-AT, CG-TA, GC-GC, GC-CG) using density functional theory (DFT). The proton affinity of the DNA intercalator daunomycin in water was computed to be 159.2 kcal/mol at BP86/TZ2P, which is in line with the experimental observation that daunomycin is protonated under physiological conditions. The intercalation interaction of protonated daunomycin with two stacked DNA base pairs was studied through a hybrid approach in which intercalation is treated at LDA/TZP while the molecular structure of daunomycin and hydrogen-bonded Watson-Crick pairs is computed at BP86/TZ2P. We find that the affinity of the drug for the six considered base pair dimers decreases in the order AT-AT > AT-TA > GC-AT > GC-TA > GC-CG > GC-GC, in excellent agreement with experimental data on the thermodynamics of the interaction between daunomycin and synthetic polynucleotides in aqueous solution. Our analyses show that the overall stability of the intercalation complexes comes mainly from pi-pi stacking but an important contribution to the computed and experimentally observed sequence specificity comes from hydrogen bonding between daunomycin and hetero atoms in the minor groove of AT base pairs.     

Theoretical study of the interaction between 5-methylcytosine and acrylamide

The hydrogen-bonded complexes between 5-methylcytosine and acrylamide have been investigated using the density function theory (DFT) method. Five stable complexes have been found with no imaginary frequencies. Complex C3 is the most stable one with interaction energies of -69.01?kJ?mol(-1) corrected for basis set superposition error (BSSE). The charge change in the process of these complexes formation has also been examined. The atoms in molecules (AIM) theory and natural bond orbital (NBO) method have been performed to investigate the hydrogen bonds involved in all the complexes. The electron density and its corresponding Laplacian at the bond and ring critical points have been analyzed. In C3 complex, there is the largest stabilization energy (18.17?kJ?mol(-1)) between N11-H12 antibonding orbital and lone electron pair of O17. It can be seen that the hydrogen bonds play a crucial role in the stability of all the complexes between 5-methylcytosine and acrylamide. The theoretical results could provide helpful information for other researchers in further work.     

Profiling cellular protein complexes by proximity ligation with dual tag microarray readout

Patterns of protein interactions provide important insights in basic biology, and their analysis plays an increasing role in drug development and diagnostics of disease. We have established a scalable technique to compare two biological samples for the levels of all pairwise interactions among a set of targeted protein molecules. The technique is a combination of the proximity ligation assay with readout via dual tag microarrays. In the proximity ligation assay protein identities are encoded as DNA sequences by attaching DNA oligonucleotides to antibodies directed against the proteins of interest. Upon binding by pairs of antibodies to proteins present in the same molecular complexes, ligation reactions give rise to reporter DNA molecules that contain the combined sequence information from the two DNA strands. The ligation reactions also serve to incorporate a sample barcode in the reporter molecules to allow for direct comparison between pairs of samples. The samples are evaluated using a dual tag microarray where information is decoded, revealing which pairs of tags that have become joined. As a proof-of-concept we demonstrate that this approach can be used to detect a set of five proteins and their pairwise interactions both in cellular lysates and in fixed tissue culture cells. This paper provides a general strategy to analyze the extent of any pairwise interactions in large sets of molecules by decoding reporter DNA strands that identify the interacting molecules.     

DFT and MD study of adsorption sensitivity of aluminium phosphide nanotube towards some air pollutant gas molecules

To investigate the adsorption behaviour of CS₂, CO₂, SO₂, H₂Se and H₂S gas molecules on the external surface of (6, 0) single-walled aluminium phosphide nanotube (AlPNT), the density functional theory (DFT) calculations at the B3LYP level of theory are performed. The partial densities of states (PDOS) for the SO₂ molecule, the S and O atoms of SO₂ molecule before and after adsorption on the surface of AlPNT have been plotted. The vibrational frequencies and physical properties such as chemical potential, chemical hardness, dipole moment and chemical electrophilicity of all studied complexes have been systematically investigated. The electron density and the Laplacian of the electron density for bond critical points have been examined by the AIM theory. Also the molecular dynamics (MD) simulations of two complexes with the minimum and maximum negative interaction energies that is: AlPNT/CO₂ and AlPNT/SO₂ complexes, respectively, have been considered.     

Experimental studies of molecular interactions between nitrogen bases of nucleic acids.

New experimental results concerning molecular interactions between the nitrogen bases of nucleic acids in the crystalline phase and in vacuo are reported. The temperature dependence of the evaporation rate is measured for solid species. The sensitivity of conventional methods of sublimation heat measurements was improved essentially using a quartz resonator serving as a precise sensor of evaporation rate. Sublimation heats were found for both canonical bases and a number of their derivatives. The in vacuo formation of base associates interacting through hydrogen bonds was observed with a field mass spectrometer. The dimer formation enthalpies, which are indicative of a stronger attraction in complementary pairs compared with noncomplementary ones, were derived from the temperature dependence of ionic currents. Hydrogen-bound complexes of more intricate associates (base trimers and aqueous molecules associates) were studied. The energy gain in the formation of trimers of identical molecules was shown to be larger (per base molecule) than that for dimers.     

AT base pairs are less stable than GC base pairs in Z-DNA: The crystal structure of d(m5CGTAm5CG)

Two hexanucleoside pentaphosphates, 5-methyl and 5-bromo cytosine derivatives of d(CpGpTp-ApCpG) have been synthesized, crystallized, and their three-dimensional structure solved. They both form left-handed Z-DNA and the methylated derivative has been refined to 1.2 Å resolution. These are the first crystal Z-DNA structures that contain AT base pairs. The overall form of the molecule is very similar to that of the unmethylated or the fully methylated (dC-dG)₃ hexamer although there are slight changes in base stacking. However, significant differences are found in the hydration of the helical groove. When GC base pairs are present, the helical groove is systematically filled with two water molecules per base pair hydrogen bonded to the bases. Both of these water molecules are not seen in the electron density map in the segments of the helix containing AT base pairs, probably because of solvent disorder. This could be one of the features that makes AT base pairs form Z-DNA less readily than GC base pairs.     

Intermolecular duplexes in heterogeneous nuclear RNA from HeLa cells

Rapidly sedimenting hnRNA complexes contain regions of stable intermolecular duplex. Disruption of such complexes, as judged by a reduction in sedimentation rate, requires conditions sufficient to denature the duplex regions. Rapidly sedimenting molecules reappear only when the complementary sequences reanneal — that is, the formation of such complexes is dependent upon time and the concentration of homologous RNA. These experiments lead us to the conclusion that rapidly sedimenting hnRNA complexes consist of two or more largely single-stranded RNA molecules held together by short duplex regions. Precisely such structures have been visualized in the electron microscope. Rapidly sedimenting fractions of native nuclear RNA from preparative sucrose gradients consist primarily of large, multi-molecular complexes interconnected by duplex regions averaging 300 base pairs in length. Exposure of the RNA to severely denaturing conditions eliminates such complexes. Reannealing of the RNA reconstitutes complexes which are indistinguishable from those observed in preparations before denaturation.     

Analysis of the conformational features of Watson–Crick duplex fragments by molecular mechanics and quantum mechanics methods

It is generally accepted that important features of the Watson–Crick duplex (WCD) originate from the molecular structure of its subunits. However, it is still unclear what properties of each subunit are responsible for significant features of the WCD structure. Computations of deoxydinucleoside monophosphate (dDMP) complexes with Na ions on the basis of the density functional theory (DFT) have shown that conformational properties of minimal single-stranded fragments of DNA play a pivotal role in the origin of the unique features of the WCD. The directionality of the sugar-phosphate backbone (SPB) and preferable ranges of its torsion angles combined with the difference in geometry between purines and pyrimidines have been found to define the nucleotide sequence dependence of the WCD three-dimensional structure. In this work, density functional theory computations were extended to minimal duplex fragments, that is, complementary dDMP (cdDMP) complexes with Na ions. Using several computational methods and various functionals, energy minima were searched for the BI conformation of cdDMP complexes with different nucleotide sequences. Two sequences were optimized using an ab initio method at the MP2/6-31++G** level of theory. An analysis of the SPB torsion angles, sugar-ring puckering, and mutual base positions in the optimized structures showed that the conformational features of cdDMP complexes with Na ions remained within the BI ranges and become more similar to the corresponding features that WCDs display in a crystal. Qualitatively, the main features of each cdDMP complex were invariant with different computational methods, although the values of certain conformational parameters could vary, but still within the limits that are typical of the corresponding family. Common functionals that are employed in DFT calculations were observed to overestimate the distance between base pairs, while MP2 computations and new complex functionals yielded structures with atom–atom contacts that are too close. Several energy minima that correspond to the BI conformation have been proven to exist for certain cdDMP complexes with Na ions, indicating that the topography of the potential energy surface is complex. This circumstance accounts for the variation of conformational parameters among duplex fragments with the same nucleotide sequence. The common AMBER and CHARMM molecular mechanics force fields reproduce many conformational characteristics of dDMPs and their complementary complexes with Na ions, but fail to reproduce certain details of the nucleotide sequence dependence of the WCD conformation.     

Collective motion in DNA and its role in drug intercalation

The effects of collective motion in DNA as reflected by resonance coupling among its intact segments have been discussed for both linear and circular DNA molecules. The results indicate that due to the effects of this kind of internal collective motion, the energy will be at times highly concentrated at some spots. As a result of the overfocus of energy, the stress built up along the direction of hydrogen bonds between complementary base pairs will be dramatically increased, rupturing a series of consecutive hydrogen bonds simultaneously and resulting in a suddenly free jerk, such that the DNA molecule will undergo a local “quake.” The “hole” formed by this kind of quake-like motion will be large enough for bulky drugs to gain entrance and intercalate into DNA. Even for smaller drugs, this local quake-like motion can also provide a significant mode of entry for intercalation. Energy minimizations carried out for DNA–drug complexes indicate that, for most drugs, a distortion or disruption of 2 to 4 base pairs occurs at the intercalation site in DNA molecules.     

Predicted mode of intercalation of doxorubicin with dinucleotide dimers

Intermolecular molecular mechanics energy calculations have been carried out for doxorubicin interacting with two dinucleotide dimer sequences. The preferred mode of intercalation is in the minor groove with the anthraquinone ring of the drug nearly perpendicular to the base pairs for the (CpG) sequence having alternate C3′ endo-C2′ endo sugar ring puckering. The preferred intercalation conformation of the drug is nearly identical to the N-bromacetyldaunomycin crystal structure. This prediction is qualitatively consistent with the recently reported crystal structure of a d(CpGpCpGpCpG) dimer-daunomycin complex. For the other dinucleotide sequence, (TpC-ApG), minor groove intercalation is also preferred, but the drug conformation can be changed.     

Molecular dynamics and quantum chemical studies of solvent effects on cyclo glycylglycine and glycylalanine dipeptides

Hydrogen bond (H-bond) interactions between the two cyclo dipeptides, cyclo(glycyl-glycine) (CGG) and cyclo(glycyl-alanine) (CGA), and water have been studied using molecular dynamics (MD) and quantum chemical methods. The MD studies have been carried out on CGG and CGA in water using fixed charge force field (AMBER ff03) for over 10 ns with a MD time step of 2 fs. The results of this study show that the solvation pattern influences the conformations of the cyclo dipeptides. Following molecular simulations, post Hartree–Fock and density functional theory methods have been used to explore the molecular properties of the cyclo dipeptides in gaseous and aqueous phase environments. The self-consistent reaction field theory has been used to optimise the cyclopeptides in diethyl ether (? = 4.3) and water (? = 78.5), and the solvent effects have been analysed. A cluster of eight water molecules leads to the formation of first solvation shell of CGG and CGA and the strong H-bonding mainly contributes to the interaction energies. The H-bond interactions have been analysed by the calculation of electron density ρ(r) and its Laplacian ▽²ρ(r) at bond critical points using atoms in molecules theory. The natural bond orbital analysis was carried out to reveal the nature of H-bond interactions. In the solvated complexes, the keto carbons registered the maximum NMR chemical shifts.     

Binding of small molecules at interface of protein–protein complex – A newer approach to rational drug design

Protein–protein interaction is a vital process which drives many important physiological processes in the cell and has also been implicated in several diseases. Though the protein–protein interaction network is quite complex but understanding its interacting partners using both in silico as well as molecular biology techniques can provide better insights for targeting such interactions. Targeting protein–protein interaction with small molecules is a challenging task because of druggability issues. Nevertheless, several studies on the kinetics as well as thermodynamic properties of protein–protein interactions have immensely contributed toward better understanding of the affinity of these complexes. But, more recent studies on hot spots and interface residues have opened up new avenues in the drug discovery process. This approach has been used in the design of hot spot based modulators targeting protein–protein interaction with the objective of normalizing such interactions.     

Distamycin-stabilized Antiparallel-Parallel Combination (APC) DNA

Abstract

The formation of Antiparallel-Parallel-Combination (APC) DNA, a liner duplex with a segment of parallel-stranded (ps) helix flanked by conventional B-DNA, was tested with a number of synthetic oligonucleotides. The groove-binding ligand distamycin A (DstA) was used to stabilize the ps segment comprising five A·T base pairs. Two drug molecules bound per APC, one in each of the two equivalent grooves characteristic of ps-DNA. APC-DNA, reference molecules and their complexes with DstA were analysed by several methods: circular dichroism and absorption spectroscopy, thermal denaturation, chemical modification, and molecular modeling. The dye binding stoichiometry differed significantly due to inherent structural differences in the groove geometries of ps-DNA (trans base pairs, similar grooves) and conventional antiparallel-stranded (aps) B-DNA (cis base pairs, distinct major and minor grooves). The data support the existence of APC folding in solution.     

Investigation of the dark interaction between furocoumarins and DNA

The complexes between some furocoumarins and DNA have been studied using various physicochemical techniques. Flow-dichroism measurements data strongly support the intercalation of the planar furocoumarin molecules between two base pairs of duplex DNA. The equilibrium dialysis and spectrophotometric data show relatively low values of the association constants of the complexes and a small number of molecules able to intercalate in DNA, thus indicating that furocoumarins have a relatively low affinity for DNA in the complex formation. The biological and photobiological consequences connected with these results are discussed.The binding curves obtained using some polynucleotides and various DNA samples having different composition with regard to base pairs, have shown that the regions of the macromolecule having alternate sequences of purine and pyrimidine represent sites useful for intercalation. No preference has been observed for A-T or G-C.     

Chromatin-like structures in polyoma virus and simian virus 10 lytic cycle.

Nucleoprotein complexes containing viral DNA and cellular histones were extracted from nuclei of permissive cells infected with polyoma virus or simian virus 40 (SV40) and examined by electron microscopy. Polyoma and SV40 nucleoprotein complexes are almost identical. They appear as relaxed circular molecules consisting of 20 to 21 globular particles interconnected by thin filaments. Their contour length in 0.02 M salt is 2.7 times shorter than that of viral DNA form I obtained after dissociation of the proteins in 1 M NaCl. The nucleosomes have an average diameter of 12.5 nm. Each nucleosome contains 175 to 205 DNA base pairs condensed fivefold in length. The nucleosomes are regularly spaced on the circular molecule. The internucleosomal filaments are made of naked DNA, and each filament contains about 55 base pairs. The partial sensitivity of the nucleoprotein complex to cleavage by EcoR1 endonuclease suggests that the nucleosomes are not formed at specific sites on the viral genome. Faster sedimenting nucleoprotein complexes containing replicative intermediates were studied. Isopycnic centrifugation in metrizamide gradients in the absence of aldehyde fixation showed that these molecules conserved the same DNA-to-protein ratio as the form I DNA-containing complexes.     

Regulation of nuclear histone acetyltransferase by nucleic acids,histone · DNA complex,and chromatin

Nuclear histone acetyltransferase is found to be inhibited by various nucleic acids and components. Of the adenosine phosphates, the order of inhibitory potency is ATP>ADP>AMP. Among the nucleoside triphosphates, GTP seems to be the best inhibitor, followed by ATP, CTP, and UTP. Deoxymononucleotides have the same order of inhibition potential as their ribonucleotide counterparts, with inhibition constants in the low millimolar range. Oligonucleotides and polynucleotides are much better inhibitors than mononucleotides. The inhibition constants of the DNA molecules are size dependent. Molecules larger than 40 base pairs have inhibition constants less than 18 µg/ml, whereas molecules with decreasing numbers of base pairs have increasing magnitudes of inhibition constants. However, acetyltransferase has a lower affinity for free DNA molecules than for DNA · histone complexes as revealed by its interaction with DNA-Sepharose and histone · DNA-Sepharose columns. Furthermore, native chromatin depleted of endogenous histone acetyltransferase activity shows no inhibitory effect on the enzyme. Yet heated chromatin not only loses substrate activity but also becomes an inhibitor for the enzyme. Since unmodified sea urchin sperm chromatin has been shown to be a potent acetyltransferase inhibitor, it seems possible that DNA · histone complexes may be the true inhibitory species and that the conformational states of such complexes may serve as a regulatory mechanism in the control of the enzyme activity.     

Binding of the antitumor drug nogalamycin and its derivatives to DNA: structural comparison

The three-dimensional molecular structures of the complexes between a novel antitumor drug nogalamycin and its derivative U-58872 with a modified DNA hexamer d[m5CGT(pS)Am5CG] have been determined at 1.7- and 1.8-A resolution, respectively, by X-ray diffraction analyses. Both structures (in space group P6(1)) have been refined with constrained refinement procedure to final R factors of 0.208 (3386 reflections) and 0.196 (2143 reflections). In both complexes, two nogalamycins bind to the DNA hexamer double helix in a 2:1 ratio with the elongated aglycon chromophore intercalated between the CpG steps at both ends of the helix. The aglycon chromophore spans across the GC Watson-Crick base pairs with its nogalose lying in the minor groove and the aminoglucose lying in the major groove of the distorted B-DNA double helix. Most of the sugars remain in the C2'-endo pucker family, except three deoxycytidine residues (terminal C1, C7, and internal C5). All nucleotides are in the anti conformation. Specific hydrogen bonds are found in the complex between the drug and guanine-cytosine bases in both grooves of the helix. One hydroxyl group of the aminoglucose donates a hydrogen bond to the N7 of guanine, while the other receives a hydrogen bond from the N4 amino group of cytosine. The orientation of these two hydrogen bonds suggests that nogalamycin prefers a GC base pair with its aglycon chromophore intercalating at the 5'-side of a guanine (between NpG), or at the 3'-side of a cytosine (between CpN) with the sugars pointing toward the GC base pair. The binding of nogalamycin to DNA requires that the base pairs in DNA open up transiently to allow the bulky sugars to go through, suggesting that nogalamycin prefers GC sequences embedded in a stretch of AT sequences.