Molecular Docking, Molecular Dynamics Simulations, and Free Energy Calculation Insights into the Binding Mechanism between VS-4718 and Focal Adhesion Kinase

Molecular Mechanism of Cynodon dactylon Phytosterols Targeting MAPK3 and PARP1 to Combat Epithelial Ovarian Cancer: A Multifaceted Computational Approach

Epithelial Ovarian Cancer (EOC) presents a global health concern, necessitating the development of innovative therapeutic strategies to combat its impact. This study was employed to investigate the unexplored therapeutic efficacy of Cynodon dactylon phytochemicals against EOC using a multifaceted computational approach. A total of 19 out of 89 rigorously curated phytochemicals were assessed as potential drug targets via ADMET profiling, while protein-protein interaction analysis scrutinized the top 20 hub genes among 264 disease targets, revealing their involvement in cancer-related pathways and underscoring their significance in EOC pathogenesis. In molecular docking, Stigmasterol acetate showed the highest binding affinity (-10.9 kcal/mol) with Poly [ADP-ribose] polymerase-1 (PDB: 1UK1), while Arundoin and Beta-Sitosterol exhibited strong affinities (-10.4 kcal/mol and -10.1 kcal/mol, respectively); additionally, Beta-Sitosterol interacting with Mitogen-activated protein kinase 3 (PDB: 4QTB) showed a binding affinity of -10.1 kcal/mol, forming 2 hydrogen bonds and a total of 10 bonds with 10 residues. Molecular dynamics simulations exhibited the significant structural stability of the Beta-Sitosterol-4QTB complex with superior binding free energy (-36.61 kcal/mol) among the three complexes. This study identified C. dactylon phytosterols, particularly Beta-Sitosterol, as effective in targeting MAPK3 and PARP1 to combat EOC, laying the groundwork for further experimental validation and drug development efforts.

Evaluating binding and interaction of selected pesticides with serum albumin proteins by Raman, 1H NMR, mass spectrometry and molecular dynamics simulation

Addressing the acute pesticide poisoning and toxicity to humans, is a global challenge of top priority. Serum albumin is the most abundant plasma protein, capable of binding with herbicide and pesticide residues. This study reports multifaceted approaches for in-depth and robust investigation of the molecular interactions of selected pesticides, including propanil (PPL), bromoxynil (BXL), metolachlor (MLR) and glyphosate (GPE) with bovine serum albumin (BSA) proteins using experimental (Raman and FTIR spectroscopy, native mass spectrometry and high field 1H NMR), molecular dynamics (MD) simulation and principal component analysis (PCA). The binding of pesticides with BSA resulted in BSA amide I and amide II Raman spectral shifts. PCA of Raman spectra of serum-pesticide complexes showed the grouping of pesticides on the score plot based on the similarities and differences in pesticides' chemical structures. Native mass spectrometry results revealed strong adduct formation of the pesticides with the protein. The observed changes in chemical shifts, peak broadening or peak disappearance of characteristic proton signals of the pesticides, indicated altered chemical environments due to binding BSA-pesticides interactions. The results of MD simulation conducted for over 500 ns revealed strong pesticides interaction with LEU197, LEU218, LEU237, TRP213, SER286 and ILE289 residues to the site I of BSA. Free energy landscapes provided insights into the conformational changes in BSA on the binding of pesticides. Overall, the experimental and computational results are in consonant and indicate the binding of pesticides into the site I and site II (sub-domain IIA) of the BSA via hydrogen bonding, non-covalent and hydrophobic interactions.Communicated by Ramaswamy H. Sarma.

Bioinformatics and immunoinformatics assisted multiepitope vaccine construct against Burkholderia anthina

