key: cord-322562-3gvsn9vf authors: Hatada, Ryo; Okuwaki, Koji; Mochizuki, Yuji; Handa, Yuma; Fukuzawa, Kaori; Komeiji, Yuto; Okiyama, Yoshio; Tanaka, Shigenori title: Fragment Molecular Orbital Based Interaction Analyses on COVID-19 Main Protease − Inhibitor N3 Complex (PDB ID: 6LU7) date: 2020-06-15 journal: J Chem Inf Model DOI: 10.1021/acs.jcim.0c00283 sha: doc_id: 322562 cord_uid: 3gvsn9vf [Image: see text] The worldwide spread of COVID-19 (new coronavirus found in 2019) is an emergent issue to be tackled. In fact, a great amount of works in various fields have been made in a rather short period. Here, we report a fragment molecular orbital (FMO) based interaction analysis on a complex between the SARS-CoV-2 main protease (Mpro) and its peptide-like inhibitor N3 (PDB ID: 6LU7). The target inhibitor molecule was segmented into five fragments in order to capture site specific interactions with amino acid residues of the protease. The interaction energies were decomposed into several contributions, and then the characteristics of hydrogen bonding and dispersion stabilization were made clear. Furthermore, the hydration effect was incorporated by the Poisson–Boltzmann (PB) scheme. From the present FMO study, His41, His163, His164, and Glu166 were found to be the most important amino acid residues of Mpro in interacting with the inhibitor, mainly due to hydrogen bonding. A guideline for optimizations of the inhibitor molecule was suggested as well based on the FMO analysis. The 2019 novel coronavirus (SARS-CoV-2 or its disease name caused the pneumonia outbreak in Wuhan, China, in late December 2019 and has rapidly spread around the world. 1 As of March 19, 2020, around 190 000 people were infected and over 7700 have succumbed to the epidemic (as of May 25, 5 200 000 and 337 000, respectively). The World Health Organization (WHO) has declared this novel coronavirus outbreak a global health emergency, and subsequently has described the outbreak as pandemic on March 11. At present, there are no targeted therapeutics, and effective treatment options remain very limited. Although a large number of researchers around the world are engaged in developing antiviral drugs against COVID-19, it is also wellknown that new drug discovery and development is a long, costly, and rigorous scientific process. Many researchers investigating the genetic and functional data of SARS-CoV-2 compare it with other coronaviruses (CoVs) to design proper infection control strategy and seek potential drugs that can prevent and/or cure this serious epidemic. CoVs infect humans and other animal species, causing a variety of highly prevalent and severe diseases, including Severe Acute Respiratory Syndrome (SARS) and Middle East Respiratory Syndrome (MERS). CoVs are positive-sense, single-stranded RNA viruses featuring the largest viral RNA genomes known to date. The SARS-CoV-2 genome is comprised of approximately 30 000 nucleotides. As for the development of relevant drugs, viral proteases are wellknown common targets in dealing with human viruses such as the HIV virus and hepatitis C virus. Protease inhibitors would thus be remarkably effective in blocking the replication of coronavirus, including the SARS and the MERS, providing a promising foundation for the development of anticoronaviral therapeutics. Based on its functions, the main protease (Mpro) or chymotrypsin-like protease (3CLpro) is suggested to be a potential drug target to combat COVID-19, which is highly conserved among coronaviruses. Sequence alignment revealed that the Mpro of SARS-CoV-2 shares 96% similarity with that of previous SARS-CoV. 2 Studies for identifying the inhibitors of COVID-19 Mpro have then quickly begun for discovering and developing drugs against the disease. For example, on the basis of homology models for SARS-CoV-2 Mpro and docking simulations of 1903 small drug compounds, Xu et al. 3 selected four promising drug candidates according to the docking score, 3D similarity of binding mode to known ligands and the binding free energy evaluation. They thus suggested that nelfinavir can be a potential inhibitor against the SARS-CoV-2 Mpro. The X-ray crystal structure of SARS-CoV-2 Mpro was released on February 5, 2020 (and revised a week later) by Liu et al. 4 It is of dimer-like structure in complex with a peptidelike inhibitor called N3, 5 and very similar to that of previous SARS-CoV 6, 7 reflecting the 96% sequence homology between them. We can thus anticipate the development of dependable structure-based design of inhibitors against SARS-CoV-2 Mpro. In order to obtain reliable strategies for this aim, it is now desirable to know accurate information on the manner of ligand-protein interactions in the pharmacophore. As illustrated in a recent book of in silico drug design, 8 the fragment molecular orbital (FMO) method 9−12 provides an efficient tool for performing ab initio quantum-chemical calculations for biomolecular systems and accurately analyzing their intermolecular interactions in terms of the interfragment interaction energies (IFIEs). Namely, the set of IFIEs can be used as numerical indices to capture the nature of interactions in a given target system. In fact, the FMO method has been successfully applied to the comprehensive analyses of protein− ligand 13−16 and protein−protein 17−19 interactions, including viral molecular systems. The aim of the present work is to apply the FMO method to the analysis of the interactions between SARS-CoV-2 Mpro and its inhibitor N3 to elucidate the essence of its pharmacophore structure and associated molecular recognition. In the following sections, after addressing the present methodology by means of the FMO-IFIE, we show the calculated results on the interactions between the amino-acid residues of Mpro and the fragments of N3. With the aid of the advanced analyses based on the Pair Interaction Energy Decomposition Analysis (PIEDA) 20 and the incorporation of the solvent effect 21,22 as well, we will provide some useful information toward a rational drug design for the inhibitors against SARS-CoV-2 Mpro. In this study, we adopted the crystal structure of the complex 4 between COVID-19 Mpro and inhibitor N3 and processed it for the FMO based interaction analyses. The setup of calculations is described as follows. 2.1. Processing of the Crystal Structure. First, the crystal structure of the MPro:N3 complex (PDB ID: 6LU7; Version 2, 2.16 Å resolution, 4 ) was downloaded. As stated in the Introduction, although Mpro usually functions as a homodimer, we calculated only the monomeric unit stored in the PDB file since the dimer interface has no direct interaction with N3 ligand. We further processed the molecular structure with the MOE program. 23 All of the 84 crystal water molecules were preserved and subjected to subsequent FMO calculation. Figure 1 illustrates the peptide-like N3 inhibitor as the ligand to bind to the protein. Note that this ligand has a covalent bond to Cys145 as the result of the Michael addition. We performed standard modeling operations such as hydrogen attachment and pK a adaptation (pH 7). The charged N-and C-terminuses were assigned to Ser1 and Gln306, respectively. Molecular mechanics based energy minimization (with AMBER10:EHT force field 24 in MOE 23 ) was carried out for all the attached hydrogen atoms, the ligand, and amino acid residues within 4.5 Å of the ligand, where a tether mask of 1.0 Å was used for the crucial pharmacophore region. The entire structure of the processed protein−ligand complex and the close-up view of pharmacophore are shown in Figure 2 . Figure 3 is a schematic illustration of the ligand−residue interactions, based on geometrical relations judged by MOE. 23 This illustration may be informative as a guide in discussing the FMO results. The first term in this equation is the fictitiously isolated monomer energy, whereas the second term corresponds to pair interaction energy (PIE) 11 or interfragment interaction energy (IFIE), 12, 25 which is quite useful for interaction analyses 10, 12 (see refs 13−19 if necessary). We therefore used IFIE for the interaction analysis in this paper. In the FMO calculation of a protein, 10−12 the protein is usually fragmented at the C α carbon atom with a sp 3 hybridization 26 as exemplified in Figure 4 , where there is no hydrogen-capping for the bond segmentation via the technique of bond detachment atom (BDA). In this study, we regarded an amino acid residue as a fragment to maintain a one-to-one correspondence between a residue and a fragment. It should be noted that care is necessary when the >CO part in a peptide bond is involved in interactions because the formal assignment of reside is shifted as in Figure 4 27 (it will be addressed later). In the present case, as the N3 inhibitor is too large to be regarded as a fragment, we divided it into 5 fragments as shown in Figure 5 . This fragmentation pattern was similar to that often used to divide peptides. To be specific, each fragment in Figure 5 was designed so as to mimic an amino acid residue fragment of a peptide. We chose this pattern because the overall molecular structure of N3 is analogous to that of a peptide. The fragmentation of both Mpro and N3 was done with the FMOe program, 28 and the total number of fragments became 395 (=306 (protein) + 5 (ligand) + 84 (water)). Our FMO program ABINIT-MP 12 was used for the present study. The electron correlation effect was incorporated by the . Example of fragmentation for mini protein consisting of five residues (total number of fragments is thus five). The bonding electron pair at segmentation point is expressed as a colon symbol. Care is necessary for fac that the segmentation is not made at the peptide bond (refer to texts). second-order Møller−Plesset perturbation (MP2) calculation 29 which is efficiently executable with ABINIT-MP. 30−32 A partial renormalization 33 was utilized to reduce the trend of overstabilizations by MP2. The 6-31G* basis set 29 was adopted as a standard choice in FMO calculations. 9, 12 The decomposition analysis of PIE (called as PIEDA 20 ) can grasp the nature of interactions in detail. We thus adopted the PIEDA option in ABINIT-MP. 34 In PIEDA, an interaction energy, PIE or IFIE, between fragments I and J is expressed as where ES, EX, CT, and DI mean "ElectroStatic", "EXchange repulsion", "Charge-Transfer and mix terms", and DIspersion", respectively. The DI energy is usually evaluated by the MP2correlated calculation. Based on our experiences, contributions from ES and CT in eq 4 take large values for the hydrogen bonding, while the so-called CH/π 35 and π/π interactions have large DI portions. As the hydrogen-bonding is potentially influenced by the hydration effect, the Poisson−Boltzmann (PB) model 21,22 was employed to incorporate the hydration effect. The PB grand iteration was iterated until the reaction field of hydration converged to 1.0 × 10 −5 au in energy. Once converged, PIEDA was performed to obtain the respective contributions including the additional "Solvation" term. To identity the type of interaction, we introduced the following index with quantities of eq 4 (inspired by ref 36) where the "(Solv)" term in the denominator is valid only for the case of PB solvation. The quantity of eq 5 is multiplied by 100% (not necessarily normalized) for discussion in the next section. Note that a large value of F ES+CT suggests the existence of hydrogen bonding. Model without N3-Cys145 Bond. We additionally considered a noncovalent bonding model in which the N3 ligand is not connected to Cys145. This model was generated from the 6LU7 structure through the following steps with the MOE program; 23 (I) The covalent bond between N3 and Cys145 was cleaved to restore the α,β-unsaturated carbonyl reactant of the Michael addition reaction, (II) constrained optimization (tether = 1) was performed for the region around N3 (4.5 Å) by fixing the main chain of residues, and (III) selected optimization for N3, Cys145, and two connected residues. Note that the point of fragmentation in N3 was shifted due to the presence of double bond by which the Michael addition takes place. Finally, the FMO calculation was carried out without PB. Since the ligand moiety was divided into five fragments as shown in Figure 5 , we describe the interaction analyses in the order of fragment index. We will directly cite the evaluated values in the tables for detailed discussions. The corresponding bar graph representations of PIEDA are given in Figures S1− S5 (in the Supporting Information) for graphical understanding. 3.1. Fragment 1. Table 1 compiles the results of IFIE and decomposed contributions (PIEDA) interacting with Fragment 1 of the ligand, where the listing threshold is set as 2.0 kcal/mol The distance between the main chain >CO of Thr190 and the N−H part of Fragment 1 is as close as 1.99 Å as illustrated in Figure 6 , suggesting that they are forming a typical hydrogen bond. As addressed already, there is a concern of assignment shift 27 that this carbonyl part belongs to theAla191 fragment in the segmentation of FMO scheme. 10−12 The IFIE of Ala191 fragment is −11.58 kcal/mol. The ES term of the IFIE is −9.15 kcal/mol, and the CT term is −4.52 kcal/ mol, indicating an electrostatic interaction with charge transfer. The F ES+CT index is 69.5% certainly. These features indicate the existence of hydrogen bond. The IFIE and PIEDA values under the PB condition are given in Table 2 . For the hydrogen bonding on Ala191 (actually Thr190 residue), the hydration effect increments the stabilization. The sums of IFIE values (of all residues in the protease) for Fragment 1 with and without hydration are −26.20 and −24.21 kcal/mol, respectively. The interaction between Fragment 1 and protease can be characterized as hydrogen bonding type, consistent with the situation in Figure 3 . 23 3.2. Fragment 2. Fragment 2 has a Val-like side chain. Table 3 lists the IFIE and PIEDA results, and Figure 7 shows a couple of important geometrical features. There is a hydrogen bond between the main chain >CO part of Glu166 and the N−H part of Fragment 2, where the distance is found to be 1.91 Å. The assignment shift due to the fragmentation 27 appears again, and thus the IFIE of Leu167 fragment should be of interest. The corresponding value is −14.50 kcal/mol, and its F ES+CT index is 87.3%. Another notable stabilization by hydrogen bonding (−10.69 kcal/mol) is associated with the side chain of Gln189, where the associated distance is 2.44 Å as indicated in Figure 7 . Note that the importance of Glu166 and Gln189 is consistently suggested in Figure 3 . The PB results of IFIE and PIEDA are compiled in Table 4 . It is notable that the hydration enhances the stabilization for Leu167 fragment (actually Glu166 residue) but provides destabilization for Gln189 although their amounts are small. The sums of IFIE values for Fragment 2 with and without hydration are −26.92 and −28.41 kcal/mol, respectively. In other words, the PB hydration leads to a certain electrostatic shielding for Fragment 2. 3.3. Fragment 3. As shown in both Figures 3 and 8 , the main chain >N−H part of Glu166 and the >CO part of Fragment 3 have a hydrogen bonding, where the distance of interest is 1.88 Å and the corresponding IFIE value is −9.18 kcal/mol found in Table 5 . There is another hydrogen bonding with the side chain of Gln189 (distance of 1.90 Å), whose IFIE value is −10.96 kcal/mol. These two interactions have sizable F ES+CT values, as expected. Interestingly, several residues such as His41 and Met165 provide vital stabilizations with small F ES+CT values, suggesting that dispersion type interactions play an additional role in binding of Fragment 3. Both are hydrophobic interactions with the Leu-like side chain of Fragment 3; they are CH/π interaction with the imidazole ring of His41 (see Figure S6 ) and hydrophobic interaction with the thioether chain of Met165. Table 6 lists the PB results. The hydrogen bonding stabilizations of Glu166 and Gln189 are increased by the hydration. The sums of IFIE values for Fragment 3 with and without hydration are −49.45 and −47.50 kcal/mol, respectively. These values are about twice larger than the corresponding values of Fragment 1 and Fragment 2. 3.4. Fragment 4. Figure 3 suggests that His164 has a hydrogen bonding with the ligand, and this bonding is actually formed between its main chain >CO part and the >N−H part of Fragment 4 as illustrated in Figure 9 (distance of 1.87 Å). Due to a concern of assignment shift, 27 the corresponding stabilization is found as the IFIE value of −15.78 kcal/mol for Met165 fragment (see Table 7 ). Notably, the neighbored residue His163 has a large stabilization of −37.25 kcal/mol as IFIE, where the protonated side chain ring of His 163 plays a vital role in interactions with the >CO part of the ring in Fragment 4 (close distance of 1.74 Å). The F ES+CT index of His163 is as large as 89.2%, just indicating the hydrogen bonding. Glu166 provides a stabilization of −8.73 kcal/mol, and a relatively small F ES+CT value (57.0%) implies a contribution from dispersion. Phe140 and Asn142 have dispersion type stabilizations as well. The PB results are listed in Table 8 . The hydration leads to a stabilization for Met165 (actually His164 fragment). In contrast, Glu166 is rather destabilized via the electrostatic shielding (positive value of "Solv" term). Because the charged carboxylic acid in the side chain of Glu166 is exposed to the solvent, opposite-signed and sizable charges should be induced on its solvent surface. The interaction with these surface charges effectively suppresses the "ES" component between Fragment 4 and Glu166 as the "Solv" one. These can also be interpreted from more microscopic changes in electron density distribution using natural population analysis (NPA) 38,39 (refer to Figure S7 ). One carboxyl O atom of Glu166 and H of Fragment 4 that form a hydrogen bonding are depolarized with the magnitudes of 0.008e and 0.018e in atomic charge by the hydration, respectively. This weakens the "ES" interaction slightly than in vacuo. Another carboxyl O atom exposed to the solvent is rather polarized with the magnitude of 0.110e while inducing the solvent charges reactively. The presence of such a solvent-charge induction increases the significance of the "Solv" term because it incorporates the electrostatic shielding. Figure 3 , His41 can interact with Fragment 5. This suggestion is confirmed by the calculated results in Table 9 . A hydrogen bonding forms (distance of 2.16 Å in Figure 10 ) between the >N−H part of side chain ring of His41 and the >CO part of Fragment 5. The IFIE value is −10.57 kcal/mol, and consistently the F ES+CT value is 83.4%. CH/π interaction between phenyl ring of Fragment 5 and Thr25 is also observed ( Figure S6 ), whose main energy component is DI (−3.17 kcal/mol). The PB hydration increases the stabilization for His41 as shown in Table 10 . The sums of IFIE values for Fragment 5 with and without hydration are −17.77 and −21.32 kcal/mol, respectively, indicating an electrostatic shielding. The role of Fragment 5 is observed to be rather minor relative to other four fragments. 3.6. Effect of Elimination of Covalent Bond between Cys145 and N3. Figure 11 compares the MOE-based optimized structures around the N3 ligand with and without the covalent bond to Cys145. As illustrated in this figure, most of the structural changes before and after the Michael addition reaction took place on the ligand, while the protein was almost unchanged except for Cys145, where the S−C bond distance of Cys145 was 1.83 Å with the bond but 3.22 Å without it. The IFIE and PIEDA results for Cys145 are compiled in Table S1 for nonbonded Cys145, where the fragmentation points are illustrated in Figure S8 . Interestingly, the interaction energies between the nonbonded Cys145 and Fragments 4 and 5 are not large before the Michael reaction. The remaining Fragments 1−3 show almost no interaction with the protein. For Fragments 4 and 5, the interaction energies are listed for important residues in Tables S2 and S3, respectively. It is notable that the conformational changes of the ligand lead to destabilization of the hydrogen bonding between Fragment 4 and the main chain > CO of His164 (ΔE of +5.86 kcal/mol, relative to the corresponding entry of Met165 in Table 7 ) and also stabilization (ΔE of −11.63 kcal/mol, similarly) of that with Gln166. About Fragment 5, there was a sizable conformational change for the terminal Phe ring, which destabilized the interaction with Thr25 (+2.00 kcal/mol, relative to the value in Table 9 ) and that with Leu27 (+3.39 kcal/mol) mainly due to the loss of DI. For Fragments 1−3, no table difference was found with and without the covalent bond. 3.7. Summary of Interaction Analyses. In this section, a series of interaction analyses have been made in fragment-byfragment fashion. Here, we would make several notes as follows. First, all five fragments of the N3 ligand have the hydrogen bonding with neighbored residues in the pharmacophore. In particular, Fragment 4 interacts strongly with His163, His164, and Glu166, and this fragment plays a leading role in stabilizing the ligand. Besides these three residues, His41, Gln189, and Thr190 have hydrogen bonds with fragments of the ligand. Second, the dispersion-type interactions assist the ligand binding (e.g., Phe140, Asn142, and Glu166 interacting with Fragment 4). Third, the phenyl ring of Fragment 5 and Journal of Chemical Information and Modeling pubs.acs.org/jcim Article isoxazole ring of Fragment 1 do not have prominent interaction with surrounding residues. This fact suggests that if some chemical modifications are made on this ring to further optimize the inhibitor, interactions with Thr25, Thr26, and Asn142 for Fragment 5, and Ala191 and Pro168 for Fragment 1 may be induced. Fourth, the role of the peptide-like side chains in the ligand is illuminated. For example, the Leu-like side chain of Fragment 3 has a CH/π interaction with His41. His41 has a CH /π bond to Fragment 3 and also a hydrogen bond to fragment 5, and thus plays an important role in ligand binding. Supplemental modifications on the groups such as Ala-like side chain (currently not interacting with the neighboring residues) in Fragment 1 may also be effective. These findings obtained by the FMO calculation were in general agreement with the consideration given by X-ray crystallography. 4 Our discussion with quantitative interaction energy values could provide additional quantitative insights. In addition to confirming and reinforcing the experimental observations, the FMO calculation has newly found that His41 plays an important role in N3 binding. In the present work, a series of FMO-based interaction analyses 10,12 with ABINIT-MP 12 have been made on the complex between the COVID-19 main protease and the peptide-like inhibitor N3, whose fundamental structure was obtained from PDB ID: 6LU7. 4 The object ligand was divided into five fragments, and their interactions were investigated by using the IFIE 12,25 and PIEDA 34 with and without the PB hydration. 21, 22 It was found that the hydrogen bonds with surrounding residues in the pharmacophore were important throughout five fragments. The fourth fragment of the ligand was identified as the most important moiety in interactions with His163, His164 and Glu166. It was also noted that dispersion interactions provided additional stabilization of the ligand, such as the CH/π interaction with His41. A possibility of further optimization of ligand binding was suggested at the The dimer-ES approximation 37 was adopted due to a long distance from this fragment. fifth fragment of the ligand for the interactions with Thr25, Thr26 and Asn142 in addition to His41. The influence of the covalent bond between the ligand and Cys145 was also investigated. Though some differences were observed for several residues, essential discussion on interaction energies were kept unchanged. It would be noted technically that these biochemical findings based on the IFIE and PIEDA results were carefully derived by avoiding the issue of assignment shift due to C α fragmentation; 27 the function of peptide bond fragmentation 40 has been recently available in ABINIT-MP, but a cross-reference testing should be required before routine usages. Virtual screening and high-throughput screening based on the crystal structure of the complex of Mpro and N3 have recently been performed. 4 In ref 4, the authors have discussed why cinanserin is a potential inhibitor for SARS-CoV-2. The present FMO analysis supports this idea because the significance of attractive interactions between the inhibitor and surrounding residues such as His41, Met49, Asn142, Met165 and Glu166 were identified quantitatively. Compounds that mimic many of these interactions are then expected to have higher inhibitory activity, indicating a direction of HIT-to-lead process. As a whole, our quantumchemical analysis brings additional quantitative information to the recent comprehensive study 4 and can provide useful insights into the rational drug discovery based on the N3 pharmacophore. Thus, we hope that the present work by the Figure 11 . Pharmacophore structures around the N3 ligand with and without covalent bond to Cys145. Green and light blue frames correspond to the ligands with and without covalent bond, respectively. FMO calculations for the SARS-CoV-2 Mpro inhibitor could contribute to the worldwide efforts toward the development of effective drugs against the COVID-19. 41 From a stand-point of open science promotion against COVID-19, we have uploaded the data set of IFIE and PIEDA on the present system to the FMO database (FMODB). 42, 43 ■ ASSOCIATED CONTENT * sı Supporting Information The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acs.jcim.0c00283. Additional experimental results as found in supplemental figures and tables (PDF) Global Spread of COVID-19 and Pandemic Potential Evolution of the Novel Coronavirus from the Ongoing Wuhan Outbreak and Modeling of Its Spike Protein for Risk of Human Transmission Nelfinavir Was Predicted to Be a Potential Inhibitor of 2019-NCov Main Protease by an Integrative Approach Combining Homology Modelling, Molecular Docking and Binding Free Energy Calculation Structure of M pro from SARS-CoV-2 and discovery of its inhibitors Design of Wide-Spectrum Inhibitors Targeting Coronavirus Main Proteases Coronavirus Main Proteinase (3CL pro ) Structure: Basis for Design of Anti-SARS Drugs The Crystal Structures of Severe Acute Respiratory Syndrome Virus Main Protease and Its Complex with an Inhibitor Quantum Mechanics in Drug Discovery The Fragment Molecular Orbital Method: Practical Applications to Large Molecular Systems Fragment Molecular Orbital Method: An Approximate Computational Method for Large Molecules Exploring Chemistry with the Fragment Molecular Orbital Method Electron-Correlated Fragment-Molecular-Orbital Calculations for Biomolecular and Nano Systems Molecular Interactions between Estrogen Receptor and Its Ligand Studied by the Ab Initio Fragment Molecular Orbital Method Theoretical Analysis of Binding Specificity of Influenza Viral Hemagglutinin to Avian and Human Receptors Based on the Fragment Molecular Orbital Method Interaction Energy Analysis on Specific Binding of Influenza Virus Hemagglutinin to Avian and Human Sialosaccharide Receptors: Importance of Mutation-Induced Structural Change Towards Good Correlation between Fragment Molecular Orbital Interaction Energies and Experimental IC(50) for Ligand Binding: A Case Study of P38 MAP Kinase Possibility of Mutation Prediction of Influenza Hemagglutinin by Combination of Hemadsorption Experiment and Quantum Chemical Calculation for Antibody Binding Prediction of Probable Mutations in Influenza Virus Hemagglutinin Protein Based on Large-Scale Ab Initio Fragment Molecular Orbital Calculations Computational Analysis of the Interaction Energies between Amino Acid Residues of the Measles Virus Hemagglutinin and Its Receptors Fragment Molecular Orbital Calculations with Implicit Solvent Based on the Poisson-Boltzmann Equation: Implementation and DNA Study Fragment Molecular Orbital Calculations with Implicit Solvent Based on the Poisson-Boltzmann Equation: II. Protein and Its Ligand-Binding System Studies A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules VISCANA: Visualized Cluster Analysis of Protein-Ligand Interaction Based on the Ab Initio Fragment Molecular Orbital Method for Virtual Ligand Screening Fragment Molecular Orbital Method: Application to Polypeptides Antigen-antibody interactions of influenza virus hemagglutinin revealed by the fragment molecular orbital calculation Modern Quantum Chemistry Large Scale MP2 Calculations with Fragment Molecular Orbital Scheme A Parallelized Integral-Direct Second-Order Møller-Plesset Perturbation Theory Method with a Fragment Molecular Orbital Scheme Large Scale FMO-MP2 Calculations on a Massively Parallel-Vector Computer Enhanced Second-Order Treatment of Electron Pair Correlation Implementation of Pair Interaction Energy Decomposition Analysis and Its Applications to Protein-Ligand Systems CH// Interaction in the Crystal Structure of Organic Compounds. A Database Study Characterising Inter-Helical Interactions of G Protein-Coupled Receptors with the Fragment Molecular Orbital Method Fragment molecular orbital method: use of approximate electrostatic potential Natural Bond Orbital Analysis of Near-Hartree−Fock Water Dimer Natural Population Analysis Fragmentation at sp 2 carbon in fragment molecular orbital (FMO) method Research and Development on Therapeutic Agents and Vaccines for COVID-19 and Related Human Coronavirus Diseases FMO database (FMODB) The database of quantum mechanical data based on the FMO method Development of an automated fragment molecular orbital (FMO) calculation protocol toward construction of quantum mechanical calculation database for large biomolecules The authors declare no competing financial interest.