Binding Site Residues Research Articles

Binding of Selected Ligands to Human Protein Disulfide Isomerase and Microsomal Triglyceride Transfer Protein Complex and the Associated Conformational Changes: A Computational Molecular Modelling Study.

Human protein disulfide isomerase (PDI) is a multifunctional protein, and also serves as the β subunit of the human microsomal triglyceride transfer protein (MTP) complex, a lipid transfer machinery. Dysfunction of the MTP complex is associated with certain disease conditions such as abetalipoproteinemia and cardiovascular diseases. It is known that the functions of PDI or the MTP complex can be regulated by the binding of a small-molecule ligand to either of these two proteins. In the present study, the conformational changes of the MTP complex upon the binding of three selected small-molecule ligands (17β-estradiol, lomitapide and a phospholipid) are investigated based on the available biochemical and structural information by using the protein-ligand docking method and molecular dynamics (MD) simulation. The ligand-binding sites, the binding poses and binding strengths, the key binding site residues, and the ligand binding-induced conformational changes in the MTP complex are analyzed based on the MD trajectories. The open-to-closed or closed-to-open transitions of PDI is found to occur in both reduced and oxidized states of PDI and also independent of the presence or absence of small-molecule ligands. It is predicted that lomitapide and 1,2-diacyl-sn-glycero-3-phosphocholine (a phospholipid) can bind inside the lipid-binding pocket in the MTP complex with high affinities, whereas 17β-estradiol interacts with the lipid-binding pocket in addition to its binding to the interface region of the MTP complex. Additionally, lomitapide can bind to the b' domain of PDI as reported earlier for E2. Key residues for the ligand-binding interactions are identified in this study. It will be of interest to further explore whether the binding of small molecules can facilitate the conformational transitions of PDI in the future. The molecular and structural insights gained from the present work are of value for understanding some of the important biological functions of PDI and the MTP complex.

Synthesis, cytotoxicity, apoptosis-inducing activity and molecular docking studies of novel isatin-podophyllotoxin hybrids.

Podophyllotoxin, along with its numerous derivatives and related compounds, is well known for its broad-spectrum pharmacological activity, especially for anticancer potential. In this study, several isatin-podophyllotoxin hybrid compounds were successfully synthesized with good yields through microwave-prompted three-component reactions of 2-amino-1,4-naphthoquinone, various substituted isatins, and tetronic acid. Their cytotoxicity was assessed against four types of human cancer cell lines, HepG2 (hepatoma carcinoma), MCF7 (breast cancer), A549 (non-small lung cancer), and KB (epidermoid carcinoma), alongside nontumorigenic HEK-293 human embryonic kidney cells. Among 14 compounds screened, 7f possessed the strongest cytotoxicity to KB and A549 cell lines, with IC50 values of 1.99 ± 0.22 and 0.90 ± 0.09 μM, respectively. Further studies revealed that product 7f could arrest the cell cycle of A549 cells at S phase and induce apoptosis of A549 cells. This compound was examined for its binding ability against cyclin-dependent kinases (CDKs) and procaspase/caspase systems. The results indicated that 7f exhibited significant interactions with the residues of the ATP binding sites of CDK2/cyclin A and CDK5/p25 and also activated procaspase 6 through stable zinc chelation. Additionally, physicochemical and pharmacokinetic properties related to drug-likeness, in parallel with toxicity, were computationally assessed to identify the main issues that need to be addressed in structural optimization. Taken together, compound 7f was identified as a potent cytotoxic agent that could be considered for anticancer drug discovery and development.

Insights into the catalytic mechanism of archaeal peptidoglycan endoisopeptidases from methanogenic phages.

