Polarization effects in protein–ligand calculations extend farther than the actual induction energy

Improving the accuracy of predicting protein-ligand binding-free energy with semiempirical quantum chemistry charge

Aim: It is a challenge to predict binding-free energy (ΔG) accurately. Methodology/results: For accurate ΔG prediction, a new strategy combining solvated interaction energy (SIE) or molecular mechanics/generalized Born surface area (MM/GBSA) approach with the Coulson charge of both protein and ligand calculated by semiempirical quantum mechanics (SQM), named SIE-SQMPC or MM/GBSA-SQMPC approach, was developed and tested on 50 protein-ligand complexes. Both approaches achieved higher correlation (R ²) between experimental and predicted ΔG than that with Amber-ff03 charge, even for ligands with highly different scaffolds. But, SIE-SQMPC is computationally much faster than MM/GBSA-SQMPC. Conclusion: SIE-SQMPC provided an effective alternative to predict ΔG of protein-ligand binding (R ² = 0.66-0.94 for SIE-AM1; R ² = 0.59-0.98 for SIE-PM7), which has the potential of high-throughput processing for molecular docking and drug design.

Thermodynamics calculation of protein-ligand interactions by QM/MM polarizable charge parameters

The calculation of protein-ligand binding free energy (ΔG) is of great importance for virtual screening and drug design. Molecular dynamics (MD) simulation has been an attractive tool to investigate this scientific problem. However, the reliability of such approach is affected by many factors including electrostatic interaction calculation. Here, we present a practical protocol using quantum mechanics/molecular mechanics (QM/MM) calculations to generate polarizable QM protein charge (QMPC). The calculated QMPC of some atoms in binding pockets was obviously different from that calculated by AMBER ff03, which might significantly affect the calculated ΔG. To evaluate the effect, the MD simulations and MM/GBSA calculation with QMPC for 10 protein-ligand complexes, and the simulation results were then compared to those with the AMBER ff03 force field and experimental results. The correlation coefficient between the calculated ΔΔG using MM/GBSA under QMPC and the experimental data is .92, while that with AMBER ff03 force field is .47 for the complexes formed by streptavidin or its mutants and biotin. Moreover, the calculated ΔΔG with QMPC for the complexes formed by ERβ and five ligands is positively related to experimental result with correlation coefficient of .61, while that with AMBER ff03 charge is negatively related to experimental data with correlation coefficient of .42. The detailed analysis shows that the electrostatic polarization introduced by QMPC affects the electrostatic contribution to the binding affinity and thus, leads to better correlation with experimental data. Therefore, this approach should be useful to virtual screening and drug design.

Net charge changes in the calculation of relative ligand-binding free energies via classical atomistic molecular dynamics simulation

The calculation of binding free energies of charged species to a target molecule is a frequently encountered problem in molecular dynamics studies of (bio-)chemical thermodynamics. Many important endogenous receptor-binding molecules, enzyme substrates, or drug molecules have a nonzero net charge. Absolute binding free energies, as well as binding free energies relative to another molecule with a different net charge will be affected by artifacts due to the used effective electrostatic interaction function and associated parameters (e.g., size of the computational box). In the present study, charging contributions to binding free energies of small oligoatomic ions to a series of model host cavities functionalized with different chemical groups are calculated with classical atomistic molecular dynamics simulation. Electrostatic interactions are treated using a lattice-summation scheme or a cutoff-truncation scheme with Barker-Watts reaction-field correction, and the simulations are conducted in boxes of different edge lengths. It is illustrated that the charging free energies of the guest molecules in water and in the host strongly depend on the applied methodology and that neglect of correction terms for the artifacts introduced by the finite size of the simulated system and the use of an effective electrostatic interaction function considerably impairs the thermodynamic interpretation of guest-host interactions. Application of correction terms for the various artifacts yields consistent results for the charging contribution to binding free energies and is thus a prerequisite for the valid interpretation or prediction of experimental data via molecular dynamics simulation. Analysis and correction of electrostatic artifacts according to the scheme proposed in the present study should therefore be considered an integral part of careful free-energy calculation studies if changes in the net charge are involved.