Next Article in Journal
A Novel Prostate-Specific Membrane-Antigen (PSMA) Targeted Micelle-Encapsulating Wogonin Inhibits Prostate Cancer Cell Proliferation via Inducing Intrinsic Apoptotic Pathway
Next Article in Special Issue
RVMAB: Using the Relevance Vector Machine Model Combined with Average Blocks to Predict the Interactions of Proteins from Protein Sequences
Previous Article in Journal
Correction: Chan, Y.-Y., et al. The Constituents of Michelia compressa var. formosana and Their Bioactivities. Int. J. Mol. Sci. 2014, 15, 10926–10935
Previous Article in Special Issue
Confocal Spectroscopy to Study Dimerization, Oligomerization and Aggregation of Proteins: A Practical Guide
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 17
Issue 5
10.3390/ijms17050742
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Recent Progress in Treating Protein–Ligand Interactions with Quantum-Mechanical Methods

by
Nusret Duygu Yilmazer
and
Martin Korth
*
Institute for Theoretical Chemistry, Ulm University, Albert-Einstein-Allee 11, 89069 Ulm, Germany
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2016, 17(5), 742; https://doi.org/10.3390/ijms17050742
Submission received: 16 March 2016 / Revised: 18 April 2016 / Accepted: 3 May 2016 / Published: 16 May 2016
(This article belongs to the Special Issue Proteins and Protein-Ligand Interactions)
Versions Notes

Abstract

:
We review the first successes and failures of a “new wave” of quantum chemistry-based approaches to the treatment of protein/ligand interactions. These approaches share the use of “enhanced”, dispersion (D), and/or hydrogen-bond (H) corrected density functional theory (DFT) or semi-empirical quantum mechanical (SQM) methods, in combination with ensemble weighting techniques of some form to capture entropic effects. Benchmark and model system calculations in comparison to high-level theoretical as well as experimental references have shown that both DFT-D (dispersion-corrected density functional theory) and SQM-DH (dispersion and hydrogen bond-corrected semi-empirical quantum mechanical) perform much more accurately than older DFT and SQM approaches and also standard docking methods. In addition, DFT-D might soon become and SQM-DH already is fast enough to compute a large number of binding modes of comparably large protein/ligand complexes, thus allowing for a more accurate assessment of entropic effects.

1. Introduction and Overview

The accurate and robust computational prediction of the relative strength of protein/ligand interactions is at the heart of virtual drug design [1,2]. A computational strategy that is able to evaluate protein/ligand binding affinities reliably and quickly over a wide range of protein targets and drug compounds is, unfortunately, still elusive [3,4,5]. This problem stems from the fact that free energies of binding for large molecular systems have to be computed to very high accuracy. Enthalpic contributions can now be covered, even at the quantum mechanical level, but the estimation of entropic effects is still complicated at all levels of theory [6,7]. A range of algorithms are available to either estimate or compute these contributions from classical molecular dynamics (MD) simulations [8], but a number of recent studies have shown that the accuracy of the underlying potential energy hypersurfaces is often at least questionable [9,10]. This is problematic, because very high accuracies are needed, and delicate enthalpy/entropy compensation effects have to be taken into account [11]. If quantum mechanical methods are applied, the problem lies in identifying the most important binding modes, as the number of modes that can be subjected to calculations is limited by computing power. In the past, researchers had to choose between either covering the dynamics of the system at hand while relying on very approximate classical (molecular mechanics, MM) energy models, or using a highly accurate energy model for a single binding mode, thus effectively ignoring the dynamics of the problem. It should be noted that thorough investigations of the role of entropic contributions are of course only possible with accurate energy models, and thus critically rely on the availability to do quantum mechanical calculations on a large number of binding modes, at least for model systems [12,13].
More advanced approaches can, on the other hand, build upon existing models. Sophisticated algorithms are available for the identification of possible binding modes—i.e., the “posing” step of standard “docking” procedures. These seem to work rather well, unlike the evaluation of binding energies—i.e., the “scoring” step [14]. Alternatively, snapshot structures can be extracted from classical molecular dynamics simulations [8,15]. A very interesting solution was put forward by Zheng et al. with their “movable type” method [16]. First, interaction pairs are constructed, then two tables with interaction strengths and interaction probabilities are constructed, and binding free energies are finally computed by applying statistical mechanics using these two tables. The same approach was developed in the past by Klamt for describing liquid phase properties with his COSMO-RS theory [17]. The most attractive features of this approach are that no explicit sampling is necessary and that pair-wise interactions can be computed at arbitrary levels of theory.
Free energies of binding can thus be calculated with quantum mechanical approaches from weighted averaging over free energies of pre-selected binding modes (or pair-wise interactions), taking into account the entropic effects of having a range of modes accessible for the system [18]. For details on ensemble averaging calculations in the context of non-covalent interactions and free energies of solvation, see recent reviews by Gilson [6] and Jensen [19]. The most important point is that entropic contributions are split into a part that can be computed for single binding modes using normal mode analysis and a part that stems from the fact that the system can change between single binding modes and which is related to the relative energies of these binding modes. While dispersion-corrected density functional theory (DFT-D) methods are still too slow for normal mode analysis of large biomolecules, dispersion and hydrogen bond-corrected semi-empirical quantum mechanical (SQM-DH) methods now allow for such calculations [20]. Free energies of the single binding modes can be computed within standard rigid-rotor, harmonic-oscillator approximations from energy and vibrational frequency calculations of the complex, protein, and ligand structures. An excellent overview of computing binding free energies of host–guest and similar systems in aqueous solution was given by Jensen [19], who additionally discusses possible error sources, including the role of explicit solvation and counter-ions. One advantage of quantum mechanical (QM) level approaches over classical methods is the correct treatment of polarization effects. This has the important side-effect that (continuum) solvation model results are usually already much improved for SQM, but especially for DFT methods [19]. “Dispersion corrected” density functional theory (DFT-D) methods have now become the de facto standard for larger molecular systems that are dominated by non-covalent interactions like protein/ligand complexes [21,22]. DFT-D was shown to perform as accurately as high-level quantum chemistry methods for such problems and is routinely applicable to systems with several hundreds of atoms. Uncorrected DFT does not account for London dispersion (dispersive van-der-Waals) interactions, which became a problem of the utmost importance when larger bio and also nano-structures came into reach during the late 1990s. The systematic application of dispersion corrections terms by Stefan Grimme [21,22] turned out to be a very efficient way of dealing with this problem. DFT methods are not yet fast enough for overall free energy calculations of more than model systems, but they can serve as a reference for faster methods, and they can give important insight into the balance of interactions, especially also concerning the interplay of enthalpic and entropic parts.
Very recently, dispersion- and also hydrogen-bond (H)-corrected semi-empirical quantum mechanical (SQM-DH) methods were developed, which perform almost as accurately as DFT-D for non-covalent interactions [23,24]. SQM-DH methods are about three orders of magnitude faster than DFT-D, which makes them very promising candidates to compute accurate free energies of binding for protein/ligand systems. In the following we will give an overview of recent developments in the field of applying SQM-DH, DFT-D, and higher-level quantum chemistry methods to protein/ligand systems or related model systems. A larger number of older (semi-empirical) quantum mechanical studies are available, but they usually apply “uncorrected” DFT or SQM methods [14,25], which are now known to perform very unreliably for non-covalent interactions, and are therefore not reviewed here. Finally, it needs to be mentioned that the overall problem of designing new drugs is of course a much more complex one. Ideas for the working mechanism of specific drugs have to be developed, target proteins have to be identified, structural information about proteins has to be collected, and input has to be prepared so that computational models can be applied. The problem of how to arrive at good model systems for calculations can easily become as complicated as actually performing calculations. However, similar to the case of enthalpic and entropic contributions, thorough investigations of error sources again have to be based on the application of trustworthy energy models.

