UV photoelectron spectroscopy and ab initio characterization of valence orbital structures and conformations of neutral phosphate esters

The HeI UV photoelectron spectrum of trimethyl phosphate (TMP) has been measured and interpreted with the aid of SCF molecular orbital calculations carried out with STO-3G, STO-3G* and 4-31G basis functions. The photoelectron spectrum of TMP is more accurately reproduced by results from 4-31G calculations than by results from STO-3G or STO-3G* calculations. However, all three basis sets yield results which predict the same assignment of the photoelectron spectrum. Results at the 4-31G level indicate that whether calculations are based on crystallographic bond angles and bond lengths or on STO-3G optimized geometries has little effect on the energetic ordering of the upper occupied orbitals. The energetic ordering of orbitals is also found to be only weakly dependent upon the torsional angle phi, describing rotation of ester groups about P-O bonds and upon the torsional angle psi, describing rotation of methyl groups about C-O bonds. For trimethyl phosphate, with C3 symmetry, the vertical ionization potentials of the upper occupied orbitals are 10.81 eV (8e), 11.4 eV (9a), 11.93 eV (7e), 12.6-12.9 eV (8a and 6e), 14.4 eV (7a) and 15.0-16.0 eV (5e and 6a). Calculations at the 4-31G level indicate that many of the highest occupied orbitals in neutral dimethyl phosphate and methyl phosphate have energies and electron distributions similar to orbitals in TMP. For TMP, a search for optimized values of phi and psi has been carried out at the STO-3G*level. In agreement with previous NMR studies and with classical potential calculations, the STO-3G* results indicate that both the gauche (phi = 53.1 degrees) and anticlinal (phi = 141.9 degrees) conformations are thermally accessible. Also in agreement with the classical potential calculations, the STO-3G* results predict that in the all gauche conformation energy is minimized when the methyl groups assume a staggered geometry (psi = 60 degrees to 80 degrees) and that an energy maximum occurs for an eclipsed geometry (phi = 0 degrees to 20 degrees). A study of the dependence of optimized values of O-P-O ester bond angles on the torsional angles, phi, was carried out at the STO-3G, STO-3G* and 4-31G levels. The results demonstrate that for C3 symmetry, the coupling of O-P-O angles to phi is influence by repulsive steric interactions.     

Spatial configuration of single-stranded polynucleotides. Calculations of average dimensions and Nmr coupling constants

The probability distributions of the torsional angles (Φ′, ω′, ω, Φ, and ψ), which fix the structure of nucleotide backbone, have been calculated using the results of energy calculations based on extended Huckel theory (EHT), complete neglect of differential overlap (CNDO), perturbative configuration interaction using localized orbitals (PCILO), and classical potential functions (CPF) methods. Statistical average values of the vicinal ¹H? ¹H, ¹H? ³¹P, and ¹³C? ³¹P nmr coupling constants 〈J〉 have been calculated from the generalized Karplus relations using the probability distribution in the Φ′, Φ, and ψ space. Experimental 〈J〉 values for polyribouridylic acid (polyU) support the theoretical predictions for these torsional angles. Using Monte Carlo technique, random coils of single-stranded polynucleotides have been simulated and the mean-square end-to-end distance 〈r²〉 has been calculated. Molecular orbital methods (EHT, CNDO, and PCILO) suggest considerable flexibility around O? P bonds, leading to fairly small values for the characteristic ratio (C_∞ ～ 4). Observed values of the unperturbed characteristic ratio for polynucleotides are quite large (C_∞ ～ 18) suggesting a relatively rigid nucleotide backbone. The results based on molecular orbital calculations can be reconciled with the experimental values by introducing an additional stabilization of ～2 kcal mol^?1 for the predicted minimum energy ragion (Φ′ ～ 240°, ω′ ～ 290°, ω 290°, Φ 180°, and ψ 60°). Such a stabilization may arise from the association of water molecules and metal ions with the phosphate group and (or) Coulomb interaction between neighboring phosphate groups. The calculations provide a semiquantitative estimate of torsional rigidity in the nucleotide backbone.     

