Towards good correlation between fragment molecular orbital interaction energies and experimental IC50 for ligand binding: A case study of p38 MAP kinase

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.

Toward Post-Hartree-Fock Accuracy for Protein-Ligand Affinities Using the Molecules-in-Molecules Fragmentation-Based Method

The complexity and size of large molecular systems, such as protein-ligand complexes, pose computational challenges for accurate post-Hartree-Fock calculations. This study delivers a thorough benchmarking of the Molecules-in-Molecules (MIM) method, presenting a clear and accessible strategy for layer/theory selections in post-Hartree-Fock computations on substantial molecular systems, notably protein-ligand complexes. An approach is articulated, enabling augmented computational efficiency by strategically canceling out common subsystem energy terms between complexes and proteins within the supermolecular equation. Employing DLPNO-based post-Hartree-Fock methods in conjunction with the three-layer MIM method (MIM3), this study demonstrates the achievement of protein-ligand binding energies with remarkable accuracy (errors <1 kcal mol-1), while significantly reducing computational costs. Furthermore, noteworthy correlations between theoretically computed interaction energies and their experimental equivalents were observed, with R2 values of approximately 0.90 and 0.78 for CDK2 and BZT-ITK sets, respectively, thus validating the efficacy of the MIM method in calculating binding energies. By highlighting the crucial role of diffuse or small Pople-style basis sets in the middle layer for reducing energy errors, this work provides valuable insights and practical methodologies for interaction energy computations in large molecular complexes and opens avenues for their application across a diverse range of molecular systems.

Identification of IL-2 inducible tyrosine kinase inhibitors by quantum mechanics and ligand based virtual screening approaches

Interleukin-2-inducible T-cell kinase (ITK) is a crucial intracellular signaling mediator in normal and malignant T-cells and natural killer cells. Selective inhibition of ITK might be useful for treating a variety of disorders including; autoimmune, inflammatory, and neoplastic disorders. Over the past two decades, the clinical management of ITK inhibitors has progressed dramatically. So far, specific inhibitor with no off-target effects against ITK is available. Herein, we aim to discover potential virtual hits to fasten the process of drug design and development against ITK. In this regard, the key chemical characteristics of ITK inhibitors were identified using ligand-based pharmacophore modeling. The validated pharmacophore comprises one hydrogen bond donor and three hydrogen bond acceptors and was utilized as a 3D query in virtual screening using ZINC, Covalent, and in-house databases. A total of 12 hit compounds were chosen on the basis of their critical interactions with the significant amino acids of ITK. The orbital energies such as HOMO and LUMO of the hit compounds were calculated to evaluate the inhibitor's potencies. Further, molecular dynamics simulation demonstrated the stability of ITK upon binding of selected virtual hits. Binding energy using the MMGBSA method showed the potential binding affinity of all the hits with ITK. The research identifies key chemical characteristics with geometric restrictions that lead to ITK inhibition.Communicated by Ramaswamy H. Sarma.

Probing RNA-Small Molecule Interactions Using Biophysical and Computational Approaches

Interest in small molecules that target RNA is flourishing, and the expectation set on them to treat diseases with unmet medical needs is high. However, several challenges remain, including difficulties in selecting suitable tools and establishing workflows for their discovery. In this context, we optimized experimental and computational approaches that were previously employed for the protein targets. Here, we demonstrate that a fluorescence-based assay can be effectively used to screen small molecule libraries for their ability to bind and stabilize an RNA stem-loop. Our screen identified several fluoroquinolones that bind to the target stem-loop. We further probed their interactions with the target using biolayer interferometry, isothermal titration calorimetry (ITC), and nuclear magnetic resonance spectroscopy. The results of these biophysical assays suggest that the fluoroquinolones bind the target in a similar manner. Armed with this knowledge, we built models for the complexes of the fluoroquinolones and the RNA target. Then, we performed fragment molecular orbital (FMO) calculations to dissect the interactions between the fluoroquinolones and the RNA. We found that the binding free energies obtained from the ITC experiments correlated strongly with the interaction energies calculated by FMO. Finally, we designed fluoroquinolone analogues and performed FMO calculations to predict their binding free energies. Taken together, the results of this study support the importance of conducting orthogonal assays in binding confirmation and compound selection and demonstrate the usefulness of FMO calculations in the rational design of RNA-targeted small molecules.

