Identification of a novel binding mechanism of Quinoline based molecules with lactate dehydrogenase of Plasmodium falciparum

Tetrahydropyridine appended 8-aminoquinoline derivatives: Design, synthesis, in silico, and in vitro antimalarial studies

Antimalarial drug resistance is a major obstacle in the ongoing quest against malaria. The disease affects half of the world's population. The majority of them are toddlers and pregnant women. Needed a potent compound to act on drug-resistant Pf at appropriate concentrations without endangering the host. Envisaged solving this issue through rational drug design by creating a novel hybrid drug possessing two pharmacophores that can act on two marvellous and independent aims within the cell. Synthesized a new series of substituted 4-phenyl-1,2,3,6-tetrahydropyridine (THP) 8-Aminoquinoline-based hybrid analogs which have been integrated with quinoline, chloroquine, pamaquine, and primaquine, which exhibited antimalarial activity against Pf. Out of thirteen 4-phenyl-1,2,3,6-THP appended 8-Aminoquinoline derivatives, the compounds 1j, 1e, 1b, and 1l have exhibited good antimalarial activity against chloroquine-sensitive (3D7) and chloroquine-resistant (RKL-9) strain with the minimum inhibitory concentration. Compound 1b was the most effective and showed consistently good potency against the drug-resistant (RKL-9) strain, although all other arrays showed good antimalarial efficacy. Additional docking and molecular dynamics studies were carried out at several targeting sites to quantify the structural parameters necessary for the activity.

Novel molecular inhibitor design for Plasmodium falciparum Lactate dehydrogenase enzyme using machine learning generated library of diverse compounds

Generative machine learning models offer a novel strategy for chemogenomics and de novo drug design, allowing researchers to streamline their exploration of the chemical space and concentrate on specific regions of interest. In cases with limited inhibitor data available for the target of interest, de novo drug design plays a crucial role. In this study, we utilized a package called 'mollib,' trained on ChEMBL data containing approximately 365,000 bioactive molecules. By leveraging transfer learning techniques with this package, we generated a series of compounds, starting from five initial compounds, which are potential Plasmodium falciparum (Pf) Lactate dehydrogenase inhibitors. The resulting compounds exhibit structural diversity and hold promise as potential novel Pf Lactate dehydrogenase inhibitors.

Antiplasmodial potential of isolated xanthones from Mesua ferrea Linn. roots: an in vitro and in silico molecular docking and pharmacokinetics study

BACKGROUND

Malaria is a major global health concern, particularly in tropical and subtropical countries. With growing resistance to first-line treatment with artemisinin, there is an urgent need to discover novel antimalarial drugs. Mesua ferrea Linn., a plant used in traditional medicine for various purposes, has previously been investigated by our research group for its cytotoxic properties. The objective of this study was to explore the compounds isolated from M. ferrea with regards to their potential antiplasmodial activity, their interaction with Plasmodium falciparum lactate dehydrogenase (PfLDH), a crucial enzyme for parasite survival, and their pharmacokinetic and toxicity profiles.

METHODS

The isolated compounds were assessed for in vitro antiplasmodial activity against a multidrug-resistant strain of P. falciparum K1 using a parasite lactate dehydrogenase (pLDH) assay. In vitro cytotoxicity against Vero cells was determined using the MTT (3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide) assay. The interactions between the isolated compounds and the target enzyme PfLDH were investigated using molecular docking. Additionally, pharmacokinetic and toxicity properties were estimated using online web tools SwissADME and ProTox-II, respectively.

RESULTS

Among the seven compounds isolated from M. ferrea roots, rheediachromenoxanthone (5), which belongs to the pyranoxanthone class, demonstrated good in vitro antiplasmodial activity, with the IC50 being 19.93 µM. Additionally, there was no toxicity towards Vero cells (CC50 = 112.34 µM) and a selectivity index (SI) of 5.64. Molecular docking analysis revealed that compound (5) exhibited a strong binding affinity of - 8.6 kcal/mol towards PfLDH and was stabilized by forming hydrogen bonds with key amino acid residues, including ASP53, TYR85, and GLU122. Pharmacokinetic predictions indicated that compound (5) possessed favorable drug-like properties and desired pharmacokinetic characteristics. These include high absorption in the gastrointestinal tract, classification as a non-substrate of permeability glycoprotein (P-gp), non-inhibition of CYP2C19, ease of synthesis, a high predicted LD50 value of 4,000 mg/kg, and importantly, non-hepatotoxic, non-carcinogenic, and non-cytotoxic effects.

CONCLUSIONS

This study demonstrated that compounds isolated from M. ferrea exhibit activity against P. falciparum. Rheediachromenoxanthone has significant potential as a scaffold for the development of potent antimalarial drugs.

In-silico identification of deleterious non-synonymous SNPs of TBX1 gene: Functional and structural impact towards 22q11.2DS

