Machine learning on ligand-residue interaction profiles to significantly improve binding affinity prediction.

Structure-based virtual screenings (SBVSs) play an important role in drug discovery projects. However, it is still a challenge to accurately predict the binding affinity of an arbitrary molecule binds to a drug target and prioritize top ligands from an SBVS. In this study, we developed a novel method, using ligand-residue interaction profiles (IPs) to construct machine learning (ML)-based prediction models, to significantly improve the screening performance in SBVSs. Such a kind of the prediction model is called an IP scoring function (IP-SF). We systematically investigated how to improve the performance of IP-SFs from many perspectives, including the sampling methods before interaction energy calculation and different ML algorithms. Using six drug targets with each having hundreds of known ligands, we conducted a critical evaluation on the developed IP-SFs. The IP-SFs employing a gradient boosting decision tree (GBDT) algorithm in conjunction with the MIN + GB simulation protocol achieved the best overall performance. Its scoring power, ranking power and screening power significantly outperformed the Glide SF. First, compared with Glide, the average values of mean absolute error and root mean square error of GBDT/MIN + GB decreased about 38 and 36%, respectively. Second, the mean values of squared correlation coefficient and predictive index increased about 225 and 73%, respectively. Third, more encouragingly, the average value of the areas under the curve of receiver operating characteristic for six targets by GBDT, 0.87, is significantly better than that by Glide, which is only 0.71. Thus, we expected IP-SFs to have broad and promising applications in SBVSs.

In silico binding profile characterization of SARS-CoV-2 spike protein and its mutants bound to human ACE2 receptor.

Severe acute respiratory syndrome coronavirus (SARS-CoV-2), a novel coronavirus, has brought an unprecedented pandemic to the world and affected over 64 million people. The virus infects human using its spike glycoprotein mediated by a crucial area, receptor-binding domain (RBD), to bind to the human ACE2 (hACE2) receptor. Mutations on RBD have been observed in different countries and classified into nine types: A435S, D364Y, G476S, N354D/D364Y, R408I, V341I, V367F, V483A and W436R. Employing molecular dynamics (MD) simulation, we investigated dynamics and structures of the complexes of the prototype and mutant types of SARS-CoV-2 spike RBDs and hACE2. We then probed binding free energies of the prototype and mutant types of RBD with hACE2 protein by using an end-point molecular mechanics Poisson Boltzmann surface area (MM-PBSA) method. According to the result of MM-PBSA binding free energy calculations, we found that V367F and N354D/D364Y mutant types showed enhanced binding affinities with hACE2 compared to the prototype. Our computational protocols were validated by the successful prediction of relative binding free energies between prototype and three mutants: N354D/D364Y, V367F and W436R. Thus, this study provides a reliable computational protocol to fast assess the existing and emerging RBD mutations. More importantly, the binding hotspots identified by using the molecular mechanics generalized Born surface area (MM-GBSA) free energy decomposition approach can guide the rational design of small molecule drugs or vaccines free of drug resistance, to interfere with or eradicate spike-hACE2 binding.

Phosphorylation at Ser289 Enhances the Oligomerization of Tau Repeat R2.

In tauopathies such as Alzheimer's disease (AD), aberrant phosphorylation causes the dissociation of tau proteins from microtubules. The dissociated tau then aggregates into sequent forms from soluble oligomers to paired helical filaments and insoluble neurofibrillary tangles (NFTs). NFTs is a hallmark of AD, while oligomers are found to be the most toxic form of the tau aggregates. Therefore, understanding tau oligomerization with regard to abnormal phosphorylation is important for the therapeutic development of AD. In this study, we investigated the impact of phosphorylated Ser289, one of the 40 aberrant phosphorylation sites of full-length tau proteins, on monomeric and dimeric structures of tau repeat R2 peptides. We carried out intensive replica exchange molecular dynamics simulation with a total simulation time of up to 0.1 ms. Our result showed that the phosphorylation significantly affected the structures of both the monomer and the dimer. For the monomer, the phosphorylation enhanced ordered-disordered structural transition and intramolecular interaction, leading to more compactness of the phosphorylated R2 compared to the wild-type one. As to the dimer, the phosphorylation increased intermolecular interaction and ß-sheet formation, which can accelerate the oligomerization of R2 peptides. This result suggests that the phosphorylation at Ser289 is likely to promote tau aggregation. We also observed a phosphorylated Ser289-Na+-phosphorylated Ser289 bridge in the phosphorylated R2 dimer, suggesting an important role of cation ions in tau aggregation. Our findings suggest that phosphorylation at Ser289 should be taken into account in the inhibitor screening of tau oligomerization.

