DFT study of alpha- and beta-D-galactopyranose at the B3LYP/6-311++G** level of theory

Forty-one conformations of alpha- and beta-d-galactopyranose were geometry optimized using the B3LYP density functional and 6-311++G** basis set. Full geometry optimization was performed on different ring geometries and different hydroxymethyl rotamers (gg/gt/tg). Analytically derived Hessians were used to calculate zero point energy, enthalpy, and entropy. The lowest energy and free-energy conformation found is the alpha-gg-(4)C(1)-c chair conformation, which is of lower electronic and free energy than the lowest energy alpha-d-glucopyranose conformer because of favorable hydrogen-bonding interactions. The in vacuo calculations showed considerable ( approximately 2.2kcal/mol) energetic preference for the alpha over the beta anomer for galactopyranose in both the (4)C(1) and (1)C(4) chair conformations. Results are compared to glucopyranose and mannopyranose calculations in vacuo. Boat and skew-boat forms were found that remained stable upon gradient optimization, although many starting conformations moved to other boat forms upon optimization. As with glucopyranose and mannopyranose, the orientation and interaction of the hydroxyl groups make the most significant contributions to the conformation-energy relationship in vacuo.     

B3LYP/6-311++G** study of alpha- and beta-D-glucopyranose and 1,5-anhydro-D-glucitol: 4C1 and 1C4 chairs, (3,O)B and B(3,O) boats, and skew-boat conformations

Geometry optimization, at the B3LYP/6-311++G** level of theory, was carried out on 4C1 and 1C4 chairs, (3,O)B and B(3,O) boats, and skew-boat conformations of alpha- and beta-D-glucopyranose. Similar calculations on 1,5-anhydro-D-glucitol allowed examination of the effect of removal of the 1-hydroxy group on the energy preference of the hydroxymethyl rotamers. Stable minimum energy boat conformers of glucose were found, as were stable skew boats, all having energies ranging from approximately 4-15 kcal/mol above the global energy 4C1 chair conformation. The 1C4 chair electronic energies were approximately 5-10 kcal/mol higher than the 4C1 chair, with the 1C4 alpha-anomers being lower in energy than the beta-anomers. Zero-point energy, enthalpy, entropy, and relative Gibbs free energies are reported at the harmonic level of theory. The alpha-anomer 4C1 chair conformations were found to be approximately 1 kcal/mol lower in electronic energy than the beta-anomers. The hydroxymethyl gt conformation was of lowest electronic energy for both the alpha- and beta-anomers. The glucose alpha/beta anomer ratio calculated from the relative free energies is 63/37%. From a numerical Hessian calculation, the tg conformations were found to be approximately 0.4-0.7 kcal/mol higher in relative free energy than the gg or gt conformers. Transition-state barriers to rotation about the C-5-C-6 bond were calculated for each glucose anomer with resulting barriers to rotation of approximately 3.7-5.8 kcal/mol. No energy barrier was found for the path between the alpha-gt and alpha-gg B(3,O) boat forms and the equivalent 4C1 chair conformations. The alpha-tg conformation has an energy minimum in the 1S3 twist form. Other boat and skew-boat forms are described. The beta-anomer boats retained their starting conformations, with the exception of the beta-tg-(3,O)B boat that moved to a skew form upon optimization.     

DFT study of alpha- and beta-D-allopyranose at the B3LYP/6-311++G ** level of theory

