Docking, molecular dynamics, binding energy-MM-PBSA studies of naphthofuran derivatives to identify potential dual inhibitors against BACE-1 and GSK-3β

Identifying potential Alzheimer's disease therapeutics through GSK-3β inhibition: A molecular docking and dynamics approach

Emerging as a promising drug target for Alzheimer's disease (AD) therapy, glycogen synthase kinase 3β (GSK-3β) has garnered attention. This study sought to rigorously scrutinize a compendium of natural compounds retrieved from the ZINC database through pharmacodynamic experiments, employing a 1 H-indazole-3-carboxamide (INDZ) scaffold, to identify compounds capable of inhibiting the GSK-3β protein. Utilizing a multi-step approach, the study involved pharmacophore analysis, followed by molecular docking to select five promising ligands for further investigation. Subsequently, ESMACS simulations were employed to assess the stability of the ligand-protein interactions. Evaluation of the binding modes and free energy of the ligands revealed that five compounds (2a-6a) exhibited crucial interactions with the active site residues. Furthermore, various methodologies, including hydrogen bond and clustering analyses, were utilized to ascertain their inhibitory potential and elucidate the factors contributing to ligand binding in the protein's active site. The findings from MMPBSA/GBSA analysis indicated that these five selected small molecules closely approached the IC50 value of the reference ligand (OH8), yielding energy values of -34.85, -32.58, -31.71, and -30.39 kcal/mol, respectively. Additionally, an assessment of the interactions using hydrogen bond and dynamic analyses delineated the effective binding of the ligands with the binding pockets in the protein. Through computational analysis, we obtained valuable insights into the molecular mechanisms of GSK-3β, aiding in the development of more potent inhibitors.

Therapeutic Potential Effect of Glycogen Synthase Kinase 3 Beta (GSK-3β) Inhibitors in Parkinson Disease: Exploring an Overlooked Avenue

Parkinson's disease (PD) is a progressive neurodegenerative disease of the brain due to degeneration of dopaminergic neurons in the substantia nigra (SN). Glycogen synthase kinase 3 beta (GSK-3β) is implicated in the pathogenesis of PD. Therefore, the purpose of the present review was to revise the mechanistic role of GSK-3β in PD neuropathology, and how GSK-3β inhibitors affect PD neuropathology. GSK-3 is a conserved threonine/serine kinase protein that is intricate in the regulation of cellular anabolic and catabolic pathways by modulating glycogen synthase. Over-expression of GSK-3β is also interconnected with the development of different neurodegenerative diseases. However, the underlying mechanism of GSK-3β in PD neuropathology is not fully clarified. Over-expression of GSK-3β induces the development of PD by triggering mitochondrial dysfunction and oxidative stress in the dopaminergic neurons of the SN. NF-κB and NLRP3 inflammasome are activated in response to dysregulated GSK-3β in PD leading to progressive neuronal injury. Higher expression of GSK-3β in the early stages of PD neuropathology might contribute to the reduction of neuroprotective brain-derived neurotrophic factor (BDNF). Thus, GSK-3β inhibitors may be effective in PD by reducing inflammatory and oxidative stress disorders which are associated with degeneration of dopaminergic in the SN.

Targeting glycogen synthase kinase-3β for Alzheimer's disease: Recent advances and future Prospects

Senile plaques induced by β-amyloid (Aβ) abnormal aggregation and neurofibrillary tangles (NFT) caused by tau hyperphosphorylation are important pathological manifestations of Alzheimer's disease (AD). Glycogen synthase kinase-3 (GSK-3) is a conserved kinase; one member GSK-3β is highly expressed in the AD brain and involved in the formation of NFT. Hence, pharmacologically inhibiting GSK-3β activity and expression is a good approach to treat AD. As summarized in this article, multiple GSK-3β inhibitors has been comprehensively summarized over recent five years. However, only lithium carbonate and Tideglusib have been studied in clinical trials of AD. Besides ATP-competitive and non-ATP-competitive inhibitors, peptide inhibitors, allosteric inhibitors and other types of inhibitors have gradually attracted more interest. Moreover, considering the close relationship between GSK-3β and other targets involved in cholinergic hypothesis, Aβ aggregation hypothesis, tau hyperphosphorylation hypothesis, oxidative stress hypothesis, neuro-inflammation hypothesis, etc., diverse multifunctional molecules and multi-target directed ligands (MTDLs) have also been disclosed. We hope that these recent advances and critical perspectives will facilitate the discovery of safe and effective GSK-3β inhibitors for AD treatment.