Burkholderia anthina is a pathogenic bacterial species belonging to the Burkholderiaceae family and it is mainly considered the etiological agent of chronic obstructive pulmonary diseases associated with cystic fibrosis, due to being intrinsic antibiotic resistant making it difficult to treat pulmonary infections. Hence increased rate of antibiotic-resistant bacterial species vaccine development is the priority to tackle this problem. In research work, we designed a multi-epitope-based vaccine construct against B. anthina using reverse vaccinology immunoinformatics and biophysical approaches. Based on the subtractive proteomic screening of core proteins we identified 3 probable antigenic proteins and good vaccine targets namely, type VI secretion system tube protein hcp Burkholderia, fimbria/pilus periplasmic chaperone and fimbrial biogenesis outer membrane usher protein. The selected 3 proteins were used for B and B cells B-derived T-cell epitopes prediction. In epitopes prediction, different epitopes were predicted with various lengths and percentile scores and subjected to further immunoinformatics analysis. In immunoinformatics screening a total number of 06, IDDGNANAL, KTVKPDPRY, SEVESGSAP, YGGDLTVEV, SVSHDTNGR, and GSKADGYQR epitopes were considered good vaccine target candidates and shortlisted for vaccine construct designing. The vaccine construct was designed by joining selected epitopes with the help of a GPGPG linker and additionally linked with cholera toxin b subunit adjuvant to increase the efficacy of the vaccine construct the sequence of the said adjuvant were retrieved from protein data bank through its (PDB ID: 5ELD). The designed vaccine construct was evaluated for its physiochemical properties analysis in which we reported that the vaccine construct comprises 216 amino acids with a molecular weight of 22.37499 kilo Dalton, 15.55 instability index (II) is computed, and this classifies that the vaccine construct is properly stable. VaxiJen v2.0 web server predicted that the vaccine construct is probable antigenic in nature with 0.6320 predicted value. Furthermore AllerTOP v. 2.0 tool predicted that the designed vaccine construct is non allergic in nature. Molecular docking analysis was done for analysis of the binding affinity of the vaccine construct with TLR-2 (PDB ID: 6NIG), the docking results predicted 799.2 kcal/mol binding energy score that represents the vaccine construct has a good binding ability with TLR-2. Moreover, molecular dynamic simulation analysis results revealed that the vaccine construct and immune cell receptor has proper binding stability over various environmental condition, i.e. change in pressure range, temperature, and motion. After each analysis, we observed that the vaccine construct is safe stable, and probably antigenic and could generate an immune response against the target pathogen but in the future, experimental analysis is still needed to verify in silico base results.

Unraveling the therapeutic potential of Senna singueana phytochemicals to attenuate pancreatic cancer using protein-protein interactions, molecular docking, and MD simulation

Pancreatic cancer (PC) presents challenges due to limited treatment options and patients seek complementary therapies alongside conventional treatments to improve well-being. This study uses computational drug discovery approaches to find potential phytochemicals from S. singueana for PC treatment. Among the 38 phytochemicals screened from S. singueana, specific inhibitors against PC were selected. Protein-protein interaction (PPI) network analysis highlighted key targets with high degrees, including PTEN (8) and PTK2 (7) genes, along with their respective proteins 5BZX and 3BZ3, which were employed for molecular docking studies. 1-methylchrysene and 3-methyl-1,8,9-anthracenetriol showed strong binding affinities of - 9.2 and - 8.1 Kcal/mol, respectively. Molecular dynamics simulations lasting 300 ns assessed structural stability and interaction energy of compound-target dockings: 1-methylchrysene-PTEN and 3-methyl-1,8,9-anthracenetriol-PTK2. In molecular dynamics simulations, the 3-methyl-1,8,9-anthracenetriol-PTK2 complex showed lower RMSD, RMSF, radius of gyration, solvent-accessible surface area, and more hydrogen bonds than the 1-methylchrysene-PTEN complex. The 3-methyl-1,8,9-anthracenetriol-PTK2 complex exhibited significantly stronger binding with a binding free energy (ΔG) of - 21.92 kcal/mol compared to the less favourable ΔG of - 10.65 kcal/mol for the 1-methylchrysene-PTEN complex. The consistent and stable binding interaction observed in the 3-methyl-1,8,9-anthracenetriol-PTK2 complex highlights its potential as a potent inhibitor of Focal Adhesion Kinase 1. Consequently, it emerges as a promising lead compound for the development of pancreatic cancer therapeutics. Conversely, the fluctuations observed in the 1-methylchrysene-PTEN complex indicate a less stable binding interaction. This indicates the potential of 3-methyl-1,8,9-anthracenetriol as a primary candidate for pancreatic cancer treatment. These findings improve our grasp of S. singueana's multi-target effects and its promise in addressing pancreatic cancer. Nevertheless, additional in-vivo and in-vitro studies are required to validate their effectiveness and therapeutic potential.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40203-023-00179-9.

