Direct Detection of Unnatural DNA Nucleotides dNaM and d5SICS using the MspA Nanopore

Malyshev et al. showed that the four-letter genetic code within a living organism could be expanded to include the unnatural DNA bases dNaM and d5SICS. However, verification and detection of these unnatural bases in DNA requires new sequencing techniques. Here we provide proof of concept detection of dNaM and d5SICS in DNA oligomers via nanopore sequencing using the nanopore MspA. We find that both phi29 DNA polymerase and Hel308 helicase are capable of controlling the motion of DNA containing dNaM and d5SICS through the pore and that single reads are sufficient to detect the presence and location of dNaM and d5SICS within single molecules.     

Exploring potentially alternative non-canonical DNA duplex structures through simulation

Hopkins proposed an alternative and chirally distinct family of double-stranded DNA (dsDNA) models that have antiparallel chains with 5′→3′ senses opposite to those of the right-handed Watson–Crick (WC) family. Termed configuration II, this family of dsDNA models contains both right-handed (II-R) and left-handed (II-L) forms, with Z-DNA as an example of the latter. Relative interstrand binding energies for six DNA duplex models, two each of configuration I-R (standard WC canonical B-DNA), II-R, and II-L for the duplex d(CGCGAATTCGCG), have been estimated under identical conditions using MM-PBSA analysis from molecular dynamics trajectories using three different AMBER force fields. These simulations support the stereo chemical soundness of configuration II dsDNA forms. Recent force fields (Barcelona Supercomputing Center [BSC] [bsc1] and Olomouc 2015 [OL15]) successfully render stable II-L structures, whereas the previous force field, bsc0, generated stable II-R structures, although with an energy difference between II-R and II-L of ～30?kcal/mol.

Communicated by Ramaswamy H. Sarma     

Systematic exploration of a class of hydrophobic unnatural base pairs yields multiple new candidates for the expansion of the genetic alphabet

We have developed a family of unnatural base pairs (UBPs), which rely on hydrophobic and packing interactions for pairing and which are well replicated and transcribed. While the pair formed between d5SICS and dNaM (d5SICS-dNaM) has received the most attention, and has been used to expand the genetic alphabet of a living organism, recent efforts have identified dTPT3-dNaM, which is replicated with even higher fidelity. These efforts also resulted in more UBPs than could be independently analyzed, and thus we now report a PCR-based screen to identify the most promising. While we found that dTPT3-dNaM is generally the most promising UBP, we identified several others that are replicated nearly as well and significantly better than d5SICS-dNaM, and are thus viable candidates for the expansion of the genetic alphabet of a living organism. Moreover, the results suggest that continued optimization should be possible, and that the putatively essential hydrogen-bond acceptor at the position ortho to the glycosidic linkage may not be required. These results clearly demonstrate the generality of hydrophobic forces for the control of base pairing within DNA, provide a wealth of new structure–activity relationship data and importantly identify multiple new candidates for in vivo evaluation and further optimization.     

Conformational features of benzo‐homologated yDNA duplexes by molecular dynamics simulation

yDNA is a base‐modified nucleic acid duplex containing size‐expanded nucleobases. Base‐modified nucleic acids could expand the genetic alphabet and thereby enhance the functional potential of DNA. Unrestrained 100 ns MD simulations were performed in explicit solvent on the yDNA NMR sequence [5′(yA T yA yA T yA T T yA T)₂] and two modeled yDNA duplexes, [5′(yC yC G yC yC G G yC G G)₂] and [(yT5′ G yT A yC yG C yA yG T3′)?(yA5′ C T C yG C G yT A yC A3′)]. The force field parameters for the yDNA bases were derived in consistent with the well‐established AMBER force field. Our results show that DNA backbone can withstand the stretched size of the bases retaining the Watson‐Crick base pairing in the duplexes. The duplexes retained their double helical structure throughout the simulations accommodating the strain due to expanded bases in the backbone torsion angles, sugar pucker and helical parameters. The effect of the benzo‐expansion is clearly reflected in the extended C1′‐C1′ distances and enlarged groove widths. The size expanded base modification leads to reduction in base pair twist resulting in larger overlapping area between the stacked bases, enhancing inter and intra strand stacking interactions in yDNA in comparison with BDNA. This geometry could favour enhanced interactions with the groove binders and DNA binding proteins., 2016. © 2015 Wiley Periodicals, Inc. Biopolymers 105: 55–64, 2016     

