On the binding of aminoalkyl adenylates to isoleucyl-tRNA synthetase from Escherichia coli MRE 600.

The binding of nine aminoalkyl adenylates to isoleucyl-tRNA synthetase from Escherichia coli MRE 600 was measured and compared with the binding of the cognate amino acids. It was found that they bind rather tightly to the enzyme, the Kd's ranging from 3.1.10(-4) M with glycinol-AMP ester to 3.7.10(-9) M with L-isoleucinol-AMP ester. The binding is not affected by magnesium. It is shown that the free energies of binding of the esters can be calculated adding a constant contribution of the AMP-moiety of about - 4.1 (- 17) kcal/mole (kJ/mole) to the free energies of binding of the cognate amino acids, which we have reported earlier (19, 25, 26).     

Energetic analysis of binding of progesterone and 5 beta-androstane-3,17-dione to anti-progesterone antibody DB3 using molecular dynamics and free energy calculations.

Molecular dynamics simulations and molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) free energy calculations were used to study the energetics of the binding of progesterone (PRG) and 5 beta-androstane-3,17-dione (5AD) to anti-PRG antibody DB3. Although the two steroids bind to DB3 in different orientations, their binding affinities are of the same magnitude, 1 nM for PRG and 8 nM for 5AD. The calculated relative binding free energy of the steroids, 8.8 kJ/mol, is in fair agreement with the experimental energy, 5.4 kJ/mol. In addition, computational alanine scanning was applied to study the role of selected amino acid residues of the ligand-binding site on the steroid cross-reactivity. The electrostatic and van der Waals components of the total binding free energies were found to favour more the binding of PRG, whereas solvation energies were more favourable for the binding of 5AD. The differences in the free energy components are due to the binding of the A rings of the steroids to different binding pockets: PRG is bound to a pocket in which electrostatic antibody-steroid interactions are dominating, whereas 5AD is bound to a pocket in which van der Waals and hydrophobic interactions dominate.     

Free energy simulations and MM-PBSA analyses on the affinity and specificity of steroid binding to antiestradiol antibody

Antiestradiol antibody 57-2 binds 17beta-estradiol (E2) with moderately high affinity (K(a) = 5 x 10(8) M(-1)). The structurally related natural estrogens estrone and estriol as well synthetic 17-deoxy-estradiol and 17alpha-estradiol are bound to the antibody with 3.7-4.9 kcal mol(-1) lower binding free energies than E2. Free energy perturbation (FEP) simulations and the molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) method were applied to investigate the factors responsible for the relatively low cross-reactivity of the antibody with these four steroids, differing from E2 by the substituents of the steroid D-ring. In addition, computational alanine scanning of the binding site residues was carried out with the MM-PBSA method. Both the FEP and MM-PBSA methods reproduced the experimental relative affinities of the five steroids in good agreement with experiment. On the basis of FEP simulations, the number of hydrogen bonds formed between the antibody and steroids, which varied from 0 to 3 in the steroids studied, determined directly the magnitude of the steroid-antibody interaction free energies. One hydrogen bond was calculated to contribute about 3 kcal mol(-1) to the interaction energy. Because the relative binding free energies of estrone (two antibody-steroid hydrogen bonds), estriol (three hydrogen bonds), 17-deoxy-estradiol (no hydrogen bonds), and 17alpha-estradiol (two hydrogen bonds) are close to each other and clearly lower than that of E2 (three hydrogen bonds), the water-steroid interactions lost upon binding to the antibody make an important contribution to the binding free energies. The MM-PBSA calculations showed that the binding of steroids to the antiestradiol antibody is driven by van der Waals interactions, whereas specificity is solely due to electrostatic interactions. In addition, binding of steroids to the antiestradiol antibody 57-2 was compared to the binding to the antiprogesterone antibody DB3 and antitestosterone antibody 3-C4F5, studied earlier with the MM-PBSA method.     

Relative binding free energy calculations of inhibitors to two mutants (Glu46----Ala/Gln) of ribonuclease T1 using molecular dynamics/free energy perturbation approaches.

