Molecular modeling of drug-DNA interactions: virtual screening to structure-based design

Virtual ligand screening (VLS) and structure-based design are strategies that have been routinely used for the development of pharmaceuticals, particularly those targeting enzymes and other protein targets. In recent years, an increased understanding of the role played by nucleic acids in biological systems made DNA an alternative candidate for the development of new drugs. This review highlights some successful applications of molecular modeling in virtual ligand screening and structure-based design of organic and inorganic molecules that target non-canonical nucleic acid structures such as G-quadruplex and triplex DNA.     

Accurate ab initio prediction of NMR chemical shifts of nucleic acids and nucleic acids/protein complexes

NMR chemical shift predictions based on empirical methods are nowadays indispensable tools during resonance assignment and 3D structure calculation of proteins. However, owing to the very limited statistical data basis, such methods are still in their infancy in the field of nucleic acids, especially when non-canonical structures and nucleic acid complexes are considered. Here, we present an ab initio approach for predicting proton chemical shifts of arbitrary nucleic acid structures based on state-of-the-art fragment-based quantum chemical calculations. We tested our prediction method on a diverse set of nucleic acid structures including double-stranded DNA, hairpins, DNA/protein complexes and chemically-modified DNA. Overall, our quantum chemical calculations yield highly/very accurate predictions with mean absolute deviations of 0.3–0.6 ppm and correlation coefficients (r²) usually above 0.9. This will allow for identifying misassignments and validating 3D structures. Furthermore, our calculations reveal that chemical shifts of protons involved in hydrogen bonding are predicted significantly less accurately. This is in part caused by insufficient inclusion of solvation effects. However, it also points toward shortcomings of current force fields used for structure determination of nucleic acids. Our quantum chemical calculations could therefore provide input for force field optimization.     

Predicting nucleic acid binding interfaces from structural models of proteins

The function of DNA‐ and RNA‐binding proteins can be inferred from the characterization and accurate prediction of their binding interfaces. However, the main pitfall of various structure‐based methods for predicting nucleic acid binding function is that they are all limited to a relatively small number of proteins for which high‐resolution three‐dimensional structures are available. In this study, we developed a pipeline for extracting functional electrostatic patches from surfaces of protein structural models, obtained using the I‐TASSER protein structure predictor. The largest positive patches are extracted from the protein surface using the patchfinder algorithm. We show that functional electrostatic patches extracted from an ensemble of structural models highly overlap the patches extracted from high‐resolution structures. Furthermore, by testing our pipeline on a set of 55 known nucleic acid binding proteins for which I‐TASSER produces high‐quality models, we show that the method accurately identifies the nucleic acids binding interface on structural models of proteins. Employing a combined patch approach we show that patches extracted from an ensemble of models better predicts the real nucleic acid binding interfaces compared with patches extracted from independent models. Overall, these results suggest that combining information from a collection of low‐resolution structural models could be a valuable approach for functional annotation. We suggest that our method will be further applicable for predicting other functional surfaces of proteins with unknown structure. Proteins 2012. © 2011 Wiley Periodicals, Inc.     

Generation of amplifiable genome-specific oligonucleotide probes and libraries

Here we describe a process for the generation of oligonucleotide libraries representative of a given nucleic acid. Starting from at random pool of DNA oligonucleotides, the technique selects only those that hybridize to the nucleic acid template. This selection yields a highly specific library that represents an oligonucleotide image of the chosen template. The novel quality of this approach is the generation of amplifiable oligonucleotide probes that are of unique length and are easily subjected to differential selection. Here we apply this technique to produce different genomic oligonucleotide libraries and show that these genomic oligonucleotide libraries do not cross-hybridize. Differential selection of these genomic oligonucleotide libraries produces oligonucleotides that can be used in the identification, characterzation, and isolation of nucleic acids.     

Single-round deoxyribozyme discovery

Artificial evolution experiments typically use libraries of ∼10¹⁵ sequences and require multiple rounds of selection to identify rare variants with a desired activity. Based on the simple structures of some aptamers and nucleic acid enzymes, we hypothesized that functional motifs could be isolated from significantly smaller libraries in a single round of selection followed by high-throughput sequencing. To test this idea, we investigated the catalytic potential of DNA architectures in which twelve or fifteen randomized positions were embedded in a scaffold present in all library members. After incubating in either the presence or absence of lead (which promotes the nonenzymatic cleavage of RNA), library members that cleaved themselves at an RNA linkage were purified by PAGE and characterized by high-throughput sequencing. These selections yielded deoxyribozymes with activities 8- to 30-fold lower than those previously isolated under similar conditions from libraries containing 10¹⁴ different sequences, indicating that the disadvantage of using a less diverse pool can be surprisingly small. It was also possible to elucidate the sequence requirements and secondary structures of deoxyribozymes without performing additional experiments. Due to its relative simplicity, we anticipate that this approach will accelerate the discovery of new catalytic DNA and RNA motifs.     

Describing RNA structure by libraries of clustered nucleotide doublets