In Silico Screening of Natural Flavonoids against 3-Chymotrypsin-like Protease of SARS-CoV-2 Using Machine Learning and Molecular Modeling.

The "Long-COVID syndrome" has posed significant challenges due to a lack of validated therapeutic options. We developed a novel multi-step virtual screening strategy to reliably identify inhibitors against 3-chymotrypsin-like protease of SARS-CoV-2 from abundant flavonoids, which represents a promising source of antiviral and immune-boosting nutrients. We identified 57 interacting residues as contributors to the protein-ligand binding pocket. Their energy interaction profiles constituted the input features for Machine Learning (ML) models. The consensus of 25 classifiers trained using various ML algorithms attained 93.9% accuracy and a 6.4% false-positive-rate. The consensus of 10 regression models for binding energy prediction also achieved a low root-mean-square error of 1.18 kcal/mol. We screened out 120 flavonoid hits first and retained 50 drug-like hits after predefined ADMET filtering to ensure bioavailability and safety profiles. Furthermore, molecular dynamics simulations prioritized nine bioactive flavonoids as promising anti-SARS-CoV-2 agents exhibiting both high structural stability (root-mean-square deviation < 5 Å for 218 ns) and low MM/PBSA binding free energy (<-6 kcal/mol). Among them, KB-2 (PubChem-CID, 14630497) and 9-O-Methylglyceofuran (PubChem-CID, 44257401) displayed excellent binding affinity and desirable pharmacokinetic capabilities. These compounds have great potential to serve as oral nutraceuticals with therapeutic and prophylactic properties as care strategies for patients with long-COVID syndrome.

In silico binding affinity prediction for metabotropic glutamate receptors using both endpoint free energy methods and a machine learning-based scoring function.

Metabotropic glutamate receptors (mGluRs) play an important role in regulating glutamate signal pathways, which are involved in neuropathy and periphery homeostasis. mGluR4, which belongs to Group III mGluRs, is most widely distributed in the periphery among all the mGluRs. It has been proved that the regulation of this receptor is involved in diabetes, colorectal carcinoma and many other diseases. However, the application of structure-based drug design to identify small molecules to regulate the mGluR4 receptor is limited due to the absence of a resolved mGluR4 protein structure. In this work, we first built a homology model of mGluR4 based on a crystal structure of mGluR8, and then conducted hierarchical virtual screening (HVS) to identify possible active ligands for mGluR4. The HVS protocol consists of three hierarchical filters including Glide docking, molecular dynamic (MD) simulation and binding free energy calculation. We successfully prioritized active ligands of mGluR4 from a set of screening compounds using HVS. The predicted active ligands based on binding affinities can almost cover all the experiment-determined active ligands, with only one ligand missed. The correlation between the measured and predicted binding affinities is significantly improved for the MM-PB/GBSA-WSAS methods compared to the Glide docking method. More importantly, we have identified hotspots for ligand binding, and we found that SER157 and GLY158 tend to contribute to the selectivity of mGluR4 ligands, while ALA154 and ALA155 could account for the ligand selectivity to mGluR8. We also recognized other 5 key residues that are critical for ligand potency. The difference of the binding profiles between mGluR4 and mGluR8 can guide us to develop more potent and selective modulators. Moreover, we evaluated the performance of IPSF, a novel type of scoring function trained by a machine learning algorithm on residue-ligand interaction profiles, in guiding drug lead optimization. The cross-validation root-mean-square errors (RMSEs) are much smaller than those by the endpoint methods, and the correlation coefficients are comparable to the best endpoint methods for both mGluRs. Thus, machine learning-based IPSF can be applied to guide lead optimization, albeit the total number of actives/inactives are not big, a typical scenario in drug discovery projects.

A multiple-step in silico screening protocol to identify allosteric inhibitors of Spike-hACE2 binding.