In silico design and identification of new peptides for mitigating hIAPP aggregation in type 2 diabetes

The aberrant misfolding and self-aggregation of human islet amyloid polypeptide (hIAPP or amylin) into cytotoxic aggregates are implicated in the pathogenesis of type 2 diabetes (T2D). Among various inhibitors, short peptides derived from the amyloidogenic regions of hIAPP have been employed as hIAPP aggregation inhibitors due to their low immunogenicity, biocompatibility, and high chemical diversity. Recently, hIAPP fragment HSSNN18-22 was identified as an amyloidogenic sequence and displayed higher antiproliferative activity to RIN-5F cells. Various hIAPP aggregation inhibitors have been designed by chemical modifications of the highly amyloidogenic sequence (NFGAIL) of hIAPP. In this work, a library of pentapeptides based on fragment HSSNN18-22 was designed and assessed for their efficacy in blocking hIAPP aggregation using an integrated computational screening approach. The binding free energy calculations by molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method identified HSSQN and HSSNQ that bind to hIAPP monomer with a binding affinity of -21.25 ± 4.90 and -19.73 ± 3.10 kcal/mol, respectively, which is notably higher as compared to HSSNN (-11.90 ± 4.12 kcal/mol). The sampling of the non aggregation-prone helical conformation was notably increased from 23.5 ± 3.0 in the hIAPP monomer to 38.1 ± 3.6, and 33.8 ± 3.0% on the incorporation of HSSQN, and HSSNQ, respectively, which indicate reduced aggregation propensity of hIAPP monomer. The pentapeptides, HSSQN and HSSNQ, identified as hIAPP aggregation inhibitors in this work can be further conjugated with various metal chelating peptides to yield more efficacious and clinically relevant multifunctional modulators for targeting various pathological hallmarks of T2D.

Structural and molecular insights into tacrine-benzofuran hybrid induced inhibition of amyloid-β peptide aggregation and BACE1 activity

Amyloid-β (Aβ) aggregation and β-amyloid precursor protein cleaving enzyme 1 (BACE1) are the potential therapeutic drug targets for Alzheimer's disease (AD). A recent study highlighted that tacrine-benzofuran hybrid C1 displayed anti-aggregation activity against Aβ42 peptide and inhibit BACE1 activity. However, the inhibition mechanism of C1 against Aβ42 aggregation and BACE1 activity remains unclear. Thus, molecular dynamics (MD) simulations of Aβ42 monomer and BACE1 with and without C1 were performed to inspect the inhibitory mechanism of C1 against Aβ42 aggregation and BACE1 activity. In addition, a ligand-based virtual screening followed by MD simulations was employed to explore potent new small-molecule dual inhibitors of Aβ42 aggregation and BACE1 activity. MD simulations highlighted that C1 promotes the non aggregating helical conformation in Aβ42 and destabilizes D23-K28 salt bridge that plays a vital role in the self-aggregation of Aβ42. C1 displays a favourable binding free energy (-50.7 ± 7.3 kcal/mol) with Aβ42 monomer and preferentially binds to the central hydrophobic core (CHC) residues. MD simulations highlighted that C1 strongly interacted with the BACE1 active site (Asp32 and Asp228) and active pockets. The scrutiny of interatomic distances among key residues of BACE1 highlighted the close flap (non-active) position in BACE1 on the incorporation of C1. The MD simulations explain the observed high inhibitory activity of C1 against Aβ aggregation and BACE1 in the in vitro studies. The ligand-based virtual screening followed by MD simulations identified CHEMBL2019027 (C2) as a promising dual inhibitor of Aβ42 aggregation and BACE1 activity.Communicated by Ramaswamy H. Sarma.

