GPU Accelerated Quantum Virtual Screening: Application for the Natural Inhibitors of New Dehli Metalloprotein (NDM-1)

Abstract

Quantum mechanical approaches for the massive computation on large biological system such as virtual screening in drug design and development have presented a challenge to computational chemists for many years. In this study, we demonstrated that by taking advantage of rapid growth of GPU-based hardware and software (i.e., teraChem), it is feasible to perform virtual screening of a refined chemical library at quantum mechanical level in order to identify the lead compounds with improved accuracy, especially for the drug targets such as metalloproteins in which significant charge transfer and polarization occur amongst the metal ions and their coordinated amino acids. Our calculations predicted four nature compounds (i.e., Curcumin, Catechin, menthol, and Ferulic acid) as the suitable inhibitors for antibiotics resistance against New Delhi Metallo-β-lactamase-1 (NDM-1). Molecular orbitals (MOs) of the QM region of metal ions and their coordinated residues indicate that the bridged hydroxide ion delocalized the electron over the Zn-OH-Zn group at HOMO, different from MOs when the OH⁻ is not presented in NDM-1. This indicates that the bridged hydroxide ion plays an important role in the design of antibiotics and other inhibitors targeting the metalloproteins.

Introduction

The drug discovery and development is an expensive and complex enterprise. The role of the computational chemist within this team has become firmly established as the computer-aided drug design impacts the development across the discovery time line, from target identification and hit discovery to hit-to-lead optimization and safety profiling (Loughney et al., 2011).

A common activity of the computational drug design is the prediction of the bound conformation of a ligand (“Docking”) or millions of ligands to a given receptor binding pocket (“virtual screening”). In general, widely used knowledge-based, empirical or force-field based scoring functions are capable of providing docked poses with reasonable accuracy. Along this line, development of virtual screening has been made in the areas of ligand-based methods including the field based methods (Jilek and Cramer, 2004), substructural descriptors (Vogt et al., 2010) etc., and of structure-based methods with many docking softwares available (see Sukumar and Das, 2011 and its references therein). More recently, the artificial intelligence (AI) and machine learning technologies have also been developed for virtual screening in drug discovery (Goh et al., 2017). However, these approaches are usually not applicable for the complicated pockets such as those with metal ions or highly polarizable sites (Wang et al., 2011). In these cases, the incorporation of the quantum mechanical calculations into the docking or virtual screening process (Illingworth et al., 2008) becomes essential for obtaining the reasonably accurate prediction which serves as the stepping-stone for further hit-to-lead optimization.

Over the last two decades, there has been a dramatic change in the computational power and accessibility of quantum mechanical software programs, enabling the application of quantum mechanical (QM) scheme to eve larger and more complex systems. It remains however that the highest level of QM can only be applied to modest system sizes of 10 s of atoms. A number of general strategies have been developed and employed to further extend the size of systems for drug-protein systems (Mucs and Bryce, 2013). The first is the QM/MM approach that describes a subregion of interest in electronic detail via QM and couples it to a larger environment modeled at the MM level. This approach has been widely used in study of enzyme activity (see review) (Senn and Thiel, 2009). The second approach is in the application of Fragment Molecular Orbital (FMO) (Fedorov et al., 2012). FMO scheme fragments the large molecules into small moieties as isolated monomers by using hybrid projection operators or adaptive frozen orbital. Non-additive terms can be recovered by considering combinations of dimers and trimers of the fragments. The third approach is the linear scaling methods (Zhou et al., 2010) which seek to reduce the scaling of computational cost order of N basis function. Several approaches have been developed by using localized molecular orbitals and divide-and-conquer method, as well as linear scaling wave plan technique. They have successfully applied to protein-ligand and protein-protein interactions at the QM level (Heady et al., 2006; Cole et al., 2010).

On the other hand, the rapid advances of computer power and technology brought up the importance of the utilization of computational architectures to accelerate the QM calculations. Distribution of the computational overheads of QM algorithms on parallel CPU-based architectures leads to tractable ab initio and density function QM calculations on thousands of basis functions (Valiev et al., 2010). Recently, work has focused on the adaption of QM algorithms to graphics processing units (GPUs) architectures. GPU-based approaches for DFT (Ufimtsev and Martinez, 2008, 2009a,b; Yasuda, 2008) and MP2 (Vogt et al., 2008) calculations have been developed and the software has also been advanced (Genovese et al., 2009; Ufimtsev and Martinez, 2009b; Asadchev and Gordon, 2012; Titov et al., 2013; Frisch et al., 2016).

