A novel peptide derived from Zingiber cassumunar rhizomes exhibits anticancer activity against the colon adenocarcinoma cells (Caco-2) via the induction of intrinsic apoptosis signaling.

This paper presents the initial exploration of the free radical scavenging capabilities of peptides derived from protein hydrolysates (PPH) obtained from Zingiber cassumunar rhizomes (Phlai). To replicate the conditions of gastrointestinal digestion, a combination of pepsin and pancreatin proteolysis was employed to generate these hydrolysates. Subsequently, the hydrolysate underwent fractionation using molecular weight cut-off membranes at 10, 5, 3, and 0.65 kDa. The fraction with a molecular weight less than 0.65 kDa exhibited the highest levels ABTS, DPPH, FRAP, and NO radical scavenging activity. Following this, RP-HPLC was used to further separate the fraction with a molecular weight less than 0.65 kDa into three sub-fractions. Among these, the F5 sub-fraction displayed the most prominent radical-scavenging properties. De novo peptide sequencing via quadrupole-time-of-flight-electron spin induction-mass spectrometry identified a pair of novel peptides: Asp-Gly-Ile-Phe-Val-Leu-Asn-Tyr (DGIFVLNY or DY-8) and Ile-Pro-Thr-Asp-Glu-Lys (IPTDEK or IK-6). Database analysis confirmed various properties, including biological activity, toxicity, hydrophilicity, solubility, and potential allergy concerns. Furthermore, when tested on the human adenocarcinoma colon (Caco-2) cell line, two synthetic peptides demonstrated cellular antioxidant activity in a concentration-dependent manner. These peptides were also assessed using the FITC Annexin V apoptosis detection kit with PI, confirming the induction of apoptosis. Notably, the DY-8 peptide induced apoptosis, upregulated mRNA levels of caspase-3, -8, and -9, and downregulated Bcl-2, as confirmed by real-time quantitative polymerase chain reaction (RT-qPCR). Western blot analysis indicated increased pro-apoptotic Bax expression and decreased anti-apoptotic Bcl-2 expression in Caco-2 cells exposed to the DY-8 peptide. Molecular docking analysis revealed that the DY-8 peptide exhibited binding affinity with Bcl-2, Bcl-xL, and Mcl-1, suggesting potential utility in combating colon cancer as functional food ingredients.

Naturally occurring quercetin and myricetin as potent inhibitors for human ectonucleotide pyrophosphatase/phosphodiesterase 1.