While the COVID-19 pandemic continues to worsen, effective medicines that target the life cycle of SARS-CoV-2 are still under development. As more highly infective and dangerous variants of the coronavirus emerge, the protective power of vaccines will decrease or vanish. Thus, the development of drugs, which are free of drug resistance is direly needed. The aim of this study is to identify allosteric binding modulators from a large compound library to inhibit the binding between the Spike protein of the SARS-CoV-2 virus and human angiotensin-converting enzyme 2 (hACE2). The binding of the Spike protein to hACE2 is the first step of the infection of host cells by the coronavirus. We first built a compound library containing 77 448 antiviral compounds. Molecular docking was then conducted to preliminarily screen compounds which can potently bind to the Spike protein at two allosteric binding sites. Next, molecular dynamics simulations were performed to accurately calculate the binding affinity between the spike protein and an identified compound from docking screening and to investigate whether the compound can interfere with the binding between the Spike protein and hACE2. We successfully identified two possible drug binding sites on the Spike protein and discovered a series of antiviral compounds which can weaken the interaction between the Spike protein and hACE2 receptor through conformational changes of the key Spike residues at the Spike-hACE2 binding interface induced by the binding of the ligand at the allosteric binding site. We also applied our screening protocol to another compound library which consists of 3407 compounds for which the inhibitory activities of Spike/hACE2 binding were measured. Encouragingly, in vitro data supports that the identified compounds can inhibit the Spike-ACE2 binding. Thus, we developed a promising computational protocol to discover allosteric inhibitors of the binding of the Spike protein of SARS-CoV-2 to the hACE2 receptor, and several promising allosteric modulators were discovered.

Anti-inflammatory effects of Aureusidin in LPS-stimulated RAW264.7 macrophages via suppressing NF-κB and activating ROS- and MAPKs-dependent Nrf2/HO-1 signaling pathways.

Aureusidin, a naturally-occurring flavonoid, is found in various plants of Cyperaceae such as Heleocharis dulcis (Burm. f.) Trin., but its pharmacological effect and active mechanism are rarely reported. This study aimed to investigate the anti-inflammatory effect and action mechanism of Aureusidin in LPS-induced mouse macrophage RAW264.7 cells. The results suggested that lipopolysaccharide (LPS)-induced nitric oxide (NO), tumor necrosis factor-α (TNF-α) and prostaglandin E2 (PGE2) production were obviously inhibited by Aureusidin. Moreover, Aureusidin also significantly decreased the mRNA expression of various inflammatory factors in LPS-stimulated RAW264.7 cells. Furthermore, mechanistic studies showed that Aureusidin significantly inhibited nuclear transfer of nuclear factor-κB (NF-κB), while increasing the nuclear translocation of nuclear factor E2-related factor 2 (Nrf2) as well as expression of Nrf2 target genes such as heme oxygenase (HO-1) and NAD(P)H:quinone oxidoreductase 1 (NQO1), but the addition of the HO-1 inhibitor Sn-protoporphyrin (Snpp) significantly abolished the anti-inflammatory effect of Aureusidin in LPS-stimulated RAW264.7 cells, confirming the view that HO-1 was involved in the anti-inflammatory effect. In addition, Aureusidin increased the levels of reactive oxygen species (ROS) and mitogen-activated protein kinase (MAPK) phosphorylation in RAW264.7 cells. Antioxidant N-acetylcysteine (NAC) or three MAPK inhibitors blocked the nuclear translocation of Nrf2 and HO-1 expression induced by Aureusidin, indicating that Aureusidin activated the Nrf2/HO-1 signaling pathway through ROS and MAPKs pathways. At the same time, co-treatment with the NAC blocked the phosphorylation of MAPKs. Results from molecular docking indicated that Aureusidin inhibited the NF-κB pathway by covalently binding to NF-κB. Thus, Aureusidin exerted the anti-inflammatory activity through blocking the NF-κB signaling pathways and activating the MAPKs and Nrf2/HO-1 signaling pathways. Based on the above results, Aureusidin may be an attractive therapeutic candidate for the inflammation-related diseases.

Geometry Optimization Algorithms in Conjunction with the Machine Learning Potential ANI-2x Facilitate the Structure-Based Virtual Screening and Binding Mode Prediction.

Structure-based virtual screening utilizes molecular docking to explore and analyze ligand-macromolecule interactions, crucial for identifying and developing potential drug candidates. Although there is availability of several widely used docking programs, the accurate prediction of binding affinity and binding mode still presents challenges. In this study, we introduced a novel protocol that combines our in-house geometry optimization algorithm, the conjugate gradient with backtracking line search (CG-BS), which is capable of restraining and constraining rotatable torsional angles and other geometric parameters with a highly accurate machine learning potential, ANI-2x, renowned for its precise molecular energy predictions reassembling the wB97X/6-31G(d) model. By integrating this protocol with binding pose prediction using the Glide, we conducted additional structural optimization and potential energy prediction on 11 small molecule-macromolecule and 12 peptide-macromolecule systems. We observed that ANI-2x/CG-BS greatly improved the docking power, not only optimizing binding poses more effectively, particularly when the RMSD of the predicted binding pose by Glide exceeded around 5 Å, but also achieving a 26% higher success rate in identifying those native-like binding poses at the top rank compared to Glide docking. As for the scoring and ranking powers, ANI-2x/CG-BS demonstrated an enhanced performance in predicting and ranking hundreds or thousands of ligands over Glide docking. For example, Pearson's and Spearman's correlation coefficients remarkedly increased from 0.24 and 0.14 with Glide docking to 0.85 and 0.69, respectively, with the addition of ANI-2x/CG-BS for optimizing and ranking small molecules binding to the bacterial ribosomal aminoacyl-tRNA receptor. These results suggest that ANI-2x/CG-BS holds considerable potential for being integrated into virtual screening pipelines due to its enhanced docking performance.

