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.

Sequence-dependent conformation of an A-DNA double helix: The crystal structure of the octamer d(G-G-T-A-T-A-C-C)

The crystal structures of the synthetic self-complementary octamer d(G-G-T-A-T-A-C-C) and its 5-bromouracil-containing analogue have been refined to R values of 20% and 14% at resolutions of 1·8 and 2·25 Å, respectively. The molecules adopt an A-DNA type double-helical conformation, which is minimally affected by crystal forces. A detailed analysis of the structure shows a considerable influence of the nucleotide sequence on the base-pair stacking patterns. In particular, the electrostatic stacking interactions between adjacent guanine and thymine bases produce symmetric bending of the double helix and a major-groove widening. The sugar-phosphate backbone appears to be only slightly affected by the base sequence. The local variations in the base-pair orientation are brought about by correlated adjustments in the backbone torsion angles and the glycosidic orientation. Sequence-dependent conformational variations of the type observed here may contribute to the specificity of certain protein-DNA interactions.

Chemistry and biology of DNA methyltransferases

Recognition of a specific DNA sequence by a protein is probably the best example of macromolecular interactions leading to various events. It is a prerequisite to understanding the basis of protein-DNA interactions to obtain a better insight into fundamental processes such as transcription, replication, repair, and recombination. DNA methyltransferases with varying sequence specificities provide an excellent model system for understanding the molecular mechanism of specific DNA recognition. Sequence comparison of cloned genes, along with mutational analyses and recent crystallographic studies, have clearly defined the functions of various conserved motifs. These enzymes access their target base in an elegant manner by flipping it out of the DNA double helix. The drastic protein-induced DNA distortion, first reported for HhaI DNA methyltransferase, appears to be a common mechanism employed by various proteins that need to act on bases. A remarkable feature of the catalytic mechanism of DNA (cytosine-5) methyltransferases is the ability of these enzymes to induce deamination of the target cytosine in the absence of S-adenosyl-L-methionine or its analogs. The enzyme-catalyzed deamination reaction is postulated to be the major cause of mutational hotspots at CpG islands responsible for various human genetic disorders. Methylation of adenine residues in Escherichia coli is known to regulate various processes such as transcription, replication, repair, recombination, transposition, and phage packaging.

New evidences on Tau-DNA interactions and relevance to neurodegeneration

Tau is mainly distributed in cytoplasm and also found to be localized in the nucleus. There is limited data on DNA binding potential of Tau.We provide novel evidence on nicking of DNA by Tau. Tau nicks the supercoiled DNA leading to open circular and linear forms. The metal ion magnesium (a co-factor for endonuclease) enhanced the Tau DNA nicking ability, while an endonuclease specific inhibitor,aurinetricarboxylic acid (ATA) inhibited the Tau DNA nicking ability Further, we also evidenced that Tau induces B-C-A mixed conformational transition in DNA and also changes DNA stability. Tau-scDNA complex is more sensitive to DNAse I digestion indicating stability changes in DNA caused by Tau. These findings indicate that Tau alters DNA helicity and integrity and also nicks the DNA. The relevance of these novel intriguing findings regarding the role Tau in neuronal dysfunction is discussed. (C) 2010 Elsevier Ltd. All rights reserved.

Interaction of synthetic analogs of distamycin and netropsin with nucleic acids. Does curvature of ligand play a role in distamycin-DNA interactions?