The TBX1 gene plays a critical role in the development of 22q11.2 deletion syndrome (22q11.2DS), a complex genetic disorder associated with various phenotypic manifestations. In this study, we performed in-silico analysis to identify potentially deleterious non-synonymous single nucleotide polymorphisms (nsSNPs) within the TBX1 gene and evaluate their functional and structural impact on 22q11.2DS. A comprehensive analysis pipeline involving multiple computational tools was employed to predict the pathogenicity of nsSNPs. This study assessed protein stability and explored potential alterations in protein-protein interactions. The results revealed the rs751339103(C>A), rs780800634(G>A), rs1936727304(T>C), rs1223320618(G>A), rs1248532217(T>C), rs1294927055 (C>T), rs1331240435 (A>G, rs1601289406 (A>C), rs1936726164 (G>A), and rs911796187(G>A) with a high-risk potential for affecting protein function and stability. These nsSNPs were further analyzed for their impact on post-translational modifications and structural characteristics, indicating their potential disruption of molecular pathways associated with TBX1 and its interacting partners. These findings provide a foundation for further experimental studies and elucidation of potential therapeutic targets and personalized treatment approaches for individuals affected by 22q11.2DS.

Designing novel anti-plasmodial quinoline-furanone hybrids: computational insights, synthesis, and biological evaluation targeting Plasmodium falciparum lactate dehydrogenase

To combat resistance against current antimalarials, modifying key pharmacophores and exploring novel parasite-specific drug targets remained one of the key drug design strategies. The resistance to quinoline-based antimalarials arises often due to the efflux of the drug. Hence, the development of newer agents containing bulkier pharmacophores will enable medicinal chemists to counteract drug resistance. In view of this, herein we designed bulkier quinoline-furanone hybrids. Initially, virtual drug-likeness and ADMET screening were conducted to optimize physicochemical properties followed by docking of the hybrids against the Plasmodium falciparum lactate dehydrogenase (PfLDH) enzyme. The most potent hybrids that emerged from the computational screening were synthesized and screened for their bioactivity against the resistant strain of Plasmodium through Schizont Maturation Inhibition assays. Among the compounds tested, 5g and 6e demonstrated the best activity, with IC50 values similar to chloroquine (CQ), and 5g exhibited superior LDH inhibition compared to CQ. Compounds 5f, 7a, and 7f showed IC50 values comparable to CQ and moderate LDH inhibition. Structure-activity relationship (SAR) analysis revealed that halogen substitutions, particularly Br and Cl, enhanced antimalarial activity, while strong electron-withdrawing (-NO2) or -donating (-OH) groups led to diminished activity. Additionally, bulkier aromatic substitutions were favoured for antimalarial activity and LDH inhibition. The investigation successfully found potent anti-plasmodial quinoline-furanone hybrids, demonstrating promising prospects for combating malaria.

Binding patterns of inhibitors to different pockets of kinesin Eg5

The kinesin-5 family member, Eg5, plays very important role in the mitosis. As a mitotic protein, Eg5 is the target of various mitotic inhibitors. There are two targeting pockets in the motor domain of Eg5, which locates in the α2/L5/α3 region and the α4/α6 region respectively. We investigated the interactions between the different inhibitors and the two binding pockets of Eg5 by using all-atom molecular dynamics method. Combined the conformational analysis with the free-energy calculation, the binding patterns of inhibitors to the two binding pockets are shown. The α2/L5/α3 pocket can be divided into 4 regions. The structures and binding conformations of inhibitors in region 1 and 2 are highly conserved. The shape of α4/α6 pocket is alterable. The space of this pocket in ADP-binding state of Eg5 is larger than that in ADP·Pi-binding state due to the limitation of a hydrogen bond formed in the ADP·Pi-binding state. The results of this investigation provide the structural basis of the inhibitor-Eg5 interaction and offer a reference for the Eg5-targeted drug design.

Screening approach by a combination of computational and in vitro experiments: identification of fluvastatin sodium as a potential SIRT2 inhibitor

BACKGROUND

Sirtuins (SIRTs) are NAD+-dependent deacetylases that play various roles in numerous pathophysiological processes, holding promise as therapeutic targets worthy of further investigation. Among them, the SIRT2 subtype is closely associated with tumorigenesis and malignancies. Dysregulation of SIRT2 activation can regulate the expression levels of related genes in cancer cells, leading to tumor occurrence and metastasis.

METHODS

In this study, we used computer simulations to screen for novel SIRT2 inhibitors from the FDA database, based on which 10 compounds with high docking scores and good interactions were selected for in vitro anti-pancreatic cancer metastasis testing and enzyme binding inhibition experiments. The results showed that fluvastatin sodium may possess inhibitory activity against SIRT2. Subsequently, fluvastatin sodium was subjected to molecular docking experiments with various SIRT isoforms, and the combined results from Western blotting experiments indicated its potential as a SIRT2 inhibitor. Next, molecular docking, molecular dynamics (MD) simulations, and binding free energy calculations were performed, revealing the binding mode of fluvastatin sodium at the SIRT2 active site, further validating the stability and interaction of the ligand-protein complex under physiological conditions.

RESULTS

Overall, this study provides a systematic virtual screening workflow for the discovery of SIRT2 activity inhibitors, identifies the potential inhibitory effect of fluvastatin sodium as a lead compound on SIRT2, and opens up a new direction for developing highly active and selectively targeted SIRT2 inhibitors.

Systematic exploration of network pharmacology, in silico modeling and pharmacokinetic profiling for vitamin E in polycystic ovarian syndrome

Aim: This study seeks to explore the possibility of using vitamin E to alleviate the symptoms of polycystic ovarian syndrome (PCOS). Methods: Various computational methods were employed, including network pharmacology utilizing a compound-target-pathway approach, Swiss ADME, OSIRIS® property explorer, pkCSM, PASS online web resource and MOLINSPIRATION® software. In addition, in silico analysis of vitamin E was performed with ten receptors. Results & discussion: Our findings highlight the diverse potential of vitamin E in alleviating PCOS. The observed influence on hormones is in line with existing PCOS theories regarding cyst development, further enhancing the therapeutic promise of vitamin E. Conclusion: In conclusion, our computational analysis indicates that vitamin E shows potential as a therapeutic agent for alleviating PCOS in adolescents.

Molecular dynamic investigation for Roco4 kinase inhibitor as treatment options for parkinsonism

CONTEXT

Parkinson's disease is a neurodegenerative condition characterized by the degeneration of dopaminergic neurons, resulting in motor disabilities such as rigidity, bradykinesia, postural instability, and resting tremors. While the exact cause of Parkinson's remains uncertain, both familial and sporadic forms are often associated with the G2019S mutation found in the kinase domain of LRRK2. Roco4 is an analogue of LRRK2 protein in Dictyostelium discoideum which is an established model organism to investigate LRRK2 inhibitors. In this study, the potential treatment of Parkinson's was explored by inhibiting the activity of the mutated LRRK2 protein using Roco4 as the base protein structure. Mongolicain-A and Bacoside-A exhibited significant selectivity towards the G2019S mutation, displaying a binding affinity of - 12.3 Kcal/mol and - 11.4 Kcal/mol respectively. Mongolicain-A demonstrated increased specificity towards Roco4, while Bacoside-A demonstrated significant binding affinity to all 34 kinases proteins alike. The Molecular Dynamics Studies (MDS) results strongly suggests that Mongolicain-A is a significant inhibitor of Roco4 kinase. ADMET and drugability analysis also suggests that among the two best ligands, Mongolicain-A demonstrates significant physicochemical properties to be suitable for best drug like molecule. Based on the in-silico molecular docking, molecular dynamic simulation, ADMET and drugability analyses, it is strongly suggested that, Mongolicain-A could be a potential candidate for treatment and management of Parkinson's disease via inhibition of LRRK2 protein. Further in-vitro and in-vivo investigations are in demand to validate these findings.

METHODS

To identify potential inhibitors, 3069 phytochemicals were screened using molecular docking via AutoDock Vina. Molecular Dynamics Simulation was carried out using GROMACS 2022.2 for a duration of 100ns per complex to study the stability and inhibition potential of the protein ligand complexes. ADMET analysis was carriedout using Molinspiration and preADMET web tool.

Molecular dynamics simulation of the effect of temperature on the conformation of ubiquitin protein

CONTENT

Ubiquitin, a ubiquitous small protein found in all living organisms, is crucial for tagging proteins earmarked for degradation and holds pivotal importance in biomedicine. Protein functionality is intricately linked to its structure. To comprehend the impact of diverse temperatures on ubiquitin protein structure, our study delved into the energy landscape, hydrogen bonding, and overall structural stability of ubiquitin protein at varying temperatures. Through meticulous analysis of root mean square deviation and root mean square fluctuation, we validated the robustness of the simulation conditions employed. Within our simulated system, the bonding energy and electrostatic potential energy exhibited linear augmentation, while the van der Waals energy demonstrated a linear decline. Additionally, our findings highlighted that the α-Helix secondary structure of the ubiquitin protein gradually transitions toward helix destabilization under high-temperature conditions. The secondary structure of ubiquitin protein experiences distinct changes under varying temperatures. The outcomes of our molecular simulations offer a theoretical framework that enhances our comprehension of how temperature impacts the structural stability of ubiquitin protein. These insights contribute not only to a deeper understanding of iniquity's behavior but also hold broader implications in the realm of biomedicine and beyond.

METHODS

All the MD simulations were performed using the GROMACS software with GROMOS96 force field and SPC for water. The ubiquitin protein was put in the center of a cubic box with a length of 8 nm, a setting that allowed > 0.8 nm in the minimal distance between the protein surface and the box wall. To remove the possible coordinate collision of the configurations, in the beginning, the steepest descent method was used until the maximum force between atoms was under 100 kJ/mol·nm with a 0.01 nm step size. Minimization was followed by 30 ps of position-restrained MD simulation. The protein was restrained to its initial position, and the solvent was freely equilibrated. The product phase was obtained with the whole system simulated for 10 ns without any restraint using an integral time step of 1 fs with different temperatures. The cutoff for short-range electronic interaction was set to 1.5 nm. The long-range interactions were treated with a particle-mesh Ewald (PME) method with a grid width of 1.2 nm.

Discovery of a novel SHP2 allosteric inhibitor using virtual screening, FMO calculation, and molecular dynamic simulation

CONTEXT

SHP2 is a non-receptor protein tyrosine phosphatase to remove tyrosine phosphorylation. Functionally, SHP2 is an essential bridge to connect numerous oncogenic cell-signaling cascades including RAS-ERK, PI3K-AKT, JAK-STAT, and PD-1/PD-L1 pathways. This study aims to discover novel and potent SHP2 inhibitors using a hierarchical structure-based virtual screening strategy that combines molecular docking and the fragment molecular orbital method (FMO) for calculating binding affinity (referred to as the Dock-FMO protocol). For the SHP2 target, the FMO method prediction has a high correlation between the binding affinity of the protein-ligand interaction and experimental values (R2 = 0.55), demonstrating a significant advantage over the MM/PBSA (R2 = 0.02) and MM/GBSA (R2 = 0.15) methods. Therefore, we employed Dock-FMO virtual screening of ChemDiv database of ∼2,990,000 compounds to identify a novel SHP2 allosteric inhibitor bearing hydroxyimino acetamide scaffold. Experimental validation demonstrated that the new compound (E)-2-(hydroxyimino)-2-phenyl-N-(piperidin-4-ylmethyl)acetamide (7188-0011) effectively inhibited SHP2 in a dose-dependent manner. Molecular dynamics (MD) simulation analysis revealed the binding stability of compound 7188-0011 and the SHP2 protein, along with the key interacting residues in the allosteric binding site. Overall, our work has identified a novel and promising allosteric inhibitor that targets SHP2, providing a new starting point for further optimization to develop more potent inhibitors.

METHODS

All the molecular docking studies were employed to identify potential leads with Maestro v10.1. The protein-ligand binding affinities of potential leads were further predicted by FMO calculations at MP2/6-31G* level using GAMESS v2020 system. MD simulations were carried out with AmberTools18 by applying the FF14SB force field. MD trajectories were analyzed using VMD v1.9.3. MM/GB(PB)SA binding free energy analysis was carried out with the mmpbsa.py tool of AmberTools18. The docking and MD simulation results were visualized through PyMOL v2.5.0.

In-silico lead identification of the pan-mutant IDH1 and IDH2 inhibitors to target glioblastoma

Glioblastoma (GBM) is an aggressive malignant type of brain tumor. Targeting one single intracellular pathway might not alleviate the disease, rather it activates the other molecular pathways that lead to the worsening of the disease condition. Therefore, in this study, we attempted to target both isocitrate dehydrogenase 1 (IDH1) and IDH2, which are one of the most commonly mutated proteins in GBM and other cancer types. Here, standard precision and extra precision docking, IFD, MM-GBSA, QikProp, and molecular dynamics (MD) simulation were performed to identify the potential dual inhibitor for IDH1 and IDH2 from the enamine database containing 59,161 ligands. Upon docking the ligands with IDH1 (PDB: 6VEI) and IDH2 (PDB: 6VFZ), the top eight ligands were selected, based on the XP Glide score. These ligands produced favourable MMGBSA scores and ADME characteristics. Finally, the top four ligands 12953, 44825, 51295, and 53210 were subjected to MD analysis. Interestingly, 53210 showed maximum interaction with Gln 277 for 99% in IDH1 and Gln 316 for 100% in IDH2, which are the crucial amino acids for the inhibitory function of IDH1 and IDH2 to target GBM. Therefore, the present study attempts to identify the novel molecules which could possess a pan-inhibitory action on both IDH1 and IDH that could be crucial in the management of GBM. Yet further evaluation involving in vitro and in vivo studies is warranted to support the data in our current study.Communicated by Ramaswamy H. Sarma.

Gallic Acid-Modified Polyethylenimine-Polypropylene Carbonate-Polyethylenimine Nanoparticles: Synthesis, Characterization, and Anti-Periodontitis Evaluation

The aim of the research was to develop novel gallic acid (GA)-modified amphiphilic nanoparticles of polyethylenimine (PEI)-polypropylene carbonate (PPC)-PEI (PEPE) and comprehensively assess its properties as an antiperiodontitis nanoparticle targeting the Toll-like receptor (TLR). The first step is to evaluate the binding potential of GA to the core trigger receptors TLR2 and TLR4/MD2 for periodontitis using molecular docking techniques. Following this, we conducted NMR, transmission electron microscopy, and dynamic light scattering analyses on the synthesized PEPE nanoparticles. As the final step, we investigated the synthetic results and in vitro antiperiodontitis properties of GA-PEPE nanoparticles. The investigation revealed that GA exhibits potential for targeted binding to TLR2 and the TLR4/MD2 complex. Furthermore, we successfully developed 91.19 nm positively charged PEPE nanoparticles. Spectroscopic analysis indicated the successful synthesis of GA-modified PEPE. Additionally, CCK8 results demonstrated that GA modification significantly reduced the biotoxicity of PEPE. The in vitro antiperiodontitis properties assessment illustrated that 6.25 μM of GA-PEPE nanoparticles significantly reduced the expression of pro-inflammatory factors TNF-α, IL-1β, and IL-6. The GA-PEPE nanoparticles, with their targeted TLR binding capabilities, were found to possess excellent biocompatibility and antiperiodontitis properties. GA-PEPE nanoparticles will provide highly innovative input into the development of anti- periodontitis nanoparticles.

Exosomes derived from human umbilical cord mesenchymal stem cells decrease neuroinflammation and facilitate the restoration of nerve function in rats suffering from intracerebral hemorrhage

Exosomes derived from human umbilical cord mesenchymal stem cells (hUCMSC-ex) have become a hopeful substitute for whole-cell therapy due to their minimal immunogenicity and tumorigenicity. The present study aimed to investigate the hypothesis that hUCMSC-ex can alleviate excessive inflammation resulting from intracerebral hemorrhage (ICH) and facilitate the rehabilitation of the nervous system in rats. In vivo, hemorrhagic stroke was induced by injecting collagenase IV into the striatum of rats using stereotactic techniques. hUCMSC-ex were injected via the tail vein at 6 h after ICH model establishment at a dosage of 200 µg. In vitro, astrocytes were pretreated with hUCMSC-ex and then stimulated with hemin (20 μmol/mL) to establish an ICH cell model. The expression of TLR4/NF-κB signaling pathway proteins and inflammatory factors, including TNF-α, IL-1β, and IL-10, was assessed both in vivo and in vitro to investigate the impact of hUCMSC-ex on inflammation. The neurological function of the ICH rats was evaluated using the corner turn test, forelimb placement test, Longa score, and Bederson score on the 1st, 3rd, and 5th day. Additionally, RT-PCR was employed to examine the mRNA expression of TLR4 following hUCMSC-ex treatment. The findings demonstrated that hUCMSC-ex downregulated the protein expression of TLR4, NF-κB/P65, and p-P65, reduced the levels of pro-inflammatory cytokines TNF-α and IL-1β, and increased the expression of the anti-inflammatory cytokine IL-10. Ultimately, the administration of hUCMSC-ex improved the behavioral performance of the ICH rats. However, the results of PT-PCR indicated that hUCMSC-ex did not affect the expression of TLR4 mRNA induced by ICH, suggesting that hUCMSCs-ex may inhibit TLR4 translation rather than transcription, thereby suppressing the TLR4/NF-κB signaling pathway. We can conclude that hUCMSC-ex mitigates hyperinflammation following ICH by inhibiting the TLR4/NF-κB signaling pathway. This study provides preclinical evidence for the potential future application of hUCMSC-ex in the treatment of cerebral injury.

A discovery of potent kaempferol derivatives as multi-target medicines against diabetes as well as bacterial infections: an in silico approach

Flavonoids demonstrate beneficial effects on human health because flavonoids contain important biological properties. Kaempferol is a flavonol, type of flavonoid found in eatable plants and in plants usually employed in ancient drugs (Moringa oleifera, Tilia spp., fern genus spp. and gingko etc.). Some medicinal studies have shown that the use of foods full of kaempferol decreases the risk of many (cancer, vascular) diseases. All the data of 50 kaempferol derivatives were collected from PubChem database. Through Schrödinger software, 3D-QSAR study was performed for 50 compounds by using method of field base. Conformer of kaempferol derivatives was docked against anti-diabetic, anti-microbial co-crystal structures and protein. To monitor the best anti-diabetic and antibacterial agent, particular kaempferol derivatives were downloaded from PubChem database. Virtual screening by molecular docking provided four lead compounds with four different proteins. These hit compounds were found to be potent inhibitor for diabetic enzymes alpha-amylase and DPP IV and had the potential to suppress DNA gyrase and dihydrofolate reductase synthesis. Molecular dynamic simulation of docked complexes evaluates the value of root mean square fluctuation by iMOD server. Kaempferol 3-O-alpha-L-(2, 3-di-Z-p-coumaroyl) rhamnoside (42) compound used as anti-diabetic and kaempferol 3-O-gentiobioside (3) as antibacterial with good results can be used for drug discovery.Communicated by Ramaswamy H. Sarma.

Genetic analysis of albinism caused by compound heterozygous mutations of the OCA2 gene in a Chinese family

BACKGROUND

Oculocutaneous albinism (OCA) is a group of rare genetic disorders characterized by a reduced or complete lack of melanin in the skin, hair, and eyes. Patients present with colorless retina, pale pink iris, and pupil, and fear of light. The skin, eyebrows, hair, and other body hair are white or yellowish-white. These conditions are caused by mutations in specific genes necessary for the production of melanin. OCA is divided into eight clinical types (OCA1-8), each with different clinical phenotypes and potential genetic factors. This study aimed to identify the genetic causes of non-syndromic OCA in a Chinese Han family.

METHODS

We performed a comprehensive clinical examination of family members, screened for mutation loci using whole exome sequencing (WES) technology, and predicted mutations using In silico tools.

RESULTS

The patient's clinical manifestations were white skin, yellow hair, a few freckles on the cheeks and bridge of the nose, decreased vision, blue iris, poorly defined optic disk borders, pigmentation of the fundus being insufficient, and significant vascular exposure. The WES test results indicate that the patient has compound heterozygous mutations in the OCA2 gene (c.1258G > A (p.G420R), c.1441G > A (p.A481T), and c.2267-2 A > C), respectively, originating from her parents. Among them, c.1258G > A (p.G420R) is a de novo mutation with pathogenic. Our analysis suggests that compound heterozygous mutations in the OCA2 gene are the primary cause of the disease in this patient.

CONCLUSIONS

The widespread application of next-generation sequencing technologies such as WES in clinical practice can effectively replace conventional detection methods and assist in the diagnosis of clinical diseases more quickly and accurately. The newly discovered c.1258G > A (p.G420R) mutation can update and expand the gene mutation spectrum of OCA2-type albinism.

Structure-based virtual screening and molecular dynamics approaches to identify new inhibitors of Staphylococcus aureus sortase A

Staphylococcus aureus is a prevalent Gram-positive bacteria leading cause of a wide range of human pathologies. Moreover, antibiotic résistance of pathogenesis bacteria is one of the worldwide health problems. In Gram-positive bacteria, the enzyme of SrtA, is responsible for the anchoring of surface-exposed proteins to the cell wall peptidoglycan. Because of its critical role in Gram-positive bacterial pathogenesis, SrtA is an attractive target for anti-virulence during drug development. To date, some SrtA inhibitors have been discovered most of them being derived from flavonoid compounds, like Myricetin. In order to provide potential hit molecules against SrtA for clinical use, we obtained a total of 293 compounds by performing in silico shape-based screening of compound libraries against Myristin as a reference structure. Employing molecular docking and scoring functions, the top 3 compounds Apigenin, Efloxate, and Compound 8261032 were screened by comparing their docking scores with Myricetin. Furthermore, MD simulations and MM-PBSA binding energy calculation studies revealed that only Compound 8261032 strongly binds to the catalytic core of the SrtA enzyme than Myricetin, and stable behavior was consistently observed in the docking complex. Compound 8261032 showed a good number of hydrogen bonds with SrtA and higher MM-PBSA binding energy when compared to all three molecules. Also, it makes strength interactions with Arg139 and His62, which are critical for SrtA biological activity. This study showed that the development of this inhibitor could be a fundamental strategy against resistant bacteria, but further studies in vitro are needed to confirm this claim.Communicated by Ramaswamy H. Sarma.

Discovery of Palbociclib as a potent c-Myc G4 stabilizer for lung cancer treatment using molecular docking, molecular dynamics simulation, and in vitro activity evaluation

Despite significant progress in lung cancer treatment, this disease remains a prevalent and serious global malignancy, leading to high rates of illness and death. Urgent research is needed to discover new or alternative therapies that can improve clinical outcomes for lung cancer patients. In our study, we successfully demonstrated the effectiveness of Palbociclib, a CDK4/6 inhibitor, in suppressing the growth of lung cancer cells. The IC50 values obtained were 11.00 μM and 11.74 μM for H1299 and A549 cells, respectively. Furthermore, our findings indicate that Palbociclib may possess strong c-Myc G4 stabilizing properties by significantly reducing both protein and mRNA expression levels of c-Myc. Additionally, Palbociclib induces apoptosis and causes cell cycle arrest at the G2/M phase in two cells. Through circular dichroism (CD), molecular docking, and molecular dynamics (MD) simulation, we have provided evidence that Palbociclib enhances the structural stability of c-Myc G4 while exhibiting a high binding affinity to its ligand's binding site on c-Myc G4. These results suggest that Palbociclib holds promise as a novel c-Myc G4 stabilizer for treating cancers associated with abnormal c-Myc activity; further optimization and development are warranted.

Cuproptosis gene-related, neural network-based prognosis prediction and drug-target prediction for KIRC

BACKGROUND

Kidney renal clear cell carcinoma (KIRC), as a common case in renal cell carcinoma (RCC), has the risk of postoperative recurrence, thus its prognosis is poor and its prognostic markers are usually based on imaging methods, which have the problem of low specificity. In addition, cuproptosis, as a novel mode of cell death, has been used as a biomarker to predict disease in many cancers in recent years, which also provides an important basis for prognostic prediction in KIRC. For postoperative patients with KIRC, an important means of preventing disease recurrence is pharmacological treatment, and thus matching the appropriate drug to the specific patient's target is also particularly important. With the development of neural networks, their predictive performance in the field of medical big data has surpassed that of traditional methods, and this also applies to the field of prognosis prediction and drug-target prediction.

OBJECTIVE

The purpose of this study is to screen for cuproptosis genes related to the prognosis of KIRC and to establish a deep neural network (DNN) model for patient risk prediction, while also developing a personalized nomogram model for predicting patient survival. In addition, sensitivity drugs for KIRC were screened, and a graph neural network (GNN) model was established to predict the targets of the drugs, in order to discover potential drug action sites and provide new treatment ideas for KIRC.

METHODS

We used the Cancer Genome Atlas (TCGA) database, International Cancer Genome Consortium (ICGC) database, and DrugBank database for our study. Differentially expressed genes (DEGs) were screened using TCGA data, and then a DNN-based risk prediction model was built and validated using ICGC data. Subsequently, the differences between high- and low-risk groups were analyzed and KIRC-sensitive drugs were screened, and finally a GNN model was trained using DrugBank data to predict the relevant targets of these drugs.

RESULTS

A prognostic model was built by screening 10 significantly different cuproptosis-related genes, the model had an AUC of 0.739 on the training set (TCGA data) and an AUC of 0.707 on the validation set (ICGC data), which demonstrated a good predictive performance. Based on the prognostic model in this paper, patients were also classified into high- and low-risk groups, and functional analyses were performed. In addition, 251 drugs were screened for sensitivity, and four drugs were ultimately found to have high sensitivity, with 5-Fluorouracil having the best inhibitory effect, and subsequently their corresponding targets were also predicted by GraphSAGE, with the most prominent targets including Cytochrome P450 2D6, UDP-glucuronosyltransferase 1A, and Proto-oncogene tyrosine-protein kinase receptor Ret. Notably, the average accuracy of GraphSAGE was 0.817 ± 0.013, which was higher than that of GAT and GTN.

CONCLUSION

Our KIRC risk prediction model, constructed using 10 cuproptosis-related genes, had good independent prognostic ability. In addition, we screened four highly sensitive drugs and predicted relevant targets for these four drugs that might treat KIRC. Finally, literature research revealed that four drug-target interactions have been demonstrated in previous studies and the remaining targets are potential sites of drug action for future research.

In-silico bioprospecting of secondary metabolites from endophytic Streptomyces spp. against Magnaporthe oryzae, a cereal killer fungus

Rice blast disease, caused by Magnaporthe oryzae, is the most devastating cereal killer worldwide. Note that melanin pigment is an essential factor of M. oryzae virulence, thus fungicides interfering with melanin biosynthesizing enzymes would reduce the pathogenicity. Scytalone dehydratase (SDH) is the key target for commercial fungicides, like carpropamid, due to its role in the dehydration reaction of the fungal melanin pathway. However, a single-point mutation (V75M) in SDH elicits resistance to carpropamid. A lack of effective fungicides against this resistant strain expedited the quest for novel bioactive inhibitors. Currently, bacterial endophytes like Streptomyces have been heralded for synthesizing bioactive metabolites to protect plants from phytopathogens. The literature search led to the identification of 21 Streptomyces spp. symbionts of paddy that can suppress M. oryzae growth. An antiSMASH server was used to explore Streptomyces spp. gene clusters and found 4463 putative metabolites. Besides, 745 unique metabolites were subjected to a series of virtual screening techniques. Ideally, this process identified five potential SDH inhibitors. The docking result highlights that the metabolite pseudopyronine A interacted hydrophobically with both Val75 of SDHWT and Met75 of SDHV75M targets. Moreover, pseudopyronine A has a higher binding free energy with SDHWT (- 89.94 kcal/mol) and SDHV75M (- 71.95 kcal/mol). Interestingly, the pyranones scaffold of pseudopyronine A was reported for antifungal activity against phytopathogens. Dynamic behavior confirms that pseudopyronine A has excellent conformational states with both SDHWT and SDHV75M. Altogether, we hope that this study creates a new avenue for the discovery of novel phytopathogen inhibitors from endophytes.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03859-7.

In Vitro Antioxidant and In Vivo Hepatoprotective Properties of Wissadula periplocifolia Extract

Wissadula periplocifolia (L.) Thwaites is a traditional medicinal plant belonging to the family Malvaceae, used in folk medicine for inflamed snake bites and bee stings. The current study was designed to investigate the in vitro antioxidant and in vivo anti-inflammatory and hepatoprotective activities of 80% ethanol extract of W. periplocifolia and its different fractions. The crude ethanolic extract (CEE) was then serially fractionated with petroleum ether fraction (PEF), chloroform fraction (CHF), and aqueous fraction (AQF). The antioxidant activity was assessed using 1,1-diphenyl-2-picrylhydrazyl (DPPH) radical assay, anti-inflammatory activity was determined in the xylene-induced ear edema model, and hepatoprotective activity was measured in the paracetamol-induced hepatic injury model. PEF showed a significant scavenging effect with an IC50 value of 33.5 μg/mL, followed by CEE (IC50 = 42.2 μg/mL), CHF (IC50 = 77 μg/mL), and AQF (IC50 = 80 μg/mL), compared to standard butylated hydroxytoluene (IC50 = 14.8 μg/mL). Both doses of CEE (250 and 500 mg/kg) could reduce ear edema by 41.3 and 50%, respectively, compared to standard diclofenac sodium (76.09%). Moreover, CEE significantly reduces the elevated liver enzymes (ALT, AST, and ALP), compared to control. Nevertheless, it elevated blood protein and reduced the blood bilirubin level (p < 0.01), compared to control. Histopathological studies also indicated significant protection of the liver from paracetamol-induced liver damage. In conclusion, W. periplocifolia could be a good source of antioxidant and hepatoprotective phytochemicals; meanwhile, toxicological and pharmacokinetic studies are recommended.

Nearly complete structure of bacteriophage DT57C reveals architecture of head-to-tail interface and lateral tail fibers

The T5 family of viruses are tailed bacteriophages characterized by a long non-contractile tail. The bacteriophage DT57C is closely related to the paradigmal T5 phage, though it recognizes a different receptor (BtuB) and features highly divergent lateral tail fibers (LTF). Considerable portions of T5-like phages remain structurally uncharacterized. Here, we present the structure of DT57C determined by cryo-EM, and an atomic model of the virus, which was further explored using all-atom molecular dynamics simulations. The structure revealed a unique way of LTF attachment assisted by a dodecameric collar protein LtfC, and an unusual composition of the phage neck constructed of three protein rings. The tape measure protein (TMP) is organized within the tail tube in a three-stranded parallel α-helical coiled coil which makes direct contact with the genomic DNA. The presence of the C-terminal fragment of the TMP that remains within the tail tip suggests that the tail tip complex returns to its original state after DNA ejection. Our results provide a complete atomic structure of a T5-like phage, provide insights into the process of DNA ejection as well as a structural basis for the design of engineered phages and future mechanistic studies.

High-throughput virtual screening of potential inhibitors of GPR52 using docking and biased sampling method for Huntington's disease therapy

Huntington's disease (HD) is a rare and progressive neurodegenerative disorder caused by polyglutamine (poly-Q) mutations of the huntingtin (HTT) gene resulting in chorea, cognitive, and psychiatric dysfunctions. Being a monogenic condition, reducing the levels of the mutated huntingtin protein (mHTT) holds promise as an effective therapeutic approach. GPR52, an orphan G-protein coupled receptor (GPCR), enriched in the striatum, is a novel target for slowing down the progression of HD by lowering the mHTT levels. Therefore, the study focuses on identifying potent small-molecule inhibitors for GPR52 using a combination of robust high-throughput virtual screening (HTVS) and pharmacokinetics profiling followed by fast pulling of ligand (FPL) and umbrella sampling (US) simulations. Initially, screening a library of 2,36,545 compounds was done against the binding pocket of GPR52. Based on binding affinity, stereochemical and non-bonded interactions, and pharmacokinetic profiling, 50 compounds were shortlisted. Selected hit compounds 1, 2, and 3 were subjected to FPL simulations with applied external bias potential to investigate their unique dissociation pathways and intermolecular interactions over time. Subsequently, the US simulations were performed on the selected hit compounds to estimate their binding free energy (ΔG). The analysis of the trajectories obtained from simulations revealed that the residues TYR34, TYR185, GLY187, ASP188, ILE189, SER299, PHE300, and THR303 within the active site of GPR52 were significant for efficient ligand binding through the formation of various hydrogen bond interactions and hydrophobic contacts. Out of the three hit compounds, compound 3 had the lowest ΔG of - 20.82 ± 0.44 kcal/mol. The study identified compounds 1, 2, and 3 as potential molecules that can be developed as GPR52 inhibitors holding promise for lowering mHTT levels.

In Silico Studies of Four Compounds of Cecropia obtusifolia against Malaria Parasite

Malaria is a disease that affects many people in the world. In Mexico, malaria remains an active disease in certain regions, particularly in the states of Chiapas and Chihuahua. While antimalarial effects have been attributed to some species of Cecropia in various countries, no such studies have been conducted in Mexico. Therefore, the objective of this study was to evaluate the in silico antimalarial activity of some active compounds identified according to the literature in the species of Cecropia obtusifolia, belonging to the Cecropiaceae family, such as ursolic acid, α-amyrin, chrysin, and isoorientin. These compounds were evaluated with specific molecular docking and molecular dynamics (MD) studies using three different malarial targets with the PDB codes 1CET, 2BL9, and 4ZL4 as well as the prediction of their pharmacokinetic (Pk) properties. Docking analysis revealed the following best binding energies (kcal/mol): isoorientin-1CET (-9.1), isoorientin-2BL9 (-8.8), and chrysin-4ZL4 (-9.6). MD simulation validated the stability of the complexes. Pharmacokinetics analysis suggested that the compounds would generally perform well if administered. Therefore, these results suggest that these compounds may be used as potential drugs for the treatment of malaria.

Molecular Drug Simulation and Experimental Validation of the CD36 Receptor Competitively Binding to Long-Chain Fatty Acids by 7-Ketocholesteryl-9-carboxynonanoate

Long-chain fatty acids (LCFAs) are one of the main energy-supplying substances in the body. LCFAs with different lengths and saturations may have contrasting biological effects that exacerbate or alleviate progress against a variety of systemic disorders of lipid metabolism in organisms. Nonalcoholic fatty liver disease is characterized by chronic inflammation and steatosis, mainly caused by the ectopic accumulation of lipids in the liver, especially LCFAs. CD36 is a scavenger receptor that recognizes and mediates the transmembrane absorption of LCFAs and is expressed in a variety of cells throughout the body. In previous studies, our group found that 7-ketocholesteryl-9-carboxynonanoate (oxLig-1) has the biological effect of targeting CD36 to inhibit oxidized low-density lipoprotein lipotoxicity-induced lipid metabolism disorder; it has an ω-carboxyl physiologically active center and is structurally similar to LCFAs. However, the biological mechanism of oxLig-1 binding to CD36 and competing for binding to different types of LCFAs is still not clear. In this study, molecular docking and molecular dynamics simulation were utilized to simulate and analyze the binding activity between oxLig-1 and different types of LCFAs to CD36 and confirmed by the enzyme-linked immunosorbent assay (ELISA) method. Absorption, distribution, metabolism, excretion, and toxicity (ADMET) platform was applied to predict the drug-forming properties of oxLig-1, and HepG2 cells model of oleic acid and nonalcoholic fatty liver disease (NAFLD) model mice were validated to verify the biological protection of oxLig-1 on lipid lowering. The results showed that there was a co-binding site of LCFAs and oxLig-1 on CD36, and the binding driving forces were mainly hydrogen bonding and hydrophobic interactions. The binding abilities of polyunsaturated LCFAs, oxLig-1, monounsaturated LCFAs, and saturated LCFAs to CD36 showed a decreasing trend in this order. There was a similar decreasing trend in the stability of the molecular dynamics simulation. ELISA results similarly confirmed that the binding activity of oxLig-1 to CD36 was significantly higher than that of typical monounsaturated and saturated LCFAs. ADMET prediction results indicated that oxLig-1 had a good drug-forming property. HepG2 cells model of oleic acid and NAFLD model mice study results demonstrated the favorable lipid-lowering biological effects of oxLig-1. Therefore, oxLig-1 may have a protective effect by targeting CD36 to inhibit the excessive influx and deposition of lipotoxicity monounsaturated LCFAs and saturated LCFAs in hepatocytes.

Improvement of binding pose prediction of the MR1 covalent ligands by inclusion of simple pharmacophore constraints and structural waters in the docking process

The major histocompatibility complex (MHC) class I-related molecule, MR1, is a key component of the immune system, presenting antigens to T-cell receptors (TCRs) and modulating the immune response against various antigens. MR1 possesses a compact ligand-binding pocket despite its ability to interact with ligands that can have either agonistic or antagonistic effects on the immune system. Agonistic ligands can stimulate the immune response, while antagonistic ligands do not elicit an immune response. In most cases, ligand binding to MR1 is mediated through a covalent bond with Lys43. However, recent studies have suggested that a variety of small molecules can interact with the MR1-binding site. In this study, we have used several approaches to improve the binding pose prediction of covalent ligands to MR1, including docking in mutated receptors, and imposing simple pharmacophore constraints and structural water molecules. The careful assignment of pharmacophore constraints and inclusion of structural water molecules in the challenging docking process of covalent docking improved the binding pose prediction and virtual screening performance. In a retrospective virtual screening, the proposed approach exhibited EF1% and EF2% values of 7.4 and 5.5, respectively. Conversely, when using the mutated receptor, both EF1% and EF2% were recorded as 0 for the conventional docking method. The performance of the pharmacophore constraints was also evaluated on other covalent docking cases, and compared to previously reported results for common covalent docking methods. The proposed approach achieved an average RMSD of 2.55, while AutoDock4, CovDock, FITTED, GOLD, ICM-Pro, and MOE exhibited average RMSD values of 3.0, 2.93, 3.04, 4.93, 2.44, and 3.36, respectively. Our results demonstrate that the inclusion of simple pharmacophore constraints and structural waters can improve the prediction of binding poses of covalent ligands to MR1, which can aid in the discovery of novel immunotherapeutic agents.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03694-w.

At the origin of the selectivity of the chlorophyll-binding sites in light harvesting complex II (LHCII)

The photosynthetic light-harvesting complexes (LHCs) are responsible for light absorption due to their pigment-binding properties. These pigments are primarily Chlorophyll (Chl) molecules of type a and b, which ensure an excellent coverage of the visible light spectrum. To date, it is unclear which factors drive the selective binding of different Chl types in the LHC binding pockets. To gain insights into this, we employed molecular dynamics simulations on LHCII binding different Chl types. From the resulting trajectories, we have calculated the binding affinities per each Chl-binding pocket using the Molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) model. To further examine the importance of the nature of the axial ligand in tuning the Chl selectivity of the binding sites, we used Density Functional Theory (DFT) calculations. The results indicate that some binding pockets have a clear Chl selectivity, and the factors governing these selectivities are identified. Other binding pockets are promiscuous, which is consistent with previous in vitro reconstitution studies. DFT calculations show that the nature of the axial ligand is not a major factor in determining the Chl binding pocket selectivity, which is instead probably controlled by the folding process.

