A computer-aided quantum chemical study of the N<Stack><Subscript>15</Subscript><Superscript>−</Superscript></Stack> cluster

Ab initio (RHF, MP2) and Density Functional Theory (DFT) methods have been used to examine six isomers of the N15m cluster with the 6-31+G* basis set. Different from the known odd-numbered anionic N7m, N9m, and N11m clusters, in which the open-chain structures are the most stable species, the most stable N15m isomer is structure 1 (C1), which may be considered as a complex between the fragments cyclic N5m (D5h) and staggered N10 (D2d). The decomposition pathways of structure 2 (CS), containing two aromatic N5 rings connected by a N5 chain, and the open-chain structure 3 (C2v) were studied at the B3LYP/6-31+G* level of theory. Relative energies were refined at the level of B3LYP/6-311+G(3df,2p)//B3LYP/6-31+G*+ZPE (B3LYP/6-31+G*). The barriers for N2 and N5m (D5h) fission reactions for structure 2 are predicted to be 18.2 and 14.2 kcal x mol(-1), respectively. The corresponding N2+N3m fission barrier for structure 3 is predicted to be 11.2 kcal x mol(-1). Supplementary material is available for this article if you access the article at http://dx.doi.org/10.1007/s00894-003-0118-0. A link in the frame on the left on that page takes you directly to the supplementary material. Figure Structure 1 of the N15m cluster, showing bond distances in A and bond angles in degrees     

QM/MM study of the active site of free papain and of the NMA-papain complex.

Hybrid quantum mechanical/molecular mechanical (QM/MM) calculations using restricted and unrestricted Hartree-Fock and B3LYP ab initio (QM) and Amber force field (MM), respectively, have been applied to study the catalytic site of papain in both free and substrate bonded forms. Ab initio geometry optimizations have been performed for the active site of papain and the N-methyl-acetamide (NMA)-papain complex within the molecular mechanical treatment of the protein environment. A covalent tetrahedral intermediate structure could be obtained only when the amide N atom of the substrate molecule was protonated through a proton transfer from the His-159 in the catalytic site. Our results support the previous assumption that a proton transfer from His-159 to the amide N atom of the substrate occurs prior to or concerted with the nucleophilic attack of the Cys-25 sulfur atom to the carbonyl group of the substrate. The electron correlation effect will reduce the proton transfer barrier. Therefore, this proton transfer can be easily observed in the B3LYP/6-31G* calculations. The HF/6-31G* method overestimates the reaction barrier against this proton transfer. The sulfur atom of Cys-25 and the imidazole ring of His-159 are found to be coplanar in the free form of the enzyme. However, the rotation of the imidazole ring of His-159 was observed during the formation of the tetrahedral intermediate. Without the papain environment, the coplanar thiolate-imidazolium ion pair RS-...ImH+ is much less stable than the neutral form of RSH....Im. Within the protein environment, however, the thiolate-imidazolium ion pair becomes more stable than its neutral form by 4.1 and 0.4 kcal/mol in HF/6-31G* and B3LYP/6-31G* calculations, respectively. The barrier of proton transfer from S-H group of Cys-25 to the imidazole ring of His-159 was reduced from 22.0 kcal/mol to 15.2 kcal/mol by the protein environment in HF/6-31G* calculations. This barrier is found to be much smaller (2.5 kcal/mol) in B3LYP/6-31G* calculations.     

Theoretical Investigation of the Molecular Structure of the πκ DNA Base Pair

Abstract

The structure of the nonclassical πκ base pair (7–methyl-oxoformycin … 2,4-diaminopyrimidine) was studied at the ab initio Hartree-Fock (HF) and MP2 levels using the 6–31G* and 6–31G** basis sets. The πκ base pair is bound by three parallel hydrogen bonds with the donor-acceptor-donor recognition pattern. Recently, these bases were proposed as an extension of the genetic alphabet from four to six letters (Piccirilli et al. Nature 343, 33(1990)). By the HF/6- 31G* method with full geometry optimization we calculated the 12 degree propeller twist for the minimum energy structure of this complex. The linearity of hydrogen bonds is preserved in the twisted structure by virtue of the pyramidal arrangement of the κ-base amino groups. The rings of both the π and κ molecules remain nearly planar. This nonplanar structure of the πκ base pair is only 0.1 kcal/mol more stable than the planar (Cs) conformation. The HF/6- 31G* level gas-phase interaction energy of πκ (—13.5 kcal/mol) calculated by us turned out to be nearly the same as the interaction energy obtained previously for the adenine-thymine base pair (—13.4 kcal/mol) at the same computational level. The inclusion of p-polarization functions on hydrogens, electron correlation effects (MP2/6–31G** level), and the correction for the basis set superposition error (BSSE) increase this energy to -14.0 kcal/mol.     

