Exploring the effect of inhibitor AKB-9778 on VE-PTP by molecular docking and molecular dynamics simulation

Engineering of human tryptophan hydroxylase 2 for efficient synthesis of 5-hydroxytryptophan

L-Tryptophan hydroxylation catalyzed by tryptophan hydroxylase (TPH) presents a promising method for synthesizing 5-hydroxytryptophan (5-HTP), yet the limited activity of wild-type human TPH2 restricts its application. A high-activity mutant, MT10 (H318E/H323E), was developed through semi-rational active site saturation testing (CAST) of wild-type TPH2, exhibiting a 2.85-fold increase in kcat/Km over the wild type, thus enhancing catalytic efficiency. Two biotransformation systems were developed, including an in vitro one-pot system and a Whole-Cell Catalysis System (WCCS). In the WCCS, MT10 achieved a conversion rate of only 31.5 % within 32 h. In the one-pot reaction, MT10 converted 50 mM L-tryptophan to 44.5 mM 5-HTP within 8 h, achieving an 89 % conversion rate, outperforming the M1 (NΔ143/CΔ26) variant. Molecular dynamics simulations indicated enhanced interactions of MT10 with the substrate, suggesting improved binding affinity and system stability. This study offers an effective approach for the efficient production of 5-HTP.

Identification of natural xanthine oxidase inhibitors: Virtual screening, anti-xanthine oxidase activity, and interaction mechanism

Xanthine oxidase (XO) is a crucial target for hyperuricemia treatment(s). Naturally occurred XO inhibitors with minimal toxicity and high efficacy have attracted researchers' attention. With the goal of quickly identifying natural XO inhibitors, an integrated computational screening strategy was constructed by molecular docking and calculating the free energy of binding. Twenty-seven hits were achieved from a database containing 19,377 natural molecules. This includes fourteen known XO inhibitors and four firstly-reported inhibitors (isolicoflavonol, 5,7-dihydroxycoumarin, parvifolol D and clauszoline M, IC50 < 40 μM). Iolicoflavonol (hit 8, IC50 = 8.45 ± 0.68 μM) and 5,7-dihydroxycoumarin (hit 25, IC50 = 10.91 ± 0.71 μM) displayed the great potency as mixed-type inhibitors. Docking study and molecular dynamics simulation revealed that both hits could interact with XO's primarily active site residues ARG880, MOS1328, and ASN768 of XO. Fluorescence spectroscopy studies showed that hit 8 bound to the active cavity region of XO, causing changes in XO's conformation and hydrophobicity. Hits 8 and 25 exhibit favorable Absorption, Distribution, Metabolism, and Excretion (ADME) properties. Additionally, no cytotoxicity against human liver cells was observed at their median inhibition concentrations against XO. Therefore, the present study offers isolicoflavonol and 5,7-dihydroxycoumarin with the potential to be disease-modifying agents for hyperuricemia.

Lead generation of UPPS inhibitors targeting MRSA: Using 3D-QSAR pharmacophore modeling, virtual screening, molecular docking, and molecular dynamic simulations

Undecaprenyl Pyrophosphate Synthase (UPPS) is a vital target enzyme in the early stages of bacterial cell wall biosynthesis. UPPS inhibitors have antibacterial activity against resistant strains such as MRSA and VRE. In this study, we used several consecutive computer-based protocols to identify novel UPPS inhibitors. The 3D QSAR pharmacophore model generation (HypoGen algorithm) protocol was used to generate a valid predictive pharmacophore model using a set of UPPS inhibitors with known reported activity. The developed model consists of four pharmacophoric features: one hydrogen bond acceptor, two hydrophobic, and one aromatic ring. It had a correlation coefficient of 0.86 and a null cost difference of 191.39, reflecting its high predictive power. Hypo1 was proven to be statistically significant using Fischer's randomization at a 95% confidence level. The validated pharmacophore model was used for the virtual screening of several databases. The resulting hits were filtered using SMART and Lipinski filters. The hits were docked into the binding site of the UPPS protein, affording 70 hits with higher docking affinities than the reference compound (6TC, - 21.17 kcal/mol). The top five hits were selected through extensive docking analysis and visual inspection based on docking affinities, fit values, and key residue interactions with the UPPS receptor. Moreover, molecular dynamic simulations of the top hits were performed to confirm the stability of the protein-ligand complexes, yielding five promising novel UPPS inhibitors.

