Senescent cells accumulate in tissues over time as part of the natural ageing process and the removal of senescent cells has shown promise for alleviating many different age-related diseases in mice. Cancer is an age-associated disease and there are numerous mechanisms driving cellular senescence in cancer that can be detrimental to recovery. Thus, it would be beneficial to develop a senolytic that acts not only on ageing cells but also senescent cancer cells to prevent cancer recurrence or progression. We used molecular modelling to develop a series of rationally designed peptides to mimic and target FOXO4 disrupting the FOXO4-TP53 interaction and releasing TP53 to induce apoptosis. We then tested these peptides as senolytic agents for the elimination of senescent cells both in cell culture and in vivo. Here we show that these peptides can act as senolytics for eliminating senescent human cancer cells both in cell culture and in orthotopic mouse models. We then further characterized one peptide, ES2, showing that it disrupts FOXO4-TP53 foci, activates TP53 mediated apoptosis and preferentially binds FOXO4 compared to TP53. Next, we show that intratumoural delivery of ES2 plus a BRAF inhibitor results in a significant increase in apoptosis and a survival advantage in mouse models of melanoma. Finally, we show that repeated systemic delivery of ES2 to older mice results in reduced senescent cell numbers in the liver with minimal toxicity. Taken together, our results reveal that peptides can be generated to specifically target and eliminate FOXO4+ senescent cancer cells, which has implications for eradicating residual disease and as a combination therapy for frontline treatment of cancer. This work was supported by the Cancer Early Detection Advanced Research Center at Oregon Health & Science University.
The docking program PLANTS, which is based on ant colony optimization (ACO) algorithm, has many advanced features for molecular docking. Among them are multiple scoring functions, the possibility to model explicit displaceable water molecules, and the inclusion of experimental constraints. Here, we add support of PLANTS to VirtualFlow (VirtualFlow Ants), which adds a valuable method for primary virtual screenings and rescoring procedures. Furthermore, we have added support of ligand libraries in the MOL2 format, as well as on the fly conversion of ligand libraries which are in the PDBQT format to the MOL2 format to endow VirtualFlow Ants with an increased flexibility regarding the ligand libraries. The on the fly conversion is carried out with Open Babel and the program SPORES. We applied VirtualFlow Ants to a test system involving KEAP1 on the Google Cloud up to 128,000 CPUs, and the observed scaling behavior is approximately linear. Furthermore, we have adjusted several central docking parameters of PLANTS (such as the speed parameter or the number of ants) and screened 10 million compounds for each of the 10 resulting docking scenarios. We analyzed their docking scores and average docking times, which are key factors in virtual screenings. The possibility of carrying out ultra-large virtual screening with PLANTS via VirtualFlow Ants opens new avenues in computational drug discovery.
MM-PB/GBSA methods represent a higher-level scoring theory than docking. This study reports an extensive testing of different MM-GBSA scoring schemes on two bromodomain (BRD) datasets. The first set is composed of 24 BRPF1 complexes, and the second one is a nonredundant set constructed from the PDBbind and composed of 28 diverse BRD complexes. A variety of MM-GBSA schemes were analyzed to evaluate the performance of four protocols with different numbers of minimization and MD steps, 10 different force fields and three different water models. Results showed that neither additional MD steps nor unfixing the receptor atoms improved scoring or ranking power. On the contrary, our results underscore the advantage of fixing receptor atoms or limiting the number of MD steps not only for a reduction in the computational costs but also for boosting the prediction accuracy. Among Amber force fields tested, ff14SB and its derivatives rather than ff94 or polarized force fields provided the most accurate scoring and ranking results. The TIP3P water model yielded the highest scoring and ranking power compared to the others. Posing power was further evaluated for the BRPF1 set. A slightly better posing power for the protocol which uses both minimization and MD steps with a fixed receptor than the one which uses only minimization with a fully flexible receptor-ligand system was observed. Overall, this study provides insights into the usage of the MM-GBSA methods for screening of BRD inhibitors, substantiating the benefits of shorter protocols and latest force fields and maintaining the crystal waters for accuracy.
This study aims to synthesise and characterise novel compounds containing 2-amino-1,3,4-thiadiazole and their acyl derivatives and to investigate antifungal activities. Similarity search, molecular dynamics and molecular docking were also studied to find out a potential target and enlighten the inhibition mechanism. As a first step, 2-amino-1,3,4-thiadiazole derivatives (compounds 3 and 4) were synthesised with high yields (81 and 84%). The target compounds (6a–n and 7a–n) were then synthesised with moderate to high yields (56–87%) by reacting 3 and 4 with various acyl chloride derivatives (5a–n). The synthesized compounds were characterized using the IR, 1H-NMR, 13C-NMR, Mass, X-ray (compound 7n) and elemental analysis techniques. Later, the in vitro antifungal activities of the synthesised compounds were determined. The inhibition zones exhibited by the compounds against the tested fungi, their minimum fungicidal activities, minimum inhibitory concentration and the lethal dose values (LD50) were determined. The compounds exhibited moderate to high levels of activity against all tested pathogens. Finally, in silico modelling was used to enlighten inhibition mechanism using ligand and structure-based methods. As an initial step, similarity search was carried out and the resulting proteins that belong to Homo sapiens were used as reference in sequence similarity search to find the corresponding amino acid sequences in target organisms. Homology modelling was used to construct the protein structure. The stabilised protein structure obtained from molecular dynamics simulation was used in molecular docking. The overall results presented here might be a good starting point for the identification of novel and more active compounds as antifungal agents.
4 February 2021 in The Journal of Physical Chemistry B