Influence of Chirality of Crizotinib on Its MTH1 Protein Inhibitory Activity: Insight from Molecular Dynamics Simulations and Binding Free Energy Calculations

Overcoming Chemoresistance in Cancer: The Promise of Crizotinib

Chemoresistance is a major obstacle in cancer treatment, often leading to disease progression and poor outcomes. It arises through various mechanisms such as genetic mutations, drug efflux pumps, enhanced DNA repair, and changes in the tumor microenvironment. These processes allow cancer cells to survive despite chemotherapy, underscoring the need for new strategies to overcome resistance and improve treatment efficacy. Crizotinib, a first-generation multi-target kinase inhibitor, is approved by the FDA for the treatment of ALK-positive or ROS1-positive non-small cell lung cancer (NSCLC), refractory inflammatory (ALK)-positive myofibroblastic tumors (IMTs) and relapsed/refractory ALK-positive anaplastic large cell lymphoma (ALCL). Crizotinib exists in two enantiomeric forms: (R)-crizotinib and its mirror image, (S)-crizotinib. It is assumed that the R-isomer is responsible for the carrying out various processes reviewed here The S-isomer, on the other hand, shows a strong inhibition of MTH1, an enzyme important for DNA repair mechanisms. Studies have shown that crizotinib is an effective multi-kinase inhibitor targeting various kinases such as c-Met, native/T315I Bcr/Abl, and JAK2. Its mechanism of action involves the competitive inhibition of ATP binding and allosteric inhibition, particularly at Bcr/Abl. Crizotinib showed synergistic effects when combined with the poly ADP ribose polymerase inhibitor (PARP), especially in ovarian cancer harboring BRCA gene mutations. In addition, crizotinib targets a critical vulnerability in many p53-mutated cancers. Unlike its wild-type counterpart, the p53 mutant promotes cancer cell survival. Crizotinib can cause the degradation of the p53 mutant, sensitizing these cancer cells to DNA-damaging substances and triggering apoptosis. Interestingly, other reports demonstrated that crizotinib exhibits anti-bacterial activity, targeting Gram-positive bacteria. Also, it is active against drug-resistant strains. In summary, crizotinib exerts anti-tumor effects through several mechanisms, including the inhibition of kinases and the restoration of drug sensitivity. The potential of crizotinib in combination therapies is emphasized, particularly in cancers with a high prevalence of the p53 mutant, such as triple-negative breast cancer (TNBC) and high-grade serous ovarian cancer (HGSOC).

Exploring MTH1 inhibitory potential of Thymoquinone and Baicalin for therapeutic targeting of breast cancer

Cancers frequently have increased ROS levels due to disrupted redox balance, leading to oxidative DNA and protein damage, mutations, and apoptosis. The MTH1 protein plays a crucial role by sanitizing the oxidized dNTP pools. Hence, cancer cells rely on MTH1 to prevent the integration of oxidized dNTPs into DNA, preventing DNA damage and allowing cancer cell proliferation. We have discovered Thymoquinone (TQ) and Baicalin (BC) as inhibitors of MTH1 using combined docking and MD simulation approaches complemented by experimental validations via assessing binding affinity and enzyme inhibition. Docking and MD simulations studies revealed an efficient binding of TQ and BC to the active site pocket of the MTH1, and the resultant complexes are appreciably stable. Fluorescence measurements estimated a strong binding affinity of TQ and BC with Ka 3.4 ×106 and 1.0 ×105, respectively. Treating breast cancer cells with TQ and BC significantly inhibited the growth and proliferation (IC50 values 28.3 µM and 34.8 µM) and induced apoptosis. TQ and BC increased the ROS production in MCF7 cells, imposing substantial oxidative stress on cancer cells and leading to cell death. Finally, TQ and BC are proven strong MTH1 inhibitors, offering promising prospects for anti-cancer therapy.

Ligand-Based Virtual Screening and Molecular Docking of Benzimidazoles as Potential Inhibitors of Triosephosphate Isomerase Identified New Trypanocidal Agents

