Search

Abstract Background Breast cancer, one of the most often diagnosed malignancies worldwide, continues to take countless women's lives. Its treatment usually involves targeting the human estrogen receptor alpha (ERα). Current research explores the potential of natural compounds to regulate ERα activity, providing a hopeful direction for breast cancer therapy. Our study utilized a comprehensive approach to identify promising natural compounds for breast cancer treatment, including quantum descriptors, molecular docking, molecular dynamics simulations, and ADMET/pharmacokinetics analysis. Results Six natural compounds derived from podophyllum medicinal plants, namely 4-demethylpodophyllotoxin (NP1), α-peltatin (NP2), podophyllotoxin (NP3), deoxypodophyllotoxin (NP4), podophyllotoxone (NP5), and β-peltatin (NP6), were investigated as potential selective estrogen receptor α (ERα) inhibiting agents for breast cancer. These compounds demonstrated the strongest binding affinity to the target enzyme, with binding energies of − 8.9 and − 8.1 kcal/mol, respectively. Further assessments of drug-likeness and ADME properties were conducted for these compounds, along with quantum calculations (HOMO–LUMO) to evaluate their reactivity. Additionally, molecular dynamics studies were performed to assess the stability of the NP1 and NP2 protein–ligand complexes. Conclusions We analyzed six natural compounds comprehensively, evaluating their ADME properties, molecular docking interactions, quantum descriptors, and dynamic simulations. Our findings demonstrate that these natural compounds are promising possibilities for treating breast cancer. Additionally, they may provide a basis for developing future compounds targeting estrogen receptor α (ERα) activity. Graphical abstract

The purpose of this study is to investigate the molecular interactions and potential therapeutic uses of Eltrombopag (EPAG), a small molecule that activates the cMPL receptor. EPAG has been found to be effective in increasing platelet levels and alleviating thrombocytopenia. We utilized computational techniques to predict and confirm the complex formed by the ligand (EPAG) and the Thrombopoietin receptor (TPO-R) cMPL, elucidating the role of RAS, JAK-2, STAT-3, and other essential elements for downstream signaling. Molecular dynamics (MD) simulations were employed to evaluate the stability of the ligand across specific proteins, showing favorable characteristics. For the first time, we examined the presence of TPO-R in human umbilical cord mesenchymal stem cells (hUCMSC) and human gingival mesenchymal stem cells (hGMSC) proliferation. Furthermore, treatment with EPAG demonstrated angiogenesis and vasculature formation of endothelial lineage derived from both MSCs. It also indicated the activation of critical factors such as RUNX-1, GFI-1b, VEGF-A, MYB, GOF-1, and FLI-1. Additional experiments confirmed that EPAG could be an ideal molecule for protecting against UVB radiation damage, as gene expression (JAK-2, ERK-2, MCL-1, NFkB, and STAT-3) and protein CD90/cMPL analysis showed TPO-R activation in both hUCMSC and hGMSC. Overall, EPAG exhibits significant potential in treating radiation damage and mitigating the side effects of radiotherapy, warranting further clinical exploration.

Allantoin represents a compound widely employed in pharmaceutical and cosmetic fields. Its safety has been acknowledged by regulatory bodies such as the US Food and Drug Administration, the European Commission for Cosmetics and Consumer and Health and European Directorate for the Quality of Medicines & HealthCare. This justifies its wide use in dermatological/cosmetic formulations and allows their safe use.Allantoin possesses an asymmetric carbon atom, resulting in two enantiomers, with the (S)-enantiomer predominating in plants, although racemization may potentially occur during manufacturing processes. Notably, literature currently lacks enantioselective LC methods for allantoin analysis.In this study, two zwitterionic Cinchona alkaloid-based chiral stationary phases (CSPs), commercially known as CHIRALPAK® ZWIX(+) (CSP1) and CHIRALPAK® ZWIX(-) (CSP2), were utilized for the enantioseparation of allantoin under polar-ionic conditions. By employing a mobile phase consisting of acetonitrile/methanol/water/acetic acid (96:2:2:0.1, v/v/v/v), nearly complete baseline separation (with α=1.08) of allantoin enantiomers was achieved in less than 15 min with both CSPs. Due to the “pseudo-enantiomeric” nature of the two chiral selectors (quinine-based in CSP1 and quinidine-based in CSP2), an inversion of the enantiomer elution order was observed with the two CSPs under identical experimental conditions. Remarkably, this represents a rare instance where these CSPs demonstrate the ability to enantioseparate a non-ionic, non-ionizable species.The application of a molecular dynamics in silico protocol proved useful in elucidating the retention mechanism in depth, casting light on the central role of the H-bond formation and the involvement of the anionic moiety of the CSP 1.

This study aims to determine the antibacterial efficacy of a bioactive phytochemical quercetin (QUR), compared to the reference antibiotic ciprofloxacin (CIP), using in vitro and in silico approaches. Following molecular docking, QUR showed strong affinity to the druggable targets from bacteria: Sortase B, Toxic Shock Syndrome Toxin-1, Multidrug Efflux Pump AdeJ and LasR, with binding energies –6.9 to –10.3 kcal/mol. CIP displayed binding energies –7.3 to –7.4 kcal/mol against the target proteins. QUR had stronger efficacy against MEPJ (binding energy: -8.1 kcal/mol) and LasR (binding energy: -10.3 kcal/mol), compared to CIP having binding energies -7.4 and -7.3 kcal/mol, respectively. Molecular dynamics simulation validated the formation of LasR complex with QUR, which was energetically more stable compared to the LasR-CIP complex. Disc diffusion method, using nutrient agar medium, demonstrated antibacterial activity of QUR against both gram-negative (Acinetobacter baumannii, Escherichia coli, Pseudomonas aeruginosa) and gram-positive (Staphylococcus aureus, Streptococcus pneumonia) bacteria with zone diameter of inhibition 10–23 mm and 7–21 mm, respectively. QUR displayed synergistic effect, with growth inhibitory index (GII) values of 0.65–1.15, combined with PIP. QUR-CEF combination displayed synergy against E. coli, P. aeruginosa, and S. pneumoniae with GIIs 0.60–0.73. Pharmacokinetics studies for QUR revealed acceptable parametric values, comparable to that of CIP. Thus, the current integrated in silico and in vitro analysis authenticated QUR as a strong lead molecule to combat bacterial infections, using QUR alone or in combination with conventionally used antibiotics (PIP and CEF). For further validation of our results, clinical experiments are suggested.

