GROMOS96 43a1 performance on the characterization of glycoprotein conformational ensembles through molecular dynamics simulations

Discovery of Substituted 2-oxoquinolinylthiazolidin-4-one Analogues as Potential EGFRK Inhibitors in Lung Cancer Treatment

PURPOSE

Cancer is the second leading cause of death globally and is responsible for an estimated 9.6 million deaths in 2018. Globally, about 1 in 6 deaths is due to cancer and the chemotherapeutic drugs available have high toxicity and have reported side effects hence, there is a need for the synthesis of novel drugs in the treatment of cancer.

METHODS

The current research work dealt with the synthesis of a series of 3-(3-acetyl-2-oxoquinolin-1-(2H)-yl-2-(substitutedphenyl)thiazolidin-4-one (Va-j) derivatives and evaluation of their in-vitro anticancer activity. All the synthesized compounds were satisfactorily characterized by IR and NMR data. Compounds were further evaluated for their in-vitro anticancer activity against A-549 (lung cancer) cell lines. The in-vitro anticancer activity was based upon the 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazoliumbromide (MTT) assay method.

RESULTS

The synthesized compounds exhibited satisfactory anticancer properties against the A-549 cell line. The compound (VH): showed the highest potency amongst the tested derivatives against the A-549 cell line with IC50 values of 100 µg/ml respectively and was also found to be more potent than Imatinib (150 µg/ml) which was used as a standard drug. Molecular docking studies of the titled compounds (Va-j) were carried out using AutoDock Vina/PyRx software. The synthesized compounds exhibited well-conserved hydrogen bonds with one or more amino acid residues in the active pocket of the EGFRK tyrosine kinase domain (PDB 1m17).

CONCLUSION

Among all the synthesized analogues, the binding affinity of the compound (Vh) was found to be higher than other synthesized derivatives and a molecular dynamics simulation study explored the stability of the docked complex system.

Mss51 protein inhibition serves as a novel target for type 2 diabetes: a molecular docking and simulation study

Myostatin is a widely recognized inhibitory factor of skeletal muscle growth and significantly influences muscle development and metabolism. In mice, myostatin inhibition improves insulin sensitivity, increases glucose uptake by skeletal muscle, and reduces body fat. Furthermore, Mss51 is downregulated in response to myostatin inhibition, and its deletion appears to improve the metabolic state of skeletal muscle and reduce adipose tissue, which makes Mss51 a potential target for the treatment of obesity and type 2 diabetes. Here, we report a computationally predicted and validated three-dimensional structure of Mss51. Computational screening was used to identify naturally occurring compounds from the Herbal and Specs chemical database that might inhibit Mss51, based on binding affinities and physiochemical and ADMET properties. ZINC00338371, ZINC95099599 and ZINC08214878 were found to bind to Mss51 with high binding affinity and specificity. In addition, 100 ns molecular dynamics simulations were conducted to assess the stabilities of the interactions between the three compounds and Mss51. MD simulation demonstrated that all three compounds bind to the active pocket site of Mss51 stably and cause conformation changes. ZINC00338371 was found to bind most stably with binding free energy -229.022 ± 13.776 kJ/mol to Mss51, suggesting that it has therapeutic potential as a treatment option for obesity and type 2 diabetes.Communicated by Ramaswamy H. Sarma.

Expression analysis and mapping of Viral-Host Protein interactions of Poxviridae suggests a lead candidate molecule targeting Mpox

BACKGROUND

Monkeypox (Mpox) is an important human pathogen without etiological treatment. A viral-host interactome study may advance our understanding of molecular pathogenesis and lead to the discovery of suitable therapeutic targets.

METHODS

GEO Expression datasets characterizing mRNA profile changes in different host responses to poxviruses were analyzed for shared pathway identification, and then, the Protein-protein interaction (PPI) maps were built. The viral gene expression datasets of Monkeypox virus (MPXV) and Vaccinia virus (VACV) were used to identify the significant viral genes and further investigated for their binding to the library of targeting molecules.

RESULTS

Infection with MPXV interferes with various cellular pathways, including interleukin and MAPK signaling. While most host differentially expressed genes (DEGs) are predominantly downregulated upon infection, marked enrichments in histone modifiers and immune-related genes were observed. PPI analysis revealed a set of novel virus-specific protein interactions for the genes in the above functional clusters. The viral DEGs exhibited variable expression patterns in three studied cell types: primary human monocytes, primary human fibroblast, and HeLa, resulting in 118 commonly deregulated proteins. Poxvirus proteins C6R derived protein K7 and K7R of MPXV and VACV were prioritized as targets for potential therapeutic interventions based on their histone-regulating and immunosuppressive properties. In the computational docking and Molecular Dynamics (MD) experiments, these proteins were shown to bind the candidate small molecule S3I-201, which was further prioritized for lead development.

RESULTS

MPXV circumvents cellular antiviral defenses by engaging histone modification and immune evasion strategies. C6R-derived protein K7 binding candidate molecule S3I-201 is a priority promising candidate for treating Mpox.

Immunoinformatics-Driven Strategies for Advancing Epitope-Based Vaccine Design for West Nile Virus

The West Nile virus (WNV) is the causative agent of West Nile disease (WND), which poses a potential risk of meningitis or encephalitis. The aim of the study was to design an epitope-based vaccine for WNV by utilizing computational analyses. The epitope-based vaccine design process encompassed WNV sequence collection, phylogenetic tree construction, and sequence alignment. Computational models identified B-cell and T-cell epitopes, followed by immunological property analysis. Epitopes were then modeled and docked with B-cell receptors, MHC I, and MHC II. Molecular dynamics simulations further explored dynamic interactions between epitopes and receptors. The findings indicated that the B-cell epitope QINHHWHKSGSSIG, along with three T-cell epitopes (FLVHREWFM for MHC I, NPFVSVATANAKVLI for MHC II, and NAYYVMTVGTKTFLV for MHC II), successfully passed the immunological evaluations. These four epitopes were further subjected to docking and molecular dynamics simulation studies. Although each demonstrated favorable affinities with their respective receptors, only NAYYVMTVGTKTFLV displayed a stable interaction with MHC II during MDS analysis, hence emerging as a potential candidate for a WNV epitope-based vaccine. This study demonstrates a comprehensive approach to epitope vaccine design, combining computational analyses, molecular modeling, and simulation techniques to identify potential vaccine candidates for WNV.

DFT and molecular simulation validation of the binding activity of PDEδ inhibitors for repression of oncogenic k-Ras

The development of effective drugs targeting the K-Ras oncogene product is a significant focus in anticancer drug development. Despite the lack of successful Ras signaling inhibitors, recent research has identified PDEδ, a KRAS transporter, as a potential target for inhibiting the oncogenic KRAS signaling pathway. This study aims to investigate the interactions between eight K-Ras inhibitors (deltarazine, deltaflexin 1 and 2, and its analogues) and PDEδ to understand their binding modes. The research will utilize computational techniques such as density functional theory (DFT) and molecular electrostatic surface potential (MESP), molecular docking, binding site analyses, molecular dynamic (MD) simulations, electronic structure computations, and predictions of the binding free energy. Molecular dynamic simulations (MD) will be used to predict the binding conformations and pharmacophoric features in the active site of PDEδ for the examined structures. The binding free energies determined using the MMPB(GB)SA method will be compared with the observed potency values of the tested compounds. This computational approach aims to enhance understanding of the PDEδ selective mechanism, which could contribute to the development of novel selective inhibitors for K-Ras signaling.

Licochalcone A and B enhance muscle proliferation and differentiation by regulating Myostatin

BACKGROUND

Myostatin (MSTN) inhibition has demonstrated promise for the treatment of diseases associated with muscle loss. In a previous study, we discovered that Glycyrrhiza uralensis (G. uralensis) crude water extract (CWE) inhibits MSTN expression while promoting myogenesis. Furthermore, three specific compounds of G. uralensis, namely liquiritigenin, tetrahydroxymethoxychalcone, and Licochalcone B (Lic B), were found to promote myoblast proliferation and differentiation, as well as accelerate the regeneration of injured muscle tissue.

PURPOSE

The purpose of this study was to build on our previous findings on G. uralensis and demonstrate the potential of its two components, Licochalcone A (Lic A) and Lic B, in muscle mass regulation (by inhibiting MSTN), aging and muscle formation.

METHODS

G. uralensis, Lic A, and Lic B were evaluated thoroughly using in silico, in vitro and in vivo approaches. In silico analyses included molecular docking, and dynamics simulations of these compounds with MSTN. Protein-protein docking was carried out for MSTN, as well as for the docked complex of MSTN-Lic with its receptor, activin type IIB receptor (ACVRIIB). Subsequent in vitro studies used C2C12 cell lines and primary mouse muscle stem cells to acess the cell proliferation and differentiation of normal and aged cells, levels of MSTN, Atrogin 1, and MuRF1, and plasma MSTN concentrations, employing techniques such as western blotting, immunohistochemistry, immunocytochemistry, cell proliferation and differentiation assays, and real-time RT-PCR. Furthermore, in vivo experiments using mouse models focused on measuring muscle fiber diameters.

RESULTS

CWE of G. uralensis and two of its components, namely Lic A and B, promote myoblast proliferation and differentiation by inhibiting MSTN and reducing Atrogin1 and MuRF1 expressions and MSTN protein concentration in serum. In silico interaction analysis revealed that Lic A (binding energy -6.9 Kcal/mol) and B (binding energy -5.9 Kcal/mol) bind to MSTN and reduce binding between it and ACVRIIB, thereby inhibiting downstream signaling. The experimental analysis, which involved both in vitro and in vivo studies, demonstrated that the levels of MSTN, Atrogin 1, and MuRF1 were decreased when G. uralensis CWE, Lic A, or Lic B were administered into mice or treated in the mouse primary muscle satellite cells (MSCs) and C2C12 myoblasts. The diameters of muscle fibers increased in orally treated mice, and the differentiation and proliferation of C2C12 cells were enhanced. G. uralensis CWE, Lic A, and Lic B also promoted cell proliferation in aged cells, suggesting that they may have anti-muslce aging properties. They also reduced the expression and phosphorylation of SMAD2 and SMAD3 (MSTN downstream effectors), adding to the evidence that MSTN is inhibited.

CONCLUSION

These findings suggest that CWE and its active constituents Lic A and Lic B have anti-mauscle aging potential. They also have the potential to be used as natural inhibitors of MSTN and as therapeutic options for disorders associated with muscle atrophy.

Conformational dynamics of α-1 acid glycoprotein (AGP) in cancer: A comparative study of glycosylated and unglycosylated AGP

α-1 acid glycoprotein (AGP) is one of the most abundant plasma proteins. It fulfills two important functions: immunomodulation, and binding to various drugs and receptors. These different functions are closely associated and modulated via changes in glycosylation and cancer missense mutations. From a structural point of view, glycans alter the local biophysical properties of the protein leading to a diverse ligand-binding spectrum. However, glycans can typically not be observed in the resolved X-ray crystallography structure of AGP due to their high flexibility and microheterogeneity, so limiting our understanding of AGP's conformational dynamics 70 years after its discovery. We here investigate how mutations and glycosylation interfere with AGP's conformational dynamics changing its biophysical behavior, by using molecular dynamics (MD) simulations and sequence-based dynamics predictions. The MD trajectories show that glycosylation decreases the local backbone flexibility of AGP and increases the flexibility of distant regions through allosteric effects. We observe that mutations near the glycosylation site affect glycan's conformational preferences. Thus, we conclude that mutations control glycan dynamics which modulates the protein's backbone flexibility directly affecting its accessibility. These findings may assist in the drug design targeting AGP's glycosylation and mutations in cancer.

Targeting myostatin using quercetin as a media supplement to improve myogenesis for cultured meat production: An in silico and in vitro study

Cultured meat (CM) is an alternative protein food and is produced by cultivating muscle satellite (stem) cells (MSCs) derived from livestock animals (bovine, chickens, and porcine) through myogenesis leading to generate muscle mass. Myostatin (MSTN) is well well-known negative regulator of myogenesis, and in the present study, in silico screening of natural compounds was performed to identify MSTN inhibitors. Interestingly, quercetin was found to inhibit MSTN (binding energy -7.40 kcal/mol), and this was further validated by a 100 ns molecular dynamics simulation. Quercetin was added to culture media to boost myogenesis, and its potent antioxidant property helped maintain media pH. Furthermore, quercetin increased the myotube thickness and length, increased MSC differentiation, and upregulated the gene and protein expressions of myoblast determination protein 1 (MYOD), Myogenin (MYOG), and Myosin heavy chains (MYH) in vitro. In addition, quercetin inhibited the activities of MSTN, activin receptor type-2B (ACVR2B), and SMAD2 and 3, and thus significantly enhanced MSC differentiation and myotube formation. Overall, this study shows that quercetin might be useful for enhancing large-scale CM production. It is hoped that this study provides a starting point for research in the CM area aimed to enhancing product quality, nutritional values, and the efficacy of large-scale production.

2D-QSAR, molecular docking and MD simulation based virtual screening of the herbal molecules against Alzheimer's disorder: an approach to predict CNS activity

Acetylcholinesterase (AChE) is one of the key enzyme targets that have been used clinically for the management of Alzheimer's Disorder (AD). Numerous reports in the literature predict and demonstrate in-vitro, and in-silico anticholinergic activity of herbal molecules, however, majority of them failed to find clinical application. To address these issues, we developed a 2D-QSAR model that could efficiently predict the AChE inhibitory activity of herbal molecules along with predicting their potential to cross the blood-brain-barrier (BBB) to exert their beneficial effects during AD. Virtual screening of the herbal molecules was performed and amentoflavone, asiaticoside, astaxanthin, bahouside, biapigenin, glycyrrhizin, hyperforin, hypericin, and tocopherol were predicted as the most promising herbal molecules for inhibiting AChE. Results were validated through molecular docking, atomistic molecular dynamics simulations and Molecular mechanics-Poisson Boltzmann surface area (MM-PBSA) studies against human AChE (PDB ID: 4EY7). To determine whether or not these molecules can cross BBB to inhibit AChE within the central nervous system (CNS) for being beneficial for the management of AD, we determined a CNS Multi-parameter Optimization (MPO) score, which was found in the range of 1 to 3.76. Overall, the best results were observed for amentoflavone and our results demonstrated a PIC50 value of 7.377 nM, molecular docking score of -11.5 kcal/mol, and CNS MPO score of 3.76. In conclusion, we successfully developed a reliable and efficient 2D-QSAR model and predicted amentoflavone to be the most promising molecule that could inhibit human AChE enzyme within the CNS and could prove beneficial for the management of AD.Communicated by Ramaswamy H. Sarma.

Virtual Insights into Natural Compounds as Potential 5α-Reductase Type II Inhibitors: A Structure-Based Screening and Molecular Dynamics Simulation Study

Androgenic alopecia (AGA) is a dermatological disease with psychosocial consequences for those who experience hair loss. AGA is linked to an increase in androgen levels caused by an excess of dihydrotestosterone in blood capillaries produced from testosterone by 5α-reductase type II (5αR2), which is expressed in scalp hair follicles; 5αR2 activity and dihydrotestosterone levels are elevated in balding scalps. The diverse health benefits of flavonoids have been widely reported in epidemiological studies, and research interest continues to increase. In this study, a virtual screening approach was used to identify compounds that interact with active site residues of 5αR2 by screening a library containing 241 flavonoid compounds. Here, we report two potent flavonoid compounds, eriocitrin and silymarin, that interacted strongly with 5αR2, with binding energies of -12.1 and -11.7 kcal/mol, respectively, which were more significant than those of the control, finasteride (-11.2 kcal/mol). Molecular dynamic simulations (200 ns) were used to optimize the interactions between compounds and 5αR2 and revealed that the interaction of eriocitrin and silymarin with 5αR2 was stable. The study shows that eriocitrin and silymarin provide developmental bases for novel 5αR2 inhibitors for the management of AGA.

Underpinning Endogeneous Damp EDA-Fibronectin in the Activation of Molecular Targets of Rheumatoid Arthritis and Identifcation of its Effective Inhibitors by Computational Methods

Rheumatoid arthritis (RA) is one of the most severe inflammatory diseases that cause swelling, stiffness and pain in the joints, which pose a significant threat worldwide. Damage-associated molecular patterns (DAMPs) are danger molecules of endogenous origin, released during cell injury or cell death, interacts with various Pattern recognition receptors (PRRs) and activates various inflammatory diseases. One of the DAMP molecules, so-called EDA-fibronectin (Fn) is also responsible for causing RA. EDA-Fn triggers RA through its interaction with TLR4. Apart from TLR4, it is divulged that certain other PRR's are also responsible for RA, but the identity and mechanism of those PRRs remain unknown until now. Hence, for the first time, we tried to reveal those PRR's interaction with EDA-Fn in RA through computational methods. Protein-protein interaction (PPI) was checked using ClusPro between EDA-Fn and certain Pattern recognition receptors (PRRs) to explore the binding affinities of the potential PRRs. Protein-protein docking unveiled that TLR5, TLR2 and RAGE has good interaction with EDA-Fn than the well-reported TLR4. Macromolecular simulation was performed for TLR5, TLR2 and RAGE complexes along with the control group TLR4 for 50 ns to further investigate the stability, leading to the identification of TLR2, TLR5 and RAGE as the stable complexes. Hence, TLR2, TLR5 and RAGE on interaction with EDA-Fn may lead to the progression of RA that may need additional validations through in vitro and in vivo animal models. Molecular docking was used to analyse the binding force of the top 33 active anti-arthritic compounds with the target protein EDA-Fn. Molecular docking study showed that withaferin A has a good binding activity with EDA-fibronectin target. Hence, it is emphasized that guggulsterone and berberine could modulate the EDA-Fn-mediated TLR5/TLR2/RAGE pathways, thereby it could inhibit the deteriorating effects of RA which needs further in vitro and in vivo experimental validations.

Identification of hits as anti-obesity agents against human pancreatic lipase via docking, drug-likeness, in-silico ADME(T), pharmacophore, DFT, molecular dynamics, and MM/PB(GB)SA analysis

Obesity, characterized by excessive fat accumulation, is a major health concern. Inhibition of human pancreatic lipase, an enzyme involved in fat digestion, offers a potential strategy for weight loss and obesity treatment. This study aimed to identify polyphenols capable of forming stable complexes with human pancreatic lipase to block its activity. Molecular docking, density functional theory (DFT), molecular dynamics (MD) simulations, and MMPBGBSA calculations were employed to evaluate ligand binding, stability, and energy profiles. Pharmacophore modeling was also performed to identify key structural features for effective inhibition. Virtual screening identified ZINC000015120539, ZINC000000899200, ZINC000001531702, and ZINC000013340267 as potential candidates, exhibiting favorable binding and stable interactions over 100 ns MD simulations. These findings provide insights into the inhibitory potential of selected polyphenols on human pancreatic lipase and support further experimental investigations for obesity treatment.Communicated by Ramaswamy H. Sarma.

Identification of novel inhibitors of Neisseria gonorrhoeae MurI using homology modeling, structure-based pharmacophore, molecular docking, and molecular dynamics simulation-based approach

MurI is one of the most significant role players in the biosynthesis of the peptidoglycan layer in Neisseria gonorrhoeae (Ng). We attempted to highlight the structural and functional relationship between Ng-MurI and D-glutamate to design novel molecules targeting this interaction. The three-dimensional (3D) model of the protein was constructed by homology modeling and the quality and consistency of generated model were assessed. The binding site of the protein was identified by molecular docking studies and a pharmacophore was identified using the interactions of the control ligand. The structure-based pharmacophore model was validated and employed for high-throughput virtual screening and molecular docking to identify novel Ng-MurI inhibitors. Finally, the model was optimized by molecular dynamics (MD) simulations and the optimized model complex with the substrate glutamate and novel molecules facilitated us to confirm the stability of the protein-ligand docked complexes. The 100 ns MD simulations of the potential lead compounds with protein confirmed that the modeled complexes were stable. This study identifies novel potential compounds with good fitness and docking scores, which made the interactions of biological significance within the protein active site. Hence, the identified compounds may act as new leads to design and develop Ng-MurI inhibitors.Communicated by Ramaswamy H. Sarma.

Reducing the assemblies of amyloid-beta multimers by sodium dodecyl sulfate surfactant at concentrations lower than critical micelle concentration: molecular dynamics simulation exploration

Amyloid-β peptide, the predominant proteinaceous component of senile plaques, is responsible for the incidence of Alzheimer's disease (AD), an age-associated neurodegenerative disorder. Specifically, the amyloid-β(1-42) (Aβ1-42) isoform, known for its high toxicity, is the predominant biomarker for the preliminary diagnosis of AD. The aggregation of the Aβ1-42 peptides can be affected by the components of the cellular medium through changing their structures and molecular interactions. In this study, we investigated the effect of sodium dodecyl sulfate (SDS) at much lower concentrations than the critical micelle concentration (CMC) on Aβ1-42 aggregation. For this purpose, we studied mono-, di-, tri- and tetramers of Aβ1-42 peptide in two different concentrations of SDS molecules (10 and 40 molecules) using a 300 ns molecular dynamics simulation for each system. The distance between the center of mass (COM) of Aβ1-42 peptides confirms that an increase in the number of SDS molecules decreases their aggregation probability due to greater interaction with SDS molecules. Besides, the less compactness parameter reveals the reduced aggregation probability of Aβ1-42 peptides. Based on the energetic FEL landscapes, SDS molecules with the concentration closer to the CMC are an effective inhibitory agent to prevent the formation of Aβ1-42 fibrils. Also, the aggregation direction of the peptide pairs can be predicted by determining the direction of the accumulation-deterrent forces.Communicated by Ramaswamy H. Sarma.

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.

Effect of co-pigments on anthocyanins of Rhododendron arboreum and insights into interaction mechanism

The impact of intermolecular copigmentation between five phenolic acids, two flavonoid and three amino acids with R. arboreum anthocyanins (ANS) and its isolated cyanidin-3-O-monoglycosides were investigated through experimental and theoretical approach. On addition of different copigments, phenolic acid induced strong hyperchromic (0.26-0.55 nm) and bathochromic shift (6.6-14.2 nm). The color intensity and stability of ANS with, storage at 4 °C & 25 °C, sunlight, oxidation and heat were evaluated by chromaticity, anthocyanin content, kinetic and structural simulation analysis. The strongest copigmentation reaction was observed with narningin (NA) and also showed high thermostability and highest half-life i.e. 3.39 h-1.24 h at 90-160 °C. The cyanidin-3-O-monoglycosides were analysed for their copigmentation effect and observations revealed that NA displayed best copigmentation effect to cyanidin-3-O-arabinoside (B) followed by cyanidin-3-O-galactoside (A), and cyanidin-3-O-rhamnoside (C). Additionally, structural simulation and steered molecular dynamics insights NA is the most favourable co-pigment involving π-π stacking and H-bonding.

A biallelic variant of DCAF13 implicated in a neuromuscular disorder in humans

Neuromuscular disorders encompass a broad range of phenotypes and genetic causes. We investigated a consanguineous family in which multiple patients had a neuromuscular disorder characterized by a waddling gait, limb deformities, muscular weakness and facial palsy. Exome sequencing was completed on the DNA of three of the four patients. We identified a novel missense variant in DCAF13, ENST00000612750.5, NM_015420.7, c.907 G > A;p.(Asp303Asn), ENST00000616836.4, NM_015420.6, c.1363 G > A:p.(Asp455Asn) (rs1209794872) segregating with this phenotype; being homozygous in all four affected patients and heterozygous in the unaffected individuals. The variant was extremely rare in the public databases (gnomAD allele frequency 0.000007081); was absent from the DNA of 300 ethnically matched controls and affected an amino acid which has been conserved across 1-2 billion years of evolution in eukaryotes. DCAF13 contains three WD40 domains and is hypothesized to have roles in both rRNA processing and in ubiquitination of proteins. Analysis of DCAF13 with the p.(Asp455Asn) variant predicted that the amino acid change is deleterious and affects a β-hairpin turn, within a WD40 domain of the protein which may decrease protein stability. Previously, a heterozygous variant of DCAF13 NM_015420.6, c.20 G > C:p.(Trp7Ser) with or without a heterozygous missense variant in CCN3, was suggested to cause inherited cortical myoclonic tremor with epilepsy. In addition, a heterozygous DCAF13 variant has been associated with autism spectrum disorder. Our study indicates a potential role of biallelic DCAF13 variants in neuromuscular disorders. Screening of additional patients with similar phenotype may broaden the allelic and phenotypic spectrum due to DCAF13 variants.

Identifying promising GSK3β inhibitors for cancer management: a computational pipeline combining virtual screening and molecular dynamics simulations

Glycogen synthase kinase-3 (GSK3β), a serine/threonine protein kinase, has been discovered as a novel target for anticancer drugs. Although GSK3β is involved in multiple pathways linked to the etiology of various cancers, no specific GSK3β inhibitor has been authorized for cancer therapy. Most of its inhibitors have toxicity effects therefore, there is a need to develop safe and more potent inhibitors. In this study, a library of 4,222 anti-cancer compounds underwent rigorous computational screening to identify potential candidates for targeting the binding pocket of GSK3β. The screening process involved various stages, including docking-based virtual screening, physicochemical and ADMET analysis, and molecular dynamics simulations. Ultimately, two hit compounds, BMS-754807 and GSK429286A, were identified as having high binding affinities to GSK3β. BMS-754807 and GSK429286A exhibited binding affinities of -11.9, and -9.8 kcal/mol, respectively, which were greater than that of the positive control (-7.6 kcal/mol). Further, molecular dynamics simulations for 100 ns were employed to optimize the interaction between the compounds and GSK3β, and the simulations demonstrated that the interaction was stable and consistent throughout the study. These hits were also anticipated to have good drug-like properties. Finally, this study suggests that BMS-754807 and GSK429286A may undergo experimental validation to evaluate their potential as cancer treatments in clinical settings.

Multi-population genome-wide association study implicates immune and non-immune factors in pediatric steroid-sensitive nephrotic syndrome

Pediatric steroid-sensitive nephrotic syndrome (pSSNS) is the most common childhood glomerular disease. Previous genome-wide association studies (GWAS) identified a risk locus in the HLA Class II region and three additional independent risk loci. But the genetic architecture of pSSNS, and its genetically driven pathobiology, is largely unknown. Here, we conduct a multi-population GWAS meta-analysis in 38,463 participants (2440 cases). We then conduct conditional analyses and population specific GWAS. We discover twelve significant associations-eight from the multi-population meta-analysis (four novel), two from the multi-population conditional analysis (one novel), and two additional novel loci from the European meta-analysis. Fine-mapping implicates specific amino acid haplotypes in HLA-DQA1 and HLA-DQB1 driving the HLA Class II risk locus. Non-HLA loci colocalize with eQTLs of monocytes and numerous T-cell subsets in independent datasets. Colocalization with kidney eQTLs is lacking but overlap with kidney cell open chromatin suggests an uncharacterized disease mechanism in kidney cells. A polygenic risk score (PRS) associates with earlier disease onset. Altogether, these discoveries expand our knowledge of pSSNS genetic architecture across populations and provide cell-specific insights into its molecular drivers. Evaluating these associations in additional cohorts will refine our understanding of population specificity, heterogeneity, and clinical and molecular associations.

Exploring biogenic chalcones as DprE1 inhibitors for antitubercular activity via in silico approach

Cases of drug-resistant tuberculosis (TB) have increased worldwide in the last few years, and it is a major threat to global TB control strategies and the human population. Mycobacterium tuberculosis is a common causative agent responsible for increasing cases of TB and as reported by WHO, approximately, 1.5 million death occurred from TB in 2020. Identification of new therapies against drug-resistant TB is an urgent need to be considered primarily. The current investigation aims to find the potential biogenic chalcone against the potential targets of drug-resistant TB via in silico approach. The ligand library of biogenic chalcones was screened against DprE1. Results of molecular docking and in silico ADMET prediction revealed that ZINC000005158606 has lead-like properties against the targeted protein. Pharmacophore modeling was done to identify the pharmacophoric features and their geometric distance present in ZINC000005158606. The binding stability study performed using molecular dynamics (MD) simulation of the DprE1-ZINC000005158606 complex revealed the conformational stability of the complex system over 100 ns with minimum deviation. Further, the in silico anti-TB sensitivity of ZINC000005158606 was found to be higher as compared to the standards against Mycobacterium tuberculosis. The overall in silico investigation indicated the potential of identified hit to act as a lead molecule against Mycobacterium tuberculosis.

In Silico Screening and Molecular Dynamics Simulation Studies in the Identification of Natural Compound Inhibitors Targeting the Human Norovirus RdRp Protein to Fight Gastroenteritis

Norovirus (HNoV) is a leading cause of gastroenteritis globally, and there are currently no treatment options or vaccines available to combat it. RNA-dependent RNA polymerase (RdRp), one of the viral proteins that direct viral replication, is a feasible target for therapeutic development. Despite the discovery of a small number of HNoV RdRp inhibitors, the majority of them have been found to possess a little effect on viral replication, owing to low cell penetrability and drug-likeness. Therefore, antiviral agents that target RdRp are in high demand. For this purpose, we used in silico screening of a library of 473 natural compounds targeting the RdRp active site. The top two compounds, ZINC66112069 and ZINC69481850, were chosen based on their binding energy (BE), physicochemical and drug-likeness properties, and molecular interactions. ZINC66112069 and ZINC69481850 interacted with key residues of RdRp with BEs of -9.7, and -9.4 kcal/mol, respectively, while the positive control had a BE of -9.0 kcal/mol with RdRp. In addition, hits interacted with key residues of RdRp and shared several residues with the PPNDS, the positive control. Furthermore, the docked complexes showed good stability during the molecular dynamic simulation of 100 ns. ZINC66112069 and ZINC69481850 could be proven as potential inhibitors of the HNoV RdRp in future antiviral medication development investigations.

Computational Exploration of Anti-cancer Potential of Flavonoids against Cyclin-Dependent Kinase 8: An In Silico Molecular Docking and Dynamic Approach

Over the centuries, cancer has been considered one of the significant health threats. It holds the position in the list of deadliest diseases over the globe. In women, breast cancer is the most common among many cancers and is the second most common cancer all over the world, while lung cancer is the first. Cyclin-dependent kinase 8 (CDK8) has been identified as a critical oncogenic driver that is found in breast cancer and associated with tumor progression. Flavonoids were virtually screened against CDK8 using molecular docking, drug-likeness, ADMET prediction, and a molecular dynamics (MD) simulation approach to determine the potential flavonoid structure against CDK8. The results indicated that ZINC000005854718 showed the highest negative binding affinity of -10.7 kcal/mol with the targeted protein and passed all the drug-likeness parameters. Performed molecular dynamics simulation showed that docked complex systems have good conformational stability over 100 ns in different temperatures (298, 300, 305, 310, and 320 K). The comparison between calculated binding free energy via MM/PB(GB)SA methods and binding affinity calculated via molecular docking suggested tight binding of ZINC000005854718 with targeted protein. The results concluded that ZINC000005854718 has drug-like properties with tight and stable binding with the targeted protein.

Identifying the natural compound Catechin from tropical mangrove plants as a potential lead candidate against 3CLpro from SARS-CoV-2: An integrated in silico approach

SARS-CoV-2, a member of beta coronaviruses, is a single-stranded, positive-sense RNA virus responsible for the COVID-19 pandemic. With global fatalities of the pandemic exceeding 4.57 million, it becomes crucial to identify effective therapeutics against the virus. A protease, 3CL^pro, is responsible for the proteolysis of viral polypeptides into functional proteins, which is essential for viral pathogenesis. This indispensable activity of 3CL^pro makes it an attractive target for inhibition studies. The current study aimed to identify potential lead molecules against 3CL^pro of SARS-CoV-2 using a manually curated in-house library of antiviral compounds from mangrove plants. This study employed the structure-based virtual screening technique to evaluate an in-house library of antiviral compounds against 3CL^pro of SARS-CoV-2. The library was comprised of thirty-three experimentally proven antiviral molecules extracted from different species of tropical mangrove plants. The molecules in the library were virtually screened using AutoDock Vina, and subsequently, the top five promising 3CL^pro-ligand complexes along with 3CL^pro-N3 (control molecule) complex were subjected to MD simulations to comprehend their dynamic behaviour and structural stabilities. Finally, the MM/PBSA approach was used to calculate the binding free energies of 3CLpro complexes. Among all the studied compounds, Catechin achieved the most significant binding free energy (-40.3 ± 3.1 kcal/mol), and was closest to the control molecule (-42.8 ± 5.1 kcal/mol), and its complex with 3CL^pro exhibited the highest structural stability. Through extensive computational investigations, we propose Catechin as a potential therapeutic agent against SARS-CoV-2. Communicated by Ramaswamy H. Sarma.

Dipeptidyl peptidase-4 inhibitory potentials of Glycyrrhiza uralensis and its bioactive compounds licochalcone A and licochalcone B: An in silico and in vitro study

Type 2 diabetes mellitus (T2DM) is a growing global public health issue, and dipeptidyl peptidase-4 (DPP-4) is a potential therapeutic target in T2DM. Several synthetic anti-DPP-4 medications can be used to treat T2DM. However, because of adverse effects, there is an unmet demand for the development of safe and effective medications. Natural medicines are receiving greater interest due to the inherent safety of natural compounds. Glycyrrhiza uralensis (licorice) is widely consumed and used as medicine. In this study, we investigated the abilities of a crude water extract (CWE) of G. uralensis and two of its constituents (licochalcone A (LicA) and licochalcone B (LicB)) to inhibit the enzymatic activity of DPP-4 in silico and in vitro. In silico studies showed that LicA and LicB bind tightly to the catalytic site of DPP-4 and have 11 amino acid residue interactions in common with the control inhibitor sitagliptin. Protein-protein interactions studies of LicA-DPP4 and LicB-DPP4 complexes with GLP1 and GIP reduced the DPP-4 to GLP1 and GIP interactions, indicated that these constituents might reduce the degradations of GLP1 and GIP. In addition, molecular dynamics simulations revealed that LicA and LicB stably bound to DPP-4 enzyme. Furthermore, DPP-4 enzyme assay showed the CWE of G. uralensis, LicA, and LicB concentration-dependently inhibited DPP-4; LicA and LicB had an estimated IC50 values of 347.93 and 797.84 μM, respectively. LicA and LicB inhibited DPP-4 at high concentrations, suggesting that these compounds could be used as functional food ingredients to manage T2DM.

Identification of compelling inhibitors of human norovirus 3CL protease to combat gastroenteritis: A structure-based virtual screening and molecular dynamics study

Human noroviruses (NV) are the most prevalent cause of sporadic and pandemic acute gastroenteritis. NV infections cause substantial morbidity and death globally, especially amongst the aged, immunocompromised individuals, and children. There are presently no authorized NV vaccines, small-molecule therapies, or prophylactics for humans. NV 3 C L protease (3CLP) has been identified as a promising therapeutic target for anti-NV drug development. Herein, we employed a structure-based virtual screening method to screen a library of 700 antiviral compounds against the active site residues of 3CLP. We report three compounds, Sorafenib, YM201636, and LDC4297, that were revealed to have a higher binding energy (BE) value with 3CLP than the control (Dipeptidyl inhibitor 7) following a sequential screening, in-depth molecular docking and visualization, physicochemical and pharmacological property analysis, and molecular dynamics (MD) study. Sorafenib, YM201636, and LDC4297 had BEs of -11.67, -10.34, and -9.78 kcal/mol with 3CLP, respectively, while control had a BE of -6.38 kcal/mol. Furthermore, MD simulations of the two best compounds and control were used to further optimize the interactions, and a 100 ns MD simulation revealed that they form stable complexes with 3CLP. The estimated physicochemical, drug-like, and ADMET properties of these hits suggest that they might be employed as 3CLP inhibitors in the management of gastroenteritis. However, wet lab tests are a prerequisite to optimize them as NV 3CLP inhibitors.

Identification of a novel farnesoid X receptor agonist, kaempferol-7-O-rhamnoside, a compound ameliorating drug-induced liver injury based on virtual screening and in vitro validation

Farnesoid X receptor (FXR), a bile acid receptor, plays an essential role in maintaining bile acid and liver homeostasis and has been recognized as an essential target for drug-induced liver injury (DILI). This study aimed to identify potential FXR agonists by virtual screening, molecular dynamics (MD) simulation, and biological assays. First, an in-house Traditional Chinese medicine compound database was screened using a virtual approach based on molecular docking to reveal potential FXR agonists. Secondly, MD was applied to analyze the process of agonist binding. Finally, the acetaminophen (APAP)-induced L02 cells model evaluated the pharmacodynamic activity of agonists treating DILI. Virtual screening results showed that kaempferol-7-O-rhamnoside was confirmed as the FXR agonist. MD results showed that kaempferol-7-O-rhamnoside could stably bind the FXR. In addition, in vitro cell-based assay showed that kaempferol-7-O-rhamnoside could promote the expression of the FXR gene and inhibit the Cyp7a1 gene expression in APAP-induced cells, significantly reducing the activities of AST, AKP and ROS, and enhancing the expression of GSH. The current study confirmed that kaempferol-7-O-rhamnoside might improve liver function by promoting proliferation, ameliorating oxidative stress, and regulating FXR target genes as observed in vitro. Therefore, in this study, discovering the FXR agonist, kaempferol-7-O-rhamnoside, provides valuable guidance for developing novel drugs against DILI.

Targeting Glutaminase by Natural Compounds: Structure-Based Virtual Screening and Molecular Dynamics Simulation Approach to Suppress Cancer Progression

Cancer cells change their glucose and glutamine (GLU) metabolism to obtain the energy required to continue growing. Glutaminase (GLS) plays a crucial role in promoting cell metabolism for cancer cell growth; targeting GLU metabolism by inhibiting GLS has attracted interest as a potential cancer management strategy. Herein, we employed a sequential screening of traditional Chinese medicine (TCM) database followed by drug-likeness and molecular dynamics simulations against the active site of GLS. We report 12 potent compounds after screening the TCM database against GLS, followed by a drug-likeness filter with Lipinski and Veber rule criteria. Among them, ZINC03978829 and ZINC32296657 were found to have higher binding energy (BE) values than the control compound 6-Diazo-5-Oxo-L-Norleucine, with BEs of -9.3 and -9.7 kcal/mol, respectively, compared to the BE of 6-Diazo-5-Oxo-L-Norleucine (-4.7 kcal/mol) with GLS. Molecular dynamics simulations were used to evaluate the results further, and a 100 ns MD simulation revealed that the hits form stable complexes with GLS and formed 2-5 hydrogen bond interactions. This study indicates that these hits might be employed as GLS inhibitors in the battle against cancer. However, more laboratory tests are a prerequisite to optimize them as GLS inhibitors.

Multifaceted Computational Modeling in Glycoscience

Glycoscience assembles all the scientific disciplines involved in studying various molecules and macromolecules containing carbohydrates and complex glycans. Such an ensemble involves one of the most extensive sets of molecules in quantity and occurrence since they occur in all microorganisms and higher organisms. Once the compositions and sequences of these molecules are established, the determination of their three-dimensional structural and dynamical features is a step toward understanding the molecular basis underlying their properties and functions. The range of the relevant computational methods capable of addressing such issues is anchored by the specificity of stereoelectronic effects from quantum chemistry to mesoscale modeling throughout molecular dynamics and mechanics and coarse-grained and docking calculations. The Review leads the reader through the detailed presentations of the applications of computational modeling. The illustrations cover carbohydrate-carbohydrate interactions, glycolipids, and N- and O-linked glycans, emphasizing their role in SARS-CoV-2. The presentation continues with the structure of polysaccharides in solution and solid-state and lipopolysaccharides in membranes. The full range of protein-carbohydrate interactions is presented, as exemplified by carbohydrate-active enzymes, transporters, lectins, antibodies, and glycosaminoglycan binding proteins. A final section features a list of 150 tools and databases to help address the many issues of structural glycobioinformatics.

Identification and Evaluation of Traditional Chinese Medicine Natural Compounds as Potential Myostatin Inhibitors: An In Silico Approach

Myostatin (MSTN), a negative regulator of muscle mass, is reported to be increased in conditions linked with muscle atrophy, sarcopenia, and other muscle-related diseases. Most pharmacologic approaches that treat muscle disorders are ineffective, emphasizing the emergence of MSTN inhibition. In this study, we used computational screening to uncover natural small bioactive inhibitors from the Traditional Chinese Medicine database (~38,000 compounds) for the MSTN protein. Potential ligands were screened, based on binding affinity (150), physicochemical (53) and ADMET properties (17). We found two hits (ZINC85592908 and ZINC85511481) with high binding affinity and specificity, and their binding patterns with MSTN protein. In addition, molecular dynamic simulations were run on each complex to better understand the interaction mechanism of MSTN with the control (curcumin) and the hit compounds (ZINC85592908 and ZINC85511481). We determined that the hits bind to the active pocket site (Helix region) and trigger conformational changes in the MSTN protein. Since the stability of the ZINC85592908 compound was greater than the MSTN control, we believe that ZINC85592908 has therapeutic potential against the MSTN protein and may hinder downstream singling by inhibiting the MSTN protein and increasing myogenesis in the skeletal muscle tissues.

Heparin-Assisted Amyloidogenesis Uncovered through Molecular Dynamics Simulations

Glycosaminoglycans (GAGs), in particular, heparan sulfate and heparin, are found colocalized with Aβ amyloid. They have been shown to enhance fibril formation, suggesting a possible pathological connection. We have investigated heparin's assembly of the KLVFFA peptide fragment using molecular dynamics simulation, to gain a molecular-level mechanistic understanding of how GAGs enhance fibril formation. The simulations reveal an exquisite process wherein heparin accelerates peptide assembly by first "gathering" the peptide molecules and then assembling them. Heparin does not act as a mere template but is tightly coupled to the peptides, yielding a composite protofilament structure. The strong intermolecular interactions suggest composite formation to be a general feature of heparin's interaction with peptides. Heparin's chain flexibility is found to be essential to its fibril promotion activity, and the need for optimal heparin chain length and concentration has been rationalized. These insights yield design rules (flexibility; chain-length) and protocol guidance (heparin:peptide molar ratio) for developing effective heparin mimetics and other functional GAGs.

Structural and allergenic properties of the fatty acid binding protein from shrimp Litopenaeus vannamei

BACKGROUND

The shrimp Litopenaeus vannamei is an important source of food allergens but its allergenic repertoire is poorly characterized. Cross-reactivity between crustacean and mites has been reported, with tropomyosin, the most relevant allergen involved. The aim of this study was to investigate the structural and immunological properties of a recombinant Fatty Acid Binding Protein (FABP) family from L. vannamei (LvFABP).

METHODS

ELISA, skin prick test (SPT) and basophil activation assays were performed to determine IgE reactivity and allergenic activity of LvFABP. LC-MS/MS and Circular Dichroism experiments were done for structural analysis. B-cell epitope mapping with overlapping peptides, and cross-inhibition studies using human sera were done to identify antigenic regions and cross-reactivity.

RESULTS

The recombinant LvFABP bound serum IgE from 27% of 36 shrimp allergic patients and showed allergenic activity when tested for basophil activation and SPT in a selected number of them. CD-spectroscopy of LvFABP revealed that the protein is folded with a secondary structure composed of mainly β-strands and a smaller fraction of α helices. This is consistent with molecular modelling results, which exhibit a typical β barrel fold with two α-helices and ten β-strands. Epitope mapping identified two IgE-binding antigenic regions and inhibition assays found high cross-reactivity between LvFABP and Blo t 13, mediated by the antigenic region involving amino acids 54 to 72.

CONCLUSIONS

Our results show that LvFABP is a shrimp allergen that cross reacts with the house dust mite allergen Blo t 13 and has allergenic activity, which suggest that it could be clinically relevant in case of shellfish allergy. This new allergen, named Lit v 13, will also help to understand basic mechanisms of sensitization to shrimp.

In silico ADMET, molecular docking and molecular simulation-based study of glabridin’s natural and semisynthetic derivatives as potential tyrosinase inhibitors

Hyper-pigmentation conditions may develop due to erroneous melanogenesis cascade which leads to excess melanin production. Recently, inhibition of tyrosinase is the main focus of investigation as it majorly contributes to melanin production. This inhibition property can be exploited in medicine, agriculture, and in cosmetics. Present study aims to find a natural and safe alternative molecule as tyrosinase inhibitor. In this study, human tyrosinase enzyme was modelled due to unavailability of its crystal structure to look into the degree of efficacy of glabridin and its 15 derivatives as tyrosinase inhibitor. Docking was performed by Autodock Vina at the catalytic core enzyme. Glabridin effects on melanoma cell lines was also elucidated by analysing cytotoxicity and effect on melanin production. Computational ADME analysis was done by SwissADME. Molecular dynamic simulation was also performed to further evaluate the interaction profile of these molecules and kojic acid (positive inhibitor) with respect to apo protein. Notably, four derivatives 5′-formylglabridin, glabridin dimer, 5′-prenyl glabridin and R-glabridin exhibited better binding affinity than glabridin. Glabridin effectively inhibited melanin production in a dose dependent manner. Among these, 5′-formylglabridin displayed highest binding affinity with docking score − 9.2 kcal/mol. Molecular properties and bioactivity analysis by Molinspiration web server and by SwissADME also presented these molecules as potential drug candidates. The study explores the understanding for the development of suitable tyrosinase inhibitor/s for the prevention of hyperpigmentation. However, a detailed in vivo study is required for glabridin derivatives to suggest these molecules as anti-melanogenic compound.

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.

Targeting Cytotoxin-Associated Antigen A, a Virulent Factor of Helicobacter pylori-Associated Gastric Cancer: Structure-Based In Silico Screening of Natural Compounds

Gastric cancer is the fifth most frequent cancer and the third major cause of mortality worldwide. Helicobacter pylori, a bacterial infection linked with GC, injects the cytotoxin-associated antigen A (CagA; an oncoprotein) into host cells. When the phosphorylated CagA protein enters the cell, it attaches to other cellular components, interfering with normal cellular signaling pathways. CagA plays an important role in the progression of GC by interacting with phosphatidylserine of the host cell membrane. Therefore, disrupting the CagA-phosphatidylserine connection using small molecules appears to be a promising therapeutic approach. In this report, we screened the natural compounds from ZINC database against the CagA protein using the bioinformatics tools. Hits were initially chosen based on their physicochemical, absorption, distribution, metabolism, excretion, and toxicity (ADMET) characteristics, as well as other drug-like characteristics. To locate safe and effective hits, the PAINS filter, binding affinities estimation, and interaction analysis were used. Three compounds with high binding affinity and specificity for the CagA binding pocket were discovered. The final hits, ZINC153731, ZINC69482055, and ZINC164387, were found to bind strongly with CagA protein, with binding energies of -11.53, -10.67, and -9.21 kcal/mol, respectively, which were higher than that of the control compound (-7.25 kcal/mol). Further, based on binding affinity and interaction pattern, two leads (ZINC153731, ZINC69482055) were chosen for molecular dynamics (MD) simulation analysis. MD results showed that they displayed stability in their vicinity at 100 ns. This study suggested that these compounds could be used as possible inhibitors of CagA protein in the fight against GC. However, additional benchwork tests are required to validate them as CagA protein inhibitors.

Holistic review of polysaccharides isolated from pumpkin: Preparation methods, structures and bioactivities

Pumpkin polysaccharides have arrested researchers' attention in fields of food supplements for healthy product and traditional Chinese medicine due to their multiple bioactivities with non-toxic and highly biocompatible. This review emphatically summarized recent progresses in the primary and spatial structural features, various bioactivities, structure-to-function associations, different preparation techniques, and absorption characteristics across intestinal epithelial and in vivo bio-distribution of pumpkin polysaccharides. Additionally, current challenges and future trends in development of pumpkin polysaccharides were pointed out. We found that pumpkin polysaccharides were primary structure (e.g. glucan, galactoglucan, galactomannan, galactan, homogalacturonan (HG), and rhamnogalacturonan-Ι (RG-Ι)) and special structure diverse (e.g. hollow helix, linear, and sphere-like) and significant functional foods or therapeutic agents (e.g. oral hypoglycemic agents). Moreover, we found that the molecular weight (Mw), uronic acid, linkage types, and modifications all could affect their bioactivities (e.g. anti-oxidant, anti-coagulant, and anti-diabetic activities), and pumpkin polysaccharides may across intestinal epithelial into the blood reaching to target organs. Collectively, the structures diversity and pharmacological values of pumpkin polysaccharides support their therapeutic potentials and sanitarian functions.

Identification of Potent Natural Resource Small Molecule Inhibitor to Control Vibrio cholera by Targeting Its Outer Membrane Protein U: An In Silico Approach

Vibrio cholerae causes the diarrheal disease cholera which affects millions of people globally. The outer membrane protein U (OmpU) is the outer membrane protein that is most prevalent in V. cholerae and has already been recognized as a critical component of pathogenicity involved in host cell contact and as being necessary for the survival of pathogenic V. cholerae in the host body. Computational approaches were used in this study to screen a total of 37,709 natural compounds from the traditional Chinese medicine (TCM) database against the active site of OmpU. Following a sequential screening of the TCM database, we report three lead compounds-ZINC06494587, ZINC85510056, and ZINC95910434-that bind strongly to OmpU, with binding affinity values of -8.92, -8.12, and -8.78 kcal/mol, which were higher than the control ligand (-7.0 kcal/mol). To optimize the interaction, several 100 ns molecular dynamics simulations were performed, and the resulting complexes were shown to be stable in their vicinity. Additionally, these compounds were predicted to have good drug-like properties based on physicochemical properties and ADMET assessments. This study suggests that further research be conducted on these compounds to determine their potential use as cholera disease treatment.

In silico study of some selective phytochemicals against a hypothetical SARS-CoV-2 spike RBD using molecular docking tools

Background

This world is currently witnessing a pandemic outbreak of ‘COVID-19′ caused by a positive-strand RNA virus ‘SARS-CoV-2’. Millions have succumbed globally to the disease, and the numbers are increasing day by day. The viral genome enters into the human host through interaction between the spike protein (S) and host angiotensin-converting enzyme-2 (ACE2) proteins. S is the common target for most recently rolled-out vaccines across regions. A recent surge in single/multiple mutations in S region is of great concern as it may escape vaccine induced immunity. So far, the treatment regime with repurposed drugs has not been too successful.

Hypothesis

Natural compounds are capable of targeting mutated spike protein by binding to its active site and destabilizing the spike-host ACE2 interaction.

Materials and methods

A hypothetical mutated spike protein was constructed by incorporating twelve different mutations from twelve geographical locations simultaneously into the receptor-binding domain (RBD) and docked with ACE2 and seven phytochemicals namely allicin, capsaicin, cinnamaldehyde, curcumin, gingerol, piperine and zingeberene. Molecular Dynamic (MD) simulation and Principal Component Analysis (PCA) were finally used for validation of the docking results.

Result

The docking results showed that curcumin and piperine were most potent to bind ACE2, mutated spike, and mutated spike-ACE2 complex, thereby restricting viral entry. ADME analysis also proved their drug candidature. The docking complexes were found to be stable by MD simulation.

Conclusion

This result provides a significant insight about the phytochemicals' role, namely curcumin and piperine, as the potential therapeutic entities against mutated spike protein of SARS-CoV-2.

The Underlying Molecular Mechanisms Involved in Traditional Chinese Medicine Smilax china L. for the Treatment of Pelvic Inflammatory Disease

Smilax china L. (SCL) is extensively used in the treatment of pelvic inflammatory disease (PID). This study aimed to clarify the potential active ingredients of SCL and mechanisms on PID. SCL was widely distributed in Japan, South Korea, and China, which was traditionally considered heat-clearing, detoxicating, and dampness-eliminating medicine. Systems pharmacology revealed that 32 compounds in SCL may interact with 19 targets for immunoenhancement, antiapoptosis, anti-inflammation, and antioxidant activity of the PID model. Molecular docking revealed that isorhamnetin, moracin M, rutin, and oxyresveratrol may have higher binding potential with prostaglandin-endoperoxide synthase 2 (PTGS2), mitogen-activated protein kinase 1 (MAPK1), siderocalin (LCN2), tumor necrosis factor (TNF), and matrix metalloprotein-9 (MMP9), respectively. Molecular dynamics simulation showed that the binding modes of moracin M-MAPK1, rutin-TNF, and oxyresveratrol-MMP9 complexes were more stable, evidenced by relatively smaller fluctuations in root mean square deviation values. Conclusively, SCL may treat PID by inhibiting inflammatory factors, antitissue fibrosis, and microbial growth.

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.

Three-Dimensional Structures of Carbohydrates and Where to Find Them

Analysis and systematization of accumulated data on carbohydrate structural diversity is a subject of great interest for structural glycobiology. Despite being a challenging task, development of computational methods for efficient treatment and management of spatial (3D) structural features of carbohydrates breaks new ground in modern glycoscience. This review is dedicated to approaches of chemo- and glyco-informatics towards 3D structural data generation, deposition and processing in regard to carbohydrates and their derivatives. Databases, molecular modeling and experimental data validation services, and structure visualization facilities developed for last five years are reviewed.

The Lazy Life of Lipid-Linked Oligosaccharides in All Life Domains

Lipid-linked oligosaccharides (LLOs) play an important role in the N-glycosylation pathway as the donor substrate of oligosaccharyltransferases (OSTs), which are responsible for the en bloc transfer of glycan chains onto a nascent polypeptide. The lipid component of LLO in both eukarya and archaea consists of a dolichol, and an undecaprenol in prokarya, whereas the number of isoprene units may change between species. Given the potential relevance of LLOs and their related enzymes to diverse biotechnological applications, obtaining reliable LLO models from distinct domains of life could support further studies on complex formation and their processing by OSTs, as well as protein engineering on such systems. In this work, molecular modeling techniques, such as quantum mechanics calculations, molecular dynamics simulations, and metadynamics were employed to study eukaryotic (Glc₃-Man₉-GlcNAc₂-PP-Dolichol), bacterial (Glc₁-GalNAc₅-Bac₁-PP-Undecaprenol), and archaeal (Glc₁-Man1-Gal₁-Man1-Glc₁-Gal₁-Glc₁-P-Dolichol) LLOs in membrane bilayers. Microsecond molecular dynamics simulations and metadynamics calculations of LLOs revealed that glycan chains are more prone to interact with the membrane lipid head groups, while the PP linkages are positioned at the lipid phosphate head groups level. The dynamics of isoprenoid chains embedded within the bilayer are described, and membrane dynamics and related properties are also investigated. Overall, there are similarities regarding the structure and dynamics of the eukaryotic, the bacterial, and the archaeal LLOs in bilayers, which can support the comprehension of their association with OSTs. These data may support future studies on the transferring mechanism of the oligosaccharide chain to an acceptor protein.

Current Status of Carbohydrates Information in the Protein Data Bank

Carbohydrates are well known for their physicochemical, biological, functional, and therapeutic characteristics. Unfortunately, their chemical nature imposes severe challenges for the structural elucidation of these phenomena, impairing not only the depth of our understanding of carbohydrates but also the development of new biotechnological and therapeutic applications based on these molecules. In the recent past, the amount of structural information, obtained mainly from X-ray crystallography, has increased progressively, as well as its quality. In this context, the current work presents a global analysis of the carbohydrate information available in the Protein Data Bank (PDB). From high quality structures, it is clear that most of the data are highly concentrated on a few sets of residue types, on their monosaccharidic forms, and connected by a small diversity of glycosidic linkages. The geometries of these linkages can be mostly associated with the types of linkages instead of residues, while the level of puckering distortion was characterized, quantified, and located in a pseudorotational equilibrium landscape, not only to local minima but also to transitional states. These qualitative and quantitative analyses offer a global picture of the carbohydrate structural content in the PDB, potentially supporting the building of new models for carbohydrate-related biological phenomena at the atomistic level, including new developments on force field parameters.

Discovery of selective farnesoid X receptor agonists for the treatment of hyperlipidemia from traditional Chinese medicine based on virtual screening and in vitro validation

Farnesoid X receptor (FXR), a bile acid receptor, has important roles in maintaining bile acid and cholesterol homeostasis, which is an attractive target for hyperlipidemia. Present study aimed to discover potential selective FXR agonists over G-protein coupled bile acid receptor 1 (GPBAR1, TGR5) from traditional Chinese medicine (TCM) by using virtual screening, in vitro studies and molecular dynamics simulation (MD). Ligand-based pharmacophore model for FXR was firstly built to screen FXR agonists from the Traditional Chinese Medicine Database (TCMD). Then, 21 FXR crystal structures were clustered in two types and two representative structures (PDB ID: 3OMM and 3P89) were, respectively, used to carry out molecular docking to refine the screened result. Moreover, the pharmacophore model for GPBAR1 was built to screen selective FXR agonists with no activity on GPBAR1. A set of 24 candidate selective FXR agonists which fitvalue of FXR pharmacophore model and docking score of 3OMM and 3P89 were in the top 100 and cannot match the pharmacophore model for GPBAR1 were obtained. By the lipid-lowering activity test in HepG2 cell lines, Arctigenin was identified to be potential selective FXR agonist with the activity of 20 μmol·L^-1. After down-regulating FXR, Arctigenin could increase the mRNA of FXR while exerted no effect on the mRNA of GPBAR1. MD was further used to interpret the mechanism of Arctigenin with the representative structures. This research provided a new screening procedure for finding selective candidate compounds and appropriate docking models of a target by considering the structure diversity of PDB structures, which was applied to discovery novel selective FXR agonists to treat hyperlipidemia.Communicated by Ramaswamy H. Sarma.

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.