Anisotropic dynamics of an interfacial enzyme active site observed using tethered substrate analogs and ultrafast 2D IR spectroscopy

Enzymes accelerate the rates of biomolecular reactions by many orders of magnitude compared to bulk solution, and it is widely understood that this catalytic effect arises from a combination of polar pre-organization and electrostatic transition state stabilization. A number of recent reports have also implicated ultrafast (femtosecond-picosecond) timescale motions in enzymatic activity. However, complications arising from spatially-distributed disorder, the occurrence of multiple substrate binding modes, and the influence of hydration dynamics on solvent-exposed active sites still confound many experimental studies. Here we use ultrafast two-dimensional infrared (2D IR) spectroscopy and covalently-tethered substrate analogs to examine dynamical properties of the promiscuous Pyrococcus horikoshii ene-reductase (PhENR) active site in two binding configurations mimicking proposed "inactive" and "reactive" Michaelis complexes. Spectral diffusion measurements of aryl-nitrile substrate analogs reveal an end-to-end tradeoff between fast (sub-ps) and slow (>5 ps) motions. Fermi resonant aryl-azide analogs that sense interactions of coupled oscillators are described. Lineshape and quantum beat analyses of these probes reveal characteristics that correlate with aryl-nitrile frequency fluctuation correlation functions parameters, demonstrating that this anisotropy is an intrinsic property of the water-exposed active site, where countervailing gradients of fast dynamics and disorder in the reactant ground state are maintained near the hydration interface. Our results suggest several plausible factors leading to state-selective rate enhancement and promiscuity in PhENR. This study also highlights a strategy to detect perturbations to vibrational modes outside the transparent window of the mid-IR spectrum, which may be extended to other macromolecular systems.

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.

Structural and Functional Characterization of One Unclassified Glutathione S-Transferase in Xenobiotic Adaptation of Leptinotarsa decemlineata

Arthropod Glutathione S-transferases (GSTs) constitute a large family of multifunctional enzymes that are mainly associated with xenobiotic or stress adaptation. GST-mediated xenobiotic adaptation takes place through direct metabolism or sequestration of xenobiotics, and/or indirectly by providing protection against oxidative stress induced by xenobiotic exposure. To date, the roles of GSTs in xenobiotic adaptation in the Colorado potato beetle (CPB), a notorious agricultural pest of plants within Solanaceae, have not been well studied. Here, we functionally expressed and characterized an unclassified-class GST, LdGSTu1. The three-dimensional structure of the LdGSTu1 was solved with a resolution up to 1.8 Å by X-ray crystallography. The signature motif VSDGPPSL was identified in the “G-site”, and it contains the catalytically active residue Ser14. Recombinant LdGSTu1 was used to determine enzyme activity and kinetic parameters using 1-chloro-2, 4-dinitrobenzene (CDNB), GSH, p-nitrophenyl acetate (PNA) as substrates. The enzyme kinetic parameters and enzyme-substrate interaction studies demonstrated that LdGSTu1 could catalyze the conjugation of GSH to both CDNB and PNA, with a higher turnover number for CDNB than PNA. The LdGSTu1 enzyme inhibition assays demonstrated that the enzymatic conjugation of GSH to CDNB was inhibited by multiple pesticides, suggesting a potential function of LdGSTu1 in xenobiotic adaptation.

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.

Interaction analyses of SARS-CoV-2 spike protein based on fragment molecular orbital calculations

At the stage of SARS-CoV-2 infection in human cells, the spike protein consisting of three chains, A, B, and C, with a total of 3300 residues plays a key role, and thus its structural properties and the binding nature of receptor proteins to host human cells or neutralizing antibodies has attracted considerable interest. Here, we report on interaction analyses of the spike protein in both closed (PDB-ID: 6VXX) and open (6VYB) structures, based on large-scale fragment molecular orbital (FMO) calculations at the level of up to the fourth-order Møller–Plesset perturbation with singles, doubles, and quadruples (MP4(SDQ)). Inter-chain interaction energies were evaluated for both structures, and a mutual comparison indicated considerable losses of stabilization energies in the open structure, especially in the receptor binding domain (RBD) of chain-B. The role of charged residues in inter-chain interactions was illuminated as well. By two separate calculations for the RBD complexes with angiotensin-converting enzyme 2 (ACE2) (6M0J) and B38 Fab antibody (7BZ5), it was found that the binding with ACE2 or antibody partially compensated for this stabilization loss of RBD.

