The selective N-acylation of 1,2-amino alcohols has been proposed to occur through the proton shuttle mechanism. However, the O-acetylation of propranolol catalyzed by Candida antarctica lipase B is an exception. We investigated the relation between the chemoselectivity of this reaction and the acyl group length. For this purpose, we compared the acyl groups: ethanoyl, butanoyl, octanoyl, and hexadecanoyl. We studied the Michaelis complexes between serine-acylated Candida antarctica lipase B and propranolol, employing a computational approach that involved sampling Michaelis complex conformations through ensemble docking plus consensus scoring and molecular dynamics simulations. The conformations were then classified as near attack conformations for acylation of the amino or hydroxy group. The relative populations of these two classes of conformations were found to be consistent with the experimentally observed chemoselective O-acetylation. We predict that increasing the length of the hydrocarbon chain of the acyl group will cause O-acylation to be unfavorable with respect to N-acylation. The nucleophilic attack of propranolol to the acylated lipase was found to be more favorable through the classical mechanism when compared with the proton shuttle mechanism.
Keywords: Amidation; Amino alcohols; Enzyme catalysis; Esterification; Molecular modeling; Near attack conformations.