Identification of defactinib derivatives targeting focal adhesion kinase using ensemble docking, molecular dynamics simulations and binding free energy calculations

Focal adhesion kinase (FAK) belongs to the nonreceptor tyrosine kinases, which selectively phosphorylate tyrosine residues on substrate proteins. FAK is associated with bladder, esophageal, gastric, neck, breast, ovarian and lung cancers. Thus, FAK has been considered as a potential target for tumor treatment. Currently, there are six adenosine triphosphate (ATP)-competitive FAK inhibitors tested in clinical trials but no approved inhibitors targeting FAK. Defactinib (VS-6063) is a second-generation FAK inhibitor with an IC50 of 0.6 nM. The binding model of VS-6063 with FAK may provide a reference model for developing new antitumor FAK-targeting drugs. In this study, the VS-6063/FAK binding model was constructed using ensemble docking and molecular dynamics simulations. Furthermore, the molecular mechanics/generalized Born (GB) surface area (MM/GBSA) method was employed to estimate the binding free energy between VS-6063 and FAK. The key residues involved in VS-6063/FAK binding were also determined using per-residue energy decomposition analysis. Based on the binding model, VS-6063 could be separated into seven regions to enhance its binding affinity with FAK. Meanwhile, 60 novel defactinib-based compounds were designed and verified using ensemble docking. Overall, the present study improves our understanding of the binding mechanism of human FAK with VS-6063 and provides new insights into future drug designs targeting FAK.Communicated by Ramaswamy H. Sarma.

Volatile Organic Compounds Adsorption Capacities of Zeolite/Activated Carbon Composites Formed by Electrostatic Self-Assembly

Activated carbon (AC) is a broad-spectrum adsorbent but is flammable and has low adsorption capacities for polar and/or high-boiling volatile organic compounds (VOCs), while zeolites exhibit high thermal stability but poor adsorption of macromolecular and nonpolar VOCs. In this study, zeolite/AC composites were synthesized with the aim of obtaining broad-spectrum, efficient, and safe adsorbents for VOCs. Dimethyldiallylammonium chloride (DDA)-modified AC was used as a carrier for an in situ hydrothermal reaction enabling assembly with zeolites due to electrostatic attraction. Interface models were constructed for their phases, which revealed the binding force and simulated the binding process. The adsorption and flame resistance of the composites were evaluated. The results showed that DDA effectively modified AC to give it a long-lasting positive charge in solutions. High-silicon and pure-silicon zeolites exhibited low negative charges or were even neutral; it was difficult to combine with the modified AC via electrostatic attractions. Instead, LTA zeolites with high aluminum contents and negative charges were used, and the seed-induction method was used. Ethanol and ultrasonic dispersion were used to prevent agglomeration of the seeds and modified AC powder, so they were self-assembled electrostatically. Moreover, the crystallization time was extended and composites with high zeolite loadings were successfully prepared. According to the model calculation, the binding energy between the zeolite and AC before and after the DDA modification were 324.97 and 1076.46 kcal mol-1, respectively, and the distance between them was shortened by 2.7 Å after DDA treatment. As a result, AC and zeolite combined more closely and exhibited a stronger binding energy. The adsorption capacity for highly polar dichloromethane was improved by zeolite loading on the AC, and the bed penetration time was doubled. However, impregnation with inorganic sodium enhanced the reactivities of the organic components in the composite, and the ignition point was slightly reduced. Furthermore, the electrostatic self-assembly method can expand to prepare the LTA zeolite/columnar AC composite from shaped AC, greatly improving its application prospects.

Identification of abemaciclib derivatives targeting cyclin-dependent kinase 4 and 6 using molecular dynamics, binding free energy calculation, synthesis, and pharmacological evaluation

