Quantum Mechanical Modeling of the Interactions between Noble Metal (Ag and Au) Nanoclusters and Water with the Effective Fragment Potential Method

Comparative assessment of QM-based and MM-based models for prediction of protein-ligand binding affinity trends

Methods which accurately predict protein-ligand binding strengths are critical for drug discovery. In the last two decades, advances in chemical modelling have enabled steadily accelerating progress in the discovery and optimization of structure-based drug design. Most computational methods currently used in this context are based on molecular mechanics force fields that often have deficiencies in describing the quantum mechanical (QM) aspects of molecular binding. In this study, we show the competitiveness of our QM-based Molecules-in-Molecules (MIM) fragmentation method for characterizing binding energy trends for seven different datasets of protein-ligand complexes. By using molecular fragmentation, the MIM method allows for accelerated QM calculations. We demonstrate that for classes of structurally similar ligands bound to a common receptor, MIM provides excellent correlation to experiment, surpassing the more popular Molecular Mechanics Poisson-Boltzmann Surface Area (MM/PBSA) and Molecular Mechanics Generalized Born Surface Area (MM/GBSA) methods. The MIM method offers a relatively simple, well-defined protocol by which binding trends can be ascertained at the QM level and is suggested as a promising option for lead optimization in structure-based drug design.

Effective Fragment Potentials for Flexible Molecules: Transferability of Parameters and Amino Acid Database

An accurate but efficient description of noncovalent interactions is a key to predictive modeling of biological and materials systems. The effective fragment potential (EFP) is an ab initio-based force field that provides a physically meaningful decomposition of noncovalent interactions of a molecular system into Coulomb, polarization, dispersion, and exchange-repulsion components. An EFP simulation protocol consists of two steps, preparing parameters for molecular fragments by a series of ab initio calculations on each individual fragment, and calculation of interaction energy and properties of a total molecular system based on the prepared parameters. As the fragment parameters (distributed multipoles, polarizabilities, localized wave function, etc.) depend on a fragment geometry, straightforward application of the EFP method requires recomputing parameters of each fragment if its geometry changes, for example, during thermal fluctuations of a molecular system. Thus, recomputing fragment parameters can easily become both computational and human bottlenecks and lead to a loss of efficiency of a simulation protocol. An alternative approach, in which fragment parameters are adjusted to different fragment geometries, referred to as "flexible EFP", is explored here. The parameter adjustment is based on translations and rotations of local coordinate frames associated with fragment atoms. The protocol is validated on extensive benchmark of amino acid dimers extracted from molecular dynamics snapshots of a cryptochrome protein. A parameter database for standard amino acids is developed to automate flexible EFP simulations in proteins. To demonstrate applicability of flexible EFP in large-scale protein simulations, binding energies and vertical electron ionization and electron attachment energies of a lumiflavin chromophore of the cryptochrome protein are computed. The results obtained with flexible EFP are in a close agreement with the standard EFP procedure but provide a significant reduction in computational cost.