Introduction: The most common primary brain tumor in adults is glioblastoma multiforme (GBM), accounting for 45.2% of all cases. The characteristics of GBM, a highly aggressive brain tumor, include rapid cell division and a propensity for necrosis. Regretfully, the prognosis is extremely poor, with only 5.5% of patients surviving after diagnosis.Methodology: To eradicate these kinds of complicated diseases, significant focus is placed on developing more effective drugs and pinpointing precise pharmacological targets. Finding appropriate biomarkers for drug discovery entails considering a variety of factors, including illness states, gene expression levels, and interactions between proteins. Using statistical techniques like p-values and false discovery rates, we identified differentially expressed genes (DEGs) as the first step in our research for identifying promising biomarkers in GBM. Of the 132 genes, 13 showed upregulation, and only 29 showed unique downregulation. No statistically significant changes in the expression of the remaining genes were observed.Results: Matrix metallopeptidase 9 (MMP9) had the greatest degree in the hub biomarker gene identification, followed by (periostin (POSTN) at 11 and Hes family BHLH transcription factor 5 (HES5) at 9. The significance of the identification of each hub biomarker gene in the initiation and advancement of glioblastoma multiforme was brought to light by the survival analysis. Many of these genes participate in signaling networks and function in extracellular areas, as demonstrated by the enrichment analysis.We also identified the transcription factors and kinases that control proteins in the proteinprotein interactions (PPIs) of the DEGs.Discussion: We discovered drugs connected to every hub biomarker. It is an appealing therapeutic target for inhibiting MMP9 involved in GBM. Molecular docking investigations indicated that the chosen complexes (carmustine, lomustine, marimastat, and temozolomide) had high binding affinities of −6.3, −7.4, −7.7, and −8.7 kcal/mol, respectively, the mean root-mean-square deviation (RMSD) value for the carmustine complex and marimastat complex was 4.2 Å and 4.9 Å, respectively, and the lomustine and temozolomide complex system showed an average RMSD of 1.2 Å and 1.6 Å, respectively. Additionally, high stability in root-mean-square fluctuation (RMSF) analysis was observed with no structural conformational changes among the atomic molecules. Thus, these in silico investigations develop a new way for experimentalists to target lethal diseases in future.

BackgroundSoil- and water-transmitted helminths are a major concern in the developing world due to their high prevalence. More than a quarter of the population were estimated to be infected with helminths in these endemic zones.Research designAn in silico approach was used to design a vaccine construct against the Taenia genus utilizing the proteomic information and evaluation of the construct using immune-informatics.ResultsOur study identified 451 conserved proteins in Taenia spp. using the existing proteome; out of these, 141 were found to be expressed in cysticerci. These proteins were screened for antigenic epitopes and a multi-subunit vaccine was constructed. The constructed vaccine was assessed for its efficacy in mounting the appropriate immune response. Our constructed vaccine showed stability and optimal performance against the TLR 4 receptor, which is reported to be upregulated in Taenia infections in hosts.ConclusionImmune-informatics tools help design vaccines for neglected diseases such as those attributed to helminths, which are known to cause widespread morbidity. Our vaccine construct holds tremendous potential in conferring protection against all Taenia spp. of clinical relevance to human.

Abstract Background Single nucleotide polymorphisms (SNPs) have been linked with prostate cancer (PCa) and have shown potential as prognostic markers for advanced stages. Loss of function mutations in PKCι have been linked with increased risk of malignancy by enhancing tumor cell motility and invasion. We have evaluated the impact of two coding region SNPs on the PKCι gene (PRKCI) and their prognostic potential. Methods Genotypic association of non-synonymous PKCι SNPs rs1197750201 and rs1199520604 with PCa was determined through tetra-ARMS PCR. PKCι was docked with interacting partner Par-6 to determine the effect of these variants on PKCι binding capabilities. Molecular dynamic simulations of PKCι docked with Par-6 were performed to determine variant effects on PKCι protein interactions. The possible impact of changes in PKCι protein interactions on epithelial cell polarity was hypothesized. Results PKCι rs1199520604 mutant genotype TT showed association with PCa (p = 0.0055), while rs1197750201 mutant genotype AA also showed significant association with PCa (P = 0.0006). The binding interaction of PKCι with Par-6 was altered for both variants, with changes in Van der Waals energy and electrostatic energy of docked structures. Conclusion Genotypic analysis of two non-synonymous PKCι variants in association with PCa prognosis was performed. Both variants in the PB1 domain showed potential as a prognostic marker for PCa. In silico analysis of the effect of the variants on PKCι protein interactions indicated they may be involved in PCa progression through aberration of epithelial cell polarity pathways.

Abstract Most of the breast cancers are estrogen receptor-positive recurring with a steady rate of up to 20 years dysregulating the normal cell cycle. Dinaciclib is still in clinical trials and considered as a research drug against such cancers targeting CDK2. The major goal of this study was to identify the potential inhibitors of CDK-2 present in Moringa oleifera for treating hormonal receptor positive breast cancers. For this purpose, in silico techniques; molecular docking, MM-GBSA and molecular dynamics simulations were employed to screen Moringa oleifera compounds and their anticancer potential was determined against CDK-2 protein targets. Among 36 compounds of Moringa oleifera reported in literature, chlorogenic acid (1), quercetin (2), ellagic acid (3), niazirin (4), and kaempferol (5) showed good affinity with the target. The interaction of the compounds was visualized using PYMOL software. The profiles of absorption, distribution, metabolism, excretion (ADME) and toxicity were determined using SWISS and ProTox II webservers. The MTT assay was performed in-vitro using MCF-7 cancer cell lines to validate the anticancer potential of Moringa oleifera leaf extract. MTT assay results revealed no significant change in proliferation of Mcf-7 cells following 24 h treatment with fraction A (petroleum ether). However, significant antiproliferative effect was observed at 200 µg/mL dose of fraction B (ethyl acetate) and cell viability was reduced to 40%. In conclusion, the data suggested that all the compounds with highest negative docking score than the reference could be the potential candidates for cyclin dependent kinase-2 (CDK-2) inhibition while ellagic acid, chlorogenic acid and quercetin being the most stable and potent inhibitors to treat estrogen receptor positive breast cancer targeting CDK-2. Moreover, the data suggested that further investigation is required to determine the optimum dose for significant antiproliferative effects using in-vivo models to validate our findings of in-silico analysis.

Abstract In this study, an environmentally friendly, solvent- and catalyst-free synthesis of 2-anilino nicotinic acids derivatives is reported. This operationally simple and green procedure was applied to a selection of primary aromatic amines giving rise to 23 derivatives of 2-anilino nicotinic acids in a very short reaction time (15–120 min) with good to excellent yield. Next, similarity searches were executed on these derivatives to find the possible biological target. These products were screened for inhibition of COX-1 and COX-2 by molecular docking and dynamic studies. In silico studies revealed that among these derivatives, the structure 10 bearing meta-chlorine substitutions could act as COX-1 and COX-2 inhibitors. These results can be used in designing important lead compounds for further development as potential anti-inflammatory drugs.

SHP2 belongs to a cytoplasmic non-receptor protein tyrosine phosphatase class. It plays a critical role in the development of various cancers, such as gastric cancer, leukemia, and breast cancer. Thus, SHP2 has gained the interest of researchers as a potential target for inhibiting tumor cell proliferation in SHP2-dependent cancers. This study employed pharmacophore-based virtual screening, molecular docking, molecular dynamic (MD) simulations, MM/PBSA, and principal component analysis (PCA), followed by ADME prediction. We selected three potential hits from a collective database of more than one million chemical compounds. The stability of these selected hit–protein complexes was analyzed using 500 ns MD simulations and binding free energy calculations. The identified hits Lig_1, Lig_6, and Lig_14 demonstrated binding free energies of −161.49 kJ/mol, −151.28 kJ/mol, and −107.13 kJ/mol, respectively, compared to the reference molecule (SHP099) with a ΔG of −71.48 kJ/mol. Our results showed that the identified compounds could be used as promising candidates for selective SHP2 allosteric inhibition in cancer.

Abstract Background COVID-19 was declared a pandemic by the World Health Organisation in 2020 after its outbreak in December 2019 in Wuhan, China. Since researchers have been working to develop specific drugs to cure COVID-19. COVID-19 is caused by the severe acute respiratory cornonavirus2 or popularly known as SARS-CoV2 attacking the ACE2 receptor in the human respiratory system. The main protease translated by the viral genome is a highly conserved protein that plays a crucial role in viral protein replication and transcription. Compounds such as Darunavir and danoprevir have been tested to show potential biological activity against the viral protein, but a high mutation rate defies a permanent solution to this problem. Results In this study, virtual screening of natural ligands (around 170,000 molecules) and FDA-approved repurposed drugs retrieved from ZINC Database was carried out against SARS-CoV2 main protease (PDB ID: 7DJR). Molecular coupling was performed for the top three ligands, where ZINC70699832 showed a significantly good binding affinity of − 11.05 kcal/mol. It has shown an interaction affinity for the residues THR-25, PHE-140, LEU-141, ASN-142, GLY-143, SER-144, CYS-145, MET-165, GLU-166, GLN-189 and GLN-192. The molecular dynamic simulation was also performed using GROMACS, for all complexes where the ZINC70699832–7DJR complex showed stability in terms of root mean square deviation. Conclusion The study recommends that ZINC70699832 has great potential to serve as a potent inhibitor of the main protease of SARS-CoV2 main protease.

Peptide-based liquid droplets (coacervates) produced by spontaneous liquid-liquid phase separation (LLPS), have emerged as a promising class of drug delivery systems due to their high entrapping efficiency and the simplicity of their formulation. However, the detailed mechanisms governing their interaction with cell membranes and cellular uptake remain poorly understood. In this study, we investigated the interactions of peptide coacervates composed of HBpep—peptide derived from the histidine-rich beak proteins (HBPs) of the Humboldt squid—with model cellular membranes in the form of supported lipid bilayers (SLBs). We employed quartz crystal microbalance with dissipation monitoring (QCM-D), neutron reflectometry (NR) and atomistic molecular dynamics (MD) simulations to reveal the nature of these interactions in the absence of fluorescent labels or tags. HBpep forms small oligomers at pH 6 whereas it forms µm-sized coacervates at physiological pH. Our findings reveal that both HBpep oligomers and HBpep-coacervates adsorb onto SLBs at pH 6 and 7.4, respectively. At pH 6, when the peptide carries a net positive charge, HBpep oligomers insert into the SLB, facilitated by the peptide’s interactions with the charged lipids and cholesterol. Importantly, however, HBpep coacervate adsorption at physiological pH, when it is largely uncharged, is fully reversible, suggesting no significant lipid bilayer rearrangement. HBpep coacervates, previously identified as efficient drug delivery vehicles, do not interact with the lipid membrane in the same manner as traditional cationic drug delivery systems or cell-penetrating peptides. Based on our findings, HBpep coacervates at physiological pH cannot cross the cell membrane by a simple passive mechanism and are thus likely to adopt a non-canonical cell entry pathway.

Indoleamine 2,3-dioxygenase 1 (IDO1) is the key enzyme that catalyzes the conversion of l-tryptophan (Trp) to N-formyl-kynurenine in the tryptophan-kynurenine (Trp-Kyn) pathway. Overexpression of IDO1 contributes to the depletion of Trp and the accumulation of Kyn, which can result in tumor immune escape. Inhibition of IDO1 can restore the host immune response to eradicate cancer cells. 5-(Pyridin-3-yl)-1H-indole-4,7-dione was developed as the scaffold for a type of IDO1 inhibitors from simplification of the structure of exiguamine A in our previous work. In the present study, we designed and synthesized a series of compounds with different side-chain substituents, and linkers of varying lengths, at the 3-position, to design compounds with an aryl motif that can occupy pocket B of the IDO1 protein. Most compounds exhibited potent IDO1 inhibitory activity with IC50 values at the micromolar level, and 3-(2-((4-fluorobenzyl)amino)ethyl)-5-(3-pyridyl)-1H-indole-4,7-dione (1d) displayed the most potent inhibition with a half-maximal inhibitory concentration (IC50) value of 0.125 μM in an enzymatic assay and a half-maximal effective concentration (EC50) value of 0.821 μM in a cellular assay. Compound 1d showed higher selectivity for IDO1 over indoleamine 2,3-dioxygenase 2 and tryptophan 2,3-dioxygenase at the effective concentration. Molecular docking studies and molecular dynamic simulations suggested that the phenyl ring of these inhibitors can enter into pocket B of IDO1 and interact with the residue Phe 226 through hydrophobic interactions.

2-phenoxyanline was thoroughly examined using several spectroscopy nuclear magnetic resonance, IR, UV–vis, and quantum mechanical methodologies. D F T with basis set is employed to determine structural optimization and distinct vibration modes. The optimum binding parameters agree closely with the observed binding characteristics. VEDA performed the duties relating to potential energy distribution (PED) satisfactorily. The G I A O approach was adopted to compute chemical shifts in the 13C, 1H- N.M.R and the outcomes were compared to experimental spectra. TD-DFT with PCM model was used to calculate UV–vis with different solvents that, when assessed against observational spectra. The FMO energy provides sufficient proof of such. Molecular docking and dynamic simulations gave a better idea of the interaction of ligands with receptors.

We provide promising computational (in silico) data on phytochemicals (compounds 1–10) from Arabian Peninsula medicinal plants as strong binders, targeting 3-chymotrypsin-like protease (3CLPro) and papain-like proteases (PLPro) of SARS-CoV-2. Compounds 1–10 followed the Lipinski rules of five (RO5) and ADMET analysis, exhibiting drug-like characters. Non-covalent (reversible) docking of compounds 1–10 demonstrated their binding with the catalytic dyad (CYS145 and HIS41) of 3CLPro and catalytic triad (CYS111, HIS272, and ASP286) of PLPro. Moreover, the implementation of the covalent (irreversible) docking protocol revealed that only compounds 7, 8, and 9 possess covalent warheads, which allowed the formation of the covalent bond with the catalytic dyad (CYS145) in 3CLPro and the catalytic triad (CYS111) in PLPro. Root-mean-square deviation (RMSD), root-mean-square fluctuation (RMSF), and radius of gyration (Rg) analysis from molecular dynamic (MD) simulations revealed that complexation between ligands (compounds 7, 8, and 9) and 3CLPro and PLPro was stable, and there was less deviation of ligands. Overall, the in silico data on the inherent properties of the above phytochemicals unravel the fact that they can act as reversible inhibitors for 3CLPro and PLPro. Moreover, compounds 7, 8, and 9 also showed their novel properties to inhibit dual targets by irreversible inhibition, indicating their effectiveness for possibly developing future drugs against SARS-CoV-2. Nonetheless, to confirm the theoretical findings here, the effectiveness of the above compounds as inhibitors of 3CLPro and PLPro warrants future investigations using suitable in vitro and in vivo tests.

Herein, network pharmacology was used to identify the active components in Ilex kudingcha and common hypertension-related targets. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were conducted, and molecular docking was performed to verify molecular dynamic simulations. Six active components in Ilex kudingcha were identified; furthermore, 123 target genes common to hypertension were identified. Topological analysis revealed the strongly associated proteins, with RELA, AKT1, JUN, TP53, TNF, and MAPK1 being the predicted targets of the studied traditional Chinese medicine. In addition, GO enrichment analysis revealed significant enrichment of biological processes such as oxidative stress, epithelial cell proliferation, cellular response to chemical stress, response to xenobiotic stimulus, and wound healing. Furthermore, KEGG enrichment analysis revealed that the genes were particularly enriched in lipid and atherosclerosis, fluid shear stress and atherosclerosis, and other pathways. Molecular docking revealed that the key components in Ilex kudingcha exhibited good binding potential to the target genes RELA, AKT1, JUN, TP53, TNF, and IL-6. Our study results suggest that Ilex kudingcha plays a role in hypertension treatment by exerting hypolipidemic, anti-inflammatory, and antioxidant effects and inhibiting the transcription of atherosclerosis-related genes.

Guaiane-type sesquiterpenoids are most prevalent in the genus Cinnamomum. Hence this study investigates the structures, anti-nociceptive and IL-6 targeted anti-inflammatory potential of three novels C-14 guaiane-type sesquiterpenoids and two new monoterpenoids, isolated from Cinnamomum migao. The structures were precisely confirmed and characterized through the modern chromatographic and spectroscopic techniques of HRESIMS, 1D NMR, 2D NMR, experimental circular dichroism (ECD), and calculated (ECD). Novel sesquiterpenoids 1 and 2 exhibited significant anti-inflammatory activities against the NO production and pro-inflammatory cytokines. Their IC50 values were determined as 9.52 and 13.42 μΜ against IL-6 mRNA, respectively. Similarly, subcutaneous injection of n-BuT and EA extracts showed a dose-dependent suppression of formalin-induced tonic biting/licking responses during the tonic antinociceptive phase. Furthermore, absorption, distribution, metabolism, excretion, and toxicity (ADMET) analysis of guaiane-type sesquiterpenoids 1 and 2 displayed that both compounds have a high level of GIT absorption, with a high zone of safety for cardiac and hepatotoxicity and no inhibition of cytochromes. In addition, molecular docking and simulation studies strengthen the anti-inflammatory potential of sesquiterpene 2 which showed a good binding affinity with IL-6 protein. Overall the inclusive results showed that the extracts and newly isolated guaiane-type sesquiterpenoids from C. migao will provide new evidence for the traditional use of this species to treat inflammation and nociception.

Objective: To investigate the renoprotective effects of luteolin on diabetes in rats. Methods: One week after administration of streptozotocin 55 mg/kg intraperitoneally, rats were given 25, 50, and 75 mg/kg/day of luteolin orally for another eight weeks. At the end of the experiment, body weight, blood glucose level, biochemical parameters for renal function (serum creatinine, blood urea nitrogen, uric acid, serum albumin, and total protein), kidney histology, matrix metalloproteinase (MMP)-2, MMP-9, and histone deacetylase 2 (HDAC-2) expression, and malondialdehyde, myeloperoxidase, and hydroxyproline content in renal tissue were evaluated. High glucose-induced damage using NRK-52E cell line was studied to evaluate cell viability and metalloenzyme expression. Additionally, in silico studies including docking and molecular dynamics simulations were conducted. Results: MMP-2, MMP-9, and HDAC-2 expressions were significantly increased in high glucose-induced NRK-52E cells and the renal tissue of diabetic rats. However, these changes were reversed by luteolin at the administered doses. Additionally, luteolin significantly reduced oxidative stress, inflammation, and fibrosis, as well as improved biochemical parameters in diabetic rats. Furthermore, luteolin at the examined doses markedly alleviated diabetes-induced histopathological changes in renal tissues. Conclusions: Luteolin effectively attenuates streptozotocin-induced diabetic nephropathy in rats by inhibiting MMP-2, MMP-9, and HDAC-2 expression, and reducing oxidative stress and inflammation.

Monkeypox is a zoonotic viral infection caused by monkeypox virus which belongs to the Poxviridae family of genus Orthopoxvirus. Usually, the virus transmission happens when the individual comes in contact with the infected person through body fluids, animal lesions, respiratory droplets or through virus contaminated materials. Clinical presentation of the monkeypox has shown significant resemblance to that of smallpox and chickenpox, belonging to the same orthopoxvirus genus but were eradicated during 1980s globally. Monkeypox may lead to a range of medical complications including clinical symptoms like fever, rashes, headaches, back pain, myodynia and swollen lymph nodes. As far as the treatment modalities are concerned, the antiviral therapeutic agents developed for the smallpox treatment, were also permitted to be used for the monkeypox treatment. However, there is no proven treatment for human monkeypox. In the current study, we have focused on designing of a best probable ligand against the target MPXVgp158 (Monkeypox virus protein). Since Tecovirimat is an FDA approved compound known as an antipoxviral drug, the study aimed to develop a Monkeypox virus protein MPXVgp158 inhibitor which is bioavailable and biocompatible as well through drug designing using computational tools. Molecular docking (MD) analysis displayed Tecovirimat with lesser binding energy, higher non-bonded interaction capability, and more stability against MPXVgp158, with efficient binding mode of interactions. Hence, Tecovirimat was adjudged to be the potential candidate against MPXVgp158 inhibition.

Abstract Background Alpha-1 antitrypsin deficiency (AATD) is characterized by reduced serum levels of the AAT protein and predisposes to liver and lung disease. The characterization at structural level of novel pathogenic SERPINA1 mutants coding for circulating AAT could provide novel insights into the mechanisms of AAT misfolding. The present study aimed to provide a practical framework for the identification and analysis of new AAT mutations, combining structural simulations and clinical data. Methods We analysed a total of five mutations (four not previously described) in a total of six subjects presenting moderate to severe AATD: Gly95Alafs*18, Val210Glu, Asn247Ser, Pi*S + Asp341His and Pi*S + Leu383Phe + Lys394Ile. Clinical data, genotyping and phenotyping assays, structural mapping, and conformational characterization through molecular dynamic (MD) simulations were developed and combined. Results Newly discovered AAT missense variants were localized both on the interaction surface and the hydrophobic core of the protein. Distribution of mutations across the structure revealed Val210Glu at the solvent exposed s4C strand and close to the “Gate” region. Asn247Ser was located on the accessible surface, which is important for glycan attachment. On the other hand, Asp341His, Leu383Phe were mapped close to the “breach” and “shutter” regions. MD analysis revealed the reshaping of local interactions around the investigated substitutions that have varying effects on AAT conformational flexibility, hydrophobic packing, and electronic surface properties. The most severe structural changes were observed in the double- and triple-mutant (Pi*S + Asp341His and Pi*S + Leu383Phe + Lys394Ile) molecular models. The two carriers presented impaired lung function. Conclusions The results characterize five variants, four of them previously unknown, of the SERPINA1 gene, which define new alleles contributing to the deficiency of AAT. Rare variants might be more frequent than expected, and therefore, in discordant cases, standardized screening of the S and Z alleles needs complementation with gene sequencing and structural approaches. The utility of computational modelling for providing supporting evidence of the pathogenicity of rare single nucleotide variations is discussed.

Background: SARS-CoV-2 is a deadly viral disease and uncounted deaths occurs since its first appearance in the year 2019. The antiviral drugs, benzylisoquinoline alkaloids, and coumarin molecules were searched using different online engines for drug repurposing with SARS-CoV-2 and to investigate the effects on main viral protease (Mpro) upon their bindings. Methods: A database composed of antiviral drugs, benzylisoquinoline alkaloids, and Coumarin molecules was screened through a molecular docking strategy to uncover the interactions of collected molecules with SARS-CoV-2 Mpro. Further, molecular dynamics simulations (MDS) were implemented for 100 ns to calculate the stability of the best complexed molecular scaffold with Mpro. The conformations of the simulated complexes were investigated by using principal component analysis (PCA) and Gibbs energy landscape (FEL) and DSSP together. Next, free binding energy (ΔGbind) was calculated using the mmpbsa method. Results: Molecular docking simulations demonstrate 17 molecules exhibited better binding affinity out of 99 molecules present in the database with the viral protease Mpro, followed ADMET properties and were documented. The Coumarin-EM04 molecular scaffold exhibited interactions with catalytical dyad HIS41, CYS145, and neighboring amino acids SER165 and GLN189 in the catalytical site. The crucial factor RMSD was calculated to determine the orientations of Coumarin-EM04. The Coumarin-EM04 complexed with Mpro was found stable in the binding site during MDS. Furthermore, the free energy binding ΔGbind of Coumarin-EM04 was found to be −187.471 ± 2.230 kJ/mol, and for Remdesivir ΔGbind was −171.926 ± 2.237 kJ/mol with SARS-CoV-2 Mpro. Conclusion: In this study, we identify potent molecules that exhibit interactions with catalytical dyad HIS41 and CYS145 amino acids and unravel Coumarin-EM04 exhibited ΔGbind higher than Remdesivir against Mpro and thus may serve better antiviral agent against SARS-CoV-2.

The investigation of black phosphorus-based (BP)-based mass sensors provides theoretical support for the development of mass-detection devices. This study examines the non-linear dynamic behavior of a rectangular single-layered BP (SLBP) with an attached mass through the utilization of molecular dynamics (MD) simulations and a nonlinear orthotropic plate model (NOPM) with a concentrated mass. The results indicate that significant deformation of an SLBP with an attached mass necessitates consideration of geometric nonlinearity, although the attached mass does not affect the deformation. Additionally, this paper discusses the impact of the attachment mass and the amplitude of harmonic force on the non-linear forced vibration of the SLBP. It is observed that as the attachment mass increases, the nonlinear vibration resonance frequency decreases, while the peak amplitude increases. Furthermore, the thermal nonlinear vibration of SLBP with an attached mass has been investigated, revealing that an increase in the attached mass leads to a decrease in the nonlinear vibration frequency but an increase in the amplitude of the nonlinear vibration of SLBP with an attached mass. Overall, comparison with MD simulation results, this investigation suggests that NOPM with a concentrated mass effectively describes the nonlinear vibration behavior of an SLBP with an attached mass, providing theoretical support for designing such devices to detect attached masses.

RuBisCO is the main photosynthetic enzyme of carbon assimilatory pathways in nature. Despite being the most abundant protein on earth, RuBisCO is still relatively underutilised in the food chain. Although there are sequence and structure details in the database, there are few instances of studies on evolutionary relationships. A bioinformatics and in silico study was conducted to check sequence and structural differences of RuBisCO among different photosynthetic organisms. RuBisCO from Oryza sativa showed an abundance of charged amino acids, salt-bridges, and intra-protein interactions and was more hydrophilic in nature compared to Nostoc sp., Chlamydomonas reinhardtii, and Nicotiana tabacum. From molecular dynamics simulations, lower root mean square deviation and root mean square fluctuation indicate that RuBisCO from Oryza sativa was more stable, followed by Nicotiana tabacum, and a lower radius of gyrations indicates their tight packing. From this study, it was clear that some specific evolutions in charged amino acids of RuBisCO of monocot, i.e., Oryza sativa, make it more stable and stronger than other plant groups. The study concludes that a more stable nature of RuBisCO is gained from monocot Oryza sativa.

Fatty acid transport protein 1 (FATP1) is an integral transmembrane protein that is involved in facilitating the translocation of long-chain fatty acids (LCFA) across the plasma membrane, thereby orchestrating the importation of LCFA into the cell. FATP1 also functions as an acyl-CoA ligase, catalyzing the ATP-dependent formation of fatty acyl-CoA using LCFA and VLCFA (very-long-chain fatty acids) as substrates. It is expressed in various types of tissues and is involved in the regulation of crucial signalling pathways, thus playing a vital role in numerous physiological and pathological conditions. Structural insight about FATP1 is, thus, extremely important for understanding the mechanism of action of this protein and developing efficient treatments against its anomalous expression and dysregulation, which are often associated with pathological conditions such as breast cancer. As of now, there has been no prior prediction or evaluation of the 3D configuration of the human FATP1 protein, hindering a comprehensive understanding of the distinct functional roles of its individual domains. In our pursuit to unravel the structure of the most commonly expressed isoforms of FATP1, we employed the cutting-edge ALPHAFOLD 2 model for an initial prediction of the entire protein’s structure. This prediction was complemented by molecular dynamics simulations, focusing on the most promising model. We predicted the structure of FATP1 in silico and thoroughly refined and validated it using coarse and molecular dynamics in the absence of the complete crystal structure. Their relative dynamics revealed the different properties of the characteristic FATP1.

Over time, the antigenic evolution of emerging variants of SARS-CoV-2 has demanded the development of potential protective vaccines. Administration of additional doses of current vaccines based on the WT spike protein may boost immunity, but their effectiveness has dwindled for patients with more recent variants. Here, we studied the neutralization activity of post-WT strain-based vaccination and a structural simulation in-silico based on the interactions of the RBD-hACE2 as the key to initiating infection among the VOCs of SARS-CoV-2. Our data display shows that WT sera showed a markedly greater reduction in Delta and Omicron, suggesting that the Wuhan-based vaccines may be more susceptible to breakthrough and new VOCs. According to the MD simulation, mutations of Omicron result in a significant change in the variant charge distribution throughout the binding interface that consequently alters the critical interface electrostatic potential in comparison to other variants. This observation provides new insights into immunization policy and next-generation vaccine development.

Many endothermic liquid–liquid transitions, occurring at a temperature Tn+ above the melting temperature Tm, are related to previous exothermic transitions, occurring at a temperature Tx after glass formation below Tg, with or without attached crystallization and predicted by the nonclassical homogenous nucleation equation. A new thermodynamic phase composed of broken bonds (configurons), driven by percolation thresholds, varying from ~0.145 to Δε, is formed at Tx, with a constant enthalpy up to Tn+. The liquid fraction Δε is a liquid glass up to Tn+. The solid phase contains glass and crystals. Molecular dynamics simulations are used to induce, in NiTi2, a reversible first-order transition by varying the temperature between 300 and 1000 K under a pressure of 1000 GPa. Cooling to 300 K, without applied pressure, shows the liquid glass presence with Δε = 0.22335 as memory effect and Tn+ = 2120 K for Tm = 1257 K.

Chitin is one of the most common polysaccharides and is abundant in the cell walls of fungi and the shells of insects and aquatic organisms as a skeleton. The mechanism of how chitin responds to pH is essential to the precise control of brewing and the design of smart chitin materials. However, this molecular mechanism remains a mystery. Results from single-molecule studies, including single-molecule force spectroscopy (SMFS), AFM imaging, and molecular dynamic (MD) simulations, have shown that the mechanical and conformational behaviors of chitin molecules show surprising pH responsiveness. This can be compared with how, in natural aqueous solutions, chitin tends to form a more relaxed spreading conformation and show considerable elasticity under low stretching forces in acidic conditions. However, its molecular chain collapses into a rigid globule in alkaline solutions. The results show that the chain state of chitin can be regulated by the proportions of inter- and intramolecular H-bonds, which are determined via the number of water bridges on the chain under different pH values. This basic study may be helpful for understanding the cellular activities of fungi under pH stress and the design of chitin-based drug carriers.

The escalating prevalence of drug-resistant strains of Mycobacterium tuberculosis has posed a significant challenge to global efforts in combating tuberculosis. To address this issue, innovative therapeutic strategies are required that target essential biochemical pathways while minimizing the potential for resistance development. The concept of dual targeting has gained prominence in drug discovery against resistance bacteria. Dual targeting recognizes the complexity of cellular processes and disrupts more than one vital pathway, simultaneously. By inhibiting more than one essential process required for bacterial growth and survival, the chances of developing resistance are substantially reduced. A previously reported study investigated the dual-targeting potential of a series of novel compounds against the folate pathway in Mycobacterium tuberculosis. Expanding on this study, we investigated the predictive pharmacokinetic profiling and the structural mechanism of inhibition of UCP1172, UCP1175, and UCP1063 on key enzymes, dihydrofolate reductase (DHFR) and 5-amino-6-ribitylamino-2,4(1H,3H)-pyrimidinedione 5′-phosphate reductase (RV2671), involved in the folate pathway. Our findings indicate that the compounds demonstrate lipophilic physiochemical properties that promote gastrointestinal absorption, and may also inhibit the drug-metabolizing enzyme, cytochrome P450 3A4, thus enhancing their biological half-life. Furthermore, key catalytic residues (Serine, Threonine, and Aspartate), conserved in both enzymes, were found to participate in vital molecular interactions with UCP1172, which demonstrated the most favorable free binding energies to both DHFR and RV2671 (−41.63 kcal/mol, −48.04 kcal/mol, respectively). The presence of characteristic loop shifts, which are similar in both enzymes, also indicates a common inhibitory mechanism by UCP1172. This elucidation advances the understanding of UCP1172’s dual inhibition mechanism against Mycobacterium tuberculosis.

Oleanolic acid (OA), as a ubiquitous compound in the plant kingdom, is studied for both its neuroprotective and neurotoxic properties. The mechanism of acetylcholinesterase (AChE) inhibitory potential of OA is investigated using molecular dynamic simulations (MD) and docking as well as biomimetic tests. Moreover, the in vitro SH-SY5Y human neuroblastoma cells and the in vivo zebrafish model were used. The inhibitory potential towards the AChE enzyme is examined using the TLC-bioautography assay (the IC50 value is 9.22 μM). The CH-π interactions between the central fragment of the ligand molecule and the aromatic cluster created by the His440, Phe288, Phe290, Phe330, Phe331, Tyr121, Tyr334, Trp84, and Trp279 side chains are observed. The results of the in vitro tests using the SH-SY5Y cells indicate that the viability rate is reduced to 71.5%, 61%, and 43% at the concentrations of 100 µg/mL, 300 µg/mL, and 1000 µg/mL, respectively, after 48 h of incubation, whereas cytotoxicity against the tested cell line with the IC50 value is 714.32 ± 32.40 µg/mL. The in vivo tests on the zebrafish prove that there is no difference between the control and experimental groups regarding the mortality rate and morphology (p > 0.05).

SARS-CoV-2 exploits the homotrimer transmembrane Spike glycoproteins (S protein) during host cell invasion. The Omicron XBB subvariant, delta, and prototype SARS-CoV-2 receptor-binding domain show similar binding strength to hACE2 (human Angiotensin-Converting Enzyme 2). Here we utilized multiligand virtual screening to identify small molecule inhibitors for their efficacy against SARS-CoV-2 virus using QPLD, pseudovirus ACE2 Inhibition -Time Resolved Forster/Fluorescence energy transfer (TR-FRET) Assay Screening, and Molecular Dynamics simulations (MDS). Three hundred and fifty thousand compounds were screened against the macrodomain of the nonstructural protein 3 of SARS-CoV-2. Using TR-FRET Assay, we filtered out two of 10 compounds that had no reported activity in in vitro screen against Spike S1: ACE2 binding assay. The percentage inhibition at 30 µM was found to be 79% for “Compound F1877-0839” and 69% for “Compound F0470-0003”. This first of its kind study identified “FILLY” pocket in macrodomains. Our 200 ns MDS revealed stable binding poses of both leads. They can be used for further development of preclinical candidates.

Muhammad Faheem,1 Yusuf S Althobaiti,2,3 Abdul Waheed Khan,4 Aman Ullah,1 Syed Hussain Ali,1 Umair Ilyas1 1Riphah Institute of Pharmaceutical Sciences, Riphah International University, Islamabad, Pakistan; 2Department of Pharmacology and Toxicology, College of Pharmacy, Taif University, Taif, 21944, Saudi Arabia; 3Addiction and Neuroscience Research Unit, Taif University, Taif, 21944, Saudi Arabia; 4Department of Pharmacy, The University of Lahore, Islamabad, PakistanCorrespondence: Muhammad Faheem; Yusuf S Althobaiti Email muhammad.faheem@riphah.edu.pk; althobaiti@tu.edu.saObjective: The study investigated the effect 5-[(naphthalen-2-yloxy) methyl]-1,3,4-oxadiaszole2-thiol (B3) in animal model of acute epileptic shock.Methods: The pharmacokinetics profile of B3 was checked through SwissADME software. The binding affinities of B3, diazepam, and flumazenil (FLZ) were obtained through Auto Dock and PyRx. Post docking analysis and interpretation of hydrogen bonds were performed through Discovery Studio Visualizer 2016. Molecular dynamics simulations of three complexes were carried out through Desmond software package. B3 was then proceeded in PTZ-induced acute seizures models. Flumazenil was used in animal studies for elucidation of possible mechanism of B3. After behavioral studies, the animals were sacrificed, and the brain samples were isolated and stored in 4% formalin for molecular investigations including H and E staining, IHC staining and Elisa etc.Results: The results demonstrate that B3 at 20 and 40 mg/kg prolonged the onset time of generalized seizures. B3 considerably increased the expression of protective glutathione S-transferase and glutathione reductase and reduced lipid peroxidation and inducible nitric oxide synthase (P < 0.001) in the cortex. B3 significantly suppressed (P < 0.01) the over expression of the inflammatory mediator tumor necrosis factor–α, whose up-regulation is reported in acute epileptic shocks.Conclusion: Hence, it is concluded from the aforementioned results that B3 provides neuroprotective effects PTZ-induced acute epileptic model. FLZ pretreatment resulted in inhibition of the anticonvulsant effect of B3. B3 possesses anticonvulsant effect which may be mediated through GABAA mediated antiepileptic pathway.Keywords: docking, molecular dynamic simulations, gamma amino butyric acid A

This work investigated the adsorption mechanism of heavy metals: silver (Ag), mercury (Hg), cadmium (Cd), palladium (Pd), and zinc (Zn) on the carbon graphite (111) surface based on density functional theory (DFT), Metropolis Monte Carlo (MMC), and the molecular dynamics (MD) simulation methods. The obtained results from the adsorption of these species showed that the process is spontaneous and exothermic in nature. The maximum adsorption capacities were obtained in neutral to low acid medium, and the interaction between Hg-carbon graphite was more favored than other systems. These findings showed that carbon graphite was more efficient in the removal of the studied metals. Moreover, this study better explains the adsorption mechanism of heavy metals onto carbon graphite and gives a theoretical basis for the wider application of graphite adsorbent in the removal of heavy metals. This paper provides theoretical support for heavy metals removal in various mediums and also provides some new ideas for the secondary utilization of quantum chemical descriptors and molecular dynamics simulation methods.

The activation of caspases is central to apoptotic process in living systems. Defects in apoptosis have been implicated with carcinogenesis. Need to develop smart agents capable of inducing apoptosis in tumor cells is obvious. With this motive, diversity oriented synthesis of 1-benzylpyrrolidin-3-ol analogues was envisaged. The multi component Ugi reaction synthesized library of electronically diverse analogues was explored for cytotoxic propensity towards a panel of human cancer cell lines at 10 µM. The lead compounds exhibit a selective cytotoxicity towards HL-60 cells as compared to cell lines derived from solid tumors. Besides, their milder cytotoxic effect on non-cancerous cell lines reaffirm their selective action towards cancer cells only. The lead molecules were tested for their ability to target caspase-3, as a vital protease triggering apoptosis. The lead compounds were observed to induce apoptosis in HL-60 cells around 10 µM concentration. The lead compounds exhibited various non-covalent supra type interactions with caspase-3 key residues around the active site. The binding ability of lead compounds with caspase-3 was studied via molecular docking and molecular dynamic (MD) simulations. MD simulations indicated the stability of compound-caspase-3 complex throughout the 50 ns simulation run. The stability and bio-availability of the lead compounds under physiological conditions was assessed by their interaction with Bovine Serum Albumin (BSA) as model protein. BSA interactions of lead compounds were studied by various bio-physical methods and further substantiated with in silico MD simulations.

Abstract Background Canine circovirus is a deadly pathogen of dogs and causes vasculitis and hemorrhagic enteritis. It causes lethal gastroenteritis in pigs, fox, and dogs. Canine circovirus genome contains two main (and opposite) transcription units which encode two open reading frames (ORFs), a replicase-associated protein (Rep) and the capsid (Cap) protein. The replicase protein and capsid protein consist of 303 amino acids and 270 amino acids respectively. Several immuno-informatics methods such as epitope screening, molecular docking, and molecular-dynamics simulations were used to craft peptide-based vaccine construct against canine circovirus. Results The vaccine construct was designed by joining the selected epitopes with adjuvants by suitable linker. The cloning and expression of the vaccine construct was also performed using in silico methods. Screening of epitopes was conducted by NetMHC server that uses ANN (Artificial neural networking) algorithm. These methods are fast and cost-effective for screening epitopes that can interact with dog leukocyte antigens (DLA) and initiate an immune response. Overall, 5 epitopes, YQHLPPFRF, YIRAKWINW, ALYRRLTLI, HLQGFVNLK, and GTMNFVARR, were selected and used to design a vaccine construct. The molecular docking and molecular dynamics simulation studies show that these epitopes can bind with DLA molecules with stability. The codon adaptation and in silico cloning studies show that the vaccine can be expressed by Escherichia coli K12 strain. Conclusion The results suggest that the vaccine construct can be useful in preventing the dogs from canine circovirus infections. However, the results need further validation by performing other in vitro and in vivo experiments.

Coronavirus disease 2019 (COVID-19) is a pandemic whose adverse effects have been felt all over the world. As of August 2022, reports indicated that over 500 million people in the world had been infected and the number of rising deaths from the disease were slightly above 6.4 million. New variants of the causative agent, SARS-CoV-2 are emanating now and then and some are more efficacious and harder to manage. SARS-CoV-2 main protease (Mpro) has essential functions in viral gene expression and replication through proteolytic cleavage of polyproteins. Search for SARS-CoV-2 Mpro inhibitors is a vital step in the treatment and management of COVID-19. In this study, we investigated whether alkaloids with antiviral and myriad other bioactivities from the genus Lycoris can act as SARS-CoV-2 Mpro inhibitors. We conducted a computer-aided drug design study through screening optimal ligands for SARS-CoV-2 Mpro from a list of over 150 Lycoris alkaloids created from online databases such as ChEMBL, PubChem, ChemSpider, and published journal papers. The In silico study involved molecular docking of Lycoris alkaloids to SARS-CoV-2 Mpro active site, absorption, distribution, metabolism, elimination and toxicity (ADMET) screening and finally molecular dynamic (MD) simulations of the most promising ligand-SARS-CoV-2 Mpro complexes. The study identified 3,11-dimethoxy-lycoramine, narwedine, O-demethyllycoramine and epilycoramine as drug-like and lead-like Lycoris alkaloids with favorable ADMET properties and are very likely to have an inhibition activity on SARS-CoV-2 Mpro and may become potential drug candidates. DOI: http://dx.doi.org/10.5281/zenodo.7041808

Janus tyrosine kinase 3 (JAK3) is primarily expressed in immune cells and is needed for signaling by the common gamma chain (γc) family of cytokines. Abnormal JAK3 signal transduction can manifest as hematological disorders, e.g., leukemia, severe combined immunodeficiency (SCID) and autoimmune disease states. While regulatory JAK3 phosphosites have been well studied, here a functional proteomics approach coupling a JAK3 autokinase assay to mass spectrometry revealed ten previously unreported autophosphorylation sites (Y105, Y190, Y238, Y399, Y633, Y637, Y738, Y762, Y824, and Y841). Of interest, Y841 was determined to be evolutionarily conserved across multiple species and JAK family members, suggesting a broader role for this residue. Phospho-substitution mutants confirmed that Y841 is also required for STAT5 tyrosine phosphorylation. The homologous JAK1 residue Y894 elicited a similar response to mutagenesis, indicating the shared importance for this site in JAK family members. Phospho-specific Y841-JAK3 antibodies recognized activated kinase from various T-cell lines and transforming JAK3 mutants. Computational biophysics analysis linked Y841 phosphorylation to enhanced JAK3 JH1 domain stability across pH environments, as well as to facilitated complementary electrostatic JH1 dimer formation. Interestingly, Y841 is not limited to tyrosine kinases, suggesting it represents a conserved ubiquitous enzymatic function that may hold therapeutic potential across multiple kinase families.

The monkeypox virus (MPXV) is an enveloped, double-stranded DNA virus belonging to the genus Orthopox viruses. In recent years, the virus has spread to countries where it was previously unknown, turning it into a worldwide emergency for public health. This study employs a structural-based drug design approach to identify potential inhibitors for the core cysteine proteinase of MPXV. During the simulations, the study identified two potential inhibitors, compound CHEMBL32926 and compound CHEMBL4861364, demonstrating strong binding affinities and drug-like properties. Their docking scores with the target protein were −10.7 and −10.9 kcal/mol, respectively. This study used ensemble-based protein–ligand docking to account for the binding site conformation variability. By examining how the identified inhibitors interact with the protein, this research sheds light on the workings of the inhibitors’ mechanisms of action. Molecular dynamic simulations of protein–ligand complexes showed fluctuations from the initial docked pose, but they confirmed their binding throughout the simulation. The MMGBSA binding free energy calculations for CHEMBL32926 showed a binding free energy range of (−9.25 to −9.65) kcal/mol, while CHEMBL4861364 exhibited a range of (−41.66 to −31.47) kcal/mol. Later, analogues were searched for these compounds with 70% similarity criteria, and their IC50 was predicted using pre-trained machine learning models. This resulted in identifying two similar compounds for each hit with comparable binding affinity for cysteine proteinase. This study’s structure-based drug design approach provides a promising strategy for identifying new drugs for treating MPXV infections.

Amphiphilic peptides, such as Aß amyloids, can adsorb at an interface between two immiscible electrolyte solutions (ITIES). Based on previous work (vide infra), a hydrophilic/hydrophobic interface is used as a simple biomimetic system for studying drug interactions. The ITIES provides a 2D interface to study ion-transfer processes associated with aggregation, as a function of Galvani potential difference. Here, the aggregation/complexation behaviour of Aβ(1-42) is studied in the presence of Cu (II) ions, together with the effect of a multifunctional peptidomimetic inhibitor (P6). Cyclic and differential pulse voltammetry proved to be particularly sensitive to the detection of the complexation and aggregation of Aβ(1-42), enabling estimations of changes in lipophilicity upon binding to Cu (II) and P6. At a 1:1 ratio of Cu (II):Aβ(1-42), fresh samples showed a single DPV (Differential Pulse Voltammetry) peak half wave transfer potential (E1/2) at 0.40 V. Upon increasing the ratio of Cu (II) two-fold, fluctuations were observed in the DPVs, indicating aggregation. The approximate stoichiometry and binding properties of Aβ(1-42) during complexation with Cu (II) were determined by performing a differential pulse voltammetry (DPV) standard addition method, which showed two binding regimes. A pKa of 8.1 was estimated, with a Cu:Aβ1-42 ratio~1:1.7. Studies using molecular dynamics simulations of peptides at the ITIES show that Aβ(1-42) strands interact through the formation of β-sheet stabilised structures. In the absence of copper, binding/unbinding is dynamic, and interactions are relatively weak, leading to the observation of parallel and anti-parallel arrangements of β-sheet stabilised aggregates. In the presence of copper ions, strong binding occurs between a copper ion and histidine residues on two peptides. This provides a convenient geometry for inducing favourable interactions between folded β-sheet structures. Circular Dichroism spectroscopy (CD spectroscopy) was used to support the aggregation behaviour of the Aβ(1-42) peptides following the addition of Cu (II) and P6 to the aqueous phase.

The newly FDA-approved drug, Axitinib, is an effective therapy against RTKs, but it possesses severe adverse effects like hypertension, stomatitis, and dose-dependent toxicity. In order to ameliorate Axitinib’s downsides, the current study is expedited to search for energetically stable and optimized pharmacophore features of 14 curcumin (1,7-bis(4-hydroxy-3-methoxyphenyl)hepta-1,6-diene-3,5-dione) derivatives. The rationale behind the selection of curcumin derivatives is their reported anti-angiogenic and anti-cancer properties. Furthermore, they possessed a low molecular weight and a low toxicity profile. In the current investigation, the pharmacophore model-based drug design, facilitates the filtering of curcumin derivatives as VEGFR2 interfacial inhibitors. Initially, the Axitinib scaffold was used to build a pharmacophore query model against which curcumin derivatives were screened. Then, top hits from pharmacophore virtual screening were subjected to in-depth computational studies such as molecular docking, density functional theory (DFT) studies, molecular dynamics (MD) simulations, and ADMET property prediction. The findings of the current investigation revealed the substantial chemical reactivity of the compounds. Specifically, compounds S8, S11, and S14 produced potential molecular interactions against all four selected protein kinases. Docking scores of −41.48 and −29.88 kJ/mol for compounds S8 against VEGFR1 and VEGFR3, respectively, were excellent. Whereas compounds S11 and S14 demonstrated the highest inhibitory potential against ERBB and VEGFR2, with docking scores of −37.92 and −38.5 kJ/mol against ERBB and −41.2 and −46.5 kJ/mol against VEGFR-2, respectively. The results of the molecular docking studies were further correlated with the molecular dynamics simulation studies. Moreover, HYDE energy was calculated through SeeSAR analysis, and the safety profile of the compounds was predicted through ADME studies.

A quantitative structure-property relationship (QSPR) in-silico study was performed to develop a mathematical model that correlates 2D and 3D descriptors of 37 antioxidant lubricant additives (compounds) with their properties. A molecular dynamics simulation study was also carried out to access these additives' binding strength on diamond-like carbon (DLC) and steel crystal surfaces. Five novel antioxidant lubricant additives were designed from the information derived from the QSPR mathematical model’s high coefficient molecular descriptors. All the novel lubricant additive’s antioxidant properties were found to be better than our previous study, with the lubricant additive (Z)-3-(4-(5-amino-1-phenyl-1H-pyrazol-3-yl)-3,5-dimethylphenyl)-2-phenyl-5-(thiophen-2-ylmethylene)-3,5-dihydro-4H-imidazol-4-one found to possessed excellent antioxidant properties of 0.850281 total acid values (T.A.V 0.1g/L) than its co-additives. Moreover, all the designed additives dynamically bind to steel crystal surfaces excellently from our dynamic simulation study than the DLC crystal surface. The molecular dynamics simulation results were found to be better than the one reported by our previous study. This investigation will help synthesize novel and excellent antioxidant lubricant additives that will hinder the base oil from undergoing a complete oxidation cycle and meet environmental requirements as these novel additives do not contain Zinc and Phosphorus, which often rendered exhaust pipes catalytic converter inactive, thereby increasing environmental pollution.

Abstract Background Lipases (EC 3.1.1.3) catalyze the hydrolysis of oil into free fatty acids and glycerol forming the 3rd largest group of commercialized enzymes. Plant lipases grab attention recently because of their specificity, less production and purified cost, and easy availability. In silico approach is the first step to identify different genes coding for lipase in a most common indigenous plant, wheat, to explore the possibility of this plant as an alternative source for commercial lipase production. As the hierarchy organization of genes reflects an ancient process of gene duplication and divergence, many of the theoretical and analytical tools of the phylogenetic systematics can be utilized for comparative genomic studies. Also, in addition to experimental identification and characterization of genes, for computational genomic analysis, Arabidopsis has become a popular strategy to identify crop genes which are economically important, as Arabidopsis genes had been well identified and characterized for lipase. A number of articles had been reported in which genes of wheat have shown strong homology with Arabidopsis. The complete genome sequences of rice and Arabidopsis constitute a valuable resource for comparative genome analysis as they are representatives of the two major evolutionary lineages within the angiosperms. Here, in this in silico approach, Arabidopsis and Oryza sativa serve as models for dicotyledonous and monocotyledonous species, respectively, and the genomic sequence data available was used to identify the lipase genes in wheat. Results In this present study, Ensembl Plants database was explored for lipase producing gene present in wheat genome and 21 genes were screened down as they contain specific domain and motif for lipase (GXSXG). According to the evolutionary analysis, it was found that the gene TraesCS5B02G157100, located in 5B chromosome, has 58.35% sequence similarity with the reported lipase gene of Arabidopsis thaliana and gene TraesCS3A02G463500 located in the 3A chromosome has 51.74% sequence similarity with the reported lipase gene of Oryza sativa. Homology modeling was performed using protein sequences coded by aforementioned genes and optimized by molecular dynamic simulations. Further with the help of molecular docking of modeled structures with tributyrin, binding efficiency was checked, and the difference in energies (DE) was −9.83 kcal/mol and −6.67 kcal/mol, respectively. Conclusions The present work provides a basic understanding of the gene-encoding lipase in wheat, which could be easily accessible and used as a potent industrial enzyme. The study enlightens another direction which can be used further to explore plant lipases.

Aromatic aldehydes are important industrial raw materials mainly synthesized by anti-Markovnikov (AM) oxidation of corresponding aromatic olefins. The AM product selectivity remains a big challenge. P450 aMOx is the first reported enzyme that could catalyze AM oxidation of aromatic olefins. Here, we reported a rational design strategy based on the “butterfly” model of the active site of P450 aMOx. Constrained molecular dynamic simulations and a binding energy analysis of key residuals combined with an experimental alanine scan were applied. As a result, the mutant A275G showed high AM selectivity of >99%. The results also proved that the “butterfly” model is an effective design strategy for enzymes.

In this work, it was proposed to employ the particle dispersion velocity (νD), pressure, heat release rate, delay time, and the numbers of Al-clusters and Al–O bonds derived from reactive molecular dynamics simulations to atomically identify the microexplosion of Al nanoparticles (ANPs), which is feasible to reduce agglomeration and promote the combustion efficiency in applications. The microexplosion is significantly different from common oxidation in terms of the evolution of morphology, reactivity, and energy, since it can be strengthened with an increase in O2 content and temperature and a decrease in particle size. That is, the microexplosion causes the more rapid increase in νD, pressure, heat release, and the number of Al-clusters and Al–O bonds, and the decrease in delay time compared to common oxidation. Meanwhile, the methods for assessing the intensity of the microexplosion and the promotion effect of the microexplosion on ANP oxidation were successfully established in this work. Importantly, the promotion efficiency was ascertained using the heat release rate (0.05 ​kJ·g−1·ps−1) and the ratio of the number of Al–O bonds to the number of Al atoms at the early stage of oxidation (0–50 ps). This work is expected to present a comprehensive perspective of the oxidation of ANPs in O2 under different O2 content, temperatures, and particle sizes, which may be instructive for other violent oxidation of various metal-based fuel nanoparticles.

Protein kinase p38γ is an attractive target against cancer because it plays a pivotal role in cancer cell proliferation by phosphorylating the retinoblastoma tumour suppressor protein. Therefore, inhibition of p38γ with active small molecules represents an attractive alternative for developing anti-cancer drugs. In this work, we present a rigorous and systematic virtual screening framework to identify potential p38γ inhibitors against cancer. We combined the use of machine learning-based quantitative structure activity relationship modelling with conventional computer-aided drug discovery techniques, namely molecular docking and ligand-based methods, to identify potential p38γ inhibitors. The hit compounds were filtered using negative design techniques and then assessed for their binding stability with p38γ through molecular dynamics simulations. To this end, we identified a promising compound that inhibits p38γ activity at nanomolar concentrations and hepatocellular carcinoma cell growth in vitro in the low micromolar range. This hit compound could serve as a potential scaffold for further development of a potent p38γ inhibitor against cancer.

A growing body of evidence suggests that only a few amino acids (“hot-spots”) at the interface contribute most of the binding energy in transient protein-protein interactions. However, experimental protocols to identify these hot-spots are highly labor-intensive and expensive. Computational methods, including evolutionary couplings, have been proposed to predict the hot-spots, but they generally fail to provide details of the interacting amino acids. Here we showed that unbiased evolutionary methods followed by biased molecular dynamic simulations could achieve this goal and reveal critical elements of protein complexes. We applied the methodology to selected G-protein coupled receptors (GPCRs), known for their therapeutic properties. We used the structure-prior-assisted direct coupling analysis (SP-DCA) to predict the binding interfaces of A2aR/D2R, CB1R/D2R, A2aR/CB1R, 5HT2AR/D2R, and 5-HT2AR/mGluR2 receptor heterodimers, which all agreed with published data. In order to highlight details of the interactions, we performed molecular dynamic (MD) simulations using the newly developed AWSEM energy model. We found that these receptors interact primarily through critical residues at the C and N terminal domains and the third intracellular loop (ICL3). The MD simulations showed that these residues are energetically necessary for dimerization and revealed their native conformational state. We subsequently applied the methodology to the 5-HT2AR/5-HTR4R heterodimer, given its implication in drug addiction and neurodegenerative pathologies such as Alzheimer’s disease (AD). Further, the SP-DCA analysis showed that 5-HT2AR and 5-HTR4R heterodimerize through the C-terminal domain of 5-HT2AR and ICL3 of 5-HT4R. However, elucidating the details of GPCR interactions would accelerate the discovery of druggable sites and improve our knowledge of the etiology of common diseases, including AD.

Candida dubliniensis is an opportunistic pathogen associated with oral and invasive fungal infections in immune-compromised individuals. Furthermore, the emergence of C. dubliniensis antifungal drug resistance could exacerbate its treatment. Hence, in this study a multi-epitope vaccine candidate has been designed using an immunoinformatics approach by targeting C. dubliniensis secreted aspartyl proteinases (SAP) proteins. In silico tools have been utilized to predict epitopes and determine their allergic potential, antigenic potential, toxicity, and potential to elicit interleukin-2 (IL2), interleukin-4 (IL4), and IFN-γ. Using the computational tools, eight epitopes have been predicted that were then linked with adjuvants for final vaccine candidate development. Computational immune simulation has depicted that the immunogen designed emerges as a strong immunogenic candidate for a vaccine. Further, molecular docking and molecular dynamics simulation analyses revealed stable interactions between the vaccine candidate and the human toll-like receptor 5 (TLR5). Finally, immune simulations corroborated the promising candidature of the designed vaccine, thus calling for further in vivo investigation.

The non-estrogenic 2,5-disubstituted tetrazole core-bearing bisphenol structures (TbB) are being researched as emerging structural congeners of Bisphenol A, an established industrial endocrine disruptor. However, there is no understanding of TbB’s adverse effects elicited via metabolic activation. Therefore, the current study aimed to investigate the metabolism of TbB ligands, with in silico results serving as a guide for in vitro studies. The Cytochrome P450 enzymes (CYP) inhibitory assay of TbB ligands on the seven human liver CYP isoforms (i.e., 1A2, 2A6, 2D6, 2C9, 2C8, 2C19, and 3A4) using human liver microsomes (HLM) revealed TbB ligand 223-3 to have a 50% inhibitory effect on all the CYP isoforms at a 10 μM concentration, except 1A2. The TbB ligand 223-10 inhibited 2B6 and 2C8, whereas the TbB ligand 223-2 inhibited only 2C9. The first-order inactivity rate constant (Kobs) studies indicated TbB ligands 223-3, 223-10 to be time-dependent (TD) inhibitors, whereas the TbB 223-2 ligand did not show such a significant effect. The 223-3 exhibited a TD inhibition for 2C9, 2C19, and 1A2 with Kobs values of 0.0748, 0.0306, and 0.0333 min−1, respectively. On the other hand, the TbB ligand 223-10 inhibited 2C9 in a TD inhibition manner with Kobs value 0.0748 min−1. However, the TbB ligand 223-2 showed no significant TD inhibition effect on the CYPs. The 223-2 ligand biotransformation pathway by in vitro studies in cryopreserved human hepatocytes suggested the clearance via glucuronidation with the predominant detection of only 223-2 derived mono glucuronide as a potential inactive metabolite. The present study demonstrated that the 223-2 ligand did not elicit any significant adverse effect via metabolic activation, thus paving the way for its in vivo drug–drug interactions (DDI) studies.

The neuropeptide vasopressin (VP) and its three G protein-coupled receptors (V1aR, V1bR and V2R) are of high interest in a wide array of drug discovery programs. V1aR is of particular importance due to its cardiovascular functions and diverse roles in the central nervous system. The structure–activity relationships underpinning ligand-receptor interactions remain however largely unclear, hindering rational drug design. This is not least due to the high structural flexibility of VP in its free as well as receptor-bound states. In this work, we developed a novel approach to reveal features of conformational selectivity upon VP-V1aR complex formation. We employed virtual screening strategies to probe VP’s conformational space for transiently adopted structures that favor binding to V1aR. To this end, we dissected the VP conformational space into three sub-ensembles, each containing distinct structural sets for VP’s three-residue C-terminal tail. We validated the computational results with experimental nuclear magnetic resonance (NMR) data and docked each sub-ensemble to V1aR. We observed that the conformation of VP’s three-residue tail significantly modulated the complex dissociation constants. Solvent-exposed and proline trans-configured VP tail conformations bound to the receptor with three-fold enhanced affinities compared to compacted or cis-configured conformations. The solvent-exposed and more flexible structures facilitated unique interaction patterns between VP and V1aR transmembrane helices 3, 4, and 6 which led to high binding energies. The presented “virtual conformational space screening” approach, integrated with NMR spectroscopy, thus enabled identification and characterization of a conformational selection-type complex formation mechanism that confers novel perspectives on targeting the VP-V1aR interactions at the level of the encounter complex – an aspect that opens novel research avenues for understanding the functionality of the evolutionary selected conformational properties of VP, as well as guidance for ligand design strategies to provide more potent and selective VP analogues.

Lysosomal acid lipase (LAL) is a serine hydrolase that hydrolyzes cholesteryl ester (CE) and TGs delivered to the lysosomes into free cholesterol and fatty acids. LAL deficiency due to mutations in the LAL gene (LIPA) results in accumulation of TGs and cholesterol esters in various tissues of the body leading to pathological conditions such as Wolman's disease and CE storage disease (CESD). Here, we present the first crystal structure of recombinant human LAL (HLAL) to 2.6 Å resolution in its closed form. The crystal structure was enabled by mutating three of the six potential glycosylation sites. The overall structure of HLAL closely resembles that of the evolutionarily related human gastric lipase (HGL). It consists of a core domain belonging to the classical α/β hydrolase-fold family with a classical catalytic triad (Ser-153, His-353, Asp-324), an oxyanion hole, and a “cap” domain, which regulates substrate entry to the catalytic site. Most significant structural differences between HLAL and HGL exist at the lid region. Deletion of the short helix, 238NLCFLLC244, at the lid region implied a possible role in regulating the highly hydrophobic substrate binding site from self-oligomerization during interfacial activation. We also performed molecular dynamic simulations of dog gastric lipase (lid-open form) and HLAL to gain insights and speculated a possible role of the human mutant, H274Y, leading to CESD.