Enhancement of energy decomposition analysis in fragment molecular orbital calculations

Energy decomposition analysis is one of the most attractive features of fragment molecular orbital (FMO) calculations from the point of view of practical applications. Here we report some enhancements for PIEDA in the ABINIT-MP program. One is a separation of the dispersion-type stabilization from the electron correlation energy, traditionally referred to as the "dispersion interaction" (DI). Another is an alternative evaluation of the electrostatic (ES) interaction using the restrained electrostatic potential (RESP) charges. The GA:CT stacked base pair and the Trp-Cage miniprotein were used as illustrative examples.

Differential molecular mechanisms of substrate recognition by selenium methyltransferases, INMT and TPMT, in selenium detoxification and excretion

It is known that the recommended dietary allowance of selenium (Se) is dangerously close to its tolerable upper intake level. Se is detoxified and excreted in urine as trimethylselenonium ion (TMSe) when the amount ingested exceeds the nutritional level. Recently, we demonstrated that the production of TMSe requires two methyltransferases: thiopurine S-methyltransferase (TPMT) and indolethylamine N-methyltransferase (INMT). In this study, we investigated the substrate recognition mechanisms of INMT and TPMT in the Se-methylation reaction. Examination of the Se-methyltransferase activities of two paralogs of INMT, namely, nicotinamide N-methyltransferase and phenylethanolamine N-methyltransferase, revealed that only INMT exhibited Se-methyltransferase activity. Consistently, molecular dynamics simulations demonstrated that dimethylselenide was preferentially associated with the active center of INMT. Using the fragment molecular orbital method, we identified hydrophobic residues involved in the binding of dimethylselenide to the active center of INMT. The INMT-L164R mutation resulted in a deficiency in Se- and N-methyltransferase activities. Similarly, TPMT-R152, which occupies the same position as INMT-L164, played a crucial role in the Se-methyltransferase activity of TPMT. Our findings suggest that TPMT recognizes negatively charged substrates, whereas INMT recognizes electrically neutral substrates in the hydrophobic active center embedded within the protein. These observations explain the sequential requirement of the two methyltransferases in producing TMSe.

Binding Free Energy Calculation Based on the Fragment Molecular Orbital Method and Its Application in Designing Novel SHP-2 Allosteric Inhibitors

The accurate prediction of binding free energy is a major challenge in structure-based drug design. Quantum mechanics (QM)-based approaches show promising potential in predicting ligand-protein binding affinity by accurately describing the behavior and structure of electrons. However, traditional QM calculations face computational limitations, hindering their practical application in drug design. Nevertheless, the fragment molecular orbital (FMO) method has gained widespread application in drug design due to its ability to reduce computational costs and achieve efficient ab initio QM calculations. Although the FMO method has demonstrated its reliability in calculating the gas phase potential energy, the binding of proteins and ligands also involves other contributing energy terms, such as solvent effects, the 'deformation energy' of a ligand's bioactive conformations, and entropy. Particularly in cases involving ionized fragments, the calculation of solvation free energy becomes particularly crucial. We conducted an evaluation of some previously reported implicit solvent methods on the same data set to assess their potential for improving the performance of the FMO method. Herein, we develop a new QM-based binding free energy calculation method called FMOScore, which enhances the performance of the FMO method. The FMOScore method incorporates linear fitting of various terms, including gas-phase potential energy, deformation energy, and solvation free energy. Compared to other widely used traditional prediction methods such as FEP+, MM/PBSA, MM/GBSA, and Autodock vina, FMOScore showed good performance in prediction accuracies. By constructing a retrospective case study, it was observed that incorporating calculations for solvation free energy and deformation energy can further enhance the precision of FMO predictions for binding affinity. Furthermore, using FMOScore-guided lead optimization against Src homology-2-containing protein tyrosine phosphatase 2 (SHP-2), we discovered a novel and potent allosteric SHP-2 inhibitor (compound 8).

Leveraging the Fragment Molecular Orbital Method to Explore the PLK1 Kinase Binding Site and Polo-Box Domain for Potent Small-Molecule Drug Design