The rapidly increasing wealth of structural information on RNA and knowledge of its varying roles in biology have facilitated the study of RNA structure using computational methods. Here, we present a new method to describe RNA structure based on nucleotide doublets, where a doublet is any two nucleotides in a structure. We restrict our search to doublets that are close together in space, but not necessarily in sequence, and obtain doublet libraries of various sizes by clustering a large set of doublets taken from a data set of high-resolution RNA structures. We demonstrate that these libraries are able to both capture structural features present in RNA and fit local RNA structure with a high level of accuracy. Libraries ranging in size from ten to 100 doublets are examined, and a detailed analysis shows that a library with as few as 30 doublets is sufficient to capture the most common structural features, while larger libraries would be more appropriate for accurate modeling. We anticipate many uses for these libraries, from annotation to structure refinement and prediction.     

Dynamics and Conformational Changes in Human NEIL2 DNA Glycosylase Analyzed by Hydrogen/Deuterium Exchange Mass Spectrometry

Base excision DNA repair (BER) is necessary for removal of damaged nucleobases from the genome and their replacement with normal nucleobases. BER is initiated by DNA glycosylases, the enzymes that cleave the N-glycosidic bonds of damaged deoxynucleotides. Human endonuclease VIII-like protein 2 (hNEIL2), belonging to the helix–two-turn–helix structural superfamily of DNA glycosylases, is an enzyme uniquely specific for oxidized pyrimidines in non-canonical DNA substrates such as bubbles and loops. The structure of hNEIL2 has not been solved; its closest homologs with known structures are NEIL2 from opossum and from giant mimivirus. Here we analyze the conformational dynamics of free hNEIL2 using a combination of hydrogen/deuterium exchange mass spectrometry, homology modeling and molecular dynamics simulations. We show that a prominent feature of vertebrate NEIL2 – a large insert in its N-terminal domain absent from other DNA glycosylases – is unstructured in solution. It was suggested that helix–two-turn–helix DNA glycosylases undergo open–close transition upon DNA binding, with the large movement of their N- and C-terminal domains, but the open conformation has been elusive to capture. Our data point to the open conformation as favorable for free hNEIL2 in solution. Overall, our results are consistent with the view of hNEIL2 as a conformationally flexible protein, which may be due to its participation in the repair of non-canonical DNA structures and/or to the involvement in functional and regulatory protein–protein interactions.     

In vitro selection of high-affinity nucleic acid ligands to parasite target molecules

The logic of using nucleic acids as pharmaceutical reagents is in part based on their capacity to interact with high affinity and specificity with other biological components. Considerable progress has been made over the past 10 years in the development of nucleic acid-based drug molecules using a variety of different technologies. One approach is a combinatorial technology that involves an iterative Darwinian-type in vitro evolution process, which has been termed SELEX for 'systematic evolution of ligands by exponential enrichment'. The procedure is a highly efficient method of identifying rare ligands from combinatorial nucleic acid libraries of very high complexity. It allows the selection of nucleic acid molecules with desired functions and it has been instrumental in the identification of a number of synthetic DNA and RNA molecules, so-called aptamers that recognise ligands of different chemical origin. The method is fast, it does not require special equipment and the selected aptamers typically bind their target with high affinity and high specificity. Here we summarise the recent examples of the SELEX technique within the context of identifying high-affinity ligands against parasite target molecules.     

The proto-Nucleic Acid Builder: a software tool for constructing nucleic acid analogs

