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.

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.