Quantum clustering and network analysis of MD simulation trajectories to probe the conformational ensembles of protein-ligand interactions

In this article, we present a novel application of a quantum clustering (QC) technique to objectively cluster the conformations, sampled by molecular dynamics simulations performed on different ligand bound structures of the protein. We further portray each conformational population in terms of dynamically stable network parameters which beautifully capture the ligand induced variations in the ensemble in atomistic detail. The conformational populations thus identified by the QC method and verified by network parameters are evaluated for different ligand bound states of the protein pyrrolysyl-tRNA synthetase (DhPylRS) from D. hafniense. The ligand/environment induced re-distribution of protein conformational ensembles forms the basis for understanding several important biological phenomena such as allostery and enzyme catalysis. The atomistic level characterization of each population in the conformational ensemble in terms of the re-orchestrated networks of amino acids is a challenging problem, especially when the changes are minimal at the backbone level. Here we demonstrate that the QC method is sensitive to such subtle changes and is able to cluster MD snapshots which are similar at the side-chain interaction level. Although we have applied these methods on simulation trajectories of a modest time scale (20 ns each), we emphasize that our methodology provides a general approach towards an objective clustering of large-scale MD simulation data and may be applied to probe multistate equilibria at higher time scales, and to problems related to protein folding for any protein or protein-protein/RNA/DNA complex of interest with a known structure.

DNA binding properties of novel dansylated distamycin analogues in which the fluorophore is directly conjugated to the N-methyl-pyrrole carboxamide backbone

Polyamides that are structural analogues of the naturally occurring DNA minor groove binding antibiotic distamycin (Dst) are promising candidates as gene modulators. Developing strategies for the large scale screening and monitoring of the cellular distribution of such ligands would aid the faster discovery of molecules, which would have eventual utility in molecular biology and medicine. Attachment of fluorescent tags would be a useful step towards this end. A fundamental question in this connection is whether the tag modifies the DNA binding affinity of the parent compounds. Towards answering this question, we have developed two oligopeptides that bear the dansyl (N, N-dimethylaminonaphthalene sulfonamido fluorophore) coupled directly to the N-terminus of the conjugated N-methylpyrrole carboxamide network, and possess three or four N-methyl pyrrole carboxamide units (abbreviated as Dn3 and Dn4 respectively). DNA binding abilities of these molecules were assessed from fluorescence titration experiments, duplex-DNA T-m analysis (employing both UV and fluorescence spectroscopy), induced circular dichroism measurements (ICD), salt dependence of ICD and apparent binding constant measurements (K-app) employing ethidium bromide (EtBr) displacement assay. Both these molecules 'reported' DNA binding in the form of an enhanced fluorescence emission. As judged from the ICD measurements, salt dependence of ICD, T-m analysis and K-app measurements, the binding affinities of the molecules that possessed dansyl group at their N-termini were lower than the ones with equivalent number of amide units, but possessed N-methylpyrrole carboxamide unit at their N-termini. These results would have implications in the future design of fluorescent polyamides.

Structure, function, and mechanism of HhaI DNA methyltransferases

A vast amount of literature has accumulated on the characterization of DNA methyltransferases. The HhaI DNA methyltransferase, a C5-cytosine methyltransferase, has been the subject of investigation for the last 2 decades. Biochemical and kinetic characterization have led to an understanding of the catalytic and kinetic mechanism of the methyltransfer reaction. The HhaI methyltransferase has also been subjected to extensive structural analysis, with the availability of 12 structures with or without a cofactor and a variety of DNA substrates. The mechanism of base flipping, first described for the HhaI methyltransferase, is conserved among all DNA methyltransferases and is also found to occur in numerous DNA repair enzymes. Studies with other methyltransferase reveal a significant structural and functional similarity among different types of methyltransferases. This review aims to summarize the available information on the HhaI DNA methyltransferase.

Dynamic and thermodynamic anomalies of water at low temperatures: from bulk water to reverse micelles and DNA hydration layer