Ab initio Investigation of the Molecular Structure of Methyl Methoxymethyl Phosphonate,a Promising Nuclease-resistant Alternative of the Phosphodiester Linkage

Abstract

Conformational flexibility of the methyl methoxymethyl phosphonate anion (CH₃-O-PO₂- CH₂-O-CH₃)^?, a nuclease resistant alternative to the phosphodiester linkage in DNA, have been investigated by ab initio quantum mechanical calculations. The potential of backbone torsional degrees of freedom of methyl methoxymethyl phosphonate anion (MMP) was determined at the Hartree-Fock (HF) 3–21G* level using the adiabatic mapping technique. Energies, geometries, and effective atomic charges of different conformers were calculated at HF/6–31G* and MP2/6–31G* levels of theory. These were compared to the results obtained for dimethyl phosphate calculated at the same level. The impact on DNA structure from inserting a methylene group between phosphorus and oxygen of the nucleoside sugar moiety was examined via distance- and angle-constrained geometry optimizations. Due to its high flexibility, MMP has been shown to be compatible with both A and B forms of DNA.     

An experimental and theoretical approach to the molecular structure of 2-{4-[3-(2,5-dimethylphenyl)-3-methylcyclobutyl]thiazol-2-yl}isoindoline-1,3-dione

The title compound, 2-{4-[3-(2,5-dimethylphenyl)-3-methylcyclobutyl]thiazol-2-yl}isoindoline-1,3-dione (C₂₄H₂₂N₂O₂S), was synthesized and characterized by IR-NMR spectroscopy and single-crystal X-ray diffraction. The compound crystallizes in the monoclinic space group P2₁/c with a?=?19.7799(13) Å, b?=?6.7473(4) Å, c?=?15.7259(9) Å and β?=?103.416(5)°. In addition, the molecular geometry, vibrational frequencies and gauge including atomic orbital (GIAO) ¹H and ¹³C chemical shift values of the title compound in the ground state have been calculated by using the Hartree-Fock (HF) and density functional method (DFT/B3LYP) with 6–31G(d), 6–31 + G(d,p) and LANL2DZ basis sets, and compared with the experimental data. To determine conformational flexibility, molecular energy profile of the title compound was obtained by semi-empirical (AM1) calculations with respect to two selected degrees of torsional freedom, which were varied from ?180° to +180° in steps of 5°. Besides, molecular electrostatic potential, frontier molecular orbitals (FMO) analysis and thermodynamic properties of the title compound were investigated by theoretical calculations.     

Ab initio molecular orbital calculations on furanose sugars: a study with the 6-31G basis set

Ab initio molecular orbital calculations were performed on 2-deoxy-beta-D-glycero-tetrofuranose (1) using the 6-31G* basis set to evaluate the effect of ring conformation on the molecular parameters (bond lengths, angles, and torsions). Geometric optimizations were conducted on the planar and ten envelope conformers of 1, and these data were compared to those obtained from previous calculations using the STO-3G and 3-21G basis sets. Conformational energy profiles derived from 3-21G and 6-31G* data were found to be qualitatively comparable. The effect of furanose ring conformation on key bond lengths (e.g., C-H, C-O), bond angles (e.g., COC), and bond torsions (e.g., the exoanomeric C-1-O-1 torsion) was examined, and a qualitative agreement was observed between the 3-21G and 6-31G* analyses. The results indicate that, for semi-quantitative ab initio studies of intact carbohydrates, the 3-21G basis set is sufficient, and that the STO-3G basis set should not be employed unless crude structural approximations are desired. The observed concerted behavior of C-O bond lengths in the vicinity of the anomeric carbon of the aldofuranose ring has suggested a possible role of C-1-O-1 bond orientation in affecting the mechanism of glycoside bond hydrolysis.     