One hundred and two conformations of alpha- and beta-D-allopyranose, the C-3 substituted epimer of glucopyranose, were geometry optimized using the density functional, B3LYP, and the basis set, 6-311++G **. Full geometry optimization was performed on different ring geometries and on the hydroxymethyl rotamers (gg/gt/tg). Analytically derived Hessians were used to calculate zero point energy, enthalpy, and entropy. The lowest energy and free energy conformation found is the alpha-tg(g-)-4C1-c conformation, which is only slightly higher in electronic (approximately 0.2 kcal/mol) and free energy than the lowest energy alpha-D-glucopyranose. The in vacuo calculations showed a small (approximately 0.3 kcal/mol) energetic preference for the alpha- over the beta-anomer for allopyranose in the 4C1 conformation, whereas in the 1C4 conformation a considerable (approximately 1.6 kcal/mol) energetic preference for the beta- over the alpha-anomer for allopyranose was encountered. The results are compared to previous aldohexose calculations in vacuo. Boat and skew forms were found that remained stable upon gradient optimization although many starting boat conformations moved to other skew forms upon optimization. As found for glucose, mannose, and galactose the orientation and interaction of the hydroxyl groups make the most significant contributions to the conformation/energy relationship in vacuo. A comparison of different basis sets and density functionals is made in the Discussion section, confirming the appropriateness of the level of theory used here.     

Noncontact dipole effects on channel permeation. II. Trp conformations and dipole potentials in gramicidin A

The four Trp dipoles in the gramicidin A (gA) channel modulate channel conductance, and their side chain conformations should therefore be important, but the energies of different conformations are unknown. A conformational search for the right-handed helix based on molecular mechanics in vacuo yielded 46 conformations within 20 kcal/mol of the lowest energy conformation. The two lowest energy conformations correspond to the solid-state and solution-state NMR conformations, suggesting that interactions within the peptide determine the conformation. For representative conformations, the electrostatic potential of the Trp side chains on the channel axis was computed. A novel application of the image-series method of. Biophys. J. 9:1160-1170) was introduced to simulate the polarization of bulk water by the Trp side chains. For the experimentally observed structures, the CHARm toph19 potential energy (PE) of a cation in the channel center is -1.65 kcal/mol without images. With images, the PE is -1.9 kcal/mol, demonstrating that the images further enhance the direct dipole effect. Nonstandard conformations yielded less favorable PEs by 0.4-1.1 kcal/mol.     

An investigation of the pyranose ring interconversion path of alpha-L-idose calculated using density functional theory

The interconversion pathways of the pyranose ring conformation of alpha-L-idose from a (4)C1 chair to other conformations were investigated using density functional calculations. From these calculations, four different ring interconversion paths and their transition state structures from the (4)C1 chair to other conformations, such as B(3,O), and (1)S3, were obtained. These four transition-state conformations cover four possible combinations of the network patterns of the hydroxyl group hydrogen bonds (clockwise and counterclockwise) and the conformations of the primary alcohol group (tg and gg). The optimized conformations, transition states, and their intrinsic reaction coordinates (IRC) were all calculated at the B3LYP/6-31G** level. The energy differences among the structures obtained were evaluated at the B3LYP/6-311++G** level. The optimized conformations indicate that the conformers of (4)C1, (2)S(O), and B(3,O) have similar energies, while (1)S3 has a higher energy than the others. The comparison of the four transition states and their ring interconversion paths, which were confirmed using the IRC calculation, suggests that the most plausible ring interconversion of the alpha-L-idopyranose ring occurs between (4)C1 and B(3,O) through the E3 envelope, which involves a 5.21 kcal/mol energy barrier.     

Ab initio computational study of beta-cellobiose conformers using B3LYP/6-311++G**

The molecular structure of 27 conformers of beta-cellobiose were studied in vacuo through gradient geometry optimization using B3LYP density functionals and the 6-311++G** basis set. The conformationally dependent geometry changes and energies were explored as well as the hydrogen-bonding network. The lowest electronic energy structures found were not those suggested from available crystallographic and NMR solution data, where the glycosidic dihedral angles fall in the region (phi, psi) approximately (40 degrees, -20 degrees ). Rather, 'flipped' conformations in which the dihedral angles are in the range (phi, psi) approximately (180 degrees, 0 degrees ) are energetically more stable by approximately 2.5 kcal/mol over the 'experimentally accepted' structure. Further, when the vibrational free energy, deltaG, obtained from the calculated frequencies, is compared throughout the series, structures with (phi, psi) in the experimentally observed range still have higher free energy ( approximately 2.0 kcal/mol) than 'flipped' forms. The range of bridging dihedral angles of the 'normal' conformers, resulting from the variance in the phi dihedral is larger than that found in the 'flipped' forms. Due to this large flat energy surface for the normal conformations, we surmise that the summation of populations of these conformations will favor the 'normal' conformations, although evidence suggests that polar solvent effects may play the dominant role in providing stability for the 'normal' forms. Even though some empirical studies previously found the 'flipped' conformations to be lowest in energy, these studies have been generally discredited because they were in disagreement with experimental results. Most of the DFT/ab initio conformations reported here have not been reported previously in the ab initio literature, in part because the use of less rigorous theoretical methods, i.e. smaller basis sets, have given results in general agreement with experimental data, that is, they energetically favored the 'normal' forms. These are the first DFT/ab initio calculations at this level of theory, apparently because of the length and difficulty of carrying out optimizations at these high levels.     