This study emphasizes on the implementation of QM/MM approach into virtual screening for chemical libraries. By taking advantage of GPU-based architectures, we have developed a combined strategy to screen a library at QM level with a promising speed. To our knowledge, this is the first time to investigate the feasibility of GPU-based virtual screening of a chemical library at the QM/MM level.

We use New Delhi Metallo-β-lactamase-1 (NDM-1) as the example. The NDM-1 protein bearing Klebsiella pneumoniae and Escherichia coli infections, which have broad spectrum against all beta-lactams tested, were first discovered in India in 2008 (Yong et al., 2009). The occurrence of NDM-1 exaggerates the problem of bacteria resistance, and its rapid spreading also makes it become one of major threats to human health in recent years. However, no clinical effective inhibitors were reported so far. NDM-1 belongs to the family of β-lactamases, which can further be divided into A, B, C, and D subgroups. Basically, subfamily A, C, and D can be grouped into serine proteases, while subfamily B is a metalloprotein containing either one or two zinc ions to keep its activity. NDM-1 belongs to B1 family, in which two zinc ions have been identified with a hydroxide ion bridged between two zinc cofactors (Thomas et al., 2011). Based on the crystal structure NDM-1 complexed with captopril (King et al., 2012), the reaction mechanism of hydrolysis of ampicillin which leads to the drug resistance of bacteria has been studied both experimentally (Zhang and Hao, 2011) and theoretically (Zheng and Xu, 2013).

Methodological details

Model construction

The structure of NDM-1 was taken from crystal structure of the complex between NDM-1 and hydrolysed ampicilin (PDB 3Q6X, Zhang and Hao, 2011). Previous QM/MM study has suggested a hydroxide ion connecting two zinc cofactors for the hydrolysis of antibiotics and therefore been included in our previous model (Thomas et al., 2011; Zheng and Xu, 2013). Overall, the zinc ions are coordinated with residues (OH⁻, His120, His122, Asp123), and with residues (OH⁻, His189, Cys208, His250), respectively. The position of OH⁻ group was taken from previous study (Zheng and Xu, 2013), while the Cys208 is deprotonated with net charge of −1 as shown in Figure 1.

Figure 1

Molecular mechanical virtual screening

Our docking program was used for virtual screening of the natural product library to the NDM-1 protein. The docking program consists of two-steps: (1) performing a multiple-copy simultaneous search “MCSS” (Zeng and Treutlein, 1999; Zeng, 2000) of common fragments such as benzene, pyridine, pyrimidine, pyrazine and phenol on the binding sites of a protein. The details have been described else (Zeng and Treutlein, 1999). In this study, we use the residue Glu123 as the center of the binding site. (2) docking a compound into the binding site of the protein. The fragments of functional groups were used as the guidance to place the compound inside the binding site and a torsional Monte Carlo minimization was performed to optimize the ligand conformation. Details of docking algorithm have been briefly described before (Stehn et al., 2013; Treutlein et al., 2017). Our software package quCBit (www.medchemsoft.com) has been developed and applied in this purpose. All-atom Charmm22 force field (MacKerell et al., 1998) was used and dielectric constant of 10 was used to screen the electrostatic interactions as implemented.

The docked conformations were selected based on its overlay criteria with the MCSS minima and the molecular mechanical interaction energies between the compound and protein. The overlay criterion is designated as following steps. Firstly, the sub-groups correspond to the fragment types within a compound. Secondly, the center of the sub-group of each docked conformation to the center of each MCSS minimum was calculated. If the distance is less than 2.0 Å, the sub-group is considered as overlay with MCSS minima. Only the docked conformations with all the sub-groups overlay were retained. By trial and error, interaction energy of 10.0 kcal/mol was used as threshold to select these conformations that will be subjected to further QM refinement.

GPU-based QM refinement