Archaeal peptidoglycan, a crucial component of the cell walls of Methanobacteria and Methanopyri, enhances the tightness of methanogenic cells and their resistance to known lytic enzymes and antibiotics. Although archaeal peptidoglycan endoisopeptidases (Pei) can reportedly degrade archaeal peptidoglycan, their biochemistry is still largely unknown. In this study, we investigated the activity and catalytic properties of the endoisopeptidases PeiW and PeiP using synthesized isopeptides identical to natural substrates. Enzymatic assays demonstrated their distinct substrate specificity and cleavage efficiency. The crystal structure of Pei revealed a catalytic mechanism resembling that of cysteine peptidases that use the 'CHD' triad to cleave isopeptide bonds. We also identified several key residues in the substrate binding site that confer recognition specificity, including Y174, V252 and C265. Based on the residues present in the active site and their influence on activity, we propose a classification of the archaeal peptidoglycan endoisopeptide family into four categories to facilitate the identification of new archaeal peptidases in the future. These insights into the structure and function of Pei suggest new strategies for use in methanogen biotechnology.

Design of bisamide inhibitors of the TASK-1 potassium channel in silico.

TWIK-related acid-sensitive potassium channel 1 (TASK-1) is expressed ubiquitously across various tissues and plays a significant role in neural activity and anesthetic modulation, making it a crucial target for pharmaceutical research. The high conservation of binding site residues within the TASK family, particularly between TASK-1 and TASK-3, necessitates the development of selective inhibitors for TASK-1. In this study, we utilized a combination of structure-based drug design (SBDD) and ligand-based drug design (LBDD) approaches. Initially, several bisamide-centered molecules were designed using the program MolAICal, which is recognized for its ability to generate selective inhibitors containing bisamide segments, and conducted preliminary screening via molecular docking. Subsequently, 3D-QSAR models were developed for 56 bisamide derivatives targeting TASK-1 and TASK-3, with the models exhibiting robust predictive capabilities (TASK-1: Q2 = 0.61, R2pred = 0.84; TASK-3: Q2 = 0.60, R2pred = 0.71). Using these models, the candidate molecules were subjected to activity prediction and subsequent filtering. Ultimately, molecular dynamics simulations, coupled with free energy calculations, pinpointed two bisamide-core molecules with favorable ADMET properties as potential selective inhibitors for TASK-1. Furthermore, molecular dynamics simulations revealed the critical role of the key residue Leu122 in conferring selectivity to bisamide compounds for TASK-1 channel proteins.

Phytochemical analysis, therapeutic and molecular docking studies for the compounds of wild type and Senkambu variants of curry leaves targeting HER2 Kinase domain a potential gastric cancer receptor

Gastric cancer is the fifth most notable health concern globally. In recent years, molecular docking, a computational technique, has emerged as tool in drug discovery. The present investigation aimed to identify the major bioactive compounds in the wild-type curry leaves found in the Shevaroy Hills and the local Senkambu variant from Karamadai. Virtual screening of 40 ligands from Curry leaves of wild type and Senkambu type was identified through GC-MS profiling. These compounds were targeted against HER2 Kinase domain which is a potential receptor for Gastric cancer. Information regarding the binding site residues for the receptor was predicted using CASTp server. Molecular docking was performed for HER2 kinase domain with the predicted compounds through GC-MS profiling. The top 3 hits reported with least binding affinity for the target protein were considered for further interaction analysis using Biovia Discovery studio visualizer. Upon analyzing the interacted compounds, the Piperine from Wild type curry leaves was found to have good interaction with HER2 Kinase domain by forming two hydrogen bonds and binding score of -8.3 k cal/mol. The current study might guide the designing of analogues of piperine in the evolution of effective broad spectrum drug development in cancer therapy.

Computational analysis of the deleterious non-synonymous single nucleotide polymorphisms (nsSNPs) in TYR gene impacting human tyrosinase protein and the protein stability.