B3LYP/6-311++G** geometry-optimization study of pentahydrates of alpha- and beta-D-glucopyranose

Five water molecules were placed in 37 different configurations around alpha- and beta-D-glucopyranose in the gt, gg, and tg conformational states, and the glucose-water complexes were geometry optimized using density functionals at the B3LYP/6-311++G** level of theory. The five water molecules were organized in space and energy minimized using an empirical potential, AMB02C, and then further geometry optimized using DFT algorithms to minimum energy positions. Electronic energy, zero point vibrational energy, enthalpy, entropy, stress energy on glucose and the water cluster, hydrogen-bond energy, and relative free energy were obtained for each configuration using thermodynamic procedures and an analytical Hessian program. The lowest energy complex was that of a clustering of water molecules around the 1- and 6-hydroxyl positions of the beta-gt anomer. Configurations in which the water molecules created a favorable network completely around and under glucose were found to have low energy for both alpha and beta anomers. Calculation of the alpha/beta anomeric ratio using the zero point corrected energy gave, approximately 32/68%, highly favoring the beta anomer in agreement with the experimental approximately 36/64% value. This ratio is better than the approximately 50/50% ratio found in our previous monohydrate study. An approximate hydroxymethyl population was obtained by noting average relative energies among the three conformational states, gt, gg, and tg. In the beta anomer complexes the gt conformation was favored over the gg state, while in the alpha anomer complexes the gg state was favored over the gt conformation, with the tg conformations all being of higher energy making little or no contribution to the rotamer population. Some geometry variances, found between glucose in vacuo and glucose after interaction with water molecules, are described and account for some observed C-5-C-6 bond length anomalies reported by us previously for the vacuum glucose structures.     

Modeling ring puckering in strained systems: application to 3,6-anhydroglycosides

Different conformations of methyl 3,6-anhydroglycosides with the beta-D-galacto, alpha-D-galacto, and beta-D-gluco configurations were studied by molecular mechanics (using the program mm3) and by quantum mechanical (QM) methods at the HF/- and B3LYP/6-31+G** levels, with and without solvent emulation. Using molecular mechanics, the energies were plotted against the phi, theta puckering coordinates of Cremer and Pople. In such strained systems, only two extreme conformations of the six-membered ring are likely: (1)C(4) and B(1,4), or any one close to either of them. Results show the preponderance of a distorted chair conformation over that of the distorted boat, though the energy difference is lower and the distortions are larger for the compound with the beta-D-galacto configuration. For derivatives of this compound, experimental data in solution indicate both chair and boat forms, depending on the compound and the solvent, whereas for the remaining compounds, experimental data always show the preponderance of the chair conformation. The more accurate DFT calculations lead to the lower energy differences, suggesting that HF and MM3 underestimate the stability of the boat-like conformations. Similar studies on model compounds depict the importance of the anomeric effect in the conformational preferences.     

Nickel(II) complexes with Schiff-base ligands derived from epimeric pyranose backbones as 2,3-chelators: modeling the coordination chemistry of chitosan

