Fragment Molecular Orbital Calculations Research Articles

Quantum chemical calculation dataset for representative protein folds by the fragment molecular orbital method

The function of a biomacromolecule is not only determined by its three-dimensional structure but also by its electronic state. Quantum chemical calculations are promising non-empirical methods available for determining the electronic state of a given structure. In this study, we used the fragment molecular orbital (FMO) method, which applies to biopolymers such as proteins, to provide physicochemical property values on representative structures in the SCOP2 database of protein families, a subset of the Protein Data Bank. Our dataset was constructed by over 5,000 protein structures, including over 200 million inter-fragment interaction energies (IFIEs) and their energy components obtained by pair interaction energy decomposition analysis (PIEDA) using FMO-MP2/6-31 G*. Moreover, three basis sets, 6-31 G*, 6-31 G**, and cc-pVDZ, were used for the FMO calculations of each structure, making it possible to compare the energies obtained with different basis functions for the same fragment pair. The total data size is approximately 6.7 GB. Our dataset will be useful for functional analyses and machine learning based on the physicochemical property values of proteins.

New Catalytic Residues and Catalytic Mechanism of the RNase T1 Family.

The ribonuclease T1 family, including RNase Po1 secreted by Pleurotus ostreatus, exhibits antitumor activity. Here, we resolved the Po1/guanosine-3'-monophosphate complex (3'GMP) structure at 1.75 Å. Structure comparison and fragment molecular orbital (FMO) calculation between the apo form and the Po1/3'GMP complex identified Phe38, Phe40, and Glu42 as the key binding residues. Two types of the RNase/3'GMP complex in RNasePo1 and RNase T1 were homologous to Po1, and FMO calculations elucidated that the biprotonated histidine on the β3 sheet (His36) on the β3 sheet and deprotonated Glu54 on the β4 sheet were advantageous to RNase activity. Moreover, tyrosine (Tyr34) on the β3 sheet was elucidated as a crucial catalytic residues. Mutation of Tyr34 with phenylalanine decreased RNase activity and diminished antitumor efficacy compared to that in the wild type. This suggests the importance of RNase activity in antitumor mechanisms.

Size-consistency and orbital-invariance issues revealed by VQE-UCCSD calculations with the FMO scheme.

The fragment molecular orbital (FMO) scheme is one of the popular fragmentation-based methods and has the potential advantage of making the circuit shallow for quantum chemical calculations on quantum computers. In this study, we used a GPU-accelerated quantum simulator (cuQuantum) to perform the electron correlation part of the FMO calculation as unitary coupled-cluster singles and doubles (UCCSD) with the variational quantum eigensolver (VQE) for hydrogen-bonded (FH) and (FH) -H O systems with the STO-3G basis set. VQE-UCCSD calculations were performed using both canonical and localized MO sets, and the results were examined from the point of view of size-consistency and orbital-invariance affected by the Trotter error. It was found that the use of localized MO leads to better results, especially for (FH) -H O. The GPU acceleration was substantial for the simulations with larger numbers of qubits, and was about a factor of 6.7-7.7 for 18 qubit systems.

Systematic improvement of empirical energy functions in the era of machine learning.

The impact of targeted replacement of individual terms in empirical force fields is quantitatively assessed for pure water, dichloromethane (CH Cl ), and solvated K and Cl ions. For the electrostatic interactions, point charges (PCs) and machine learning (ML)-based minimally distributed charges (MDCM) fitted to the molecular electrostatic potential are evaluated together with electrostatics based on the Coulomb integral. The impact of explicitly including second-order terms is investigated by adding a fragment molecular orbital (FMO)-derived polarization energy to an existing force field, in this case CHARMM. It is demonstrated that anisotropic electrostatics reduce the RMSE for water (by 1.4 kcal/mol), CH Cl (by 0.8 kcal/mol) and for solvated Cl clusters (by 0.4 kcal/mol). An additional polarization term can be neglected for CH Cl but further improves the models for pure water (by 1.0 kcal/mol) and hydrated Cl (by 0.4 kcal/mol), and is key for solvated K , reducing the RMSE by 2.3 kcal/mol. A 12-6 Lennard-Jones functional form performs satisfactorily with PC and MDCM electrostatics, but is not appropriate for descriptions that account for the electrostatic penetration energy. The importance of many-body contributions is assessed by comparing a strictly 2-body approach with self-consistent reference data. Two-body interactions suffice for CH Cl whereas water and solvated K and Cl ions require explicit many-body corrections. Finally, a many-body-corrected dimer potential energy surface exceeds the accuracy attained using a conventional empirical force field, potentially reaching that of an FMO calculation. The present work systematically quantifies which terms improve the performance of an existing force field and what reference data to use for parametrizing these terms in a tractable fashion for ML fitting of pure and heterogeneous systems.