Ecto-nucleotide pyrophosphatases/phosphodiesterases 1 (ENPP1) is a key enzyme in purinergic signaling pathways responsible for cell-to-cell communications and regulation of several fundamental pathophysiological processes. In this study, Kyoto Green, a rapid chemical sensor of pyrophosphate, was employed to screen for effective ENPP1 inhibitors among five representative flavonoids (quercetin, myricetin, morin, kaempferol, and quercetin-3-glucoside), five nucleosides (adenosine, guanosine, inosine, uridine, and cytidine), and five deoxynucleosides (2'- and 3'-deoxyadenosine, 2'-deoxyguanosine, 2'-deoxyinosine, and 2'-deoxyuridine). Conventional colorimetric, fluorescence, and bioluminescence assays revealed that ENPP1 was effectively inhibited by quercetin (Ki ~ 4 nM) and myricetin (Ki ~ 32 nM) when ATP was used as a substrate at pH 7.4. In silico analysis indicated that the presence of a chromone scaffold, particularly one containing a hydroxyl group at the 3' position on the B ring, may promote binding to the active site pocket of ENPP1 and enhance inhibition. This study demonstrated that the naturally derived quercetin and myricetin could effectively inhibit ENPP1 enzymatic activity and may offer health benefits in arthritis management.

Structural Dynamics of the Precatalytic State of Human Cytochrome c upon T28C, G34C, and A50C Mutations: A Molecular Dynamics Simulation Perspective.

The native structure of cytochrome c (cytc) contains hexacoordinate heme iron with His18 and Met80 residues ligated at the axial sites. Mutations of cytc at Ω-loops have been investigated in modulating the peroxidase activity and, hence, related to the initiation of the apoptotic pathway. Our previous experimental data reported on the peroxidase activity of the cysteine-directed mutants at different parts of the Ω-loop of human cytc (hCytc), that is, T28C, G34C, and A50C. In this work, we performed 1 µs molecular dynamics (MD) simulations to elucidate the detailed structural and dynamic changes upon these mutations, particularly at the proximal Ω-loop. The structures of hCytc were modeled in the hexacoordinated form, which was referred to as the "precatalytic state". The results showed that the structural features of the G34C mutant were more distinctive than those of other mutants. G34C mutation caused local destabilization and flexibility at the proximal Ω-loop (residues 12-28) and an extended distance between this Ω-loop region and heme iron. Besides, analysis of the orientation of the Arg38 side chain of the G34C mutant revealed the Arg38 conformer facing away from the heme iron. The obtained MD results also suggested structural diversity of the precatalytic states for the three hCytc mutants, specifically the effect of G34C mutation on the flexibility of the proximal Ω-loops. Therefore, our MD simulations combined with previous experimental data provide detailed insights into the structural basis of hCytc that could contribute to its pro-apoptotic function.

In Vitro and In Silico Study on the Molecular Encapsulation of α-Tocopherol in a Large-Ring Cyclodextrin.

α-tocopherol is the physiologically most active form of vitamin E, with numerous biological activities, such as significant antioxidant activity, anticancer capabilities, and anti-aging properties. However, its low water solubility has limited its potential use in the food, cosmetic, and pharmaceutical industries. One possible strategy for addressing this issue is the use of a supramolecular complex with large-ring cyclodextrins (LR-CDs). In this study, the phase solubility of the CD26/α-tocopherol complex was investigated to assess the possible ratios between host and guest in the solution phase. Next, the host-guest association of the CD26/α-tocopherol complex at different ratios of 1:2, 1:4, 1:6, 2:1, 4:1, and 6:1 was studied by all-atom molecular dynamics (MD) simulations. At 1:2 ratio, two α-tocopherol units interact spontaneously with CD26, forming an inclusion complex, as supported by the experimental data. In the 2:1 ratio, a single α-tocopherol unit was encapsulated by two CD26 molecules. In comparison, increasing the number of α-tocopherol or CD26 molecules above two led to self-aggregation and consequently limited the solubility of α-tocopherol. The computational and experimental results indicate that a 1:2 ratio could be the most suitable stoichiometry to use in the CD26/α-tocopherol complex to improve α-tocopherol solubility and stability in inclusion complex formation.

Molecular insights on the conformational dynamics of a P76C mutant of human cytochrome c and the enhancement on its peroxidase activity.

In apoptotic pathway, the interaction of Cytochrome c (Cytc) with cardiolipin in vivo is a key process to induce peroxidase activity of Cytc and trigger the release of Cytc in the inner mitochondria into cytosol. The peroxidase active form of Cytc occurs due to local conformational changes that support the opening of the heme crevice and the loss of an axial ligand between Met80 and heme Fe. Structural adjustments at the Ω-loop segments of Cytc are required for such process. To study the role of the distal Ω-loop segments comprising residues 71-85 in human Cytc (hCytc), we investigated a cysteine mutation at Pro76, one of the highly conserved residues in this loop. The effect of P76C mutant was explored by the combination of experimental characterizations and molecular dynamics (MD) simulations. The peroxidase activity of the P76C mutant was found to be significantly increased by ∼13 folds relative to the wild type. Experimental data on global denaturation, alkaline transition, heme bleaching, and spin-labeling Electron Spin Resonance were in good agreement with the enhancement of peroxidase activity. The MD results of hCytc in the hexacoordinate form suggest the important changes in P76C mutant occurred due to the unfolding at the central Ω-loop (residues 40-57), and the weakening of H-bond between Tyr67 and Met80. Whereas the experimental data implied that the P76C mutant tend to be in equilibrium between the pentacoordinate and hexacoordinate forms, the MD and experimental information are complementary and were used to support the mechanisms of peroxidase active form of hCytc.

Titania Nanosheet Generates Peroxynitrite-Dependent S-Nitrosylation and Enhances p53 Function in Lung Cancer Cells.

Metal nanomaterials can enhance the efficacy of current cancer therapies. Here, we show that Ti0.8O2 nanosheets cause cytotoxicity in several lung cancer cells but not in normal cells. The nanosheet-treated cells showed certain apoptosis characteristics. Protein analysis further indicated the activation of the p53-dependent death mechanism. Transmission electron microscopy (TEM) and scanning electron microscopy (SEM) analyses revealed the cellular uptake of the nanosheets and the induction of cell morphological change. The nanosheets also exhibited a substantial apoptosis effect on drug-resistant metastatic primary lung cancer cells, and it was found that the potency of the nanosheets was dramatically higher than standard drugs. Ti0.8O2 nanosheets induce apoptosis through a molecular mechanism involving peroxynitrite (ONOO-) generation. As peroxynitrite is known to be a potent inducer of S-nitrosylation, we further found that the nanosheets mediated the S-nitrosylation of p53 at C182, resulting in higher protein-protein complex stability, and this was likely to induce the surrounding residues, located in the interface region, to bind more strongly to each other. Molecular dynamics analysis revealed that S-nitrosylation stabilized the p53 dimer with a ΔGbindresidue of <-1.5 kcal/mol. These results provide novel insight on the apoptosis induction effect of the nanosheets via a molecular mechanism involving S-nitrosylation of the p53 protein, emphasizing the mechanism of action of nanomaterials for cancer therapy.

Bibenzyl analogue DS-1 inhibits MDM2-mediated p53 degradation and sensitizes apoptosis in lung cancer cells.

BACKGROUND: Lung cancer is a leading fatal malignancy due to the high incidence of treatment failure. Dysfunction of the tumor suppressor p53 contributes to cancer initiation, progression, and therapeutic resistance. Targeting MDM2, a negative regulator of p53, has recently attracted interest in cancer drug research as it may restore tumor suppressive function. PURPOSE: The present study aimed to investigate the effect of 3,4-dihydroxy-5,4'-dimethoxybibenzyl (DS-1) on targeting MDM2 and restoring p53 function in lung cancer cells. METHODS: The efficacy of DS-1 alone or in combination with cisplatin in lung cancer cells was determined by MTT, nuclear staining, and annexin V/PI assay. The expression of apoptosis-related proteins was determined by western blot analysis. To evaluate the role of DS-1 on the stabilization and degradation of p53, cycloheximide chasing assay and immunoprecipitation were conducted, and the active form of p53 was investigated by immunofluorescent staining assay. To confirm and demonstrate the site interaction between DS-1 and the MDM2 protein, in silico computational analysis was performed. RESULTS: DS-1 exhibited a cytotoxic effect and sensitized lung cancer cells to cisplatin-induced apoptosis. DS-1 caused a significant increase in the cellular level of p53 protein, while the active form of p53 (phosphorylation at Ser15) was unaltered. DS-1 treatment in combination with cisplatin could enhance activated p-p53 (Ser15) and p53 downstream signaling (Bax, Bcl-2, and Akt), leading to a higher level of apoptosis. Immunoprecipitation analysis revealed that DS-1 decreased the p53-ubiquitin complex, a prerequisite step in p53 proteasomal degradation. Molecular docking simulation further evidenced that DS-1 interacts with MDM2 within the p53-binding domain by carbon-hydrogen bond interaction at Lys27, π-alkyl interactions at Ile37 and Leu30, and van der Waals interactions at Ile75, Val51, Val69, Phe67, Met38, Tyr43, Gly34, and Phe31. Treatment by DS-1 and cisplatin in patient-derivated primary lung cancer cells showed consistent effects by increasing cisplatin sensitivity. CONCLUSIONS: Our findings provide evidence that DS-1 is an MDM2 inhibitor and its underlying mechanism involves MDM2 binding and p53 induction, which may benefit the development of this compound for lung cancer treatment.

Delivery of Alpha-Mangostin Using Cyclodextrins through a Biological Membrane: Molecular Dynamics Simulation.

α-Mangostin (MGS) exhibits various pharmacological activities, including antioxidant, anticancer, antibacterial, and anti-inflammatory properties. However, its low water solubility is the major obstacle for its use in pharmaceutical applications. To increase the water solubility of MGS, complex formation with beta-cyclodextrins (ßCDs), particularly with the native ßCD and/or its derivative 2,6-dimethyl-ß-CD (DMßCD) is a promising technique. Although there have been several reports on the adsorption of ßCDs on the lipid bilayer, the release of the MGS/ßCDs inclusion complex through the biological membrane remains unclear. In this present study, the release the MGS from the two different ßCDs (ßCD and DMßCD) across the lipid bilayer was investigated. Firstly, the adsorption of the free MGS, free ßCDs, and inclusion complex formation was studied by conventional molecular dynamics simulation. The MGS in complex with those two ßCDs was able to spontaneously release free MGS into the inner membrane. However, both MGS and DMßCD molecules potentially permeated into the deeper region of the interior membrane, whereas ßCD only adsorbed at the outer membrane surface. The interaction between secondary rim of ßCD and the 1-palmitoeyl-2-oleoyl-glycero-3-phosphocholine (POPC) phosphate groups showed the highest number of hydrogen bonds (up to 14) corresponding to the favorable location of ßCD on the POPC membrane. Additionally, the findings suggested that electrostatic energy was the main driving force for ßCD adsorption on the POPC membrane, while van der Waals interactions played a predominant role in DMßCD adsorption. The release profile of MGS from the ßCDs pocket across the lipid bilayer exhibited two energy minima along the reaction coordinate associated with the permeation of the MGS molecule into the deeper region of the POPC membrane.

Why Are Lopinavir and Ritonavir Effective against the Newly Emerged Coronavirus 2019? Atomistic Insights into the Inhibitory Mechanisms.

Since the emergence of a novel coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) was first reported from Wuhan, China, neither a specific vaccine nor an antiviral drug against SARS-CoV-2 has become available. However, a combination of two HIV-1 protease inhibitors, lopinavir and ritonavir, has been found to be effective against SARS-CoV, and both drugs could bind well to the SARS-CoV 3C-like protease (SARS-CoV 3CLpro). In this work, molecular complexation between each inhibitor and SARS-CoV-2 3CLpro was studied using all-atom molecular dynamics simulations, free energy calculations, and pair interaction energy analyses based on MM/PB(GB)SA and FMO-MP2/PCM/6-31G* methods. Both anti-HIV drugs interacted well with the residues at the active site of SARS-CoV-2 3CLpro. Ritonavir showed a somewhat higher number atomic contacts, a somewhat higher binding efficiency, and a somewhat higher number of key binding residues compared to lopinavir, which correspond with the slightly lower water accessibility at the 3CLpro active site. In addition, only ritonavir could interact with the oxyanion hole residues N142 and G143 via the formation of two hydrogen bonds. The interactions in terms of electrostatics, dispersion, and charge transfer played an important role in the drug binding. The obtained results demonstrated how repurposed anti-HIV drugs could be used to combat COVID-19.