Theoretical investigations of the hydrolysis pathway of verdoheme to biliverdin

Conversion of iron(II) verdoheme to iron biliverdin in the presence of OH(-) was investigated using B3LYP method. Both 3-21G and 6-31G* basis sets were employed for geometry optimization calculation as well as energy stabilization estimation. Calculation at 6-31G* level was found necessary for a correct spin state estimation of the iron complexes. Two possible pathways for the conversion of iron verdoheme to iron biliverdin were considered. In one path the iron was six-coordinate while in the other it was considered to be five-coordinate. In the six-coordinated pathway, the ground state of bis imidazole iron verdoheme is singlet while that for open chain iron biliverdin it is triplet state with 4.86 kcal/mol more stable than the singlet state. The potential energy surface suggests that a spin inversion take place during the course of reaction after TS. The ring opening process in the six-coordinated pathway is in overall -2.26 kcal/mol exothermic with a kinetic barrier of 9.76 kcal/mol. In the five-coordinated pathway the reactant and product are in the ground triplet state. In this path, hydroxyl ion attacks the iron center to produce a complex, which is only 1.59 kcal/mol more stable than when OH(-) directly attacks the macrocycle. The activation barrier for the conversion of iron hydroxy species to the iron biliverdin complex by a rebound mechanism is estimated to be 32.68 kcal/mol. Large barrier for rebound mechanism, small barrier of 4.18 kcal/mol for ring opening process of the hydroxylated macrocycle, and relatively same stabilities for complexes resulted by the attack of nucleophile to the iron and macrocycle indicate that five-coordinated pathway with direct attack of nucleophile to the 5-oxo position of macrocycle might be the path for the conversion of verdoheme to biliverdin.     

The investigation of hydrocarbon cracking reaction energetics with composite energy methods

Hydrocarbo cracking reactions are one of the most commonly encountered reactions in the petroleum industry, and the energetics of the reactions are crucial in understanding the reaction mechanisms and predicting reaction rates. In this work, a modified composite energy method (CBS-RAD(MP2)) is created as a version of the CBS-RAD method which gives accurate energetics for hydrocarbon free radical reactions. It replaces the time consuming QCISD(fc)/6-31g* method in the geometry optimization and frequency calculation steps with MP2(full)/6-31g* level calculations. The accuracy of the new CBS-RAD(MP2) method is compared with the widely used G2, G3 and CBS-QB3 composite methods for predicting heats of reaction and activation barriers of 14 hydrocarbon cracking reactions. We find that the new CBS-RAD(MP2) method has the second least RMS error of 1.22 kcal/mol for heats of reaction calculations. For activation energy calculations, the new CBS-RAD(MP2) method has the least RMS error of 1.37 kcal/mol. Moreover, the CBS-RAD(MP2) method was found to require only 81% of the computational time required compared to the CBS-QB3 method, 32% of G3 and 15% of the G2 method, making it an attractive alternative for predicting hydrocarbon cracking reaction energetics.     

Study of the hydrogen bond in different orientations of adenine-thymine base pairs: An <Emphasis Type="Italic">ab initio</Emphasis> study

In order to gain deeper insight into structure, charge distribution, and energies of A-T base pairs, we have performed quantum chemical ab initio and density functional calculations at the HF (Hartree-Fock) and B3LYP levels with 3-21G*, 6-31G*, 6-31G**, and 6-31++G** basis sets. The calculated donor-acceptor atom distances in the Watson-Crick A-T base pair are in good agreement with the experimental mean values obtained from an analysis of 21 high resolution DNA structures. In addition, for further correction of interaction energies between adenine and thymine, the basis set superposition error (BSSE) associated with the hydrogen bond energy has been computed via the counterpoise method using the individual bases as fragments. In the Watson-Crick A-T base pair there is a good agreement between theory and experimental results. The distances for (N2...H23-N19), (N8-H13...O24), and (C1...O18) are 2.84, 2.94, and 3.63 A, respectively, at B3LYP/6-31G** level, which is in good agreement with experimental results (2.82, 2.98, and 3.52 A). Interaction energy of the Watson-Crick A-T base pair is -13.90 and -10.24 kcal/mol at B3LYP/6-31G** and HF/6-31G** levels, respectively. The interaction energy of model (9) A-T base pair is larger than others, -18.28 and -17.26 kcal/mol, and for model (2) is the smallest value, -13.53 and -13.03 kcal/mol, at B3LYP/6-31G** and B3LYP/6-31++G** levels, respectively. The computed B3LYP/6-31G** bond enthalpies for Watson-Crick A-T pairs of -14.4 kcal/mol agree well with the experimental results of -12.1 kcal/mol deviating by as little as -2.3 kcal/mol. The BSSE of some cases is large (9.85 kcal/mol) and some is quite small (0.6 kcal/mol).     

