Computational evaluation of intermolecular interactions of a universal base 3-nitropyrrole in stacked dimers and DNA duplexes

The stacking interactions between a universal base of 3-nitropyrrole (3NP) and four canonical nucleobases were studied by means of ab initio molecular orbital calculations. The stabilities of the complexes are comparable to those of the stacked dimers of canonical bases reported previously. The detailed analysis of the interaction energies revealed the importance of the dipole-dipole interaction included in the Hartree-Fock terms to determine the geometry dependence of the stacking energies. It was also clarified that the dispersion energies included in the electron-correlation terms were essential to obtain adequate stabilities. The contribution of the nitro group was evaluated by the comparative studies of pyrrole and 3NP. The increased molecular dipole moment and surface are expected to account for the enhancement of the stability of the stacked dimers containing 3NP. The force field parameters required for calculation of the molecular mechanics of 3NP were obtained for 3NP on the basis of these molecular orbital calculations. The energy-minimized structures obtained by the molecular mechanics calculations of 3NP accorded with those obtained by the molecular orbital calculations described above. A DNA duplex structure containing 3NP-A, 3NP-T, or 3NP-C was calculated by use of these force field parameters. In the case of 3NP-A, the computationally calculated structure was in good agreement with that previously determined by use of (1)H-NMR except for the orientation of the nitro group.     

Potential of mean force calculations of the stacking-unstacking process in single-stranded deoxyribodinucleoside monophosphates.

The free energy of the stacking-unstacking process of deoxyribodinucleoside monophosphates in aqueous solution has been investigated by potential of mean force calculations along a reaction coordinate, defined by the distance between the glycosidic nitrogen atoms of the bases. The stacking-unstacking process of a ribodinucleoside monophosphate was observed to be well characterized by this coordinate, which has the advantage that it allows for a dynamical backbone and flexible bases. All 16 naturally occurring DNA dimers composed of the adenine, cytosine, guanine, or thymine bases in both the 5' and the 3' positions were studied. From the free-energy profiles we observed the deepest minima for the stacked states of the purine-purine dimers, but good stacking was also observed for the purine-pyrimidine and pyrimidine-purine dimers. Substantial stacking ability was found for the dimers composed of a thymine base and a purine base and also for the deoxythymidylyl-3',5'-deoxythymidine dimer. Very poor stacking was observed for the dCpdC dimer. Conformational properties and solvent accessibility are discussed for the stacked and unstacked dimers. The potential of mean force profiles of the stacking-unstacking process for the DNA dimers are compared with the RNA dimers.     

The physical basis of nucleic acid base stacking in water

It has been argued that the stacking of adenyl groups in water must be driven primarily by electrostatic interactions, based upon NMR data showing stacking for two adenyl groups joined by a 3-atom linker but not for two naphthyl groups joined by the same linker. In contrast, theoretical work has suggested that adenine stacking is driven primarily by nonelectrostatic forces, and that electrostatic interactions actually produce a net repulsion between adenines stacking in water. The present study provides evidence that the experimental data for the 3-atom-linked bis-adenyl and bis-naphthyl compounds are consistent with the theory indicating that nonelectrostatic interactions drive adenine stacking. First, a theoretical conformational analysis is found to reproduce the observed ranking of the stacking tendencies of the compounds studied experimentally. A geometric analysis identifies two possible reasons, other than stronger electrostatic interactions, why the 3-atom-linked bis-adenyl compounds should stack more than the bis-naphthyl compounds. First, stacked naphthyl groups tend to lie further apart than stacked adenyl groups, based upon both quantum calculations and crystal structures. This may prevent the bis-naphthyl compound from stacking as extensively as the bis-adenyl compound. Second, geometric analysis shows that more stacked conformations are sterically accessible to the bis-adenyl compound than to the bis-naphthyl compound because the linker is attached to the sides of the adenyl groups, but to the ends of the naphthyl groups. Finally, ab initio quantum mechanics calculations and energy decompositions for relevant conformations of adenine and naphthalene dimers support the view that stacking in these compounds is driven primarily by nonelectrostatic interactions. The present analysis illustrates the importance of considering all aspects of a molecular system when interpreting experimental data, and the value of computer models as an adjunct to chemical intuition.     

Energy gaps of graphene clusters: the first-principles calculations based on high-throughput screening

Abstract

With the enumeration of the triangular lattice fragments, we have systematically investigated the graphene clusters (C_nH_m, n = 14 – 24) with various sizes and shapes, whose structural stabilities and electronic properties are studied by the Hückel molecular orbital (HMO) method and the first-principles calculation. According to the formation energies, we show the structural stabilities of the clusters are closely related to the shape and size, as well as the chemical potential of hydrogen. The energy gaps obtained from the HMO method are in the same trend with the ones calculated by the first-principles calculations, indicating the effective screening of the gap minimum and maximum in a fast speed. There is a general decreasing of the energy gaps with the size increment due to the quantum confinement, meanwhile, the gaps are also highly dependent on the shape of the clusters for those with the same number of carbon atom.     

Conformational Studies on Nucleosides with Furanose Ring Modifications. 1.