DPD simulation to reproduce lipid membrane microdomains based on fragment molecular orbital calculations

Lipid domains play a critical role in signal transduction and transport across cell membranes. The formation of domains in “HLC” ternary lipid bilayers composed of high transition temperature (high-Tm) lipids, low-Tm lipids, and cholesterol (Chol) has been extensively studied as a raft-like system. Recently, experiments were performed to control the formation of submicron domains in LLC lipid bilayers containing low-Tm phosphatidylethanolamine (PE), low-Tm phosphatidylcholine (PC), and Chol by manipulating the presence or absence of Chol. The formation of microdomains in this LLC mixture was replicated by dissipative particle dynamics simulation. The results show that domain formation can be replicated.

Prediction of Binding Pose and Affinity of Nelfinavir, a SARS-CoV-2 Main Protease Repositioned Drug, by Combining Docking, Molecular Dynamics, and Fragment Molecular Orbital Calculations.

A novel in silico drug design procedure is described targeting the Main protease (Mpro) of the SARS-CoV-2 virus. The procedure combines molecular docking, molecular dynamics (MD), and fragment molecular orbital (FMO) calculations. The binding structure and properties of Mpro were predicted for Nelfinavir (NFV), which had been identified as a candidate compound through drug repositioning, targeting Mpro. Several poses of the Mpro and NFV complexes were generated by docking, from which four docking poses were selected by scoring with FMO energy. Then, each pose was subjected to MD simulation, 100 snapshot structures were sampled from each of the generated MD trajectories, and the structures were evaluated by FMO calculations to rank the pose based on binding energy. Several residues were found to be important in ligand recognition, including Glu47, Asp48, Glu166, Asp187, and Gln189, all of which interacted strongly with NFV. Asn142 is presumably regarded to form hydrogen bonds or CH/π interaction with NFV; however, in the present calculation, their interactions were transient. Moreover, the tert-butyl group of NFV had no interaction with Mpro. Identifying such strong and weak interactions provides candidates for maintaining and substituting ligand functional groups and important suggestions for drug discovery using drug repositioning. Besides the interaction between NFV and the amino acid residues of Mpro, the desolvation effect of the binding pocket also affected the ranking order. A similar procedure of drug design was applied to Lopinavir, and the calculated interaction energy and experimental inhibitory activity value trends were consistent. Our approach provides a new guideline for structure-based drug design starting from a candidate compound whose complex crystal structure has not been obtained.

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.

Synthesis and biological evaluation of 2′-hydroxychalcone derivatives as AMPK activators

A series of 2′-hydroxychalcone derivatives with various substituents on B-ring were synthesized and evaluated for AMP-activated protein kinase (AMPK) activation activity in podocyte cells. The results displayed that hydroxy, methoxy and methylenedioxy groups on B-ring could enhance the activitiy better than O-saturated alkyl, O-unsaturated alkyl or other alkoxy groups. Compounds 27 and 29 possess the highest fold change of 2.48 and 2.73, respectively, which were higher than those of reference compound (8) (1.28) and metformin (1.88). Compounds 27 and 29 were then subjected to a concentration–response study to obtain the EC50 values of 2.0 and 4.8 µM, respectively and MTT assays also showed that cell viability was not influenced by the exposure of podocytes to compounds 27 and 29 at concentrations up to 50 μM. In addition, compound 27 was proved to activate AMPK via calcium/calmodulin-dependent protein kinase kinase β (CaMKKβ)-dependent pathway without affecting intracellular calcium levels. The computational study showed that the potent compounds exhibited stronger ligand-binding strength to CaMKKβ, particularly compounds 27 (−8.4 kcal/mol) and 29 (−8.0 kcal/mol), compared to compound 8 (-7.5 kcal/mol). Fragment molecular orbital (FMO) calculation demonstrated that compound 27 was superior to compound 29 due to the presence of methyl group, which amplified the binding by hydrophobic interactions. Therefore, compound 27 would represent a promising AMPK activator for further investigation of the treatment of diabetes and diabetic nephropathy.