Structure restoration and aggregate inhibition of V30M mutant transthyretin protein by potential quinoline molecules

Transthyretin (TTR) is a tetrameric protein found in human plasma and cerebrospinal fluid that functions as a transporter of thyroxine (T4) and retinol. A mutation resulting in the substitution of valine to methionine at position 30 (V30M) is the most common mutation that destabilizes the tetramer structure of TTR protein resulting in a fatal neuropathy known as TTR amyloidosis. The V30M TTR-induced neuropathy can be inhibited through stabilization of the TTR tetramer by the binding of small molecules. We accessed the potential of in-house synthesized quinoline molecules to stabilize the V30M TTR structure and analyzed the impact of protein-ligand interactions through molecular docking, molecular dynamics (MD) simulations, steered MD, and umbrella sampling simulations. This study revealed that the binding of quinoline molecules reverted back the structural changes including the residual flexibility, changes in secondary structural elements, and also restored the alterations in the electrostatic surface potential induced by the V30M mutation. Further, the top-most 4G and 4R molecules were compared with an FDA-approved drug (Tafamidis) and a reference quinoline molecule 14C. Here, we intend to suggest that the quinoline molecules could revert the structural changes, cease tetramer dissociation, prevent abnormal oligomerization and therefore could be developed as an effective therapeutics against TTR amyloidosis.