Repurposing of phyto-ligand molecules from the honey bee products for Alzheimer's disease as novel inhibitors of BACE-1: small molecule bioinformatics strategies as amyloid-based therapy

Alzheimer's disease (AD) is one of the neurodegenerative diseases, manifesting dementia, spatial disorientation, language, cognitive, and functional impairment, mainly affects the elderly population with a growing concern about the financial burden on society. Repurposing can improve the traditional progress of drug design applications and could speed up the identification of innovative remedies for AD. The pursuit of potent anti-BACE-1 drugs for AD treatment has become a pot boiler topic in the recent past and to instigate the design of novel improved inhibitors from the bee products. Drug-likeness characteristics (ADMET: absorption, distribution, metabolism, excretion, and toxicity), docking (AutoDock Vina), simulation (GROMACS), and free energy interaction (MM-PBSA, molecular mechanics Poisson-Boltzmann surface area) analyses were performed to identify the lead candidates from the bee products (500 bioactives from the honey, royal jelly, propolis, bee bread, bee wax, and bee venom) for Alzheimer's disease as novel inhibitors of BACE-1 (beta-site amyloid precursor protein cleaving enzyme (1) receptor using appropriate bioinformatics tools. Forty-four bioactive lead compounds were screened from the bee products through high throughput virtual screening on the basis of their pharmacokinetic and pharmacodynamics characteristics, showing favorable intestinal and oral absorption, bioavailability, blood brain barrier penetration, less skin permeability, and no inhibition of cytochrome P450 inhibitors. The docking score of the forty-four ligand molecules was found to be between -4 and -10.3 kcal/mol, respectively, exhibiting strong binding affinity to BACE1 receptor. The highest binding affinity was observed in the rutin (-10.3 kcal/mol), 3,4-dicaffeoylquinic acid (-9.5 kcal/mol), nemorosone (-9.5 kcal/mol), and luteolin (-8.9 kcal/mol). Furthermore, these compounds demonstrated high total binding energy -73.20 to -105.85 kJ/mol), and low root mean square deviation (0.194-0.202 nm), root mean square fluctuation (0.0985-0.1136 nm), radius of gyration (2.12 nm), number of H-bonds (0.778-5.436), and eigenvector values (2.39-3.54 nm²) in the molecular dynamic simulation, signifying restricted motion of Cα atoms, proper folding and flexibility, and highly stable with compact of the BACE1 receptor with the ligands. Docking and simulation studies concluded that rutin, 3,4-dicaffeoylquinic acid, nemorosone, and luteolin are plausibly used as novel inhibitors of BACE1 to combat AD, but further in-depth experimental investigations are warranted to prove these in silico findings.

Effect Produced by a Cyclooctyne Derivative on Both Infarct Area and Left Ventricular Pressure via Calcium Channel Activation

BACKGROUND

There are reports which indicate that some cyclooctyne derivatives may exert changes in cardiovascular system; however, its molecular mechanism is not very clear.

OBJECTIVE

The aim of this study was to evaluate the biological activity of four cyclooctyne derivatives (compounds 1: to 4: ) produced on infarct area and left ventricular pressure.

METHODS

Biological activity produced by cyclooctyne derivatives on infarct area was determinate using an ischemia/reperfusion injury model. In addition, to characterize the molecular mechanism of this effect, the following strategies were carried out as follows; i) biological activity produced by cyclooctyne derivative (compound 4: ) on either perfusion pressure or left ventricular pressure was evaluated using an isolated rat heart; ii) theoretical interaction of cyclooctyne derivative with calcium channel (1t0j protein surface) using a docking model.

RESULTS

The results showed that cyclooctyne derivative (compound 4: ) decrease infarct area of in a dose-dependent manner compared with compound 1: to 3: . Besides, this cyclooctyne derivative increase both perfusion pressure and left ventricular pressure which was inhibited by nifedipine. Other theoretical data suggests that cyclooctyne derivative could interact with some aminoacid residues (Met₈₃, Ile₈₅, Ser₈₆, Leu₁₀₈, Glu₁₁₄) involved in 1t0j protein surface.