Trinucleotide’s quadruplet symmetries and natural symmetry law of DNA creation ensuing Chargaff’s second parity rule

For almost 50 years the conclusive explanation of Chargaff’s second parity rule (CSPR), the equality of frequencies of nucleotides A=T and C=G or the equality of direct and reverse complement trinucleotides in the same DNA strand, has not been determined yet. Here, we relate CSPR to the interstrand mirror symmetry in 20 symbolic quadruplets of trinucleotides (direct, reverse complement, complement, and reverse) mapped to double-stranded genome. The symmetries of Q-box corresponding to quadruplets can be obtained as a consequence of Watson–Crick base pairing and CSPR together. Alternatively, assuming Natural symmetry law for DNA creation that each trinucleotide in one strand of DNA must simultaneously appear also in the opposite strand automatically leads to Q-box direct-reverse mirror symmetry which in conjunction with Watson–Crick base pairing generates CSPR. We demonstrate quadruplet’s symmetries in chromosomes of wide range of organisms, from Escherichia coli to Neanderthal and human genomes, introducing novel quadruplet-frequency histograms and 3D-diagrams with combined interstrand frequencies. These “landscapes” are mutually similar in all mammals, including extinct Neanderthals, and somewhat different in most of older species. In human chromosomes 1–12, and X, Y the “landscapes” are almost identical and slightly different in the remaining smaller and telocentric chromosomes. Quadruplet frequencies could provide a new robust tool for characterization and classification of genomes and their evolutionary trajectories.     

FF12MC: A revised AMBER forcefield and new protein simulation protocol

Specialized to simulate proteins in molecular dynamics (MD) simulations with explicit solvation, FF12MC is a combination of a new protein simulation protocol employing uniformly reduced atomic masses by tenfold and a revised AMBER forcefield FF99 with (i) shortened C? H bonds, (ii) removal of torsions involving a nonperipheral sp³ atom, and (iii) reduced 1–4 interaction scaling factors of torsions ? and ψ. This article reports that in multiple, distinct, independent, unrestricted, unbiased, isobaric–isothermal, and classical MD simulations FF12MC can (i) simulate the experimentally observed flipping between left‐ and right‐handed configurations for C14–C38 of BPTI in solution, (ii) autonomously fold chignolin, CLN025, and Trp‐cage with folding times that agree with the experimental values, (iii) simulate subsequent unfolding and refolding of these miniproteins, and (iv) achieve a robust Z score of 1.33 for refining protein models TMR01, TMR04, and TMR07. By comparison, the latest general‐purpose AMBER forcefield FF14SB locks the C14–C38 bond to the right‐handed configuration in solution under the same protein simulation conditions. Statistical survival analysis shows that FF12MC folds chignolin and CLN025 in isobaric–isothermal MD simulations 2–4 times faster than FF14SB under the same protein simulation conditions. These results suggest that FF12MC may be used for protein simulations to study kinetics and thermodynamics of miniprotein folding as well as protein structure and dynamics. Proteins 2016; 84:1490–1516. © 2016 The Authors Proteins: Structure, Function, and Bioinformatics Published by Wiley Periodicals, Inc.     

Structural basis of error‐prone replication and stalling at a thymine base by human DNA polymerase ι