A QM/MM scheme (Warshel and Levitt, 1976) was used for the QM refinement of all the docked conformations for each compound. In the scheme, the QM region contains residues two zinc ions, hydroxide ion, 6 residues (i.e., His120, His122, Asp124, His189, Cys208, and His250), as well as the compound. The remaining part is the MM region. Saturated hydrogen was added for the boundary between QM and MM regions (Field et al., 1990). The total binding energies of protein-compound in solutions were calculated as

where and should be calculated using QM/MM for the complex and protein with Polarizable continuum model (PCM) to mimic the solvent. However, current version of teraChem cannot perform this type of calculations so that we use gas-phase calculation for complex and protein. As the protein atoms are fixed in both complex and protein during the calculations, the errors between the results of gas-phase and PCM could be minor. Therefore, the Equation (1) thus becomes

in which the E_complex and E_protein are the total energies of complex between protein and compound, and the total energy of the protein only, respectively. is the energy of each compound in aqueous solution obtained using polarizable continuum model (PCM) (Barone et al., 1997) via Gaussian09 (Frisch et al., 2009) with dielectric constant of 80.0.

The structures of complex and NDM-1 protein were minimized using the combined QM/MM approach in which the QM is calculated using teraChem (Ufimtsev and Martinez, 2009b) and the MM region calculated using Amber12 (Case et al., 2012). Two hundred minimization steps were performed for each docked conformation obtained from section Molecular Mechanical Virtual Screening. During the minimization, DFT method with B3LYP functional (Lee et al., 1988; Becke, 1993) and 6-31G basis set (Ditchfield et al., 1971) was used for the QM region. We choose the standard basis set of 6-31G only as it significantly reduces the computational cost comparing with the 6-31G(d) with the minimum effect on the energy and structures of the QM region. For the MM region, Amber ff14SB force field (Maier et al., 2015) was used for all amino acids, and general atom force field (GAFF) (Wang et al., 2004) was applied for the compounds. All the teraChem/Amber calculations were carried out on GPU computers (GTX780) with two E3-1200 CPUs and all the Guassian09/Amber calculations were performed on two E5-2620 CPUs with 12 threads.

Results and discussions

Outcome of virtual screening

Our docking protocol was firstly used to predict the binding mode of hydrolysed ampicillin to the NDM-1, in comparison with crystal structure (PDB 3Q6X, Zhang and Hao, 2011). At first, our MM docking calculations give 21 conformations based on the score of protein-ligand interaction energies less than 10.00 kcal/mol and overlay with the MCSS minima as described in Method section Molecular Mechanical Virtual Screening. Using the GPU-accelerated QM refinement, these conformations were further optimized in the complex with NDM-1. Figure 2A plots the correlation between binding energies of docked conformations after QM refinement and its heavy atom RMSD to the crystal structure. Overall, the best conformation (shown in Figure 2B) of hydrolysed ampicillin with NDM-1 have the lowest binding energies of −79.82 kcal/mol with RMSD of 1.42 Å from the crystal structure. Based the QM protein-ligand binding energies ΔE, our combined approach of MM docking and GPU-accelerated QM refinement is capable to predict the native structure of ligand binding to the protein.

Figure 2

By applying our combined approach to screen the library of 34 nature compounds (Thakur et al., 2013), only 29 compounds were docked into the binding sites. Table 1 list the best binding energies of these compounds with NDM-1. Using the hydrolysed ampicillin as the reference (E = −79.82 kcal/mol), four compounds are considered as the potential inhibitors because they have value of E lower than the hydrolysed ampicillin. These hits include three compounds (Curcumin-Keto, Catechin, and Menthol) with strong binding energies of −84.25 kcal/mol, −83.76 kcal/mol, and −81.08 kcal/mol, and one compound (Ferulic acid) with equivalent binding energy of −79.12 kcal/mol, respectively.

The complex structures of these inhibitors to NDM-1 are shown in Figure 3. The keto form of Curcumin has the strongest binding energy to NDM1 (−84.25 kcal/mol). The binding is dominated by interactions between the Zn-OH-Zn and the hydroxyl group of Curcumin. Similar feature is also found for Catechin and Menthol (Figures 3B,C). As shown in Figure 3D), Ferulic acid has the acidic group coordinated to Zinc ion directly; in addition, it also form electrostatic interactions to residue Lys393. This reminisces the binding structures of β-lactam antibiotics (e.g., ampicillin) to NDM-1 (Figure 2).

