Computational study of the enantioselectivity of the O-acetylation of (R,S)-propranolol catalyzed by Candida antarctica lipase B

Application of molecular dynamics simulation in the field of food enzymes: improving the thermal-stability and catalytic ability

Enzymes can produce high-quality food with low pollution, high function, high acceptability, and medical aid. However, most enzymes, in their native form, do not meet the industrial requirements. Sequence-based and structure-based methods are the two main strategies used for enzyme modification. Molecular Dynamics (MD) simulation is a sufficiently comprehensive technology, from a molecular perspective, which has been widely used for structure information analysis and enzyme modification. In this review, we summarize the progress and development of MD simulation, particularly for software, force fields, and a standard procedure. Subsequently, we review the application of MD simulation in various food enzymes for thermostability and catalytic improvement was reviewed in depth. Finally, the limitations and prospects of MD simulation in food enzyme modification research are discussed. This review highlights the significance of MD simulation and its prospects in food enzyme modification.

Hydrogen Bonds and n → π* Interactions in the Acetylation of Propranolol Catalyzed by Candida antarctica Lipase B: A QTAIM Study

Enzyme-substrate interactions play a crucial role in enzymatic catalysis. Quantum theory of atoms in molecules (QTAIM) calculations are extremely useful in computational studies of these interactions because they provide very detailed information about the strengths and types of molecular interactions. QTAIM also provides information about the intramolecular changes that occur in the catalytic reaction. Here, we analyze the enzyme-substrate interactions and the topological properties of the electron density in the enantioselective step of the acylation of (R,S)-propranolol, an aminoalcohol with therapeutic applications, catalyzed by Candida antarctica lipase B. Eight reaction paths (four for each enantiomer) are investigated and the energies, atomic charges, hydrogen bonds, and n → π* interactions of propranolol, the catalytic triad (composed of D187, H224, and S105), and the oxyanion hole are analyzed. It is found that D187 acts as an electron density reservoir for H224, and H224 acts as an electron density reservoir for the active site of the protein. It releases electron density when the tetrahedral intermediate is formed from the Michaelis complex and receives it when the enzyme-product complex is formed. Hydrogen bonds can be grouped into noncovalent and covalent hydrogen bonds. The latter are stronger and more important for the reaction than the former. We also found weak n → π* interactions, which are characterized by QTAIM and the natural bond orbital (NBO) analysis.

Finite-temperature effect in the O-acylation of (R,S)-propranolol catalyzed by Candida antarctica lipase B

CalB is a triacylglycerol hydrolase (E.C.3.1.1.3) used in the O-acylation of the beta-adrenergic blocking agent (R,S)-propranolol. The catalytic mechanism involves two steps: enzyme acylation and enzyme deacylation. The enantioselectivity of the O-acylation of (R,S)-propranolol originates from the second step, where the acyl-enzyme transfers the acyl group to the racemic substrate. This step occurs via an initial Michaelis complex (MCC) and a tetrahedral intermediate (TI-2). To gain more insight into the molecular basis of this reaction, we performed an exhaustive conformational sampling along the reaction coordinate of the enantioselective step of the reaction (MCC→TI-2→EPC) applying a QM/MM MD protocol (SCC-DFTB/CHARMM) in combination with umbrella sampling and the weighted histogram analysis method. To identify finite temperature effects we compare the PMF and the potential energy pathway. It is found that the effect of the finite temperature in this reaction is a destabilization of the tetrahedral intermediate and an increase of the barrier height of its formation. This increase is higher for the S-enantiomer.

Chemoenzymatic synthesis of enantiomerically enriched diprophylline and xanthinol nicotinate

