Coupling of the guanosine glycosidic bond conformation and the ribonucleotide cleavage reaction: implications for barnase catalysis

To examine the possible relationship of guanine-dependent GpA conformations with ribonucleotide cleavage, two potential of mean force (PMF) calculations were performed in aqueous solution. In the first calculation, the guanosine glycosidic (Gchi) angle was used as the reaction coordinate, and computations were performed on two GpA ionic species: protonated (neutral) or deprotonated (negatively charged) guanosine ribose O2 '. Similar energetic profiles featuring two minima corresponding to the anti and syn Gchi regions were obtained for both ionic forms. For both forms the anti conformation was more stable than the syn, and barriers of approximately 4 kcal/mol were obtained for the anti --> syn transition. Structural analysis showed a remarkable sensitivity of the phosphate moiety to the conformation of the Gchi angle, suggesting a possible connection between this conformation and the mechanism of ribonucleotide cleavage. This hypothesis was confirmed by the second PMF calculations, for which the O2 '--P distance for the deprotonated GpA was used as reaction coordinate. The computations were performed from two selected starting points: the anti and syn minima determined in the first PMF study of the deprotonated guanosine ribose O2'. The simulations revealed that the O2 ' attack along the syn Gchi was more favorable than that along the anti Gchi: energetically, significantly lower barriers were obtained in the syn than in the anti conformation for the O--P bond formation; structurally, a lesser O2 '--P initial distance, and a better suited orientation for an in-line attack was observed in the syn relative to the anti conformation. These results are consistent with the catalytically competent conformation of barnase-ribonucleotide complex, which requires a guanine syn conformation of the substrate to enable abstraction of the ribose H2 ' proton by the general base Glu73, thereby suggesting a coupling between the reactive substrate conformation and enzyme structure and mechanism.     

Conformation of adenosine 3',5'-monophosphate in solution as studied by the NMR-desert method. II. Self-association and temperature-dependent glycosidic isomerization at pH 7.