Ab initio Models for Six-Center Multiple Proton Exchange and Ion Pair Formation Assisted by Lewis Acids

High level ab initio and density functional calculations, extrapolated to QCISD(T)/6-311+G(3df,2p)//MP2/6-31+G^**+ZPE, reveal that cyclic ion pairs can form in the hydrogen bonded complexes of haloboric acids BH_nX_3-n–HX, X=F, Cl, with Lewis bases HX, H₂O, CH₃OH, and NH₃, even in isolation (e.g., in the gas phase). The intrinsic acidities (deprotonation energies) required for protonation of these bases with formation of gas phase ion pairs are calculated to be <295 kcal/mol for water, <301 kcal/mol for methanol, and <306 kcal/mol for ammonia; such values are common for acidic sites in zeolites. All gas phase ion pairs prefer symmetric bidentate or tridentate structures. In the other cases where hydrogen bonded complexes prevail, symmetric ion pair-like transition structures for multiple hydrogen exchange are computed.Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1007/s0089400060563     

Quantum-chemical investigation of tautomerization ways of Watson-Crick DNA base pair guanine-cytosine

A novel physico-chemical mechanism of the Watson-Crick DNA base pair Gua.Cyt tautomerization Gua.Cyt*<---->Gua.Cyt<---->Gua*.Cyt (mutagenic tautomers of bases are marked by asterisks) have been revealed and realized in a pathway of single proton transfer through two mutual isoenergetic transition states with Gibbs free energy of activation 30.4 and 30.6 kcal/mol and they are ion pairs stabilized by three (N2H...N3, N1H...N4- and O6+H...N4-) and five (N2H...O2, N1H...O2, N1H...N3, O6+H...N4- and 06+H...N4-) H-bonds accordingly. Stable base pairs Gua-Cyt* and Gua*.Cyt which dissociate comparably easy into monomers have acceptable relative Gibbs energies--12.9 and 14.3 kcal/mol--for the explanation of the nature of the spontaneous transitions of DNA replication. Results are obtained at the MP2/6-311++G(2df,pd)//B3LYP/6-31 1++G(d,p) level of theory in vacuum approach.     

Theoretical investigation of (-)1-(benzofuran-2-yl)-2-propylaminopentane[(-)-BPAP] as a hydroxyl radical scavenger

The chemical reactions between (-)-BPAP and .OH were studied using molecular orbital theory, with several simplified models. (-)-BPAP was proved to be a good radical scavenger. It was found that every atom of the benzofuran ring, except carbon 3, was capable of easily trapping the radical, where the most active site was carbon 1 on the furan part. The activation energies for the trappings were ca. 10kcal/mol by the calculations at UHF/6-31G(d)//UHF/3-21G level. Since the single radical trapping products were still radicals, these could trap further radicals by way of cascade without any activation energy. Thus, the double radical trapping products were very stable, of which the lowest energy level was about 80kcal/mol lower than the starting reactants.     

Examining the reaction between antioxidant compounds and 2,2-diphenyl-1-picrylhydrazyl (DPPH) through a computational investigation

