Pseudorotational Angle Research Articles

A local mode study of ring puckering effects in the infrared spectra of cyclopentane.

We report and interpret recently recorded high-resolution infrared spectra for the fundamentals of the CH2 scissors and CH stretches of gas phase cyclopentane at -26.1 and -50 °C, respectively. We extend previous theoretical studies of this molecule, which is known to undergo barrierless pseudorotation due to ring puckering, by constructing local mode Hamiltonians of the stretching and scissor vibrations for which the frequencies, couplings, and linear dipoles are calculated as functions of the pseudorotation angle using B3LYP/6-311++(d,p) and MP2/cc-pVTZ levels of theory. Symmetrization (D5h) of the vibrational basis sets leads to simple vibration/pseudorotation Hamiltonians whose solutions lead to good agreement with the experiment at medium resolution, but which miss interesting line fractionation when compared to the high-resolution spectra. In contrast to the scissor motion, pseudorotation leads to significant state mixing of the CH stretches, which themselves are Fermi coupled to the scissor overtones.

Variability of the DNA Backbone Geometry in DNA-Protein Complexes: Experimental Data Analysis.

We have analyzed and compared the available experimental data (PDB) on the backbone geometry of the DNA in solution (NMR), in crystals (X-rays), and in complexes with proteins (X-rays and cryo-electron microscopy). The deoxyribose (pseudorotational angle τ0) and ε/ζ (BI-BII transition in phosphates) flexibilities are practically the same in the four samples. The α/γ mobility is minimal in crystalline DNA: on the histograms, there is one canonical and one noncanonical t/t peak. The α/γ mobility increases in DNA solutions (three more noncanonical peaks) and is maximal in DNA-protein complexes (another additional peak). On a large amount of data, we have confirmed that the three main degrees of freedom of the sugar-phosphate backbone are "orthogonal": changes in any of the angles τ0, (ζ-ε), and (γ-α) occur, as a rule, at a constant (usually canonical) value of any other. In the DNA-protein complexes, none of the geometrical parameters commonly used to distinguish the A and B forms of DNA, except for Zp and its simpler analog Zp', show an unambiguous correlation with τ0. Proteins, binding to DNA, in 59% of cases change the local shape of the helix up to the characteristic of the A-form without switching the deoxyribose conformation from south to north. However, we have found simple local characteristics of one nucleotide that correlate with the angles τ0 and (ζ-ε). These are the angles C3'C1'N* and C4'C3'P(2), respectively. They are orthogonal in DNA-protein complexes exactly as the pair τ0 and (ζ-ε). Most characteristics of DNA in complexes with proteins are the same in X-ray and in cryo-EM data, except for the histogram for the angle τ0. We offer a possible explanation for this difference. We also discuss the artifacts on the ε/ζ histogram for DNA in solutions caused by the currently used NMR refinement protocols.

Interplay of Rotational and Pseudorotational Motions in Flexible Cyclic Molecules.

Solutions to the time-independent nuclear Schrödinger equation associated with the pseudorotational motion of three flexible cyclic molecules are presented and discussed. Structural relaxations related to the pseudorotational motion are described as functions of a pseudorotation angle ϕ which is formulated according to the definition of ring-puckering coordinates originally proposed by Cremer and Pople (J. Am. Chem. Soc.1975, 97 ( (6), ), 1354−1358). In order to take into account the interplay between pseudorotational and rotational motions, the rovibrational Hamiltonian matrices are formulated for the rotational quantum numbers J = 0 and J = 1. The rovibrational Hamiltonian matrices are constructed and diagonalized using a Python program developed by the authors. Suitable algorithms for (i) the construction of one-dimensional cuts of potential energy surfaces along the pseudorotation angle ϕ and (ii) the assignment of the vibrorotational wave functions (which are needed for the automatic calculation of rotational transition energies J = 0 → J = 1) are described and discussed.

Structural and functional computational analysis of nicotine analogs as potential neuroprotective compounds in Parkinson disease

