Quantum chemical geometry optimizations in proteins using crystallographic raw data

A method is developed for the combination of quantum chemical geometry optimizations and crystallographic structure refinement. The method is implemented by integrating the quantum chemical software Turbomole with the crystallographic software Crystallography and NMR System (CNS), using three small procedures transferring information between the two programs. The program (COMQUM-X)is used to study the binding of the inhibitor N-methylmesoporphyrin to ferrochelatase, and we show that the method behaves properly and leads to an improvement of the structure of the inhibitor. It allows us to directly quantify in energy terms how much the protein distort the structure of the bound inhibitor compared to the optimum vacuum structure (4-6 kJ/mol). The approach improves the standard combined quantum chemical and molecular mechanics (QC/MM) approach by guaranteeing that the final structure is in accordance with experimental data (the reflections) and avoiding the risk of propagating errors in the crystal coordinates. The program can also be seen as an improvement of standard crystallographic refinement, providing an accurate empirical potential function for any group of interest. The results can be directly interpreted in standard crystallographic terms (e.g., R factors or electron density maps). The method can be used to interpret crystal structures (e.g., the protonation status of metal-bound water molecules) and even to locally improve them.