Abstract
The fragment molecular orbital (FMO) method is an efficient quantum chemical calculation technique for large biomolecules, dividing each into smaller fragments and providing inter-fragment interaction energies (IFIEs) that support our understanding of molecular recognition. The ab initio fragment MO method program (ABINIT-MP), an FMO processing software, can automatically divide typical proteins and nucleic acids. In contrast, small molecules such as ligands and hetero systems must be manually divided. Thus, we developed a graphical user interface to easily handle such manual fragmentation as a library for Molecular Operating Environment (MOE) that preprocesses and visualizes FMO calculations. We demonstrated fragmentation with IFIE analyses for the two following cases: 1) covalent cysteine–ligand bonding inside the SARS-CoV-2 main protease (Mpro) and nirmatrelvir (Paxlovid) complex, and 2) the metal coordination inside a zinc-bound cyclic peptide. IFIE analysis successfully identified the key amino acid residues for the molecular recognition of nirmatrelrvir with Mpro and the details of their interactions (e.g., hydrogen bonds and CH/π interactions) via ligand fragmentation of functional group units. In metalloproteins, we found an efficient and accurate scheme for the fragmentation of Zn2+ ions with four histidines coordinated to the ion. FMOe simplifies manual fragmentation, allowing users to experiment with various fragmentation patterns and perform in-depth IFIE analysis with high accuracy. In the future, our findings will provide valuable insight into complicated cases, such as ligand fragmentation in modality drug discovery, especially for medium-sized molecules and metalloprotein fragmentation around metals.
Supplementary materials
Title
Supporting Information
Description
The interaction energy analysis between nirmatrelvir (9I7) and each amino acid residue of SARS-CoV-2 Mpro using ligand fragmentation is displayed in Figures S1 and S2. Interaction energy analysis of Zn2+ ion fragment, including the four histidines for the cyclic peptide containing Zn2+ ion, is shown in Figures S3 and S4 based on the main chain and the main/side chain fragmentation schema, respectively. Atomic charges of His6, His18, D-His12 (DHI12), D-His24 (DHI24), and Zn2+ ion using Mulliken, natural population analysis, and Merz–Kollman approaches at the FMO-HF/6-31G and FMO-HF/6-31G* levels are listed in Tables S1 and S2, respectively. Mulliken, natural population analysis, and Merz–Kollman approaches for each fragmentation scheme at the HF/6-31G* level are displayed in Figures S5, S6, and S7, respectively. FMO-based electron density analysis of a Zn2+ ion and its coordinated histidines is demonstrated in Figure S8 (PDF).
Actions