Recent years have seen tremendous developments in methods for computational modeling of (bio-) molecular systems. Ever larger reactive systems are being studied with high accuracy approaches, and high-level QM/MM calculations are being routinely performed. However, applying high-accuracy methods to large biological systems is computationally expensive and becomes problematic when conformational sampling is needed. To address this challenge, classical force field based approaches such as free energy perturbation (FEP) and empirical valence bond calculations (EVB) have been employed in this work. Specifically: Force-field independent metal parameters have been developed for a range of alkaline earth and transition metal ions, which successfully reproduce experimental solvation free energies, metal-oxygen distances, and coordination numbers. These are valuable for the computational study of biological systems. Experimental studies have shown that the epoxide hydrolase from Solanum tuberosum (StEH1) is not only an enantioselective enzyme, but for smaller substrates, displays enantioconvergent behavior. For StEH1, two detailed studies, involving combined experimental and computational efforts have been performed: We first used trans-stilbene oxide to establish the basic reaction mechanism of this enzyme. Importantly, a highly conserved and earlier ignored histidine was identified to be important for catalysis. Following from this, EVB and experiment have been used to investigate the enantioconvergence of the StEH1-catalyzed hydrolysis of styrene oxide. This combined approach involved wildtype StEH1 and an engineered enzyme variant, and established a molecular understanding of enantioconvergent behavior of StEH1. A novel framework was developed for the Computer-Aided Directed Evolution of Enzymes (CADEE), in order to be able to quickly prepare, simulate, and analyze hundreds of enzyme variants. CADEE’s easy applicability is demonstrated in the form of an educational example. In conclusion, classical approaches are a computationally economical means to achieve extensive conformational sampling. Using the EVB approach has enabled me to obtain a molecular understanding of complex enzymatic systems. I have also increased the reach of the EVB approach, through the implementation of CADEE, which enables efficient and highly parallel in silico testing of hundreds-to-thousands of individual enzyme variants.