2. Dispersion-Corrected Density Functional Theory (DFT-D) and Higher-Level Quantum Chemistry Investigations

2.1. Host/Guest and Other Model Systems

Most DFT-D and higher-level wave function theory (WFT)-based work in the field has been done on model systems and especially host/guest complexes, which are like protein/ligand complexes dominated by non-covalent interactions. The Grimme group has published not only a large number of studies on the role of dispersion interactions in chemistry (see for instance [26]), but has also produced a number of very thorough studies on how to compute host/guest interaction energies in a reliable fashion. Concerning electronic energies, Antony et al. [27] benchmarked generalized gradient approximations (GGA) and hybrid DFT-D3 methods against approximate Coupled Cluster data at the LPNO-CEPA-1/CBS (local pair natural orbital coupled electron pair approximation/complete basis set) level for protein/ligand model systems and found typical errors of 1–2 kcal (or 5%) with respect to theoretical reference data and 2–4 kcal/mol (or 10%) with respect to the experimental references. Notably, second order Møller–Plesset perturbation method (MP2) was found to perform less accurately and slower than DFT-D methods. Antony and Grimme [28] also showed how DFT-D3 can be applied to complexes of up to a few thousand atoms via molecular fragmentation with conjugated caps (MFCC). For MFCC (further discussed in references [29,30,31]), the protein is decomposed into fragments, which are then capped with typical end groups. MFCC is found to work well for the artificial case of a neutral peptide, but convergence problems are found for fragments with multiple charges. Grimme [12] then turned to go from enthalpies to free energies of binding for host/guest systems, and found similar overall errors when additionally applying implicit solvation models and lower-level estimates for the dynamical terms. Further studies from the Grimme group followed, an overview is given in Antony et al. [32]. They found that free energies of binding for host/guest systems can typically be computed up to an accuracy of about 2 kcal/mol if a single binding mode is dominating, and emphasize the critical role of dispersion interactions, so that either DFT-D with triple-zeta basis sets or Coupled Cluster (rather than MP2) with very large basis set need to be employed.
Skylaris and co-workers follow a somewhat complementary route: they perform “linear-scaling” plain-wave (periodic boundary condition) DFT calculations, which in principle allows them to go to very large system sizes, but limits the applicability of many existing techniques from computational chemistry [33,34,35,36]. However, they also work on using electrostatic embedding techniques, for instance [33]. Cole et al. [34] investigate protein/protein interactions, which are of fundamental importance to biological processes. Protein/protein interactions are also a very good example of when periodic DFT calculations might become advantageous, as the relevant systems cannot be easily broken down into smaller “cluster-based” calculations. Looking at the interaction of the breast cancer suppressor BRCA2 with RAD51, they want to understand how the former binds and delivers the latter to sites of DNA damage for repair. The lack of experimental structural information is addressed with model system calculations at the force field and DFT-D level. Periodic DFT calculations are of course no remedy for the problem of how to take into account entropic contributions, and poor agreement with experiment is accordingly found in cases where these contributions matter [35]. Fox et al. [36] ran into the same problem when investigating a protein/model system of more than 2500 atoms for the interaction of eight small aromatic ligands with the T4-lysozyme double mutant L99A/M102Q protein polar binding site.
The Hobza group has published a large body of work on the accurate calculation of non-covalent interactions and the benchmarking of SQM vs. DFT vs. WFT methods for this purpose (see for instance reference [20]). For a number of model systems, theoretical references are now available up to a full configuration interaction level [37]. Other benchmark work includes the extensive sampling of the conformational space of nine HIV-1 protease inhibitors [38], where 1600 conformations of each ligand are investigated at the SQM-DH level, and errors from using single conformations and short molecular dynamics runs are about 3 kcal/mol. Further studies address the effects of using implicit solvent models and conformational ensembles [39], where the multiconfigurational COSMO-RS is found to perform best. For modeling protein/ligand interactions, Hobza and co-workers most often turn to SQM-DH methods, which is why this work is discussed in the following section.
Ryde and co-workers have published extensively on using QM/MM type approaches for treating protein/ligand interactions, but (S)QM methods in the older works “uncorrected” QM methods and these studies will accordingly not be referenced here. One older investigation was using MP2 [40], which again cannot be recommended for non-covalent interactions. Newer work is partly based on “enhanced” methods and also very high-level WFT approaches. Mikulskis et al. [41] for instance computed free energies of binding for nine cyclic carboxylate guest molecules to an octa-acid host. The best agreement with experiment was found at the classical level. Results were worse with DFT, and again worse with their linearized coupled-cluster theory with single and double excitations LCCSD (T0) approximate Coupled Cluster approach. Andrejic et al. arrived at LCCSD (T0)/CBS calculations for a 200 atom host/guest model system [42] by calculating short-range interactions pairwise at the QM level and treating long-range interactions with a polarizable multipole approach. This approach gives an accuracy of 2–3 kcal/mol; i.e., comparable to a proper DFT-D treatment, but almost no correlation with experiment is observed, as only a single binding mode could be treated.
Protein/ligand-screenings beyond proof-of-concept studies on model systems [43,44,45,46,47,48,49] are not yet done with DFT or higher-level WFT methods, but most often SQM-DH methods—see the next section for details. It was nevertheless already shown that QM methods are advantageous for protein structure refinement [50,51,52,53]. Christensen et al. [54] were even able to use backbone amide proton chemical shifts calculated at the QM-level for this purpose [29,55].

2.2. Local Wave Function Theory Methods

