| 1. |
- Ciancetta, Antonella, et al.
(författare)
-
A QM/MM study of the binding of RAPTA ligands to cathepsin B
- 2011
-
Ingår i: Journal of Computer-Aided Molecular Design. - Springer. - 0920-654X. ; 25:8, s. 729-742
-
Tidskriftsartikel (refereegranskat)abstract
- We have carried out quantum mechanical (QM) and QM/MM (combined QM and molecular mechanics) calculations, as well as molecular dynamics (MD) simulations to study the binding of a series of six RAPTA (Ru(II)-arene-1,3,5-triaza-7-phosphatricyclo-3.3.1.1 decane) complexes with different arene substituents to cathepsin B. The recently developed QM/MM-PBSA approach (QM/MM combined with Poisson-Boltzmann solvent-accessible surface area solvation) has been used to estimate binding affinities. The QM calculations reproduce the antitumour activities of the complexes with a correlation coefficient (r(2)) of 0.35-0.86 after a conformational search. The QM/MM-PBSA method gave a better correlation (r(2) = 0.59) when the protein was fixed to the crystal structure, but more reasonable ligand structures and absolute binding energies were obtained if the protein was allowed to relax, indicating that the ligands are strained when the protein is kept fixed. In addition, the best correlation (r(2) = 0.80) was obtained when only the QM energies were used, which suggests that the MM and continuum solvation energies are not accurate enough to predict the binding of a charged metal complex to a charged protein. Taking into account the protein flexibility by means of MD simulations slightly improves the correlation (r(2) = 0.91), but the absolute energies are still too large and the results are sensitive to the details in the calculations, illustrating that it is hard to obtain stable predictions when full flexible protein is included in the calculations.
|
|
| 2. |
- Diehl, Carl, et al.
(författare)
-
Conformational entropy changes upon lactose binding to the carbohydrate recognition domain of galectin-3.
- 2009
-
Ingår i: Journal of biomolecular NMR. - Springer Netherlands. - 1573-5001. ; 45:1-2, s. 157-169
-
Tidskriftsartikel (refereegranskat)abstract
- The conformational entropy of proteins can make significant contributions to the free energy of ligand binding. NMR spin relaxation enables site-specific investigation of conformational entropy, via order parameters that parameterize local reorientational fluctuations of rank-2 tensors. Here we have probed the conformational entropy of lactose binding to the carbohydrate recognition domain of galectin-3 (Gal3), a protein that plays an important role in cell growth, cell differentiation, cell cycle regulation, and apoptosis, making it a potential target for therapeutic intervention in inflammation and cancer. We used (15)N spin relaxation experiments and molecular dynamics simulations to monitor the backbone amides and secondary amines of the tryptophan and arginine side chains in the ligand-free and lactose-bound states of Gal3. Overall, we observe good agreement between the experimental and computed order parameters of the ligand-free and lactose-bound states. Thus, the (15)N spin relaxation data indicate that the molecular dynamics simulations provide reliable information on the conformational entropy of the binding process. The molecular dynamics simulations reveal a correlation between the simulated order parameters and residue-specific backbone entropy, re-emphasizing that order parameters provide useful estimates of local conformational entropy. The present results show that the protein backbone exhibits an increase in conformational entropy upon binding lactose, without any accompanying structural changes.
|
|
| 3. |
- Diehl, Carl, et al.
(författare)
-
Protein Flexibility and Conformational Entropy in Ligand Design Targeting the Carbohydrate Recognition Domain of Galectin-3.
- 2010
-
Ingår i: Journal of the American Chemical Society. - American Chemical Society. - 1520-5126. ; 132, s. 14577-14589
-
Tidskriftsartikel (refereegranskat)abstract
- Rational drug design is predicated on knowledge of the three-dimensional structure of the protein-ligand complex and the thermodynamics of ligand binding. Despite the fundamental importance of both enthalpy and entropy in driving ligand binding, the role of conformational entropy is rarely addressed in drug design. In this work, we have probed the conformational entropy and its relative contribution to the free energy of ligand binding to the carbohydrate recognition domain of galectin-3. Using a combination of NMR spectroscopy, isothermal titration calorimetry, and X-ray crystallography, we characterized the binding of three ligands with dissociation constants ranging over 2 orders of magnitude. (15)N and (2)H spin relaxation measurements showed that the protein backbone and side chains respond to ligand binding by increased conformational fluctuations, on average, that differ among the three ligand-bound states. Variability in the response to ligand binding is prominent in the hydrophobic core, where a distal cluster of methyl groups becomes more rigid, whereas methyl groups closer to the binding site become more flexible. The results reveal an intricate interplay between structure and conformational fluctuations in the different complexes that fine-tunes the affinity. The estimated change in conformational entropy is comparable in magnitude to the binding enthalpy, demonstrating that it contributes favorably and significantly to ligand binding. We speculate that the relatively weak inherent protein-carbohydrate interactions and limited hydrophobic effect associated with oligosaccharide binding might have exerted evolutionary pressure on carbohydrate-binding proteins to increase the affinity by means of conformational entropy.
|
|
| 4. |
- Genheden, Samuel, et al.
(författare)
-
A comparison of different initialization protocols to obtain statistically independent molecular dynamics simulations.
- 2011
-
Ingår i: Journal of computational chemistry. - Wiley Periodicals. - 1096-987X. ; 32:2, s. 187-195
-
Tidskriftsartikel (refereegranskat)abstract
- We study how the results of molecular dynamics (MD) simulations are affected by various choices during the setup, e.g., the starting velocities, the solvation, the location of protons, the conformation of His, Asn, and Gln residues, the protonation and titration of His residues, and the treatment of alternative conformations. We estimate the binding affinity of ligands to four proteins calculated with the MM/GBSA method (molecular mechanics combined with a generalized Born and surface area solvation energy). For avidin and T4 lysozyme, all variations gave similar results within 2 kJ/mol. For factor Xa, differences in the solvation or in the selection of alternative conformations gave results that are significantly different from those of the other approaches by 4-6 kJ/mol, whereas for galectin-3, changes in the conformations, rotations, and protonation gave results that differed by 10 kJ/mol, but only if residues close to the binding site were modified. This shows that the results of MM/GBSA calculations are reasonably reproducible even if the MD simulations are set up with different software. Moreover, we show that the sampling of phase space can be enhanced by solvating the systems with different equilibrated water boxes, in addition to the common use of different starting velocities. If different conformations are available in the crystal structure, they should also be employed to enhance the sampling. Protonation, ionization, and conformations of Asn, Gln, and His may also be used to enhance sampling, but great effort should be spent to obtain as reliable predictions as possible close to the active site.
|
|
| 5. |
- Genheden, Samuel, et al.
(författare)
-
Accurate Predictions of Nonpolar Solvation Free Energies Require Explicit Consideration of Binding-Site Hydration
- 2011
-
Ingår i: Journal of the American Chemical Society. - American Chemical Society. - 0002-7863. ; 133:33, s. 13081-13092
-
Tidskriftsartikel (refereegranskat)abstract
- Continuum solvation methods are frequently used to increase the efficiency of computational methods to estimate free energies. In this paper, we have evaluated how well such methods estimate the nonpolar solvation free-energy change when a ligand binds to a protein. Three different continuum methods at various levels of approximation were considered, viz., the polarized continuum model (PCM), a method based on cavity and dispersion terms (CD), and a method based on a linear relation to the solvent-accessible surface area (SASA). Formally rigorous double-decoupling thermodynamic integration was used as a benchmark for the continuum methods. We have studied four protein-ligand complexes with binding sites of varying solvent exposure, namely the binding of phenol to ferritin, a biotin analogue to avidin, 2-aminobenzimidazole to trypsin, and a substituted galactoside to galectin-3. For ferritin and avidin, which have relatively hidden binding sites, rather accurate nonpolar solvation free energies could be obtained with the continuum methods if the binding site is prohibited to be filled by continuum water in the unbound state, even though the simulations and experiments show that the ligand replaces several water molecules upon binding. For the more solvent exposed binding sites of trypsin and galectin-3, no accurate continuum estimates could be obtained, even if the binding site was allowed or prohibited to be filled by continuum water. This shows that continuum methods fail to give accurate free energies on a wide range of systems with varying solvent exposure because they lack a microscopic picture of binding-site hydration as well as information about the entropy of water molecules that are in the binding site before the ligand binds. Consequently, binding affinity estimates based upon continuum solvation methods will give absolute binding energies that may differ by up to 200 kJ/mol depending on the method used. Moreover, even relative energies between ligands with the same scaffold may differ by up to 75 kJ/mol. We have tried to improve the continuum solvation methods by adding information about the solvent exposure of the binding site or the hydration of the binding site, and the results are promising at least for this small set of complexes.
|
|
| 6. |
- Genheden, Samuel, et al.
(författare)
-
An MM/3D-RISM approach for ligand binding affinities.
- 2010
-
Ingår i: Journal of physical chemistry. B. - American Chemical Society. - 1520-5207. ; 114:25, s. 8505-8516
-
Tidskriftsartikel (refereegranskat)abstract
- We have modified the popular MM/PBSA or MM/GBSA approaches (molecular mechanics for a biomolecule, combined with a Poisson-Boltzmann or generalized Born electrostatic and surface area nonelectrostatic solvation energy) by employing instead the statistical-mechanical, three-dimensional molecular theory of solvation (also known as 3D reference interaction site model, or 3D-RISM-KH) coupled with molecular mechanics or molecular dynamics ( Blinov , N. ; et al. Biophys. J. 2010 ; Luchko , T. ; et al. J. Chem. Theory Comput. 2010 ). Unlike the PBSA or GBSA semiempirical approaches, the 3D-RISM-KH theory yields a full molecular picture of the solvation structure and thermodynamics from the first principles, with proper account of chemical specificities of both solvent and biomolecules, such as hydrogen bonding, hydrophobic interactions, salt bridges, etc. We test the method on the binding of seven biotin analogues to avidin in aqueous solution and show it to work well in predicting the ligand-binding affinities. We have compared the results of 3D-RISM-KH with four different generalized Born and two Poisson-Boltzmann methods. They give absolute binding energies that differ by up to 208 kJ/mol and mean absolute deviations in the relative affinities of 10-43 kJ/mol.
|
|
| 7. |
- Genheden, Samuel
(författare)
-
Are homology models sufficiently good for free-energy simulations?
- 2012
-
Ingår i: Journal of chemical information and modeling. - American Chemical Society. - 1549-960X. ; 52:11, s. 3013-3021
-
Tidskriftsartikel (refereegranskat)abstract
- In this paper, I evaluate the usefulness of protein homology models in rigorous free-energy simulations to determine ligand affinities. Two templates were used to create models of the factor Xa protein and one template was used for dihydrofolate reductase from Plasmodium falciparum. Then, the relative free energies for several pairs of ligands were estimated using thermodynamic integration with the homology models as starting point of the simulation. These binding affinities were compared to affinities obtained when using published crystal structures as starting point of the simulations. Encouragingly, the differences between the affinities obtained when starting from either homology models or crystal structure were not statistical significant for a majority of the considered pairs of ligands. Differences between 1 and 2 kJ/mol were observed for the dihydrofolate reductase ligands and differences between 0 and 8 kJ/mol were observed for the factor Xa ligands. The largest difference for factor Xa was caused by an erroneous modeling of a loop region close to two of the ligands, and it was only observed when using one of the templates. Therefore, it is advisible to always use more than one template when creating homology models if they should be used in free-energy simulations.
|
|
| 8. |
- Genheden, Samuel, et al.
(författare)
-
Binding Affinities of Factor Xa Inhibitors Estimated by Thermodynamic Integration and MM/GBSA.
- 2011
-
Ingår i: Journal of Chemical Information and Modeling. - American Chemical Society. - 1549-960X. ; 51:Online March 18, 2011, s. 947-958
-
Tidskriftsartikel (refereegranskat)abstract
- We present free energy estimates of nine 3-amidinobenzyl-1H-indole-2-carboxamide inhibitors of factor Xa. Using alchemical thermodynamic integration (TI) calculations, we estimate the difference in binding free energies with high accuracy and precision, except for mutations involving one of the amidinobenzyl rings. Crystal studies show that the inhibitors may bind in two distinct conformations, and using TI, we show that the two conformations give a similar binding affinity. Furthermore, we show that we can reduce the computational demand, while still retaining a high accuracy and precision, by using fewer integration points and shorter protein-ligand simulations. Finally, we have compared the TI results to those obtained with the simpler MM/GBSA method (molecular-mechanics with generalized Born surface-area solvation). MM/GBSA gives better results for the mutations that involve a change of net charge, but if a precision similar to that of the TI method is required, the MM/GBSA method is actually slightly more expensive. Thus, we have shown that TI could be a valuable tool in drug design.
|
|
| 9. |
- Genheden, Samuel, et al.
(författare)
-
Comparison of end-point continuum-solvation methods for the calculation of protein-ligand binding free energies.
- 2012
-
Ingår i: Proteins. - Wiley-Blackwell Publishing, Inc. - 1097-0134. ; 80:5, s. 1326-1342
-
Tidskriftsartikel (refereegranskat)abstract
- We have compared the predictions of ligand-binding affinities from several methods based on end-point molecular dynamics simulations and continuum solvation, i.e. methods related to MM/PBSA (molecular mechanics combined with Poisson-Boltzmann and surface area solvation). Two continuum-solvation models were considered, viz. the Poisson-Boltzmann (PB) and generalised Born (GB) approaches. The non-electrostatic energies were also obtained in two different ways, viz. either from the sum of the bonded, van der Waals, non-polar solvation energies, and entropy terms (as in MM/PBSA), or from the scaled protein-ligand van der Waals interaction energy (as in the linear interaction energy approach, LIE). Three different approaches to calculate electrostatic energies were tested, viz. the sum of electrostatic interaction energies and polar solvation energies, obtained either from a single simulation of the complex or from three independent simulations of the complex, the free protein, and the free ligand, or the linear-response approximation (LRA). Moreover, we investigated the effect of scaling the electrostatic interactions by an effective internal dielectric constant of the protein (ε(int) ). All these methods were tested on the binding of seven biotin analogues to avidin and nine 3-amidinobenzyl-1H-indole-2-carboxamide inhibitors to factor Xa. For avidin, the best results were obtained with a combination of the LIE non-electrostatic energies with the MM+GB electrostatic energies from a single simulation, using ε(int) = 4. For fXa, standard MM/GBSA, based on one simulation and using ε(int) = 4-10 gave the best result. The optimum internal dielectric constant seems to be slightly higher with PB than with GB solvation. Proteins 2012. © 2012 Wiley-Liss, Inc.
|
|
| 10. |
- Genheden, Samuel, et al.
(författare)
-
Comparison of the Efficiency of the LIE and MM/GBSA Methods to Calculate Ligand-Binding Energies
- 2011
-
Ingår i: Journal of Chemical Theory and Computation. - American Chemical Society. - 1549-9618. ; 7:11, s. 3768-3778
-
Tidskriftsartikel (refereegranskat)abstract
- We have evaluated the efficiency of two popular end-point methods to calculate ligand-binding free energies, LIE (linear interaction energy) and MM/GBSA (molecular mechanics with generalized Born surface-area solvation), i.e. the computational effort needed to obtain estimates of a similar precision. As a test case, we use the binding of seven biotin analogues to avidin. The energy terms used by MM/GBSA and LIE exhibit a similar correlation time (similar to 5 ps), and the equilibration time seems also to be similar, although it varies much between the various ligands. The results show that the LIE method is more effective than MM/GBSA, by a factor of 2-7 for a truncated spherical system with a radius of 26 angstrom and by a factor of 1.0-2.4 for the full avidin tetramer (radius 47 angstrom). The reason for this is the cost for the MM/GBSA entropy calculations, which more than compensates for the extra simulation of the free ligand in LIE. On the other hand, LIE requires that the protein is neutralized, whereas MM/GBSA has no such requirements. Our results indicate that both the truncation and neutralization of the proteins may slow the convergence and emphasize small differences in the calculations, e.g., differences between the four subunits in avidin. Moreover, LIE cannot take advantage of the fact that avidin is a tetramer. For this test case, LIE gives poor results with the standard parametrization, but after optimizing the scaling factor of the van der Waals terms, reasonable binding affinities can be obtained, although MM/GBSA still gives a significantly better predictive index and correlation to the experimental affinities.
|
|
