Application of singular value decomposition to the inter-fragment interaction energy analysis for ligand screening

Effect of Direct and Water-Mediated Interactions on the Identification of Hotspots in Biomolecular Complexes with Multiple Subsystems

Identification of important residues in biochemical complexes is often a crucial step for many problems in molecular biology and biochemistry. A method is proposed to identify hotspots in biomolecular complexes based on a new metric, derived from networks representing molecular subunits (residues, bridging solvent molecules, ligands etc.) connected by interactions. A singular value decomposition of the weighted adjacency matrix is used to construct a scalar rank for each subunit that reflects its importance in the residue interaction network. This metric is called the singular value centrality. In addition, a new formalism is proposed to account for water-mediated interactions in the ranking of residues. Interactions for a residue network can be provided by various computational methods. In this work interactions are obtained from full quantum-mechanical calculations of protein-protein complexes using the fragment molecular orbital method. The ranking results are shown to be in good agreement with earlier computational and experimental studies. The developed method can be used to gain a deeper insight into the role of subunits in complex biomolecular systems.

FMO-guided design of darunavir analogs as HIV-1 protease inhibitors

The prevalence of HIV-1 infection continues to pose a significant global public health issue, highlighting the need for antiretroviral drugs that target viral proteins to reduce viral replication. One such target is HIV-1 protease (PR), responsible for cleaving viral polyproteins, leading to the maturation of viral proteins. While darunavir (DRV) is a potent HIV-1 PR inhibitor, drug resistance can arise due to mutations in HIV-1 PR. To address this issue, we developed a novel approach using the fragment molecular orbital (FMO) method and structure-based drug design to create DRV analogs. Using combinatorial programming, we generated novel analogs freely accessible via an on-the-cloud mode implemented in Google Colab, Combined Analog generator Tool (CAT). The designed analogs underwent cascade screening through molecular docking with HIV-1 PR wild-type and major mutations at the active site. Molecular dynamics (MD) simulations confirmed the assess ligand binding and susceptibility of screened designed analogs. Our findings indicate that the three designed analogs guided by FMO, 19-0-14-3, 19-8-10-0, and 19-8-14-3, are superior to DRV and have the potential to serve as efficient PR inhibitors. These findings demonstrate the effectiveness of our approach and its potential to be used in further studies for developing new antiretroviral drugs.

System truncation accelerates binding affinity calculations with the fragment molecular orbital method: A benchmark study

The fragment molecular orbital (FMO) method is a fast quantum-mechanical method that divides systems into pieces of fragments and performs ab initio calculations. The system truncation enables further speed improvement. In this article, we systematically study the effects of system truncations on binding affinity calculations obtained with FMO in combination with either the polarizable continuum model (FMO/PCM) or in combination with the Møller-Plesset method (FMO-MP2). We have used five protein complexes with ligands of several charged states. The calculated binding energies of the size variants of the truncated system, including only a restricted number of atoms around the ligand, are compared to the energy obtained from a full system. The result shows that the systems could be truncated to a radius of 8 Å from neutral ligands within an error of 0.7 kcal/mol, and 12 Å from charged ligands within an error of 1.1 kcal/mol for calculating the binding energy in solution.

Protein-Protein Interaction Modelling with the Fragment Molecular Orbital Method

Fragment molecular orbital (FMO) method enables ab initio quantum-chemical calculations for biomolecular systems with high accuracy and moderate computational cost. Through this analysis we can evaluate the inter-fragment interaction energies (IFIEs) that provide useful measures for effective interactions between the fragments representing amino-acid residues and ligand molecules. Here I describe how to prepare the input structures and perform the FMO calculations for protein-protein complex system. In addition to the pre-processing, some useful tools for the post-processing analysis are also illustrated.

Interspecies Comparison of Interaction Energies between Photosynthetic Protein RuBisCO and 2CABP Ligand

Ribulose 1,5-bisphosphate carboxylase/oxygenase (RuBisCO) functions as the initial enzyme in the dark reactions of photosynthesis, catalyzing reactions that extract CO2 from the atmosphere and fix CO2 into organic compounds. RuBisCO is classified into four types (isoforms I–IV) according to sequence-based phylogenetic trees. Given its size, the computational cost of accurate quantum-chemical calculations for functional analysis of RuBisCO is high; however, recent advances in hardware performance and the use of the fragment molecular orbital (FMO) method have enabled the ab initio analyses of RuBisCO. Here, we performed FMO calculations on multiple structural datasets for various complexes with the 2′-carboxylarabinitol 1,5-bisphosphate (2CABP) ligand as a substrate analog and investigated whether phylogenetic relationships based on sequence information are physicochemically relevant as well as whether novel information unobtainable from sequence information can be revealed. We extracted features similar to the phylogenetic relationships found in sequence analysis, and in terms of singular value decomposition, we identified residues that strongly interacted with the ligand and the characteristics of the isoforms for each principal component. These results identified a strong correlation between phylogenetic relationships obtained by sequence analysis and residue interaction energies with the ligand. Notably, some important residues were located far from the ligand, making comparisons among species using only residues proximal to the ligand insufficient.