Polo-like kinase 1 (PLK1) plays a pivotal role in cell division regulation and emerges as a promising therapeutic target for cancer treatment. Consequently, the development of small-molecule inhibitors targeting PLK1 has become a focal point in contemporary research. The adenosine triphosphate (ATP)-binding site and the polo-box domain in PLK1 present crucial interaction sites for these inhibitors, aiming to disrupt the protein's function. However, designing potent and selective small-molecule inhibitors can be challenging, requiring a deep understanding of protein-ligand interaction mechanisms at these binding sites. In this context, our study leverages the fragment molecular orbital (FMO) method to explore these site-specific interactions in depth. Using the FMO approach, we used the FMO method to elucidate the molecular mechanisms of small-molecule drugs binding to these sites to design PLK1 inhibitors that are both potent and selective. Our investigation further entailed a comparative analysis of various PLK1 inhibitors, each characterized by distinct structural attributes, helping us gain a better understanding of the relationship between molecular structure and biological activity. The FMO method was particularly effective in identifying key binding features and predicting binding modes for small-molecule ligands. Our research also highlighted specific "hot spot" residues that played a critical role in the selective and robust binding of PLK1. These findings provide valuable insights that can be used to design new and effective PLK1 inhibitors, which can have significant implications for developing anticancer therapeutics.

Protein-ligand binding affinity prediction of cyclin-dependent kinase-2 inhibitors by dynamically averaged fragment molecular orbital-based interaction energy

Fragment molecular orbital (FMO) method is a powerful computational tool for structure-based drug design, in which protein-ligand interactions can be described by the inter-fragment interaction energy (IFIE) and its pair interaction energy decomposition analysis (PIEDA). Here, we introduced a dynamically averaged (DA) FMO-based approach in which molecular dynamics simulations were used to generate multiple protein-ligand complex structures for FMO calculations. To assess this approach, we examined the correlation between the experimental binding free energies and DA-IFIEs of six CDK2 inhibitors whose net charges are zero. The correlation between the experimental binding free energies and snapshot IFIEs for X-ray crystal structures was R²  = 0.75. Using the DA-IFIEs, the correlation significantly improved to 0.99. When an additional CDK2 inhibitor with net charge of -1 was added, the DA FMO-based scheme with the dispersion energies still achieved R²  = 0.99, whereas R² decreased to 0.32 employing all the energy terms of PIEDA.

Sulfated Hyaluronan Binds to Heparanase and Blocks Its Enzymatic and Cellular Actions in Carcinoma Cells

We examined whether sulfated hyaluronan exerts inhibitory effects on enzymatic and biological actions of heparanase, a sole endo-beta-glucuronidase implicated in cancer malignancy and inflammation. Degradation of heparan sulfate by human and mouse heparanase was inhibited by sulfated hyaluronan. In particular, high-sulfated hyaluronan modified with approximately 2.5 sulfate groups per disaccharide unit effectively inhibited the enzymatic activity at a lower concentration than heparin. Human and mouse heparanase bound to immobilized sulfated hyaluronan. Invasion of heparanase-positive colon-26 cells and 4T1 cells under 3D culture conditions was significantly suppressed in the presence of high-sulfated hyaluronan. Heparanase-induced release of CCL2 from colon-26 cells was suppressed in the presence of sulfated hyaluronan via blocking of cell surface binding and subsequent intracellular NF-κB-dependent signaling. The inhibitory effect of sulfated hyaluronan is likely due to competitive binding to the heparanase molecule, which antagonizes the heparanase-substrate interaction. Fragment molecular orbital calculation revealed a strong binding of sulfated hyaluronan tetrasaccharide to the heparanase molecule based on electrostatic interactions, particularly characterized by interactions of (−1)- and (−2)-positioned sulfated sugar residues with basic amino acid residues composing the heparin-binding domain-1 of heparanase. These results propose a relevance for sulfated hyaluronan in the blocking of heparanase-mediated enzymatic and cellular actions.