Identification of natural products against enoyl-acyl-carrier-protein reductase in malaria via combined pharmacophore modeling, molecular docking and simulations studies

Plasmodium falciparum is counted as one of the deadly species causing malaria. In that respect, enoyl acyl carrier protein reductase is recognized as one of the attractive druggable targets for the identification of antimalarials. Thus, from the structural proteome of ENR, common feature pharmacophores were constructed. To identify the representative models, all the hypotheses were subjected to validation methods, like, test set, enrichment factor, and Güner-Henry method, and the selected representative hypotheses were used to screen out the drug-like natural products. Further, the screened candidates were advanced to molecular docking calculations. Based on the docking score criteria and presence of essential interaction with Tyr277, seven candidates were shortlisted to conduct the HYDE and QSAR assessment. Further, the stability of these complexes was evaluated by employing molecular dynamics simulations, molecular mechanics-generalized born surface area approach-based free binding energy calculations with the residue-wise contribution of PfENR to the total binding free energy of the complex. On comparing the root mean square deviation, and fluctuation plots of the docked candidates with the reference, all the candidates displayed stable behavior, and the same outcome was depicted from the secondary structure element. However, from the free energy calculations, and residue-wise contribution conducted after dynamics, it was observed that out of seven, only five candidates sustain the binding with Tyr277 and cofactor of PfENR. Therefore, in the current work, the hybrid study of screening and stability lead to the identification of five structurally diverse candidates that can be employed for the design of novel antimalarials.Communicated by Ramaswamy H. Sarma.