Two mononuclear nickel(II) complexes with Schiff-base ligands derived from the epimeric sugars glucosamine and mannosamine have been synthesized. The X-ray crystal structure reveals a distorted octahedral geometry at the nickel(II) ions with an N4O2 donor set and the rare 2,3-chelation of the donor atoms of the carbohydrate backbone. Upon complexation only the glucopyranose ring maintains the 4C1 chair conformation, whereas the mannopyranose ring adopts the OS5 screw-boat conformation. Dimeric units of complex cations are formed by intermolecular hydrogen bonding which are further assembled by pi-stacking affording one-dimensional chains with a twofold screw symmetry.     

Effects of electric field on alamethicin bound at the lipid-water interface: a molecular mechanics study.

A systematic molecular mechanics study of the alamethicin molecule was made to determine a set of low-energy conformers in vacuo and in aqueous environment. The behavior of these conformers was investigated at the phase boundary which was modeled as a plane dividing two compartments with solvation properties of water and octanol with a constant electric field applied normal to the boundary. The calculations were performed with a molecular mechanics program for calculation of stable conformations at the phase boundary utilizing the Empiric Conformational Energy Program for Peptides force field and the Hopfinger-Scheraga solvation model. 371 minimum energy conformers of alamethicin, determined in vacuo with the build-up procedure, were used as starting conformations for energy minimization in aqueous environment and at the phase boundary. Only 49 interphase-bound structures were within 12 kcal/mol of the minima which was found. No helical structures having values close to the canonical parameters for an alpha- or 3(10)-helix were found despite the presence of eight alpha-methylalanine residues which favor the formation of these helices; four helix-like structures were found, having all negative phi, psi values. All the helical conformers have very high energies in water (approximately 14 kcal/mol), but are quite stable at the phase boundary (3.7-6.8 kcal/mol above the lowest minima found). The implications of these results for proposed mechanisms for membrane-binding and voltage-dependent gating are considered.     

Flexibility of the pyranose ring in α- and β-D-glucoses

Conformational energies of α- and β-D -glucopyranoses were computed by varying all the ring bond angles and torsional angles using semiempirical potential functions. Solvent accessibility calculations were also performed to obtain a measure of solvent interaction. The results indicate that the ⁴C₁ (D ) chair is the most favored conformation, both by potential energy and solvent accessibility criteria. The ⁴C₁ (D ) chair conformation is also found to be somewhat flexible, being able to accommodate variations up to 10° in the ring torsional angles without appreciable change in energy. Observed solid-state conformations of these sugars and their derivatives lie in the minimum-energy region, suggesting that the substituents and crystal field forces play a minor role in influencing the pyranose ring conformation. Theory also predicts the variations in the ring torsional angles, i.e., CCCC < CCCO < CCOC, in agreement with the experimental results. The boat and twist-boat conformations are found to be at least 5 kcal mol^?1 higher in energy compared to the ⁴C₁ (D ) chair, suggesting that these forms are unlikely to be present in a polysaccharide chain. The ¹C₄ (D ) chair has energy intermediate between that of the ⁴C₁ (D ) chair and that of the twist-boat conformation. The calculated energy barrier between ⁴C₁ (D ) and ¹C₄ (D ) conformations is high—about 11 kcal mol^?1.     

Conformational analysis of the 20-residue membrane-bound portion of melittin by conformational space annealing

The conformational space of the 20-residue membrane-bound portion of melittin has been investigated extensively with the conformational space annealing (CSA) method and the ECEPP/3 (Empirical Conformational Energy Program for Peptides) algorithm. Starting from random conformations, the CSA method finds that there are at least five different classes of conformations, within 4 kcal/mol, which have distinct backbone structures. We find that the lowest energy conformation of this peptide from previous investigations is not the global minimum-energy conformation (GMEC); but it belongs to the second lowest energy class of the five classes found here. In four independent runs, one conformation is found repeatedly as the lowest energy conformation of the peptide (two of the four lowest energy conformations are identical; the other two have essentially identical backbone conformations but slightly different side-chain conformations). We propose this conformation, whose energy is lower than that found previously by 1.9 kcal/mol, as the GMEC of the ECEPP/3 force field. The structure of the proposed GMEC is less helical and more compact than the previous one. It appears that the CSA method can find several classes of conformations of a 20-residue peptide starting from random conformations utilizing only its amino acid sequence information. The proposed GMEC has also been found with a modified electrostatically driven Monte Carlo method [D. R. Ripoll, A. Liwo, and H.A. Scheraga (1998) “New Developments of the Electrostatically Driven Monte Carlo Method: Test on the Membrane-Bound Portion of Melittin,” Biopolymers, Vol. 46, pp. 117–126]. © 1998 John Wiley & Sons, Inc. Biopoly 46: 103–115, 1998     