In Silico Design of Natural Inhibitors of ApoE4 from the Plant Moringa oleifera: Molecular Docking and Ab Initio Fragment Molecular Orbital Calculations.

Alzheimer's disease (AD) is a neurological disease, and its signs and symptoms appear slowly over time. Although current Alzheimer's disease treatments can alleviate symptoms, they cannot prevent the disease from progressing. To accurately diagnose and treat Alzheimer's disease, it is therefore necessary to establish effective methods for diagnosis. Apolipoprotein E4 (ApoE4), the most frequent genetic risk factor for AD, is expressed in more than half of patients with AD, making it an attractive target for AD therapy. We used molecular docking simulations, classical molecular mechanics optimizations, and ab initio fragment molecular orbital (FMO) calculations to investigate the specific interactions between ApoE4 and the naturally occurring compounds found in the plant Moringa Oleifera. According to the FMO calculations, quercetin had the highest binding affinity to ApoE4 among the sixteen compounds because its hydroxyl groups generated strong hydrogen bonds with the ApoE4 residues Trp11, Asp12, Arg15, and Asp130. As a result, we proposed various quercetin derivatives by introducing a hydroxyl group into quercetin and studied their ApoE4 binding properties. The FMO data clearly showed that adding a hydroxyl group to quercetin improved its binding capacity to ApoE4. Furthermore, ApoE4 Trp11, Asp12, Arg15, and Asp130 residues were discovered to be required for significant interactions between ApoE4 and quercetin derivatives. They had a higher ApoE4 binding affinity than our previously proposed epicatechin derivatives. Accordingly, the current results evaluated using the ab initio FMO method will be useful for designing potent ApoE4 inhibitors that can be used as a candidate agent for AD treatment.

Elucidating specific interactions for designing novel pyrrolamide derivatives as potential GyrB inhibitors based on ab initio fragment molecular orbital calculations

Tuberculosis (TB), the second leading infectious killer, causes serious public health problems worldwide. To develop novel anti-TB agents, many biochemical studies have targeted the subunit B of DNA gyrase (GyrB), which captures a second DNA segment and responses for ATP hydrolysis. Here, we investigated specific interactions between GyrB residues and existing pyrrolamide derivatives at an electronic level using ab initio fragment molecular orbital (FMO) calculations and designed potent inhibitors against GyrB. The evaluated binding affinities between GyrB and pyrrolamides were confirmed to be consistent with the IC50 values obtained from previous experiments. Thus, we employed the most potent pyrrolamide (compound 1) as a lead compound and proposed novel pyrrolamide derivatives. The specific interactions between GyrB and these derivatives were investigated using molecular mechanic optimizations and FMO calculations. The results revealed that our proposed derivatives had strong hydrogen bonds with Asp79 and Arg141 and exhibited electrostatic interactions with Glu56 and Ile84 of GyrB. In addition, the binding affinity between GyrB and compound 1 was enhanced significantly by the replacement at the R3 site of compound 1. The present results may provide structural concepts for the rational design of potent GyrB inhibitors as anti-TB agents.Communicated by Ramaswamy H. Sarma

Development of reverse mapping system bridging dissipative particle dynamics and fragment molecular orbital calculation

We have developed a reverse mapping system to convert mesoscale structures generated by dissipative particle dynamics (DPD) simulations into nanoscale structures. This system is called DSRMS (DPD-based structure reverse mapping system) and is controlled by Python3 scripts using OCTA’s COGNAC program for DPD and molecular dynamics. The restored structures can be subjected to fragment molecular orbital calculations using the ABINIT-MP program for detailed nanoscale interaction analysis. Polyelectrolyte and lipid membranes have been used as illustrative example.

Moving out of CF3 -Land: Synthesis, Receptor Affinity, and in silico Studies of NK1 Receptor Ligands Containing a Pentafluorosulfanyl (SF5 ) Group.