Trypanosoma cruzi (T. cruzi) is a parasite that affects humans and other mammals. T. cruzi depends on glycolysis as a source of adenosine triphosphate (ATP) supply, and triosephosphate isomerase (TIM) plays a key role in this metabolic pathway. This enzyme is an attractive target for the design of new trypanocidal drugs. In this study, a ligand-based virtual screening (LBVS) from the ZINC15 database using benzimidazole as a scaffold was accomplished. Later, a molecular docking on the interface of T. cruzi TIM (TcTIM) was performed and the compounds were grouped by interaction profiles. Subsequently, a selection of compounds was made based on cost and availability for in vitro evaluation against blood trypomastigotes. Finally, the compounds were analyzed by molecular dynamics simulation, and physicochemical and pharmacokinetic properties were determined using SwissADME software. A total of 1604 molecules were obtained as potential TcTIM inhibitors. BP2 and BP5 showed trypanocidal activity with half-maximal lytic concentration (LC50) values of 155.86 and 226.30 µM, respectively. Molecular docking and molecular dynamics simulation analyzes showed a favorable docking score of BP5 compound on TcTIM. Additionally, BP5 showed a low docking score (-5.9 Kcal/mol) on human TIM compared to the control ligand (-7.2 Kcal/mol). Both compounds BP2 and BP5 showed good physicochemical and pharmacokinetic properties as new anti-T. cruzi agents.

Identification of a novel immune signature for optimizing prognosis and treatment prediction in colorectal cancer

BACKGROUND

Globally, colorectal cancer (CRC) is one of the most lethal malignant diseases. However, the currently approved therapeutic options for CRC failed to acquire satisfactory treatment efficacy. Tailoring therapeutic strategies for CRC individuals can provide new insights into personalized prediction approaches and thus maximize clinical benefits.

METHODS

In this study, a multi-step process was used to construct an immune-related genes (IRGs) based signature leveraging the expression profiles and clinical characteristics of CRC from the Gene Expression Omnibus (GEO) database and the Cancer Genome Atlas (TCGA) database. An integrated immunogenomic analysis was performed to determine the association between IRGs with prognostic significance and cancer genotypes in the tumor immune microenvironment (TIME). Moreover, we performed a comprehensive in silico therapeutics screening to identify agents with subclass-specific efficacy.

RESULTS

The established signature was shown to be a promising biomarker for evaluating clinical outcomes in CRC. The immune risk score as calculated by this classifier was significantly correlated with over-riding malignant phenotypes and immunophenotypes. Further analyses demonstrated that CRCs with low immune risk scores achieved better therapeutic benefits from immunotherapy, while AZD4547, Cytochalasin B and S-crizotinib might have potential therapeutic implications in the immune risk score-high CRCs.

CONCLUSIONS

Overall, this IRGs-based signature not only afforded a useful tool for determining the prognosis and evaluating the TIME features of CRCs, but also shed new light on tailoring CRCs with precise treatment.

Inhibitors of DNA Glycosylases as Prospective Drugs

DNA glycosylases are enzymes that initiate the base excision repair pathway, a major biochemical process that protects the genomes of all living organisms from intrinsically and environmentally inflicted damage. Recently, base excision repair inhibition proved to be a viable strategy for the therapy of tumors that have lost alternative repair pathways, such as BRCA-deficient cancers sensitive to poly(ADP-ribose)polymerase inhibition. However, drugs targeting DNA glycosylases are still in development and so far have not advanced to clinical trials. In this review, we cover the attempts to validate DNA glycosylases as suitable targets for inhibition in the pharmacological treatment of cancer, neurodegenerative diseases, chronic inflammation, bacterial and viral infections. We discuss the glycosylase inhibitors described so far and survey the advances in the assays for DNA glycosylase reactions that may be used to screen pharmacological libraries for new active compounds.

