Atomic charges derived from semiempirical electrostatic potentials; an interaction energy method

Extension of the MST continuum solvation model to the RM1 semiempirical Hamiltonian

The need to simulate large-sized molecules or to deal with large series of compounds is a challenging topic in computational chemistry, which has stimulated the development of accurate semiempirical methods, such as the recently reported Recife Model 1 (RM1; J Comput Chem 2006, 27, 1101). Even though RM1 may prove to be of value simply due to the enhanced quantitative accuracy in gas phase, it is unclear how the new parameters optimized for RM1 affect the suitability of this semiempirical Hamiltonian to study chemical processes in condensed phases. To address this question, we report the parametrization of the MST/RM1 continuum model for neutral solutes in water, octanol, chloroform and carbon tetrachloride, and for ions in water. The results are used to discuss the transferability of the solvation parameters implemented in previous MST/AM1 and MST/PM3 models.

Multipole electrostatic potential derived atomic charges in NDDO-methods with spd-basis sets

The recently introduced multipole approach for computing the molecular electrostatic potential (MEP) within the semiempirical neglect of diatomic differential overlap (NDDO) framework [Horn AHC, Lin Jr-H., Clark T (2005) Theor Chem Acc 114:159-168] has been used to obtain atomic charges of nearly ab initio quality by scaling the semiempirical MEP. The parameterization set comprised a total of 797 compounds and included not only the newly parameterized AM1* elements Al, Si, P, S, Cl, Ti, Zr, and Mo but also the standard AM1 elements H, C, N, O and F. For comparison, the ZDO-approximated MEP was also calculated analytically in the spd-basis. For the AM1*-optimized structures, single-point calculations at the B3LYP, HF and MP2 levels with the 6-31G(d) and LanL2DZP basis sets were performed to obtain the MEP. The regression analysis of all 12 combinations of semiempirical and ab initio MEP data yielded correlation coefficients of at least 0.99 in all cases. Scaling the analytical and multipole-derived semiempirical MEP by the regression coefficients yielded mean unsigned errors below 2.6 and 1.9 kcal mol(-1), respectively. Subsequently, for 22 drug molecules from the World Drug Index, atomic charges were computed according to the RESP procedure using XX/6-31G(d) (XX=B3LYP, HF, MP2) and scaled AM1* multipole MEP; the correlation coefficients obtained are 0.83, 0.85 and 0.83, respectively.

Transferability of fragmental contributions to the octanol/water partition coefficient: an NDDO-based MST study

This study examines the transferability of fragmental contributions to the octanol/water partition coefficient. As a previous step, we report the parameterization of the AM1 and PM3 versions of the MST model for n-octanol. The final AM1 and PM3 MST models reproduce the experimental free energy of solvation and the octanol/water partition coefficient (log P(ow)) with a root-mean-square deviation of around 0.7 kcal/mol and 0.5 (in units of log P), respectively. Based on this parameterization, an NNDO-based procedure is presented to dissect the free energy of transfer between octanol and water in contributions directly associated with specific atoms or functional groups. The application of this procedure to a set of representative molecular systems illustrates the dependence of the log P(ow) fragmental contribution due to electronic, hydrogen bonding, and steric effects, which cannot be easily accounted for in simple additive-based empirical schemes. The results point out the potential use of theoretical methods to refine the fragmental contributions in empirical methods.

Computer-aided drug design: a free energy perturbation study on the binding of methyl-substituted pterins and N5-deazapterins to dihydrofolate reductase

Molecular dynamics simulation and free energy perturbation techniques have been used to study the relative binding free energies of 8-methylpterins and 8-methyl-N5-deazapterins to dihydrofolate reductase (DHFR). Methyl-substitution at the 5, 6 and 7 positions in the N-heterocyclic ring gives rise to a variety of ring substituent patterns and biological activity: several of these methyl derivatives of the 8-methyl parent compounds (8-methylpterin and 8-methyl-N5-deazapterin) have been identified as substrates or inhibitors of vertebrate DHFR in previous work. The calculated free energy differences reveal that the methyl-substituted compounds are thermodynamically more stable than the primary compounds (8-methylpterin and 8-methyl-N5-deazapterin) when bound to the enzyme, due largely to hydrophobic hydration phenomena. Methyl substitution at the 5 and/or 7 positions in the 6-methyl-substituted compounds has only a small effect on the stability of ligand binding. Furthermore, repulsive interactions between the 6-methyl substituent and DHFR are minimal, suggesting that the 6-methyl position is optimal for binding. The results also show that similarly substituted 8-methylpterins and 8-methyl-N5-deazapterins have very similar affinities for binding to DHFR. The computer simulation predictions are in broad agreement with experimental data obtained from kinetic studies, i.e. 6,8-dimethylpterin is a more efficient substrate than 8-methylpterin and 6,8-dimethyl-N5-deazapterin is a better inhibitor than 8-methyl-N5-deazapterin.

Novel mechanism-based substrates of dihydrofolate reductase and the thermodynamics of ligand binding: a comparison of theory and experiment for 8-methylpterin and 6,8-dimethylpterin

Molecular dynamics simulation and free energy perturbation techniques have been used to study the relative binding free energies of the designed mechanism-based pterins, 8-methylpterin and 6,8-dimethylpterin, to dihydrofolate reductase (DHFR), with cofactor nicotinamide adenine dinucleotide phosphate (NADPH). The calculated free energy differences suggest that DHFR.NADPH.6,8-dimethylpterin is thermodynamically more stable than DHFR.NADPH.8-methylpterin by 2.4 kcal/mol when the substrates are protonated and by 1.3 kcal/mol when neutral. The greater binding strength of 6,8-dimethylpterin may be attributed largely to hydration effects. In terms of an appropriate model for the pH-dependent kinetic mechanism, these differences can be interpreted consistently with experimental data obtained from previous kinetic studies, i.e., 6,8-dimethylpterin is a more efficient substrate of vertebrate DHFRs than 8-methylpterin. The kinetic data suggest a value of 6.6 +/- 0.2 for the pKa of the active site Glu-30 in DHFR.NADPH. We have also used experimental data to estimate absolute values for thermodynamic dissociation constants of the active (i.e., protonated) forms of the substrates: these are of the same order as for the binding of folate (0.1-10 microM). The relative binding free energy calculated from the empirically derived dissociation constants for the protonated forms of 8-methylpterin and 6,8-dimethylpterin is 1.4 kcal/mol, a value which compares reasonably well with the theoretical value of 2.4 kcal/mol.