The NK1 receptor (NK1R) is a molecular target for both approved and experimental drugs intended for a variety of conditions, including emesis, pain, and cancers. While contemplating modifications to the typical NK1R pharmacophore, we wondered whether the CF3 groups common for many NK1R ligands, could be replaced with some other moiety. Our attention was drawn by the SF5 group, and so we designed, synthesized, and tested ten novel SF5 -containing compounds for NK1R affinity. All analogues exhibit detectable NK1R binding, with the best of them, compound 5 a, (3-bromo-5-(pentafluoro-λ6 -sulfanyl)benzyl acetyl-L-tryptophanate) binding only slightly worse (IC50 =34.3 nM) than the approved NK1R-targeting drug, aprepitant (IC50 =27.7 nM). Molecular docking provided structural explanation of SAR. According to our analysis, the SF5 group in our compounds occupies a position similar to that of one of the CF3 groups of aprepitant as found in the crystal structure. Additionally, we checked whether the docking scoring function or energies derived from Fragment Molecular Orbital quantum chemical calculations may be helpful in explaining and predicting the experimental receptor affinities for our analogues. Both these methods produce moderately good results. Overall, this is the first demonstration of the utility of the SF5 group in the design of NK1R ligands.

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.

The 8-bromobaicalein alleviated chikungunya-induced musculoskeletal inflammation and reduced the viral load in healthy adult mice

ABSTRACT Chikungunya virus is a re-emerging arbovirus that has caused epidemic outbreaks in recent decades. Patients in older age groups with high viral load and severe immunologic response during acute infection are likely to develop chronic arthritis and severe joint pain. Currently, no antiviral drug is available. Previous studies suggested that a flavone derivative, 8-bromobaicalein, was a potential dengue and Zika replication inhibitor in a cell-based system targeting flaviviral polymerase. Here we characterized that 8-bromobaicalein inhibited chikungunya virus replication with EC50 of 0.49 ± 0.11 µM in Vero cells. The molecular target predicted at viral nsP1 methyltransferase using molecular binding and fragment molecular orbital calculation. Additionally, oral administration of 250 mg/kg twice daily treatment alleviated chikungunya-induced musculoskeletal inflammation and reduced viral load in healthy adult mice. Pharmacokinetic analysis indicated that the 250 mg/kg administration maintained the compound level above EC99.9 for 12 h. Therefore, 8-bromobaicalein should be a potential candidate for further development as a pan-arboviral drug.

Structure-based lead optimization to improve potency and selectivity of a novel tetrahydroimidazo[1,2-a]pyridine-5-carboxylic acid series of heparanase-1 inhibitor

Heparanase-1 (HPSE1) is an endo-β-d-glucuronidase that is the only mammalian enzyme known to cleave heparan sulfate (HS) of heparan sulfate proteoglycans (HSPG), a key component of the glycocalyx layer of the vascular endothelium matrix. Inhibition of HPSE1 has therapeutic potential for cancer and proteinuric kidney diseases. We previously reported that 2 showed a moderate potency as an HPSE1 inhibitor and an issue of selectivity against exo-β-d-glucuronidase (GUSβ) and glucocerebrosidase (GBA) remained. A structure-based lead optimization of 2 using X-ray co-crystal structure analysis and fragment molecular orbital calculation resulted in 4e, which showed a more than 7-fold increase in HPSE1 inhibitory activity. The subsequent introduction of a methyl group into the 6-hydroxy group of 4e resulted in 18 with reduced inhibitory activities against GUSβ and GBA while maintaining the inhibitory activity against HPSE1. The inhibitory activities of 18 against serum HPSE1 in mice were significant and lasted for 4 h at doses of 3, 30, and 100 mg/kg. Compound 18 could be a novel lead compound for HPSE1 inhibitors with improved inhibitory activity against HPSE1 and increased HPSE1 selectivity over GUSβ and GBA.

Identification of Promising Sulfonamide Chalcones as Inhibitors of SARS-CoV-2 3CLpro through Structure-Based Virtual Screening and Experimental Approaches.