Comparison and summary of in silico prediction tools for CYP450-mediated drug metabolism.

The cytochrome P450 (CYP450) enzyme system is responsible for the metabolism of more than two-thirds of xenobiotics. This review summarizes reports of a series of in silico tools for CYP450 enzyme-drug interaction predictions, including the prediction of sites of metabolism (SOM) of a drug and the identification of inhibitor/substrates for CYP subtypes. We also evaluated four prediction tools to identify CYP inhibitors utilizing 52 of the most frequently prescribed drugs. ADMET Predictor and CYPlebrity demonstrated the best performance. We hope that this review provides guidance for choosing appropriate enzyme prediction tools from a variety of in silico platforms to meet individual needs. Such predictions are useful for medicinal chemists to prioritize their designed compounds for further drug discovery.

Physiologically-Based Pharmacokinetics Modeling for Hydroxychloroquine as a Treatment for Malaria and Optimized Dosing Regimens for Different Populations.

Malaria is a severe parasite infectious disease with high fatality. As one of the approved treatments of this disease, hydroxychloroquine (HCQ) lacks clinical administration guidelines for patients with special health conditions and co-morbidities. This may result in improper dosing for different populations and lead them to suffer from severe side effects. One of the most important toxicities of HCQ overdose is cardiotoxicity. In this study, we built and validated a physiologically based pharmacokinetic modeling (PBPK) model for HCQ. With the full-PBPK model, we predicted the pharmacokinetic (PK) profile for malaria patients without other co-morbidities under the HCQ dosing regimen suggested by Food and Drug Administration (FDA) guidance. The PK profiles for different special populations were also predicted and compared to the normal population. Moreover, we proposed a series of adjusted dosing regimens for different populations with special health conditions and predicted the concentration-time (C-T) curve of the drug plasma concentration in these populations which include the pregnant population, elderly population, RA patients, and renal impairment populations. The recommended special population-dependent dosage regimens can maintain the similar drug levels observed in the virtual healthy population under the original dosing regimen provided by FDA. Last, we developed mathematic formulas for predicting dosage based on a patient's body measurements and two indexes of renal function (glomerular filtration rate and serum creatine level) for the pediatric and morbidly obese populations. Those formulas can facilitate personalized treatment of this disease. We hope to provide some advice to clinical practice when taking HCQ as a treatment for malaria patients with special health conditions or co-morbidities so that they will not suffer from severe side effects due to higher drug plasma concentration, especially cardiotoxicity.

In Silico Prediction of Pharmacokinetic Profile for Human Oral Drug Candidates Which Lack Clinical Pharmacokinetic Experiment Data.