Design, synthesis, and molecular dynamics simulation studies of quinoline derivatives as protease inhibitors against SARS-CoV-2

A new series of quinoline derivatives has been designed and synthesized as probable protease inhibitors (PIs) against severe acute respiratory syndrome coronavirus 2. In silico studies using DS v20.1.0.19295 software have shown that these compounds behaved as PIs while interacting at the allosteric site of target Mpro enzyme (6LU7). The designed compounds have shown promising docking results, which revealed that all compounds formed hydrogen bonds with His41, His164, Glu166, Tyr54, Asp187, and showed π-interaction with His41, the highly conserved amino acids in the target protein. Toxicity Prediction by Komputer Assisted Technology results confirmed that the compounds were found to be less toxic than the reference drug. Further, molecular dynamics simulations were performed on compound 5 and remdesivir with protease enzyme. Analysis of conformational stability, residue flexibility, compactness, hydrogen bonding, solvent accessible surface area (SASA), and binding free energy revealed comparable stability of protease:5 complex to the protease: remdesivir complex. The result of hydrogen bonding showed a large number of intermolecular hydrogen bonds formed between protein residues (Glu166 and Gln189) and ligand 5, indicating strong interaction, which validated the docking result. Further, compactness analysis, SASA and interactions like hydrogen-bonding demonstrated inhibitory properties of compound 5 similar to the existing reference drug. Thus, the designed compound 5 might act as a potential inhibitor against the protease enzyme.Communicated by Ramaswamy H. SarmaHighlightsQuinoline derivatives have been designed as protease inhibitors against SARS-CoV-2.The compounds were docked at the allosteric site of SARS-CoV-2-Mpro enzyme (PDB ID: 6LU7) to study the stability of protein-ligand complex.Docking studies indicated the stable ligand-protein complexes for all designed compounds.The Toxicity Prediction by Komputer Assisted Technology protocol in DS v20.1.0.19295 software was used to evaluate the toxicity of the designed quinoline derivatives.Molecular dynamics studies indicated the formation of stable ligand-Mpro complexes.