Differential anti-tumour effects of MTH1 inhibitors in patient-derived 3D colorectal cancer cultures

Reactive oxygen species (ROS) function as second messengers in signal transduction, but high ROS levels can also cause cell death. MTH1 dephosphorylates oxidized nucleotides, thereby preventing their incorporation into DNA and protecting tumour cells from oxidative DNA damage. Inhibitors of MTH1 (TH588 and (S)-crizotinib) were shown to reduce cancer cell viability. However, the MTH1-dependency of the anti-cancer effects of these drugs has recently been questioned. Here, we have assessed anti-tumour effects of TH588 and (S)-crizotinib in patient-derived 3D colorectal cancer cultures. Hypoxia and reoxygenation – conditions that increase intracellular ROS levels – increased sensitivity to (S)-crizotinib, but not to TH588. (S)-crizotinib reduced tyrosine phosphorylation of c-MET and ErbB3 whereas TH588 induced a mitotic cell cycle arrest, which was not affected by adding ROS-modulating compounds. Furthermore, we show that both compounds induced DNA damage that could not be prevented by adding the ROS inhibitor N-acetyl-L-cysteine. Moreover, adding ROS-modulating compounds did not alter the reduction in viability in response to TH588 and (S)-crizotinib. We conclude that TH588 and (S)-crizotinib have very clear and distinct anti-tumour effects in 3D colorectal cancer cultures, but that these effects most likely occur through distinct and ROS-independent mechanisms.

The presumed MTH1-inhibitor TH588 sensitizes colorectal carcinoma cells to ionizing radiation in hypoxia

Background

The nudix family member enzyme MutT homologue-1 (MTH1) hydrolyses the oxidized nucleotides 8-oxo-dGTP and 2-hydroxy-dATP and thus prevents the incorporation of damaged nucleotides into nuclear and mitochondrial DNA. Therefore MTH1 was proposed to protect cancer cells from oxidative DNA lesions and subsequent cell death. We investigated whether the bona fide MTH1 inhibitor TH588 affects responses of cultured colorectal tumor cells to ionizing radiation (IR) in normoxia and in moderate or severe hypoxia.

Methods

TH588 was tested in cell viability and survival assays (tetrazolium dye (MTT), propidium iodide staining, caspase-3 activity, and colony formation assays (CFA)) in colorectal carcinoma cells (HCT116 and SW480) in combination with IR in normoxia and in hypoxia. Additionally, MTH1 was targeted by lentiviral shRNA expression. Human umbilical vein endothelial cells (HUVEC) were assessed in MTT assays.

Results

In all cell lines tested, TH588 dose-dependently impaired cell survival. In CFAs, TH588 and IR effects on carcinoma cells were additive in normoxia and in hypoxia. Using 3 different shRNAs, the lentiviral approach was detrimental to SW480, but not to HCT116.

Conclusions

TH588 has cytotoxic effects on transformed and untransformed cells and synergizes with IR in normoxia and in hypoxia. TH588 toxicity is not fully explained by MTH1 inhibition as HCT116 were unaffected by lentiviral suppression of MTH1 expression. TH588 should be explored further because it has radiosensitizing effects in hypoxia.

Molecular dynamics provides insight into how N251A and N251Y mutations in the active site of Bacillus licheniformis RN-01 levansucrase disrupt production of long-chain levan