As the mechanism of interaction between nicotinic receptors with nicotine analogs is not yet fully understood, information at molecular level obtained from computational calculations is needed. In this sense, this work is a computational study of eight nicotine analogs, all with pyrrolidine ring modifications over a nicotine-based backbone optimized with B3LYP-D3/aug-cc-pVDZ. A molecular characterization was performed focusing on geometrical parameters such as pseudo-rotation angles, atomic charges, HOMO and LUMO orbitals, reactivity indexes and intermolecular interactions. Three analogs, A2 (3-(1,3-dimethyl-4,5-dihydro-1h-pirazole-5-yl) pyridine), A3 (3-(3-methyl-4,5-dihydro-1H-pyrazol-5-yl)-pyridine) and A8 (5-methyl-3-(pyridine-3-yl)-4,5-dihydroisoxazole), were filtered suggesting putative neuroprotective activity taking into account different reactivity values, such as their lowest hardness: 2.37 eV (A8), 2.43 eV (A2) and 2.56 eV (A3), compared to the highest hardness value found: 2.71 eV for A5 (3-((2S,4R)-4-(fluoromethyl)-1-methylpyrrolidine-2-il) pyridine), similar to the value of nicotine (2.70 eV). Additionally, molecular docking of all 8 nicotine analogs with the α 7 nicotinic acetylcholine receptor (α 7 nAChR) was performed. High values of interaction between the receptor and the three nicotine analogs were obtained: A3 (-7.1 kcal/mol), A2 (-6.9 kcal/mol) and A8 (-6.8 kcal/mol); whereas the affinity energy of nicotine was -6.4 kcal/mol. Leu116 and Trp145 are key residues in the binding site of α 7 nAChR interacting with nicotine analogs. Therefore, based upon these results, possible application of these nicotine analogs as neuroprotective compounds and potential implication at the design of novel Parkinson's treatments is evidenced.

Structure of an RNA/DNA dodecamer corresponding to the HIV-1 polypurine tract at 1.6 Å resolution

The crystal structure of an RNA/DNA hybrid dodecamer, r(5'-uaaaagaaaagg):d(5'-CCTTTTCTTTTA), which contains three-quarters of the polypurine tract (PPT) sequence of theHIV RNA genome is reported. The hybrid structure was determined at 1.6 Å resolution and was found to have the A-form conformation. However, the presence of alternate conformations along the RNA template strand indicated increased flexibility of the PPT sequence. Two segments (atnucleotides 1-2 and 6-8) of the RNA chain have two conformations exhibiting differences in torsion and pseudorotation angles. For conformation I((1-2), (6-8)), 25% of the RNA sugars have the C2'-exo pucker and the rest have the expected C3'-endo pucker. The II(1-2) and II(6-8) conformations of the RNA strand have one sugar with the C2'-exo pucker. None of theribose rings exist in the C2'-endo form, in contrast to a previous report which postulated a C2'-endo ribose as a key structural element of the PPT. The widths of the minor groove for conformations I((1-2), (6-8)) and II((1-2), (6-8)) of the RNA strand are 9.2-10.5 and 9.4-10.7 Å, respectively. Both ranges are very close to the intervals accepted for A-form RNA duplexes. On the opposing DNA primer strand most of the sugars are C3'-endo, except for the 3'-terminal sugars, which are C2'-endo (T22) or O4'-endo (T23 and A24). The duplex includes a noncanonical u1(anti)·A24(syn) base interaction with only one hydrogen bond between the bases. This noncanonical base interaction at the 5'-end of the template distorts the values of the helical parameters of the adjacent base pair.

Geometry of the C-/N-Terminal Connection of Tropomyosin on F-Actin Based on Molecular Dynamics Simulations of the Thin Filament