Tyrosinase, a copper-containing oxidase, plays a vital role in the melanin biosynthesis pathway. Mutations in the tyrosinase gene can disrupt the hydroxylation of tyrosine, leading to decreased production of 3,4-dihydroxyphenylalanine (DOPA). Consequently, this impairs the subsequent formation of dopaquinone, a key precursor in melanin pigment synthesis. This study aimed to identify the deleterious non-synonymous single nucleotide polymorphisms (nsSNPs) within the TYR gene that exert an influence on the human TYR protein. Additionally, we evaluated the impact of 10 FDA-approved drugs on the protein stability of mutated structures, exploring the potential for inhibitory pharmaceutical interventions. Through various bioinformatics tools, we detected 47900 nsSNPs, particularly K142M, I151N, M179R, S184L, L189P, and C321R, which were found to be the most deleterious variants, decreasing the protein stability. These drugs (Sapropterin, Azelaic Acid, Menobenzone, Levodopda, Mequinol, Arbutin, Hexylresorcinol, Artenimol, Alloin and Curcumin) interacted with the binding sites in four mutant models K142M, I151N, M179R, and S184L proving that these ligands directly bind with the active site of mutant tyrosinase protein to inhibit it's working. On the other hand, two mutant models L189P and C321R did not show any binding site residue interaction with any ligands. In conclusion, this in-silico analysis of deleterious nsSNPs in the TYR gene, coupled with the evaluation of ligands/drugs on mutated tyrosinase structures not only advances our understanding of molecular variations but also highlights promising pathways for targeted inhibitory interventions in the intricate network of melanin biosynthesis.

Computational and theoretical insights into the cytotoxic prospects of compounds isolated from Elaeodendron buchananii against Leukemia

The present study investigated the cytotoxic prospects of isolated compounds from Elaeodendron buchananii against leukemia, using computational tools. Comprehensive literature searches revealed only buchaninoside, mutangin, methyl 3β-acetoxy-11α, 19α, 28-trihydroxyurs-12-en-23-oic acid, 3β, 11α, 19α-trihydroxyurs-12-en-23, 28-dioic acid, 3β-acetoxy-19α, 24, 28-trihydroxyurs-12-ene, 3-oxo-19α,28-dihydroxyurs-12-en-24-oic acid, and elabunin have been isolated from E. buchananii. The compounds were subjected to Density Functional Theory (DFT) and Molecular Dynamics (MD) analyses, with Fms-like tyrosine kinase (FLT3) and catalytic binding sites of Murine Leukemia Virus (MLV) as the target proteins in lukemia. Following DFT analysis, the structures of the compounds were optimized at the PW6B95D3/Def2-TZVP level of theory; their UV-Visible peaks were in the UV region, with mutangin, 3-oxo-19α,28-dihydroxyurs-12-en-24-oic acid and elabunin exhibiting one single peak. The potent Root-Mean-Square Deviation, Root-Mean-Square Fluctuation, solvent-accessible surface area and radius of gyration values indicated a strong and stable molecular interaction between the compounds and the proteins. These were further supported by high ∆G values, with MLV showing the best interaction. Per-residue decomposition plots also revealed high energy contributions in the interactions’ binding sites residues. These results indicate that the cytotoxic prospects of the isolated compounds against leukemia as indicated by its molecular interactions with FLT3 and MLV.

Irreversible Inhibition of DNMT3A by an N-Mustard Analog of S-Adenosyl-L-Methionine.

DNA methylation, an important epigenetic modification, is catalyzed by DNA methyltransferases and is essential in the regulation of gene expression. Here, the utility of an N-mustard analog designed to mimic the native methyl donor, S-adenosyl-L-methionine (SAM), was explored with the DNA methyltransferase 3A catalytic domain (DNMT3AC). In lieu of the expected analog transfer to DNA, methyltransferase activity was instead inhibited in a concentration dependent manner. Further investigation into the mechanism of analog inhibition did not reveal a typical competitive mechanism. Instead, mass spectrometry analysis provided direct evidence of two cysteine residues in the SAM binding site covalently modified by the SAM analog and confirmed its' function as an irreversible inhibitor of DNMT3AC.

Phytoconstituents of Withania somnifera (L.) Dunal (Ashwagandha) unveiled potential cerebroside sulfotransferase inhibitors: insight through virtual screening, molecular dynamics, toxicity, and reverse pharmacophore analysis