Produced by levansucrase, levan and levan oligosaccharides (GF_n) have potential applications in food and pharmaceutical industries such as prebiotics, anti-tumor and anti-inflammatory agents. Previous study reported that Bacillus licheniformis RN-01 levansucrase could produce levan oligosaccharides and long-chain levan. However, its N251A and N251Y mutants could effectively produce short-chain oligosaccharides upto GF_3, but they could not produce long-chain levan. We hypothesized that these mutations probably reduced GF₃ binding affinity in levansucrase active site that contains fructosyl-Asp93 intermediate and caused GF₃ to be in an unfavorable orientation for transfructosylation; therefore, levansucrase could not effectively extend GF₃ by one fructosyl residue to produce GF₄ and subsequently long-chain levan. However, these mutations probably did not significantly reduce binding affinity or drastically change orientation of GF₂; therefore, levansucrase could still extend GF₂ to produce GF₃. Using this hypothesis, we employed molecular dynamics to investigate effects of these mutations on GF₂/GF₃ binding in levansucrase active site. Our results reasonably support this hypothesis as N251A and N251Y mutations did not significantly reduce GF₂ binding affinity, as calculated by MM-GBSA technique and hydrogen bond occupations, or drastically change orientation of GF₂ in levansucrase active site, as measured by distance between atoms necessary for transfructosylation. However, these mutations drastically decreased GF₃ binding affinity and caused GF₃ to be in an unfavorable orientation for transfructosylation. Furthermore, the free energy decomposition and hydrogen bond occupation results suggest the importance of Arg255 in GF₂/GF₃ binding in levansucrase active site. This study provides important and novel insight into the effects of N251A and N251Y mutations on GF₂/GF₃ binding in levansucrase active site and how they may disrupt production of long-chain levan. This knowledge could be beneficial in designing levansucrase to efficiently produce levan oligosaccharides with desired length.

Resistance mechanism of the oncogenic β3-αC deletion mutation in BRAF kinase to dabrafenib and vemurafenib revealed by molecular dynamics simulations and binding free energy calculations

BRAF kinase is an essential target for anti-cancer drug development. Emergence of the β3-αC loop deletion mutation (ΔNVTAP) in BRAF kinase frequently occurred in human cancers seriously compromises the therapeutic efficacy of some BRAF kinase inhibitors, such as dabrafenib and vemurafenib. However, the mechanism of this resistance is still not well understood. In this study, the influence of the β3-αC deletion mutation on the binding profiles of three BRAF kinase inhibitors (AZ628, dabrafenib, and vemurafenib) with BRAF^V600E or BRAF^ΔNVTAP was explored by conventional molecular dynamics (MD) simulations and binding free energy calculations. The simulation results indicated that the β3-αC deletion mutation enhances the flexibility of the αC helix and alters their conformations, which amplify the conformational entropy change (-TΔS) and weaken the interactions between the inhibitors and BRAF. The further per-residue binding free energy decomposition analysis revealed that the ΔNVTAP mutation changed the contributions of a few key residues to the bindings of dabrafenib or vemurafenib, such as L57, L66, W83, C84, F135, G145, and F147, but did not have obvious impact on the contributions of these residues to AZ628. Our results provide valuable clues to understand the mechanisms of drug resistance conferred by the β3-αC deletion mutation.

The influence of residue in the position of 116 on the inhibitory potency of TH588 for MTH1

As one of the first-in-class inhibitor, TH588 was found to be efficient in the suppression of MutT homolog1 (MTH1). A recent work shows that the inhibitory potency of TH588 against human MTH1 (hsMTH1) is approximately 20-fold over that of mouse MTH1 (mmMTH1) and identifies residue in position 116 in MTH1 has an important contribution to TH588 affinity. But the effect of residue Leu or Met in position 116 on the binding affinity remains unclear. In this study, molecular dynamics (MD) simulations and free energy calculations were used to elucidate the mechanism about the effect of residue 116 to the different inhibitory potency of TH588 against MTH1. The binding free energy of TH588 in M116 complexes predicated by the Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) is much lower than that in L116 complexes, which is consistent with the experiment results. The analysis of the individual energy terms suggests that the non-polar interactions are important for distinguishing the binding of TH588. The MD results show that the Leu116 disrupts the interactions between Asn33 and TH588, thus induces the conformational changes of Asn33 as well as TH588. The altered interactions between TH588 and mmMTH1 change the flexibility of TH588, which could induce the remarkable conformational fluctuation of mmMTH1. The conformations of the two loops covering the binding pocket have obvious influence on the opening or closure of the active site. The more open binding site may explain the lower inhibitor potency of TH588 against mmMTH1 than hsMTH1. Our results provide mechanistic insight into the effect of different residue Leu or Met in position 116 on the binding affinity of TH588 for MTH1, which is expected to contribute to the further rational design of more potent inhibitors.