A very interesting and important development for the future of high-level quantum mechanical calculations on large molecular systems (and therefore also protein/ligand complexes) is the advent of robust local coupled cluster (CC) methods, and especially local pair natural orbital (LPNO) approaches [56]. LPNO-CC methods lose very little (usally substantially less than 1.0 kcal/mol) upon canonical CC in terms of accuracy, but can be applied to systems that are normally only in the reach of DFT. Results furthermore converge systematically to the canonical CC ones if tighter thresholds are chosen, which of course also increases computational costs. LPNO-CC methods furthermore show an improved basis set convergence, leading to higher gains with larger basis sets, thus bringing several hundred atoms in reach of highly accurate quantum chemistry methods. The most imporant feature is that LNPO-CC methods can (unlike other local CC approaches) be applied to a wide range of problems in a black box fashion, which makes it in principle suitable for re-screening purposes. LPNO-CCSD [57] is based on earlier work on LNPO coupled electron pair approximation (CEPA) methods [58], and was then re-designed using pre-screening, truncation, and domain techniques, which lead to the DLPNO-CCSD method. The advantage of this ansatz to exploit locality for CC is that PNOs allow for very compact wave functions, which in turn allow for large and therefore reliable domain approximations. First benchmark calculations showed that DLPNO-CCSD performs very similarly to LNPO-CCSD and thus canonical CCSD, but becomes substantially faster for medium to large systems, as now CPU time, as well as RAM and disk space resources scale nearly linearly with system size.
Liakos et al. [59] thoroughly benchmarked DLPNO-CCSD(T) for a broad range of chemical applications. They found that default thresholds (three parameters overall) allow for an accuracy of 1 kcal/mol or better, while tighter thresholds can reduce errors to about 0.25 kcal/mol. Most importantly for the case of protein/ligand interactions, default thresholds already lead to average errors as small as 0.3 kcal/mol for non-covalent interactions. Excellent results are also achieved for linear versus folded conformations of unbranched alkanes [60], which could be considered as the simplest conceivable atomic-scale model system for protein folding. As in many cases, the CC step in DLNPO-CCSD is not more costly than the Hartree-Fock (HF) step; Liakos and Neese [61] benchmarked DLPNO-CCSD(T) against DFT methods with respect to both accuracy and computational cost. Even with loose thresholds, DLPNO-CCSD(T) is more accurate than any DFT approach, while computational costs for aug-cc-pVTZ basis sets are roughly similar to B3LYP. Errors below 1 kcal/mol are achieved with default thresholds, computation times become two times higher. CCSD(T)/CBS level accuracy is achieved with tight thresholds, while computation times are around only four times higher than B3LYP. The scaling of the new approach is furthermore like HF or DFT dominated by the self-consistent field equations, which is quadratic to cubic, but can be made linear scaling for large systems. Convergence with basis set size is nevertheless somewhat slower for DLNP-CCSD(T) than DFT. As several of the above-listed benchmark studies included model systems for biomolecular interactions, it seems very likely that the availability of fast and extremely accurate WFT methods will have a strong impact on the quantum-level re-scoring of protein/ligand interactions in the near future.

3. Dispersion and Hydrogen Bond-Corrected Semi-Empirical Quantum Mechanical (SQM-DH) Applications

3.1. Model Systems and Related Studies

A large number of studies have now shown that SQM-DH methods perform almost as accurately as DFT-D for a wide variety of systems that are dominated by non-covalent interactions [62,63,64,65,66,67,68,69,70,71,72,73,74,75]. Furthermore, detailed benchmark studies, in comparison to classical approaches, DFT-D, and highly accurate WFT data is available in the original development papers as well as in separate benchmark studies [76]. This includes investigations into hydrogen-bond cooperativity effects [77] and halogen-bonded systems [78]. Several SQM-DH models exist, but it is important to note that the main difference for non-covalent interactions comes from including DH correction, not the actual model one is using [79]. Methods include AM1-FS1 [80,81], PM6-DH2 [82], and PM6-DH+ [83]. PM6-DH2 and PM6-DH+ were further improved for stable molecular dynamics simulations, leading to PM6-D3H4 [84] and arguably the most advanced model, PM6-D3H+ [85]. See references [23,24] for details on SQM-DH methods. For complexes of hundreds of atoms, SQM-DH methods are fast enough for the calculation of hundreds of millions of configurations. PM6-DH+ for instance was applied for the calculation of 320 million interaction energies for a peptide/ligand system [86]. Excellent correlation with experiment was observed, with an R2 value of 0.9, in comparison to 0.7 for an classical approach. Like DFT-D, SQM-DH methods were shown to be valuable for structure refinement purposes [75,87,88,89,90]. Again, like DFT-D, SQM-DH proved to be an intesting tool for applications in materials science [71,91,92,93,94,95], giving some hope for tackling problems at the interface of life and materials science.

3.2. Host/Guest Systems and Protein/Ligand Interactions

3.2.1. Host/Guest Systems

The most important model systems with relation to protein/ligand interactions are host/guest systems, for which SQM-DH methods are now routinely applied [12,96,97,98,99]. Muddana et al. calculated free energies of binding with their minima mining approach; i.e., the weighted averaging over dominant binding modes from classical dynamics, for CB7 host/guest systems in very good agreement with experiment [13]. This study clearly shows the power of SQM-DH-based approaches, although massive upscaling is necessary to address real-life virtual drug design problems.

3.2.2. Protein/Ligand Interactions

Such an upscaling as mentioned above is indeed possible, and is shown in studies from the Korth group for protein/ligand interactions [10,100]. The complete PDBbind2007 benchmark set is subjected to a systematic benchmark of MM, DFT, and SQM methods, and a set of smaller model systems is additionally checked with high-level WFT calculations. It is found that the accuracy first of all depends on the level of theory: SQM-DH methods perform very similarly to DFT-D and substantially better than MM methods. The performance for different methods at the same level of theory turned out to be rather similar. QM calculations can be restricted to smaller model systems of the full complexes without losing predictive capability, which brings all benchmark entries well into the reach of SQM-DH methods. Automatic model system preparation and the treatment of entropic effects are listed as major remaining challenges, though a more detailed analysis also revealed options for the further improvement of SQM-DH methods. An exhaustive comparison of WFT, DFT, and SQM methods with high-level reference data for the smaller complexes led the authors to recommend TPSS-D3(abc)/def-TZVPP as DFT method and PM6-DH+ as SQM-DH approach.
A series of docking studies from the Hobza group rely on SQM-DH methods with implicit solvation models and force field based entropic terms [63,101,102,103,104,105]. In these studies, SQM-DH agrees very well with DFT-D results and is found to be consistently better than classical methods, especially also for the identification of decoy structures.
A more systematic approach for the calculation of free energies of binding for protein/ligand systems was put forward by Stigliani et al. [106]. They apply SQM-DH methods on selected binding modes based on Autodock Vina results, and indeed find a good performance of Vina for the posing step, but the correct score only after SQM-DH re-optimization. Another systematic approach was published by Ucisik et al. [107]. Here, a special type of weighted averaging is applied, which the authors call “blurring”, and which arrives at ensembles of binding modes via a local sampling algorithm. Correlation between calculated and experimental data increases substantially with the use of SQM-DH methods. Further studies were published with SQM-DH level QM or QM/MM treatment of protein/ligand interactions [25,108,109,110,111,112,113].

4. Summary and Outlook

Recent years have seen increased activity in the field of improved virtual drug design models based on the QM-rescoring of classically-obtained binding mode ensembles. “Enhanced”, dispersion (D) and hydrogen-bond (H) corrected QM methods like DFT-D and SQM-DH are at the core of this “new wave” of QM-based scoring schemes. Unlike older DFT and SQM approaches, these methods perform very accurately for the non-covalent interactions that dominate protein/ligand binding. DFT-D shows an accuracy close to the most sophisticated (and much more costly) quantum chemistry approaches. SQM-DH methods are almost as accurate as DFT-D, but are additionally extremely fast, though potential hypersurface errors can lead to problems when optimizing highly-symmetric systems. DFT-D is now the de facto standard approach for large, non-covalent interaction-dominated systems, but SQM-DH and probably also other minimal (slower, but more robust) QM models like HF-3c [114] are an excellent choice for systems outside of the reach of DFT calculations. Such systems are, for instance, protein/protein interfaces as well as functional biomacromolecules and hybrid bio/nano-systems, at least when it comes to in silico design and optimization. Less attention needs to be paid to the actual DFT-D or SQM-DH method that is used. Fortunately, DFT-D methods are now readily available in standard quantum chemistry codes and PM6-DH2 as well as PM6-DH+ are available in MOPPAC [115]. Free D and DH codes are available for the implementation of new research ideas [116]. It seems very likely that DFT-D and SQM-DH methods will have a great impact on the development of more accurate approaches for treating protein/ligand interactions.