Figure 3

MOs of divalent zinc ions with bridged hydroxide ion

The metalloproteins involved in bacteria infection have been a common target for the antibiotics. As a metalloprotein, NDM-1 can hydrolyse the antibiotics, and therefore, reduces or diminish its therapeutic effects. The hydrolysis reaction involves hydroxide ion between the two divalent Zinc ions (Thomas et al., 2011; Zheng and Xu, 2013). Figure 4 shows the Molecular Orbitals (MOs) of the QM region of the NDM-1 system with/without hydroxide ion. The MOs are obtained from the Gaussian09 calculations. While the LUMO are similar between two systems with the electron clouds concentrated around Zn2 ion, the OH⁻ bridge makes the HOMO delocalized over the Zn-OH-Zn cluster with a through-space π-π interactions between OH⁻ and S(–) of the deprotonated Cys208. Therefore, the inhibitors targeting on the NDM-1 will be of conjugate nature and interact Zn2 ion between deprotonated Cys208 and residue Gly211. This is consistent to the docked conformation predicted for the potential inhibitors of Curcumin-Keto, Catechin, Menthol, and Ferulic Acid as described above (Figure 3). In fact, the bridged hydroxide ions have been found in many metalloproteins, including oxygen transporters hemocyanin (Alzuet et al., 1997) and hemerythrin (Thomann et al., 1993), methane monoxygenase hydroxylse (Rosenzweig et al., 1993), and aminopeptidases (Bzymek et al., 2005) etc., and therefore, play critical role in the design of compounds that alter biological functions of metalloproteins.

Figure 4

GPU vs. CPU performance

Figure 5A shows the difference of QM calculations of a single docked conformation from each nature compound using GPU and CPU (12 threads) via teraChem and Guassian09, respectively. Overall, the GPU could significantly reduce the time cost up to 4-fold. The computational enhancement also depends on the size of molecules. As shown in Figure 5B, the significant reduction (up to 4-folds) of computational costs was achieved when the calculations were performed for the system with a QM region of more than 120 atoms.

Figure 5

GPU hardware also plays an important role in the performance of calculations. As shown in Figure 5C), GTX1080/1080ti speeded up the calculation up to 5-fold from GTX670. With two units of the most recent hardware GTX1080ti, maximal computer time cost for QM computation of a docked conformation can be reduced to less than 30 min, making the QM virtual screening of chemical library feasible.

It is noticeable that the computational time cost using GPU are consistent for each compound (Figure 5A), while the cost using CPU varies significantly according to its size. This is because the difference of the algorithm used for the computation of integral matrix between the teraChem and others. In order to make the full utilization of GPU architecture, tetraChem has optimized the computations of the integrals onto GPU units (Ufimtsev and Martinez, 2009b), and reduced certain of computation into single precision (Luehr et al., 2011).

Conclusions

QM method has been used for analysis of protein-ligand interaction in drug design and discovery for many years. However, its expensive cost of computer power has hindered its application for virtual screening in drug design and development. Recent years, utilization of computational architecture to accelerate QM calculations has attracted significant attention (Furukawa et al., 2012; Ratcliff et al., 2017), and GPU-based approaches to perform QM calculations (Ufimtsev and Martinez, 2008) have been developed. Here, we have developed a combined strategy with docking the protein-ligand at MM level and re-rank the binding interaction at QM level. We choose the New Delhi metalloprotein 1(NDM-1) as the target due to its important metal binding sites in drug resistance and its significant threat to the world health (Johnson and Woodford, 2013).

With the GPU-accelerated QM refinement, our calculations have demonstrated that it is essential to re-rank the docked conformations at the Quantum mechanical level in order to identify the native binding conformation of a ligand close to the crystal structure, as demonstrated from the calculations of the hydrolysed antibiotic ampicillin to NDM-1 (Figure 2).

