The ONIOM Method and Its Applications

A simple and efficient dispersion correction to the Hartree-Fock theory (2): Incorporation of a geometrical correction for the basis set superposition error

One of the most challenging problems in computer-aided drug discovery is the accurate prediction of the binding energy between a ligand and a protein. For accurate estimation of net binding energy ΔEbind in the framework of the Hartree-Fock (HF) theory, it is necessary to estimate two additional energy terms; the dispersion interaction energy (Edisp) and the basis set superposition error (BSSE). We previously reported a simple and efficient dispersion correction, Edisp, to the Hartree-Fock theory (HF-Dtq). In the present study, an approximation procedure for estimating BSSE proposed by Kruse and Grimme, a geometrical counterpoise correction (gCP), was incorporated into HF-Dtq (HF-Dtq-gCP). The relative weights of the Edisp (Dtq) and BSSE (gCP) terms were determined to reproduce ΔEbind calculated with CCSD(T)/CBS or /aug-cc-pVTZ (HF-Dtq-gCP (scaled)). The performance of HF-Dtq-gCP (scaled) was compared with that of B3LYP-D3(BJ)-bCP (dispersion corrected B3LYP with the Boys and Bernadi counterpoise correction (bCP)), by taking ΔEbind (CCSD(T)-bCP) of small non-covalent complexes as 'a golden standard'. As a critical test, HF-Dtq-gCP (scaled)/6-31G(d) and B3LYP-D3(BJ)-bCP/6-31G(d) were applied to the complex model for HIV-1 protease and its potent inhibitor, KNI-10033. The present results demonstrate that HF-Dtq-gCP (scaled) is a useful and powerful remedy for accurately and promptly predicting ΔEbind between a ligand and a protein, albeit it is a simple correction procedure.

Insights into the Nature of Anesthetic-Protein Interactions: An ONIOM Study

Anesthetics have been employed widely to relieve surgical suffering, but their mechanism of action is not yet clear. For over a century, the mechanism of anesthesia was previously thought to be via lipid bilayer interactions. In the present work, a rigorous three-layer ONIOM(M06-2X/6-31+G*:PM6:AMBER) method was utilized to investigate the nature of interactions between several anesthetics and actual protein binding sites. According to the calculated structural features, interaction energies, atomic charges, and electrostatic potential surfaces, the amphiphilic nature of anesthetic-protein interactions was demonstrated for both inhalational and injectable anesthetics. The existence of hydrogen and halogen bonding interactions between anesthetics and proteins was clearly identified, and these interactions served to assist ligand recognition and binding by the protein. Within all complexes of inhalational or injectable anesthetics, the polarization effects play a dominant role over the steric effects and induce a significant asymmetry in the otherwise symmetric atomic charge distributions of the free ligands in vacuo. This study provides new insight into the mechanism of action of general anesthetics in a more rigorous way than previously described. Future rational design of safer anesthetics for an aging and more physiologically vulnerable population will be predicated on this greater understanding of such specific interactions.

In silico prediction and screening of modular crystal structures via a high-throughput genomic approach

High-throughput computational methods capable of predicting, evaluating and identifying promising synthetic candidates with desired properties are highly appealing to today's scientists. Despite some successes, in silico design of crystalline materials with complex three-dimensionally extended structures remains challenging. Here we demonstrate the application of a new genomic approach to ABC-6 zeolites, a family of industrially important catalysts whose structures are built from the stacking of modular six-ring layers. The sequences of layer stacking, which we deem the genes of this family, determine the structures and the properties of ABC-6 zeolites. By enumerating these gene-like stacking sequences, we have identified 1,127 most realizable new ABC-6 structures out of 78 groups of 84,292 theoretical ones, and experimentally realized 2 of them. Our genomic approach can extract crucial structural information directly from these gene-like stacking sequences, enabling high-throughput identification of synthetic targets with desired properties among a large number of candidate structures.

High-throughput computation aids design of new functional materials. Here, Yu et al. develop a high-throughput screening method for a group of zeolites with crystalline modular structures which are viewed as having gene-like stacking codes, and identify the most promising structures with desired properties.

Requirement for transient metal ions revealed through computational analysis for DNA polymerase going in reverse

DNA polymerases facilitate faithful insertion of nucleotides, a central reaction occurring during DNA replication and repair. DNA synthesis (forward reaction) is "balanced," as dictated by the chemical equilibrium by the reverse reaction of pyrophosphorolysis. Two closely spaced divalent metal ions (catalytic and nucleotide-binding metals) provide the scaffold for these reactions. The catalytic metal lowers the pKa of O3' of the growing primer terminus, and the nucleotide-binding metal facilitates substrate binding. Recent time-lapse crystallographic studies of DNA polymerases have identified an additional metal ion (product metal) associated with pyrophosphate formation, leading to the suggestion of its possible involvement in the reverse reaction. Here, we establish a rationale for a role of the product metal using quantum mechanical/molecular mechanical calculations of the reverse reaction in the confines of the DNA polymerase β active site. Additionally, site-directed mutagenesis identifies essential residues and metal-binding sites necessary for pyrophosphorolysis. The results indicate that the catalytic metal site must be occupied by a magnesium ion for pyrophosphorolysis to occur. Critically, the product metal site is occupied by a magnesium ion early in the pyrophosphorolysis reaction path but must be removed later. The proposed dynamic nature of the active site metal ions is consistent with crystallographic structures. The transition barrier for pyrophosphorolysis was estimated to be significantly higher than that for the forward reaction, consistent with kinetic activity measurements of the respective reactions. These observations provide a framework to understand how ions and active site changes could modulate the internal chemical equilibrium of a reaction that is central to genome stability.

A modelling approach for MOF-encapsulated metal catalysts and application to n-butane oxidation

A surrogate pore mimics the steric constraints provided by a MOF on n-butane in a MOF-encapsulated metal catalyst.

Computational study on the mechanism and enantioselectivity of Rh2(S-PTAD)4 catalyzed asymmetric [4+3] cycloaddition between vinylcarbenoids and dienes

In this paper, the mechanism of asymmetric [4+3] cycloaddition between a vinylcarbenoid and a diene to form cycloheptadiene has been studied using a two-layer ONIOM methodology consisting of density functional theory and semiempirical PM6.

An explicit account of solvation is essential for modeling Suzuki–Miyaura coupling in protic solvents

We compared explicit and implicit solvation approaches in modeling the free energy profile of the final step of Suzuki–Miyaura coupling.