CONCLUSIONS

All these data indicate that cyclooctyne derivative increase left ventricular pressure via calcium channel activation and this phenomenon could be translated as a decrease of infarct area.

Designing Tailored Thiosemicarbazones with Bespoke Properties: The Styrene Moiety Imparts Potent Activity, Inhibits Heme Center Oxidation, and Results in a Novel "Stealth Zinc(II) Complex"

A novel, potent, and selective antitumor agent, namely (E)-3-phenyl-1-(2-pyridinyl)-2-propen-1-one 4,4-dimethyl-3-thiosemicarbazone (PPP44mT), and its analogues were synthesized and characterized and displayed strikingly distinctive properties. This activity was mediated by the inclusion of a styrene moiety, which through steric and electrochemical mechanisms prevented deleterious oxy-myoglobin or oxy-hemoglobin oxidation relative to other potent thiosemicarbazones, i.e., di-2-pyridylketone-4-cyclohexyl-4-methyl-3-thiosemicarbazone (DpC) or di-2-pyridylketone-4,4-dimethyl-3-thiosemicarbazone (Dp44mT). Structure-activity relationship analysis demonstrated specific tuning of PPP44mT electrochemistry further inhibited oxy-myoglobin or oxy-hemoglobin oxidation. Both PPP44mT and its Cu(II) complexes showed conspicuous almost immediate cytotoxicity against SK-N-MC tumor cells (within 3 h). In contrast, [Zn(PPP44mT)₂] demonstrated a pronounced delay in activity, taking 48 h before marked antiproliferative efficacy was apparent. As such, [Zn(PPP44mT)₂] was designated as a "stealth Zn(II) complex" that overcomes the near immediate cytotoxicity of PPP44mT or its copper complexes. Upon examination of the suppression of oncogenic signaling, [Zn(PPP44mT)₂] was superior at inhibiting cyclin D1 expression compared to DpC or Dp44mT.

The investigation into the effect of the length of RGD peptides and temperature on the interaction with the αIIbβ3 integrin: a molecular dynamic study

The tripeptide Arg-Gly-Asp acid (RGD) is a protein sequence in the binding of proteins to cell surfaces, and is involved in various biological processes such as cell adhesion to the extracellular matrix, platelet activation, hemostasis, etc. The C2 domain of the Von Willebrand Factor (VWF), containing the RGD motif, plays an important role in the initial homeostasis process. It binds to the αIIbβ3 integrin and stimulates platelet aggregation. We have investigated, using the molecular Dynamic (MD) simulation method, the effect of the RGD-peptide length, and temperature variation, on the binding to the αIIbβ3 integrin receptor. We examined 10 different structural modes of the αIIbβ3 at three different temperatures; 237 K, 310 K and 318 K. Our findings show that the amino acids that form a binding pocket include Asp₂₂₄, Tyr₂₃₄, Ser₂₂₆, Tyr₁₉₀, Tyr₁₈₉, Trp₂₆₀, Trp₂₆₂, Asp_259, Lys₂₅₃, Arg₂₁₄, Asp₂₁₇, Ser₁₆₁ and Ala₂₁₈ and that the ligand-receptor interaction was increased at higher temperatures. It was also found that the increase in the number of ligands' amino acids and their types (% glycine) plays an important role in the stability, conformation, and ligand-receptor interaction.Communicated by Ramaswamy H. Sarma.

Mass Spectrometric Behavior and Molecular Mechanisms of Fermented Deoxyanthocyanidins to Alleviate Ulcerative Colitis Based on Network Pharmacology