Our calculations predict potential inhibitors different from the previous docking study in which Isomargololone and Nimbolide were estimated to be the inhibitors of NDM-1 (Thakur et al., 2013). The discrepancy could be due to the factors that they used the crystal structure of NDM-1 (PDB 3Q6X, Zhang and Hao, 2011) as target without the bridged hydroxide ion and employed empirical docking and scoring function (Goodsell et al., 1996; Morris et al., 1998) only to estimate the binding affinities of the compounds. From our QM calculations, MOs of the QM region shows that the HOMO with the bridged OH⁻ between zinc ions has electron cloud delocalized over the Zn-OH-Zn cluster, while the LUMOs are similar with electron cloud concentrated at Zn2, independent of the hydroxide ion. As a result, the inhibitors of antibiotics resistance against the NDM-1 will be of conjugate nature, by binding over the Zn-OH-Zn, such as Curcumin, Catechin, and Feruclic acid (Figure 3).

The benchmark calculations presented here also provide some insight on the feasibility of virtual screening at quantum mechanical level. For a large QM region of more than 120 atoms which are very common for the metalloproteins, the GPU calculations can achieve the speedup of 4-fold comparing to the CPU method (Figure 5A). Previous study has demonstrated that 200-300 atoms would be needed as the minimum size of QM region in order to obtain accurate activation energies of enzyme reactions in proteins (Kulik et al., 2016). Moreover, the performance based on the GPU hardware infer that two units of the most advanced GPU architecture (GTX1080ti) could achieve the QM computation of single conformation within less than 30 min (Figure 5C). Furthermore, the speedup performance can be scaled up ca. twice by using two GPU processes, as demonstrated using two GTX1080ti units (Figure 5C).

Overall, our calculations demonstrated that it is feasible to use GPU-accelerated quantum mechanical method for virtual screening of chemical library. With the rapid advance of GPU hardware, quantum mechanical calculation in virtual screening will become more sophisticated and accurate at higher level. Here, we also showed that keto form of Curcumin, Catechin, Menthol, and Ferulic Acid may be considered as inhibitors for the antibiotics resistance against NDM-1. As natural product library is significantly smaller library than the synthetic chemicals, the approach proposed here will be of significance for accurate prediction of natural product hits as the lead compounds for the drug development.

References

• 1

  AlzuetG.BubaccoL.CasellaL.RoccoG. P.SalvatoB.BeltraminiM. (1997). The binding of azide to copper-containing and cobalt-containing forms of hemocyanin from the Mediterranean crab Carcinus aestuarii. Eur. J. Biochem.247, 688–694. 10.1111/j.1432-1033.1997.00688.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 2

  AsadchevA.GordonM. S. (2012). New multithreaded hybrid CPU/GPU approach to Hartree-Fock. J. Chem. Theory Comput.8, 4166–4176. 10.1021/ct300526w

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 3

  BaroneV.CossiM.TomasiJ. (1997). A new definition of cavities for the computation of solvation free energies by the polarizable continuum model. J. Chem. Phys.107, 3210–3221. 10.1063/1.474671

  • CrossRef
  • Google Scholar

• 4

  BeckeA. D. (1993). Density-functional thermochemistry. 3. The role of exact exchange. J. Chem. Phys.98, 5648–5652. 10.1063/1.464913

  • CrossRef
  • Google Scholar

• 5

  BzymekK. P.SwierczekS. I.BennettB.HolzR. C. (2005). Spectroscopic and thermodynamic characterization of the E151D and E151A altered leucine aminopeptidases from Aeromonas proteolytica. Inorg. Chem.44, 8574–8580. 10.1021/ic051034g

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 6

  CaseD. A.DardenT.CheathamT.III.SimmerlingC.WangJ.DukeR. E.et al. (2012). AMBER 12. San Francisco, CA: University of California.

  • Google Scholar

• 7

  ColeD. J.SkylarisC. K.RajendraE.VenkitaramanA. R.PayneM. C. (2010). Protein-protein interactions from linear-scaling first-principles quantum-mechanical calculations. Epl91:6. 10.1209/0295-5075/91/37004

  • CrossRef
  • Google Scholar

• 8

  DenningtonR.KeithT. A.MillamJ. M. (2016). GaussView, Version 6, Shawnee, KS: Semichem Inc.

  • Google Scholar

• 9

  DitchfieldR.HehreW. J.PopleJ. A. (1971). Self-consistent molecular orbital methods. 9. Extended Gaussian-type basis for molecular-orbital studies of organic molecules. J. Chem. Phys.54:724. 10.1063/1.1674902

  • CrossRef
  • Google Scholar

