FMOe: Preprocessing and Visualizing Package of the Fragment Molecular Orbital Method for Molecular Operating Environment and Its Applications in Covalent Ligand and Metalloprotein Analyses.

The fragment molecular orbital (FMO) method is an efficient quantum chemical calculation technique for large biomolecules, dividing each into smaller fragments and providing interfragment interaction energies (IFIEs) that support our understanding of molecular recognition. The ab initio fragment MO method (ABINIT-MP), an FMO processing program, can automatically divide typical proteins and nucleic acids. In contrast, small molecules such as ligands and heterosystems must be manually divided. Thus, we developed a graphical user interface to easily handle such manual fragmentation as a library for the 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 nirmatrelvir 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.

New antimalarials identified by a cell-based phenotypic approach: Structure-activity relationships of 2,3,4,9-tetrahydro-1H-ß-carboline derivatives possessing a 2-((coumarin-5-yl)oxy)alkanoyl moiety.

The identification, structure-activity relationships (SARs), and biological effects of new antimalarials consisting of a 2,3,4,9-tetrahydro-1H-ß-carboline core, a coumarin ring, and an oxyalkanoyl linker are described. A cell-based phenotypic approach was employed in this search for novel antimalarial drugs with unique modes of action. Our screening campaign of the RIKEN compound library succeeded in the identification of the known tetrahydro-ß-carboline derivative (4e) as a hit compound showing significant in vitro activity. SAR studies on this chemical series led to the discovery of compound 4h having a (R)-methyl group on the oxyacetyl linker with potent inhibition of parasite growth (IC50 = 2.0 nM). Compound 4h was also found to exhibit significant in vivo antimalarial effects in mouse models. Furthermore, molecular modeling studies on 4e, 4h, and its diastereomer (4j) suggested that the (R)-methyl group of 4h forces the preferential adoption of a specific conformer which is considered to be an active conformer.

Special Features of COVID-19 in the FMODB: Fragment Molecular Orbital Calculations and Interaction Energy Analysis of SARS-CoV-2-Related Proteins.

SARS-CoV-2 is the causative agent of coronavirus (known as COVID-19), the virus causing the current pandemic. There are ongoing research studies to develop effective therapeutics and vaccines against COVID-19 using various methods and many results have been published. The structure-based drug design of SARS-CoV-2-related proteins is promising, however, reliable information regarding the structural and intra- and intermolecular interactions is required. We have conducted studies based on the fragment molecular orbital (FMO) method for calculating the electronic structures of protein complexes and analyzing their quantitative molecular interactions. This enables us to extensively analyze the molecular interactions in residues or functional group units acting inside the protein complexes. Such precise interaction data are available in the FMO database (FMODB) (https://drugdesign.riken.jp/FMODB/). Since April 2020, we have performed several FMO calculations on the structures of SARS-CoV-2-related proteins registered in the Protein Data Bank. We have published the results of 681 structures, including three structural proteins and 11 nonstructural proteins, on the COVID-19 special page (as of June 8, 2021). In this paper, we describe the entire COVID-19 special page of the FMODB and discuss the calculation results for various proteins. These data not only aid the interpretation of experimentally determined structures but also the understanding of protein functions, which is useful for rational drug design for COVID-19.

FMODB: The World's First Database of Quantum Mechanical Calculations for Biomacromolecules Based on the Fragment Molecular Orbital Method.

We developed the world's first web-based public database for the storage, management, and sharing of fragment molecular orbital (FMO) calculation data sets describing the complex interactions between biomacromolecules, named FMO Database (https://drugdesign.riken.jp/FMODB/). Each entry in the database contains relevant background information on how the data was compiled as well as the total energy of each molecular system and interfragment interaction energy (IFIE) and pair interaction energy decomposition analysis (PIEDA) values. Currently, the database contains more than 13 600 FMO calculation data sets, and a comprehensive search function implemented at the front-end. The procedure for selecting target proteins, preprocessing the experimental structures, construction of the database, and details of the database front-end were described. Then, we demonstrated a use of the FMODB by comparing IFIE value distributions of hydrogen bond, ion-pair, and XH/π interactions obtained by FMO method to those by molecular mechanics approach. From the comparison, the statistical analysis of the data provided standard reference values for the three types of interactions that will be useful for determining whether each interaction in a given system is relatively strong or weak compared to the interactions contained within the data in the FMODB. In the final part, we demonstrate the use of the database to examine the contribution of halogen atoms to the binding affinity between human cathepsin L and its inhibitors. We found that the electrostatic term derived by PIEDA greatly correlated with the binding affinities of the halogen containing cathepsin L inhibitors, indicating the importance of QM calculation for quantitative analysis of halogen interactions. Thus, the FMO calculation data in FMODB will be useful for conducting statistical analyses to drug discovery, for conducting molecular recognition studies in structural biology, and for other studies involving quantum mechanics-based interactions.

High-Precision Atomic Charge Prediction for Protein Systems Using Fragment Molecular Orbital Calculation and Machine Learning.

Here, we have constructed neural network-based models that predict atomic partial charges with high accuracy at low computational cost. The models were trained using high-quality data acquired from quantum mechanics calculations using the fragment molecular orbital method. We have succeeded in obtaining highly accurate atomic partial charges for three representative molecular systems of proteins, including one large biomolecule (approx. 2000 atoms). The novelty of our approach is the ability to take into account the electronic polarization in the system, which is a system-dependent phenomenon, being important in the field of drug design. Our high-precision models are useful for the prediction of atomic partial charges and expected to be widely applicable in structure-based drug designs such as structural optimization, high-speed and high-precision docking, and molecular dynamics calculations.

Quantum Chemical Interaction Analysis between SARS-CoV-2 Main Protease and Ensitrelvir Compared with Its Initial Screening Hit.

A non-covalent oral drug targeting SARS-CoV-2 main protease (Mpro), ensitrelvir (Xocova), has been developed using structure-based drug design (SBDD). To elucidate the factors responsible for enhanced inhibitory activities from an in silico screening hit compound to ensitrelvir, we analyzed the interaction energies of the inhibitors with each residue of Mpro using fragment molecular orbital (FMO) calculations. This analysis reveals that functional group conversion for P1' and P1 parts in the inhibitors increases the strength of existing interactions with Mpro and also provides novel interactions for ensitrelvir; the associated changes in the conformation of Mpro induce further interactions for ensitrelvir in other parts, including hydrogen bonds, a halogen bond, and π-orbital interactions. Thus, we illuminate the promising strategies of SBDD for leading ensitrelvir to get higher activity against Mpro by elucidating microscopic interactions through FMO-based analysis. These detailed mechanism findings, including water cross-linkings, will help to design novel inhibitors in SBDD.