The helical structures of DNA and RNA were originally revealed by experimental data. Likewise, the development of programs for modeling these natural polymers was guided by known structures. These nucleic acid polymers represent only two members of a potentially vast class of polymers with similar structural features, but that differ from DNA and RNA in the backbone or nucleobases. Xeno nucleic acids (XNAs) incorporate alternative backbones that affect the conformational, chemical, and thermodynamic properties of XNAs. Given the vast chemical space of possible XNAs, computational modeling of alternative nucleic acids can accelerate the search for plausible nucleic acid analogs and guide their rational design. Additionally, a tool for the modeling of nucleic acids could help reveal what nucleic acid polymers may have existed before RNA in the early evolution of life. To aid the development of novel XNA polymers and the search for possible pre-RNA candidates, this article presents the proto-Nucleic Acid Builder (https://github.com/GT-NucleicAcids/pnab), an open-source program for modeling nucleic acid analogs with alternative backbones and nucleobases. The torsion-driven conformation search procedure implemented here predicts structures with good accuracy compared to experimental structures, and correctly demonstrates the correlation between the helical structure and the backbone conformation in DNA and RNA.     

Modeling nucleic acids

Nucleic acids are an important class of biological macromolecules that carry out a variety of cellular roles. For many functions, naturally occurring DNA and RNA molecules need to fold into precise three-dimensional structures. Due to their self-assembling characteristics, nucleic acids have also been widely studied in the field of nanotechnology, and a diverse range of intricate three-dimensional nanostructures have been designed and synthesized. Different physical terms such as base-pairing and stacking interactions, tertiary contacts, electrostatic interactions and entropy all affect nucleic acid folding and structure. Here we review general computational approaches developed to model nucleic acid systems. We focus on four key areas of nucleic acid modeling: molecular representation, potential energy function, degrees of freedom and sampling algorithm. Appropriate choices in each of these key areas in nucleic acid modeling can effectively combine to aid interpretation of experimental data and facilitate prediction of nucleic acid structure.     

Evaluation of single-stranded nucleic acids as carriers in the DNA-directed assembly of macromolecules

Current developments in nanosciences indicate that the self-assembly of macromolecules, such as proteins or metallic nanoclusters, can be conveniently achieved by means of nucleic acid hybridization. Within this context, we here report on the evaluation of single-stranded nucleic acids to be utilized as carrier backbones in DNA-directed self-assembly. A microplate solid-phase hybridization assay is described which allows rapid experimental determination of the hybridization efficiencies of various sequence stretches within a given nucleic acid carrier strand. As demonstrated for two DNA fragments of different sequence, the binding efficiencies of several oligonucleotides depend on the formation of specific secondary structure elements within the carrier molecule. A correlation of sequence-specific hybridization capability with modeled secondary structure is also obvious from experiments using the fluorescence gel-shift analysis. Electrophoretic studies on the employment of helper oligonucleotides in the formation of supramolecular conjugates of several oligonucleotide-tagged proteins indicate, that structural constraints can be minimized by disruption of intramolecular secondary structures of the carrier molecule. To estimate the influences of the chemical nature of the carrier, gel-shift experiments are carried out to compare a 170mer RNA molecule with its DNA analogue. Ternary aggregates, containing two protein components bound to the carrier, are formed with a greater efficiency on the DNA instead of the RNA carrier backbone.     

DNA Ligases: Progress and Prospects

DNA ligases seal 5′-PO₄ and 3′-OH polynucleotide ends via three nucleotidyl transfer steps involving ligase-adenylate and DNA-adenylate intermediates. DNA ligases are essential guardians of genomic integrity, and ligase dysfunction underlies human genetic disease syndromes. Crystal structures of DNA ligases bound to nucleotide and nucleic acid substrates have illuminated how ligase reaction chemistry is catalyzed, how ligases recognize damaged DNA ends, and how protein domain movements and active-site remodeling are used to choreograph the end-joining pathway. Although a shared feature of DNA ligases is their envelopment of the nicked duplex as a C-shaped protein clamp, they accomplish this feat by using remarkably different accessory structural modules and domain topologies. As structural, biochemical, and phylogenetic insights coalesce, we can expect advances on several fronts, including (i) pharmacological targeting of ligases for antibacterial and anticancer therapies and (ii) the discovery and design of new strand-sealing enzymes with unique substrate specificities.     

Structural insights into regulation and action of SWI2/SNF2 ATPases

This review focuses on recent structural insights into regulation and nucleic acid binding of Superfamily 2 (SF2)-type helicases as they relate to chromatin remodelers. We review structural features of the Chd1 chromatin remodeler regarding regulation of the ATPase motor, and discuss related strategies observed for other SF2 ATPases. Since no SWI2/SNF2 ATPases have yet been captured bound to DNA in a state competent for ATP hydrolysis, we turn to structural examples from the DEAD-box RNA helicase family, and suggest that SWI2/SNF2-specific inserts may be poised to alter canonical duplex DNA structure.     

Difference in conformational diversity between nucleic acids with a six-membered 'sugar' unit and natural 'furanose' nucleic acids

Natural nucleic acids duplexes formed by Watson-Crick base pairing fold into right-handed helices that are classified in two families of secondary structures, i.e. the A- and B-form. For a long time, these A and B allomorphic nucleic acids have been considered as the 'non plus ultra' of double-stranded nucleic acids geometries with the only exception of Z-DNA, a left-handed helix that can be adopted by some DNA sequences. The five-membered furanose ring in the sugar-phosphate backbone of DNA and RNA is the underlying cause of this restriction in conformational diversity. A collection of new Watson-Crick duplexes have joined the 'original' nucleic acid double helixes at the moment the furanose sugar was replaced by different types of six-membered ring systems. The increase in this structural and conformational diversity originates from the rigid chair conformation of a saturated six-membered ring that determines the orientation of the ring substituents with respect to each other. The original A- and B-form oligonucleotide duplexes have expanded into a whole family of new structures with the potential for selective cross-communication in a parallel or antiparallel orientation, opening up a new world for information storage and for molecular recognition-directed self-organization.     

Recent advances in the in vitro evolution of nucleic acids

Molecular evolution allows chemists and biologists to generate nucleic acids with tailor-made binding or catalytic activities. Recent examples of nucleic acid evolution in vitro provide insights into natural ribozyme evolution and also demonstrate potential applications of evolved DNA and RNA molecules. Efforts to expand the scope of nucleic acid evolution are also underway, including the development of novel methods for exploring nucleic acid sequence-space and the incorporation of non-natural chemical functionality into nucleic acid libraries.