A study has been made of the association and the temperature-dependent conformation of adenosine 3',5'-monophosphate (cyclic AMP) in a neutral aqueous (2H2O) solution by means of proton magnetic resonance chemical shift and relaxation. The concentration and temperature-dependent chemical shifts of H(1'), H(2), and H(8), have enabled us to estimate the self-association constant, Ka = 1.1 +/- 0.3 M-1 at 25 degrees C and thermodynamic parameters delta H = -5.8 +/- 1.5 kcal/mol and delta S (25 degrees C) = -19.0 +/- 3 cal/mol per degree. The NMR-DESERT (Deuterium Substitution Effect on Relaxation Times) method has been utilized for the determination of the syn-anti conformational equilibrium in the monomeric state and for the determination of the mutual orientation of the two adenine rings in the dimeric state of cyclic AMP. The molecules were found to coexist with nearly equimolarity or syn-anti conformers and thermal activation of the molecules perturbs the syn-anti conformational equilibrium to comprise the syn form in preference at higher temperature. The glycosidic isomerization (from anti to syn) was found to be characterized both by a positive enthalpy change and by a positive entropy change. The cyclic AMP molecules prefer to take a 'trans-stacking' conformation in the dimeric state where the two molecules are arranged in such a way that the H(2) of one molecule is close to the H(8) of the other.     

Intramolecular CH···O hydrogen bonds in the AI and BI DNA-like conformers of canonical nucleosides and their Watson-Crick pairs. Quantum chemical and AIM analysis

The aim of this work is to cast some light on the H-bonds in double-stranded DNA in its AI and BI forms. For this purpose, we have performed the MP2 and DFT quantum chemical calculations of the canonical nucleoside conformers, relative to the AI and BI DNA forms, and their Watson-Crick pairs, which were regarded as the simplest models of the double-stranded DNA. Based on the atoms-in-molecules analysis (AIM), five types of the CH···O hydrogen bonds, involving bases and sugar, were detected numerically from 1 to 3 per a conformer: C2'H···O5', C1'H···O2, C6H···O5', C8H···O5', and C6H···O4'. The energy values of H-bonds occupy the range of 2.3-5.6 kcal/mol, surely exceeding the kT value (0.62 kcal/mol). The nucleoside CH···O hydrogen bonds appeared to "survive" turns of bases against the sugar, sometimes in rather large ranges of the angle values, pertinent to certain conformations, which points out to the source of the DNA lability, necessary for the conformational adaptation in processes of its functioning. The calculation of the interactions in the dA·T nucleoside pair gives evidence, that additionally to the N6H···O4 and N1···N3H canonical H-bonds, between the bases adenine and thymine the third one (C2H···O2) is formed, which, though being rather weak (about 1 kcal/mol), satisfies the AIM criteria of H-bonding and may be classified as a true H-bond. The total energy of all the CH···O nontraditional intramolecular H-bonds in DNA nucleoside pairs appeared to be commensurable with the energy of H-bonds between the bases in Watson-Crick pairs, which implies their possible important role in the DNA shaping.     

Conformational analysis of perezone and dihydroperezone using vibrational circular dichroism

The vibrational circular dichroism (VCD) spectra of perezone and dihydroperezone measured from CDCl₃ solutions were quite similar, suggesting analogous conformations for both molecules. Their absolute configurations were confirmed by comparison of the experimental VCD spectrum of each compound with curves generated from theoretical calculations using density functional theory (DFT) at the B3LYP/DGDZVP level of theory taking into account their conformational mobility. Conformational analysis of the 8-(R) enantiomer showed 19 low energy conformers in a 2.4 kcal/mol energy range, while for 8-(R), with the saturated side alkyl chain, 34 conformers were considered in the first 2 kcal/mol. Initial analyses were carried out using a Monte Carlo searching with the MMFF94 molecular mechanics force field, all MMFF94 conformers were geometrically optimized using DFT at the B3LYP/6-31G(d) level of theory, followed by reoptimization and calculations of their vibrational frequencies at the B3LYP/DGDZVP level. Good agreement between the theoretical 8-(R) enantiomers and experimental VCD curves were observed for both.     

Conformational response of tartaric acid to derivatization: role of 1,3-dipole-dipole interactions

The four-carbon chain in (R,R)-tartaric acid derivatives is predominantly antiperiplanar (trans) in the acid, its salts, esters, and NH-amides, while (-)-synclinal (gauche) conformer is the most abundant in N,N'-tetraalkyltartramides. Trialkylsilylation or tert-butylation of the hydroxy groups at C2 and C3 does not appear to affect the conformational preference of NH-tartramides, but it does change the conformational equilibrium in the case of tartrates (toward (-)-gauche) and N,N'-tetraalkyltartramides (toward trans), as judged from the NMR data. X-ray diffraction data point to the stabilizing role of antiparallel dipole-dipole interactions due to the 1,3-CO/CH bonds. These interactions can be found in the trans and (-)-gauche conformers but are not possible for the (+)-gauche conformers of (R,R)-tartaric acid derivatives. This rationalizes small proportion of (+)-gauche conformers in tartaric acid derivatives and points to a significance of 1,3-dipole-dipole interactions. The conformation around the C1-C2 (and C3-C4) bond is different in tartrates (O-C-C=O, syn) and tartramides (O-C-C=O, anti); the CD data (n-pi* band) show that O-silylation or O-tert-butylation brings about conformational changes around the C1-C2 bond in the case of N,N'-tetraalkyldiamides only.     

Effect of carboxylate and Na ions on the tautomeric status of 9-methylguanine

Interactions of 9-methylguanine (m9Gua) with carboxylate ion of acetic acid (CH3COO-) and Na+ were studied by 1H NMR spectroscopy and ab initio quantum chemical calculations of the B3LYP/6-31++G(d,p) and B3LYP/6-311++G(d,p) levels of theory. Changes in the m9Gua 1H NMR spectrum in the presence of the equimolar amount of sodium acetate (NaAc), which in anhydrous DMSO dissociates to CH3COO- and Na+, were interpreted as a consequence of a complex formation of m9Gua in the amino-keto-N1H tautomeric form (m9GuaN1H) with carboxylate ions via two H-bonds involving amino and N1H-imino protons. Quantum chemical calculations of interactions of the m9GuaN1H ground-state tautomer and the m9GuaN3H high energy one with relative energy 20.01 kcal/mol show that the ground state tautomer forms the ground-state complex CH3COO-:mgGuaN1H, by 5.57 kcal/mol more stable than the CH3COO-:m9GuaN3H complex, and coordination of Na+ with the O6 and N7 atoms reduces this energy difference to 2.57 kcal/mol. Such a coordination of Na+ with tautomer m9GuaN3H therewith decreases its relative energy only to 13.31 kcal/mol. Non-additivity of the two ligands contributions to the 8-times reduction of the relative energy of the high energy tautomer in the CH3COO-:m9GuaN3H:Na+(O6,N7) triple complex was concluded, the role of CH3COO- being dominant. Besides, coordination with Na+ resulted in an iminoproton transfer from the base to CH3COO- in the triple complexes of both tautomers, according to calculations in vacuum. Biological significance of the results is noticed.     

An electron spin resonance investigation and molecular orbital calculation of the anion radical intermediate in the enzymatic cis-trans isomerization of furylfuramide, a nitrofuran derivative of ethylene

The enzymatic cis-trans isomerization of nitrofuran derivatives has been proposed to occur via the formation of a radical anion intermediate. ESR investigations, in conjunction with intermediate neglect of differential overlap (INDO) molecular orbital calculations, support this concept by demonstrating the enzymatic generation of cis and trans radical anions of 3-(5-nitro-2-furyl)-2-(2-furyl) acrylamide. The INDO calculations further indicate that the rotational barrier between the cis and trans anion radicals of this compound is only 5--10 kcal/mol, whereas a 70 kcal/mol barrier exists for the parent geometric isomers. Hyperfine splitting constants for the cis-trans conformers have been assigned on the basis of INDO calculations. Surprisingly, only the nitrogen hyperfine splitting of the nitro group is distinguishably different in the two conformers, a result which is not inconsistent with the INDO calculations.     

The anomeric effect revisited. A possible role of the CH/n hydrogen bond

Ab initio MO calculations were carried out at the MP2/6-311++G(d,p) level to investigate the conformational energy of 2-substituted oxanes and 1,3-dioxanes. It has been found that the Gibbs free energies of the axial conformers are smaller than those of the corresponding equatorial conformers in every case when the 2-substituent Z is electron withdrawing (OCH(3), F, Cl, Br). The difference in Gibbs energy between the equatorial and axial conformers DeltaG(eq-ax) increases from Z=OCH(3) to F, Cl, and then to Br. In the axial conformers, the interatomic distance between Z and the axial C-H, separated by four covalent bonds, has been found to be appreciably shorter than the van der Waals distance, suggesting the importance of the five-membered CH/n (CH/O or CH/halogen) hydrogen bond in stabilizing these conformations. Natural bonding orbital (NBO) charges of the relevant atoms have been shown to be different between the two conformers: more positive for H and more negative for C in the axial conformers than in the corresponding equatorial conformers. In view of the above findings, we suggest that the CH/n hydrogen bond plays an important role in stabilizing the axial conformation in 2-substituted oxanes and 1,3-dioxanes, and by implication, in the anomeric effect in carbohydrate chemistry.     

The Raman spectra of left-handed DNA oligomers incorporating adenine-thymine base pairs+.

Raman spectra were obtained from single crystals of [d(CGCATGCG)]2 and [d(m5CGTAm5CG)]2, both of which incorporate A-T pairs into Z-DNA structures and contain C2'-endo/syn conformers of deoxyguanosine at the oligonucleotide ends. Correlation with x-ray results permits the following Raman assignments for nucleoside conformers: C3'-endo/syn G, 623 +/- 1; C2'-endo/syn G, 671 +/- 2; C2'-endo/anti C, 782 +/- 1; C2'endo/anti T, 650 +/- 5 and ca. 750; C3'-endo/syn A, 729 +/- 1 cm-1. These results show that (i) the 670 cm-1 line of syn G is highly sensitive to the change from C3'-endo to C2'-endo pucker, (ii) the 729 cm-1 line of A is affected neither by furanose pucker nor glycosidic bond orientation and (iii) the 1200-1500 cm-1 region of the Raman spectrum of the A-T double helix is greatly altered by the B-to-Z transition. Conformation sensitive Raman frequencies in the 850-1700 cm-1 region are identified for both octamer and hexamer, and the Z-to-B transition of each is monitored by spectral changes which occur upon dissolving the crystal in H2O solution.     

Conformational stability and normal coordinate analyses for 1-halovinyl azides CH<Subscript>2</Subscript>=CX–NNN (X is F,Cl and Br)

The conformational behavior of 1-halovinyl azides CH2=CX-NNN (X=F, Cl and Br) were investigated by DFT-B3LYP and ab initio MP2 calculations with the 6-311++G** basis set. The molecules were predicted to exist predominantly in the trans (the vinyl CH2=CH- and the azide -NNN groups are trans to each other) conformation. The relative energy between cis and trans were calculated to decrease in order: bromide>chloride>fluoride. Full optimization was performed at the ground and transition states in the molecule at both MP2 and B3LYP levels. The barrier to internal rotation around the C-N single bond in the three molecules was calculated to be about 4-5 kcal mol(-1). The vibrational frequencies were computed at the DFT-B3LYP level and the calculated infrared and Raman spectra of the cis- trans mixture of the three molecules were plotted. Complete vibrational assignments were made on the basis of normal coordinate calculations for both stable conformers of the three molecules.     

Conformational analysis of canonical 2-deoxyribonucleotides. 2. Purine nucleotides

The molecular structure of different conformers of isolated canonical purine 2'-deoxyribonucleotides 2-deoxyadenosine-5'-phosphate (pdA) and 2'-deoxyguanosine-5'-phosphate (pdG) was optimized using the B3LYP/6-31G(d) method. The results of the calculations reveal that the geometrical parameters and relative stability of the conformers significantly depend on the nature of the nucleobase, its orientation, the conformation of the furanose ring, the charge of the phosphate group and the character of the intramolecular hydrogen bonds. Analysis of the electron density distribution in purine nucleotides reveals the existence of a number of intramolecular hydrogen bonds. In general, the south conformer has a lower energy at the anti orientation of the base, and both conformers occur as the most stable for the syn orientation of the nucleobases. The results of the calculations reveal that the geometry and relative energy of the conformers of purine DNTs may be easily tuned by the charge of the phosphate group.     

Calculation of conformational energies and optical rotation of the most simple chiral alkane

Quantum chemical calculations have been performed to investigate the conformer distribution of 4-ethyl-4-methyloctane and its optical rotation. With the reference methods MP2 and SCS-MP2, the energies of seven conformers are found within a range of about 1.5 kcal mol(-1). It is demonstrated that the relative energies cannot be reliably predicted with conventional GGA or hybrid density functionals, Hartree-Fock, semiempirical AM1, and classical force field (MM3) calculations. An empirical dispersion correction to GGA (PBE-D), hybrid (B3LYP-D), or double hybrid (B2PLYP-D) functionals corrects these errors and results in very good agreement with the reference energies. Optical rotations have been calculated for all seven conformers at the TDDFT(BHLYP/aTZV2P) level. The computed macroscopic rotation is derived from a classical Boltzmann average. The result (1.9-3.2 deg dm(-1) (g/mL)(-1)) is very close to the experimental value of 0.2 deg dm(-1) (g/mL)(-1) for the (R)-enantiomer and has the right sign. Because six conformers are significantly populated at room temperature and the rotations of individual conformers differ in sign and magnitude, the calculated average rotation is rather sensitive to the conformer population used. From the electronic structure point of view, this example emphasizes the effect of long-range dispersion effects for the evaluation of population averaged quantities in large molecules. Computations based on free enthalpies are in worse agreement with experiment that is attributed to artefacts of the harmonic approximation used to compute the vibrational entropy terms.     

Structural properties of the [d(G3T4G3)]2 quadruplex: evidence for sequential syn-syn deoxyguanosines.

Two-dimensional 1H NMR studies on the dimeric hairpin quadruplex formed by d(G3T4G3) in the presence of either NaCl or KCl are presented. In the presence of either salt, the quadruplex structure is characterized by half the guanine nucleosides in the syn conformation about the glycosidic bond, the other half in the anti conformation, as reported for other similar sequences. However, 1H NOESY and 1H-31P heteronuclear correlation experiments demonstrate that the deoxyguanosines do not strictly alternate between syn and anti along individual strands. Thus we find the following sequences with regard to glycosidic bond conformation: 5'-G1SG2SG3AT4AT5A-T6AT7AG8SG9AG10A-3' and 5'-G11SG12AG13AT14AT1 5AT16AT17AG18SG19SG20A-3', where S and A denote syn and anti, respectively. This represents the first experimental evidence for a nucleic acid structure containing two sequential nucleosides in the syn conformation. The stacking interactions of the resulting quadruplex quartets and their component bases have been evaluated using unrestrained molecular dynamics calculations and energy component analysis. These calculations suggest that the sequential syn-syn/anti-anti and syn-anti quartet stacks are almost equal in energy, whereas the anti-syn stack, which is not present in our structure, is energetically less favorable by about 4 kcal/mol. Possible reasons for this energy difference and its implications for the stability of quadruplex structures are discussed.     

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.     

Intrinsic conformational energetics associated with the glycosyl torsion in DNA: a quantum mechanical study

The glycosyl torsion (chi) in nucleic acids has long been recognized to be a major determinant of their conformational properties. chi torsional energetics were systematically mapped in deoxyribonucleosides using high-level quantum mechanical methods, for north and south sugar puckers and with gamma in the g(+) and trans conformations. In all cases, the syn conformation is found higher in energy than the anti. When gamma is changed from g(+) to trans, the anti orientation of the base is strongly destabilized, and the energy difference and barrier between anti and syn are significantly decreased. The barrier between anti and syn in deoxyribonucleosides is found to be less than 10 kcal/mol and tends to be lower with purines than with pyrimidines. With gamma = g(+)/chi = anti, a south sugar yields a significantly broader energy well than a north sugar with no energy barrier between chi values typical of A or B DNA. Contrary to the prevailing view, the syn orientation is not more stable with south puckers than with north puckers. The syn conformation is significantly more energetically accessible with guanine than with adenine in 5-nucleotides but not in nucleosides. Analysis of nucleic acid crystal structures shows that gamma = trans/chi = anti is a minor but not negligible conformation. Overall, chi appears to be a very malleable structural parameter with the experimental chi distributions reflecting, to a large extent, the associated intrinsic torsional energetics.     

Kinetics and thermodynamics of the interaction of 5-fluoro-2'-deoxyuridylate with thymidylate synthase

Thymidylate synthase (TS), 5-fluorodeoxyuridylate (FdUMP), and 5,10-methylenetetrahydrofolate (CH2-H4folate) form a covalent complex in which a Cys thiol of TS is attached to the 6-position of FdUMP and the one-carbon unit of the cofactor is attached to the 5-position. The kinetics of formation of this covalent complex have been determined at several temperatures by semirapid quench methods. Together with previously reported data the results permit calculation of every rate and equilibrium constant in the interaction. Conversion of the noncovalent ternary complex to the corresponding covalent complex proceeds at a rate of 0.6 s-1 at 25 degrees C, and the dissociation constant for loss of CH2-H4folate from the noncovalent ternary complex is approximately 1 microM. Activation parameters for the formation of the covalent complex were shown to be Ea = 20 kcal/mol, delta G+ = 17.9 kcal/mol, delta H+ = 19.3 kcal/mol, and delta S+ = 0.005 kcal/(mol.deg). The equilibrium constant between the noncovalent and covalent ternary complexes is approximately 2 X 10(4), and the overall dissociation constant of CH2-H4folate from the covalent complex is approximately 10(-11) M. The conversion of the noncovalent ternary complex to the covalent adduct is about 12-fold slower than kcat in the normal enzymic reaction. However, because the dissociation constant for CH2-H4folate from the noncovalent ternary complex is about 10-fold lower than that from the TS-dUMP-CH2-H4folate Michaelis complex, the terms corresponding to kcat/Km are nearly equal. We propose that some of the intrinsic binding energy of CH2-H4folate may be used to facilitate formation of a 5-iminium ion intermediate.(ABSTRACT TRUNCATED AT 250 WORDS)     

A density functional investigation of the extradiol cleavage mechanism in non-heme iron catechol dioxygenases

The mechanism for extradiol cleavage in non-heme iron catechol dioxygenase was modelled theoretically via density functional theory. Based on the Fe(II)-His,His,Glu motif observed in enzymes, an active site model complex, [Fe(acetate)(imidazole)(2)(catecholate)(O(2))](-), was optimized for states with six, four and two unpaired electrons (U6, U4 and U2, respectively). The transfer of the terminal atom of the coordinated dioxygen leading to "ferryl" Fe=O intermediates spontaneously generates an extradiol epoxide. The computed barriers range from 19 kcal mol(-1) on the U6 surface to approximately 25 kcal mol(-1) on the U4 surface, with overall reaction energies of +11.6, 6.3 and 7.1 kcal mol(-1) for U6, U4 and U2, respectively. The calculations for a protonated process reveal the terminal oxygen of O(2) to be the thermodynamically favoured site but subsequent oxygen transfer to the catechol has a barrier of approximately 30-40 kcal mol(-1), depending on the spin state. Instead, protonating the acetate group gives a slightly higher energy species but a subsequent barrier on the U4 surface of only 7 kcal mol(-1) relative to the hydroperoxide complex. The overall exoergicity increases to 13 kcal mol(-1). The favoured proton-assisted pathway does not involve significant radical character and has features reminiscent of a Criegee rearrangement which involves the participation of the aromatic ring pi-orbitals in the formation of the new carbon-oxygen bond. The subsequent collapse of the epoxide, attack by the coordinated hydroxide and final product formation proceeds with an overall exoergicity of approximately 75 kcal mol(-1) on the U4 surface.     

C-H ... O hydrogen bonding in solutions of methylated nucleic acid base analogs as revealed by NMR

Formation and thermodynamic characteristics of C-H ... O hydrogen bonding of methylated uracils and caffeine have been studied by nmr along two lines. 1. The concentration and temperature dependencies of the PMR spectra of 1,3-dimethyluracil (m2 1,3Ura), 1,3-dimethylthymine (m2 1,3Thy), and 1,3,6-trimethyluracil (m3 1,3,6Ura) in chloroform at high concentrations of base analogs indicated the self-association of m2 1,3Ura and m2 1,3Thy via C(6)H ... O hydrogen bonding and the competitive formation of C-H ... O bonds between carbonyl oxygens and chloroform. The intermolecular interaction energy and the arrangement of molecules in the local minima of various m2 1,3Ura dimers were calculated by the method of atom-atom potentials. The deepest minimum for the m2 1,3Ura coplanar dimer corresponds to a C(6)-H ... O hydrogen-bond formation. 2. At low concentration of m2 1,3Ura and caffeine in CCl4, C(6)-H ... O bonding for m2 1,3Ura and C(8)-H ... O bonding for caffeine with oxygens of dimethyl sulfoxide (DMSO) and acetone were observed. The association constants of these complexes were obtained at different temperatures. The enthalpies delta H, of the m2 1,3Ura-DMSO, m2 1,3Ura-accetone, caffeine-DMSO, and caffeine-acetone complexes were -2 +/- 0.1 kcal/mol. The calculations showed that the deepest minimum of the caffeine-acetone coplanar complex corresponds to C(8)-H ... O bonding with energy of -3.5 kcal/mol and that of the m2 1,3Ura-acetone complexes corresponds to C(6)-H ... O bonding with energy of -3.4 kcal/mol. The approximate correction for the solvent effect provides good agreement of the experimental data with the calculations.     

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.     

Light induces destabilization of photoactive yellow protein

To understand the effect of visible light on the stability of photoactive yellow protein (PYP), urea denaturation experiments were performed with PYP in the dark and with PYP(M) under continuous illumination. The urea concentrations at the midpoint of denaturation were 5.26 +/- 0.29 and 3.77 +/- 0.19 M for PYP and PYP(M), respectively, in 100 mM acetate buffer, and 5.26 +/- 0.24 and 4.11 +/- 0.12 M for PYP and PYP(M), respectively, in 100 mM citrate buffer. The free energy change upon denaturation (DeltaG(D)(H2O)), obtained from the denaturation curve, was 11.0 +/- 0.4 and 7.6 +/- 0.2 kcal/mol for PYP and PYP(M), respectively, in acetate buffer, and 11.5 +/- 0.3 and 7.8 +/- 0.1 kcal/mol for PYP and PYP(M), respectively, in citrate buffer. Even though the DeltaG(D)(H2O) value for PYP(M) is almost identical in the two buffer systems, the urea concentration at the midpoint of denaturation is lower in acetate buffer than in citrate buffer. Although their CD spectra indicate that the protein conformations of the denatured states of PYP and PYP(M) are indistinguishable, the configurations of the chromophores in their denatured structures are not necessarily identical. Both denatured states are interconvertible through PYP and PYP(M). Therefore, the free energy difference between PYP and PYP(M) is 3.4-3.7 kcal/mol for the protein moiety, plus the additional contribution from the difference in configuration of the chromophore.