Biomedical subjects

Art E Cho

Publications and source records attributed to Art E Cho.

4 recordsLinked to original sources

PubMed2005

Importance of accurate charges in molecular docking: quantum mechanical/molecular mechanical (QM/MM) approach.

The extent to which accuracy of electric charges plays a role in protein-ligand docking is investigated through development of a docking algorithm, which incorporates quantum mechanical/molecular mechanical (QM/MM) calculations. In this algorithm, fixed charges of ligands obtained from force field parameterization are replaced by QM/MM calculations in the protein environment, treating only the ligands as the quantum region. The algorithm is tested on a set of 40 cocrystallized structures taken from the Protein Data Bank (PDB) and provides strong evidence that use of nonfixed charges is important. An algorithm, dubbed "Survival of the Fittest" (SOF) algorithm, is implemented to incorporate QM/MM charge calculations without any prior knowledge of native structures of the complexes. Using an iterative protocol, this algorithm is able in many cases to converge to a nativelike structure in systems where redocking of the ligand using a standard fixed charge force field exhibits nontrivial errors. The results demonstrate that polarization effects can play a significant role in determining the structures of protein-ligand complexes, and provide a promising start towards the development of more accurate docking methods for lead optimization applications.

Algorithms↗

PubMed2005

The MPSim-Dock hierarchical docking algorithm: application to the eight trypsin inhibitor cocrystals.

To help improve the accuracy of protein-ligand docking as a useful tool for drug discovery, we developed MPSim-Dock, which ensures a comprehensive sampling of diverse families of ligand conformations in the binding region followed by an enrichment of the good energy scoring families so that the energy scores of the sampled conformations can be reliably used to select the best conformation of the ligand. This combines elements of DOCK4.0 with molecular dynamics (MD) methods available in the software, MPSim. We test here the efficacy of MPSim-Dock to predict the 64 protein-ligand combinations formed by starting with eight trypsin cocrystals, and crossdocking the other seven ligands to each protein conformation. We consider this as a model for how well the method would work for one given target protein structure. Using as a criterion that the structures within 2 kcal/mol of the top scoring include a conformation within a coordinate root mean square (CRMS) of 1 A of the crystal structure, we find that 100% of the 64 cases are predicted correctly. This indicates that MPSim-Dock can be used reliably to identify strongly binding ligands, making it useful for virtual ligand screening.

Algorithms↗

PubMed2005

Integrated Modeling Program, Applied Chemical Theory (IMPACT).

We provide an overview of the IMPACT molecular mechanics program with an emphasis on recent developments and a description of its current functionality. With respect to core molecular mechanics technologies we include a status report for the fixed charge and polarizable force fields that can be used with the program and illustrate how the force fields, when used together with new atom typing and parameter assignment modules, have greatly expanded the coverage of organic compounds and medicinally relevant ligands. As we discuss in this review, explicit solvent simulations have been used to guide our design of implicit solvent models based on the generalized Born framework and a novel nonpolar estimator that have recently been incorporated into the program. With IMPACT it is possible to use several different advanced conformational sampling algorithms based on combining features of molecular dynamics and Monte Carlo simulations. The program includes two specialized molecular mechanics modules: Glide, a high-throughput docking program, and QSite, a mixed quantum mechanics/molecular mechanics module. These modules employ the IMPACT infrastructure as a starting point for the construction of the protein model and assignment of molecular mechanics parameters, but have then been developed to meet specialized objectives with respect to sampling and the energy function.

Computer Simulation↗

PubMed2004

Broken-symmetry unrestricted hybrid density functional calculations on nickel dimer and nickel hydride.

In the present work we investigate the adequacy of broken-symmetry unrestricted density functional theory for constructing the potential energy curve of nickel dimer and nickel hydride, as a model for larger bare and hydrogenated nickel cluster calculations. We use three hybrid functionals: the popular B3LYP, Becke's newest optimized functional Becke98, and the simple FSLYP functional (50% Hartree-Fock and 50% Slater exchange and LYP gradient-corrected correlation functional) with two basis sets: all-electron (AE) Wachters+f basis set and Stuttgart RSC effective core potential (ECP) and basis set. We find that, overall, the best agreement with experiment, comparable to that of the high-level CASPT2, is obtained with B3LYP/AE, closely followed by Becke98/AE and Becke98/ECP. FSLYP/AE and B3LYP/ECP give slightly worse agreement with experiment, and FSLYP/ECP is the only method among the ones we studied that gives an unacceptably large error, underestimating the dissociation energy of Ni(2) by 28%, and being in the largest disagreement with the experiment and the other theoretical predictions. We also find that for Ni(2), the spin projection for the broken-symmetry unrestricted singlet states changes the ordering of the states, but the splittings are less than 10 meV. All our calculations predict a deltadelta-hole ground state for Ni(2) and delta-hole ground state for NiH. Upon spin projection of the singlet state of Ni(2), almost all of our calculations: Becke98 and FSLYP both AE and ECP and B3LYP/AE predict (1)(d(A)(x(2)-y(2)d(B)(x(2)-y(2)) or (1)(d(A)(xy) (d)(B)(xy)) ground state, which is a mixture of (1)Sigma(g) (+) and (1)Gamma(g). B3LYP/ECP predicts a (3)(d(A)(x(2)-y(2))d(B)(xy) (mixture of (3)Sigma(g) (-) and (3)Gamma(u)) ground state virtually degenerate with the (1)(d(A)(x(2)-y(2)d(B)(x)(2)-y(2)/(1)(d(A)(xy)D(B)(xy) state. The doublet delta-hole ground state of NiH predicted by all our calculations is in agreement with the experimentally predicted (2)Delta ground state. For Ni(2), all our results are consistent with the experimentally predicted ground state of 0(g) (+) (a mixture of (1)Sigma(g) (+) and (3)Sigma(g) (-)) or 0(u) (-) (a mixture of (1)Sigma(u) (-) and (3)Sigma(u) (+)).

Journal Article↗