Structure-Activity Relationships and Molecular Docking Analysis of Mcl-1 Targeting Renieramycin T Analogues in Patient-derived Lung Cancer Cells.

Myeloid cell leukemia 1 (Mcl-1) and B-cell lymphoma 2 (Bcl-2) proteins are promising targets for cancer therapy. Here, we investigated the structure-activity relationships (SARs) and performed molecular docking analysis of renieramycin T (RT) and its analogues and identified the critical functional groups of Mcl-1 targeting. RT have a potent anti-cancer activity against several lung cancer cells and drug-resistant primary cancer cells. RT mediated apoptosis through Mcl-1 suppression and it also reduced the level of Bcl-2 in primary cells. For SAR study, five analogues of RT were synthesized and tested for their anti-cancer and Mcl-1- and Bcl-2-targeting effects. Only two of them (TM-(-)-18 and TM-(-)-4a) exerted anti-cancer activities with the loss of Mcl-1 and partly reduced Bcl-2, while the other analogues had no such effects. Specific cyanide and benzene ring parts of RT's structure were identified to be critical for its Mcl-1-targeting activity. Computational molecular docking indicated that RT, TM-(-)-18, and TM-(-)-4a bound to Mcl-1 with high affinity, whereas TM-(-)-45, a compound with a benzene ring but no cyanide for comparison, showed the lowest binding affinity. As Mcl-1 helps cancer cells evading apoptosis, these data encourage further development of RT compounds as well as the design of novel drugs for treating Mcl-1-driven cancers.