Liquid water is known to exhibit remarkable thermodynamic and dynamic anomalies, ranging from solvation properties in supercritical state to an apparent divergence of the linear response functions at a low temperature. Anomalies in various dynamic properties of water have also been observed in the hydration layer of proteins, DNA grooves and inside the nanocavity, such as reverse micelles and nanotubes. Here we report studies on the molecular origin of these anomalies in supercooled water, in the grooves of DNA double helix and reverse micelles. The anomalies have been discussed in terms of growing correlation length and intermittent population fluctuation of 4- and 5-coordinated species. We establish correlation between thermodynamic response functions and mean squared species number fluctuation. Lifetime analysis of 4- and 5-coordinated species reveals interesting differences between the role of the two species in supercooled and constrained water. The nature and manifestations of the apparent and much discussed liquid-liquid transition under confinement are found to be markedly different from that in the bulk. We find an interesting `faster than bulk' relaxation in reverse micelles which we attribute to frustration effects created by competition between the correlations imposed by surface interactions and that imposed by hydrogen bond network of water.

Structure of tandem and its implication for bifunctional intercalation into dna

Quinoxaline antibiotics (Fig. 1a, b) form a useful group of compounds for the study of drug–nucleic acid interactions1,2. They consist of a cross-bridged cyclic octadepsipeptide, variously modified, bearing two quinoxaline chromophores. These antibiotics intercalate bifunctionally into DNA2,3 probably via the narrow groove, forming a complex in which, most probably, two base pairs are sandwiched between the chromophores4,5. Depending on the nature of their sulphur-containing cross-bridge and modifications to their amino acid side chains, they display characteristic patterns of nucleotide sequence selectivity when binding to DNAs of different base composition and to synthetic polydeoxynucleotides4,6,7. This specificity has been tentatively ascribed to specific hydrogen-bonding interactions between functional groups in the DNA and complementary moieties on the peptide ring2,4,5. Variations in selectivity have been attributed both to changes in the conformation of the peptide backbone6 and no modifications of the cross-bridge7. These suggestions were made, however, in the absence of firm knowledge about the three-dimensional structure and conformation of the antibiotic molecules. We now report the X-ray structure analysis of the synthetic analogue of the antibiotic triostin A, TANDEM (des-N-tetramethyl triostin A) (Fig. 1c), which binds preferentially to alternating adenine-thymine sequences7. The X-ray structure provides a starting point for exploring the origin of this specificity and suggests possible models for the binding of other members of the quinoxaline series.

Synthesis and DNA cleavage studies of novel quinoline oxime esters

New 2-chloro-3-formyl quinoline oxime esters were synthesized by the reaction of 2-chloro-3-formyl quinoline oximes with various benzoyl chlorides in the presence of triethyl amine and dichloromethane at 0 degrees C. The DNA photo cleavage studies of some new oxime esters were investigated by neutral agarose gel electrophoresis at different concentrations (40 mu M and 80 mu M). Analysis of the cleavage products in agarose gel indicated that few of quinoline oxime esters (3d-i) converted into supercoiled pUC19 plasmid DNA to its nicked or linear form. (C) 2011 Elsevier Ltd. All rights reserved.

Photo-induced DNA cleavage activity and remarkable photocytotoxicity of lanthanide(III) complexes of a polypyridyl ligand

Lanthanide(III) complexes [Ln(pyphen)(acac)(2)(NO3)] (1, 2), [Ln(pydppz)(acac)(2)(NO3)] (3, 4) and [La(pydppz)(anacac)(2)(NO3)] (5), where Ln is La(III) (in 1, 3, 5) and Gd(III) (in 2, 4), pyphen is 6-(2-pyridyl)-1,10-phenanthroline, pydppz is 6-(2-pyridyl)-dipyrido[3,2-a:2',3'-c] phenazine, anacac is anthracenylacetylacetonate and acac is acetylacetonate, were prepared, characterized and their DNA photocleavage activity and photocytotoxicity studied. The crystal structure of complex 2 displays a GdO6N3 coordination. The pydppz complexes 3-5 show an electronic spectral band at similar to 390 nm in DMF. The La(III) complexes are diamagnetic, while the Gd(III) complexes are paramagnetic with seven unpaired electrons. The molar conductivity data suggest 1 : 1 electrolytic nature of the complexes in aqueous DMF. They are avid binders to calf thymus DNA giving K-b in the range of 5.4 10(4)-1.2 x 10(6) M-1. Complexes 3-5 efficiently cleave supercoiled DNA to its nicked circular form in UV-A light of 365 nm via formation of singlet oxygen (O-1(2)) and hydroxyl radical (HO center dot) species. Complexes 3-5 also exhibit significant photocytotoxic effect in HeLa cancer cells giving respective IC50 value of 0.16(+/- 0.01), 0.15(+/- 0.01) and 0.26 +/-(0.02) mu M in UV-A light of 365 nm, while they are less toxic in dark with an IC50 value of >3 mu M. The presence of an additional pyridyl group makes the pydppz complexes more photocytotoxic than their dppz analogues. FACS analysis of the HeLa cells treated with complex 4 shows apoptosis as the major pathway of cell death. Nuclear localization of complex 5 having an anthracenyl moiety as a fluorophore is evidenced from the confocal microscopic studies.

IN-VITRO CYTOTOXICITY AND CELL CYCLE ANALYSIS OF TWO NOVEL BIS-1, 2, 4-TRIAZOLE DERIVATIVES: 1,4-BIS[5-(5-MERCAPTO-1,3,4-OXADIAZOL-2-yl-METHYL)-THIO-4-(p-TOLYL)-1, 2,4-TRIAZOL-3-yl]-BUTANE (MNP-14) AND 1,4-BIS[5-(CARBETHOXY-METHYL)-THIO-4-(p-ETHOXY PHENYL)-1,2,4-TRIAZOL-3-yl]-BUTANE (MNP-16)

In the present study, we have tested the cytotoxic and DNA damage activity of two novel bis-1,2,4 triazole derivatives, namely 1,4-bis[5-(5-mercapto-1,3,4-oxadiazol-2-yl-methyl)-thio4-(p-tolyl)-1,2 ,4-triazol-3-yl]-butane (MNP-14) and 1,4-bis[5-(carbethoxy-methyl)-thio-4-(p-ethoxy phenyl) -1,2,4-triazol-3-yl]-butane (MNP-16). The effect of these molecules on cellular apoptosis was also determined. The in-vitro cytotoxicity was evaluated by a 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyl tetrazolium bromide (MTT) assay as well as Trypan blue dye exclusion methods against human acute lymphoblastic leukemia (MOLT4) and lung cancer cells (A549). Our results showed that MNP-16 induced significant cytotoxicity (IC50 of 3-5 mu M) compared with MNP-14. The cytotoxicity induced by MNP-16 was time and concentration dependent. The cell cycle analysis by flow cytometry (fluorescence-activated cell sorting [FACS]) revealed that though there was a significant increase in the apoptotic population (sub-G1 phase) with an increased concentration of MNP-14 and 16, there was no cell cycle arrest. Further, the comet assay results indicated considerable DNA

Tuning DNA-amphiphile condensate architecture with strongly binding counterions

Electrostatic self-assembly of colloidal and nanoparticles has attracted a lot of attention in recent years, since it offers the possibility of producing novel crystalline structures that have the potential to be used as advanced materials for photonic and other applications. The stoichiometry of these crystals is not constrained by charge neutrality of the two types of particles due to the presence of counterions, and hence a variety of three-dimensional structures have been observed depending on the relative sizes of the particles and their charge. Here we report structural polymorphism of two-dimensional crystals of oppositely charged linear macroions, namely DNA and self-assembled cylindrical micelles of cationic amphiphiles. Our system differs from those studied earlier in terms of the presence of a strongly binding counterion that competes with DNA to bind to the micelle. The presence of these counterions leads to novel structures of these crystals, such as a square lattice and a root 3 x root 3 superlattice of an underlying hexagonal lattice, determined from a detailed analysis of the small-angle diffraction data. These lower-dimensional equilibrium systems can play an important role in developing a deeper theoretical understanding of the stability of crystals of oppositely charged particles. Further, it should be possible to use the same design principles to fabricate structures on a longer length-scale by an appropriate choice of the two macroions.

Evaluation of DNA binding, antioxidant and cytotoxic activity of mononuclear Co(III) complexes of 2-oxo-1,2-dihydrobenzoh]quinoline-3-carbaldehyde thiosemicarbazones

Four new 2-oxo-1,2-dihydrobenzoh]quinoline-3-carbaldehyde N-substituted thiosemicarbazone ligands (H-2-LR, where R = H, Me, Et or Ph) and their corresponding new cobalt(III) complexes have been synthesized and characterized. The structures of the complexes 2 and 3 were determined by single crystal X-ray diffraction analysis. The interactions of the new complexes with DNA were investigated by absorption, emission and viscosity studies which indicated that the complexes bind to DNA via intercalation. Antioxidant studies of the new complexes showed that the significant antioxidant activity against DPPH radical. In addition, the in vitro cytotoxicity of complexes 1-4 against A549 cell line was assayed which showed higher cytotoxic activity with lower IC50 values indicating their efficiency in killing the cancer cells even at very low concentrations. (C) 2012 Elsevier Masson SAS. All rights reserved.

Distinct mechanisms of DNA repair in mycobacteria and their implications in attenuation of the pathogen growth

About a third of the human population is estimated to be infected with Mycobacterium tuberculosis. Emergence of drug resistant strains and the protracted treatment strategies have compelled the scientific community to identify newer drug targets, and to develop newer vaccines. In the host macrophages, the bacterium survives within an environment rich in reactive nitrogen and oxygen species capable of damaging its genome. Therefore, for its successful persistence in the host, the pathogen must need robust DNA repair mechanisms. Analysis of M. tuberculosis genome sequence revealed that it lacks mismatch repair pathway suggesting a greater role for other DNA repair pathways such as the nucleotide excision repair, and base excision repair pathways. In this article, we summarize the outcome of research involving these two repair pathways in mycobacteria focusing primarily on our own efforts. Our findings, using Mycobacterium smegmatis model, suggest that deficiency of various DNA repair functions in single or in combinations severely compromises their DNA repair capacity and attenuates their growth under conditions typically encountered in macrophages. (C) 2011 Elsevier Ireland Ltd. All rights reserved.

Genetic analysis of an Indian family with members affected with Waardenburg syndrome and Duchenne muscular dystrophy

Purpose: Waardenburg syndrome (WS) is characterized by sensorineural hearing loss and pigmentation defects of the eye, skin, and hair. It is caused by mutations in one of the following genes: PAX3 (paired box 3), MITF (microphthalmia-associated transcription factor), EDNRB (endothelin receptor type B), EDN3 (endothelin 3), SNAI2 (snail homolog 2, Drosophila) and SOX10 (SRY-box containing gene 10). Duchenne muscular dystrophy (DMD) is an X-linked recessive disorder caused by mutations in the DMD gene. The purpose of this study was to identify the genetic causes of WS and DMD in an Indian family with two patients: one affected with WS and DMD, and another one affected with only WS. Methods: Blood samples were collected from individuals for genomic DNA isolation. To determine the linkage of this family to the eight known WS loci, microsatellite markers were selected from the candidate regions and used to genotype the family. Exon-specific intronic primers for EDN3 were used to amplify and sequence DNA samples from affected individuals to detect mutations. A mutation in DMD was identified by multiplex PCR and multiplex ligation-dependent probe amplification method using exon-specific probes. Results: Pedigree analysis suggested segregation of WS as an autosomal recessive trait in the family. Haplotype analysis suggested linkage of the family to the WS4B (EDN3) locus. DNA sequencing identified a novel missense mutation p.T98M in EDN3. A deletion mutation was identified in DMD. Conclusions: This study reports a novel missense mutation in EDN3 and a deletion mutation in DMD in the same Indian family. The present study will be helpful in genetic diagnosis of this family and increases the mutation spectrum of EDN3.

Photocytotoxic Oxidovanadium(IV) Complexes of Polypyridyl Ligands Showing DNA-Cleavage Activity in Near-IR Light

Oxidovanadium(IV) complexes VO(pyphen)(L)]Cl2 (1, 2) and VO(pydppz)(L)]Cl2 (3, 4), where L is 1,10-phenanthroline (phen in 1 and 3) and dipyrido3,2-a:2',3'-c]phenazine (dppz in 2 and 4) are prepared and characterized. The crystal structure of VO(pyphen)(phen)](ClO4)2 (1a) shows a six-coordinate VN5O geometry with a VO2+ moiety in which the polypyridyl ligand binds in a meridional fashion and the phen ligand displays a chelating binding mode with an N-donor site trans to the oxidovanadyl group. The complexes show a dd band within 720-750 nm in DMF. The one-electron paramagnetic complexes are 1:2 electrolytes in DMF. The complexes exhibit an irreversible VIV/VIII redox response near -0.85 V vs. SCE in DMF/0.1 M TBAP. The complexes bind to calf thymus (ct) DNA giving Kb values within 7.5 x 104 to 1.1 x 106 M1. The complexes show poor chemical nuclease activity in the dark and exhibit significant DNA-photocleaving activity in near-IR light of 705 and 785 nm forming .OH radicals. Complexes 2-4 show remarkable photocytotoxicity in HeLa cancer cells. FACS analysis of the HeLa cells treated with complex 4 shows cell death as highlighted by the sub G1 peak. Propidium iodide staining data indicate apoptosis as the primary mode of cell death.

Under-representation of intrinsic terminators across bacterial genomic islands: Rho as a principal regulator of xenogenic DNA expression

Two transcription termination mechanisms - intrinsic and Rho-dependent - have evolved in bacteria. The Rho factor occurs in most bacterial lineages, and has been hypothesized to play a global regulatory role. Genome-wide studies using microarray, 2D-gel electrophoresis and ChIP-chip provided evidence that Rho serves to silence transcription from horizontally acquired genes and prophages in Escherichia coli K-12, implicating the factor to be a part of the ``cellular immune mechanism'' protecting against deleterious phages and aberrant gene expression from acquired xenogenic DNA. We have investigated this model by adopting an alternate in silico approach and have extended the study to other species. Our analysis shows that several genomic islands across diverse phyla have under-representation of intrinsic terminators, similar to that experimentally observed in E. coli K-12. This implies that Rho-dependent termination is the predominant process operational in these islands and that silencing of foreign DNA is a conserved function of Rho. From the present analysis, it is evident that horizontally acquired islands have lost intrinsic terminators to facilitate Rho-dependent termination. These results underscore the importance of Rho as a conserved, genome-wide sentinel that regulates potentially toxic xenogenic DNA. (C) 2012 Elsevier B.V. All rights reserved.

Role of Specific Cations and Water Entropy on the Stability of Branched DNA Motif Structures

DNA three-way junctions (TWJs) are important intermediates in various cellular processes and are the simplest of a family of branched nucleic acids being considered as scaffolds for biomolecular nanotechnology. Branched nucleic acids are stabilized by divalent cations such as Mg2+, presumably due to condensation and neutralization of the negatively charged DNA backbone. However, electrostatic screening effects point to more complex solvation dynamics and a large role of interfacial waters in thermodynamic stability. Here, we report extensive computer simulations in explicit water and salt on a model TWJ and use free energy calculations to quantify the role of ionic character and strength on stability. We find that enthalpic stabilization of the first and second hydration shells by Mg2+ accounts for 1/3 and all of the free energy gain in 50% and pure MgCl2 solutions, respectively. The more distorted DNA molecule is actually destabilized in pure MgCl2 compared to pure NaCl. Notably, the first shell, interfacial waters have very low translational and rotational entropy (i.e., mobility) compared to the bulk, an entropic loss that is overcompensated by increased enthalpy from additional electrostatic interactions with Mg2+. In contrast, the second hydration shell has anomalously high entropy as it is trapped between an immobile and bulklike layer. The nonmonotonic entropic signature and long-range perturbations of the hydration shells to Mg2+ may have implications in the molecular recognition of these motifs. For example, we find that low salt stabilizes the parallel configuration of the three-way junction, whereas at normal salt we find antiparallel configurations deduced from the NMR. We use the 2PT analysis to follow the thermodynamics of this transition and find that the free energy barrier is dominated by entropic effects that result from the decreased surface area of the antiparallel form which has a smaller number of low entropy waters in the first monolayer.