Naturally Occurring 8ß,13ß-kaur-15-en-17-al and Anti-Malarial Activity from Podocarpus polystachyus Leaves

Despite much interest and studies toward the genus Podocarpus, the anti-malarial evaluation of Podocarpus polystachyus’s phytoconstituents remains lacking. Herein, the phytoconstituents of P. polystachyus leaves and their anti-malarial effect against Plasmodium falciparum were investigated for the first time. One new natural product, 8ß,13ß-kaur-15-en-17-al (1), along with three known compounds, 8ß,13ß-kaur-15-en-17-ol (2) and 13ß-kaur-16-ene (3), and α-tocopherol hydroquinone (4) were isolated via HR-ESI-MS and NMR analyses. Compounds 1 and 2 inhibited P. falciparum growth at 12 and 52 µM of IC₅₀, respectively. Their anti-malarial activity was associated with the in silico P. falciparum lactate dehydrogenase (PfLDH) inhibition. Molecular docking of ligands 1 and 2 with the putative target PfLDH revealed ~−2 kcal/mol of binding energies more negative than the control. Molecular dynamic simulations (100 ns) showed equal or smaller deviation values (RMSD, RMSF, Rg) and stronger interactions of PfLDH-1 and PfLDH-2 complexes via at least one consistent H-bond than the control. Additionally, a slightly increased PfLDH H-bond profile in their interactions improved the PfLDH dynamic and structural stabilities. Overall, this study supports the relevance of 1 and 2 as plasmodial growth inhibitors with their putative anti-PfLDH activity, which could be a potential scaffold for developing anti-malarial drugs.

Structure-based study to identify alkaloids as promising cytochrome P450 (CYP1A1) inhibitors: An in silico approach using virtual screening, molecular dynamic simulations, and binding free energy calculation

Carcinogens present in smokeless tobacco (SLT) like tobacco-specific nitrosamines can be metabolized by the cytochrome P450 (CYP450) enzyme. Functionally, the CYP450 enzyme resides in a heme pigment to perform the catalytic activity. The CYP1A1 is one of the main extrahepatic CYP450 enzymes known to detoxify toxic substances and activate carcinogens. The CYP1A1 inhibition by potential inhibitors reduce the chance of oral cancer. The current study aimed to explore more about the inhibitor binding site and identification of lead alkaloids, that could work as putative inhibitors against target CYP1A1. In respect, we have performed docking studies, virtual screening of alkaloids, and natural product libraries against CYP1A1 followed by molecular dynamic simulations and binding free energy calculations. Docking studies of tobacco-specific nitrosamine (TSNA) products and their similar carcinogen analogs revealed that the heme group is bound to the floor of the bowl-shaped cavity whereas carcinogens are bound to the roof of the rounded shape cavity. Furthermore, virtual screening and binding free energy calculations revealed Tomatidine as a putative inhibitor against CYP1A1. On the basis of altogether outcomes of the current study, we have concluded that the addition of lead-hit alkaloid Tomatidine and others in SLT products may be working as a supplement that could be able to reduce the expression of human CYP1A1 and suppresses carcinogenic by-products formations.

Integrating Network Pharmacology and Experimental Validation to Explore the Key Mechanism of Gubitong Recipe in the Treatment of Osteoarthritis

Background

Gubitong Recipe (GBT) is a prescription based on the Traditional Chinese Medicine (TCM) theory of tonifying the kidney yang and strengthening the bone. A previous multicentral randomized clinical trial has shown that GBT can effectively relieve joint pain and improve quality of life with a high safety in treating osteoarthritis (OA). This study is aimed at elucidating the active compounds, potential targets, and mechanisms of GBT for treating OA.

Method

The network pharmacology method was used to predict the key active compounds, targets, and mechanisms of GBT in treating OA. An OA rat model was established with Hulth surgery, and the pathological changes of articular cartilage were observed to evaluate the effects of GBT. Chondrocytes were stimulated with LPS to establish in vitro models, and key targets and mechanisms predicted by network pharmacology were verified via qRT-PCR, ELISA, western blot, and immunofluorescence. The Contribution Index Model and molecular docking were used to determine the key active compounds of GBT and the major nodes affecting predicted pathways.

Result

A total of 500 compounds were acquired from related databases, where 87 active compounds and their 254 corresponding targets were identified. 2979 OA-related genes were collected from three databases, 150 of which were GBT-regulating OA genes. The compound-target network weight analysis and PPI results showed that IL-6 and PGE2 are key targets of GBT in treating OA. KEGG results showed that PI3K/AKT, Toll-like receptor, NFκB, TNF, and HIF-1 are the key signaling pathways. An in vivo experiment showed that GBT could effectively suppress cartilage degradation of OA rats. In vitro experiments demonstrated that GBT can inhibit the key targets of KEGG-related pathways. Molecular-docking results suggested that luteolin, licochalcone A, and β-carotene were key targets of GBT, and the mechanisms may be associated with the NFκB signaling pathway. Blockage experiments showed that the NFκB pathway is the key pathway of GBT in treating OA.

Conclusion

This study verified that GBT can effectively protect articular cartilage through multitarget and multipathway, and its inhibitory effect on the NFκB pathway is the most key mechanism in treating OA.

Molecular Dynamics Simulation and Essential Dynamics of Deleterious Proline 12 Alanine Single-Nucleotide Polymorphism in PPARγ2 Associated with Type 2 Diabetes, Cardiovascular Disease, and Nonalcoholic Fatty Liver Disease