We present free energy perturbation calculations on the complexes of Glu46----Ala46 (E46A) and Glu46----Gln46 (E46Q) mutants of ribonuclease T1 (RNaseT1) with inhibitors 2'-guanosine monophosphate (GMP) and 2'-adenosine monophosphate (AMP) by a thermodynamic perturbation method implemented with molecular dynamics (MD). Using the available crystal structure of the RNaseT1-GMP complex, the structures of E46A-GMP and E46Q-GMP were model built and equilibrated with MD simulations. The structures of E46A-AMP and E46Q-AMP were obtained as a final structure of the GMP----AMP perturbation calculation respectively. The calculated difference in the free energy of binding (delta delta Gbind) was 0.31 kcal/mol for the E46A system and -1.04 kcal/mol for the E46Q system. The resultant free energies are much smaller than the experimental and calculated value of approximately 3 kcal/mol for the native RNaseT1, which suggests that both mutants have greater relative adenine affinities than native RNaseT1. Especially E46Q is calculated to have a larger affinity for adenine than guanine, as we suggested previously from the calculation on the native RNaseT1. Thus, the molecular dynamics/free energy perturbation method may be helpful in protein engineering, directed toward increasing or changing the substrate specificity of enzymes.     

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 co-factor 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 pK_a 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 μM). 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. © 1993 Wiley-Liss, Inc.     

PSI-DOCK: towards highly efficient and accurate flexible ligand docking

We have developed a new docking method, Pose-Sensitive Inclined (PSI)-DOCK, for flexible ligand docking. An improved SCORE function has been developed and used in PSI-DOCK for binding free energy evaluation. The improved SCORE function was able to reproduce the absolute binding free energies of a training set of 200 protein-ligand complexes with a correlation coefficient of 0.788 and a standard error of 8.13 kJ/mol. For ligand binding pose exploration, a unique searching strategy was designed in PSI-DOCK. In the first step, a tabu-enhanced genetic algorithm with a rapid shape-complementary scoring function is used to roughly explore and store potential binding poses of the ligand. Then, these predicted binding poses are optimized and compete against each other by using a genetic algorithm with the accurate SCORE function to determine the binding pose with the lowest docking energy. The PSI-DOCK 1.0 program is highly efficient in identifying the experimental binding pose. For a test dataset of 194 complexes, PSI-DOCK 1.0 achieved a 67% success rate (RMSD < 2.0 A) for only one run and a 74% success rate for 10 runs. PSI-DOCK can also predict the docking binding free energy with high accuracy. For a test set of 64 complexes, the correlation between the experimentally observed binding free energies and the docking binding free energies for 64 complexes is r = 0.777 with a standard deviation of 7.96 kJ/mol. Moreover, compared with other docking methods, PSI-DOCK 1.0 is extremely easy to use and requires minimum docking preparations. There is no requirement for the users to add hydrogen atoms to proteins because all protein hydrogen atoms and the flexibility of the terminal protein atoms are intrinsically taken into account in PSI-DOCK. There is also no requirement for the users to calculate partial atomic charges because PSI-DOCK does not calculate an electrostatic energy term. These features are not only convenient for the users but also help to avoid the influence of different preparation methods.     

Analysis of the binding energies of testosterone, 5alpha-dihydrotestosterone,androstenedione and dehydroepiandrosterone sulfate with an antitestosterone antibody

Molecular dynamics simulations and molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) free energy calculations were used to study the binding of testosterone (TES), 5alpha-dihydrotestosterone (5ADHT), androstenedione (AND), and dehydroepiandrosterone sulfate (DHEAS) to the monoclonal antitestosterone antibody 3-C(4)F(5). The relative binding free energy of TES and AND was also calculated with free energy perturbation (FEP) simulations. The antibody 3-C(4)F(5) has a relatively high affinity (3 x 10(8) M(-1)) and on overall good binding profile for testosterone but its cross-reactivity with DHEAS has been the main reason for the failure to use this antibody in clinical immunoassays. The relative binding free energies obtained with the MM-PBSA method were 1.5 kcal/mol for 5ADHT, 3.8 kcal/mol for AND, and 4.3 kcal/mol for DHEAS, as compared to TES. When a water molecule of the ligand binding site, observed in the antibody-TES crystal structure, was explicitly included in MM-PBSA calculations, the relative binding energies were 3.4, 4.9, and 5.4 kcal/mol for 5ADHT, AND, and DHEAS, respectively. The calculated numbers are in correct order but larger than the corresponding experimental energies of 1.3, 1.5, and 2.6 kcal/mol, respectively. The fact that the MM-PBSA method reproduced the relative binding free energies of DHEAS, a steroid having a negatively charged sulfate group, and the neutrally charged TES, 5ADHT, and AND in satisfactory agreement with experiment shows the robustness of the method in predicting relative binding affinities. The 800-ps FEP simulations predicted that the antibody 3-C(4)F(5) binds TES 1.3 kcal/mol tighter than AND. Computational mutagenesis of selected amino acid residues of the ligand binding site revealed that the lower affinities of AND and DHEAS as compared to TES are due to a combined effect of several residues, each contributing a small fraction to the tighter binding of TES. An exception to this is Tyr99H, whose mutation to Ala lowered the binding of DHEAS 0.7 kcal/mol more than the binding of TES. This is probably due to the hydrogen bonding interaction formed between the OH group of Tyr99H and the sulfate group of DHEAS. Computational mutagensis data also showed that the affinity of the steroids to the antitestosterone antibody 3-C(4)F(5) would be enhanced if Trp47H were repositioned so that it would make more extensive contacts with the bound ligands. In addition, the binding of steroids to antitestosterone, antiprogesterone, and antiestradiol antibodies is discussed.     