Aims. Ulcerative colitis (UC) is a type of chronic idiopathic inflammatory bowel disease with a multifactorial pathogenesis and limited treatment options. The aim of the present study is to investigate the hydrogen deuterium exchange mass spectrometry (HDX-MS) behaviors of fermented deoxyanthocyanidins and their molecular mechanisms to alleviate UC by using quantum chemistry and network pharmacology. Methods. Tandem MS indicated at least two fragmentation pathways through which deuterated vinylphenol-deoxyanthocyanidins could generate different product ions. Quantum calculations were conducted to determine the transition states of the relevant molecules and analyze their optimized configuration, vibrational characteristics, intrinsic reaction coordinates, and corresponding energies. The potential targets of deoxyanthocyanidins in UC were screened from a public database. The R package was used for Gene Ontology (GO) and KEGG pathway analyses, and the protein–protein interactions (PPIs) of the targets were assessed using Search Tool for the Retrieval of Interacting Genes (STRING). Finally, molecular docking was implemented to analyze the binding energies and action modes of the target compounds through the online tool CB-Dock. Results. Quantum calculations indicated two potential fragmentation pathways involving the six-membered ring and dihydrogen cooperative transfer reactions of the vinylphenol-deoxyanthocyanidins. A total of 146 and 57 intersecting targets of natural and fermented deoxyanthocyanidins were separately screened out from the UC database and significant overlaps in GO terms and KEGG pathways were noted. Three shared hub targets (i.e., PTGS2, ESR1, and EGFR) were selected from the two PPI networks by STRING. Molecular docking results showed that all deoxyanthocyanidins have a good binding potential with the hub target proteins and that fermented deoxyanthocyanidins have lower binding energies and more stable conformations compared with natural ones. Conclusions. Deoxyanthocyanidins may provide anti-inflammatory, antioxidative, and immune system regulatory effects to suppress UC progression. It is proposed for the first time that fermentation of deoxyanthocyanidins can help adjust the structure of the intestinal microbiota and increase the biological activity of the natural compounds against UC. Furthermore, HDX-MS is a helpful strategy to analyze deoxyanthocyanidin metabolites with unknown structures.

A Multi-target Drug Designing for BTK, MMP9, Proteasome and TAK1 for the Clinical Treatment of Mantle Cell Lymphoma

BACKGROUND

Mantle cell lymphoma (MCL) is a type of non-Hodgkin lymphoma characterized by the mutation and overexpression of the cyclin D1 protein by the reciprocal chromosomal translocation t(11;14)(q13:q32).

AIM

The present study aims to identify potential inhibition of MMP9, Proteasome, BTK, and TAK1 and determine the most suitable and effective protein target for the MCL.

METHODOLOGY

Nine known inhibitors for MMP9, 24 for proteasome, 15 for BTK and 14 for TAK1 were screened. SB-3CT (PubChem ID: 9883002), oprozomib (PubChem ID: 25067547), zanubrutinib (PubChem ID: 135565884) and TAK1 inhibitor (PubChem ID: 66760355) were recognized as drugs with high binding capacity with their respective protein receptors. 41, 72, 102 and 3 virtual screened compounds were obtained after the similarity search with compound (PubChem ID:102173753), PubChem compound SCHEMBL15569297 (PubChem ID:72374403), PubChem compound SCHEMBL17075298 (PubChem ID:136970120) and compound CID: 71814473 with best virtual screened compounds.

RESULT

MMP9 inhibitors show commendable affinity and good interaction profile of compound holding PubChem ID:102173753 over the most effective established inhibitor SB-3CT. The pharmacophore study of the best virtual screened compound reveals its high efficacy based on various interactions. The virtual screened compound's better affinity with the target MMP9 protein was deduced using toxicity and integration profile studies.

CONCLUSION

Based on the ADMET profile, the compound (PubChem ID: 102173753) could be a potent drug for MCL treatment. Similar to the established SB-3CT, the compound was non-toxic with LD50 values for both the compounds lying in the same range.

Glycogen Synthase Kinase 3 (GSK3): Its Role and Inhibitors

Glycogen Synthase Kinase 3 (GSK3) is one of the Serine/Threonine protein kinases, which has gained a lot of attention for its role in a variety of pathways. It has two isoforms, GSK3α and GSK3β. However, GSK3β is highly expressed in different areas of the brain and has been implicated in Alzheimer's disease as it is involved in tau phosphorylation. Due to its high specificity concerning substrate recognition, GSK3 has been considered as an important target. In the last decade, several GSK3 inhibitors have been reported and two molecules are in clinical trials. This review collates the information published in the last decade about the role of GSK3 in Alzheimer's disease and progress in the development of its inhibitors. Using this collated information, medicinal chemists can strategize and design novel GSK3 inhibitors that could be useful in the treatment of Alzheimer's disease.