BACKGROUNDS AND OBJECTIVES: In silico methods which can generate high-quality physiologically based pharmacokinetic (PBPK) models for arbitrary drug candidates are greatly needed to select developable drug candidates that escape drug attrition because of the poor pharmacokinetic profile. The purpose of this study is to develop a novel protocol to preliminarily predict the concentration profile of a target drug based on the PBPK model of a structurally similar template drug by combining two software platforms for PBPK modeling, the SimCYP simulator and ADMET Predictor. METHODS: The method was evaluated by utilizing 13 drug pairs from 18 drugs in the built-in database of the SimCYP software. All drug pairs have Tanimoto scores (TS) no less than 0.5. As each drug in a drug pair can serve as both target and template, 26 sets were studied in this work. Three versions (V1, V2 and V3) of models for the target drug were constructed by replacing the corresponding parameters of the template drug step by step with those predicted by ADMET Predictor for the target drug. V1 represents the replacement of molecular weight (MW), V2 includes the replacement of parameter MW, fraction unbound in plasma (fu), blood-to-plasma partition ratio (B/P), logarithm of the octanol-buffer partition coefficient (log Po:w) and acid dissociation constant (pKa). In V3, all above-mentioned parameters as well as human jejunum effective permeability (Peff), Vd and cytochrome P450 (CYP) metabolism parameters (Km, Vmax or CLint) are modified. Normalized root mean square error (NRMSE) was used for the evaluation of the model performance. RESULTS: We found that the performance of the three versions of the models depends on structural similarity of the drug pairs. For Group I drug pairs (TS ≤ 0.7), V2 and V3 performed better than V1 in terms of NRMSE; for Group II drug pairs (0.7 < TS ≤ 0.9), 8 out of 10 V3 models had NRMSE < 0.2, the cutoff we applied to judge whether the simulated concentration-time (C-T) curve was satisfactory or not. V3 outperformed the V1 and V2 versions. For the two drug pairs belonging to Group III (TS > 0.9), V2 outperformed V1 and V3, suggesting more unnecessary replacement can lower the performance of PBPK models. We also investigated how the prediction accuracy of ADMET Predictor as well as its collaboration with SimCYP influences the quality of PBPK models constructed using SimCYP. CONCLUSION: In conclusion, we generated practical guidance on applying two mainstream software packages, ADMET Predictor and SimCYP, to construct PBPK models for drugs or drug candidates that lack ADME parameters in model construction.

Incorporating structural similarity into a scoring function to enhance the prediction of binding affinities.

In this study, we developed a novel algorithm to improve the screening performance of an arbitrary docking scoring function by recalibrating the docking score of a query compound based on its structure similarity with a set of training compounds, while the extra computational cost is neglectable. Two popular docking methods, Glide and AutoDock Vina were adopted as the original scoring functions to be processed with our new algorithm and similar improvement performance was achieved. Predicted binding affinities were compared against experimental data from ChEMBL and DUD-E databases. 11 representative drug receptors from diverse drug target categories were applied to evaluate the hybrid scoring function. The effects of four different fingerprints (FP2, FP3, FP4, and MACCS) and the four different compound similarity effect (CSE) functions were explored. Encouragingly, the screening performance was significantly improved for all 11 drug targets especially when CSE = S4 (S is the Tanimoto structural similarity) and FP2 fingerprint were applied. The average predictive index (PI) values increased from 0.34 to 0.66 and 0.39 to 0.71 for the Glide and AutoDock vina scoring functions, respectively. To evaluate the performance of the calibration algorithm in drug lead identification, we also imposed an upper limit on the structural similarity to mimic the real scenario of screening diverse libraries for which query ligands are general-purpose screening compounds and they are not necessarily structurally similar to reference ligands. Encouragingly, we found our hybrid scoring function still outperformed the original docking scoring function. The hybrid scoring function was further evaluated using external datasets for two systems and we found the PI values increased from 0.24 to 0.46 and 0.14 to 0.42 for A2AR and CFX systems, respectively. In a conclusion, our calibration algorithm can significantly improve the virtual screening performance in both drug lead optimization and identification phases with neglectable computational cost.

Gambogic acid induces heme oxygenase-1 through Nrf2 signaling pathway and inhibits NF-κB and MAPK activation to reduce inflammation in LPS-activated RAW264.7 cells.

Gambogic acid (GA), a natural product with a xanthone structure, was previously demonstrated to exert anti-inflammatory effects. The aim of this study was to evaluate the anti-inflammatory activity of GA on LPS-stimulated mouse macrophage RAW264.7 and its anti-inflammatory mechanism. Pretreatment with GA inhibited LPS-induced production of nitric oxide (NO) and prostaglandin E2 (PGE2) through reducing the expressions of inducible nitric oxide synthase (iNOS) and cyclooxygenase-2 (COX-2). GA also decreased the expressions of pro-inflammatory cytokines, including TNF-α, IL-6 and IL-1ß. The activation of nuclear factor-κB (NF-κB) and the Mitogen activated phosphokinases (MAPKs) regulates pro-inflammatory factors. Further experiments demonstrated that the nuclear translocation of NF-κB, a promoting regulator in inflammation, was blocked via inhibiting the phosphorylation event of IκBα by GA. Meanwhile, the Mitogen activated phosphokinase (MAPK) signaling pathways were also suppressed. However, activation of nuclear factor erythroid 2-related factor (Nrf2) can inhibit inflammation. GA could activate the nucleus translocation of Nrf2 and up-regulated the expression of heme oxygenase-1 (HO-1). Taken together, GA exhibited its anti-inflammatory activities through Nrf2 activation and NF-κB depression, thus could be a candidate for the prevention and treatment of diseases that involve excessive inflammation.