The binding of tropomyosin molecules to actin filaments is extremely weak and only is effective because of cooperative interactions resulting from the head-to-tail polymerization of tropomyosin on thin filaments. Experiments show that interfering with the head/tail interactions prevents actin-binding (Johnson and Smillie, 1977; Heald and Hitchcock-DeGregori, 1988). Crystal- and NMR-structures of mixtures of head/tail-fragments stabilized with various foreign adducts, yet unconstrained by interactions with F-actin, have suggested that head-to-tail linkage results from a splayed C-terminal coiled-coil region enveloping a more compact N-terminus to form a four-helix nexus with varying degree of overlap. However, it is possible that structures observed previously may have been influenced by the adducts and the absence of F-actin. In fact, N-/C-terminal sequences vary greatly between tropomyosin isoforms and are not strictly tetra-coiled-coiled sequences, complicating interpretation. Here, based on static and dynamic models of αα-striated muscle tropomyosin on F-actin, we derive the overlap and the relative pseudo-rotation angle (which determines the face of tropomyosin which is seen by F-actin) between the C- and N-terminal ends. Our recently described atomic model of tropomyosin on F-actin (Li et al., 2011) suggests that, when constrained by interaction with the F-actin surface, the tropomyosin ends overlap by five amino acids, slightly shorter than in the adduct-structures. MD simulations provide a measure of the variability in end-to-end distance of tropomyosin on F-actin (and therefore of the possible overlap length) and in the degree of pseudo-rotation between N- and C-tropomyosin termini (and therefore C-/N-nexus geometry). These simulations demonstrate local C-terminal coiled-coil splaying, relatively short overlap (7.9 ± 2.2A, i.e. ∼four to seven residues) and predict a rotation of 17.5 ± 13° at the end of one tropomyosin relative to the beginning of the next.

How Flexible are DNA Constituents? The Quantum-Mechanical Study

Relaxed force constants (RFCs) and vibrational root-mean-square deviations have been evaluated by the original calculation method for conformational parameters of the DNA structural units and their constituents: nucleic acid bases (uracile, thymine, cytosine, adenine and guanine) and their ‘building blocks’ (benzene, pyrimidine, imidazole and purine molecules), as well as the DNA backbone structural units: tetrahydrofuran, 1,2-dideoxyribose, methanol and orthophosphoric acid. It has been found that the RFCs for nomenclature torsions β, γ, ε and sugar pseudorotation angle P in 1,2-dideoxyribose are sensible to the molecule conformation and their values are in the range of 1–25 kcal/(mole-rad2) obeying the inequality K γ > K ε > KP > K β. The RFCs values for endocyclic torsions of nucleic acid bases six-member rings lie within 15–45 kcal/(mole·rad2) in pyrimidines and within 20–60 kcal/(mole·rad2) in purines. It is shown that the quantum zero-point motion effectively neglects the amino group non-planarity in cytosine, adenine and partially in guanine.

Molecular structure and character of bonding of mono and divalent metal cations (Li+, Na+, K+, Mg2+, Ca2+, Zn2+, and Cu+) with guanosine: AIM and NBO analysis

The B3LYP/6-311++G (d,p) density functional approach was used to study the gas-phase metal affinities of Guanosine (ribonucleoside) for the Li+, Na+, K+, Mg2+, Ca2+, Zn2+, and Cu+ cations. In this study we determine coordination geometries, binding strength, absolute metal ion affinities, and free energies for the most stable products. We have also compared the results for Guanosine, with our previously reported results for 2′-Deoxyguanosine. Based on the results, it is obvious that MIA is strongly dependent on the charge-to-size ratio of the cation. Guanosine interacts more strongly with Zn2+ than do with Mg2+, Ca2+, and Cu+ and therefore stronger interactions lead to higher MIA. In both free molecules and their complexes, the Syn orientation of the base is stabilized by an intramolecular O5′–H···N3 hydrogen bond and the anti orientation of the base is stabilized by an intramolecular C–H···O hydrogen bond formed between the (C8-H8) and the O5′ atom of the sugar moiety. It is also interesting to mention that linear correlation between calculated MIA values and the atomic numbers (Z) of the metal ions of Li+, Na+, and K+ were found. Furthermore, the influences of metal cationization on the strength of the N-glycosidic bond, torsion angles, angle of pseudorotation (P), and intramolecular C–H···O and O–H···O hydrogen bonds have been studied. Natural bond orbital (NBO) analysis was performed to calculate the charge transfer and natural population analysis of the complexes. Quantum theory of atoms in molecules (QTAIM) was also applied to determine the nature of interactions.