Human DNA polymerase ι (polι) is a unique member of Y‐family polymerases, which preferentially misincorporates nucleotides opposite thymines (T) and halts replication at T bases. The structural basis of the high error rates remains elusive. We present three crystal structures of polι complexed with DNA containing a thymine base, paired with correct or incorrect incoming nucleotides. A narrowed active site supports a pyrimidine to pyrimidine mismatch and excludes Watson–Crick base pairing by polι. The template thymine remains in an anti conformation irrespective of incoming nucleotides. Incoming ddATP adopts a syn conformation with reduced base stacking, whereas incorrect dGTP and dTTP maintain anti conformations with normal base stacking. Further stabilization of dGTP by H‐bonding with Gln59 of the finger domain explains the preferential T to G mismatch. A template ‘U‐turn’ is stabilized by polι and the methyl group of the thymine template, revealing the structural basis of T stalling. Our structural and domain‐swapping experiments indicate that the finger domain is responsible for polι's high error rates on pyrimidines and determines the incorporation specificity.     

Comparison of force fields for Alzheimer's A : A case study for intrinsically disordered proteins

Intrinsically disordered proteins are essential for biological processes such as cell signalling, but are also associated to devastating diseases including Alzheimer's disease, Parkinson's disease or type II diabetes. Because of their lack of a stable three‐dimensional structure, molecular dynamics simulations are often used to obtain atomistic details that cannot be observed experimentally. The applicability of molecular dynamics simulations depends on the accuracy of the force field chosen to represent the underlying free energy surface of the system. Here, we use replica exchange molecular dynamics simulations to test five modern force fields, OPLS, AMBER99SB, AMBER99SB*ILDN, AMBER99SBILDN‐NMR and CHARMM22*, in their ability to model Aβ₄₂, an intrinsically disordered peptide associated with Alzheimer's disease, and compare our results to nuclear magnetic resonance (NMR) experimental data. We observe that all force fields except AMBER99SBILDN‐NMR successfully reproduce local NMR observables, with CHARMM22* being slightly better than the other force fields.     

Solvated and generalised Born calculations differences using GPU CUDA and multi-CPU simulations of an antifreeze protein with AMBER

While there has been an increase in the number of biomolecular computational studies employing graphics processing units (GPU), results describing their use with the molecular dynamics package AMBER with the CUDA implementation are scarce. No information is available comparing MD methodologies pmemd.cuda, pmemd.mpi or sander.mpi, available in AMBER, for generalised Born (GB) simulations or with solvated systems. As part of our current studies with antifreeze proteins (AFP), and for the previous reasons, we present details of our experience comparing performance of MD simulations at varied temperatures between multi-CPU runs using sander.mpi, pmemd.mpi and pmemd.cuda with the AFP from the fish ocean pout (1KDF). We found extremely small differences in total energies between multi-CPU and GPU CUDA implementations of AMBER12 in 1ns production simulations of the solvated system using the TIP3P water model. Additionally, GPU computations achieved typical one order of magnitude speedups when using mixed precision but were similar to CPU speeds when computing with double precision. However, we found that GB calculations were highly sensitive to the choice of initial GB parametrisation regardless of the type of methodology, with substantial differences in total energies.     

Hybrid simulation approach incorporating microscopic interaction along with rigid body degrees of freedom for stacking between base pairs

Stacking interaction between the aromatic heterocyclic bases plays an important role in the double helical structures of nucleic acids. Considering the base as rigid body, there are total of 18 degrees of freedom of a dinucleotide step. Some of these parameters show sequence preferences, indicating that the detailed atomic interactions are important in the stacking. Large variants of non‐canonical base pairs have been seen in the crystallographic structures of RNA. However, their stacking preferences are not thoroughly deciphered yet from experimental results. The current theoretical approaches use either the rigid body degrees of freedom where the atomic information are lost or computationally expensive all atom simulations. We have used a hybrid simulation approach incorporating Monte‐Carlo Metropolis sampling in the hyperspace of 18 stacking parameters where the interaction energies using AMBER‐parm99bsc0 and CHARMM‐36 force‐fields were calculated from atomic positions. We have also performed stacking energy calculations for structures from Monte‐Carlo ensemble by Dispersion corrected density functional theory. The available experimental data with Watson–Crick base pairs are compared to establish the validity of the method. Stacking interaction involving A:U and G:C base pairs with non‐canonical G:U base pairs also were calculated and showed that these structures were also sequence dependent. This approach could be useful to generate multiscale modeling of nucleic acids in terms of coarse‐grained parameters where the atomic interactions are preserved. This method would also be useful to predict structure and dynamics of different base pair steps containing non Watson–Crick base pairs, as found often in the non‐coding RNA structures. © 2015 Wiley Periodicals, Inc. Biopolymers 105: 212–226, 2016.     

