Enzymes are extraordinarily effective catalysts evolved to perform well-defined and highly specific chemical transformations. Studying the nature of rate enhancements and the mechanistic strategies in enzymes is very important, both from a basic scientific point of view, as well as in order to improve rational design of biomimetics. Kinetic isotope effect (KIE) is a very important tool in the study of chemical reactions and has been used extensively in the field of enzymology. Theoretically, the prediction of KIEs in condensed phase environments such as enzymes is challenging due to the need to include nuclear quantum effects (NQEs). Herein we describe recent progress in our group in the development of multiscale simulation methods for the calculation of NQEs and accurate computation of KIEs. We also describe their application to several enzyme systems. In particular we describe the use of combined quantum mechanics/molecular mechanics (QM/MM) methods in classical and quantum simulations. The development of various novel path-integral methods is reviewed. These methods are tailor suited to enzyme systems, where only a few degrees of freedom involved in the chemistry need to be quantized. The application of the hybrid QM/MM quantum-classical simulation approach to three case studies is presented. The first case involves the proton transfer in alanine racemase. The second case presented involves orotidine 5'-monophosphate decarboxylase where multidimensional free energy simulations together with kinetic isotope effects are combined in the study of the reaction mechanism. Finally, we discuss the proton transfer in nitroalkane oxidase, where the enzyme employs tunneling as a catalytic fine-tuning tool. (C) 2015 Elsevier Inc. All rights reserved.