Methyl 2-methoxy-7-(4-methylbenzoyl)-4-oxo-6-p-tolyl-4H-furo[3,2-c]pyran-3-carboxylate:A combined experimental and theoretical investigation

The title compound, methyl 2-methoxy-7-(4-methylbenzoyl)-4-oxo-6-p-tolyl-4H-furo[3,2-c]pyran-3-carboxylate (C₂₅H₂₀O₇), was prepared and characterized by IR and single-crystal X-ray diffraction (XRD). The compound crystallizes in the triclinic space group P ?1 with a?=?8.9554(9) Å, b?=?10.0018(10) Å, c?=?12.7454(13) Å, α?=?67.678(7)°, β?=?89.359(8)° and γ?=?88.961(8)°. In addition to the molecular geometry from X-ray experiment, the molecular geometry and vibrational frequencies of the title compound in the ground state have been calculated using semiempirical AM1 and PM3 methods, as well as Hartree-Fock (HF) and density functional (B3LYP) levels of theory with 6–31G(d) basis set. To determine conformational flexibility, molecular energy profile of the title compound was obtained by semi-empirical (AM1) calculations with respect to two selected degrees of torsional freedom, which were varied from ?180° to +180° in steps of 10°. Besides, frontier molecular orbitals (FMO) analysis and thermodynamic properties of the title compound were performed by the B3LYP/6–31G(d) method.     

An experimental Ramachandran plot for retropeptide derivatives: Conformational features of derivatives of GEM-diamino and malonyl amino acids

Malonic and diaminomethane residues, equivalent to the two possible retro modifications of a glycine unit, with an inverted peptide group, present particular conformations that differ from those found in glycine and, in general, in α-amino acids. In both cases the φ_i and ψ_i torsional angles have restricted values as deduced from inspection of the Cambridge Structural Data Bank and from compounds studied by us. Thus, both ψ_i angles tend to be equal to 115° (or −115°) in the malonyl residues, whereas the φ_i angles tend to be equal to 88° (or −88°) in the diaminomethane residues. These results are in agreement with previous experimental data on polymers, but in the case of malonyl residues they differ from theoretical calculations on isolated molecules. The experimental data for both residues can be represented in a way similar to the usual Ramachandran plot, which will be useful in analyzing the incorporation of these residues into proteins. When side chains are present in either type of residue, they are similar to conventional α-amino acids, although the orientation of the peptide groups is different. In such cases they acquire conformations similar to those found in peptide residues in the α-helix and β-sheet conformations, although other conformations are also possible. © 1998 John Wiley & Sons, Inc. Biopoly 45: 149–155, 1998     

Experimental and theoretical investigation of the molecular and electronic structure of 5-(4-aminophenyl)-4-(3-methyl-3-phenylcyclobutyl)thiazol-2-amine

The title molecule, 5-(4-aminophenyl)-4-(3-methyl-3-phenylcyclobutyl)thiazol-2-amine (C₂₀H₂₁N₃S), was prepared and characterized by ¹H-NMR, ¹³C-NMR, IR and single-crystal X-ray diffraction. The compound crystallizes in the monoclinic space group P2₁/c with a?=?9.4350(5) Å, b?=?11.2796(6) Å, c?=?18.4170(8) Å and β?=?113.378(3)°. In addition to the molecular geometry from X-ray experiment, the molecular geometry, vibrational frequencies, gauge including atomic orbital (GIAO) ¹H- and ¹³C-NMR chemical shift values and atomic charges distribution of the title compound in the ground state have been calculated using the Hartree–Fock (HF) and density functional method (DFT) (B3LYP) with 6-31G(d) basis set. To determine conformational flexibility, molecular energy profile of the title compound was obtained by semi-empirical (AM1) calculations with respect to two selected degrees of torsional freedom, which were varied from ?180° to +180° in steps of 10°. Besides, frontier molecular orbitals (FMO) analysis was performed by the B3LYP/6-31G(d) method.     