Binding recognition of substrates in NS2B/NS3 serine protease of Zika virus revealed by molecular dynamics simulations.

Zika virus (ZIKV) has become a global public health concern. The recent epidemiological data has revealed a possible association of ZIKV infection with more serious complications, particularly for Guillain-Barré syndrome in adults and microcephaly in newborn children. Till now, there is no vaccine or effective drug commercially available to combat with ZIKV infection. An attractive drug target for the ZIKV treatment is the NS2B/NS3 serine protease, which is essential for viral polyprotein processing. Herein, classical molecular dynamics (MD) simulations were performed on the ZIKV NS2B/NS3 serine protease in complex with four peptide substrates to investigate the binding recognition and protein-substrate interactions. The obtained results indicate that the P1 and P2 positions of the substrate play a significant role in binding with the protease enzyme, while the P3 and P4 positions show a minor contribution in binding interaction. Moreover, the binding free energy calculation based on the MM/PBSA method suggests that among the four similar peptide substrates, the peptide Ac-D-RKOR-ACC displays the strongest binding affinity towards the ZIKV protease due to the high energy contribution at the S2 subsite particularly for the NS3 residue D75 with the P2(O) residue of this substrate, which is in line with the experimental data. Thus, the information derived from MD simulations presented here would be useful for the design of potent protease inhibitors.