• 10

  FedorovD. G.NagataT.KitauraK. (2012). Exploring chemistry with the fragment molecular orbital method. Phys. Chem. Chem. Phys.14, 7562–7577. 10.1039/c2cp23784a

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 11

  FieldM. J.BashP. A.KarplusM. (1990). A combined quantum mechanical and molecular mechanical potential for molecular dynamics simulations. J. Comput. Chem.11, 700–733. 10.1002/jcc.540110605

  • CrossRef
  • Google Scholar

• 12

  FrischM. J.TrucksG. W.SchlegelH. B.ScuseriaG. E.RobbM. A.CheesemanJ. R.et al. (2009). Gaussian 09 Rev. A01. Wallingford, CT.

  • Google Scholar

• 13

  FrischM. J.TrucksG. W.SchlegelH. B.ScuseriaG. E.RobbM. A.CheesemanJ. R.et al. (2016). Gaussian 16 Rev. B.01. Wallingford, CT.

  • Google Scholar

• 14

  FurukawaY.KogaR.YasudaK. (2012). Acceleration of computational quantum chemistry by AVX plus CUDA on heterogeneous computer architectures. Abstr. Papers Am. Chem. Soc.243:1.

  • Google Scholar

• 15

  GenoveseL.OspiciM.DeutschT.MehautJ. F.NeelovA.GoedeckerS. (2009). Density functional theory calculation on many-cores hybrid central processing unit-graphic processing unit architectures. J. Chem. Phys.131:8. 10.1063/1.3166140

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 16

  GohG. B.HodasN. O.VishnuA. (2017). Deep learning for computational chemistry. J. Comput. Chem.38, 1291–1307. 10.1002/jcc.24764

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 17

  GoodsellD. S.MorrisG. M.OlsonA. J. (1996). Automated docking of flexible ligands: applications of AutoDock. J. Mol. Recogn.9, 1–5. 10.1002/(sici)1099-1352(199601)9:1<1::Aid-jmr241>3.0.Co;2-6

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 18

  HeadyL.Fernandez-SerraM.ManceraR. L.JoyceS.VenkitaramanA. R.ArtachoE.et al. (2006). Novel structural features of CDK inhibition revealed by an ab initio computational method combined with dynamic simulations. J. Med. Chem.49, 5141–5153. 10.1021/jm060190+

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 19

  IllingworthC. J. R.MorrisG. M.ParkesK. E. B.SnellC. R.ReynoldsC. A. (2008). Assessing the role of polarization in docking. J. Phys. Chem. A112, 12157–12163. 10.1021/jp710169m

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 20

  JilekR. J.CramerR. D. (2004). Topomers: a validated protocol for their self-consistent generation. J. Chem. Inf. Comput. Sci.44, 1221–1227. 10.1021/ci049961d

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 21

  JohnsonA. P.WoodfordN. (2013). Global spread of antibiotic resistance: the example of New Delhi metallo-beta-lactamase (NDM)-mediated carbapenem resistance. J. Med. Microbiol.62, 499–513. 10.1099/jmm.0.052555-0

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 22

  KingD. T.WorrallL. J.GruningerR.StrynadkaN. C. J. (2012). New Delhi Metallo-beta-lactamase: structural insights into beta-lactam recognition and inhibition. J. Am. Chem. Soc.134, 11362–11365. 10.1021/ja303579d

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 23

  KulikH. J.ZhangJ. Y.KlinmanJ. P.MartinezT. J. (2016). How large should the QM region be in QM/MM calculations? The case of catechol O-methyltransferase. J. Phys. Chem. B120, 11381–11394. 10.1021/acs.jpcb.6b07814

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 24

  LeeC. T.YangW. T.ParrR. G. (1988). Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. Phys. Rev. B37, 785–789. 10.1103/PhysRevB.37.785

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 25

  LoughneyD.ClausB. L.JohnsonS. R. (2011). To measure is to know: an approach to CADD performance metrics. Drug Discov. Today16, 548–554. 10.1016/j.drudis.2011.05.003

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 26

  LuehrN.UfimtsevI. S.MartinezT. J. (2011). Dynamic precision for electron repulsion integral evaluation on graphical processing units (GPUs). J. Chem. Theory Comput.7, 949–954. 10.1021/ct100701w

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 27

  MacKerellA. D.BashfordD.BellottM.DunbrackR. L.EvanseckJ. D.FieldM. J.et al. (1998). All-atom empirical potential for molecular modeling and dynamics studies of proteins. J. Phys. Chem. B102, 3586–3616. 10.1021/jp973084f

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 28

  MaierJ. A.MartinezC.KasavajhalaK.WickstromL.HauserK. E.SimmerlingC. (2015). ff14SB: Improving the accuracy of protein side chain and backbone parameters from ff99SB. J. Chem. Theory Comput.11, 3696–3713. 10.1021/acs.jctc.5b00255

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 29

  MorrisG. M.GoodsellD. S.HallidayR. S.HueyR.HartW. E.BelewR. K.et al. (1998). Automated docking using a Lamarckian genetic algorithm and an empirical binding free energy function. J. Comput. Chem.19, 1639–1662. 10.1002/(SICI)1096-987X(19981115)19:14<1639::AID-JCC10>3.0.CO;2-B

  • CrossRef
  • Google Scholar