Theoretical study of the geometrical and electronic structures and thermochemistry of spherophanes

A set of supramolecular cage-structures—spherophanes—was studied at the density functional B3LYP level. Full geometrical structure optimisations were made with 6–31G and 6–31G(d) basis sets followed by frequency calculations, and electronic energies were evaluated at B3LYP/6–31++G(d,p). Three different symmetries were considered: C1, Ci, and Oh. It was found that the bonds between the benzene rings are very long to allow π-electron delocalisation between them. These spherophanes show portal openings of 2.596 Å in Spher1, 4.000 Å in Meth2, 3.659 Å in Oxa3, and 4.412 Å in Thia4. From the point of view of potential host–guest interaction studies, it should also be noted that the atoms nearest to the centre of the cavities are carbons bonded to X groups. These supramolecules seem to exhibit relatively large gap HOMO?LUMO: 2.89 eV(Spher1), 5.26 eV(Meth2), 5.73 eV(Oxa3), and 4.82 eV(Thia4). The calculated ΔH°_f (298.15 K) values at B3LYP/6–31G(d) are (in kcal mol^?1) 750.98, 229.78, ?10.97, and 482.49 for Spher1, Meth2, Oxa3, and Thia4, respectively. Using homodesmotic reactions, relative to Spher1, the spherophanes Meth2, Oxa3, and Thia4 are less strained by ?399.13 kcal mol^?1, ?390.40 kcal mol^?1, and ?411.38 kcal mol^?1, respectively. Their infrared and ¹³C NMR calculated spectra are reported.     

Quantum chemical studies on the conformational structure of nucleic acids. IV. Calculation of backbone structure by CNDO method

Quantum chemical calculations using the CNDO/2 method, have been carried out to determine the energetically favoured ranges of the torsional angles (φ′, ω′, ω, φ, ψ) which fix the conformational structure of nucleic acid backbone. The two dimensional isoenergy maps have been constructed in the (ω′, ω) and (φ, ψ) hyperspaces. The variation of total energy with respect to φ′ has also been studied. The results show that the non-bonding interactions play a major role in the conformational stability of nucleic acids and polynucleotides. The theoretical predictions show good correspondence with the experimental data (X-ray and ¹³C NMR) as well as the other reported theoretical calculations (EHT, PCILO and classical potential functions). The most favoured structure has the conformational angles close to 240, 290, 290, 180 and 60° and these values lead to a helical structure with a pitch of 34 Å and about ten nucleotide units per turn of the helix. The proposed models of Watson &; Crick, DNA-B and DNA-C lie in high energy regions.     

Density functional theory study of the potassium complexation of an unsymmetrical 1,3-alternate calix[4]-crown-5-N-azacrown-5 bearing two different crown rings

Theoretical studies of an unsymmetrical calix[4]-crown-5-N-azacrown-5 (1) in a fixed 1,3-alternate conformation and the complexes 1·K⁺(a), 1·K⁺(b), 1·K⁺(c) and 1·K⁺K⁺ were performed using density functional theory (DFT) at the B3LYP/6-31G* level. The fully optimized geometric structures of the free macroligand and its 1:1 and 1:2 complexes, as obtained from DFT calculations, were used to perform natural bond orbital (NBO) analysis. The two main types of driving force metal–ligand and cation–π interactions were investigated. NBO analysis indicated that the stabilization interaction energies (E ₂) for O…K⁺ and N…K⁺ are larger than the other intermolecular interactions in each complex. The significant increase in electron density in the RY* or LP* orbitals of K⁺ results in strong host–guest interactions. In addition, the intermolecular interaction thermal energies (ΔE, ΔH, ΔG) were calculated by frequency analysis at the B3LYP/6-31G* level. For all structures, the most pronounced changes in the geometric parameters upon interaction are observed in the calix[4]arene molecule. The results indicate that both the intermolecular electrostatic interactions and the cation–π interactions between the metal ion and π orbitals of the two pairs that face the inverted benzene rings play a significant role.     

