A Full QM Computational Study of the Catalytic Mechanism of α-1,4-Glucan Lyases.

We investigated the catalytic mechanism of α-1,4-glucan lyases using a full QM DFT approach based on the M06-2X functional. The reaction profile of the whole catalytic process can be divided into three phases: glycosylation, deglycosylation-elimination and tautomerization. Glycosylation is a highly asynchronous SN 1-like process with an energy barrier of 10.2 kcal mol-1 . A proton moves from the Asp665 residue to the glycosidic oxygen. Asp553 acts as a nucleophile and attacks the anomeric carbon causing the cleavage of the glycosidic bond. Deglycosilation-elimination is the rate-determining step of the entire process with an overall barrier of 18.3 kcal mol-1 . The final step (restoring the catalyst and tautomerization) occurs rather easily, since the Asp553 carboxylate group "assists" the proton transfer in the tautomerization process. Our computations clearly indicate that tautomerization must occur inside the enzyme before leaving the active site rather than in the aqueous solution. Outside of the protein environment the enol-AF→keto-AF process "assisted" by a water molecule has a barrier of 35.8 kcal mol-1 .

The Reaction Pathway of Cellulose Pyrolysis to a Multifunctional Chiral Building Block: The Role of Water Unveiled by a DFT Computational Investigation.

LAC (hydroxylactone (1R,5S)-1-hydroxy-3,6-dioxabicyclo[3.2.1]octan-2-one) is one of the most interesting products of the pyrolysis of cellulose and represents a useful chiral building block in organic synthesis. A computational investigation at the DFT level on the mechanism of formation of LAC shows that this species can be obtained following two reaction paths, path A and path B, starting from a well-known pyrolysis product (ascopyrone P). A series of internal rearrangements involving in all cases a proton transfer leads directly to LAC (path B). An alternative path (path A) can be also followed. From this path, via a "gate" connecting the two reaction channels, it is possible to reach path B and form LAC. In both cases, the rate-determining step of the process is the initial keto-enol isomerization. We found that water, which is present in the reaction mixture, "catalyzes" the reaction by assisting the proton transfers present in all the steps of the process. In particular, water lowers the barrier of the rate-determining step that becomes 40.9 kcal mol-1 (79.4 kcal mol-1 in the absence of water). The corresponding computed rate constant is 4.3×10 s-1 at 500 °C, a value which is consistent with the presence of LAC in the absence of metal catalysts. The results of this study on the non-catalyzed process underpin the important role played by water in the formation of pyrolysis products of cellulose where proton transfer is a key mechanistic step.

Computational Evidence for the Catalytic Mechanism of Tyrosylprotein Sulfotransferases: A Density Functional Theory Investigation.

In this paper we have examined the mechanism of tyrosine O-sulfonation catalyzed by human TPST-2. Our computations, in agreement with Teramoto's hypothesis, indicate a concerted SN2-like reaction (with an activation barrier of 18.2 kcal mol(-1)) where the tyrosine oxygen is deprotonated by Glu(99) (base catalyst) and simultaneously attacks as a nucleophile the sulfuryl group. For the first time, using a quantum mechanics protocol of alanine scanning, we identified unequivocally the role of the amino acids involved in the catalysis. Arg(78) acts as a shuttle that "assists" the sulfuryl group moving from the 3'-phosphoadenosine-5'-phosphosulfate molecule to threonine and stabilizes the transition state (TS) by electrostatic interactions. The residue Lys(158) keeps close the residues participating in the overall H-bond network, while Ser(285), Thr(81), and Thr(82) stabilize the TS via strong hydrogen interactions and contribute to lower the activation barrier.

Gold(I)-Catalyzed Dearomative [2+2]-Cycloaddition of Indoles with Activated Allenes: A Combined Experimental-Computational Study.

The gold-catalyzed synthesis of methylidene 2,3-cyclobutane-indoles is documented through a combined experimental/computational investigation. Besides optimizing the racemic synthesis of the tricyclic indole compounds, the enantioselective variant is presented to its full extent. In particular, the scope of the reaction encompasses both aryloxyallenes and allenamides as electrophilic partners providing high yields and excellent stereochemical controls in the desired cycloadducts. The computational (DFT) investigation has fully elucidated the reaction mechanism providing clear evidence for a two-step reaction. Two parallel reaction pathways explain the regioisomeric products obtained under kinetic and thermodynamic conditions. In both cases, the dearomative CC bond-forming event turned out to be the rate-determining step.

Gold(I)-Assisted α-Allylation of Enals and Enones with Alcohols.

The intermolecular α-allylation of enals and enones occurs by the condensation of variously substituted allenamides with allylic alcohols. Cooperative catalysis by [Au(ItBu)NTf2] and AgNTf2 enables the synthesis of a range of densely functionalized α-allylated enals, enones, and acyl silanes in good yield under mild reaction conditions. DFT calculations support the role of an α-gold(I) enal/enone as the active nucleophilic species.

CNT-Confinement Effects on the Menshutkin SN2 Reaction: The Role of Nonbonded Interactions.

We investigated the effects of CNT confinement ((6,6) tube) on the model Menshutkin reaction H3N + H3CCl = H3NCH3((+)) + Cl((-)), which is representative of chemical processes involving developing of charge separation along the reaction pathway. We used either a full QM approach or a hybrid QM/MM approach. We found that the CNT significantly lowers the activation barrier with respect to the hypothetical gas-phase reaction: The activation barrier Ea varies from 34.6 to 25.7 kcal mol(-1) (a value similar to that found in a nonpolar solvent) and the endothermicity ΔE from 31.2 to 13.5 kcal mol(-1). A complex interplay between C-H···π, N-H···π, and Cl···π nonbonded interactions of the endohedral system with the CNT wall explains the lower barrier and lower endothermicity. The hybrid QM/MM approach (MM = UFF force field) does not reproduce satisfactorily the QM energy ΔE (18.1 vs 13.5 kcal mol(-1)), while optimum agreement is found in the barrier Ea (25.8 vs 25.7 kcal mol(-1)). These results suggest that the simple Qeq formalism (included in the MM potential) does not describe properly the effect of CNT polarization in the presence of the net charge separation featuring the final product. A more accurate estimate of the tube polarization was obtained with single-point QM/MM computations including PCM corrections (using the benzene dielectric constant) on the QM/MM optimized structures. After PCM corrections, Ea changes slightly (from 25.8 to 24.5 kcal mol(-1)), but a more significant variation is observed for ΔE that becomes 13.1 kcal mol(-1), in rather good agreement with the full QM. This level of theory (QM/MM with PCM correction, MM = UFF) represents a more general approach suitable for describing CNT-confined chemical processes involving significant charge separation. QM/MM computations were extended to CNTs of different radii: (4,4), (5,5), (7,7), (8,8), (9,9), (10,10), (12,12), (14,14) CNTs and, as a limit case, a graphene sheet. The lack of space available in the small tube (4,4) causes a strong structural distortion and a consequent increase in Ea and ΔE (40.8 and 44.0 kcal mol(-1), respectively). These quantities suddenly decrease with the augmented volume inside the (5,5) tube. For larger tubes, different structural arrangements of the endohedral system are possible, and Ea and ΔE remain almost constant until the limiting case of graphene.