Computational prediction method to decipher receptor-glycoligand interactions in plant immunity

Microbial and plant cell walls have been selected by the plant immune system as a source of microbe- and plant damage-associated molecular patterns (MAMPs/DAMPs) that are perceived by extracellular ectodomains (ECDs) of plant pattern recognition receptors (PRRs) triggering immune responses. From the vast number of ligands that PRRs can bind, those composed of carbohydrate moieties are poorly studied, and only a handful of PRR/glycan pairs have been determined. Here we present a computational screening method, based on the first step of molecular dynamics simulation, that is able to predict putative ECD-PRR/glycan interactions. This method has been developed and optimized with Arabidopsis LysM-PRR members CERK1 and LYK4, which are involved in the perception of fungal MAMPs, chitohexaose (1,4-β-d-(GlcNAc)6 ) and laminarihexaose (1,3-β-d-(Glc)6 ). Our in silico results predicted CERK1 interactions with 1,4-β-d-(GlcNAc)6 whilst discarding its direct binding by LYK4. In contrast, no direct interaction between CERK1/laminarihexaose was predicted by the model despite CERK1 being required for laminarihexaose immune activation, suggesting that CERK1 may act as a co-receptor for its recognition. These in silico results were validated by isothermal titration calorimetry binding assays between these MAMPs and recombinant ECDs-LysM-PRRs. The robustness of the developed computational screening method was further validated by predicting that CERK1 does not bind the DAMP 1,4-β-d-(Glc)6 (cellohexaose), and then probing that immune responses triggered by this DAMP were not impaired in the Arabidopsis cerk1 mutant. The computational predictive glycan/PRR binding method developed here might accelerate the discovery of protein-glycan interactions and provide information on immune responses activated by glycoligands.

Enfuvirtide, an HIV-1 fusion inhibitor peptide, can act as a potent SARS-CoV-2 fusion inhibitor: an in silico drug repurposing study

Regarding the urgency of therapeutic measures for coronavirus disease 2019 (COVID-19) pandemic, the use of available drugs with FDA approval is preferred because of the less time and cost required for their development. In silico drug repurposing is an accurate way to speed up the screening of the existing FDA-approved drugs to find a therapeutic option for COVID-19. The similarity in SARS-CoV-2 and HIV-1 fusion mechanism to host cells can be a key point for Inhibit SARS-CoV-2 entry into host cells by HIV fusion inhibitors. Accordingly, in this study, an HIV-1 fusion inhibitor called Enfuvirtide (Enf) was selected. The affinity and essential residues involving in the Enf binding to the S2 protein of SARS-CoV-2, HIV-1 gp41 protein and angiotensin-converting enzyme 2 (ACE-2) as a negative control, was evaluated using molecular docking. Eventually, Enf-S2 and Enf-gp41 protein complexes were simulated by molecular dynamics (MD) in terms of binding affinity and stability. Based on the most important criteria such as docking score, cluster size, energy and dissociation constant, the strongest interaction was observed between Enf with the S2 protein. In addition, MD results confirmed that Enf-S2 protein interaction was remarkably stable and caused the S2 protein residues to undergo the fewest fluctuations. In conclusion, it can be stated that Enf can act as a strong SARS-CoV-2 fusion inhibitor and demonstrates the potential to enter the clinical trial phase of COVID-19.

Communicated by Ramaswamy H. Sarma

Designing high affinity target-binding peptides to HLA-E: a key membrane antigen of multiple myeloma