Favored and less favored codon–anticodon duplexes arising from the GC codon family box encoding for alanine: some computational perspectives

Alanine is encoded by the four codons of the GC box (GCA, GCG, GCU, and GCC). Known alanine anticodons include the UGC, IGC, and VGC triplets (I = inosine; V = uridine-5-oxyacetic acid). The energy-minimized structures of all possible codon–anticodon combinations involving all the alanine codons GCA, GCG, GCU, and GCC with the alanine anticodons UGC, IGC, and VGC are studied using the AMBER software. Fifteen H-bonded duplex structures arising out of these combinations are studied here, all having Watson–Crick-type base pairs at the first and second codon positions, and a variety of base pairing possibilities at the third (or wobble) position. Structural and stability considerations suggest that some codon–anticodon duplexes would be more favored than others for accommodation during the translation process. The UGC anticodon is predicted to favor the GCA codon for reading, while the GCC codon is least favored. The IGC anticodon would prefer to read the GCC codon, the GCG codon being least favored, while a syn conformer for A in the GCA codon could allow for it to be read. For the VGC anticodon, the GCA codon is predicted to be read most favorably, and the GCC codon least favorably, while a syn conformer for V in the anticodon would allow for the codon GCU to be read through a wobble pair which involves the exocyclic 5-oxyacetate group of V in H-bonding.     

Strike‐induced chemosensory searching (SICS) in northern Pacific rattlesnakes (Crotalus viridis oreganus,Holbrook, 1840) rescued from substandard husbandry conditions. II. Complete recovery of function after two years

Six northern Pacific rattlesnakes (Crotalus viridis oreganus) rescued from substandard husbandry conditions were shipped from Woodland Park Zoological Gardens to the University of Colorado in December 1994. The snakes accepted prey and rapidly gained weight, but their searching behavior after predatory strikes was seriously depressed. We administered repeated tests of strike‐induced chemosensory searching (SICS) finding that the snakes improved steadily with each succeeding test, until performance eventually reached the level characteristic of this taxon. Although the snakes reached normal body weight within 4 months of exposure to proper husbandry conditions, normal SICS was not seen until 2 years had elapsed. Release of these snakes into natural habitat prior to recovery of SICS may well have resulted in the demise of the animals, suggesting that rehabilitation efforts must consider recovery of essential behavior patterns that can take far longer than recovery of body weight. Zoo Biol 18:141–146, 1999. © 1999 Wiley‐Liss, Inc.     

Simple molecular engineering of glycol nucleic acid: Progression from self-pairing to cross-pairing with cDNA and RNA

The acyclic chiral nucleic acid analogue, Glycol Nucleic Acid (GNA), displayed exceptional structural simplicity and atom economy while forming self-paired duplexes, using canonical Watson–Crick base pairing. We disclose here that the replacement of phosphodiester linker in GNA with somewhat rigid and shorter carbamate linker in Glycol Carbamate Nucleic Acid (GCNA) backbone allows unprecedented stability to the antiparallel self-paired duplexes. The R-GCNA oligomers were further found to form cross-paired antiparallel duplexes with cDNA and RNA following Watson–Crick base pairing. The stability of cross-paired GCNA:DNA and GCNA:RNA duplexes was higher than the corresponding DNA:DNA and DNA:RNA duplexes. The chiral (R) and (S) precursors were easily accessible from naturally occurring l-serine.     

Thermodynamics calculation of protein–ligand interactions by QM/MM polarizable charge parameters