Insight into the substrate specificity change caused by the Y227H mutation of α-glucosidase III from the European honeybee (Apis mellifera) through molecular dynamics simulations

Honey from the European honeybee, Apis mellifera, is produced by α-glucosidases (HBGases) and is widely used in food, pharmaceutical, and cosmetic industries. Categorized by their substrate specificities, HBGases have three isoforms: HBGase I, II and III. Previous experimental investigations showed that wild-type HBGase III from Apis mellifera (WT) preferred sucrose to maltose as a substrate, while the Y227H mutant (MT) preferred maltose to sucrose. This mutant can potentially be used for malt hydrolysis because it can efficiently hydrolyze maltose. In this work, to elucidate important factors contributing to substrate specificity of this enzyme and gain insight into how the Y227H mutation causes substrate specificity change, WT and MT homology models were constructed, and sucrose/maltose was docked into active sites of the WT and MT. AMBER14 was employed to perform three independent molecular dynamics runs for these four complexes. Based on the relative binding free energies calculated by the MM-GBSA method, sucrose is better than maltose for WT binding, while maltose is better than sucrose for MT binding. These rankings support the experimentally observed substrate specificity that WT preferred sucrose to maltose as a substrate, while MT preferred maltose to sucrose, suggesting the importance of binding affinity for substrate specificity. We also found that the Y227H mutation caused changes in the proximities between the atoms necessary for sucrose/maltose hydrolysis that may affect enzyme efficiency in the hydrolysis of sucrose/maltose. Moreover, the per-residue binding free energy decomposition results show that Y227/H227 may be a key residue for preference binding of sucrose/maltose in the WT/MT active site. Our study provides important and novel insight into the binding of sucrose/maltose in the active site of Apis mellifera HBGase III and into how the Y227H mutation leads to the substrate specificity change at the molecular level. This knowledge could be beneficial in the design of this enzyme for increased production of desired products.

(S)-crizotinib reduces gastric cancer growth through oxidative DNA damage and triggers pro-survival akt signal

Gastric cancer (GC), a common gastrointestinal malignancy worldwide, has poor prognosis and frequent recurrence. There is a great need to identify effective therapy for GC. Crizotinib is a multi-targeted, clinically available oral tyrosine kinase inhibitor approved for lung cancer, but its use for the highly heterogeneous disease of GC is unknown. The goal of this study was to investigate the anti-cancer mechanisms of the (S)-crizotinib in inhibiting GC growth. Human GC cell lines (SGC-7901 and BGC-823) and the (S)-crizotinib-resistant BGC-823/R were cultured for determining the effects of (S)-crizotinib on cell viability, apoptosis, oxidant generation, and cell cycle progression. Involvement of ROS, Akt signaling, MTH1, and DNA damage was tested with respective pharmacological blockade. The in vivo anti-tumor effects of (S)-crizotinib were determined using xenograft tumor mice. Results indicated that (S)-crizotinib decreased GC cell viability, induced growth arrest and apoptosis, and increased levels of γH2AX and Ser1981-phosphorylated ATM, which were inhibited by NAC. The anti-cancer mechanism of (S)-crizotinib was independent of MTH1. Moreover, ATM-activated Akt, a pro-survival signal, whose inhibition further enhanced (S)-crizotinib-induced inhibition of GC cell growth and tumor growth in xenograft mice, and re-sensitized resistant GC cells to (S)-crizotinib. (S)-crizotinib reduced GC cell and tumor growth through oxidative DNA damage mechanism and triggered pro-survival Akt signaling. We conclude that inclusion of Akt inhibition (to block the survival signaling) with (S)-crizotinib may provide an effective and novel combination therapy for GC in the clinical setting.