Quantum Mechanics/Molecular Mechanics (QM/MM) Calculations Support a Concerted Reaction Mechanism for the Zika Virus NS2B/NS3 Serine Protease with Its Substrate.

Zika virus (ZIKV) is mainly transmitted to humans by Aedes species mosquitoes and is associated with serious pathological disorders including microcephaly in newborns and Guillain-Barré syndrome in adults. Currently, there is no vaccine or anti-ZIKV drug available for preventing or controlling ZIKV infection. An attractive drug target for ZIKV treatment is a two-compartment (NS2B/NS3) serine protease that processes viral polyprotein during infection. Here, conventional molecular dynamics simulations of the ZIKV protease in complex with peptide substrate (TGKRS) sequence at the C-terminus of NS2B show that the substrate is in the active conformation for the cleavage reaction by ZIKV protease. Hybrid quantum mechanics/molecular mechanics (QM/MM) umbrella sampling simulations (PM6/ff14SB) of acylation results reveal that proton transfer from S135 to H51 and nucleophilic attack on the substrate by S135 are concerted. The rate-limiting step involves the formation of a tetrahedral intermediate. In addition, the single-point energy QM/MM calculations, precisely at the level of coupled cluster theory (LCCSD(T)/(aug)-cc-pVTZ), were performed to correct the potential energy profiles for the first step of the acylation process. The average computed activation barrier at this level of theory is 16.3 kcal mol-1. Therefore, the computational approaches presented here are helpful for further designing of NS2B/NS3 inhibitors based on transition-state analogues.