Background. Peroxisome proliferator-activated receptor-γ (PPARγ) gene is located at 3p25 position. PPARγ functions as the master regulator of glucose homeostasis and lipoprotein metabolism, and recent studies have reported that it is involved in various metabolic diseases such as diabetes mellitus, hyperlipidemia, coronary artery disease (CAD), and nonalcoholic fatty liver disease (NAFLD). PPARγ1 and PPARγ2 are necessary for the development of adipose tissue and insulin sensitivity regulation. But PPARγ2 is the isoform that was controlled in response to nutrient intake and obesity. Methodology. In this study, we used computational techniques to show the impact of Pro12Ala polymorphism on PPARγ2. The 3-D structure of PPARγ2 was modeled using I-TASSER server. The modeled structure was validated with the ZLab server, and the mutation was created with SPDB viewer. Stability prediction tools were used. Molecular dynamics simulation (MDS) was used to understand the structural and functional behavior of the wild type and mutant. Essential dynamics was also applied. Results and Conclusions. Stability prediction tools were showed that this mutation has a destabilizing effect on the PPARγ2 structure. The RMSD, RMSF, Rg, SASA, and DSSP were in line with H-bond results that showed less flexibility in the mutant structure. Essential dynamics was used to verify MDS results. Pro12Ala polymorphism leads to rigidity of the PPARγ2 protein and might disturb the conformational changes and interactions of PPARγ2 and results in type 2 diabetes mellitus (T2DM), CAD, and NAFLD. This study can help scientists to develop a drug therapy against these diseases.

Interaction analyses of hTAAR1 and mTAAR1 with antagonist EPPTB

Trace amine-associated receptor 1 (TAAR1) plays a critical role in regulating monoaminergic activity. EPPTB is the only known selective potent antagonist of the mouse (m) TAAR1 presently, while it was shown to be weak at antagonizing human (h) TAAR1. The lack of high-resolution structure of TAAR1 hinders the understanding of the differences in the interaction modes between EPPTB and m/hTARR1. The purpose of this study is to probe these interaction modes using homology modeling, molecular docking, molecular dynamics (MD) simulations, and molecular mechanics-generalized Born surface area (MM-GBSA) binding energy calculations. Eight populated conformers of hTAAR1-EPPTB complex were observed during the MD simulations and could be used in structure-based virtual screening in future. The MM-GBSA binding energy of hTAAR1-EPPTB complex (-96.5 kcal/mol) is larger than that of mTAAR1-EPPTB complex (-106.7 kcal/mol), which is consistent with the experimental finding that EPPTB has weaker binding affinity to hTAAR1. The several residues in binding site of hTAAR1 (F154^4.56, T194^5.42 and I290^7.39) are different from these of mTAAR1 (Y153^4.56, A193^5.42 and Y287^7.39), which might contribute to the binding affinity difference. Our docking analysis on another hTAAR1 antagonist Compound 3 has found that 1). this compound binds in different pockets of our mTAAR1 and hTAAR1 homology models with a slightly stronger binding affinity to hTAAR1; 2). both antagonists bind to a very similar pocket of hTAAR1.

Designing Novel Compounds for the Treatment and Management of RET-Positive Non-Small Cell Lung Cancer-Fragment Based Drug Design Strategy

Rearranged during transfection (RET) is an oncogenic driver receptor that is overexpressed in several cancer types, including non-small cell lung cancer. To date, only multiple kinase inhibitors are widely used to treat RET-positive cancer patients. These inhibitors exhibit high toxicity, less efficacy, and specificity against RET. The development of drug-resistant mutations in RET protein further deteriorates this situation. Hence, in the present study, we aimed to design novel drug-like compounds using a fragment-based drug designing strategy to overcome these issues. About 18 known inhibitors from diverse chemical classes were fragmented and bred to form novel compounds against RET proteins. The inhibitory activity of the resultant 115 hybrid molecules was evaluated using molecular docking and RF-Score analysis. The binding free energy and chemical reactivity of the compounds were computed using MM-GBSA and density functional theory analysis, respectively. The results from our study revealed that the developed hybrid molecules except for LF21 and LF27 showed higher reactivity and stability than Pralsetinib. Ultimately, the process resulted in three hybrid molecules namely LF1, LF2, and LF88 having potent inhibitory activity against RET proteins. The scrutinized molecules were then subjected to molecular dynamics simulation for 200 ns and MM-PBSA analysis to eliminate a false positive design. The results from our analysis hypothesized that the designed compounds exhibited significant inhibitory activity against multiple RET variants. Thus, these could be considered as potential leads for further experimental studies.

The molecular docking and molecular dynamics study of flavonol synthase and flavonoid 3'-monooxygenase enzymes involved for the enrichment of kaempferol

Kaempferol is a natural flavonol that shows many pharmacological properties including anti-inflammatory, antioxidant, anticancer, antidiabetic activities etc. It has been reported in many vegetables, fruits, herbs and medicinal plants. The enzyme flavonol synthase (FLS, EC 1.14.20.6) catalyses the conversion of dihydroflavonols to flavonols. Whereas flavonoid 3'-monooxygenase (F3'H, EC 1.14.14.82) catalyses the hydroxylation of dihydroflavonol, and flavonol. FLS is involved in the synthesis of the kaempferol whereas F3'H causes degradation of kaempferol. The present study aimed to analyse the binding affinity, stability and activating activity of enzyme FLS as well as inhibitory activity of enzyme F3'H involved in the enrichment of the kaempferol using the in-silico approaches. Computational study for physico-chemical properties, conserved domain identification, 3-D structure prediction and its validation, conservation analysis, molecular docking followed by molecular dynamics analysis of FLS and F3'H, protein-activator (FLS-LIG Complex) and protein-inhibitor (F3'H-LIG Complex) complexes have been performed. Other structural analyses like root mean square fluctuation (RMSF), root mean square deviation (RMSD), surface area solvent accessibility (SASA), radius of gyration (Rg), hydrogen bond analysis, principal component analysis (PCA), Poisson-Boltzmann analysis (MM_PBSA) and the dynamic cross correlation map (DCCM) analysis to explore the structural, functional and thermodynamic stability of the proteins and the complexes were also studied. The molecular docking result showed that FLS binds strongly with the activator ascorbate (CID _54670067) while F3'H binds with the inhibitor ketoconazole (CID_456201). The most powerful inhibitor (ketoconazole for F3'H) and activator (ascorbate for FLS) is determined by computing the thermodynamic binding free energy through MM_PBSA analysis. The current work provides wide-ranging structural and functional information about FLS and F3'H enzymes showing detailed molecular mechanism of kaempferol biosynthesis and its degradation and hence kaempferol enrichment. Finding of the present work opens up new possibilities for future research towards enrichment of kaempferol by using activator (ascorbate) for FLS and inhibitor (ketoconazole) for F3'H as well as for its large-scale production using in vitro approaches.Communicated by Ramaswamy H. Sarma.

Hexaconazole exposure ravages biosynthesis pathway of steroid hormones: revealed by molecular dynamics and interaction

Widespread application of hexaconazole for agriculture purpose poses a threat to human health by disrupting normal endocrine homeostasis. To avoid adverse health effects on human, it is crucial to identify the effects of hexaconazole on key enzymes responsible for steroidal hormone synthesis. In view of this, present study was conducted to investigate the interaction mechanisms of hexaconazole with key enzymes in comparison with their food drug administration (FDA) approved inhibitor by molecular docking and molecular dynamics simulations. Results indicate that hexaconazole contacts with the active site of the key enzymes required for steroidal hormonal synthesis. Results pertaining to root-mean-square deviation, root-mean-square calculation, radius of gyration, hydrogen bonding and solvent accessible surface area exhibited that the interaction pattern and stability of interaction of hexaconazole was similar to enzyme specific inhibitor. In addition, ligand and enzyme complex interaction energy of hexaconazole was almost similar to key enzyme and FDA-approved enzyme specific inhibitor complex. This study offers a molecular level of understanding of hexaconazole with different enzymes required for steroidal hormonal synthesis. Findings of the study clearly suggest that hexaconazole has efficacy to stably interact with various enzyme required to progress the pathway of hormonal synthesis. If incessant exposure of hexaconazole occurs during agricultural work it may lead to ravage hormonal synthesis or potent endocrine disruption. The result of binding energy and complex interaction energy is depicted in the graphical abstract.

Influence of extraction technology on rapeseed oil functional quality: a study on rapeseed polyphenols

Extraction technology can influence the vegetable oil functional quality. Polyphenols in rapeseed oil have been proved to be beneficial for cardiovascular health. In this study, we evaluated the effect of extraction methods on the functional quality of rapeseed oil from the perspective of phenolic compounds. The results showed that hot pressing produces the highest amount of phenolic compounds in rapeseed oil. Its most abundant phenolic compound, sinapine (9.18 μg g^-1), showed the highest activity in inhibiting anaerobic choline metabolism with an EC₅₀ value of 1.9 mM, whose downstream products are related to cardiovascular diseases. Molecular docking and molecular dynamics (MD) simulations revealed that sinapine exhibits good binding affinity toward CutC, and CutC-sinapine is a stable complex with fewer conformational fluctuations and similar tightness. Taken together, hot pressing can be considered the best extraction method for rapeseed oil from the perspective of phenolic compounds.