B3LYP/6-311++G** study of monohydrates of alpha- and beta-D-glucopyranose: hydrogen bonding, stress energies, and effect of hydration on internal coordinates

Twenty-six monohydrates of alpha- and beta-D-glucopyranose were studied using gradient methods at the B3LYP/6-311++G** level of theory. Geometry optimization was carried out with the water molecules at different configurations around the glucose molecule. A new nomenclature for hydrated carbohydrates was developed to describe the water configurations. Zero-point vibrational energy, enthalpy, entropy, and relative free energy were obtained using the harmonic approximation. Hydrogen-bond energies for the monohydrates range from approximately -5 to -12 kcal/mol, and the average relative free energy is approximately 5 kcal/mol. The 1-hydroxy position is the most energetically favored site for hydration, and the region between the two and three positions is the next-most favored site. A water molecule approaching alpha-D-glucose between the 1- and 2-hydroxy positions pulls the 2-hydroxyl hydrogen atom away from the 1-hydroxy oxygen atom, thus increasing the hydrogen-bond length and also increasing the alpha-D-glucose energy. The increase in energy that occurs with a similar interaction on the beta-anomer is much less effective since the hydrogen bond is much longer. Using the calculated free energies of all 26 configurations, the anomer population (alpha/beta) increases in the beta-anomer population relative to the in vacuo case by approximately 10% at the expense of the alpha-anomer, giving an (alpha/beta) ratio of approximately 50/50. This result arises from entropy contributions favoring the beta-anomer more than the alpha-anomer. From analysis of donor and acceptor hydrogen-bond lengths, excellent correlation is found between the DFT calculated distances and those taken from carbohydrate structures in the Cambridge Crystallographic Data Bank.     

DFT conformation and energies of amylose fragments at atomic resolution. Part 2: ‘band-flip’ and ‘kink’ forms of α-maltotetraose

In Part 2 of this series of DFT optimization studies of α-maltotetraose, we present results at the B3LYP/6-311++G∗∗ level of theory for conformations denoted ‘band-flips’ and ‘kinks’. Recent experimental X-ray studies have found examples of amylose fragments with conformations distorted from the usual syn forms, and it was of interest to examine these novel structural motifs by the same high-level DFT methods used in Part 1. As in Part 1, we have examined numerous hydroxymethyl rotamers (gg, gt, and tg) at different locations in the residue sequence, and include the two hydroxyl rotamers, the clockwise ‘c’ and counterclockwise ‘r’ forms. A total of fifty conformations were calculated and energy differences were found to attempt to identify those sources of electronic energy that dictate stressed amylose conformations. Most stressed conformations were found to have relative energies considerably greater (i.e., ∼4 to 12 kcal/mol) than the lowest energy syn forms. Relative energy differences between ‘c’ and ‘r’ forms are somewhat mixed with some stressed conformations being ‘c’ favored and some ‘r’ favored, with the lowest energy ‘kink’ form being an all-gg-r conformation with the ‘kink’ in the bc glycosidic dihedral angles. Comparison of our calculated structures with experimental results shows very close correspondence in dihedral angles.     

New developments of the electrostatically driven monte carlo method: Test on the membrane-bound portion of melittin