Distamycin and netropsin, a class of minor groove binding nonintercalating agents, are characterized by their B-DNA and A-T basespecific interactions. To understand the CQI I ~OIT~ ~ I ~ ~aOnMd ~c hemical basis of the above specificities, the DNA-binding characteristics of a novel synthetic analogue of distamycin have been studied. The analogue, mPD derivative, has the requisite charged end groups and a number of potential hydrogen-bonding loci equal to those of distamycin. The difference in the backbone curvatures of the ligands, distamycin, the mPD derivative, and NSC 101327 (another structurally analogous compound),is a major difference between these ligands. UV and CD spectrosoopic studies reported here show the following salient features: The mPD derivative recognizes only B-DNA, to which it binds via the minor groove. On the other hand, unlike distamycin, it binds with comparable affinities to A-T and G-C base pairs in a natural DNA. These DNA-binding properties are compared with those reported earlier for distamycin and NSC 101327 [Zimmer, Ch., & Wahnert, U. (1986) Prog. Biophys. Mol. Biol. 47, 31-1121. The backbone structures of these three ligands were compared to show the progressive decrease in curvatures in the order distamycin, mPD derivative, and NSC 101327. The plausible significance of the backbone curvature vis-&vis the characteristic B-DNA and AT-specific binding of distamycin is discussed. To our knowledge, this is the first attempt (with a model synthetic analogue) to probe the possible influence of backbone curvature upon the specificity of interactions of the distamycin class of groove-binding ligands with DNA.

Theoretical studies on alpha-helix--DNA interactions

The interaction energies between (Ala)10 and alpha-helix fragment and different nucleotide sequences in right-handed B-form have been optimized using semi-empirical potential energy functions. The energies are calculated for two different orientations of the alpha-helix, viz., when the alpha-helix axis taken in the N----C direction is (i) parallel and (ii) antiparallel to the 5'-3' ascending strand of DNA, proximal to it. When both the DNA molecule as well as the alpha-helix are treated as rigid molecules it is found that a polyalanine alpha-helix has slightly more favourable contacts when it is in the proximity of a four nucleotide sequence of 5'-(N-A-T-N)-3' type, where N is either a purine or a pyrimidine. However, when the two interacting molecules are allowed to undergo local structural variations then the interaction energy appears to be independent of the base sequence confirming the non-specific nature of these interactions.

Photoresist functionalisation method for high-density protein microarrays using photolithography

Since the last decade, there is a growing need for patterned biomolecules for various applications ranging from diagnostic devices to enabling fundamental biological studies with high throughput. Protein arrays facilitate the study of protein-protein, protein-drug or protein-DNA interactions as well as highly multiplexed immunosensors based on antibody-antigen recognition. Protein microarrays are typically fabricated using piezoelectric inkjet printing with resolution limit of similar to 70-100 mu m limiting the array density. A considerable amount of research has been done on patterning biomolecules using customised biocompatible photoresists. Here, a simple photolithographic process for fabricating protein microarrays on a commercially available diazo-naphthoquinone-novolac-positive tone photoresist functionalised with 3-aminopropyltriethoxysilane is presented. The authors demonstrate that proteins immobilised using this procedure retain their activity and therefore form functional microarrays with the array density limited only by the resolution of lithography, which is more than an order of magnitude compared with inkjet printing. The process described here may be useful in the integration of conventional semiconductor manufacturing processes with biomaterials relevant for the creation of next-generation bio-chips.

Sensitivity of polarization fluctuations to the nature of protein-water interactions: Study of biological water in four different protein-water systems

Since the time of Kirkwood, observed deviations in magnitude of the dielectric constant of aqueous protein solution from that of neat water (similar to 80) and slower decay of polarization have been subjects of enormous interest, controversy, and debate. Most of the common proteins have large permanent dipole moments (often more than 100 D) that can influence structure and dynamics of even distant water molecules, thereby affecting collective polarization fluctuation of the solution, which in turn can significantly alter solution's dielectric constant. Therefore, distance dependence of polarization fluctuation can provide important insight into the nature of biological water. We explore these aspects by studying aqueous solutions of four different proteins of different characteristics and varying sizes, chicken villin headpiece subdomain (HP-36), immunoglobulin binding domain protein G (GB1), hen-egg white lysozyme (LYS), and Myoglobin (MYO). We simulate fairly large systems consisting of single protein molecule and 20000-30000 water molecules (varied according to the protein size), providing a concentration in the range of similar to 2-3 mM. We find that the calculated dielectric constant of the system shows a noticeable increment in all the cases compared to that of neat water. Total dipole moment auto time correlation function of water < dM(W) (0)delta M-W (t) > is found to be sensitive to the nature of the protein. Surprisingly, dipole moment of the protein and total dipole moment of the water molecules are found to be only weakly coupled. Shellwise decomposition of water molecules around protein reveals higher density of first layer compared to the succeeding ones. We also calculate heuristic effective dielectric constant of successive layers and find that the layer adjacent to protein has much lower value (similar to 50). However, progressive layers exhibit successive increment of dielectric constant, finally reaching a value close to that of bulk 4-5 layers away. We also calculate shellwise orientational correlation function and tetrahedral order parameter to understand the local dynamics and structural re-arrangement of water. Theoretical analysis providing simple method for calculation of shellwise local dielectric constant and implication of these findings are elaborately discussed in the present work. (C) 2014 AIP Publishing LLC.

Kinetics of interaction of Cotton Leaf Curl Kokhran Virus-Dabawali (CLCuKV-Dab) coat protein and its mutants with ssDNA

Gemini viral assembly and transport of viral DNA into nucleus for replication, ssentially involve DNA-coat protein interactions. The kinetics of interaction of Cotton LeafCtirl Kokhran Virus-Dabawali recombinant coat protein (rCP) with DNA was studied by electrophoretic mobility shift assay (EMSA) and Surface plasmon resonance (SPR). The rCP interacted with ssDNA with a K-A, of 2.6 +/- 0.29 x 10(8) M-1 in a sequence non-specific manner. The CP has a conserved C2H2 type zinc finger motif composed of residues C68, C72, H81 and H85. Mutation of these residues to alanine resulted in reduced binding to DNA probes. The H85A mutant rCP showed the least binding with approximately 756 fold loss in the association rate and a three order magnitude decrease in the binding affinity as compared to rCP. The CP-DNA interactions via the zinc finger motif could play a Crucial role ill Virus assembly and in nuclear transport. (C) 2009 Elsevier Inc.

Use of structure-directed DNA ligands to probe the binding of recA protein to narrow and wide grooves of DNA and on its ability to promote homologous pairing

We have used circular dichroism and structure-directed drugs to identify the role of structural features, wide and narrow grooves in particular, required for the cooperative polymerization, recognition of homologous sequences, and the formation of joint molecules promoted by recA protein. The path of cooperative polymerization of recA protein was deduced by its ability to cause quantitative displacement of distamycin from the narrow groove of duplex DNA. By contrast, methyl green bound to the wide groove was retained by the nucleoprotein filaments comprised of recA protein-DNA. Further, the mode of binding of these ligands and recA protein to DNA was confirmed by DNaseI digestion. More importantly, the formation of joint molecules was prevented by distamycin in the narrow groove while methyl green in the wide groove had no adverse effect. Intriguingly, distamycin interfered with the production of coaggregates between nucleoprotein filaments of recA protein-M13 ssDNA and naked linear M13 duplex DNA, but not with linear phi X174 duplex DNA. Thus, these data, in conjunction with molecular modeling, suggest that the narrow grooves of duplex DNA provide the fundamental framework required for the cooperative polymerization of recA protein and alignment of homologous sequences. These findings and their significance are discussed in relation to models of homologous pairing between two intertwined DNA molecules.

Single-Molecule DNA Analysis Reveals That Yeast Hop1 Protein Promotes DNA Folding and Synapsis: Implications for Condensation of Meiotic Chromosomes

During meiosis, long-range interaction between homologous chromosomes is thought to be crucial for homology recognition, exchange of DNA strands, and production of normal haploid gametes. However, little is known about the identity of the proteins involved and the actual molecular mechanism(s) by which chromosomes recognize and recombine with their appropriate homologous partners. Single-molecule analyses have the potential to provide insights into our understanding of this fascinating and long-standing question. Using atomic force microscopy and magnetic tweezers techniques, we discovered that Hop1 protein, a key structural component of Saccharomyces cerevisiae synaptonemal complex, exhibits the ability to bridge noncontiguous DNA segments into intramolecular stem-loop structures in which the DNA segments appear to be fully synapsed within the filamentous protein stems. Additional evidence suggests that Hop1 folds DNA into rigid protein DNA filaments and higher-order nucleoprotein structures. Importantly, Hop1 promotes robust intra- and intermolecular synapsis between double-stranded DNA molecules, suggesting that juxtaposition of DNA sequences may assist in strand exchange between homologues by recombination-associated proteins. Finally, the evidence from ensemble experiments is consistent with the notion that Hop1 causes rigidification of DNA molecules. These results provide the first direct evidence for long-range protein-mediated DNA DNA synapsis, independent of crossover recombination, which is presumed to occur during meiotic recombination.

DNA looping and translocation provide an optimal cleavage mechanism for the type III restriction enzymes

EcoP15I is a type III restriction enzyme that requires two recognition sites in a defined orientation separated by up to 3.5 kbp to efficiently cleave DNA. The mechanism through which site- bound EcoP15I enzymes communicate between the two sites is unclear. Here, we use atomic force microscopy to study EcoP15I-DNA pre-cleavage complexes. From the number and size distribution of loops formed, we conclude that the loops observed do not result from translocation, but are instead formed by a contact between site- bound EcoP15I and a nonspecific region of DNA. This conclusion is confirmed by a theoretical polymer model. It is further shown that translocation must play some role, because when translocation is blocked by a Lac repressor protein, DNA cleavage is similarly blocked. On the basis of these results, we present a model for restriction by type III restriction enzymes and highlight the similarities between this and other classes of restriction enzymes.