Computational approach to elucidate the formation and stabilization mechanism of amorphous formulation using molecular dynamics simulation and fragment molecular orbital calculation

α-Glycosyl rutin (Rutin-G) consists of a flavonol skeleton and sugar groups and is a promising additive for amorphous formulations. In our previous study, experimental approaches suggested an interaction between the model drug carbamazepine (CBZ) and flavonol skeleton of Rutin-G that stabilizes amorphous formulations. In the present study, the formation and stabilization mechanisms of CBZ/Rutin-G amorphous formulation were investigated using a computational approach. The CBZ/Rutin-G amorphous formulation was obtained via molecular dynamics (MD) simulation, which mimicked the melt-quenching method. Root mean square deviation analysis revealed that the translational motion of CBZ during the cooling process was suppressed by adding Rutin-G. Monitoring the atomic distance during the cooling process revealed that hydrogen bonds via carboxamide oxygen of CBZ with hydroxyl hydrogen of Rutin-G were preferentially formed with flavonol skeletons than sugar groups. The simulated amorphous formulation was then calculated using fragment molecular orbital (FMO) method. The quantitative evaluation of multiple interactions revealed that the hydrogen bond energy was higher in CBZ-sugar groups than in CBZ-flavonol skeleton, while the π-type of interaction energy was higher in CBZ-flavonol skeleton than in CBZ-sugar groups. The computational approach combining MD simulation and FMO calculation provides information on various interactions that are difficult to detect using experimental approaches, which helps understand the formation and stabilization mechanism of amorphous formulations.

Hot Spot Analysis of YAP-TEAD Protein-Protein Interaction Using the Fragment Molecular Orbital Method and Its Application for Inhibitor Discovery

The Hippo pathway is an important signaling pathway modulating growth control and cancer cell proliferation. Dysregulation of the Hippo pathway is a common feature of several types of cancer cells. The modulation of the interaction between yes-associated protein (YAP) and transcriptional enhancer associated domain (TEAD) in the Hippo pathway is considered an attractive target for cancer therapeutic development, although the inhibition of PPI is a challenging task. In order to investigate the hot spots of the YAP and TEAD1 interacting complex, an ab initio Fragment Molecular Orbital (FMO) method was introduced. With the hot spots, pharmacophores for the inhibitor design were constructed, then virtual screening was performed to an in-house library. Next, we performed molecular docking simulations and FMO calculations for screening results to study the binding modes and affinities between PPI inhibitors and TEAD1. As a result of the virtual screening, three compounds were selected as virtual hit compounds. In order to confirm their biological activities, cellular (luciferase activity, proximity ligation assay and wound healing assay in A375 cells, qRT-PCR in HEK 293T cells) and biophysical assays (surface plasmon resonance assays) were performed. Based on the findings of the study, we propose a novel PPI inhibitor BY03 and demonstrate a profitable strategy to analyze YAP-TEAD PPI and discover novel PPI inhibitors.

Revisiting the Charge-Transfer States at Pentacene/C60 Interfaces with the GW/Bethe-Salpeter Equation Approach

Molecular orientations and interfacial morphologies have critical effects on the electronic states of donor/acceptor interfaces and thus on the performance of organic photovoltaic devices. In this study, we explore the energy levels and charge-transfer states at the organic donor/acceptor interfaces on the basis of the fragment-based GW and Bethe-Salpeter equation approach. The face-on and edge-on orientations of pentacene/C60 bilayer heterojunctions have employed as model systems. GW+Bethe-Salpeter equation calculations were performed for the local interface structures in the face-on and edge-on bilayer heterojunctions, which contain approximately 2000 atoms. Calculated energy levels and charge-transfer state absorption spectra are in reasonable agreements with those obtained from experimental measurements. We found that the dependence of the energy levels on interfacial morphology is predominantly determined by the electrostatic contribution of polarization energy, while the effects of induction contribution in the edge-on interface are similar to those in the face-on. Moreover, the delocalized charge-transfer states contribute to the main absorption peak in the edge-on interface, while the face-on interface features relatively localized charge-transfer states in the main absorption peak. The impact of the interfacial morphologies on the polarization and charge delocalization effects is analyzed in detail.

