Journal of Physical Chemistry A, Vol.102, No.30, 6152-6160, 1998
Ab initio investigation of the methylimidazole-indole complexes as models of the histidine-tryptophan pair
Several histidine-tryptophan complexes, derived from the crystal structures available in the Brookhaven Protein Data Bank, have been examined with ab initio theoretical methods (using as model systems 5-methylimidazole and indole, respectively), in order to identify the most favorable arrangements of the two side chains, elucidating also the strength and the nature of the intermolecular interaction established between them. The equilibrium geometries of the isolated partners were optimized at the HF/6-31G* level and the interaction energy of the adducts was computed, employing the 6-31G* basis set with the d exponents reduced to 0.25, thus named 6-31G*(0.25), at the HF and MP2 (frozen-core approximation) levels. For a few typical orientations, the dependence of the interaction energy upon the intermolecular distance, as measured from the ring centroids, was then examined while keeping fixed reciprocal orientations and internal geometries of the partners. There is a fair linear correlation between the equilibrium distances (R-eq) at the MP2 level and the experimental (R-exp) ones and between the MP2 interaction energies at R-eq and those computed at R-exp. For three arrangements with a shallow or even repulsive HF interaction energy, the counterpoise correction to the basis set superposition error (BSSE) was introduced both at the HF and MP2 levels, using Pople's 6-31G* standard, 6-31G*(0.25), and Dunning's DZP basis sets, to test the reliability of the results obtained along the whole approaching path. This is made necessary by the noticeable displacement in the equilibrium distances usually found at the various levels. The DZP HF interaction energies turn out to be less affected by BSSEs than the 6-31G* and the 6-31G*(0.25) ones and are located in an intermediate position between them. As a general rule for these complexes, the counterpoise correction is larger at the correlated level; therefore the addition of the correlation effect to the counterpoise-corrected SCF energy produces a curve fairly close to the MP2 one that seems to represent a lower bound to the true interaction energy. Kitaura and Morokuma's decomposition analysis of the interaction energies was also carried out on these typical complexes.
Keywords:SET SUPERPOSITION ERRORS;HYDROGEN-BONDED DIMERS;DNA-BASE PAIRS;MOLECULAR-INTERACTIONS;COUNTERPOISE CORRECTIONS;INTERACTION ENERGY;ELECTRON CORRELATION;CYTOSINE DIMER;WATER;DECOMPOSITION