The electrostatically driven Monte Carlo (EDMC) method has been greatly improved by adding a series of new features, including a procedure for cluster analysis of the accepted conformations. This information is used to guide the search for the global energy minimum. Alternative procedures for generating perturbed conformations to sample the conformational space were also included. These procedures enhance the efficiency of the method by generating a larger number of low-energy conformations. The improved EDMC method has been used to explore the conformational space of a 20-residue polypeptide chain whose sequence corresponds to the membrane-bound portion of melittin. The ECEPP/3 (Empirical Conformational Energy Program for Peptides) algorithm was used to describe the conformational energy of the chain. After an exhaustive search involving 14 independent runs, the lowest energy conformation (LEC) (−91.0 kcal/mol) of the entire study was encountered in four of the runs, while conformations higher in energy by no more than 1.8 kcal/mol were found in the remaining runs with the exception of one of them (run 8). The LEC is identical to the conformation found recently by J. Lee, H.A. Scheraga, and S. Rackovsky [(1998) “Conformational Analysis of the 20-Residue Membrane-Bound Portion of Melittin by Conformational Space Annealing,” Biopolymers, Vol. 46, pp. 103–115] as the lowest energy conformation obtained in their study using the conformational space annealing method. These results suggest that this conformation corresponds to the global energy minimum of the ECEPP/3 potential function for this specific sequence; it also appears to be the conformation of lowest free energy. © 1998 John Wiley & Sons, Inc. Biopoly 46: 117–126, 1998     

Structural studies on the bacterial cell wall peptidoglycan pseudomurein. I. Conformational energy calculations on the glycan strands in C1 conformation and comparison with murein

Conformational energy calculations have been used to explore the conformations which may be realized for the sugar moiety of murein and pseudomurein. For the building blocks of the pseudomurein sugar strands, i.e. for the monosaccharides beta-D-N-acetylglucosamine (NAG) and alpha-L-N-acetyltalosaminuronic acid (NAT), both in C1 ring conformation, as well as for their 1,3 and 1,4 linked disaccharides, the favoured conformations were obtained. The helical parameters of sugar strands of both linkage types, which describe the regular structure of the corresponding polysaccharides, poly-(1,3-NAT-NAG) and poly-(1,4-NAT-NAG), were calculated. Both types of polysaccharides poly-(NAG-NAT) considered in this study favoured extended conformations, which in the case of 1,3 linked polymers showed less gain of length per saccharide unit compared to 1,4 linked poly-(NAG-NAT) residues. For a 1,3 linked sugar moiety of pseudomurein every pair of neighbouring peptides attached to glycan chain pointed in favoured conformations approximately to opposite sides of the strands, whereas in a 1,4 linked poly-(NAG-NAT) the peptides protruded approximately to the same side of the glycan moiety. A comparison between pseudomurein and murein revealed that the sugar moieties of both peptidoglycans have similar features in respect to their overall structure, i.e. both favoured more or less extended structures. In contrast to these data the shapes of the resulting polysaccharide moieties were remarkably different. In poly-(1,3-NAG-NAT) the glycan chains possessed a zig-zag-like arrangement, whereas for glycan chains of the murein type relatively flat structures were preferred. These remaining differences in the conformational arrangement between both peptidoglycans depend strongly on the C1 chair conformation of NAT. It is, therefore, attractive to speculate about an hypothetical pseudomurein sugar chain configuration comprising beta-L-N-acetyltalosaminuronic acid in its 1C conformation.     

Several transition states from C1 to skew conformations of β-d-glucopyranose

The interconversion pathways in the ring distortion of β-d-glucopyranose were investigated using density functional calculations. We examined the energies of several conformers of β-d-glucopyranose and tried to obtain the transition-state conformation and determine the pathway between a ⁴C₁ chair and some distorted ring conformers. The results showed that two E₃/²H₃ conformations and one E₃/⁴H₃ conformation were transition states in such ring puckering. The transition state with the lowest energy conformation is the E₃/²H₃ ring conformation with the side-chain conformation of r-ggG⁺. Intrinsic reaction coordinate calculations indicated that the E₃/²H₃ conformation with the lowest conformational energy is a transition state of the ring interconversion path between the conformations of ⁴C₁ and ²S_O/B_3,O. The energy barrier of this interconversion was 6.13 kcal/mol. As far as we know, this is the first example of finding pathways for an interconversion of glucopyranose ring puckering at the level of a quantum chemical calculation.     