CDK4/6 plays a crucial role in various cancers and is an effective anticancer drug target. However, the gap between clinical requirements and approved CDK4/6 drugs is unresolved. Thus, there is an urgent need to develop selective and oral CDK4/6 inhibitors, particularly for monotherapy. Here, we studied the interaction between abemaciclib and human CDK6 using molecular dynamics simulations, binding free energy calculations, and energy decomposition. V101 and H100 formed stable hydrogen bonds with the amine-pyrimidine group, and K43 interacted with the imidazole ring via an unstable hydrogen bond. Meanwhile, I19, V27, A41, and L152 interacted with abemaciclib through π-alkyl interactions. Based on the binding model, abemaciclib was divided into four regions. With one region modification, 43 compounds were designed and evaluated using molecular docking. From each region, three favorable groups were selected and combined with each other to obtain 81 compounds. Among them, C2231-A, which was obtained by removing the methylene group from C2231, showed better inhibition than C2231. Kinase profiling revealed that C2231-A showed inhibitory activity similar to that of abemaciclib; additionally, C2231-A inhibited the growth of MDA-MB-231 cells to a greater extent than did abemaciclib. Based on molecular dynamics simulation, C2231-A was identified as a promising candidate compound with considerable inhibitory effects on human breast cancer cell lines.

Decellularized brain extracellular matrix slice glioblastoma culture model recapitulates the interaction between cells and the extracellular matrix without a nutrient-oxygen gradient interference

Decellularized extracellular matrix (dECM) is a valuable tool for generating three-dimensional in vitro tumor models that effectively recapitulate tumor-extracellular matrix (ECM) interactions. However, in current culture models, the components and structures of dECM are enzymatically disrupted to form hydrogels, making it difficult to recapitulate the native ECM. Additionally, when studying ECM-cell interactions, large-volume tumor culture models are incompatible with traditional experimental techniques and the nutrient-oxygen concentration gradient, which is a significant confounding factor. To address these issues, we developed a decellularized brain extracellular matrix slice (dBECMS) glioblastoma (GBM) culture model. This model possesses good light transmittance and substance diffusivity, making it compatible with traditional experimental techniques without forming nutrient-oxygen concentration gradients. Through transcriptomic analysis, we found that native brain ECM has a broad impact on glioma cells; the impact involves the ECM-ECM receptor interactions and the ECM and metabolic reprogramming. Further experiments demonstrated that dBECMS promoted glucose consumption and lactate production in GBM cells. Silver staining experiments revealed abundant proteins in the media of dBECMS, suggesting the degradation of the brain ECM by GBM cells. Transcriptome analysis also showed that the dBECMS-GBM culture model more accurately recapitulated the transcriptional profile of GBM than the two-dimensional culture. We experimentally demonstrated that the dBECMS-GBM model enhanced the resistance of GBM cells to temozolomide and increased the stemness of GBM cells. Additionally, we demonstrated the feasibility of the dBECMS-GBM model as a platform for drug response modeling. STATEMENT OF SIGNIFICANCE: The decellularized brain extracellular matrix (ECM) slice glioblastoma culture model mimics the interaction between native brain ECM and glioblastoma when glioblastoma infiltrates the brain and reveals the effects of native brain ECM on glioblastoma metabolism, ECM reprogramming, drug responsiveness, and stemness.

Interactions between curcumin and human salt-induced kinase 3 elucidated from computational tools and experimental methods

Natural products are widely used for treating mitochondrial dysfunction-related diseases and cancers. Curcumin, a well-known natural product, can be potentially used to treat cancer. Human salt-induced kinase 3 (SIK3) is one of the target proteins for curcumin. However, the interactions between curcumin and human SIK3 have not yet been investigated in detail. In this study, we studied the binding models for the interactions between curcumin and human SIK3 using computational tools such as homology modeling, molecular docking, molecular dynamics simulations, and binding free energy calculations. The open activity loop conformation of SIK3 with the ketoenol form of curcumin was the optimal binding model. The I72, V80, A93, Y144, A145, and L195 residues played a key role for curcumin binding with human SIK3. The interactions between curcumin and human SIK3 were also investigated using the kinase assay. Moreover, curcumin exhibited an IC₅₀ (half-maximal inhibitory concentration) value of 131 nM, and it showed significant antiproliferative activities of 9.62 ± 0.33 µM and 72.37 ± 0.37 µM against the MCF-7 and MDA-MB-23 cell lines, respectively. This study provides detailed information on the binding of curcumin with human SIK3 and may facilitate the design of novel salt-inducible kinases inhibitors.