• 30

  MucsD.BryceR. A. (2013). The application of quantum mechanics in structure-based drug design. Expert Opin. Drug Discov.8, 263–276. 10.1517/17460441.2013.752812

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 31

  RatcliffL. E.MohrS.HuhsG.DeutschT.MasellaM.GenoveseL. (2017). Challenges in large scale quantum mechanical calculations. Wiley Interdiscipl. Rev. Comput. Mol. Sci.7:24. 10.1002/wcms.1290

  • CrossRef
  • Google Scholar

• 32

  RosenzweigA. C.FrederickC. A.LippardS. J.NordlundP. (1993). Crystal structure of a bacterial non-haem iron hydroxylase thatcatalyses the biological oxidation of methane. Nature366, 537–543. 10.1038/366537a0

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 33

  Schrodinger (2012). PyMol: The PyMol Molecular Graphics System, Version 1.5.0.4. Schrodinger LLC.

  • Google Scholar

• 34

  SennH. M.ThielW. (2009). QM/MM methods for biomolecular systems. Angew. Chem. Int. Edn.48, 1198–1229. 10.1002/anie.200802019

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 35

  StehnJ. R.HaassN. K.BonelloT.DesouzaM.KottyanG.TreutleinH.et al. (2013). A novel class of anticancer compounds targets the actin cytoskeleton in tumor cells. Cancer Res.73, 5169–5182. 10.1158/0008-5472.Can-12-4501

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 36

  SukumarN.DasS. (2011). Current trends in virtual high throughput screening using ligand-based and structure-based methods. Comb. Chem. High Throughput Screen.14, 872–888. 10.2174/138620711797537120

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 37

  ThakurP. K.KumarJ.RayD.AnjumF.HassanM. I. (2013). Search of potential inhibitor against New Delhi metallo-beta-lactamase 1 from a series of antibacterial natural compounds. J. Nat. Sci. Biol. Med.4, 51–56. 10.4103/0976-9668.107260

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 38

  ThomannH.BernardoM.McCormickJ. M.PulverS.AnderssonK. K.LipscombJ. D.et al. (1993). Pulsed EPR studies of mixed valent [Fe(II)Fe(III)] forms of hemerythrin and methane monooxygenase: evidence for a hydroxide bridge. J. Am. Chem. Soc.115, 8881–8882. 10.1021/ja00072a068

  • CrossRef
  • Google Scholar

• 39

  ThomasP. W.ZhengM.WuS. S.GuoH.LiuD. L.XuD. G.et al. (2011). Characterization of purified New Delhi Metallo-beta-lactamase-1. Biochemistry50, 10102–10113. 10.1021/bi201449r

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 40

  TitovA. V.UfimtsevI. S.LuehrN.MartinezT. J. (2013). Generating efficient quantum chemistry codes for novel architectures. J. Chem. Theory Comput.9, 213–221. 10.1021/ct300321a

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 41

  TreutleinH. R.StylesM.ZengJ. (2017). Successful Computational Prediction of the Structure-Activity Relationship of a Potent JAK2 Inhibitor. 10.4225/03/5a1372ea64301

  • CrossRef
  • Google Scholar