Design, synthesis, activity and molecular dynamics studies of 1,3,4-thiadiazole derivatives as selective allosteric inhibitors of SHP2 for the treatment of cancer

Overexpression or gene mutation of SHP2 is closely linked with a variety of cancers and has been identified as a crucial anticancer target. In the study, we took SHP2 allosteric inhibitor SHP099 as the lead compound, and 32 1,3,4-thiadiazole derivatives were identified as selective allosteric inhibitors of SHP2. In vitro enzyme activity test showed that some compounds had high inhibition on full length SHP2, and almost no activity on homologous protein SHP1, exhibiting high selectivity. Compound YF704 (4w) had the best inhibition activity, with IC₅₀ value of 0.25 ± 0.02 μM, and also showed strong inhibitory activity on SHP2-E76K and SHP2-E76A, with IC₅₀ values of 6.88 ± 0.69 μM and 1.38 ± 0.12 μM, respectively. CCK8 proliferation test found that multiple compounds would effectively inhibit the proliferation of a variety of cancer cells. Among them, the IC₅₀ values of compound YF704 on MV4-11 and NCI-H358 cells were 3.85 ± 0.34 μM and 12.01 ± 0.62 μM, respectively. Specially, these compounds were sensitive to NCI-H358 cells containing KRAS^G12C mutation, thus overcoming the problem that SHP099 was insensitive to such cells. Apoptosis experiment showed that compound YF704 would effectively induce apoptosis of MV4-11 cells. Western blot showed that compound YF704 would downregulate the phosphorylation levels of Erk1/2 and Akt in MV4-11 and NCI-H358 cells. Molecular docking study show that compound YF704 would effectively bind to the allosteric region of SHP2 and form hydrogen bond interactions with key residues Thr108, Arg111 and Phe113. Molecular dynamics study further revealed the binding mechanism of SHP2 and compound YF704. In conclusion, we hope to provide potential SHP2 selective inhibitors and provide valuable clues for cancer treatment.

Identification of natural product inhibitors of PTP1B based on high-throughput virtual screening strategy: In silico, in vitro and in vivo studies

Protein tyrosine phosphatase 1B (PTP1B) is a key negative regulator of the insulin signaling pathway, which is a potential therapeutic target for the treatment of type 2 diabetes mellitus (T2DM). In this study, we identified several PTP1B inhibitors with high activity by using high-throughput virtual screening and in vitro enzyme inhibition activity verification strategies. Among them, baicalin was first reported as a selective mixed inhibitor of PTP1B, with IC₅₀ value of 3.87 ± 0.45 μM, and its inhibitory activity against homologous proteins TCPTP, SHP2, and SHP1 exceeded 50 μM. Molecular docking study found that baicalin and PTP1B could bind stably, and revealed that baicalin had a dual inhibitory effect. Cell experiments showed that baicalin was almost non-toxic and could significantly enhance the phosphorylation of IRS-1 in C2C12 myotube cells. Animal experiments showed that baicalin could significantly reduce the blood sugar of STZ-induced diabetic mice models, and had a liver protective effect. In conclusion, this study can provide new ideas for the development of PTP1B selective inhibitors.

Physicochemical properties of soybean β-conglycinin-based films affected by linoleic acid

To explore the interaction between lipids and proteins during emulsion film formation, the linoleic acid concentration effects on the physicochemical properties of soybean β-conglycinin (7S) films were studied. The viscosity and size of oil droplets in the film-forming solution gradually increased with the increasing linoleic acid concentration. As the linoleic acid concentration increased, the number of oil droplets on the film surfaces and elongation at break of films gradually increased, whereas the tensile strength decreased. The films containing 20% linoleic acid had the highest water vapor permeability value, which was decreased by increasing or decreasing the linoleic acid concentration. According to the molecular dynamics simulation and chemical interactions, 7S could be adsorbed at the linoleic acid interface and bind stably, resulting in the decreased ionic and hydrogen bonds but the increased hydrophobic interactions and covalent bonds among proteins in the films.

