An iterative optimization protocol implemented in LICHEM on the QM region and all atoms within 30 Å of the O3' on the five chosen representatives of each system, the representative with the lowest energy was selected as a starting structure for product formation. The product structures were generated from the reactant structures using GaussView and optimized using the same iterative optimization approach.