The possible formation of nano-clusters of succinic acid and maleic acid in tetrahydrofuran due to incomplete solvation

In continuation of our work on the conformational analysis of succinic acid (SA) and maleic acid (MA) in different solvents, we present here the experimental dielectric and IR and also the ab initio Hartree–Fock calculations of the two dicarboxylic acids in tetrahydrofuran (THF). The dielectric measurements are carried out at microwave X-band frequency of 9.7 GHz and the calculations are performed at STO-3G and 6-31G(d) basis sets. The dielectric data and the dipole moment determined experimentally are compared with the dipole moment determined from the conformal analysis. It is seen that the dielectric properties of SA/MA in THF are much different from that of SA/MA in 1-4, dioxane (1-4D) that we had reported previously. The IR spectra of SA–THF system is also reported here. The present study indicates the possible formation of nano-clusters of SA/MA in THF due to incomplete solvation by THF.     

Electronic and molecular structure of M-DNA fragments

To study M-DNA molecular structure (such DNA with transition metal ions placed between the nucleic bases is able to conduct the electric current) and its conductivity mechanisms, we carried out ab initio quantum-mechanical calculations of electronic and spatial structures, thermodynamic characteristics of adenine-thymine (АТ) and guanine-cytosine (GC) base pair complexes with Zn²⁺ and Ni²⁺. To take into account the influence of the alkaline environment, calculations for these complexes were also carried out with hydroxyl and two water molecules. Computations were performed at MP2 level of theory using 6–31+G* basis set. Analogous calculations were carried out for (AC)(TG) stacking dimer of nucleic acid base pairs with two Zn²⁺. The calculation of the interaction energy in complexes has shown the preference of locating the metal ion (instead of the imino proton) between bases in M-DNA. The electronic transition energy calculation has revealed the reduction of the first singlet transition energy in АТ and GC complexes with Ni²⁺ from 4.5 eV to 0.4 - 0.6 eV. Ni²⁺ orbitals take part in the formation of HOMO and LUMO on the complexes investigated. It was shown that charges of metal ions incorporated into complexes with nucleic bases and in dimer decrease significantly.     

DFT studies on the intrinsic conformational properties of non-ionic pyrrolysine in gas phase

B3LYP/6-31G(d,p) level of theory is used to carry out a detailed gas phase conformational analysis of non-ionized (neutral) pyrrolysine molecule about its nine internal back-bone torsional angles. A total of 13 minima are detected from potential energy surface exploration corresponding to the nine internal back-bone torsional angles. These minima are then subjected to full geometry optimization and vibrational frequency calculations at B3LYP/6-31++G(d,p) level. Characteristic intramolecular hydrogen bonds present in each conformer, their relative energies, theoretically predicted vibrational spectra, rotational constants and dipole moments are systematically reported. Single point calculations are carried out at B3LYP/6-311++G(d,p) and MP2/6-31++G(d,p) levels. Six types of intramolecular H-bonds, viz. O…H–O, N…H-O, O…H–N, N…H–N, O…H–C and N…H–C, are found to exist in the pyrrolysine conformers; all of which contribute to the stability of the conformers. The vibrational frequencies are found to shift invariably toward the lower side of frequency scale corresponding to the presence of intramolecular H-bond interactions in the conformers.     

Theoretical Studies Using an Ab Initio and Molecular Modelling Combination Method on the Binding of Sequence Recognition Altered Bis-Benzimidazoles to the Minor Groove of DNA

Abstract

Ab initio calculations (Hartree-Fock) using the 3–21G and the STO-3G Gaussian basis sets were performed on synthetic analogues of the minor groove binding bis-benzimidazole Hoechst 33258 designed to exhibit altered sequence recognition. Geometry optimized conformations, energies and distribution of electrostatic charges within the molecule were derived. The binding of the optimized conformations of the drug to both alternating and non-alternating (AT)_n and (GC)_n sequences were studied.     