On the binding of Mg2+, Ca2+, Zn2+ and Cu+ metal cations to 2 ′ -deoxyguanosine: Changes on sugar puckering and strength of the [formula omitted]-glycosidic bond

The binding of Mg2+, Ca2+, Zn2+ and Cu+ metal cations to 2′-deoxyguanosine has been analyzed, using the hybrid B3LYP, Density Functional Theory (DFT) method and 6–311++G(d,p) orbital basis sets. Coordination geometries, absolute metal ion affinities, and free energies for the most stable complexes formed by Mg2+, Ca2+, Zn2+ and Cu+ with the nucleoside, 2′-deoxyguanosine, have been determined. Furthermore, the influences of metal cationization on the strength of the N-glycosidic bond, torsion angles and angle of Pseudorotation (P) have been studied. With respect to the results, it has been found that metal binding significantly changes the values of the phase angle of Pseudorotation (P) in the sugar unit of these nucleosides. In all modified forms, the length of the C1′–N9 bond increases. The Mulliken population analysis on atomic charges has been carried out on the optimized geometries. Natural Bond Orbital (NBO) analysis was also performed to calculate the charge transfer and natural population analysis of the complexes. Quantum Theory of Atoms In Molecules (QTAIM) was also applied to determine the nature of interactions.

DFT STUDY ON METAL CATIONIZATION AND O6-PROTONATION ON 2′-DEOXYGUANOSINE CONFIGURATION: CHANGES ON SUGAR PUCKERING AND STRENGTH OF THE N-GLYCOSIDIC BOND

Density functional theory (DFT) calculations were performed at the B3LYP/6-311++G (d, p) level to determine coordination geometries, absolute metal ion affinities, and free energies for all possible complexation stable products formed by monovalent metal cations including Li+ , Na+ , K+ with the nucleoside 2′-Deoxyguanosine. All computations indicate that the metal ion affinity (MIA) decreases on going from Li+ to Na+ and K+ for 2′-Deoxyguanosine. For example, the affinities for the metal ions described above are 75.2, 57.3, 43.4 kcal/mol, respectively. Furthermore, the influence of metal cationization and O6 protonation on the strength of the N-glycosidic bond, torsion angles and angle of pseudorotation (P) have been studied. With respect to the results, it has been found that metal binding significantly change the values of the phase angle of pseudorotation (P) in the sugar unit of these nucleosides but, O6 -protonation do not significantly change the values of the torsion angles and angle of pseudorotation (P). In all modified forms, the length of the C1'–N9 bond increases. The Mulliken population analysis and natural bond orbital (NBO) analysis on atomic charges have been carried out on the optimized geometries.

Free-Energy Landscape of RNA Hairpins Constructed via Dihedral Angle Principal Component Analysis

To systematically construct a low-dimensional free-energy landscape of RNA systems from a classical molecular dynamics simulation, various versions of the principal component analysis (PCA) are compared: the cPCA using the Cartesian coordinates of all atoms, the dPCA using the sine/cosine-transformed six backbone dihedral angles as well as the glycosidic torsional angle chi and the pseudorotational angle P, the aPCA which ignores the circularity of the 6 + 2 dihedral angles of the RNA, and the dPCA(etatheta), which approximates the 6 backbone dihedral angles by 2 pseudotorsional angles eta and theta. As representative examples, a 10-nucleotide UUCG hairpin and the 36-nucleotide segment SL1 of the Psi site of HIV-1 are studied by classical molecular dynamics simulation, using the Amber all-atom force field and explicit solvent. It is shown that the conformational heterogeneity of the RNA hairpins can only be resolved by an angular PCA such as the dPCA but not by the cPCA using Cartesian coordinates. Apart from possible artifacts due to the coupling of overall and internal motion, this is because the details of hydrogen bonding and stacking interactions but also of global structural rearrangements of the RNA are better discriminated by dihedral angles. In line with recent experiments, it is found that the free energy landscape of RNA hairpins is quite rugged and contains various metastable conformational states which may serve as an intermediate for unfolding.