Novel Detection Method for Evaluating the Activity of an Alkaline Serine Protease from Bacillus clausii

Until now, the detection methods for serine proteases have been quite time-consuming or cannot indicate the "real" protease activity. Here, a rapid and simple method for determining the "real" activity of serine proteases toward AAPX (a kind of mixed polypeptide substrates, with X representing 20 standard amino acids) was developed. This AAPX method has high reliability, sensitivity, and repeatability and can be used for detecting the serine protease activity spectrophotometrically. Additionally, the site-directed saturation mutagenesis library of alkaline serine protease PRO (BcPRO) from Bacillus clausii was screened with this AAPX method. Three beneficial mutants S99R, S99H, and S99W were identified, and S99W displayed the highest activity. In comparison to wild-type BcPRO, S99W exhibited enhanced catalytic performance toward eight AAPX monomers, and the molecular dynamics simulation revealed the mechanism responsible for its improved activity toward AAPM. Consequently, this work provides an efficient method for detecting, characterizing, mining, and high-throughput screening of serine proteases.

Identification of potential target endoribonuclease NSP15 inhibitors of SARS-COV-2 from natural products through high-throughput virtual screening and molecular dynamics simulation

SARS‐CoV‐2 wreaks havoc around the world, triggering the COVID‐19 pandemic. It has been confirmed that the endoribonuclease NSP15 is crucial to the viral replication, and thus identified as a potential drug target against COVID‐19. The NSP15 protein was used as the target to conduct high‐throughput virtual screening on 30,926 natural products from the NPASS database to identify potential NSP15 inhibitors. And 100 ns molecular dynamics simulations were performed on the NSP15 and NSP15‐NPC198199 system. In all, 10 natural products with high docking scores with NSP15 protein were obtained, among which compound NPC198199 scored the highest. The analysis of the binding mode between NPC198199 and NSP15 found that NPC198199 would form H‐bond interactions with multiple key residues at the catalytic site. Subsequently, a series of post‐dynamics simulation analyses (including RMSD, RMSF, PCA, DCCM, RIN, binding free energy, and H‐bond occupancy) were performed to further explore inhibitory mechanism of compound NPC198199 on NSP15 protein at the molecular level. The research strongly indicates that the 10 natural compounds screened can be used as potential inhibitors of NSP15, and provides valuable information for the subsequent drug discovery of anti‐SARS‐CoV‐2.

An Insight of RuBisCO Evolution through a Multilevel Approach

RuBisCO is the most abundant enzyme on earth; it regulates the organic carbon cycle in the biosphere. Studying its structural evolution will help to develop new strategies of genetic improvement in order to increase food production and mitigate CO₂ emissions. In the present work, we evaluate how the evolution of sequence and structure among isoforms I, II and III of RuBisCO defines their intrinsic flexibility and residue-residue interactions. To do this, we used a multilevel approach based on phylogenetic inferences, multiple sequence alignment, normal mode analysis, and molecular dynamics. Our results show that the three isoforms exhibit greater fluctuation in the loop between αB and βC, and also present a positive correlation with loop 6, an important region for enzymatic activity because it regulates RuBisCO conformational states. Likewise, an increase in the flexibility of the loop structure between αB and βC, as well as Lys330 (form II) and Lys322 (form III) of loop 6, is important to increase photosynthetic efficiency. Thus, the cross-correlation dynamics analysis showed changes in the direction of movement of the secondary structures in the three isoforms. Finally, key amino acid residues related to the flexibility of the RuBisCO structure were indicated, providing important information for its enzymatic engineering.

Thermodynamic stability, in-vitro permeability, and in-silico molecular modeling of the optimal Elaeis guineensis leaves extract water-in-oil nanoemulsion