The calculation of protein–ligand binding free energy (ΔG) is of great importance for virtual screening and drug design. Molecular dynamics (MD) simulation has been an attractive tool to investigate this scientific problem. However, the reliability of such approach is affected by many factors including electrostatic interaction calculation. Here, we present a practical protocol using quantum mechanics/molecular mechanics (QM/MM) calculations to generate polarizable QM protein charge (QMPC). The calculated QMPC of some atoms in binding pockets was obviously different from that calculated by AMBER ff03, which might significantly affect the calculated ΔG. To evaluate the effect, the MD simulations and MM/GBSA calculation with QMPC for 10 protein–ligand complexes, and the simulation results were then compared to those with the AMBER ff03 force field and experimental results. The correlation coefficient between the calculated ΔΔG using MM/GBSA under QMPC and the experimental data is .92, while that with AMBER ff03 force field is .47 for the complexes formed by streptavidin or its mutants and biotin. Moreover, the calculated ΔΔG with QMPC for the complexes formed by ERβ and five ligands is positively related to experimental result with correlation coefficient of .61, while that with AMBER ff03 charge is negatively related to experimental data with correlation coefficient of .42. The detailed analysis shows that the electrostatic polarization introduced by QMPC affects the electrostatic contribution to the binding affinity and thus, leads to better correlation with experimental data. Therefore, this approach should be useful to virtual screening and drug design.     

Molecular dynamics simulations of xDNA

xDNA is a modified DNA, which contains natural as well as expanded bases. Expanded bases are generated by the addition of a benzene spacer to the natural bases. A set of AMBER force‐field parameters were derived for the expanded bases and the structural dynamics of the xDNA decamer ( xT5 ′ G xT A xC xG C xA xG T3′ ) · ( xA5′ C T xG C G xT A xC A3′) was explored using a 22 ns molecular dynamics simulation in explicit solvent. During the simulation, the duplex retained its Watson‐Crick base‐pairing and double helical structure, with deviations from the starting B‐form geometry towards A‐form; the deviations are mainly in the backbone torsion angles and in the helical parameters. The sugar pucker of the residues were distributed among a variety of modes; C2′ endo, C1′ exo, O4′ endo, C4′ exo, C2′ exo, and C3′ endo. The enhanced stacking interactions on account of the modification in the bases could help to retain the duplex nature of the helix with minor deviations from the ideal geometry. In our simulation, the xDNA showed a reduced minor groove width and an enlarged major groove width in comparison with the NMR structure. Both the grooves are larger than that of standard B‐DNA, but major groove width is larger than that of A‐DNA with almost equal minor groove width. The enlarged groove widths and the possibility of additional hydration in the grooves makes xDNA a potential molecule for various applications. © 2009 Wiley Periodicals, Inc. Biopolymers 91: 351–360, 2009. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Molecular mechanics modeling of oligonucleotide adducts of the antitumor drug cis-diamminedichloroplatinum(II)

In order to assess the geometric changes caused when the antitumor drug cis-diammine-dichloroplatinum(II) (cis-DDP) binds to DNA, molecular mechanics calculations were performed on two double-stranded and two single-stranded oligonucleotides and their adducts with cis-{Pt(NH₃)₂}²⁺. For the platinated duplexes, three model structures have been derived, one involving only local disruption of base pairing with retention of the helix directionality, and two models showing pronounced kinking of the double helix. One of the kinked models is stabilized by bridging sodium ions. The other kinked duplex model shows retention of all Watson–Crick base pairing, including that of the coordinated guanines. All models exhibit hydrogen bonds connecting one ammine ligand of platinum with one or two phosphate groups located at the 5′ side of the platinated strand.     

NMR structure of duplex DNA containing the α‐OH‐PdG·dA base pair: A mutagenic intermediate of acrolein