The binding mechanism between azoles and FgCYP51B, sterol 14α-demethylase of Fusarium graminearum

BACKGROUND

Fusarium graminearum is the main pathogen of Fusarium head blight (FHB), a worldwide plant disease and a major disease of wheat in China. Control of FHB is mainly dependent on the application of demethylase inhibitor (DMI) fungicides. Fungal sterol 14α-demethylase enzymes (CYP51) are the main target for DMI fungicides. A molecular modeling study and biological evaluation were performed to investigate the binding mechanism between azoles and CYP51B in F. graminearum.

RESULTS

A homology model based on the crystal structure of Aspergillus fumigatus was built. Molecular docking and molecular dynamics (MD) simulations were then used to identify the optimum binding mode of propiconazole (PRP), diniconazole (DIN), triadimenol (TRL), tebuconazole (TEC) and triadimefon (TRN) with FgCYP51B. Furthermore, the binding free energy of the five protein-inhibitor complexes was calculated using molecular mechanics generalized Born surface area and Poisson-Boltzmann surface area (MM-GB/PBSA) methods. Key residues in the selective binding of azoles to FgCYP51B were recognized by per-residue free energy decomposition analysis. The five ligands have a similar binding mode in the active pocket. The binding free energy to the enzyme for inhibitors PRP and TEC is more favorable than that of TRN, TRL and DIN. Furthermore, the amino acid residues Phe511, Val136, Ile374, Ala308, Ser312 and Try137 of FgCYP51B are key residues interacting with azoles fungicides. From the experimental evaluation, the 50% effective concentration (EC₅₀ ) values for PRP, TEC, DIN, TRL and TRN are 0.024, 0.047, 0.148, 0.154 and 0.474 mg L^-1 , respectively. These five molecules exhibit potential inhibitory activity against CYP51B protein from F. graminearum.

CONCLUSION

Azole fungicides for FgCYP51B should possess more hydrophobic groups interacting with residues Phe511, Val136, Ile374, Ala308, Ser312 and Tyr137. PRP and TEC are preferable for the control of FHB than DIN, TRL and TRN. © 2017 Society of Chemical Industry.

Theoretical rationalisation of the photophysics of a TICT excited state of cinnamoyl–coumarin derivatives in homogeneous and biological membrane models

Analysis of the potential energy barriers and structural dynamics of a new TICT-probe for monitoring biological environments.

Association of MTH1 expression with the tumor malignant potential and poor prognosis in patients with resected lung cancer

OBJECTIVES

The oxidized purine nucleoside triphosphatase, mutT homolog 1 (MTH1), physiologically sanitizes 8-oxo-dGTP in the nucleotide pool. Previous studies indicated that MTH1 is associated with tumor proliferation and invasion in non-small cell lung cancer (NSCLC) cell lines; however, the role of MTH1 in patients with NSCLC remains unclear.

MATERIALS AND METHODS

Two patient cohorts that underwent surgery for NSCLC in our institution were investigated retrospectively. In one cohort consisting of 197 patients, the associations between MTH1 expression and clinicopathological factors or prognosis were analyzed. In another cohort consisting of 41 patients, the relationship between MTH1 expression in the tumors and serum oxidative stress levels (evaluated by the diacron-reactive oxygen metabolites [d-ROMs] test) or antioxidant capacity in the patients (evaluated by the biological antioxidant potential (BAP) test) was analyzed. A total of 238 patients were assessed for MTH1 protein levels using immunohistochemistry.

RESULTS

Among the 197 patients in the former cohort, 111 (56.3%) exhibited high MTH1 expression, while 86 (43.7%) exhibited low MTH1 expression. Male sex, smoking habit of ≥20 pack-years, squamous cell carcinoma, pathological stage ≥ II, tumor diameter ≥30mm, lymph node metastases, pleural invasion, lymphatic permeation and vascular infiltration were significantly associated with high MTH1 expression (p<0.05). The high MTH1 expression group had a significantly worse prognosis than that of the low MTH1 expression group (5-year overall survival: 81.6% vs. 92.3%, p=0.0011; 5-year disease-free survival: 55.0% vs. 83.7%, p=0.0002). d-ROMs and BAP test values were significantly higher in the high than in the low MTH1 expression group (p<0.05).

