Quantum Chemical Studies on Stability and Chemical Activities in Calcium Ion Bound Calmodulin Loops

Quantum chemical (QC) calculations for macromolecules require truncation of the molecule, highlighting the portion of interest due to heavy computation cost. As a result, an estimation of the effects of truncation is important to interpret the energy spectrum of such calculations. We perform density functional theory based QC calculations on calcium ion bound EF-hand loops of Calmodulin isolated from the crystal structure in an implicit solvent. We find that the terminal contributions of neutral capping are negligible across the entire ground-state energy spectrum. The coordination energy range and the nature of hybridization of the coordination state molecular orbitals remain qualitatively similar across these loops. While the HOMO and LUMO of loops in the N-terminal domain are dominated by the acidic aspartates, and the polar/hydrophobic residues, respectively, these levels of the C-terminal domain loops show strong localized electron density on the phenyl rings of the tyrosines. The Fukui index calculation identifies the hydroxyl oxygen in the phenyl ring of Y99 as a potent nucleophile. Our analysis indicates a general way of interpreting the electronic energy spectra to understand stability and functions of large biomolecules where the truncation of the molecule and, hence, the terminal capping effects are inevitable.