In this work, we present a computational investigation on the reactions between two well-known antioxidants (quercetin and morin) and 2,2-diphenyl-1-picrylhydrazyl (DPPH). A density functional theory (DFT) approach with the B3LYP functional and the 6-31G(d,p) basis set was used for the simulations. The structural and energetic parameters (Gibbs free-energy, ΔG, and Gibbs free-energy of activation, ΔG⁺⁺) were determined to provide information on the antioxidant activity as well as to evaluate the contributions of each hydroxyl group to the referred property. According to the results obtained, quercetin presented three hydroxyls as being thermodynamically spontaneous in the reaction with DPPH (4\(^{\prime }\)-ArOH, 3\(^{\prime }\)-ArOH, and 3-ArOH, with ΔG = -4.93 kcal/mol, -2.89 kcal/mol, and -1.87 kcal/mol, respectively) against only two in the case of morin (2\(^{\prime }\)-ArOH and 3-ArOH, with ΔG = -7.56 kcal/mol and -4.57 kcal/mol, respectively). Hence, quercetin was found to be a more efficient antioxidant, which is in agreement with different experimental and computational investigations of bond dissociation enthalpies (BDEs). However, the order of contribution of the OH groups of each compound to the antioxidant potential present some differences when compared to what was seen in the previous investigations, especially for morin. These findings are in contrast to what was observed in studies based on the determinations of BDEs. Therefore, experimental investigations on the hydrogen-atom transfer mechanism (HAT) for both compounds are encouraged in order to clarify these observations.     

Hydration of carbon dioxide by carbonic anhydrase: internal proton transfer of Zn2+-bound HCO3-

Proton transfer within HCO3- has been examined under various conditions through molecular orbital methods: partial retention of diatomic differential overlap and 4-31G self-consistent field programs. These conditions include the absence or presence of Zn2+, Zn2+(NH3)3, or a water ligand on Zn2+. In addition, 4-31G+ and some MP2/4-31G results are obtained. The use of Be2+ to simulate Zn2+ reproduces reaction pathways and energy barriers, except for marginal cases. The barrier of 35.6 kcal/mol for direct internal proton transfer is reduced to 3.5 kcal/mol when one water molecule, not bound to Zn2+, is included for proton relay and to 1.4 kcal/mol when two such water molecules are included. In the enzyme, either Thr-199 or solvent molecules could perform this relay function. Our results favor this facilitated proton transfer over a mechanism in which Zn2+-bound OH- attacks CO2, a bidentate intermediate forms, and the OH moiety of the resulting HCO3- dissociates from Zn2+, thus leaving one of the oxygens of the original CO2 as a ligand to Zn2+.     

Ab initio and DFT study on the electrophilic addition reaction of bromine to tetracyclo[5.3.0.02,6.03,10]deca–4,8–diene