Interaction of Mg 2+, Ca 2+, Zn 2+ and Cu + with cytosine nucleosides: Influence of metal on sugar puckering and stability of N-Glycosidic bond, a DFT study

The coordination geometries, electronic features, absolute metal ion affinities, geometrical multiplicity and binding strength for the complexes formed by Mg 2+, Ca 2+, Zn 2+ and Cu + with cytidine and deoxycytidine have been determined theoretically employing the hybrid B3LYP exchange-correlation functional and using 6-311++G(d,p) orbital basis sets. Our study shows that the most stable coordination modes of cytosine nucleosides with these cations are rather similar. In the isolated state, the most stable form corresponds to a tridentate complex in which the cation interacts with three oxygen atoms: one from the carbonyl oxygen of cytosine moiety and two from the sugar unit (i.e., O5′ and O4′). The influences of metal cationization on the strength of the N-glycosidic bond, the values of the phase angle of pseudorotation P in the sugar unit and sugar ring conformation have been studied. In all cases, the N1 C1′ bond distance changes upon introducing positive charge in the nucleosides.

Studies on the torsions of nucleic acids using ABEEMσπ/MM method

The biological functions of DNA are related to its conformation. In this work, potential energy profiles of nucleic acid torsions associated with the sugar-phosphate backbone and the glycosyl linkage were studied using the ab initio and ABEEMσπ/MM methods. The ABEEMσπ/MM well reproduced the entire potential energy profiles (not just minima) of the torsion of model compounds. The ABEEMσπ/MM method also gives reasonable results about the structures and dipole moments of DNA bases and interactions between bases and one water molecule. The torsions of four deoxyribonucleosides were investigated in aqueous solution through molecular dynamics simulations based on the ABEEMσπ/MM. Reasonable distributions of the torsional angles were observed with the pseudorotation angles in south, glycosyl torsion χ in anti conformation, and g + conformations for torsion γ. The reasonable results on simulating DNA fragments in gas and aqueous solution suggest that the newly constructed ABEEMσπ/MM method would be robust for simulating nucleic acids.

DFT study of the interaction of cytidine and 2′-deoxycytidine with Li +, Na +, and K +: effects of metal cationization on sugar puckering and stability of the N-glycosidic bond

Density functional theory (DFT) calculations were performed at the B3LYP level with a 6-311++G(d,p) basis set to systematically explore the geometrical multiplicity and binding strength for complexes formed by Li +, Na +, and K + with cytidine and 2′-deoxycytidine. All computational studies indicate that the metal ion affinity (MIA) decreases from Li + to Na + and K + for cytosine nucleosides. For example, for cytidine the affinity for the above metal ions are 79.5, 55.2, and 41.8 and for 2′-deoxycytidine, 82.8, 57.4, and 42.2 kcal/mol, respectively. It is also interesting to mention that linear correlations between calculated MIA values and the atomic numbers ( Z) of the above metal ions were found. The influence of metal cationization on the coordination modes and the strength of the N-glycosidic bond in cytosine nucleosides have been studied. In all cases, the N1–C1′ bond distance changes upon introducing a positive charge in the nucleosides. It has been found that metal binding significantly changes the values of the phase angle of pseudorotation P in the sugar unit of these nucleosides. With respect to the sugar ring, metal binding changes the values of the glycosyl torsion angle and sugar ring conformation. The present calculations in the gas phase provide the first clues on the intrinsic chemistry of these systems and may be of value for studies of the influence of metal cations on the conformational behavior and function of nucleic acids.

Experimental and Theoretical Electron Momentum Spectroscopic Study of the Valence Electronic Structure of Tetrahydrofuran under Pseudorotation