Nanoemulsion is a delivery system used to enhance bioavailability of plant-based compounds across the stratum corneum. Elaeis guineensis leaves are rich source of polyphenolic antioxidants, viz. gallic acid and catechin. The optimal E. guineensis leaves extract water-in-oil nanoemulsion was stable against coalescence, but it was under significant influence of Ostwald ripening over 90 days at 25 °C. The in-vitro permeability revealed a controlled and sustained release of the total phenolic compounds (TPC) of EgLE with a cumulative amount of 1935.0 ± 45.7 µgcm-2 after 8 h. The steady-state flux and permeation coefficient values were 241.9 ± 5.7 µgcm-2 h-1 and 1.15 ± 0.03 cm.h-1, respectively. The kinetic release mechanism for TPC of EgLE was best described by the Korsmeyer-Peppas model due to the highest linearity of R2 = 0.9961, indicating super case II transport mechanism. The in-silico molecular modelling predicted that the aquaporin-3 protein in the stratum corneum bonded preferably to catechin over gallic acid through hydrogen bonds due to the lowest binding energies of - 57.514 kcal/mol and - 8.553 kcal/mol, respectively. Thus, the in-silico study further verified that catechin could improve skin hydration. Therefore, the optimal nanoemulsion could be used topically as moisturizer to enhance skin hydration based on the in-silico prediction.

Exploration of Novel Xanthine Oxidase Inhibitors Based on 1,6-Dihydropyrimidine-5-Carboxylic Acids by an Integrated in Silico Study

Xanthine oxidase (XO) is an important target for the effective treatment of hyperuricemia-associated diseases. A series of novel 2-substituted 6-oxo-1,6-dihydropyrimidine-5-carboxylic acids (ODCs) as XO inhibitors (XOIs) with remarkable activities have been reported recently. To better understand the key pharmacological characteristics of these XOIs and explore more hit compounds, in the present study, the three-dimensional quantitative structure-activity relationship (3D-QSAR), molecular docking, pharmacophore modeling, and molecular dynamics (MD) studies were performed on 46 ODCs. The constructed 3D-QSAR models exhibited reliable predictability with satisfactory validation parameters, including q2 = 0.897, R2 = 0.983, rpred2 = 0.948 in a CoMFA model, and q2 = 0.922, R2 = 0.990, rpred2 = 0.840 in a CoMSIA model. Docking and MD simulations further gave insights into the binding modes of these ODCs with the XO protein. The results indicated that key residues Glu802, Arg880, Asn768, Thr1010, Phe914, and Phe1009 could interact with ODCs by hydrogen bonds, π-π stackings, or hydrophobic interactions, which might be significant for the activity of these XOIs. Four potential hits were virtually screened out using the constructed pharmacophore model in combination with molecular dockings and ADME predictions. The four hits were also found to be relatively stable in the binding pocket by MD simulations. The results in this study might provide effective information for the design and development of novel XOIs.

Screening potential FDA-approved inhibitors of the SARS-CoV-2 major protease 3CLpro through high-throughput virtual screening and molecular dynamics simulation

It has been confirmed that the new coronavirus SARS-CoV-2 caused the global pandemic of coronavirus disease 2019 (COVID-19). Studies have found that 3-chymotrypsin-like protease (3CL^pro) is an essential enzyme for virus replication, and could be used as a potential target to inhibit SARS-CoV-2. In this work, 3CL^pro was used as the target to complete the high-throughput virtual screening of the FDA-approved drugs, and Indinavir and other 10 drugs with high docking scores for 3CL^pro were obtained. Studies on the binding pattern of 3CL^pro and Indinavir found that Indinavir could form the stable hydrogen bond (H-bond) interactions with the catalytic dyad residues His41-Cys145. Binding free energy study found that Indinavir had high binding affinity with 3CL^pro. Subsequently, molecular dynamics simulations were performed on the 3CL^pro and 3CL^pro-Indinavir systems, respectively. The post-dynamic analyses showed that the conformational state of the 3CLpro-Indinavir system transformed significantly and the system tended to be more stable. Moreover, analyses of the residue interaction network (RIN) and H-bond occupancy revealed that the residue-residue interaction at the catalytic site of 3CL^pro was significantly enhanced after binding with Indinavir, which in turn inactivated the protein. In short, through this research, we hope to provide more valuable clues against COVID-19.