Abstract

Conformational energy calculations have been presented on adenine and thymine nucleosides in which the furanose ring is replaced by 2′,3′-dideoxy-2′,3′-didehydrofuran using molecular mechanics and conformational analysis. Conformational energies have been evaluated using the MM2 and AMBER94 force field parameters at two different dielectric constants. The results are presented in terms of isoenergy contours in the conformational space of the glycosidic (χ) and C4′-C5′ (γ) bonds torsions. In general, the χ-γ interrelationships exhibit similarities with the corresponding plots for unmodified nucleosides and nucleotides, reported previously. Consistency of the calculated preferred conformations with the X-ray data is sensitive to the force field employed.     

Structural and energetic heterogeneities of canonical and oxidized central guanine triad of B-DNA telomeric fragments

The intermolecular interaction energies in central guanine triad of telomeric B-DNA were estimated based on ab initio quantum chemistry calculations on the MP2/aDZ level of theory. The source of structural information was molecular dynamics simulation of both canonical (AGGGTT) and oxidized (AG8oxoGGTT) telomere units. Our calculations demonstrate that significant stiffness of central triad occurs if 8oxoG is present. The origin of such feature is mainly due to the increase of stacking interactions of 8oxoG with neighbouring guanine molecules and stronger hydrogen bonding formation of 8oxoG with cytosine if compared with canonical guanine. Another interesting observation is the context independence of stacking interactions of 8oxoG. Unlike to 5′-G₂/G₃-3′ and 5′-G₃/G₄-3′ sequences which are energetically different, 5′-G₂/8oxoG₃-3′ and 5′-8oxoG₃/G₄-3′ sequences are almost iso-energetic.     

Structure-Property Relationships in Liquid Crystals: Can Modelling do Better than Empiricism?

Abstract

This paper reviews the relationships observed experimentally between the physical properties of liquid crystals and the molecular structures of the constituent molecules, and reports molecular mechanics calculations designed to provide a predictive and interpretative basis for structure/property relationships in liquid crystals. The calculations are of the minimum energy configurations of dimers of interacting liquid crystal molecules, and the geometry of the dimers, relative orientation of molecular dipoles and the extent of parallel correlation are related to liquid crystal properties. A large number of structural types are examined and the results discussed in terms of dipole correlation, apolar angular correlation, chiral twist sense and transition temperatures of liquid crystals.     

Comparative theoretical study of the structures and stabilities of four typical gadolinium carboxylates in different scintillator solvents

The structural properties and stabilities of four typical gadolinium carboxylates (Gd-CBX) in toluene, linear alkyl benzene (LAB), and phenyl xylyl ethane (PXE) solvents were theoretically studied using density functional theory (DFT/B3LYP with the basis sets 6-311G(d) and MWB54) and the polarizable continuum model (PCM). The average Gd–ligand interaction energies (E _int, corrected for dispersion) and the values of the energy gap between the highest occupied molecular orbital and lowest unoccupied molecular orbital (\( \varDelta \) _HL) for the gadolinium complexes were calculated to compare the relative stabilities of the four Gd-CBX molecules in the three liquid scintillator solvents. According to the calculations, the values of E _int and \( \varDelta \) _HL for Gd-CBX in LAB are larger than the corresponding values in PXE and toluene. Gd-CBX may therefore be more compatible with LAB than with PXE and toluene. It was also found that, in the three scintillator solvents, the stabilities of the four Gd-CBX molecules increase in the order Gd-2EHA?<?Gd-2MVA?<?Gd-pivalate?<?Gd-TMHA.     

Nature of the Stacking of Nucleic Acid Bases in Water: a Monte Carlo Simulation

Abstract

The results of a Monte Carlo simulation of the hydration of uracil and thymine molecules, their stacked dimers and hydrogen-bonded base pairs are presented. Simulations have been performed in a cluster approximation. The semiempirical atom-atom potential functions have been used (cluster consisting of 200 water molecules). It has been shown that the stacking interactions of uracil and thymine molecules in water arise mainly due to the increase in the water-water interaction during the transition from monomers to dimer. It has been found out that stacked base associates are more preferable than base pairs in water. This preference is mainly due to the energetically more favourable structure of water around the stack.     

An ab initio molecular orbital study on the sequence-dependency of DNA conformation: an evaluation of intra- and inter-strand stacking interaction energy

The various nearest neighbor stacking interaction energies of stacked base pairs in the DNA double helix are calculated for both A- and B-type conformations using an ab initio molecular orbital method. It is demonstrated that the sequence-dependent conformational preference for A- or B-type results from the stacking interaction. In particular, the base sequence showing the highest preference for an A-type conformation is revealed as GC/GC, and the one with the next highest preference, AT/AT; for a B-type conformation, the respective sequences are CG/CG and CA/TG. The overall conformation of a DNA fragment is not determined by these particular sequences only but is influenced by all base pair steps. An intrinsically favorable conformation is predicted from the constituent stacking interaction.     

Simulation Studies on the Stabilities of Aggregates Formed by Fibril-Forming Segments of α-Synuclein

Abstract