Acrolein, a cell metabolic product and main component of cigarette smoke, reacts with DNA generating α‐OH‐PdG lesions, which have the ability to pair with dATP during replication thereby causing G to T transversions. We describe the solution structure of an 11‐mer DNA duplex containing the mutagenic α‐OH‐PdG·dA base pair intermediate, as determined by solution nuclear magnetic resonance (NMR) spectroscopy and retrained molecular dynamics (MD) simulations. The NMR data support a mostly regular right‐handed helix that is only perturbed at its center by the presence of the lesion. Undamaged residues of the duplex are in anti orientation, forming standard Watson‐Crick base pairs alignments. Duplication of proton signals at and near the damaged base pair reveals the presence of two enantiomeric duplexes, thus establishing the exocyclic nature of the lesion. The α‐OH‐PdG adduct assumes a syn conformation pairing to its partner dA base that is protonated at pH 6.6. The three‐dimensional structure obtained by restrained molecular dynamics simulations show hydrogen bond interactions that stabilize α‐OH‐PdG in a syn conformation and across the lesion containing base pair. We discuss the implications of the structures for the mutagenic bypass of acrolein lesions. © 2010 Wiley Periodicals, Inc. Biopolymers 93: 391–401, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Hx-amides: DNA sequence recognition by the fluorescent Hx (p-anisylbenzimidazole)•pyrrole and Hx•imidazole pairings

Hx-amides are fluorescent hybrids of imidazole (I)- and pyrrole (P)-containing polyamides and Hoechst 33258, and they bind in the minor groove of specific DNA sequences. Synthesis and DNA binding studies of HxII (5) complete our studies on the first set of Hx-amides: Hx–I/P–I/P. HxPP (2), HxIP (3) and HxPI (4) were reported earlier. Results from DNase I footprinting, biosensor-SPR, CD and ΔT_M studies showed that Hx-amides interacted with DNA via the anti-parallel and stacked, side-by-side motif. Hx was found to mimic the DNA recognition properties of two consecutive pyrrole units (PP) in polyamides. Accordingly, the stacked Hx/PP pairing binds preferentially to two consecutive AT base pairs, A/T–A/T; Hx/IP prefers C–A/T; Hx/PI prefers A/T–C; and Hx/II prefers C–C. The results also showed that Hx-amides bound their cognate sequence at a higher affinity than their formamido-triamide counterparts.     

Structural equilibrium of DNA represented with different force fields.

We have recently indicated preliminary evidence of different equilibrium average structures with the CHARMM and AMBER force fields in explicit solvent molecular dynamics simulations on the DNA duplex d(C5T5) . d(A5G5) (Feig, M. and B.M. Pettitt, 1997, Experiment vs. Force Fields: DNA conformation from molecular dynamics simulations. J. Phys. Chem. B. (101:7361-7363). This paper presents a detailed comparison of DNA structure and dynamics for both force fields from extended simulation times of 10 ns each. Average structures display an A-DNA base geometry with the CHARMM force field and a base geometry that is intermediate between A- and B-DNA with the AMBER force field. The backbone assumes B form on both strands with the AMBER force field, while the CHARMM force field produces heterogeneous structures with the purine strand in A form and the pyrimidine strand in dynamical equilibrium between A and B conformations. The results compare well with experimental data for the cytosine/guanine part but fail to fully reproduce an overall B conformation in the thymine/adenine tract expected from crystallographic data, particularly with the CHARMM force field. Fluctuations between A and B conformations are observed on the nanosecond time scale in both simulations, particularly with the AMBER force field. Different dynamical behavior during the first 4 ns indicates that convergence times of several nanoseconds are necessary to fully establish a dynamical equilibrium in all structural quantities on the time scale of the simulations presented here.     

Strike-induced chemosensory searching in old world vipers and new world pit vipers at San Diego Zoo

Rattlesnakes (Crotalus unicolor, C vegrandis) and Russell's vipers (Vipera russelli) at San Diego Zoo exhibited sustained high rates of tongue flicking after striking mice (Mus musculus) but not after seeing, smelling, and/or detecting thermal cues arising from mice. Called strike-induced chemosensory searching (SICS), this phenomenon contributes to poststrike trailing of envenomated prey. Because these zoo-raised snakes had always been offered dead rodents, and because these prey were usually ingested without first being struck, the present snakes had rarely exercised their innate predatory repertoire (ie, ambush tactics including striking, releasing, and trailing). Indeed, most specimens had never before struck a mouse and, hence, had never exhibited SICS. The occurrence of SICS in the present study clearly indicates that this important aspect of the predatory repertoire had not been degraded as a consequence of long-term captive husbandry.