pyProGA-A PyMOL plugin for protein residue network analysis

The field of protein residue network (PRN) research has brought several useful methods and techniques for structural analysis of proteins and protein complexes. Many of these are ripe and ready to be used by the proteomics community outside of the PRN specialists. In this paper we present software which collects an ensemble of (network) methods tailored towards the analysis of protein-protein interactions (PPI) and/or interactions of proteins with ligands of other type, e.g. nucleic acids, oligosaccharides etc. In parallel, we propose the use of the network differential analysis as a method to identify residues mediating key interactions between proteins. We use a model system, to show that in combination with other, already published methods, also included in pyProGA, it can be used to make such predictions. Such extended repertoire of methods allows to cross-check predictions with other methods as well, as we show here. In addition, the possibility to construct PRN models from various kinds of input is so far a unique asset of our code. One can use structural data as defined in PDB files and/or from data on residue pair interaction energies, either from force-field parameters or fragment molecular orbital (FMO) calculations. pyProGA is a free open-source software available from https://gitlab.com/Vlado_S/pyproga.

Dynamic Cooperativity of Ligand-Residue Interactions Evaluated with the Fragment Molecular Orbital Method

By the splendid advance in computation power realized with the Fugaku supercomputer, it has become possible to perform ab initio fragment molecular orbital (FMO) calculations for thousands of dynamic structures of protein-ligand complexes in a parallel way. We thus carried out electron-correlated FMO calculations for a complex of the 3C-like (3CL) main protease (M^pro) of the new coronavirus (SARS-CoV-2) and its inhibitor N3 incorporating the structural fluctuations sampled by classical molecular dynamics (MD) simulation in hydrated conditions. Along with a statistical evaluation of the interfragment interaction energies (IFIEs) between the N3 ligand and the surrounding amino-acid residues for 1000 dynamic structure samples, in this study we applied a novel approach based on principal component analysis (PCA) and singular value decomposition (SVD) to the analysis of IFIE data in order to extract the dynamically cooperative interactions between the ligand and the residues. We found that the relative importance of each residue is modified via the structural fluctuations and that the ligand is bound in the pharmacophore in a dynamic manner through collective interactions formed by multiple residues, thus providing new insight into structure-based drug discovery.

Interaction of 8-anilinonaphthalene-1-sulfonate with SARS-CoV-2 main protease and its application as a fluorescent probe for inhibitor identification

The 3C-like main protease of SARS-CoV-2 (3CLPro) is responsible for the cleavage of the viral polyprotein. This process is essential for the viral life cycle. Therefore, 3CLPro is a promising target to develop antiviral drugs for COVID-19 prevention and treatment. Traditional enzymatic assays for the identification of 3CLPro inhibitors rely on peptide-based colorimetric or fluorogenic substrates. However, the COVID-19 pandemic has limit or delay access to these substrates, especially for researchers in developing countries attempting to screen natural product libraries. We explored the use of the fluorescent probe 8-anilinonaphthalene-1-sulfonate (ANS) as an alternative assay for inhibitor identification. Fluorescence enhancement upon binding of ANS to 3CLPro was observed, and this interaction was competitive with a peptide substrate. The utility of ANS-based competitive binding assay to identify 3CLPro inhibitors was demonstrated with the flavonoid natural products baicalein and rutin. The molecular nature of ANS and rutin interaction with 3CLPro was explored with molecular modeling. Our results suggested that ANS could be employed in a competitive binding assay to facilitate the identification of novel SARS-CoV-2 antiviral compounds.

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.

Identification of correlated inter-residue interactions in protein complex based on the fragment molecular orbital method