Thermodynamics of water mediating protein-ligand interactions in cytochrome P450cam: a molecular dynamics study.

Ordered water molecules are observed by crystallography and nuclear magnetic resonance to mediate protein-ligand interactions. Here, we examine the energetics of hydrating cavities formed at protein-ligand interfaces using molecular dynamics simulations. The free energies of hydrating two cavities in the active site of two liganded complexes of cytochrome P450cam were calculated by multiconfigurational thermodynamic integration. The complex of cytochrome P450cam with 2-phenyl-imidazole contains a crystallographically well defined water molecule mediating hydrogen bonds between the protein and the inhibitor. We calculate that this water molecule is stabilized by a binding free energy of -11.6 +/- kJ/mol. The complex of cytochrome P450cam with its natural substrate, camphor, contains a cavity that is empty in the crystal structure although a water molecule in it could make a hydrogen bond to camphor. Here, solvation of this cavity is calculated to be unfavorable by +15.8 +/- 5.0 kJ/mol. The molecular dynamics simulations can thus distinguish a hydrated interfacial cavity from an empty one. They also provide support for the notion that protein-ligand complexes can accommodate empty interfacial cavities and that such cavities are likely to be unhydrated unless more than one hydrogen bond can be made to a water molecule in the cavity.     

Quantification and characterization of P-glycoprotein-substrate interactions

It is generally accepted that P-glycoprotein binds its substrates in the lipid phase of the membrane. Quantification and characterization of the lipid-transporter binding step are, however, still a matter of debate. We therefore selected 15 structurally diverse drugs and measured the binding constants from water to the activating (inhibitory) binding region of P-glycoprotein, K(tw(1)) (K(tw(2))), as well as the lipid-water partition coefficients, K(lw). The former were obtained by measuring the concentrations of half-maximum activation (inhibition), K(1) (K(2)), in living NIH-MDR-G185 mouse embryo fibroblasts using a Cytosensor microphysiometer, and the latter were derived from surface activity measurements. This allowed determination of the membrane concentration of drugs at half-maximum P-glycoprotein activation (C(b(1)) = (0.02 to 67) mmol/L lipid), which is much higher than the corresponding aqueous concentration (K(1) = (0.02 to 376) microM). Moreover we determined the free energy of drug binding from water to the activating binding region of the transporter (DeltaG degrees (tw(1)) = (-30 to -54) kJ/mol), the free energy of drug partitioning into the lipid membrane (DeltaG degrees (lw) = (-23 to -34) kJ/mol), and, as the difference of the two, the free energy of drug binding from the lipid membrane to the activating binding region of the transporter (DeltaG degrees (tl(1)) = (-7 to -27) kJ/mol). For the compounds tested DeltaG degrees (tl(1)) was less negative than DeltaG degrees (lw) but varied more strongly. The free energies of substrate binding to the transporter within the lipid phase, DeltaG degrees (tl(1)), are consistent with a modular binding concept, where the energetically most efficient binding module comprises two hydrogen bond acceptor groups.     

Semiquantitative calculations of catalytic free energies in genetically modified enzymes

The catalytic free energy and binding free energies of the native and the Asn-155----Thr, Asn-155----Leu, and Asn-155----Ala mutants of subtilisin are calculated by the empirical valence bond method and a free energy perturbation method. Two simple procedures are used; one "mutates" the substrate, and the other "mutates" the enzyme. The calculated changes in free energies (delta delta G not equal to cat and delta delta Gbind) between the mutant and native enzymes are within 1 kcal/mol of the corresponding observed values. This indicates that we are approaching a quantitative structure-function correlation. The calculated changes in catalytic free energies are almost entirely due to the electrostatic interaction between the enzyme-water system and the charges of the reacting system. This supports the idea that the electrostatic free energy associated with the changes of charges of the reacting system is the key factor in enzyme catalysis.     