The most populated structure of tetrahydrofuran (THF) has been investigated in our previous study using electron momentum spectroscopy (EMS). Because of the relatively low impact energy (600 eV) and low energy resolution (DeltaE = 1.20 eV) in the previous experiment, only the highest occupied molecular orbital (HOMO) of THF was investigated. The present study reports the most recent high-resolution EMS of THF in the valence space for the first time. The binding energy spectra of THF are measured at 1200 and 2400 eV plus the binding energies, respectively, for a series of azimuthal angles. The experimentally obtained binding energy spectra and orbital momentum distributions (MDs) are employed to study the orbital responses of the pseudorotation motion of THF. The outer valence Greens function (OVGF), the OVGF/6-311++G** model, and density function theory (DFT)-based SAOP/et-pVQZ model are employed to simulate the binding energy spectra. The orbital momentum distributions (MDs) are produced using the DFT-based B3LYP/aug-cc-pVTZ model, incorporating thermodynamic population analysis. Good agreement between theory and experiment is achieved. Orbital MDs of valence orbitals exhibit only slight differences with respect to the impact energies at 1200 and 2400 eV, indicating validation of the plane wave impulse approximation (PWIA). The present study has further discovered that the orbital MDs of the HOMO in the low-momentum region (p < 0.70 a.u) change significantly with the pseudorotation angle, phi, giving a v-shaped cross section, whereas the innermost valence orbital of THF does not vary with pseudorotation, revealing a very different bonding mechanism from the HOMO. The present study explores an innovative approach to study pseudorotation of sugar puckering, which sheds a light to study other biological systems with low energy barriers among ring-puckering conformations.

Comprehensive Conformational Analysis of the Nucleoside Analogue 2‘-β-Deoxy-6-azacytidine by DFT and MP2 Calculations

A comprehensive conformational analysis of isolated 2'-beta-deoxy-6-azacytidine (d6AC), an analogue of therapeutically active 6-azacytidine (6AC), has been performed by means of ab initio calculations at the MP2/6-311++G(2df,pd)//DFT B3LYP/6-31G(d,p) level of theory. Among the 81 conformers located within a 7.83 kcal/mol Gibbs energy range at T = 298.15 K, 38 contain syn-oriented bases with respect to 2'-deoxyribose; the other conformers include anti-oriented bases. Energetic analysis of these conformers shows that conformational equilibrium of isolated d6AC at T = 298.15 K is shifted to syn conformation with a syn/anti ratio estimated as 61.4%:38.6%. As far as the sugar conformation is concerned, 40 conformers contain north (N) (with 0.3 degrees < or = P < or = 40.1 degrees), and the rest possess south (S) (with 157.1 degrees < or = P < or = 207.0 degrees) puckers, where P is the pseudorotational angle of the furanose ring. The S/N occupancy ratio is estimated as 80.2%:19.8% (T = 298.15 K). The two most stable conformers are energetically quasidegenerate and correspond to both C2'-endo/syn conformers differing only by orientation of the O3'H hydroxyl group. They are both stabilized by means of similar intramolecular H-bonds, i.e., O5'H...O2, C2'H2...O2, and C2'H2...O5'. As examined by AIM criteria, from 1 to 3 H-bonds per conformer were identified among 13 possible interactions: O5'H...O2, O5'H...N6, O3'H...O5', O5'H...O3', C1'H...O2, C2'H2...O2, C2'H2...O5', C3'H...O2, C3'H...N6, C5'H1...O2, C5'H2...O2, C5'H1...N6, and C5'H2...N6. The biological effect of d6AC is conceived as an inhibition of replicative DNA polymerase caused by an unusual orientation of the sugar residue against the base in the only A form DNA-like conformer.

Effects of Ionization, Metal Cationization and Protonation on 2‘-Deoxyguanosine: Changes on Sugar Puckering and Stability of the N-Glycosidic Bond