A Small Molecule Inhibitor of VE-PTP Activates Tie2 in Schlemm's Canal Increasing Outflow Facility and Reducing Intraocular Pressure

Purpose

Tyrosine kinase with immunoglobulin-like and EGF-like domains 2 (Tie2) activation in Schlemm's canal (SC) endothelium is required for the maintenance of IOP, making the angiopoietin/Tie2 pathway a target for new and potentially disease modifying glaucoma therapies. The goal of the present study was to examine the effects of a Tie2 activator, AKB-9778, on IOP and outflow function.

Methods

AKB-9778 effects on IOP was evaluated in humans, rabbits, and mice. Localization studies of vascular endothelial protein tyrosine phosphatase (VE-PTP), the target of AKB-9778 and a negative regulator of Tie2, were performed in human and mouse eyes. Mechanistic studies were carried out in mice, monitoring AKB-9778 effects on outflow facility, Tie2 phosphorylation, and filtration area of SC.

Results

AKB-9778 lowered IOP in patients treated subcutaneously for diabetic eye disease. In addition to efficacious, dose-dependent IOP lowering in rabbit eyes, topical ocular AKB-9778 increased Tie2 activation in SC endothelium, reduced IOP, and increased outflow facility in mouse eyes. VE-PTP was localized to SC endothelial cells in human and mouse eyes. Mechanistically, AKB-9778 increased the filtration area of SC for aqueous humor efflux in both wild type and in Tie2+/- mice.

Conclusions

This is the first report of IOP lowering in humans with a Tie2 activator and functional demonstration of its action in remodeling SC to increase outflow facility and lower IOP in fully developed mice. Based on these studies, a phase II clinical trial is in progress to advance topical ocular AKB-9778 as a first in class, Tie2 activator for treatment for ocular hypertension and glaucoma.

Design, synthesis, biological evaluation and molecular dynamics studies of 4-thiazolinone derivatives as protein tyrosine phosphatase 1B (PTP1B) inhibitors

Protein tyrosine phosphatase 1B (PTP1B) is a key negative regulator of insulin signaling pathway, and more and more studies have shown that it is a potential target for the treatment of type 2 diabetes mellitus (T2DM). In this study, 17 new 4-thiazolinone derivatives were designed and synthesized as novel PTP1B inhibitors, and ADMET prediction confirmed that these compounds were to be drug-like. In vitro enzyme activity experiments were performed on these compounds, and it was found that a plurality of compounds had good inhibitory activity and high selectivity against PTP1B protein. Among them, compound 7p exhibited the best inhibitory activity with an IC₅₀ of 0.92 μM. The binding mode of compound 7p and PTP1B protein was explored, revealing the reason for its high efficiency. In addition, molecular dynamics simulations for the PTP1B^WT and PTP1B^comp#7p systems revealed the effects of compound 7p on PTP1B protein at the molecular level. In summary, the study reported for the first time that 4-thiazolinone derivatives as a novel PTP1B inhibitor had good inhibitory activity and selectivity for the treatment of T2DM, providing more options for the development of PTP1B inhibitors. AbbreviationsBBBblood-brain barrierCDC25Bcell division cycle 25 homolog BCYP2D6Cytochrome P450 2D6 bindingDCCMdynamic cross-correlation mapDSDiscovery StudioH bondhydrogen bondHIAhuman intestinal absorptionLARleukocyte antigen-related phosphataseMDmolecular dynamicsMEG-2maternal-effect germ-cell defective 2MM-PBSAmolecular mechanics Poisson Boltzmann surface area)PCAprincipal component analysisPDBProtein Data BankpNPPp-nitrophenyl phosphatePPBplasma protein bindingPTP1Bprotein tyrosine phosphotase 1BRMSDroot mean square deviationRMSFroot mean square fluctuationSHP-1src homologous phosphatase-1SHP-2src homologous phosphatase-2SPCsingle-point chargeTCPTPT cell protein tyrosine phosphataseT2DMType 2 diabetes mellitusVDWvan der WaalsCommunicated by Ramaswamy H. Sarma.