Hydrogen bonding of nucleotide base pairs: Application of the PM3 method

The ability of the pm3 semiempirical quantum mechanical method to reproduce hydrogen bonding in nucleotide base pairs was assessed. Results of pm3 calculations on the nucleotides 2′-deoxyadenosine 5′-monophosphate (pdA), 2′-deoxyguanosine 5′-monophosphate (pdG), 2′-deoxycytidine 5′-monophosphate (pdC), and 2′-deoxythymidine 5′-monophosphate (pdT) and the base pairs pdA–pdT, pdG–pdC, and pdG(syn)–pdC are presented and discussed. The pm3 method is the first of the parameterized nddo quantum mechanical models with any ability to reproduce hydrogen bonding between nucleotide base pairs. Intermolecular hydrogen bond lengths between nucleotides displaying Watson–Crick base pairing are 0.1–0.2 Å less than experimental results. Nucleotide bond distances, bond angles, and torsion angles about the glycosyl bond (χ), the C4′C5′ bond (γ), and the C5′O5′ bond (β) agree with experimental results. There are many possible conformations of nucleotides. pm3 calculations reveal that many of the most stable conformations are stabilized by intramolecular CHO hydrogen bonds. These interactions disrupt the usual sugar puckering. The stacking interactions of a dT–pdA duplex are examined at different levels of gradient optimization. The intramolecular hydrogen bonds found in the nucleotide base pairs disappear in the duplex, as a result of the additional constraints on the phosphate group when part of a DNA backbone. Sugar puckering is reproduced by the pm3 method for the four bases in the dT–pdA duplex. pm3 underestimates the attractive stacking interactions of base pairs in a B-DNA helical conformation. The performance of the pm3 method implemented in SPARTAN is contrasted with that implemented in MOPAC. At present, accurate ab initio calculations are too timeconsuming to be of practical use, and molecular mechanics methods cannot be used to determine quantum mechanical properties such as reaction-path calculations, transition-state structures, and activation energies. The pm3 method should be used with extreme caution for examination of small DNA systems. Future parameterizations of semiempirical methods should incorporate base stacking interactions into the parameterization data set to enhance the ability of these methods.

Further Quantum Mechanical Evidence that Difluorotoluene does not Hydrogen Bond

Calculations were run on the methylated DNA base pairs adenine:thymine and adenine:difluorotoluene to further investigate the hydrogen-bonding properties of difluorotoluene (F). Geometries were optimized using hybrid density functional theory. Single-point calculations at the MP2(full) level were performed to obtain more rigorous energies. The functional counterpoise method was used to correct for the basis set superposition error (BSSE), and the interaction energies were also corrected for fragment relaxation. These corrections brought the B3LYP and MP2 interaction energies into excellent agreement. In the gas phase, the Gibbs free energies calculated at the B3LYP and MP2 levels of theory predict that A and T will spontaneously form an A:T pair while A:F spontaneously dissociates into A and F. Solvation effects on the pairing of the bases were explored using implicit solvent models for water and chloroform. In aqueous solution, both A:T and A:F are predicted to dissociate into their component monomers. Semiempirical calculations were performed on small sections of B-form DNA containing the two pairs, and the results provide support for the concept that base stacking is more important than hydrogen bonding for the stability of the A:F pair within a DNA helix.

Quantum Mechnical Investion of Cyclic 3',5'- Adenosine Monophosphate, the Second Hormonal Messenger

Full geometry optimizations using the PM3, AM1, 3-21G∗/HF and 6-31G∗/HF levels of theory were conducted on the syn and anti conformations of cyclic3′,5′-adenosine monophosphate (cAMP). Comparison of the anti crystal structures with the semiempirical and ab initio results revealed that the ab initio results agree well with the experimental results. The results of semiempirical calculations are in qualitative agreement with experimental and ab initio values, with the exception of the glycosyl torsion angle for the anti conformer. Sugar puckering, which is not handled properly by semiempirical methods for unconstrained sugars, nucleosides, nucleotides and nucleotide base pairs, is modeled reasonably well by the semiempirical methods for cAMP. This improvement results from the constraints introduced by the cyclization of AMP to form the phosphodiester.

Crystal structure of a CAP-DNA complex: the DNA is bent by 90 degrees

The 3 angstrom resolution crystal structure of the Escherichia coli catabolite gene activator protein (CAP) complexed with a 30-base pair DNA sequence shows that the DNA is bent by 900. This bend results almost entirely from two 400 kinks that occur between TG/CA base pairs at positions 5 and 6 on each side of the dyad axis of the complex. DNA sequence discrimination by CAP derives both from sequence-dependent distortion of the DNA helix and from direct hydrogen-bonding interactions between three protein side chains and the exposed edges of three base pairs in the major groove of the DNA. The structure of this transcription factor-DNA complex provides insights into possible mechanisms of transcription activation

Crystallization of Escherichia coli Catabolite Gene Activator Protein with its DNA Binding Site: The use of Modular DNA

To obtain crystals of the Escherichia coli catabolite gene activator protein (CAP) complexed with its DNA-binding site, we have searched for crystallization conditions with 26 different DNA segments ≥28 base-pairs in length that explore a variety of nucleotide sequences, lengths, and extended 5′ or 3′ termini. In addition to utilizing uninterrupted asymmetric lac site sequences, we devised a novel approach of synthesizing half-sites that allowed us to efficiently generate symmetric DNA segments with a wide variety of extended termini and lengths in the large size range (≥28 bp) required by this protein. We report three crystal forms that are suitable for X-ray analysis, one of which (crystal form III) gives measurable diffraction amplitudes to 3 Å resolution. Additives such as calcium, n-octyl-β-d-glucopyranoside and spermine produce modest improvements in the quality of diffraction from crystal form III. Adequate stabilization of crystal form III is unexpectedly complex, requiring a greater than tenfold reduction in the salt concentration followed by addition of 2-methyl-2,4-pentanediol and then an increase in the concentration of polyethylene glycol.

Utilization of Probabilistic Models in Short Read Assembly from Second-Generation Sequencing

With the advent of cheaper and faster DNA sequencing technologies, assembly methods have greatly changed. Instead of outputting reads that are thousands of base pairs long, new sequencers parallelize the task by producing read lengths between 35 and 400 base pairs. Reconstructing an organism’s genome from these millions of reads is a computationally expensive task. Our algorithm solves this problem by organizing and indexing the reads using n-grams, which are short, fixed-length DNA sequences of length n. These n-grams are used to efficiently locate putative read joins, thereby eliminating the need to perform an exhaustive search over all possible read pairs. Our goal was develop a novel n-gram method for the assembly of genomes from next-generation sequencers. Specifically, a probabilistic, iterative approach was utilized to determine the most likely reads to join through development of a new metric that models the probability of any two arbitrary reads being joined together. Tests were run using simulated short read data based on randomly created genomes ranging in lengths from 10,000 to 100,000 nucleotides with 16 to 20x coverage. We were able to successfully re-assemble entire genomes up to 100,000 nucleotides in length.