CONCLUSION

This study showed that MTH1 protein expression was closely related to factors associated with a high malignant potential and poor patient survival. MTH1 may be a novel therapeutic target for NSCLC.

Revealing inhibition difference between PFI-2 enantiomers against SETD7 by molecular dynamics simulations, binding free energy calculations and unbinding pathway analysis

SETD7 is associated with multiple diseases related signaling pathways. (R)-PFI-2 is the first SETD7 inhibitor with nanomolar inhibitory potency. The activity of (R)-PFI-2 is about 500 times over that of (S)-PFI-2. Understanding the mechanism behind this difference will be helpful to discovery and design of more potent SETD7 inhibitors. A computational study combining molecular dynamics simulation, binding free energy calculations, and residue interaction network (RIN) was performed on the (S)-PFI-2/SETD7 and (R)-PFI-2/SETD7 complexes to explore the molecular mechanism behind the different inhibition activity. The results from Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) calculation show (R)-PFI-2 has lower binding free energy. Residues H252, D256, L267, Y335, G336 and H339 are responsible for the binding of SETD7 to the (R)-PFI-2. RIN analysis indicates van der Waals interaction is critical for the binding of (R)-PFI-2. The results from adaptive basing force (ABF) simulation confirm that the free energy barrier of (R)-PFI-2 dissociating from the SETD7 is larger than that of (S)-PFI-2. (S)-PFI-2 and (R)-PFI-2 dissociate from the SETD7 binding site along different reaction coordinate and have potential mean of force (PMF) depth. Our simulations results will be useful to understand molecular mechanism of activity difference between PFI-2 enantiomers against SETD7.

Understanding the molecular mechanism for the differential inhibitory activities of compounds against MTH1

MTH1 can hydrolyze oxidized nucleotides and is required for cancer survival. The IC₅₀ values were 0.8 nM for TH287 with a methyl substitution, 5.0 nM for TH588 with a cyclopropyl substitution, and 2.1 μM for TH650 with an oxetanyl substitution. Thus, it is very significant to understand inhibitory mechanisms of these structurally similar compounds against MTH1 and influences of the substituent on the bioactivities. Our MD researches indicate that TH287 maintains significant hydrogen bonds with Asn33 and Asp119, stabilizes the binding site, and induces MTH1 adopt a closed motion, leading to a high inhibitory activity. When bound with TH588, the binding site can be partially stabilized and take a semi-closed state, which is because the cyclopropyl group in TH588 has larger steric hindrance than a methyl group in TH287. So TH588 has a slightly reduced inhibitory activity compared to TH287. TH650 induces greater conformation fluctuations than TH588 and the binding site adopts an opening state, which is caused by the large bulk of oxetanyl group and the interference of solvent on the oxetanyl substituent, leading to the lowest inhibitory activity. Thus, the inhibitory activity follows a TH287 > TH588 > TH650 trend, which well matches with the experimental finding.

Peptide interactions with zigzag edges in graphene

Recognition and manipulation of graphene edges enable the control of physical properties of graphene-based devices. Recently, the authors have identified a peptide that preferentially binds to graphene edges from a combinatorial peptide library. In this study, the authors examine the functional basis for the edge binding peptide using experimental and computational methods. The effect of amino acid substitution, sequence context, and solution pH value on the binding of the peptide to graphene has been investigated. The N-terminus glutamic acid residue plays a key role in recognizing and binding to graphene edges. The protonation, substitution, and positional context of the glutamic acid residue impact graphene edge-binding. Our findings provide insights into the binding mechanisms and the design of peptides for recognizing and functionalizing graphene edges.