Visualized IFIE results seen from chain-B of spike protein.

Intermolecular interaction among Remdesivir, RNA and RNA-dependent RNA polymerase of SARS-CoV-2 analyzed by fragment molecular orbital calculation

COVID-19, a disease caused by a new strain of coronavirus (SARS-CoV-2) originating from Wuhan, China, has now spread around the world, triggering a global pandemic, leaving the public eagerly awaiting the development of a specific medicine and vaccine. In response, aggressive efforts are underway around the world to overcome COVID-19. In this study, referencing the data published on the Protein Data Bank (PDB ID: 7BV2) on April 22, we conducted a detailed analysis of the interaction between the complex structures of the RNA-dependent RNA polymerase (RdRp) of SARS-CoV-2 and Remdesivir, an antiviral drug, from the quantum chemical perspective based on the fragment molecular orbital (FMO) method. In addition to the hydrogen bonding and intra-strand stacking between complementary strands as seen in normal base pairs, Remdesivir bound to the terminus of an primer-RNA strand was further stabilized by diagonal π-π stacking with the -1A' base of the complementary strand and an additional hydrogen bond with an intra-strand base, due to the effect of chemically modified functional group. Moreover, stable OH/π interaction is also formed with Thr687 of the RdRp. We quantitatively revealed the exhaustive interaction within the complex among Remdesivir, template-primer-RNA, RdRp and co-factors, and published the results in the FMODB database.

Fragment Molecular Orbital Based Interaction Analyses on COVID-19 Main Protease - Inhibitor N3 Complex (PDB ID: 6LU7)

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.

Characterizing Interhelical Interactions of G-Protein Coupled Receptors with the Fragment Molecular Orbital Method

G-protein coupled receptors (GPCRs) are the largest superfamily of membrane proteins, regulating almost every aspect of cellular activity and serving as key targets for drug discovery. We have identified an accurate and reliable computational method to characterize the strength and chemical nature of the interhelical interactions between the residues of transmembrane (TM) domains during different receptor activation states, something that cannot be characterized solely by visual inspection of structural information. Using the fragment molecular orbital (FMO) quantum mechanics method to analyze 35 crystal structures representing different branches of the class A GPCR family, we have identified 69 topologically equivalent TM residues that form a consensus network of 51 inter-TM interactions, providing novel results that are consistent with and help to rationalize experimental data. This discovery establishes a comprehensive picture of how defined molecular forces govern specific interhelical interactions which, in turn, support the structural stability, ligand binding, and activation of GPCRs.

Conformational Searching with Quantum Mechanics

Estimating the range of three-dimensional structures (conformations) that are available to a molecule is a key component of computer-aided drug design. Quantum mechanical simulation offers improved accuracy over forcefield methods, but at a high computational cost. The question is whether this increased cost can be justified in a context in which high-throughput analysis of large numbers of molecules is often key. This chapter discusses the application of quantum mechanics to conformational searching, with a focus on three key challenges: (1) the generation of ensembles that include a good approximation to a molecule's bioactive conformation at as prominent a ranking as possible; (2) rational analysis and modification of a pre-established bioactive conformation in terms of its energetics; and (3) approximation of real solution-phase conformational ensembles in tandem with NMR data. The impact of QM on the high-throughput application (1) is debatable, meaning that for the moment its primary application is still lower-throughput applications such as (2) and (3). The optimal choice of QM method is also discussed. Rigorous benchmarking suggests that DFT methods are only acceptable when used with large basis sets, but a trickle of papers continue to obtain useful results with relatively low-cost methods, leading to a dilemma that the literature has yet to fully resolve.

Characterising GPCR-ligand interactions using a fragment molecular orbital-based approach

There has been fantastic progress in solving GPCR crystal structures. However, the ability of X-ray crystallography to guide the drug discovery process for GPCR targets is limited by the availability of accurate tools to explore receptor-ligand interactions. Visual inspection and molecular mechanics approaches cannot explain the full complexity of molecular interactions. Quantum mechanical approaches (QM) are often too computationally expensive, but the fragment molecular orbital (FMO) method offers an excellent solution that combines accuracy, speed and the ability to reveal key interactions that would otherwise be hard to detect. Integration of GPCR crystallography or homology modelling with FMO reveals atomistic details of the individual contributions of each residue and water molecule towards ligand binding, including an analysis of their chemical nature.