A concise chemoenzymatic route toward enantiomerically enriched active pharmaceutical ingredients (API) - diprophylline and xanthinol nicotinate - is reported for the first time. The decisive step is an enantioselective lipase-mediated methanolysis of racemic chlorohydrin-synthon acetate, namely 1-chloro-3-(1,3-dimethyl-2,6-dioxo-1,2,3,6-tetrahydro-7H-purin-7-yl)propan-2-yl acetate, performed under kinetically-controlled conditions on a preparative 500 mg-scale. The best results in terms of reaction enantioselectivity (E = 14) were obtained for the enantiomers resolution performed with lipase type B from Candida antarctica immobilized on acrylic resin (CAL-B, Novozym 435) suspended in homophasic acetonitrile-methanol mixture. The elaborated biocatalytic system furnished the key chlorohydrin intermediate (in 71% ee and 38% yield), which was then smoothly converted into enantioenriched active agents: (R)-(-)-diprophylline (57% ee) and (S)-(+)-xanthinol nicotinate (65% ee). To support the assignment of absolute configurations of EKR-products as well as to confirm the stereochemical outcome of the remaining reaction steps, docking studies toward the prediction of enantiomers binding selectivity in CAL-B active site as well as the respective chemical correlations with enantiomerically enriched analytical standards obtained from commercially available (R)-(-)-epichlorohydrin, were applied. In addition, single-crystal X-ray diffraction (XRD) analyses were performed for the synthesized optically active APIs furnishing by this manner a first crystal structures of nicotinic acid salt of xanthinol.

Molecular dynamics study of taxadiene synthase catalysis

Molecular dynamics (MD) simulations have been performed to study the dynamic behavior of noncovalent enzyme carbocation complexes involved in the cyclization of geranylgeranyl diphosphate to taxadiene catalyzed by taxadiene synthase (TXS). Taxadiene and the observed four side products originate from the deprotonation of carbocation intermediates. The MD simulations of the TXS carbocation complexes provide insights into potential deprotonation mechanisms of such carbocations. The MD results do not support a previous hypothesis that carbocation tumbling is a key factor in the deprotonation of the carbocations by pyrophosphate. Instead water bridges are identified which may allow the formation of side products via multiple proton transfer reactions. A novel reaction path for taxadiene formation is proposed on the basis of the simulations. © 2018 Wiley Periodicals, Inc.

Función de los confórmeros de ataque cercano en la acilación enantioselectiva del (R,S)-propranolol catalizada por lipasa B de Candida antarctica

La lipasa B de Candida antarctica (CalB) se ha utilizado en la acilación quimio- y enantioselectiva del racemato (R,S)-propranolol. CalB tiene enantioselectividad moderada (E=63) por el R-propranolol. La enantioselectividad, se origina en la reacción de transferencia del grupo acilo desde la serina catalítica, acilada, al propranolol. La fase inicial de esta reacción involucra la formación de complejos de Michaelis y posteriormente conformaciones de ataque cercano. El análisis de las conformaciones de ataque cercano ha permitido en varios casos explicar el origen de la catálisis o reproducir el efecto catalítico. En este trabajo se profundiza en la comprensión la función de las conformaciones de ataque cercano en la enantioselectividad de la acilación del (R,S)-propranolol catalizada por CalB. Para lo anterior se realizó un estudio detallado de los complejos de Michaelis y de las conformaciones de ataque cercano del paso enantioselectivo de la reacción de acilación del (R,S)-propranolol utilizando un protocolo de dinámica molecular QM/MM (SCCDFTB/CHARMM) utilizando 6 distribuciones de velocidades iniciales y simulaciones de 2,5 ns. Se estudiaron 7 complejos CalB-propranolol. Los enlaces de hidrógeno del sitio activo de CalB acilada relevantes para la actividad catalítica fueron estables en todas las simulaciones. Las poblaciones de los complejos de Michaelis y de las conformaciones de ataque cercano son dependientes de la distribución de las velocidades iniciales de la dinámica molecular. La enantioselectividad moderada de CalB acilada, encontrada experimentalmente, puede ser parcialmente atribuida a la alta población de conformaciones de ataque cercano observada para el S-propranolol.