Cerebroside sulfotransferase (CST) is considered as therapeutic target for substrate reduction therapy (SRT) for metachromatic leukodystrophy (MLD). The present study evaluates the therapeutic potential of 57 phytoconstituents of Withania somnifera against CST. Using binding score cutoff ≤-7.0 kcal/mol, top 10 compounds were screened and after ADME and toxicity-based screening, Withasomidienone, 2,4-methylene-cholesterol, and 2,3-Didehydrosomnifericin were identified as safe and potent drug candidates for CST inhibition. Key substrate binding site residues involved in interaction were LYS82, LYS85, SER89, TYR176, PHE170, PHE177. Four steroidal Lactone-based withanolide backbone of these compounds played a critical role in stabilizing their position in the active site pocket. 100 ns molecular dynamics simulation and subsequent trajectory analysis through structural deviation and compactness, principal components, free energy landscape and correlation matrix confirmed the stability of CST-2,3-Didehydrosomnifericin complex throughout the simulation and therefore is considered as the most potent drug candidate for CST inhibition and Withasomidienone as the second most potent drug candidate. The reverse pharmacophore analysis further confirmed the specificity of these two compounds towards CST as no major cross targets were identified. Thus, identified compounds in this study strongly present their candidature for oral drug and provide route for further development of more specific CST inhibitors.

Investigating the Effect of GLU283 Protonation State on the Conformational Heterogeneity of CCR5 by Molecular Dynamics Simulations.

CCR5 is a class A GPCR and serves as one of the coreceptors facilitating HIV-1 entry into host cells. This receptor has vital roles in the immune system and is involved in the pathogenesis of different diseases. Various studies were conducted to understand its activation mechanism, including structural studies in which inactive and active states of the receptor were determined in complex with various binding partners. These determined structures provided opportunities to perform molecular dynamics (MD) simulations and to analyze conformational changes observed in the protein structures. The atomic-level dynamic studies allow us to explore the effects of ionizable residues on the receptor. Here, our aim was to investigate the conformational changes in CCR5 when it forms a complex with either the inhibitor maraviroc (MRV), an approved anti-HIV drug, or HIV-1 envelope protein GP120, and compare these changes to the receptor's apo form. In our simulations, we considered both ionized and protonated states of ionizable binding site residue GLU2837.39 in CCR5 as the protonation state of this residue was considered ambiguously in previous studies. Our molecular simulations results suggested that in fact, the change in the protonation state of GLU2837.39 caused interaction profiles to be different between CCR5 and its binding partners, GP120 or MRV. We observed that when the protonated state of GLU2837.39 was considered in complex with the envelope protein GP120, there were substantial structural changes in CCR5, indicating that it adopts a more active-like conformation. On the other hand, CCR5 in complex with MRV always adopted an inactive conformation regardless of the protonation state. Hence, the CCR5 coreceptor displays conformational heterogeneity not only depending on its binding partner but also influenced by the protonation state of the binding site binding site residue GLU2837.39. This outcome is also in accordance with some studies showing that GP120 binding could activate signaling pathways. This outcome could also have significant implications for discovering novel CCR5 inhibitors as anti-HIV drugs using in silico methods such as molecular docking, as it may be necessary to consider the protonated state of GLU2837.39.

Unveiling MurM inhibitors in Enterococcus faecalis V583: a promising approach to tackle antibiotic resistance

Enterococcus faecalis is commonly found in the GI tract of humans and animals. It causes various infections, especially in hospital environments, and shows growing antibiotic resistance. This study utilized a subtractive proteomics approach to find out the potential drug targets in E. faecalis. Unique metabolic pathways were analysed and compared to the host to minimize adverse effects. Among twenty nine pathogenic specific and seventy three host-pathogen common pathways identified using the KEGG database, sixty seven essential proteins were found through the DEG BLAST search. PSORTB predicted that forty cytoplasmic proteins could be suitable as druggable targets. Further analysis identified fourteen proteins with virulence properties using the VFDB BLAST. Among these, seven proteins with more than ten antigenic sites were subjected to DrugBank BLAST, identifying three novel and four existing drug targets. One of the crucial drug targets, MurM, was selected due to its critical role in peptidoglycan biosynthesis. The reason for selecting MurM is crucial for addressing antibiotic resistance, disrupting bacterial cell wall synthesis, and attaining selective antimicrobial activity. MurM belongs to the mixed αβ class with two functional domains. The possible binding site residues of MurM are Trp31, Lys35, Trp38, Arg215, and Tyr219. Virtual screening identified potential lead candidates for MurM, and four were selected based on their physiochemical, pharmacokinetic, and structural properties. This study provides valuable insights into identifying and analysing a potential drug target, the MurM protein, and its inhibitors in E. faecalis V583.