A computational analysis of the binding affinities of FKBP12 inhibitors using the MM-PB/SA method

The FK506-binding proteins have been targets of pharmaceutical interests over years. We have studied the binding of a set of 12 nonimmunosuppressive small-molecule inhibitors to FKBP12 through molecular dynamics simulations. Each complex was subjected to 1-ns MD simulation conducted in an explicit solvent environment under constant temperature and pressure. The binding free energy of each complex was then computed by the MM-PB/SA method in the AMBER program. Our MM-PB/SA computation produced a good correlation between the experimentally determined and the computed binding free energies with a correlation coefficient (R(2)) of 0.93 and a standard deviation as low as 0.30 kcal/mol. The vibrational entropy term given by the normal mode analysis was found to be helpful for achieving this correlation. Moreover, an adjustment to one weight factor in the PB/SA model was essential to correct the absolute values of the final binding free energies to a reasonable range. A head-to-head comparison of our MM-PB/SA model with a previously reported Linear Response Approximation (LRA) model suggested that the MM-PB/SA method is more robust in binding affinity prediction for this class of compounds.     

Evaluation of protein-protein association energies by free energy perturbation calculations

The association energy upon binding of different amino acids in the specificity pocket of trypsin was evaluated by free energy perturbation calculations on complexes between bovine trypsin (BT) and bovine pancreatic trypsin inhibitor (BPTI). Three simulations of mutations of the primary binding residue (P(1)) were performed (P(1)-Ala to Gly, P(1)-Met to Gly and P(1)-Met to Ala) and the resulting differences in association energy (DeltaDeltaG(a)) are 2. 28, 5.08 and 2.93 kcal/mol for P(1)-Ala to Gly, P(1)-Met to Gly and to Ala with experimental values of 1.71, 4.62 and 2.91 kcal/mol, respectively. The calculated binding free energy differences are hence in excellent agreement with the experimental binding free energies. The binding free energies, however, were shown to be highly dependent on water molecules at the protein-protein interface and could only be quantitatively estimated if the correct number of such water molecules was included. Furthermore, the cavities that were formed when a large amino acid side-chain is perturbed to a smaller one seem to create instabilities in the systems and had to be refilled with water molecules in order to obtain reliable results. In addition, if the protein atoms that were perturbed away were not replaced by water molecules, the simulations dramatically overestimated the initial state of the free energy perturbations.     

Structure-activity relationship for extracellular block of K+ channels by tetraalkylammonium ions

External tetraalkylammonium ion binding to potassium channels is studied using microscopic molecular modelling methods and the experimental structure of the KcsA channel. Relative binding free energies of the KcsA complexes with Me4N+, Et4N+, and n-Pr4N+ are calculated with the molecular dynamics free energy perturbation approach together with automated ligand docking. The four-fold symmetry of the entrance cavity formed by the Tyr82 residues is found to provide stronger binding for the D2d than for the S4 conformation of the ligands. In agreement with experiment the Et4N+ blocker shows several kcal/mol better binding than the other tetraalkylammonium ions.     

Probing carbohydrate product expulsion from a processive cellulase with multiple absolute binding free energy methods

Understanding the enzymatic mechanism that cellulases employ to degrade cellulose is critical to efforts to efficiently utilize plant biomass as a sustainable energy resource. A key component of cellulase action on cellulose is product inhibition from monosaccharide and disaccharides in the product site of cellulase tunnel. The absolute binding free energy of cellobiose and glucose to the product site of the catalytic tunnel of the Family 7 cellobiohydrolase (Cel7A) of Trichoderma reesei (Hypocrea jecorina) was calculated using two different approaches: steered molecular dynamics (SMD) simulations and alchemical free energy perturbation molecular dynamics (FEP/MD) simulations. For the SMD approach, three methods based on Jarzynski's equality were used to construct the potential of mean force from multiple pulling trajectories. The calculated binding free energies, -14.4 kcal/mol using SMD and -11.2 kcal/mol using FEP/MD, are in good qualitative agreement. Analysis of the SMD pulling trajectories suggests that several protein residues (Arg-251, Asp-259, Asp-262, Trp-376, and Tyr-381) play key roles in cellobiose and glucose binding to the catalytic tunnel. Five mutations (R251A, D259A, D262A, W376A, and Y381A) were made computationally to measure the changes in free energy during the product expulsion process. The absolute binding free energies of cellobiose to the catalytic tunnel of these five mutants are -13.1, -6.0, -11.5, -7.5, and -8.8 kcal/mol, respectively. The results demonstrated that all of the mutants tested can lower the binding free energy of cellobiose, which provides potential applications in engineering the enzyme to accelerate the product expulsion process and improve the efficiency of biomass conversion.     

