Comparative studies for evaluation of CO2 fixation in the cavity of the Rubisco enzyme using QM, QM/MM and linear-scaling DFT methods

Abstract

We evaluate the minimum energy configuration (MM) and binding free energy (QM/MM and QM) of CO2 to Rubisco, of fundamental importance to the carboxylation step of the reaction. Two structural motifs have been used to achieve this goal, one of which starts from the initial X-ray Protein Data Bank structure of Rubisco's active centre (671 atoms), and the other is a simplified, smaller model (77 atoms) which has been used most successfully, thus far, for study. The small model is subjected to quantum chemical density functional theory (DFT) studies, both in vacuo and using implicit solvation. The effects of the protein environment are also included by means of a hybrid quantum mechanical/molecular mechanical (QM/MM) approach, using PM6/AMBER and B3LYP/AMBER schemes. Finally, linear-scaling DFT methods have also been applied to evaluate energetic features of the large motif, and the result obtained for the binding free energy of the CO2 underlines the importance of the accurate modelling of the surrounding protein milieu using a full DFT description.