Mechanistic insights into β-glucuronidase inhibition by isoprenylated flavonoids from Centaurea scoparia: Bridging experimental and computational approaches

Exploring the intricate mechanisms underlying the inhibition of β-glucuronidase, particularly the enzyme produced by gastrointestinal bacteria, is essential for the development of novel therapeutic interventions and pharmaceutical advancements. Inhibiting bacterial β-glucuronidase can prevent the reactivation of harmful drug metabolites in the gut, thereby reducing associated toxicities. The inhibitory potential of four isoprenylated flavonoids from Centaurea scoparia against β-glucuronidase was intensively investigated by combined in vitro and in silico tools. The results of in vitro inhibitory activity showed that Compounds 1 and 3 exhibited the highest inhibitory activity, as evidenced by their low IC50 values of 3.65 ± 0.28 and 3.97 ± 0.11 μM, respectively. The enzyme kinetics analysis revealed that Compound 1 and the positive control drug (EGCG) follow a mixed inhibition mechanism. In contrast, compounds 3 and 4 exhibited a noncompetitive inhibition mode, as suggested by the intersection patterns observed in the Lineweaver-Burk plots. The docking studies corroborated the in vitro inhibitory assay results, with Compounds 1 and 3 demonstrating the lowest binding affinities and the highest extent of polar and hydrophobic interactions with residues in the enzyme's binding site. Utilizing a 200 ns molecular dynamics (MD) simulation, we examined the interaction dynamics between isolated flavonoids and β-glucuronidase. The analysis of multiple MD parameters revealed that compounds 1 and 3 maintained consistent trajectories and demonstrated significant energy stabilization when interacting with β-glucuronidase. These computational findings align with experimental results, highlighting the potential of compounds 1 and 3 as potent inhibitors for β-glucuronidase.

Identification of steroidal cardenolides from Calotropis procera as novel HIV-1 PR inhibitors: A molecular docking & molecular dynamics simulation study.

Background & objectives Despite advancements in antiretroviral therapy, drug-resistant strains of HIV (human immunodeficiency virus) remain a global health concern. Natural compounds from medicinal plants offer a promising avenue for developing new HIV-1 PR (protease) inhibitors. This study aimed to explore the potential of compounds derived from Calotropis procera, a medicinal plant, as inhibitors of HIV-1 PR. Methods This in silico study utilized natural compound information and the crystal structure of HIV-1 PR. Molecular docking of 17 steroidal cardenolides from Calotropis procera against HIV-1 PR was performed using AutoDock 4.2 to identify compounds with higher antiviral potential. A dynamic simulation study was performed to provide insights into the stability, binding dynamics, and potential efficacy of the top potential antiviral compound as an HIV-1 therapeutic. Results We found that all tested cardenolides had higher binding affinities than Amprenavir, indicating their potential as potent HIV-1 PR inhibitors. Voruscharin and uscharidin displayed the strongest interactions, forming hydrogen bonds and hydrophobic interactions with HIV-1 PR. Voruscharin showed improved stability with lower RMSD (Root Mean Square Deviation) values and reduced fluctuations in binding site residues but increased flexibility in certain regions. The radius of gyration analysis confirmed a stable binding pose between HIV-1 PR and Voruscharin. Interpretation & conclusions These findings suggest that Calotropis procera could potentially be a source of compounds for developing novel HIV-1 PR inhibitors, contributing to the efforts to combat HIV. Further studies and clinical trials are needed to evaluate the safety and efficacy of these compounds as potential drug candidates for the treatment of HIV-1 infection.

Repurposing FDA-approved drugs for combating tigecycline resistance in Acinetobacter baumannii: in silico screening against BaeR protein.