The electronic and geometric structures of tetracyclo[5.3.0.0^2,6.0^3,10]deca-4,8-diene (hypostrophene) have been investigated by ab initio and DFT/B3LYP methods using the 6-31G* and 6-311G* basis sets. The double bonds of hypostrophene are endo-pyramidalized. The cationic intermediates and products formed in the addition reaction have been investigated using the HF/6-311G*, HF/6-311G**, and B3LYP/6-311G* methods. The bridged bromonium cation was more stable than the U-type cation. Considering that the bridged cation does not isomerize to the less stable U-type cation, it is not possible for the U-type product to be obtained in the reaction. The bridged bromonium cation transformed into the more stable N-type cation and the N-type product was obtained via this cation. The thermodynamic stability of the exo, exo and exo, endo isomers of the N-type dibromide molecule were almost identical. The N-type product was 16.6 kcal mol⁻¹ more stable than the U-type product. Figure General energy diagram of the hypostrophene–bromine (HS–Br₂) system (kcal mol⁻¹) (MP2/6-311G*//HF/6-311G*)     

Exploiting diverse stereochemistry of β-amino acids: toward a rational design of sheet-forming β-peptide systems

Due to the two methylene groups in their backbone, β-amino acids can adopt numerous secondary structures, including helices, sheets and nanotubes. Chirality introduced by the additional side chains can significantly influence the folding preference of β-peptides composed of chiral β-amino acids. However, only conceptual suggestions are present in the literature about the effect of chirality on folding preferences. Summarizing both the experimental and computational results, Seebach (Chem Biodivers 1:1111-1240, 2004) has proposed the first selection rule on the effect of side chain chirality, on the folding preference of β-peptides. In order to extend and fine-tune the aforementioned predictions of Seebach, we have investigated its validity to the novel type of apolar sheet proposed recently (Pohl et al. in J Phys Chem B 114:9338-9348, 2010). In order to facilitate the rational design of sheet-like structures, a systematic study on the effect of chirality on "apolar" sheet stability is presented on disubstituted [HCO-β-Ala-β(2,3)-hAla-β-Ala-NH(2)](2) model peptides calculated at the M05-2X/6-311++G(d,p)//M05-2X/6-31G(d) and B3LYP/6-311++G(d,p)//B3LYP/6-31G(d) levels of theory both in vacuum and in polar and apolar solvents. In addition, both types of "apolar" sheets were investigated; the one with two strands of identical (AA) and enantiomeric (AB) backbone structure. Our results show that heterochirally disubstituted sheets have the greatest preference for sheet formation (ΔG ~ -11 kcal mol(-1)). However, in contrast to Seebach's predictions, "homochiral disubstitution" itself does not necessarily disrupt the sheet structure, rather it could result stable fold (ΔG ~ -5 kcal mol(-1)). Results indicate that both the methyl group orientation and the local conformational effect of substitution affects sheet stability, as point chirality was found to have influence only on the backbone torsional angles. These results enabled us to extend and generalize Seebach's predictions and to propose a more general and accurate "rule of thumb" describing the effect of chirality on sheet stability. This offers an easy-to-use summary on how to design β-peptide sheet structures. We conclude that heterochirally disubstituted models are the best candidates for sheet formation, if the two strands are substituted in a way to create identical torsional angle sets on the two backbones for ideal hydrogen-bonding pattern. With adequately selected side chains, homochirally disubtituted derivatives may also form sheet structures, and the position of methyl groups would prevent assembly of more than two strands making it ideal to create hairpins.     

Energy barriers and rates of tautomeric transitions in DNA bases: ab initio quantum chemical study

Tautomeric transitions of DNA bases are proton transfer reactions, which are important in biology. These reactions are involved in spontaneous point mutations of the genetic material. In the present study, intrinsic reaction coordinates (IRC) analyses through ab initio quantum chemical calculations have been carried out for the individual DNA bases A, T, G, C and also A:T and G:C base pairs to estimate the kinetic and thermodynamic barriers using MP2/6-31G** method for tautomeric transitions. Relatively higher values of kinetic barriers (about 50-60 kcal/mol) have been observed for the single bases, indicating that tautomeric alterations of isolated single bases are quite unlikely. On the other hand, relatively lower values of the kinetic barriers (about 20-25 kcal/mol) for the DNA base pairs A:T and G:C clearly suggest that the tautomeric shifts are much more favorable in DNA base pairs than in isolated single bases. The unusual base pairing A':C, T':G, C':A or G':T in the daughter DNA molecule, resulting from a parent DNA molecule with tautomeric shifts, is found to be stable enough to result in a mutation. The transition rate constants for the single DNA bases in addition to the base pairs are also calculated by computing the free energy differences between the transition states and the reactants.     

DFT study of nano zinc/copper voltaic cells

To facilitate the development of new materials for use in batteries, it is necessary to develop ab initio full-electron computational techniques for modeling potential new battery materials. Here, we tested density functional theory procedures that are accurate enough to obtain the energetics of a zinc/copper voltaic cell. We found the magnitude of the zero-point energy correction to be 0.01–0.2 kcal/mol per atom or molecule and the magnitude of the dispersion correction to be 0.1–0.6 kcal/mol per atom or molecule for Zn _n , (H₂O) _n , \( \mathrm{Zn}{\left({\mathrm{H}}_2\mathrm{O}\right)}_n^{2+} \), \( \mathrm{Cu}{\left({\mathrm{H}}_2\mathrm{O}\right)}_n^{2+} \), and Cu _n . Counterpoise correction significantly affected the values of ?\( {E}_n^{\mathrm{abs}} \), ?\( {E}_n^{\mathrm{coh}} \), and ?E^solv by 1.0–3.1 kcal/mol per atom or molecule at the B3PW91/6-31G(d) level of theory, but by only 0.04–0.4 kcal/mol per atom or molecule at the B3PW91/cc-pVTZ level of theory. The application of B3PW91/6-31G(d) yielded results that differed from macroscopic experimental values by 0.1–7.1 kcal/mol per atom or molecule, whereas applying B3PW91/cc-pVTZ produced results that differed from macroscopic experimental values by 0.1–4.8 kcal/mol per atom or molecule, with the smallest differences occurring for reactions with a small macroscopic experimental ?E and the largest differences occurring for reactions with a large macroscopic experimental ?E, implying size consistency.     

Probing the antioxidant potential of phloretin and phlorizin through a computational investigation

The structures and energetics of two dihydrochalcones (phloretin and its glycoside phlorizin) were examined with density functional theory, using the B3LYP, M06-2X, and LC-ω PBE functionals with both the 6-311G(d,p) and 6-311 + G(d,p) basis sets. Properties connected to antioxidant activity, i.e., bond dissociation enthalpies (BDEs) for OH groups and ionization potentials (IPs), were computed in a variety of environments including the gas-phase, n-hexane, ethanol, methanol, and water. The smallest BDEs among the four OH groups for phloretin (three for phlorizin) were determined (using B3LYP/6-311 + G(d,p) in water) to be 79.36 kcal/mol for phloretin and 79.98 kcal/mol for phlorizin while the IPs (at the same level of theory) were obtained as 139.48 and 138.98 kcal/mol, respectively. By comparing with known antioxidants, these values for the BDEs indicate both phloretin and phlorizin show promise for antioxidant activity. In addition, the presence of the sugar moiety has a moderate (0-6 kcal/mol depending on functional) effect on the BDEs for all OH groups. Interestingly, the BDEs suggest that (depending on the functional chosen) the sugar moiety can lead to an increase, decrease, or no change in the antioxidant activity. Therefore, further experimental tests are encouraged to understand the substituent effect on the BDEs for phloretin and to help determine the most appropriate functional to probe BDEs for dihydrochalcones.     

Racemization of the gastrointestinal antisecretory chiral drug esomeprazole magnesium via the pyramidal inversion mechanism: A theoretical study

The pyramidal inversion mechanisms of the 6‐methoxy and the 5‐methoxy tautomers of (S)‐omeprazole were studied, employing ab initio and DFT methods. The conformational space of the model molecule (S)‐2‐[(3‐methyl‐2‐pyridinyl)methyl]sulfinyl‐1H‐benzimidazole was calculated, with respect to rotations around single bonds, at the B3LYP/6‐311G(d,p) level. All of the resulting conformations were used as starting points for full optimizations of (S)‐omeprazole, at B3LYP/6‐31G(d), B3LYP/6‐311G(d,p), B3LYP/6‐311++G(d,p), B3LYP/6‐311G(2df,2pd), MP2/6‐31G(d), and MP2/6‐311G(d,p) levels. Four distinct pathways were found for enantiomerization via the pyramidal inversion mechanism for each of the tautomers of (S)‐omeprazole. Each transition state, in which the sulfur, the oxygen and the two carbon atoms connected directly to the sulfur are in one plane, connects two diastereomeric minima. The enantiomerization is completed by free rotation around the sulfur–methylene bond, and around the methylene–pyridine ring bond. The effective Gibbs' free energy barrier for racemization ΔG of the two tautomers of (S)‐omeprazole are 39.8 kcal/mol (5‐methoxy tautomer) and 40.0 kcal/mol (6‐methoxy tautomer), indicating that the enantiomers of omeprazole are stable at room temperature (in the gas phase). The 5‐methoxy tautomer of (S)‐omeprazole was found to be slightly more stable than the 6‐methoxy tautomer, in the gas phase. The energy barrier (ΔG^?) for the(S,M) (S,P) diastereomerization of (S)‐omeprazole due to the rotation around the pyridine chiral axis was very low, 5.8 kcal/mole at B3LYP/6‐311G(d,p). Chirality 2010. © 2010 Wiley‐Liss, Inc.     

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.     

Rational design of a lipase to accommodate catalysis of Baeyer–Villiger oxidation with hydrogen peroxide

The mechanism and potential energy surface for the Baeyer-Villiger oxidation of acetone with hydrogen peroxide catalyzed by a Ser105-Ala mutant of Candida antarctica Lipase B has been determined using ab initio and density functional theories. Initial substrate binding has been studied using an automated docking procedure and molecular dynamics simulations. Substrates were found to bind to the active site of the mutant. The activation energy for the first step of the reaction, the nucleophilic attack of hydrogen peroxide on the carbonyl carbon of hydrogen peroxide, was calculated to be 4.4 kcal x mol(-1) at the B3LYP/6-31+G* level. The second step, involving the migration of the alkyl group, was found to be the rate-determining step with a computed activation energy of 19.9 kcal x mol(-1) relative the reactant complex. Both steps were found to be lowered considerably in the reaction catalyzed by the mutated lipase, compared to the uncatalyzed reaction. The first step was lowered by 36.0 kcal x mol(-1) and the second step by 19.5 kcal x mol(-1). The second step of the reaction, the rearrangement step, has a high barrier of 27.7 kcal x mol(-1) relative to the Criegee intermediate. This could lead to an accumulation of the intermediate. It is not clear whether this result is an artifact of the computational procedure, or an indication that further mutations of the active site are required. Figure Second TS (18TS) in the Baeyer-Villiger oxidation in a mutant of CALB. Distances in A