• 42

  UfimtsevI. S.MartinezT. J. (2008). Quantum chemistry on graphical processing units. 1. Strategies for two-electron integral evaluation. J. Chem. Theory Comput.4, 222–231. 10.1021/ct700268q

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 43

  UfimtsevI. S.MartinezT. J. (2009a). Quantum chemistry on graphical processing units. 2. Direct self-consistent-field implementation. J. Chem. Theory Comput.5, 1004–1015. 10.1021/ct800526s

  • CrossRef
  • Google Scholar

• 44

  UfimtsevI. S.MartinezT. J. (2009b). Quantum chemistry on graphical processing units. 3. Analytical energy gradients, geometry optimization, and first principles molecular dynamics. J. Chem. Theory Comput.5, 2619–2628. 10.1021/ct9003004

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 45

  ValievM.BylaskaE. J.GovindN.KowalskiK.StraatsmaT. P.Van DamH. J. J.et al. (2010). NWChem: a comprehensive and scalable open-source solution for large scale molecular simulations. Comput. Phys. Commun.181, 1477–1489. 10.1016/j.cpc.2010.04.018

  • CrossRef
  • Google Scholar

• 46

  VogtL.Olivares-AmayaR.KermesS.ShaoY.Amador-BedollaC.Aspuru-GuzikA. (2008). Accelerating resolution-of-the-identity second-order Moller-Plesset quantum chemistry calculations with graphical processing units. J. Phys. Chem. A112, 2049–2057. 10.1021/jp0776762

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 47

  VogtM.StumpfeD.GeppertH.BajorathJ. (2010). Scaffold hopping using two-dimensional fingerprints: true potential, black magic, or a hopeless endeavor? Guidelines for virtual screening. J. Med. Chem.53, 5707–5715. 10.1021/jm100492z

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 48

  WangJ. C.LinJ. H.ChenC. M.PerrymanA. L.OlsonA. J. (2011). Robust scoring functions for protein-ligand interactions with quantum chemical charge models. J. Chem. Inf. Model.51, 2528–2537. 10.1021/ci200220v

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 49

  WangJ. M.WolfR. M.CaldwellJ. W.KollmanP. A.CaseD. A. (2004). Development and testing of a general amber force field. J. Comput. Chem.25, 1157–1174. 10.1002/jcc.20035

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 50

  WarshelA.LevittM. (1976). Theoretical studies of enzymic reactions: dielectric, electrostatic and steric stabilization of the carbonium ion in the reaction of lysozyme. J. Mol. Biol.103, 227–249. 10.1016/0022-2836(76)90311-9

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 51

  YasudaK. (2008). Accelerating density functional calculations with graphics processing unit. J. Chem. Theory Comput.4, 1230–1236. 10.1021/ct8001046

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 52

  YongD.TolemanM. A.GiskeC. G.ChoH. S.SundmanK.LeeK.et al. (2009). Characterization of a New Metallo-beta-Lactamase gene, bla(NDM-1), and a novel erythromycin esterase gene carried on a unique genetic structure in Klebsiella pneumoniae sequence type 14 from India. Antimicrob. Agents Chemother.53, 5046–5054. 10.1128/aac.00774-09

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 53

  ZengJ. (2000). Mini-review: computational structure-based design of inhibitors that target protein surfaces. Comb. Chem. High Throughput Screen.3, 355–362. 10.2174/1386207003331490

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 54

  ZengJ.TreutleinH. R. (1999). A method for computational combinatorial peptide design of inhibitors of Ras protein. Protein Eng.12, 457–468. 10.1093/protein/12.6.457

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 55

  ZhangH. M.HaoQ. (2011). Crystal structure of NDM-1 reveals a common beta-lactam hydrolysis mechanism. Faseb J.25, 2574–2582. 10.1096/fj.11-184036

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 56

  ZhengM.XuD. G. (2013). New Delhi Metallo-beta-lactamase I: substrate binding and catalytic mechanism. J. Phys. Chem. B117, 11596–11607. 10.1021/jp4065906

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 57

  ZhouT.HuangD. Z.CaflischA. (2010). Quantum mechanical methods for drug design. Curr. Top. Med. Chem.10, 33–45. 10.2174/156802610790232242

  • Pubmed Abstract
  • CrossRef
  • Google Scholar