Fragmentation at sp2 carbon atoms in fragment molecular orbital method

In the fragment molecular orbital (FMO) method, a given molecular system is usually fragmented at sp³ carbon atoms. However, fragmentation at different sites sometimes becomes necessary. Hence, we propose fragmentation at sp² carbon atoms in the FMO method. Projection operators are constructed using sp² local orbitals. To maintain practical accuracy, it is essential to consider the three-body effect. In order to suppress the corresponding increase of computational cost, we propose approximate models considering local trimers. Numerical verification shows that the present models are as accurate as or better than the standard FMO2 method in total energy with fragmentation at sp³ carbon atoms.

RI-MP3 calculations of biomolecules based on the fragment molecular orbital method

In this study, the third-order Møller-Plesset perturbation (MP3) theory using the resolution of the identity (RI) approximation was combined with the fragment molecular orbital (FMO) method to efficiently calculate a high-order electron correlation energy of biomolecular systems. We developed a new algorithm for the RI-MP3 calculation, which can be used with the FMO scheme. After test calculations using a small molecule, the FMO-RI-MP3 calculations were performed for two biomolecular systems comprising a protein and a ligand. The computational cost of these calculations was only around 5 and 4 times higher than those of the FMO-RHF calculations. The error associated with the RI approximation was around 2.0% of the third-order correlation contribution to the total energy. However, the RI approximation error in the interaction energy between the protein and ligand molecule was insignificantly small, which reflected the negligible error in the inter fragment interaction energy. © 2018 Wiley Periodicals, Inc.

A study of CDR3 loop dynamics reveals distinct mechanisms of peptide recognition by T-cell receptors exhibiting different levels of cross-reactivity

T-cell receptors (TCRs) can productively interact with many different peptides bound within the MHC binding groove. This property varies with the level of cross-reactivity of TCRs; some TCRs are particularly hyper cross-reactive while others exhibit greater specificity. To elucidate the mechanism behind these differences, we studied five TCRs in complex with the same class II MHC (1A^b )-peptide (3K), that are known to exhibit different levels of cross-reactivity. Although these complexes have similar binding affinities, the interface areas between the TCR and the peptide-MHC (pMHC) differ significantly. We investigated static and dynamic structural features of the TCR-pMHC complexes and of TCRs in a free state, as well as the relationship between binding affinity and interface area. It was found that the TCRs known to exhibit lower levels of cross-reactivity bound to pMHC using an induced-fitting mechanism, forming large and tight interfaces rich in specific hydrogen bonds. In contrast, TCRs known to exhibit high levels of cross-reactivity used a more rigid binding mechanism where non-specific π-interactions involving the bulky Trp residue in CDR3β dominated. As entropy loss upon binding in these highly degenerate and rigid TCRs is smaller than that in less degenerate TCRs, they can better tolerate changes in residues distal from the major contacts with MHC-bound peptide. Hence, our dynamics study revealed that differences in the peptide recognition mechanisms by TCRs appear to correlate with the levels of T-cell cross-reactivity.

RI-MP2 Gradient Calculation of Large Molecules Using the Fragment Molecular Orbital Method

The second-order Møller-Plesset perturbation theory (MP2) gradient using resolution of the identity approximation (RI-MP2 gradient) was combined with the fragment molecular orbital (FMO) method to evaluate the gradient including electron correlation for large molecules. In this study, we adopted a direct implementation of the RI-MP2 gradient, in which a characteristic feature of the FMO scheme was utilized. Test calculations with a small peptide presented a computational advantage of the RI-MP2 gradient over the canonical MP2 gradient. In addition, it was shown that the error of the RI-MP2 gradient, caused by RI approximation, was negligible. As an illustrative example, we performed gradient calculations for two biomolecules-a prion protein with GN8 and a human immunodeficiency virus type 1 (HIV1) protease with lopinavir (LPV). These calculations demonstrated that the gradient including the correlation effect could be evaluated with only about twice the computational effort of the Hartree-Fock (HF) gradient.