Spectroscopic investigation on glutamic acid by Coulomb-attenuating and double hybrid density functional theory methods

This study deals with the identification of glutamic acid by means of quantum chemical approach. FT-IR, FT-Raman and UV–vis spectra were recorded in the region 4000–400, 4000–50 cm^? 1 and 200–600 nm, respectively. CAM-B3LYP/6-31G(d,p) and B2PLYP/6-31G(d,p) calculations were performed to obtain the optimised molecular structures, vibrational frequencies and corresponding vibrational assignment, thermodynamic properties and natural bonding orbital (NBO) analysis. The results show that the obtained optimised geometric parameters (bond lengths, bond angles and bond dihedrals) and vibrational frequencies were found to be in good agreement with the experimental results. The calculations of the electronic spectra were compared with the experimental ones. Furthermore, highest occupied molecular orbital and lowest unoccupied molecular orbital analyses and UV–vis spectral analysis were also performed to determine the energy band gaps and transition states. NBO analysis, calculated using density functional theory methods (CAM-B3LYP/6-31G(d,p) and B2PLYP/6-31G(d,p)), was induced to find inter-molecular atoms. ¹³C and ¹H NMR isotropic chemical shifts were calculated and the assignments made were compared with the ChemDraw Ultra values.     

Proton and phosphorus NMR studies of d-CpG(pCpG)n duplexes in solution. Helix-coil transition and complex formation with actinomycin-D.

The Watson–Crick imino and amino exchangeable protons, the nonexchangeable base and sugar protons, and the backbone phosphates for d-CpG(pCpG)_n, n = 1 and 2, have been monitored by high-resolution nmr spectroscopy in aqueous solution over the temperature range 0°–90°C. The temperature dependence of the chemical shifts of the tetramer and hexamer resonances is consistent with the formation of stable duplexes at low temperature in solution. Comparison of the spectral characteristics of the tetranucleotide with those of the hexanucleotide with temperature permits the differentiation and assignment of the cytosine proton resonances on base pairs located at the end of the helix from those in an interior position. There is fraying at the terminal base pairs in the tetranucleotide and hexanucleotide duplexes. The Watson–Crick ring imino protons exchange at a faster rate than the Watson–Crick side-chain amino protons, with exchange occurring by transient opening of the double helix. The structure of the d-CpG(pCpG)_n double helices has been probed by proton relaxation time measurements, sugar proton coupling constants, and the proton chemical shift changes associated with the helix–coil transition. The experimental data support a structural model in solution, which incorporates an anti conformation about the glycosyl bonds, C(3) exo sugar ring pucker, and base overlap geometries similar to the B-DNA helix. Rotational correlation times of 1.7 and 0.9 × 10^?9 sec have been computed for the hexanucleotide and tetranucleotide duplexes in 0.1 M salt, D₂O, pH 6.25 at 27°C. The well-resolved ³¹P resonances for the internucleotide phosphates of the tetramer and hexamer sequences at superconducting fields shift upfield by 0.2–0.5 ppm on helix formation. These shifts reflect a conformational change about the ω,ω′ phosphodiester bonds from gauche-gauche in the duplex structure to a distribution of gauche-trans states in the coil structure. Significant differences are observed in the transition width and midpoint of the chemical shift versus temperature profiles plotted in differentiated form for the various base and sugar proton and internucleotide phosphorous resonances monitoring the d-CpG(pCpG)_n helix–coil transition. The twofold symmetry of the d-CpGpCpG duplex is removed on complex formation with the antibiotic actinomycin-D. Two phosphorous resonances are shifted downfield by ～2.6 ppm and ～1.6 ppm on formation of the 1:2 Act-D:d-CpGpCpG complex in solution. Model studies on binding of the antibiotic to dinucleotides of varying sequence indicate that intercalation of the actinomycin-D occurs at the GpC site in the d-CpGpCpG duplex and that the magnitude of the downfield shifts reflects strain at the O-P-O backbone angles and hydrogen bonding between the phenoxazone and the phosphate oxygens. Actinomycin-D is known to bind to nucleic acids that exhibit a B-DNA conformation; this suggests that the d-CpG(pCpG)_n duplexes exhibit a B-DNA conformation in solution.     