Antiplasmodial activity and cytotoxicity of plant extracts from the Asteraceae and Rubiaceae families

The increasing resistance of parasites to antimalarial drugs and the limited number of effective drugs are the greatest challenges in the treatment of malaria. It is necessary to search for an alternative medicine for use as a new, more effective antimalarial drug. Therefore, this study aimed to evaluate the in vitro antimalarial activity and cytotoxicity of extracts from plants belonging to the Asteraceae and Rubiaceae families. The phytoconstituents of one hundred ten ethanolic and aqueous extracts from different parts of twenty-three plant species were analyzed. Evaluation of their antimalarial activities against the chloroquine (CQ)-resistant Plasmodium falciparum (K1) strain was carried out using the lactate dehydrogenase (pLDH) assay, and their cytotoxicity in Vero cells was assessed using the 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyl tetrazolium bromide (MTT) colorimetric method. A total of 40.91% of the extracts were active antimalarial agents. Three extracts (2.73%) exhibited high antiplasmodial activity (IC50 < 10 μg/ml), twenty-four extracts (21.82%) were moderately active with IC50 values ranging from 10-50 μg/ml, and eighteen extracts (16.36%) were mildly active with IC50 values ranging from 50-100 μg/ml. The ethanolic leaf extract of Mussaenda erythrophylla (Dona Trining; Rubiaceae) exhibited the highest activity against P. falciparum, with an IC50 value of 3.73 μg/ml and a selectivity index (SI) of 30.74, followed by the ethanolic leaf extract of Mussaenda philippica Dona Luz x M. flava (Dona Marmalade; Rubiaceae) and the ethanolic leaf extract of Blumea balsamifera (Camphor Tree; Asteraceae), with IC50 values of 5.94 and 9.66 μg/ml and SI values of 25.36 and >20.70, respectively. GC-MS analysis of these three plant species revealed the presence of various compounds, such as squalene, oleic acid amide, β-sitosterol, quinic acid, phytol, oleamide, α-amyrin, sakuranin, quercetin and pillion. In conclusion, the ethanolic leaf extract of M. erythrophylla, the leaf extract of M. philippica Dona Luz x M. flava and the leaf extract of B. balsamifera had strong antimalarial properties with minimal toxicity, indicating that compounds from these plant species have the potential to be developed into new antiplasmodial agents.

Discovery of a Potent Candidate for RET-Specific Non-Small-Cell Lung Cancer-A Combined In Silico and In Vitro Strategy

Rearranged during transfection (RET) is a tyrosine kinase oncogenic receptor, activated in several cancers including non-small-cell lung cancer (NSCLC). Multiple kinase inhibitors vandetanib and cabozantinib are commonly used in the treatment of RET-positive NSCLC. However, specificity, toxicity, and reduced efficacy limit the usage of multiple kinase inhibitors in targeting RET protein. Thus, in the present investigation, we aimed to figure out novel and potent candidates for the inhibition of RET protein using combined in silico and in vitro strategies. In the present study, screening of 11,808 compounds from the DrugBank repository was accomplished by different hypotheses such as pharmacophore, e-pharmacophore, and receptor cavity-based models in the initial stage. The results from the different hypotheses were then integrated to eliminate the false positive prediction. The inhibitory activities of the screened compounds were tested by the glide docking algorithm. Moreover, RF score, Tanimoto coefficient, prime-MM/GBSA, and density functional theory calculations were utilized to re-score the binding free energy of the docked complexes with high precision. This procedure resulted in three lead molecules, namely DB07194, DB03496, and DB11982, against the RET protein. The screened lead molecules together with reference compounds were then subjected to a long molecular dynamics simulation with a 200 ns time duration to validate the inhibitory activity. Further analysis of compounds using MM-PBSA and mutation studies resulted in the identification of potent compound DB07194. In essence, a cell viability assay with RET-specific lung cancer cell line LC-2/ad was also carried out to confirm the in vitro biological activity of the resultant compound, DB07194. Indeed, the results from our study conclude that DB07194 can be effectively translated for this new therapeutic purpose, in contrast to the properties for which it was originally designed and synthesized.

A fragment-based approach identifies an allosteric pocket that impacts malate dehydrogenase activity

Malate dehydrogenases (MDHs) sustain tumor growth and carbon metabolism by pathogens including Plasmodium falciparum. However, clinical success of MDH inhibitors is absent, as current small molecule approaches targeting the active site are unselective. The presence of an allosteric binding site at oligomeric interface allows the development of more specific inhibitors. To this end we performed a differential NMR-based screening of 1500 fragments to identify fragments that bind at the oligomeric interface. Subsequent biophysical and biochemical experiments of an identified fragment indicate an allosteric mechanism of 4-(3,4-difluorophenyl) thiazol-2-amine (4DT) inhibition by impacting the formation of the active site loop, located >30 Å from the 4DT binding site. Further characterization of the more tractable homolog 4-phenylthiazol-2-amine (4PA) and 16 other derivatives are also reported. These data pave the way for downstream development of more selective molecules by utilizing the oligomeric interfaces showing higher species sequence divergence than the MDH active site.

Romero et al. perform NMR-based screening of 1500 fragments to identify fragments that bind at the oligomeric interface of malate dehydrogenase (MDH). Their study indicates an allosteric mechanism impacting enzymatic activity, paving the way for development of more selective molecules and a starting point for the future development of specific MDH inhibitors.

Antimalarial properties and molecular docking analysis of compounds from Dioscorea bulbifera L. as new antimalarial agent candidates

Background

At present, the emergence and spread of antimalarial drug resistance has become a significant problem worldwide. There has been a challenge in searching for natural products for the development of novel antimalarial drugs. Therefore, this study aims to evaluate compounds from Dioscorea bulbifera responsible for antimalarial properties and investigate potential interactions of the compounds with Plasmodium falciparum lactate dehydrogenase (PfLDH), an essential glycolytic enzyme in the parasite’s life cycle.

Methods

An in vitro study of antimalarial activity against chloroquine (CQ)-resistant Plasmodium falciparum (K1 strain) and CQ-sensitive P. falciparum (3D7 strain) was performed using the ³H-hypoxanthine uptake inhibition method. The cytotoxic effects of the pure compounds were tested against Vero cells using a 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) assay. The interactions of the compounds with the PfLDH active site were additionally investigated using a molecular docking method.

Results

Quercetin (6) exhibited the highest antimalarial activity against the P. falciparum K1 and 3D7 strains, with IC₅₀ values of 28.47 and 50.99 μM, respectively. 2,4,3′,5′-Tetrahydroxybibenzyl (9), 3,5-dimethoxyquercetin (4) and quercetin-3-O-β-D-galactopyranoside (14) also possessed antimalarial effects against these two strains of P. falciparum. Most pure compounds were nontoxic against Vero cells at a concentration of 80 μg/ml, except for compound 9, which had a cytotoxic effect with a CC₅₀ value of 16.71 μM. The molecular docking results indicated that 9 exhibited the best binding affinity to the PfLDH enzyme in terms of low binding energy (− 8.91 kcal/mol) and formed strong hydrogen bond interactions with GLY29, GLY32, THR97, GLY99, PHE100, THR101 and ASN140, amino acids as active sites. In addition, 6 also possessed remarkable binding affinity (− 8.53 kcal/mol) to PfLDH by interacting with GLY29, ILE31, ASP53, ILE54, THR97 and THR101.

Conclusion

Quercetin is a major active compound responsible for the antimalarial activity of D. bulbifera and is an inhibitor of PfLDH. These findings provide more evidence to support the traditional use of D. bulbifera for malaria treatment. Structural models of its interactions at the PfLDH active site are plausibly useful for the future design of antimalarial agents.

Repurposing of approved drug molecules for viral infectious diseases: a molecular modelling approach

The identification of new viral drugs has become a task of paramount significance due to the frequent occurrence of viral infections and especially during the current pandemic. Despite the recent advancements, the development of antiviral drugs has not made parallel progress. Reduction of time frame and cost of the drug development process is the major advantage of drug repurposing. Therefore, in this study, a drug repurposing strategy using molecular modelling techniques, i.e. biological activity prediction, virtual screening, and molecular dynamics simulation was employed to find promising repurposing candidates for viral infectious diseases. The biological activities of non-redundant (4171) drug molecules were predicted using PASS analysis, and 1401 drug molecules were selected which showed antiviral activities in the analysis. These drug molecules were subjected to virtual screening against the selected non-structural viral proteins. A series of filters, i.e. top 10 drug molecules based on binding affinity, mean value of binding affinity, visual inspection of protein-drug complexes, and number of H-bond between protein and drug molecules were used to narrow down the drug molecules. Molecular dynamics simulation analysis was carried out to validate the intrinsic atomic interactions and binding conformations of protein-drug complexes. The binding free energies of drug molecules were assessed by employing MMPBSA analysis. Finally, nine drug molecules were prioritized, as promising repurposing candidates with the potential to inhibit the selected non-structural viral proteins.Communicated by Ramaswamy H. Sarma.