Ab initio conformational analysis of nucleic acid components: intrinsic energetic contributions to nucleic acid structure and dynamics

In recent years, the use of high-level ab initio calculations has allowed for the intrinsic conformational properties of nucleic acid building blocks to be revisited. This has provided new insights into the intrinsic conformational energetics of these compounds and its relationship to nucleic acids structure and dynamics. In this article we review recent developments and present new results. New data include comparison of various levels of theory on conformational properties of nucleic acid building blocks, calculations on the abasic sugar, known to occur in vivo in DNA, on the TA conformation of DNA observed in the complex with the TATA box binding protein, and on inosine. Tests of the Hartree-Fock (HF), second-order M?ller-Plesset (MP2), and Density Functional Theory/Becke3, Lee, Yang and Par (DFT/B3LYP) levels of theory show the overall shape of backbone torsional energy profiles (for gamma, epsilon, and chi) to be similar for the different levels, though some systematic differences are identified between the MP2 and DFT/B3LYP profiles. The east pseudorotation energy barrier in deoxyribonucleosides is also sensitive to the level of theory, with the HF and DFT/B3LYP east barriers being significantly lower (approximately 2.5 kcal/mol) than the MP2 counterpart (approximately 4.0 kcal/mol). Additional calculations at various levels of theory suggest that the east barrier in deoxyribonucleosides is between 3.0 and 4.0 kcal/mol. In the abasic sugar, the west pseudorotation energy barrier is found to be slightly lower than the east barrier and the south pucker is favored more than in standard nucleosides. Results on the TA conformation suggest that, at the nucleoside level, this conformation is significantly destabilized relative to the global energy minimum, or relative to the A- and B-DNA conformations. Deoxyribocytosine would destabilize the TA conformation more than other bases relative to the A-DNA conformation, but not relative to the B-DNA conformation.     

Conformational analysis of d-tubocurarine: Implications for minimal dimensions of its binding site within ion channels

All the minimum-energy conformations of d-tubocurarine were calculated by the method of molecular mechanics. The energy was minimized from 413 closed forms of the 18-member ring. The set of minimum-energy conformations includes 10 forms with energies less than 6 kcal/mol from the most stable one. Among the four lowest minimum-energy conformations, two forms correspond to those known from X-ray studies, whereas two conformations were not detected experimentally earlier. The flexibility of d-tubocurarine was estimated by calculating six paths of interconversion between the four lowest minimum-energy conformations. Using a molecular graphics technique, it was found that the most extended minimum-energy conformation of d-tubocurarine may fit in an ion channel of a rectangular profile of 8.7 × 11.2 Å, while one tetrahydroisoquinoline head may fit a profile as small as 6.9 × 11.0 Å. A possible model of d-tubocurarine location within the ion channel of the neuronal nicotinic acetylcholine receptor is suggested.     

Conformational preferences of 1,4,7-trithiacyclononane: a molecular mechanics and density functional theory study

Conformational preferences of 1,4,7-trithiacyclononane were studied using a highly efficient sampling technique based on local nonstochastic deformations and the MM2(91) force field. The results show that conformers that the molecule adopts in the crystal state were found to be low-energy conformers (LECs) within 5 kcal mol(-1) of the global minimum. A conformation with C1 symmetry was the global minimum and the C3 and C2 conformations were calculated to be 0.03 and 1.78 kcal mol(-1) higher in energy, respectively. The structures were further minimized using Density Functional Theory (DFT) calculations with two different functionals. The C2 and the C1 conformations were found to be LECs with the C3 conformation more than 4.0 kcal mol(-1) above the global minimum. The relative energies and structural ordering obtained using the BP86 functional are in agreement with the previously reported relative energies calculated using second-order Moller-Plesset (MP2) ab initio calculations. With the energy ordering being dependent on the molecular mechanics force field used, the approach of MM-->DFT (searching exhaustively the available conformational space at the MM level followed by generating the energy ordering through DFT calculations) appears to be appropriate for thiacrown ethers.