Low susceptibility of asunaprevir towards R155K and D168A point mutations in HCV NS3/4A protease: A molecular dynamics simulation.

Hepatitis C has become an important health problem that requires expensive treatment and leads to liver tumorigenesis. Hepatitis C virus (HCV), which is the main cause of hepatitis C, has a high mutation rate due to the lack of proofreading activity of the RNA polymerase enzyme. The NS3/4A serine protease is an important target for anti-HCV drug discovery and development because of its crucial role in the cleavage of the polypeptides involved in viral replication. In the present study, all-atom molecular dynamics simulation was performed to elucidate the effect of the single point mutations R155K and D168A in the HCV genotype 1 NS3/4A protease on the structural dynamics, molecular interactions and susceptibility of asunaprevir (ASV), a second-generation NS3/4A protease inhibitor. Principal component analysis indicated that these two mutations converted the direction of motion of residues 123, 155 and 168 in the binding pocket to significantly point outwards from ASV, resulting in a loss of the hydrogen bond network of residues R123···R155···D168. The free energy calculations based on different semiempirical QM/MM-GBSA methods revealed that the binding affinity of ASV with the two mutant forms of the NS3/4A protease was significantly decreased in the order of wild-type < R155K < D168A. This work provided useful structural information regarding the atomistic understanding of acquired drug resistance against ASV caused by the R155K and D168A mutations.

Computational screening of fatty acid synthase inhibitors against thioesterase domain.

Thioesterase (TE) domain of fatty acid synthase (FAS) is an attractive therapeutic target for design and development of anticancer drugs. In this present work, we search for the potential FAS inhibitors of TE domain from the ZINC database based on similarity search using three natural compounds as templates, including flavonoids, terpenoids, and phenylpropanoids. Molecular docking was used to predict the interaction energy of each screened ligand compared to the reference compound, which is methyl γ-linolenylfluorophosphonate (MGLFP). Based on this computational technique, rosmarinic acid and its eight analogs were observed as a new series of potential FAS inhibitors, which showed a stronger binding affinity than MGLFP. Afterward, nine docked complexes were studied by molecular dynamics simulations for investigating protein-ligand interactions and binding free energies using MM-PB(GB)SA, MM-3DRISM-KH, and QM/MM-GBSA methods. The binding free energy calculation indicated that the ZINC85948835 (R34) displayed the strongest binding efficiency against the TE domain of FAS. There are eight residues (S2308, I2250, E2251, Y2347, Y2351, F2370, L2427, and E2431) mainly contributed for the R34 binding. Moreover, R34 could directly form hydrogen bonds with S2308, which is one of the catalytic triad of TE domain. Therefore, our finding suggested that R34 could be a potential candidate as a novel FAS-TE inhibitor for further drug design.

Screening of hepatitis C NS5B polymerase inhibitors containing benzothiadiazine core: a steered molecular dynamics.

Hepatic C virus (HCV) is a global health problem, resulting in liver cirrhosis and inflammation that can develop to hepatocellular carcinoma and fatality. The NS5B polymerase of HCV plays an important role in viral RNA replication process, making it an attractive therapeutic target for design and development of anti-HCV drugs. To search new potent compounds against the HCV NS5B polymerase, the molecular docking and the steered molecular dynamics (SMD) simulation techniques were performed. The potential potent inhibitors of the NS5B polymerase were screened out from the ZINC database using structural similarity search and molecular docking technique. Five top-hit compounds (the ZINC compounds 49888724, 49054741, 49777239, 49793673, and 49780355) were then studied by the SMD simulations based on the hypothesis that a high rupture force relates to a high binding efficiency. The results demonstrated that the ZINC compound 49888724 had a greater maximum rupture force, reflecting a good binding strength and inhibitory potency than known inhibitors and the rest four ZINC compounds. Therefore, our finding indicated that the ZINC compound 49888724 is a potential candidate to be a novel NS5B inhibitor for further design. Besides, the van der Waals interaction could be considered as the main contribution for stabilizing the NS5B-ligand complex.