Multiple myeloma (MM) is a plasma cell malignancy that is currently incurable. Finding new targets and designing drugs are crucial for the treatment of MM. The two datasets (GSE6691 and GSE39754) are used to screen highly expressed antigen on MM cells. HLA-E was an ideal target for it was a hub gene, and also located in one of the key clusters. Highly expression of HLA-E mRNA on MM cells was also confirmed by real-time qPCR testing the MM patients' samples in Shengjing hospital. Crystal structure of HLA-E was obtained from Protein Data Bank (PDB ID: 3CDG) which was used to design targeting peptides with Molecular Operating Environment software. By analyzing interaction between CD94/NKG2A and HLA-E, a peptide with twelve amino acids was screened as a model peptide. Peptides library was constructed by randomly replaced non-key amino acid. Peptide-protein docking method was used to identify high affinity peptides. PEPTIDE 1-3 and model peptide were synthesized and identified the affinity to HLA-E by flow cytometer and confocal laser microscopy. At last, PEPTIDE3 (NALDEYCEDKNR) was found with the highest affinity. Taking all, HLA-E is a new treatment target, and PEPTIDE 3 is an ideal high affinity target-binding peptide candidate.

Effects of unsaturated fatty acids (Arachidonic/Oleic Acids) on stability and structural properties of Calprotectin using molecular docking and molecular dynamics simulation approach

Calprotectin is a heterodimeric protein complex with two subunits called S100A8/A9. The protein has an essential role in inflammation process and various human diseases. It has the ability to bind to unsaturated fatty acids including Arachidonic acid, Oleic acid and etc., which could be considered as a major carrier for fatty acids. In this study we aimed to appraise the thermodynamics and structural changes of Calprotectin in presence of Arachidonic acid/Oleic acid) using docking and molecular dynami simulation method. To create the best conformation of Calprotectin-Oleic acid/Arachidonic acid complexes, the docking process was performed. The complexes with the best binding energy were selected as the models for molecular dynamics simulation process. Furthermore, the structural and thermodynamics properties of the complexes were evaluated too. The Root Mean Square Deviation and Root Mean Square Fluctuation results showed that the binding of Arachidonic acid/Oleic acid to Calprotectin can cause the protein structural changes which was confirmed by Define Secondary Structure of Proteins results. Accordingly, the binding free energy results verified that binding of Oleic acid to Calprotectin leads to instability of S100A8/A9 subunits in the protein. Moreover, the electrostatic energy contribution of the complexes (Calprotectin-Oleic acid/Arachidonic acid) was remarkably higher than van der Waals energy. Thus, the outcome of this study confirm that Oleic acid has a stronger interaction with Calprotectin in comparison with Arachidonic acid. Our findings indicated that binding of unsaturated fatty acids to Calprotectin leads to structural changes of the S100A8/A9 subunits which could be beneficial to play a biological role in inflammation process.

In silico assessment of human Calprotectin subunits (S100A8/A9) in presence of sodium and calcium ions using Molecular Dynamics simulation approach

Calprotectin is a heterodimeric protein complex which consists of two subunits including S100A8 and S100A9. This protein has a major role in different inflammatory disease and various types of cancers. In current study we aimed to evaluate the structural and thermodynamic changes of the subunits and the complex in presence of sodium and calcium ions using molecular dynamics (MD) simulation. Therefore, the residue interaction network (RIN) was visualized in Cytoscape program. In next step, to measure the binding free energy, the potential of mean force (PMF) method was performed. Finally, the molecular mechanics Poisson-Boltzmann surface area (MMPBSA) method was applied as an effective tool to calculate the molecular model affinities. The MD simulation results of the subunits represented their structural changes in presence of Ca2+. Moreover, the RIN and Hydrogen bond analysis demonstrated that cluster interactions between Calprotectin subunits in presence of Ca2+ were greater in comparison with Na+. Our findings indicated that the binding free energy of the subunits in presence of Ca2+ was significantly greater than Na+. The results revealed that Ca2+ has the ability to induce structural changes in subunits in comparison with Na+ which lead to create stronger interactions between. Hence, studying the physical characteristics of the human proteins could be considered as a powerful tool in theranostics and drug design purposes.

In Silico Studies on Compounds Derived from Calceolaria: Phenylethanoid Glycosides as Potential Multitarget Inhibitors for the Development of Pesticides