3CLpro is a viable target for developing antiviral therapies against the coronavirus. With the urgent need to find new possible inhibitors, a structure-based virtual screening approach was developed. This study recognized 75 pharmacologically bioactive compounds from our in-house library of 1052 natural product-based compounds that satisfied drug-likeness criteria and exhibited good bioavailability and membrane permeability. Among these compounds, three promising sulfonamide chalcones were identified by combined theoretical and experimental approaches, with SWC423 being the most suitable representative compound due to its competitive inhibition and low cytotoxicity in Vero E6 cells (EC50 = 0.89 ± 0.32 μM; CC50 = 25.54 ± 1.38 μM; SI = 28.70). The binding and stability of SWC423 in the 3CLpro active site were investigated through all-atom molecular dynamics simulation and fragment molecular orbital calculation, indicating its potential as a 3CLpro inhibitor for further SARS-CoV-2 therapeutic research. These findings suggested that inhibiting 3CLpro with a sulfonamide chalcone such as SWC423 may pave the effective way for developing COVID-19 treatments.

Low entropic cost of binding confers high selectivity on an allosteric ERK2 inhibitor

Extracellular signal-regulated kinase 2 (ERK2), a mitogen-activated protein kinase (MAPK), plays an essential role in physiological cellular processes and is a drug target for treating cancers and type 2 diabetes. A previous in silico screening study focusing on an allosteric site that plays a crucial role in substrate anchoring conferred an ERK2 inhibitor (compound 1). In this report, compound 1 was found to show high selectivity toward ERK2 compared with the nearest off-target p38α MAPK, and the crystal structure revealed that compound 1 binds to the allosteric site of ERK2. Fragment molecular orbital calculations based upon this crystal structure provided the structural basis to improve potency of compound 1 derivatives. Further computational studies uncovered that the low entropic cost of binding conferred the high selectivity of compound 1 toward ERK2 over p38α MAPK. These findings demonstrate the feasibility of developing potent and selective ERK2 inhibitors.

Structural and Computational Analyses of the Unique Interactions of Opicapone in the Binding Pocket of Catechol O-Methyltransferase: A Crystallographic Study and Fragment Molecular Orbital Analyses.

A third-generation inhibitor of catechol O-methyltransferase (COMT), opicapone (1), has the 3-nitrocatechol scaffold as do the second-generation inhibitors such as entacapone (2) and tolcapone (3), but only 1 can sustainably inhibit COMT activity making it suitable for a once-daily regimen. These improvements should be attributed to the optimized sidechain moiety (oxidopyridyloxadiazolyl group) of 1 substituted at the 5-position of the 3-nitrocatechol ring. We analyzed the role of the sidechain moiety by solving the crystal structures of COMT/S-adenosylmethionine (SAM)/Mg/1 and COMT/S-adenosylhomocysteine (SAH)/Mg/1 complexes. Fragment molecular orbital (FMO) calculations elucidated that the dispersion interaction between the sidechains of Leu 198 and Met 201 on the β6β7-loop and the oxidopyridine ring of 1 were unique and important in both complexes. In contrast, the catechol binding site made a remarkable difference in the sidechain conformation of Lys 144. The ε-amino group of Lys 144 was outside of the catalytic pocket and was replaced by a water molecule in the COMT/SAH/Mg/1 complex. No nitrocatechol inhibitor has ever been reported to make a complex with COMT and SAH. Thus, the conformational change of Lys 144 found in the COMT/SAH/Mg/1 complex is the first crystallographic evidence that supports the role of Lys 144 as a catalytic base to take out a proton ion from the reaction site to the outside of the enzyme. The fact that 1 generated a complex with SAH and COMT also suggests that 1 could inhibit COMT twofold, as a typical substrate mimic competitive inhibitor and as a product-inhibition enhancer.

Bifurcated Hydrogen Bonds in a Peptide Crystal Unveiled by X-ray Diffraction and Polarized Raman Spectroscopy

The strength of hydrogen bonds in molecular crystals has been investigated through X-ray structural analysis and vibrational spectroscopy, and the interatomic distances and the stretching vibrational frequency involving hydrogen bonds have been known to exhibit a linear relationship. However, these studies are limited to relatively small molecular compounds, such as amino acids. In this study, we employed tetrapeptide as a model system to investigate the hydrogen bond network in the peptide crystal. Single crystal X-ray diffraction revealed that the peptide crystal structure had bifurcated hydrogen bonds. Raman spectroscopy exhibited that the vibrational frequency of amide bonds correlated linearly to their interatomic distance, but the rate of the change was significantly low compared with previous works. Fragment molecular orbital calculations also revealed that bifurcated hydrogen bonds affect the strength of the hydrogen bonds. This study provides a valuable model system for discussing the effects of bifurcated hydrogen bonds in various crystals of peptides.

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.