Acinetobacter baumannii is becoming a gravely threatening nosocomial infection with a higher mortality rate. The present study targets the BaeR protein that mediates resistance to tigecycline antibiotics. The BaeR protein, along with the aid of BaeS, senses the incoming antibiotics and stimulates the expression of resistance proteins. These resistance proteins efflux the antibiotics and protect the cells from its effect. The main goal of the current study is to determine potential inhibitors from already existing FDA-approved drugs that could mitigate the BaeR protein. A range of in silico approaches, including molecular dynamics, virtual screening, SIFT analysis, ADMET, DFT, MM/GBSA, MMPBSA and per residue interaction analysis, were performed to identify inhibitors against this protein. The screening of FDA-approved compounds against the BaeR protein yielded 620 compounds. These compounds were clustered by SIFT to distinguish related compounds, it resulted in 20 different clusters. The top five clusters that can accommodate the binding site with better interaction and score by fulfilling all criteria were selected. The DFT analysis showed a smaller energy gap among all the compounds, indicating the ability of the compound to form firm interactions. All the compounds showed less binding free energy in both MM/GBSA and MM/PBSA analyses. The compounds were observed to be stable throughout the simulation. The per-residue interaction analysis confirmed that interactions with binding site residues were stable throughout the simulation. As a result of the study, four compounds, namely ZINC000003801919, DB01203, DB11217 and ZINC0000000056652, were identified as efficient candidates to deal with antimicrobial resistance in A. baumannii.

Characterisation of seryl tRNA synthetase (srs-2) in Haemonchus contortus and Teladorsagia circumcincta

The aim of the study was to purify and characterise recombinant proteins with the potential as an anti-parasite vaccine. Full-length cDNAs encoding seryl-tRNA synthetase (srs-2) were cloned from Haemonchus contortus (HcSRS-2) and Teladorsagia circumcincta (TcSRS-2). TcSRS-2 and HcSRS-2 cDNA (1458bp) encoded proteins of 486 amino acids, each of which was present as a single band of about 55 kDa on SDS-PAGE. Multiple alignments of the protein sequences showed homology of 94% between TcSRS-2 and HcSRS-2, 76–93% with SRS-2s of eight nematodes and 68% with Mus musculus SRS-2. The predicted three-dimensional structures revealed an overall structural homology of TcSRS-2 and HcSRS-2, highly conserved binding and catalytic sites, and minor differences in the tautomerase binding site residues in other nematode SRS-2 homologues. A phylogenetic tree was constructed using helminth and mammalian SRS-2 sequences. Soluble C-terminal SRS-2 proteins were expressed in Escherichia coli strain AY2.4 and purified. Recombinant HcSRS-2 assay shows that the recombinant enzyme was active and stable. The Km and Vmax for ATP were 3.9 ± 1.0 μM and 2.7 ± 0.1 μmol min−1 mg−1 protein, respectively. Antibodies in serum and saliva from field-immune, but not nematode-naïve, sheep recognised recombinant HcSRS-2 and TcSRS-2 in enzyme-linked immunosorbent assays. Recognition of the recombinant proteins by antibodies generated by exposure of sheep to the native enzyme indicates similar antigenicity of the two proteins.

Bioinformatics and computational studies of chabamide F and chabamide G for breast cancer and their probable mechanisms of action

Globally, the prevalence of breast cancer (BC) is increasing at an alarming level, despite early detection and technological improvements. Alkaloids are diverse chemical groups, and many within this class have been reported as potential anticancer compounds. Chabamide F (F) and chabamide G (G) are two dimeric amide alkaloids found in a traditional medicinal plant, Piper chaba, and possess significant cytotoxic effects. However, their scientific rationalization in BC remains unknown. Here, we aimed to investigate their potential and molecular mechanisms for BC through in silico approaches. From network pharmacology, we identified 64 BC-related genes as targets. GO and KEGG studies showed that they were involved in various biological processes and mostly expressed in BC-related pathways such as RAS, PI3K-AKT, estrogen, MAPK, and FoxO pathways. However, PPI analysis revealed SRC and AKT1 as hub genes, which play key roles in BC tumorigenesis and metastasis. Molecular docking revealed the strong binding affinity of F (− 10.7 kcal/mol) and G (− 9.4 and − 11.7 kcal/mol) for SRC and AKT1, respectively, as well as the acquisition of vital residues to inhibit them. Their long-term stability was evaluated using 200 ns molecular dynamics simulation. The RMSD, RMSF, Rg, and SASA analyses showed that the G-SRC and G-AKT1 complexes were excellently stable compared to the control, dasatinib, and capivasertib, respectively. Additionally, the PCA and DCCM analyses revealed a significant reduction in the residual correlation and motions. By contrast, the stability of the F-SRC complex was greater than that of the control, whereas it was moderately stable in complex with AKT1. The MMPBSA analysis demonstrated higher binding energies for both compounds than the controls. In particular, the binding energy of G for SRC and AKT1 was − 120.671 ± 16.997 and − 130.437 ± 19.111 kJ/mol, respectively, which was approximately twice as high as the control molecules. Van der Waal and polar solvation energies significantly contributed to this energy. Furthermore, both of them exhibited significant interactions with the binding site residues of both proteins. In summary, this study indicates that these two molecules could be a potential ATP-competitive inhibitor of SRC and an allosteric inhibitor of AKT1.