The influence of oxidation, protonation, and metal cationization with Cu(+) and Cu(2+) on the strength of the N-glycosidic bond in 2'-deoxyguanosine has been studied by means of quantum chemical calculations. In all cases, the N9-C1' bond distance increases (0.03-0.06 A) upon introducing positive charge in the guanine moiety, the observed variations being more important for the dicationic systems. Binding energies show that the effect of X(n)(+) in guanine hinders the homolytic dissociation, whereas it largely favors the heterolytic process. With respect to the deoxyribose ring, it has been found that metal binding, oxidation, and protonation do not significantly change the values of the phase angle of pseudorotation P. However, the glycosyl torsion angle chi varies considerably (from 242.0 degrees to 189.8 degrees) as a consequence of a stabilizing guanine-sugar (H8-O4') interaction due to the increase of acidity of guanine C8-H8 upon cationization.

Free Energy Landscape of A-DNA to B-DNA Conversion in Aqueous Solution

The interconversion between the well-characterized A- and B-forms of DNA is a structural transition for which the intermediate states and the free energy difference between the two endpoints are not known precisely. In the present study, the difference between the Root Mean Square Distance (RMSD) from canonical A-form and B-form DNA is used as an order parameter to characterize this free energy difference using umbrella sampling molecular dynamics (MD) simulations with explicit solvent. The constraint imposed along this order parameter allows relatively unrestricted evolution of the intermediate structures away from both canonical A- and B-forms. The free energy difference between the A- and B-forms for the hexamer DNA sequence CTCGAG in aqueous solution is conservatively estimated to be at least 2.8 kcal/mol. A continuum of intermediate structures with no well-defined local minima links the two forms. The absence of any major barriers in the free energy surface is consistent with spontaneous conversion of the A-form DNA to B-form DNA in unconstrained simulations. The extensive sampling in the MD simulations (>0.1 mus) also allowed quantitative energetic characterization of local backbone conformational variables such as sugar pseudorotation angles and BI/BII state equilibria and their dependence on base identity. The absolute minimum in the calculated free energy profile corresponds closely to the crystal structure of the hexamer sequence, indicating that the present method has the potential to identify the most stable state for an arbitrary DNA sequence in water.

Applying the stochastic difference equation to DNA conformational transitions: A study of B–Z and B–A DNA transitions

Despite the existence of numerous models to account for the B-Z DNA transition, experimenters have not yet arrived at a conclusive answer to the structural and dynamical features of the B-Z transition. By applying the stochastic difference equation to simulate the B-Z DNA transition, we have shown that the stretched intermediate model of the B-Z transition is more probable than other B-Z transition models such as the Harvey model. This is accomplished by comparing potential energy profiles of various B-Z DNA transition models and calculating relative probabilities based on the stochastic difference equation with respect to length (SDEL) formalism. The results garnered in this article allow for new approaches in determining the structural transition of B-DNA to Z-DNA experimentally. We have also simulated the B-A DNA transition using the stochastic difference equation. Unlike the B-Z DNA transition, the mechanism for the B-A DNA transition is well established. The variation in the pseudorotation angle during the transition is in good agreement with experimental results. Qualitative features of the simulated B-A transition also agree well with experimental data. The SDEL approach is thus a suitable numerical technique to compute long-time molecular dynamics trajectory for DNA molecules.

Theoretical Study of 9-β-d-Erythrofuranosyladenine and Corresponding Carbocyclic Analogues. Evidence for a Base-Activated Conformational Lock

The conformational surfaces of three nucleoside analogues have been investigated computationally, where an adenine is attached to a diol of tetrahydrofuran, a diol of cyclopentane, and a diol of cyclopentene. In each system, the lowest-energy conformer displays a conformational lock into the south position by an internal hydrogen bond between O2‘H of the five-membered ring and the N3 nitrogen of adenine. When aqueous solvation is accounted for by the PCM method, the preference for the locked conformer is diminished. A pseudorotation angle of 9-(trans-2‘,trans-3‘dihydroxycyclopentyl)adenine has been determined to be 176.8° by fitting the measured 3JHH values using PSEUROT which is in good agreement with the calculated value of 169.3°.