Dipole-dipole interaction stabilizes the transition state of apurinic/apyrimidinic endonuclease--abasic site interaction

Human apurinic/apyrimidinic (AP) endonuclease (hAPE) initiates the repair of an abasic site (AP site). To gain insight into the mechanisms of damage recognition of hAPE, we conducted surface plasmon resonance spectroscopy to study the thermodynamics and kinetics of its interaction with substrate DNA containing an abasic site (AP DNA). The affinity of hAPE binding toward DNA increased as much as 6-fold after replacing a single adenine (equilibrium dissociation constant, K(D), 5.3 nm) with an AP site (K(D), 0.87 nm). The enzyme-substrate complex formation appears to be thermodynamically stabilized and favored by a large change in Gibbs free energy, DeltaG degrees (-50 kJ/mol). The latter is supported by a high negative change in enthalpy, DeltaH degrees (-43 kJ/mol) and also positive change in entropy, DeltaS degrees (24 J/(K mol)), and thus the binding process is spontaneous at all temperatures. Analysis of kinetic parameters reveals small enthalpy of activation for association, DeltaH degrees++(ass) (-17 kJ/mol), and activation energy for association (E(a), -14 kJ/mol) when compared with the enthalpy of activation for dissociation, DeltaH degrees++(diss) (26 kJ/mol), and activation energy in the reverse direction (E(d), 28 kJ/mol). Furthermore, varying concentration of KCl showed an increase in binding affinity at low concentration but complete abrogation of the binding at higher concentration, implying the importance of hydrophobic, but predominantly ionic, forces in the Michaelis-Menten complex formation. Thus, low activation energy and the enthalpy of activation, which are perhaps a result of dipole-dipole interactions, play critical roles in AP site binding of APE.     

Absolute binding free energy calculations: On the accuracy of computational scoring of protein–ligand interactions

Calculating the absolute binding free energies is a challenging task. Reliable estimates of binding free energies should provide a guide for rational drug design. It should also provide us with deeper understanding of the correlation between protein structure and its function. Further applications may include identifying novel molecular scaffolds and optimizing lead compounds in computer‐aided drug design. Available options to evaluate the absolute binding free energies range from the rigorous but expensive free energy perturbation to the microscopic linear response approximation (LRA/β version) and related approaches including the linear interaction energy (LIE) to the more approximated and considerably faster scaled protein dipoles Langevin dipoles (PDLD/S‐LRA version) as well as the less rigorous molecular mechanics Poisson–Boltzmann/surface area (MM/PBSA) and generalized born/surface area (MM/GBSA) to the less accurate scoring functions. There is a need for an assessment of the performance of different approaches in terms of computer time and reliability. We present a comparative study of the LRA/β, the LIE, the PDLD/S‐LRA/β, and the more widely used MM/PBSA and assess their abilities to estimate the absolute binding energies. The LRA and LIE methods perform reasonably well but require specialized parameterization for the nonelectrostatic term. The PDLD/S‐LRA/β performs effectively without the need of reparameterization. Our assessment of the MM/PBSA is less optimistic. This approach appears to provide erroneous estimates of the absolute binding energies because of its incorrect entropies and the problematic treatment of electrostatic energies. Overall, the PDLD/S‐LRA/β appears to offer an appealing option for the final stages of massive screening approaches. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Energy of binding of Aspergillus oryzae beta-glucosidase with the substrate, and the mechanism of its enzymic action