Innovative Mamba and graph transformer framework for superior protein-ligand affinity prediction

Accurate measurement of the affinity between drug targets is of great importance in the field of drug discovery, as it provides significant information about drug action. However, the diverse representations of compounds and proteins increase the difficulty of this task. Recently, numerous studies have explored the application of various deep learning methods and architectures in Drug-Target Affinity (DTA) prediction tasks, continuously improving performance. In this work, we present a multimodal late-stage fusion prediction model called MGDTA that combines the pioneering Mamba architecture with a hybrid model architecture based on graph transformers and transformers (GTT). For molecular compounds, structural characterisation is extracted using GTT layers. Additionally, a feature integration module is utilized to obtain molecular features from derived molecular fingerprints, culminating in a unified representation of drug molecules through a feature integration module. For proteins, we adopt the state-space model (SSM) module Mamba to learn feature representations from target protein sequences. This represents a preliminary exploration of the Mamba module in protein sequence feature extraction in DTA tasks. Finally, we predict the results by integrating the feature representations of both molecules and proteins. Our results indicate that the proposed model achieves superior performance on benchmark datasets, validating the feasibility of using the Mamba architecture for DTA tasks. Furthermore, sequence feature visualization demonstrates the capability of the Mamba block to focus on binding site residues within drug targets.

CC chemokine receptor 2 is allosterically modulated by sodium ions and amiloride derivatives through a distinct sodium ion binding site

CC chemokine receptor 2 and CCL2 are highly involved in cancer growth and metastasis, and immune escape. Raised sodium ion concentrations in solid tumours have also been correlated to metastasis and immune modulation. Sodium ions can modulate class A G protein-coupled receptors through the sodium ion binding site characterized by a highly conserved aspartic acid residue (D2.50), also present in CCR2. Hence, we further explored this binding site in CCR2 by radioligand binding studies and mutagenesis. Modulation of three distinctly binding radioligands by sodium ions and amiloride derivates was investigated. Sodium ions were observed to be relatively weak modulators of antagonist binding, but substantially increased 125I-CCL2 dissociation from CCR2. 6-Substituted Hexamethylene Amiloride (HMA) modulated all tested radioligands. Induced-fit docking of HMA in the presumed sodium ion binding site of CCR2 confirmed its binding site. Finally, investigation of (cancer-associated) mutations in the sodium ion binding site showed a markedly decreased expression compared to wild type. Only two mutants, G123A3.35 and G127K3.39, were able to be bound by [3H]INCB3344 and [3H]CCR2-RA-[R]. Thus, mutagenesis showed that the sodium ion binding site residues, which are distinct from other class A GPCRs and related to chemokine receptor evolution, are crucial for receptor integrity. Moreover, the tested mutations appeared to have no effect on modulation observed by HMA or a minor effect on sodium chloride modulation on the tested radioligands. All in all, these results invite further exploration of the CCR2 sodium ion binding site in (cancer) biology, and potentially as a third druggable binding site.

Phytochemical characterization, antimicrobial properties and in silico modeling perspectives of Anacyclus pyrethrum essential oil