Lanthanide complex coordination polyhedron geometry prediction accuracies of ab initio effective core potential calculations

The ability to efficiently and accurately predict solid-state geometries of lanthanide coordination compounds efficiently and accurately is central for the design of new ligands capable of forming stable and highly luminescent complexes. Accordingly, we present in this paper a report on the capability of various ab initio effective core potential calculations in reproducing the coordination polyhedron geometries of lanthanide complexes. Starting with all combinations of HF, B3LYP and MP2(Full) with STO-3G, 3-21G, 6-31G, 6-31G* and 6-31+G basis sets for [Eu(H₂O)₉]³⁺ and closing with more manageable calculations for the larger complexes, we computed the fully predicted ab initio geometries for a total of 80 calculations on 52 complexes of Sm(III), Eu(III), Gd(III), Tb(III), Dy(III), Ho(III), Er(III) and Tm(III), the largest containing 164 atoms. Our results indicate that RHF/STO-3G/ECP appears to be the most efficient model chemistry in terms of coordination polyhedron crystallographic geometry predictions from isolated lanthanide complex ion calculations. Moreover, both augmenting the basis set and/or including electron correlation generally enlarged the deviations and aggravated the quality of the predicted coordination polyhedron crystallographic geometry. Our results further indicate that Cosentino et al.’s suggestion of using RHF/3-21G/ECP geometries appears to be indeed a more robust, but not necessarily, more accurate recommendation to be adopted for the general lanthanide complex case. Figure Graphical visualization of unsigned mean errors, UME(Eu-L)s, involving only the interatomic distances between the europium central ion and the oxygen atoms of the coordination polyhedron of the cation nona-aqua-europium(III) for various model chemistries, all compared to the “Cambridge Structural Database 2004” crystallographic geometry     

Alternate g−g−. conformation for dinucleoside phosphates in solution

An alternative g^?g^? conformation (conformer I') for dinucleosides in solution has been deduced, based on potential energy calculations and nuclear magnetic resonance spectroscopy. This conformation is characterized by larger glycosidic torsional angles (X=94–111°) than those of conformer 1 (X=8–35°), although the other torsional angles are similar. There are thus four stable confonners (I, I', II and III) for dinucleosides in equilibrium with the open forms. The structure of conformer I' supports that of the ‘vertical’ double helix constructed by Olson (W.K. Olson. Proc. Natl. Acad. Sci. U.S.A. 74, (1977) 1775). Our data may suggest the possibility of interconversion between the vertical double helix and the regular double helix of A-form DNA, RNA or A'-form RNA.     

Conformational studies on polynucleotide chains. V. A comparison between the quantum-mechanical and the consistent force field approach

Consistent force field (CFF) calculations were performed for the sugar–phosphate–sugar fragment, taken as a model of the polynucleotide backbone. The potential-energy-function is the sum of four contributions, accounting for bond and angle deformation, torsional motions, and nonbonded interactions. Both deoxyribose and ribose systems, with either C(2′)-endo or C(3′)-endo puckering in the starting geometry of ribose rings, were considered. A fair number of minima of the conformational-energy hypersurface were found. Although the numerical method employed in the CFF context cannot solve the problem of finding the global minimum in a definite way, one of the final conformations has a total energy much more attractive than the others, and may be regarded as the most stable conformation attainable with our potential-energy function. The energy-minimization affects the puckering of the first ribose ring differently from that of the second: in general, for the C(2′)-endo system the second ring retains its starting conformation (Ψ′ = 152°), while in the first the Ψ′ is modified by up to 70°; the opposite occurs for the C(3′)-endo system. This is explained by the different positions of the two rings relative to the phosphate group.