We performed molecular dynamics simulations for various oligomers with different β-sheet conformations consisting of α-Synuclein 71–82 residues using an all atom force field and explicit water model. Tetramers of antiparallel β-sheet are shown to be stable, whereas parallel sheets are highly unstable due to the repulsive interactions between bulky and polar side chains as well as the weaker backbone hydrogen bonds. We also investigated the stabilities of double antiparallel β-sheets stacked with asymmetric and symmetric geometries. Our results show that this 12 amino acid residue peptide can form stable β-sheet conformers at 320K and higher temperatures. The backbone hydrogen bonds in β-sheet and the steric packing between hydrophobic side chains between β-sheets are shown to give conformational stabilities.     

Synthesis and properties of ApU analogues containing 2'-halo-2'-deoxyadenosines. Effects of 2' substituents on oligonucleotide conformation

Five A-U analogues containing deoxyadenosine or 2'-halo-2'-deoxyadenosines, which are known to have widely different C3'-endo conformer populations according to their electronegativities of the halogen substituents, dAfl-U, dAcl-U, dAbr-U, dAio-U, and dA-U, were synthesized chemically. Characterization of these dimers has been performed by UV absorption, circular dichroism, and proton nuclear magnetic resonance spectroscopy. The results show that the dimers containing 2'-halo-2'-deoxyadenosines have stacked conformations with a geometry similar to that of A-U and the degree of stacking decreases in the order dAfl-U greater than dAcl-U greater than dAbr-U greater than dAio-U. dAcl-U is assumed to have the same degree of stacking as A-U. dA-U takes a more stacked conformation than does dAio-U, but the mode of stacking is different from those of the other dimers. The effects of the 2' substituents on dimer conformation are discussed in terms of electronegativity, molecular size, and hydrophobicity.     

Molecular design of aminopolynitroazole-based high-energy materials

The density functional theory (DFT) was employed to calculate the energetic properties of several aminopolynitroazoles. The calculations were performed to study the effect of amino and nitro substituents on the heats of formation, densities, detonation performances, thermal stabilities, and sensitivity characteristics of azoles. DFT-B3LYP, DFT-B3PW91, and MP2 methods utilizing the basis sets 6-31 G* and 6-311 G (2df, 3p) were adopted to predict HOFs via designed isodesmic reactions. All of the designed aminopolynitroazoles had heats of formation of >220 kJ mol(-1). The crystal densities of the aminopolynitroazoles were predicted with the cvff force field. All of the energetic azoles had densities of >1.83 g/cm(3). The detonation velocities and pressures were evaluated using the Kamlet-Jacobs equations, utilizing the predicted densities and heats of formation. It was found that aminopolynitroazoles have a detonation velocity of about 9.1 km/s and detonation pressure of 36 GPa. The bond dissociation energies for the C-NO(2) and N-NO(2) bonds were analyzed to investigate the stabilities of the designed molecules. The charge on the nitro group was used to assess impact sensitivity in the present study. The results obtained imply that the designed molecules are stable and are expected to be candidates for high-energy materials (HEMs).     

Vibrational Study of a Nucleoside Analogue with Antitumoral and Antiviral Activity, 5-Fluoro-2′-deoxyuridine,FdU

Abstract

The advantages of different methods for obtaining a reliable assignment of the vibrational spectra of the antitumoral and antiviral nucleoside analogue, 5-fluoro-2′-deoxyuridine. FdU, are evaluated as a basis for the study of FdU-containing polymers and drug-target interactions. The experimental FT-IR and FT-Raman spectra, are compared with theoretical frequencies obtained by a classical mechanics method and a semiempirical molecular orbital (MO) method, PM3.     

Hydrogen Bonding and Stacking of DNA Bases: A Review of Quantum-chemical ab initio Studies

Abstract

Ab initio quantum-chemical calculations with inclusion of electron correlation made since 1994 (such reliable calculations were not feasible before) significantly modified our view on interactions of nucleic acid bases. These calculations allowed to perform the first reliable comparison of the strength of stacked and hydrogen bonded pairs of nucleic acid bases, and to characterize the nature of the base-base interactions. Although hydrogen-bonded complexes of nucleobases are primarily stabilized by the electrostatic interaction, the dispersion attraction is also important. The stacked pairs are stabilized by dispersion attraction, however, the mutual orientation of stacked bases is determined rather by the electrostatic energy. Some popular theories of stacking were ruled out: The theory based on attractive interactions of polar exocyclic groups of bases with delocalized electrons of the aromatic rings (Bugg et al., Biopolymers 10, 175 (1971).), and the II-II interactions model (C.A. Hunter, J. Mol. Biol. 230, 1025 (1993)). The calculations demonstrated that amino groups of nucleobases are very flexible and intrinsically nonplanar, allowing hydrogen-bond-like interactions which are oriented out of the plane of the nucleobase. Many H-bonded DNA base pairs are intrinsically nonplanar. Higher-level ab initio calculations provide a unique set of reliable and consistent data for parametrization and verification of empirical potentials. In this article, we present a short survey of the recent calculations, and discuss their significance and limitations. This summary is written for readers which are not experts in computational quantum chemistry.