A theoretical scheme to systematically describe correlated (network-like) interactions between molecular fragments is proposed within the framework of the fragment molecular orbital (FMO) method. The method is mathematically based on the singular value decomposition (SVD) of the inter-fragment interaction energy (IFIE) matrix obtained by the FMO calculation, and can be applied to a comprehensive description of protein-protein interactions in the context of molecular recognition. In the present study we apply the proposed method to a complex of measles virus hemagglutinin and human SLAM receptor, thus finding a usefulness for efficiently eliciting the correlated interactions among the amino-acid residues involved in the two proteins. Additionally, collective interaction networks by amino-acid residues important for mutation experiments can be clearly visualized.

Why Are Lopinavir and Ritonavir Effective against the Newly Emerged Coronavirus 2019? Atomistic Insights into the Inhibitory Mechanisms

Since the emergence of a novel coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) was first reported from Wuhan, China, neither a specific vaccine nor an antiviral drug against SARS-CoV-2 has become available. However, a combination of two HIV-1 protease inhibitors, lopinavir and ritonavir, has been found to be effective against SARS-CoV, and both drugs could bind well to the SARS-CoV 3C-like protease (SARS-CoV 3CLpro). In this work, molecular complexation between each inhibitor and SARS-CoV-2 3CLpro was studied using all-atom molecular dynamics simulations, free energy calculations, and pair interaction energy analyses based on MM/PB(GB)SA and FMO-MP2/PCM/6-31G* methods. Both anti-HIV drugs interacted well with the residues at the active site of SARS-CoV-2 3CLpro. Ritonavir showed a somewhat higher number atomic contacts, a somewhat higher binding efficiency, and a somewhat higher number of key binding residues compared to lopinavir, which correspond with the slightly lower water accessibility at the 3CLpro active site. In addition, only ritonavir could interact with the oxyanion hole residues N142 and G143 via the formation of two hydrogen bonds. The interactions in terms of electrostatics, dispersion, and charge transfer played an important role in the drug binding. The obtained results demonstrated how repurposed anti-HIV drugs could be used to combat COVID-19.

Analyzing Interactions with the Fragment Molecular Orbital Method

Basic concepts in the analysis of binding using the fragment molecular orbital method are discussed at length: polarization, desolvation, and interaction. The components in the pair interaction energy decomposition analysis are introduced, and the analysis is illustrated for a water dimer and a protein-ligand complex.

Multiple Virtual Screening Strategies for the Discovery of Novel Compounds Active Against Dengue Virus: A Hit Identification Study

Dengue infection is caused by a mosquito-borne virus, particularly in children, which may even cause death. No effective prevention or therapeutic agents to cure this disease are available up to now. The dengue viral envelope (E) protein was discovered to be a promising target for inhibition in several steps of viral infection. Structure-based virtual screening has become an important technique to identify first hits in a drug screening process, as it is possible to reduce the number of compounds to be assayed, allowing to save resources. In the present study, pharmacophore models were generated using the common hits approach (CHA), starting from trajectories obtained from molecular dynamics (MD) simulations of the E protein complexed with the active inhibitor, flavanone (FN5Y). Subsequently, compounds presented in various drug databases were screened using the LigandScout 4.2 program. The obtained hits were analyzed in more detail by molecular docking, followed by extensive MD simulations of the complexes. The highest-ranked compound from this procedure was then synthesized and tested on its inhibitory efficiency by experimental assays.

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

We describe several procedures for the preprocessing of fragment molecular orbital (FMO) calculations on p38 mitogen-activated protein (MAP) kinase and discuss the influence of the procedures on the protein–ligand interaction energies represented by inter-fragment interaction energies (IFIEs). The correlation between the summation of IFIEs for a ligand and amino acid residues of protein (IFIE-sum) and experimental affinity values (IC₅₀) was poor when considered for the whole set of protein–ligand complexes. To improve the correlation for prediction of ligand binding affinity, we carefully classified data set by the ligand charge, the DFG-loop state (DFG-in/out loop), which is characteristic of kinase, and the scaffold of ligand. The correlation between IFIE-sums and the activity values was examined using the classified data set. As a result, it was confirmed that there was a selected data set that showed good correlation between IFIE-sum and activity value by appropriate classification. In addition, we found that the differences in protonation and hydrogen orientation caused by subtle differences in preprocessing led to a relatively large difference in IFIE values. Further, we also examined the effect of structure optimization with different force fields. It was confirmed that the difference in the force field had no significant effect on IFIE-sum. From the viewpoint of drug design using FMO calculations, various investigations on IFIE-sum in this research, such as those regarding several classifications of data set and the different procedures of structural preparation, would be expected to provide useful knowledge for improvement of prediction ability about the ligand binding affinity.