Medicinal plants are used widely in the treatment of various infectious diseases. One of these medical plants is Moroccan plants such as Anacyclus pyrethrum. In this study, the essential oil isolated from the leaves of Anacyclus pyrethrum (APEO) by the hydrodistillation method was analyzed using (GC/MS) analysis. A total of forty-four compounds were identified form the oil and the oxygenated monoterpenes were the most abundant class of compounds. The major identified compound is santolina alcohol (40.7 %), followed by germacrene-D (8.9 %). The in-vitro assessment of the antimicrobial efficacy of APEO encompassed an investigation involving six microbial strains, including two Gram-positive bacteria, four Gram-negative bacteria, and three fungal strains. The findings revealed noteworthy antibacterial and antifungal properties against all examined microorganisms, with inhibitory zone diameters ranging from 25.67 ± 0.06 mm to 25.19 ± 0.03 mm for Gram-positive bacteria and from 22.34 ± 0.01 mm to 14.43 ± 0.02 mm for Gram-negative bacteria, as determined through the disc-diffusion assay. In the case of antifungal activity, inhibitory zones ranged from 24.57 ± 0.04 mm to 18.37 ± 0.06 mm. Further evaluation revealed that the MIC values of Gram-positive bacteria were at the concentration 0.25 % v/v, while MBC values ranged from 0.25 % to 1.0 % v/v. The Gram-negative bacteria exhibited MIC values spanning from 0.5 % to 2.0 % v/v, with MBC values in the range of 0.5 %–2.0 % v/v. For the fungal strains, MIC values ranged from 0.5 % to 1.0 % v/v, while the MFC consistently remained at 1.0 % for all tested fungal strains. The assessment of the MBC/MIC and MFC/MIC ratios collectively indicates that A. pyrethrum EO possesses bactericidal and fungicidal attributes. The in silico study of bioavailability predictions for compounds in APEO based on six physicochemical properties show optimal physiochemical properties including size, lipophilicity, solubility, flexibility, and saturation. α-Pinene, limonene, germacrene D, and (E)-β-farnesene are non-polar due to their lack of polar groups, and the ADME profile indicates desirable properties for considering these compounds in drug development. Molecular docking investigation indicates that all the compounds of APEO reside well into the binding site of the DNA gyrase B enzyme of Staphylococcus aureus by mediating a number of significant interactions with the binding site residues. The ADME analysis suggested that the major compounds APEO possess desirable properties for further consideration in drug development. In light of these findings, APEO could serve as a natural source for the elaboration of new and active antimicrobial drugs.

Direct measurements of neurosteroid binding to specific sites on GABAA receptors.

Neurosteroids are allosteric modulators of GABAA currents, acting through several functional binding sites although their affinity and specificity for each site are unknown. The goal of this study was to measure steady-state binding affinities of various neurosteroids for specific sites on the GABAA receptor. Two methods were developed to measure neurosteroid binding affinity: (1) quenching of specific tryptophan residues in neurosteroid binding sites by the neurosteroid 17-methylketone group, and (2) FRET between MQ290 (an intrinsically fluorescent neurosteroid) and tryptophan residues in the binding sites. The assays were developed using ELIC-α1GABAAR, a chimeric receptor containing transmembrane domains of the α1-GABAA receptor. Tryptophan mutagenesis was used to identify specific interactions. Allopregnanolone (3α-OH neurosteroid) was shown to bind at intersubunit and intrasubunit sites with equal affinity, whereas epi-allopregnanolone (3β-OH neurosteroid) binds at the intrasubunit site. MQ290 formed a strong FRET pair with W246, acting as a site-specific probe for the intersubunit site. The affinity and site-specificity of several neurosteroid agonists and inverse agonists was measured using the MQ290 binding assay. The FRET assay distinguishes between competitive and allosteric inhibition of MQ290 binding and demonstrated an allosteric interaction between the two neurosteroid binding sites. The affinity and specificity of neurosteroid binding to two sites in the ELIC-α1GABAAR were directly measured and an allosteric interaction between the sites was revealed. Adaptation of the MQ290 FRET assay to a plate-reader format will enable screening for high affinity agonists and antagonists for neurosteroid binding sites.