An increasing occurrence of resistance in insect pests and high mammal toxicity exhibited by common pesticides increase the need for new alternative molecules. Among these alternatives, bioinsecticides are considered to be environmentally friendly and safer than synthetic insecticides. Particularly, plant extracts have shown great potential in laboratory conditions. However, the lack of studies that confirm their mechanisms of action diminishes their potential applications on a large scale. Previously, we have reported the insect growth regulator and insecticidal activities of secondary metabolites isolated from plants of the Calceolaria genus. Herein, we report an in silico study of compounds isolated from Calceolaria against acetylcholinesterase, prophenoloxidase, and ecdysone receptor. The molecular docking results are consistent with the previously reported experimental results, which were obtained during the bioevaluation of Calceolaria extracts. Among the compounds, phenylethanoid glycosides, such as verbascoside, exhibited good theoretical affinity to all the analyzed targets. In light of these results, we developed an index to evaluate potential multitarget insecticides based on docking scores.

Identification of MMP9 as a novel key gene in mantle cell lymphoma based on bioinformatic analysis and design of cyclic peptides as MMP9 inhibitors based on molecular docking

Mantle cell lymphoma (MCL) is an aggressive disease. MCL is associated with poor patient prognosis and limited survival. To identify key genes and explore targeting cyclic peptide inhibitors for the treatment of MCL, we downloaded two gene expression profiles (GSE32018 and GSE9327) from the Gene Expression Omnibus (GEO) database. We screened 84 differentially expressed genes (DEGs). Pathway analysis showed that DEMs were mainly enriched in the ‘Pathway in cancer’, ‘PI3K-Akt signaling pathway’, ‘Cytokine-cytokine receptor interaction’, ‘Rap1 signaling pathway’, ‘NF-κB signaling pathway’ and ‘Leukocyte trans-endothelial migration’. We subsequently constructed a protein-protein interaction (PPI) network of DEGs. In addition, matrix metalloproteinase 9 (MMP9) with a high degree in the PPI network was identified as a hub gene in MCL. Meanwhile in the Molecular Complex Detection (MCODE) analysis, MMP9 was located in the important cluster. Thus, MMP9 can be used as a therapeutic target for MCL and we designed cyclic peptides as MMP9 inhibitors. MMP9 protein structure was gathered from the Protein Data Bank (PDB), with a PDB ID: 1L6J. MMP9 and cyclic peptides were docked using Molecular Operating Environment (MOE) software after structural optimization. It was revealed that cyclic peptide 2 bound deeply in the binding pocket of MMP9 and had interaction with the active-site Zn²⁺ ion in the catalytic domain. Cyclic peptides 1, 2, 4–6 also displayed potential interaction with active residues of MMP9; thus, these cyclic peptides can serve as potential drug candidates to block MMP9 activity and future studies are warranted to confirm their efficacy.

In silico identification and design of potent peptide inhibitors against PDZ-3 domain of Postsynaptic Density Protein (PSD-95)

Unique intrinsic properties of peptides like low toxicity, high biological activity, and specificity make them attractive therapeutic agents. PDZ-binding peptide inhibitors have been demonstrated for curing of Alzheimer, Parkinson, Dementia, and other central nervous system ailments. In this article, we report the successful use of an integrated computational protocol to analyze the structural basis of how peptides bind to the shallow groove of the third PDZ domain (PDZ-3) from the postsynaptic density (PSD-95) protein. This protocol employs careful and precise computational techniques for design of new strategy for predicting novel and potent peptides against PDZ protein. We attempted to generate a pharmacophore model using crystal structure of peptide inhibitor bound to the PDZ-3. A highly specific and sensitive generated pharmacophore model was used for screening virtual database generated using different combination of amino acid substitutions as well as decoy peptide database for its sensitivity and specificity. Identified hit peptides were further analyzed by docking studies, and their stability analyzed using solvated molecular dynamics. Quantum Mechanics/Molecular Mechanics (QM/MM) interaction energy and GMX-PBSA scoring schemes were used for ranking of stable peptides. Computational approach applied here generated encouraging results for identifying peptides against PDZ interaction model. The workflow can be further exercised as a virtual screening technique for reducing the search space for candidate target peptides against PDZ domains.