Acknowledgments

Nusret Duygu Yilmazer and Martin Korth gratefully acknowledge financial support from the Barbara-Mez-Starck foundation.

Author Contributions

Nusret Duygu Yilmazer and Martin Korth did literature work. Martin Korth wrote the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Leach, A.R.; Shoichet, B.K.; Peishoff, C.E. Prediction of protein-ligand interactions. Docking and scoring: Successes and gaps. J. Med. Chem. 2006, 49, 5851–5855. [Google Scholar] [CrossRef] [PubMed]
  2. Klebe, G. Virtual ligand screening: Strategies, perspectives and limitations. Drug Discov. Today 2006, 11, 580–594. [Google Scholar] [CrossRef] [PubMed]
  3. Jorgensen, W.L. Efficient Drug Lead Discovery and Optimization. Acc. Chem. Res. 2009, 42, 724–733. [Google Scholar] [CrossRef] [PubMed]
  4. Gilson, M.K.; Zhou, H.-X. Calculation of protein-ligand binding affinities. Annu. Rev. Biophys. Biomol. Struct. 2007, 36, 21–42. [Google Scholar] [CrossRef] [PubMed]
  5. Gallicchio, E.; Levy, R.M. Recent Theoretical and Computational Advances for Modeling Protein–Ligand Binding Affinities. Adv. Protein Chem. Struct. Biol. 2011, 85, 27–80. [Google Scholar] [PubMed]
  6. Zhou, H.-X.; Gilson, M.K. Theory of Free Energy and Entropy in Noncovalent Binding. Chem. Rev. 2009, 109, 4092–4107. [Google Scholar] [CrossRef] [PubMed]
  7. Korth, M. A quantum chemical view of enthalpy–entropy compensation. Med. Chem. Commun. 2013, 4, 1025–1033. [Google Scholar] [CrossRef]
  8. Mobley, D.L. Let’s get honest about sampling. J. Comput. Aided Mol. Des. 2012, 26, 93–95. [Google Scholar] [CrossRef] [PubMed]
  9. Paton, R.S.; Goodman, J.M. Hydrogen Bonding and π-Stacking: How Reliable are Force Fields? A Critical Evaluation of Force Field Descriptions of Nonbonded Interactions. J. Chem. Inf. Model. 2009, 49, 944–955. [Google Scholar] [CrossRef]
  10. Yilmazer, N.D.; Heitel, P.; Schwabe, T.; Korth, M. Benchmark of electronic structure methods for protein–ligand interactions based on high-level reference data. J. Theor. Comput. Chem. 2015, 14, 1540001. [Google Scholar] [CrossRef]
  11. Liu, L.; Guo, Q.X. Isokinetic relationship, isoequilibrium relationship, and enthalpy–entropy compensation. Chem. Rev. 2001, 101, 673–695. [Google Scholar] [CrossRef] [PubMed]
  12. Grimme, S. Supramolecular Binding Thermodynamics by Dispersion-Corrected Density Functional Theory. Chem. A Eur. J. 2012, 18, 9955–9964. [Google Scholar] [CrossRef] [PubMed]
  13. Muddana, H.S.; Gilson, M.K. Calculation of Host-Guest Binding Affinities Using a Quantum-Mechanical Energy Model. J. Chem. Theory Comput. 2012, 8, 2023–2033. [Google Scholar] [CrossRef] [PubMed]
  14. Merz, K.M. Using Quantum Mechanical Approaches to Study Biological Systems. Acc. Chem. Res. 2014, 47, 2804–2811. [Google Scholar] [CrossRef] [PubMed]
  15. Klimovich, P.V.; Shirts, M.R.; Mobley, D.L. Guidelines for the analysis of free energy calculations. J. Comput. Aided Mol. Des. 2015, 29, 397–411. [Google Scholar] [CrossRef] [PubMed]
  16. Zheng, Z.; Ucisik, M.N.; Merz, K.M. The Movable Type Method Applied to Protein–Ligand Binding. J. Chem. Theory Comput. 2013, 9, 5526–5538. [Google Scholar] [CrossRef] [PubMed]
  17. Klamt, A. The COSMO and COSMO–RS solvation models. WIREs Comput. Mol. Sci. 2011, 1, 699–709. [Google Scholar] [CrossRef]
  18. Moghaddam, S.; Inoue, Y.; Gilson, M.K. Host–Guest Complexes with Protein–Ligand-Like Affinities: Computational Analysis and Design. J. Am. Chem. Soc. 2009, 131, 4012–4021. [Google Scholar] [CrossRef] [PubMed]
  19. Jensen, J.H. Predicting accurate absolute binding energies in aqueous solution: Thermodynamic considerations for electronic structure methods. Phys. Chem. Chem. Phys. 2015, 17, 12441–12451. [Google Scholar] [CrossRef] [PubMed]
  20. Hobza, P. Calculations on Noncovalent Interactions and Databases of Benchmark Interaction Energies. Acc. Chem. Res. 2012, 45, 663–672. [Google Scholar] [CrossRef] [PubMed]
  21. Grimme, S. Density functional theory with London dispersion correction. WIREs Comput. Mol. Sci. 2011, 1, 211–228. [Google Scholar] [CrossRef]
  22. Korth, M. Error estimates for (semi-)empirical dispersion terms and large biomacromolecules. Org. Biomol. Chem. 2013, 11, 6515–6519. [Google Scholar] [CrossRef] [PubMed]
  23. Korth, M. Empirical Hydrogen-Bond Potential Functions—An Old Hat Reconditioned. Chem. Phys. Chem. 2011, 12, 3131–3142. [Google Scholar] [CrossRef] [PubMed]
  24. Yilmazer, N.D.; Korth, M. Enhanced semiempirical QM methods for biomolecular interactions. Comp. Struct. Biotech. J. 2015, 13, 169–175. [Google Scholar] [CrossRef] [PubMed]
  25. Raha, K.; Merz, K.M. Large-Scale Validation of a Quantum Mechanics Based Scoring Function: Predicting the Binding Affinity and the Binding Mode of a Diverse Set of Protein–Ligand Complexes. J. Med. Chem. 2005, 48, 4558–4575. [Google Scholar] [CrossRef] [PubMed]
  26. Risthaus, T.; Grimme, S. Benchmarking of London Dispersion-Accounting Density Functional Theory Methods on Very Large Molecular Complexes. J. Chem. Theory Comput. 2013, 9, 1580–1591. [Google Scholar] [CrossRef] [PubMed]
  27. Antony, J.; Grimme, S.; Liakos, D.G.; Neese, F. Protein–Ligand Interaction Energies with Dispersion Corrected Density Functional Theory and High-Level Wave Function Based Methods. J. Phys. Chem. A 2011, 115, 11210–11220. [Google Scholar] [CrossRef] [PubMed]
  28. Antony, J.; Grimme, S. Fully ab initio Protein–Ligand Interaction Energies with Dispersion Corrected Density Functional Theory. J. Comput. Chem. 2012, 33, 1730–1739. [Google Scholar] [CrossRef] [PubMed]
  29. He, X.; Zhu, T.; Wang, X.; Liu, J.; Zhang, J.Z.H. Fragment Quantum Mechanical Calculation of Proteins and Its Applications. Acc. Chem. Res. 2014, 47, 2748–2757. [Google Scholar] [CrossRef] [PubMed]
  30. Jia, X.; Wang, X.; Liu, J.; Zhang, J.Z.H.; Mei, Y.; He, X. An Improved Fragment-Based Quantum Mechanical Method for Calculation of Electrostatic Solvation Energy of Proteins. J. Chem. Phys. 2013, 139, 214104. [Google Scholar] [CrossRef] [PubMed]
  31. Liu, J.; Wang, X.; Zhang, J.Z.H.; He, X. Calculation of Protein–Ligand Binding Affinities Based on a Fragment Quantum Mechanical Method. RSC Adv. 2015, 5, 107020–107030. [Google Scholar] [CrossRef]
  32. Antony, J.; Sure, R.; Grimme, S. Using dispersion-corrected density functional theory to understand supramolecular binding thermodynamics. Chem. Commun. 2015, 51, 1764–1774. [Google Scholar] [CrossRef] [PubMed]
  33. Fox, S.J.; Pittock, C.; Fox, T.; Tautermann, C.S.; Malcolm, N.; Skylaris, C.-K. Electrostatic embedding in large-scale first principles quantum mechanical calculations on biomolecules. J. Chem. Phys. 2011, 135, 224107. [Google Scholar] [CrossRef] [PubMed]
  34. Cole, D.J.; Rajendra, E.; Roberts-Thomson, M.; Hardwick, B.; McKenzie, G.J.; Payne, M.C.; Venkitaraman, A.R.; Skylaris, C.K. Interrogation of the Protein–Protein Interactions between Human BRCA2 BRC Repeats and RAD51 Reveals Atomistic Determinants of Affinity. PLoS Comput. Biol. 2011, 7, e1002096. [Google Scholar] [CrossRef] [PubMed]
  35. Fox, S.; Wallnoefer, H.G.; Fox, T.; Tautermann, C.S.; Skylaris, C.K. First Principles-Based Calculations of Free Energy of Binding: Application to Ligand Binding in a Self-Assembling Superstructure. J. Chem. Theory Comput. 2011, 7, 1102–1108. [Google Scholar] [CrossRef] [PubMed]
  36. Fox, S.J.; Dziedzic, J.; Fox, T.; Tautermann, C.S.; Skylaris, C.K. Density functional theory calculations on entire proteins for free energies of binding: Application to a model polar binding site. Proteins Struct. Funct. Bioinf. 2014, 82, 3335–3346. [Google Scholar] [CrossRef]
  37. Simova, L.; Rezac, J.; Hobza, P. Convergence of the Interaction Energies in Noncovalent Complexes in the Coupled-Cluster Methods Up to Full Configuration Interaction. J. Chem. Theory Comput. 2013, 9, 3420–3428. [Google Scholar] [CrossRef] [PubMed]
  38. Kolar, M.; Fanfrlik, J.; Hobza, P. Ligand Conformational and Solvation/Desolvation Free Energy in Protein–Ligand Complex Formation. J. Phys. Chem. B 2011, 115, 4718–4724. [Google Scholar]
  39. Kolar, M.; Fanfrlik, J.; Lepsik, M.; Forti, F.; Luque, F.J.; Hobza, P. Assessing the Accuracy and Performance of Implicit Solvent Models for Drug Molecules: Conformational Ensemble Approaches. J. Phys. Chem. B 2013, 117, 5950–5962. [Google Scholar] [CrossRef]
  40. Söderhjelm, P.; Kongsted, J.; Ryde, U. Ligand Affinities Estimated by Quantum Chemical Calculations. J. Chem. Theory Comput. 2010, 6, 1726–1737. [Google Scholar] [CrossRef] [PubMed]
  41. Mikulskis, P.; Cioloboc, D.; Andrejic, M.; Khare, S.; Brorsson, J.; Genheden, S.; Mata, R.A.; Söderhjelm, P.; Ryde, U. Free-energy perturbation and quantum mechanical study of SAMPL4 octa-acid host–guest binding energies. J. Comput. Aided Mol. Des. 2014, 28, 375–400. [Google Scholar] [CrossRef] [PubMed]
  42. Andrejic, M.; Ryde, U.; Mata, R.A.; Söderhjelm, P. Coupled-Cluster Interaction Energies for 200-Atom Hos–Guest Systems. ChemPhysChem 2014, 15, 3270–3281. [Google Scholar] [CrossRef] [PubMed]
  43. Rao, L.; Zhang, I.Y.; Guo, W.; Feng, L.; Meggers, E.; Xu, X. Nonfitting protein-ligand interaction scoring function based on first-principles theoretical chemistry methods: Development and application on kinase inhibitors. J. Comput. Chem. 2013, 34, 1636–1646. [Google Scholar] [CrossRef] [PubMed]
  44. Gleeson, D.; Tehan, B.; Gleeson, M.P.; Limtrakul, J. Evaluating the enthalpic contribution to ligand binding using QM calculations: Effect of methodology on geometries and interaction energies. Org. Biomol. Chem. 2013, 10, 7053–7061. [Google Scholar] [CrossRef] [PubMed]
  45. Kantardjiev, A.A. Quantum.Ligand.Dock: Protein–ligand docking with quantum entanglement refinement on a GPU system. Nucleic Acids Res. 2012, 40, W415–W422. [Google Scholar] [CrossRef] [PubMed]
  46. Bryce, R.A. Physics-based scoring of protein-ligand interactions: Explicit polarizability, quantum mechanics and free energies. Future Med. Chem. 2011, 3, 683–698. [Google Scholar] [CrossRef] [PubMed]
  47. Gleeson, M.P.; Hannongbua, S.; Gleeson, D. QM methods in structure based design: Utility in probing protein–ligand interactions. J. Mol. Graph. Model. 2010, 29, 507–517. [Google Scholar] [CrossRef] [PubMed]
  48. Zhou, T.; Caflisch, A. High-Throughput Virtual Screening Using Quantum Mechanical Probes: Discovery of Selective Kinase Inhibitors. Chem. Med. Chem. 2010, 5, 1007–1014. [Google Scholar] [CrossRef] [PubMed]
  49. Zhang, X.; Gibbs, A.C.; Reynolds, C.H.; Peters, M.B.; Westerhoff, L.M. Quantum Mechanical Pairwise Decomposition Analysis of Protein Kinase B Inhibitors: Validating a New Tool for Guiding Drug Design. J. Chem. Inf. Model. 2010, 50, 651–661. [Google Scholar] [CrossRef] [PubMed]
  50. Li, X.; Fu, Z.; Merz, K.M. QM/MM refinement and analysis of protein bound retinoic acid. J. Comput. Chem. 2012, 33, 301–310. [Google Scholar] [CrossRef] [PubMed]
  51. Fu, Z.; Li, X.; Miao, Y.; Merz, K.M. Conformational Analysis and Parallel QM/MM X-ray Refinement of Protein Bound Anti-Alzheimer Drug Donepezil. J. Chem. Theory Comput. 2013, 9, 1686–1693. [Google Scholar] [CrossRef] [PubMed]
  52. Temelso, B.; Alser, K.A.; Gauthier, A.; Palmer, A.K.; Shields, G.C. Structural Analysis of α-Fetoprotein (AFP)-like Peptides with Anti-Breast-Cancer Properties. J. Phys. Chem. B 2014, 118, 4514–4526. [Google Scholar] [CrossRef] [PubMed]
  53. Li, X.; He, X.; Wang, B.; Merz, K., Jr. Conformational Variability of Benzamidinium-Based Inhibitors. J. Am. Chem. Soc. 2009, 131, 7742–7754. [Google Scholar] [CrossRef] [PubMed]
  54. Christensen, A.S.; Linnet, T.E.; Borg, M.; Boomsma, W.; Lindorff-Larsen, K.; Hamelryck, T.; Jensen, J.H. Protein Structure Validation and Refinement Using Amide Proton Chemical Shifts Derived from Quantum Mechanics. PLoS ONE 2013, 8, e84123. [Google Scholar] [CrossRef] [PubMed]
  55. Zhu, T.; Zhang, J.Z.H.; He, X. Correction of Erroneously Packed Protein’s Side Chains in the NMR Structure Based on Ab Initio Chemical Shift Calculations. Phys. Chem. Chem. Phys. 2014, 16, 18163–18169. [Google Scholar] [CrossRef] [PubMed]
  56. Sparta, M.; Neese, F. Chemical applications carried out by local pair natural orbital based coupled-cluster methods. Chem. Soc. Rev. 2014, 43, 5032–5041. [Google Scholar] [CrossRef] [PubMed]
  57. Neese, F.; Hansen, A.; Liakos, D.G. Efficient and accurate approximations to the local coupled cluster singles doubles method using a truncated pair natural orbital basis. J. Chem. Phys. 2009, 131, 064103. [Google Scholar]
  58. Neese, F.; Wennmohs, F.; Hansen, A. Efficient and accurate local approximations to coupled-electron pair approaches: An attempt to revive the pair natural orbital method. J. Chem. Phys. 2009, 130, 114108. [Google Scholar] [CrossRef] [PubMed]
  59. Liakos, D.G.; Sparta, M.; Kesharwani, M.K.; Martin, J.M.L.; Neese, F. Exploring the Accuracy Limits of Local Pair Natural Orbital Coupled-Cluster Theory. J. Chem. Theory Comput. 2015, 11, 1525–1539. [Google Scholar] [CrossRef] [PubMed]
  60. Liakos, D.G.; Neese, F. Domain Based Pair Natural Orbital Coupled Cluster Studies on Linear and Folded Alkane Chains. J. Chem. Theory Comput. 2015, 11, 2137–2143. [Google Scholar] [CrossRef] [PubMed]
  61. Liakos, D.G.; Neese, F. Is It Possible to Obtain Coupled Cluster Quality Energies at near Density Functional Theory Cost? Domain-Based Local Pair Natural Orbital Coupled Cluster vs Modern Density Functional Theory. J. Chem. Theory Comput. 2015, 11, 4054–4063. [Google Scholar] [CrossRef] [PubMed]
  62. Foster, M.E.; Sohlberg, K. Computational Investigation of the Role of Counterions and Reorganization Energy in a Switchable Bistable [2]Rotaxane. J. Phys. Chem. A 2011, 115, 7773–7777. [Google Scholar] [CrossRef] [PubMed]
  63. Fanfrlik, J.; Brahmkshatriya, P.S.; Rezac, J.; Jilkova, A.; Horn, M.; Mares, M.; Hobza, P.; Lepsik, M. Quantum Mechanics-Based Scoring Rationalizes the Irreversible Inactivation of Parasitic Schistosoma mansoni Cysteine Peptidase by Vinyl Sulfone Inhibitors. J. Phys. Chem. B 2013, 117, 14973–14982. [Google Scholar] [PubMed]
  64. Gordeev, E.G.; Polynski, M.V.; Ananikov, V.P. Fast and accurate computational modeling of adsorption on graphene: A dispersion interaction challenge. Phys. Chem. Chem. Phys. 2013, 15, 18815–18821. [Google Scholar] [CrossRef] [PubMed]
  65. Margiotta, N.; Marzano, C.; Gandin, V.; Osella, D.; Ravera, M.; Gabano, E.; Platts, J.A.; Petruzzella, E.; Hoeschele, J.D.; Natile, G. Revisiting [PtCl2(cis-1,4-DACH)]: An Underestimated Antitumor Drug with Potential Application to the Treatment of Oxaliplatin-Refractory Colorectal Cancer. J. Med. Chem. 2012, 55, 7182–7192. [Google Scholar] [CrossRef] [PubMed]
  66. Molla, M.R.; Gehrig, D.; Roy, L.; Kamm, V.; Paul, A.; Laquai, F.; Ghosh, S. Self-Assembly of Carboxylic Acid Appended Naphthalene Diimide Derivatives with Tunable Luminescent Color and Electrical Conductivity. Chem. Eur. J. 2014, 20, 760–771. [Google Scholar] [CrossRef] [PubMed]
  67. Raju, R.K.; Hillier, I.H.; Burton, N.A.; Vincent, M.A.; Doudou, S.; Bryce, R.A. The effects of perfluorination on carbohydrate-pi interactions: Computational studies of the interaction of benzene and hexafluorobenzene with fucose and cyclodextrin. Phys. Chem. Chem. Phys. 2010, 12, 7959–7967. [Google Scholar] [CrossRef] [PubMed]
  68. Ramraj, A.; Hillier, I.H.; Vincent, M.A.; Burton, N.A. Assessment of approximate quantum chemical methods for calculating the interaction energy of nucleic acid bases with graphene and carbon nanotubes. Chem. Phys. Lett. 2010, 484, 295–298. [Google Scholar] [CrossRef]
  69. Ramraj, A.; Raju, R.K.; Wang, Q.; Hillier, I.H.; Bryce, R.A.; Vincent, M.A. An evaluation of the GLYCAM06 and MM3 force fields, and the PM3-D* molecular orbital method for modelling prototype carbohydrate–aromatic interactions. J. Mol. Graph. Model. 2010, 29, 321–325. [Google Scholar] [CrossRef] [PubMed]
  70. Sharma, R.; McNamara, J.P.; Raju, R.K.; Vincent, M.A.; Hillier, I.H.; Morgado, C.A. The interaction of carbohydrates and amino acids with aromatic systems studied by density functional and semi-empirical molecular orbital calculations with dispersion corrections. Phys. Chem. Chem. Phys. 2008, 10, 2767–2774. [Google Scholar] [CrossRef] [PubMed]
  71. Vincent, M.A.; Hillier, I.H. Accurate Prediction of Adsorption Energies on Graphene, Using a Dispersion-Corrected Semiempirical Method Including Solvation. J. Chem. Inf. Model. 2014, 54, 2255–2260. [Google Scholar] [CrossRef] [PubMed]
  72. Wang, Q.; Suzuki, K.; Nagashima, U.; Tachikawa, M.; Yan, S. Geometric isotope effects on small chloride ion water clusters with path integral molecular dynamics simulations. Chem. Phys. 2013, 426, 38–47. [Google Scholar] [CrossRef]
  73. Wang, Q.; Suzuki, K.; Nagashima, U.; Tachikawa, M.; Yan, S. Path integral molecular dynamic study of nuclear quantum effect on small chloride water clusters of Cl−(H2O)1–4. Chem. Phys. 2013, 419, 229–236. [Google Scholar] [CrossRef]
  74. Pietrusiewicz, K.M.; Flis, A.; Ujj, V.; Koertvelyesi, T.; Drahos, L.; Pongracz, P.; Kollar, L.; Keglevich, G. Synthesis and Use of Borane and Platinum(II) Complexes of 3-Diphenylphosphino-1-phenylphospholane (LuPhos). Heteroatom Chem. 2011, 22, 730–736. [Google Scholar] [CrossRef]
  75. Lupan, A.; Kun, A.-Z.; Carrascoza, F.; Silaghi-Dumitrescu, R. Performance comparison of computational methods for modeling α-helical structures. J. Mol. Model. 2013, 19, 193–203. [Google Scholar] [CrossRef] [PubMed]
  76. Sedlak, R.; Janowski, T.; Pitonak, M.; Rezac, J.; Pulay, P.; Hobza, P. Accuracy of Quantum Chemical Methods for Large Noncovalent Complexes. J. Chem. Theory Comput. 2013, 9, 3364–3374. [Google Scholar] [CrossRef] [PubMed]
  77. Prenosil, O.; Pitonak, M.; Sedlak, R.; Kabelac, M.; Hobza, P. H-Bonding Cooperativity Effects in Amyloids: Quantum Mechanical and Molecular Mechanics Study. Z. Phys. Chem. 2011, 225, 553–574. [Google Scholar] [CrossRef]
  78. Rezac, J.; Riley, K.E.; Hobza, P. Benchmark Calculations of Noncovalent Interactions of Halogenated Molecules. J. Chem. Theory Comput. 2012, 8, 4285–4292. [Google Scholar] [CrossRef] [PubMed]
  79. Li, A.; Muddana, H.S.; Gilson, M.K. Quantum Mechanical Calculation of Noncovalent Interactions: A Large-Scale Evaluation of PMx, DFT, and SAPT Approaches. J. Chem. Theory Comput. 2014, 10, 1563–1575. [Google Scholar] [CrossRef] [PubMed]
  80. Foster, M.E.; Sohlberg, K. A New Empirical Correction to the AM1 Method for Macromolecular Complexes. J. Chem. Theory Comput. 2010, 6, 2153–2166. [Google Scholar] [CrossRef] [PubMed]
  81. Foster, M.E.; Sohlberg, K. Self-consistent addition of an atomic charge dependent hydrogen-bonding correction function. Comp. Theor. Chem. 2012, 984, 9–12. [Google Scholar] [CrossRef]
  82. Korth, M.; Pitonak, M.; Rezac, J.; Hobza, P. A Transferable H-Bonding Correction for Semiempirical Quantum-Chemical Methods. J. Chem. Theory Comput. 2010, 6, 344–352. [Google Scholar] [CrossRef] [PubMed]
  83. Korth, M. Third-Generation Hydrogen-Bonding Corrections for Semiempirical QM Methods and Force Fields. J. Chem. Theory Comput. 2010, 6, 3808–3816. [Google Scholar] [CrossRef]
  84. Rezac, J.; Hobza, P. Advanced Corrections of Hydrogen Bonding and Dispersion for Semiempirical Quantum Mechanical Methods. J. Chem. Theory Comput. 2012, 8, 141–151. [Google Scholar] [CrossRef] [PubMed]
  85. Kromann, J.C.; Christensen, A.S.; Steinmann, C.; Korth, M.; Jensen, J.H. A third-generation dispersion and third-generation hydrogen bonding corrected PM6 method: PM6-D3H+. PeerJ 2014, 2, e449. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  86. Avila-Salas, F.; Sandoval, C.; Caballero, J.; Guinez-Molinos, S.; Santos, L.S.; Cachau, R.E.; Gonzalez-Nilo, F.D. Study of Interaction Energies between the PAMAM Dendrimer and Nonsteroidal Anti-Inflammatory Drug Using a Distributed Computational Strategy and Experimental Analysis by ESI-MS/MS. J. Phys. Chem. B 2012, 116, 2031–2039. [Google Scholar] [CrossRef] [PubMed]
  87. Borbulevych, O.Y.; Plumley, J.A.; Martin, R.I.; Merz, Jr.K.M.; Westerhoff, L.M. Accurate macromolecular crystallographic refinement: Incorporation of the linear scaling, semiempirical quantum-mechanics program DivCon into the PHENIX refinement package. Acta Crystallogr. Sect. D Biol. Crystallogr. 2014, 70, 1233–1247. [Google Scholar] [CrossRef] [PubMed]
  88. Muddana, H.S.; Sapra, N.V.; Fenley, A.T.; Gilson, M.K. The electrostatic response of water to neutral polar solutes: Implications for continuum solvent modeling. J. Chem. Phys. 2013, 138, 224504. [Google Scholar] [CrossRef] [PubMed]
  89. CarvalhoMaia, J.D.; Urquiza Carvalho, G.A.; Mangueira, C.P., Jr.; Santana, S.R.; Formiga Cabral, L.A.; Rocha, G.B. GPU Linear Algebra Libraries and GPGPU Programming for Accelerating MOPAC Semiempirical Quantum Chemistry Calculations. J. Chem. Theory Comput. 2012, 8, 3072–3081. [Google Scholar]
  90. Brandon, C.J.; Martin, B.P.; McGee, K.J.; Stewart, J.J.P.; Braun-Sand, S.B. An approach to creating a more realistic working model from a protein data bank entry. J. Mol. Model. 2015, 21, 1–11. [Google Scholar] [CrossRef] [PubMed]
  91. Conti, S.; Ceccini, M. Accurate and Efficient Calculation of the Desorption Energy of Small Molecules from Graphene. J. Phys. Chem. C 2015, 119, 1867–1879. [Google Scholar] [CrossRef]
  92. Korth, M. Large-scale virtual high-throughput screening for the identification of new battery electrolyte solvents: Evaluation of electronic structure theory methods. Phys. Chem. Chem. Phys. 2014, 16, 7919–7926. [Google Scholar] [CrossRef] [PubMed]
  93. Husch, T.; Yilmazer, N.D.; Balducci, A.; Korth, M. Large-scale virtual high-throughput screening for the identification of new battery electrolyte solvents: Computing infrastructure and collective properties. Phys. Chem. Chem. Phys. 2015, 17, 3394–3401. [Google Scholar] [CrossRef] [PubMed]
  94. Husch, T.; Korth, M. How to estimate solid-electrolyte-interphase features when screening electrolyte materials. Phys. Chem. Chem. Phys. 2015, 17, 22799–22808. [Google Scholar] [CrossRef] [PubMed]
  95. Husch, T.; Korth, M. Charting the known chemical space for non-aqueous Lithium-air battery electrolyte solvents. Phys. Chem. Chem. Phys. 2015, 17, 22596–22603. [Google Scholar] [CrossRef] [PubMed]
  96. Kubota, D.; Lopes Macedo, O.F.; Santana Andrade, G.R.; Conegero, L.S.; Almeida, L.E.; Costa, N.B., Jr.; Gimenez, I.F. Structural and theoretical-experimental physicochemical study of trimethoprim/randomly methylated-β-cyclodextrin binary system. Carbohydr. Res. 2011, 346, 2746–2751. [Google Scholar] [CrossRef] [PubMed]
  97. Leclercq, L.; Suisse, I.; Roussel, P.; Agbossou-Niedercorn, F. Inclusion of tetrabutylammonium cations in a chiral thiazolium/triflate network: Solid state and solution structural investigation. J. Mol. Struct. 2012, 1010, 152–157. [Google Scholar] [CrossRef]
  98. Xia, Y.; Wang, X.; Zhang, Y.; Luo, B. Theoretical Study on Interactions of β-cyclodextrin with Trans-dichloro(dipyridine) platinum(II). Comp. Theor. Chem. 2011, 967, 213–218. [Google Scholar] [CrossRef]
  99. Gavvala, K.; Sengupta, A.; Hazra, P. Modulation of Photophysics and pKa Shift of the Anti-cancer Drug Camptothecin in the Nanocavities of Supramolecular Hosts. ChemPhysChem 2013, 14, 532–542. [Google Scholar] [CrossRef] [PubMed]
  100. Yilmazer, N.D.; Korth, M. Comparison of Molecular Mechanics, Semi-Empirical Quantum Mechanical, and Density Functional Theory Methods for Scoring Protein-Ligand Interactions. J. Phys. Chem. B 2013, 117, 8075–8084. [Google Scholar] [CrossRef] [PubMed]
  101. Fanfrlik, J.; Bronowska, A.K.; Rezac, J.; Prenosil, O.; Konvalinka, J.; Hobza, P. A Reliable Docking/Scoring Scheme Based on the Semiempirical Quantum Mechanical PM6-DH2 Method Accurately Covering Dispersion and H-Bonding: HIV-1 Protease with 22 Ligands. J. Phys. Chem. B 2010, 114, 12666–12678. [Google Scholar] [CrossRef] [PubMed]
  102. Dobes, P.; Fanfrlik, J.; Rezac, J.; Otyepka, M.; Hobza, P. Transferable scoring function based on semiempirical quantum mechanical PM6-DH2 method: CDK2 with 15 structurally diverse inhibitors. J. Comput. Aided Mol. Des. 2011, 25, 223–235. [Google Scholar] [CrossRef] [PubMed]
  103. Jilkova, A.; Rezacova, P.; Lepsik, M.; Horn, M.; Vachova, J.; Fanfrlik, J.; Brynda, J.; McKerrow, J.H.; Caffrey, C.R.; Mares, M. Structural Basis for Inhibition of Cathepsin B Drug Target from the Human Blood Fluke, Schistosoma mansoni. J. Biol. Chem. 2011, 286, 35770–35781. [Google Scholar] [CrossRef] [PubMed]
  104. Vorlova, B.; Nachtigallova, D.; Jiraskova-Vanickova, J.; Ajani, H.; Jansa, P.; Rezac, J.; Fanfrlik, J.; Otyepka, M.; Hobza, P.; Konvalinka, J.; et al. Malonate-based inhibitors of mammalian serine racemase: Kinetic characterization and structure-based computational study. Eur. J. Med. Chem. 2015, 89, 189–197. [Google Scholar] [CrossRef] [PubMed]
  105. Sulimov, V.B.; Katkova, E.V.; Oferkin, I.V.; Sulimov, A.V.; Romanov, A.N.; Roschin, A.I.; Beloglazova, I.B.; Plekhanova, O.S.; Tkachuk, V.A.; Sadovnichiy, V.A. Application of Molecular Modeling to Urokinase Inhibitors Development. BioMed Res. Int. 2014, 2014, 625176. [Google Scholar] [CrossRef] [PubMed]
  106. Stigliani, J.-L.; Bernardes-Genisson, V.; Bernadou, J.; Pratviel, G. Cross-docking study on InhA inhibitors: A combination of Autodock Vina and PM6-DH2 simulations to retrieve bio-active conformations. Org. Biomol. Chem. 2012, 10, 6341–6349. [Google Scholar] [CrossRef] [PubMed]
  107. Ucisik, M.N.; Zheng, Z.; Faver, J.C.; Merz, K.M. Bringing Clarity to the Prediction of Protein–Ligand Binding Free Energies via “Blurring”. J. Chem. Theory Comput. 2014, 10, 1314–1325. [Google Scholar] [CrossRef] [PubMed]
  108. Kamel, K.; Kolinski, A. Computational study of binding of epothilone A to β-tubulin. Acta Biochim. Pol. 2011, 58, 255–260. [Google Scholar] [PubMed]
  109. Kamel, K.; Kolinski, A. Assessment of the free binding energy of 1,25-dihydroxyvitamin D-3 and its analogs with the human VDR receptor model. Acta Biochim. Pol. 2012, 59, 653–660. [Google Scholar] [PubMed]
  110. Ahmed, M.; Sadek, M.M.; Abouzid, K.A.; Wang, F. In silico design: Extended molecular dynamic simulations of a new series of dually acting inhibitors against EGFR and HER2. J. Mol. Graph. Model. 2013, 44, 220–231. [Google Scholar] [CrossRef] [PubMed]
  111. Fong, P.; McNamara, J.P.; Hillier, I.H.; Bryce, R.A. Assessment of QM/MM Scoring Functions for Molecular Docking to HIV-1 Protease. J. Chem. Inf. Model. 2009, 49, 913–924. [Google Scholar] [CrossRef] [PubMed]
  112. Pan, X.-L.; Liu, W.; Liu, J.-Y. Mechanism of the Glycosylation Step Catalyzed by Human α-Galactosidase: A QM/MM Metadynamics Study. J. Phys. Chem. B 2013, 117, 484–489. [Google Scholar] [CrossRef] [PubMed]
  113. Fanfrlik, J.; Kolar, M.; Kamlar, M.; Hurny, D.; Ruiz, F.X.; Cousido-Siah, A.; Mitschler, A.; Rezac, J.; Munusamy, E.; Lepsik, M.; et al. Modulation of Aldose Reductase Inhibition by Halogen Bond Tuning. ACS Chem. Biol. 2013, 8, 2484–2492. [Google Scholar] [CrossRef] [PubMed]
  114. Sure, R.; Grimme, S. Corrected small basis set Hartree-Fock method for large systems. J. Comput. Chem. 2013, 34, 1672–1685. [Google Scholar] [CrossRef] [PubMed]
  115. OPENMOPAC. Available online: http://www.openmopac.net (accessed on 13 January 2015).
  116. Hydrogen-bond-correction-f3. Available online: https://github.com/jensengroup/hydrogen-bond-correction-f3 (accessed on 13 January 2015).

Share and Cite

MDPI and ACS Style

Yilmazer, N.D.; Korth, M. Recent Progress in Treating Protein–Ligand Interactions with Quantum-Mechanical Methods. Int. J. Mol. Sci. 2016, 17, 742. https://doi.org/10.3390/ijms17050742

AMA Style

Yilmazer ND, Korth M. Recent Progress in Treating Protein–Ligand Interactions with Quantum-Mechanical Methods. International Journal of Molecular Sciences. 2016; 17(5):742. https://doi.org/10.3390/ijms17050742

Chicago/Turabian Style

Yilmazer, Nusret Duygu, and Martin Korth. 2016. "Recent Progress in Treating Protein–Ligand Interactions with Quantum-Mechanical Methods" International Journal of Molecular Sciences 17, no. 5: 742. https://doi.org/10.3390/ijms17050742

APA Style

Yilmazer, N. D., & Korth, M. (2016). Recent Progress in Treating Protein–Ligand Interactions with Quantum-Mechanical Methods. International Journal of Molecular Sciences, 17(5), 742. https://doi.org/10.3390/ijms17050742

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop