Theoretical approach to drug design. 2. Relative thermodynamics of inhibitor binding by chicken dihydrofolate reductase to ethyl derivatives of trimethoprim substituted at 3'-, 4'-, and 5'-positions

Modeling Protein-Ligand Binding by Mining Minima

We present the first application of the mining minima algorithm to protein-small molecule binding. This end-point approach use an empirical force field and implicit solvent models, treats the protein binding-site as fully flexible and estimates free energies as sums over local energy wells. The calculations are found to yield encouraging agreement with experiment for three sets of HIV-1protease inhibitors and a set of phosphodiesterase 10a inhibitors. The contributions of various aspects of the model to its accuracy are examined, and the Poisson-Boltzmann correction is found to be the most critical. Interestingly, the computed changes in configurational entropy upon binding fall roughly along the same entropy-energy correlation previously observed for smaller host-guest systems. Strengths and weaknesses of the method are discussed, as are the prospects for enhancing accuracy and speed.

CHARMM: the biomolecular simulation program

CHARMM (Chemistry at HARvard Molecular Mechanics) is a highly versatile and widely used molecular simulation program. It has been developed over the last three decades with a primary focus on molecules of biological interest, including proteins, peptides, lipids, nucleic acids, carbohydrates, and small molecule ligands, as they occur in solution, crystals, and membrane environments. For the study of such systems, the program provides a large suite of computational tools that include numerous conformational and path sampling methods, free energy estimators, molecular minimization, dynamics, and analysis techniques, and model-building capabilities. The CHARMM program is applicable to problems involving a much broader class of many-particle systems. Calculations with CHARMM can be performed using a number of different energy functions and models, from mixed quantum mechanical-molecular mechanical force fields, to all-atom classical potential energy functions with explicit solvent and various boundary conditions, to implicit solvent and membrane models. The program has been ported to numerous platforms in both serial and parallel architectures. This article provides an overview of the program as it exists today with an emphasis on developments since the publication of the original CHARMM article in 1983.

Domain motions in dihydrofolate reductase: a molecular dynamics study

Molecular dynamics simulations have been carried out on the enzyme dihydrofolate reductase from Lactobacillus casei complexed with methotrexate, NADPH and 264 crystallographic water molecules. Analysis of correlations in atomic fluctuations reveal the presence of highly correlated motion (correlation coefficient > 0.6) in the region between residues 30 to 35 and 85 to 90 leading to the identification of two domains, an "adenosine-binding domain" and a "large domain", which rotate by 3 to 4 degrees with respect to each other. The strongest correlation (> 0.6) within the large domain involves a coupling between the motions of the "teen-loop", and the spatially contiguous loops linking beta 6-beta 7 and beta 7-beta 8. Moreover, there is a significant correlation (approximately 0.5) between the adenosine fragment of NADPH and the pteridine and p-aminobenzoyl fragments of methotrexate, which are separated by approximately 17 A, and is lost on removal of "rigid-body" motion from the original trajectory. This provides support for the idea that the relative motion of the two domains is a means by which the occupation of the binding site for the adenosine end of the coenzyme can affect methotrexate binding and vice versa. Quasiharmonic vibrational analysis of the trajectory reveals that the overall dynamics of the system are governed by domain motions whose contributions are dominant at low frequencies. In addition, different low-frequency modes are responsible for separately coupling the adenosine-binding site and parts of methotrexate.

Model assembly study of the ligand binding by p-hydroxybenzoate hydroxylase: correlation between the calculated binding energies and the experimental dissociation constants

The energies of binding of seven ligands by p-hydroxybenzoate hydroxylase (PHBH) were calculated theoretically. Direct enzyme-ligand interaction energies were calculated using the ab initio quantum mechanical model assembly of the active site at the 3-21G level. Solvation energies of the ligands needed in the evaluation of the binding energies were calculated with the semiempirical AM1-SM2 method and the long-range electrostatic interaction energies between the ligands and the protein matrix classically using the static charge distributions of the ligands and the protein. Energies for proton-transfer between the ligands' OH or SH substituent at position 4 and the active-site tyrosine within the ab initio model assemblies were calculated and compared to the corresponding pKas in aqueous solution. Excluding 3,4-dihydroxybenzoate, the natural product of PHBH, a linear relationship between the calculated binding energies and the experimental binding free energies was found with a correlation coefficient of 0.90. Contributions of the direct enzyme-ligand interaction energies, solvation energies and the long-range electrostatic interaction energies to the calculated binding energies were analyzed. The proton-transfer energies of the ligands with substituents ortho to the ionized OH were found to be perturbed less in the model calculations than the energies of their meta isomers as deduced from the corresponding pKas.

Molecular recognition in proteins. Simulation analysis of substrate binding by a tyrosyl-tRNA synthetase mutant

Alchemical molecular dynamics simulations are performed to determine the difference in the free energy of binding of the tyrosine substrate between the wild type of tyrosyl-tRNA synthetase (TyrRS) from Bacillus stearothermophilus and the mutant Tyr169-->Phe. The results are of general interest because the Tyr169 hydroxyl group interacts with the ammonium group of the substrate in a manner corresponding to that found in other amino acid binding proteins (e.g. the Asp receptor of the chemotactic bacterium Salmonella typhimurium and class I major histocompatibility complex molecules). The calculated free-energy change due to the Tyr169-->Phe mutation is 3.4 kcal/mol (the statistical error is +/- 0.5 kcal/mol) in satisfactory agreement with the experimental value of 3(+/- 0.5) kcal/mol. By use of thermodynamic integration, the contribution of the different terms to the free energy change are estimated. The path dependence of such a decomposition is discussed and it is suggested that the alchemical choice is of primary interest for understanding the interactions involved. There are large protein contributions to the alchemical free energy difference of the bound and free enzyme that cancel in the overall result. Due to this cancellation, the essential interactions contributing to the free-energy change are those between the OH group of Tyr169 and water in the free enzyme and those between the OH group of Tyr169 and the ammonium group of the substrate in the bound system. The results thus support simple models based on a balance of hydrogen bonding interactions.

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.