Several beta-D-glucopyranosides (p-nitrophenyl, phenyl, and ethyl), 1-thio-beta-D-glucopyranosides, and phenyl 2-deoxy, 3-deoxy, 4-deoxy, and 6-deoxy beta-D-glucopyranosides were synthesized and used to study the mechanism of the enzymatic action of Taka-beta-glucosidase [EC 3.2.1.21 Aspergillus oryzae]. Kinetic constants of the enzyme for these glycosides were determined from S/V-S or 1/V-1/S plots, and the hydrolysis rates of these compounds with the enzyme, acid (3 N HCl) and alkali (3 N NaOH) were compared. Inhibition of the enzyme by 1,5-anhydroglucitol, glucal, dihydroglucal, and 1,6-anhydroglucopyranose was also examined. Glucal and 1,5-anhydroglucitol showed strong competitive inhibition. Free energy of binding of each hydroxyl group of glucosidic glucose with the enzyme was estimated from Kms of phenyl beta-glucoside and its deoxy analogues, and also Ki values of some inhibitors. The free energies of binding of 2-OH, 3-OH, 4-OH, and 6-OH were calculated to be 1.1, 2.4, 0.7, and 1.8 kcal/mol, respectively. The free energy of binding of phenoxide at C-1 (0.3 kcal/mol) was calculated from the Km of Ph-beta-Glc and Ki of 1,5-anhydroglucitol. The energy of binding of 5-CH2OH (2.3 kcal/mol) was obtained from the Km of Ph-beta-Glc and that of Ph-beta-Xyl. The sum (6.8 kcal/mol) of each partial binding free energy was close to the value of binding free energy of Ph-beta-Glc (7.0 kcal/mol) calculated by the equation; -delta Gbind = -RT ln Km-T delta Smix, showing that the methods of estimation of each binding energy used in the present study seemed reasonable. Glucal, having a pyranose form distorted slightly, showed strong competitive inhibition and the Ki of this inhibitor was smaller than the Km of Ph-beta-Glc, suggesting that the sugar ring bound to the active site was distored to a half chair form which is labile to acid hydrolysis.     

Empirical calculation of the relative free energies of peptide binding to the molecular chaperone DnaK

We describe a methodology to calculate the relative free energies of protein-peptide complex formation. The interaction energy was decomposed into nonpolar, electrostatic and entropic contributions. A free energy-surface area relationship served to calculate the nonpolar free energy term. The electrostatic free energy was calculated with the finite difference Poisson-Boltzmann method and the entropic contribution was estimated from the loss in the conformational entropy of the peptide side chains. We applied this methodology to a series of DnaK*peptide complexes. On the basis of the single known crystal structure of the peptide-binding domain of DnaK with a bound heptapeptide, we modeled ten other DnaK*heptapeptide complexes with experimentally measured K(d) values from 0.06 microM to 11 microM, using molecular dynamics to refine the structures of the complexes. Molecular dynamic trajectories, after equilibration, were used for calculating the energies with greater accuracy. The calculated relative binding free energies were compared with the experimentally determined free energies. Linear scaling of the calculated terms was applied to fit them to the experimental values. The calculated binding free energies were between -7.1 kcal/mol and - 9.4 kcal/mol with a correlation coefficient of 0.86. The calculated nonpolar contributions are mainly due to the central hydrophobic binding pocket of DnaK for three amino acid residues. Negative electrostatic fields generated by the protein increase the binding affinity for basic residues flanking the hydrophobic core of the peptide ligand. Analysis of the individual energy contributions indicated that the nonpolar contributions are predominant compared to the other energy terms even for peptides with low affinity and that inclusion of the change in conformational entropy of the peptide side chains does not improve the discriminative power of the calculation. The method seems to be useful for predicting relative binding energies of peptide ligands of DnaK and might be applicable to other protein-peptide systems, particularly if only the structure of one protein-ligand complex is available.     

Binding of 1-benzopyran-4-one derivatives to aldose reductase: a free energy perturbation study

The relative binding affinities to human aldose reductase (ALR2) of three new 7-hydroxy-2-benzyl-4H-1-benzopyran-4-one inhibitors were predicted by free energy perturbation (FEP) simulations. Molecular substitutions were specifically designed to investigate the role of hydrogen bonding at the active site of ALR2. Starting from the lead inhibitor 7-hydroxy-2-(4'-hydroxybenzyl)-4H-1-benzopyran-4-one, the 4'-hydroxyl was mutated to methyl and to trifluoromethyl, and an hydroxyl at position 8 was additionally introduced. Once synthesized and tested as inhibitors of ALR2, the compounds displayed variations of K(i) that were in qualitative to quantitative agreement with the calculated relative free energies of binding. The results, discussed in terms of balance between free energies of solvation and free energies of binding to ALR2, elucidate the importance of hydrogen bonding with Thr113 and with Trp111 and cofactor, and provide a rationale to the observed differences in binding affinities.