Small molecule induction of MSH2-dependent cell death suggests a vital role of mismatch repair proteins in cell death

Oestrogen receptor positive breast cancer and its embedded mechanism: breast cancer resistance to conventional drugs and related therapies, a review

Traditional medication and alternative therapies have long been used to treat breast cancer. One of the main problems with current treatments is that there is an increase in drug resistance in the cancer cells owing to genetic differences such as mutational changes, epigenetic changes and miRNA (microRNA) alterations such as miR-1246, miR-298, miR-27b and miR-33a, along with epigenetic modifications, such as Histone3 acetylation and CCCTC-Binding Factor (CTCF) hypermethylation for drug resistance in breast cancer cell lines. Certain forms of conventional drug resistance have been linked to genetic changes in genes such as ABCB1, AKT, S100A8/A9, TAGLN2 and NPM. This review aims to explore the current approaches to counter breast cancer, the action mechanism, along with novel therapeutic methods endowing potential drug resistance. The investigation of novel therapeutic approaches sheds light on the phenomenon of drug resistance including genetic variations that impact distinct forms of oestrogen receptor (ER) cancer, genetic changes, epigenetics-reported resistance and their identification in patients. Long-term effective therapy for breast cancer includes selective oestrogen receptor modulators, selective oestrogen receptor degraders and genetic variations, such as mutations in nuclear genes, epigenetic modifications and miRNA alterations in target proteins. Novel research addressing combinational therapies including maytansine, photodynamic therapy, guajadiol, talazoparib, COX2 inhibitors and miRNA 1246 inhibitors have been developed to improve patient survival rates.

Anti-angiogenic potential of bioactive phytochemicals from Helicteres isora targeting VEGFR-2 to fight cancer through molecular docking and molecular dynamics simulation

Cancer is one of the leading causes of death due to its very high rate of morbidity and mortality, and there is a constant demand of effective drugs for cancer therapy. Vascular endothelial growth factor receptor-2 (VEGFR-2) plays a significant role as central modulator of angiogenesis and is targeted frequently for developing anti-angiogenic agents to fight cancer. Helicteres isora L. (Malvaceae) is reported to possess diverse medicinal properties including anti-cancer potentials. In the current investigation, 38 bioactive phytochemicals of H. isora were screened virtually to evaluate their anti-angiogenic potentials targeting VEGFR-2. The study unveiled three potential candidates such as, Diosgenin (-9.8 Kcal//mol), Trifolin (-8.4 Kcal/mol) and Yohimbine (-8.1 Kcal/mol) that showed favorable pharmacokinetic, pharmacodynamic and toxicity properties with no significant side effects. Molecular dynamics simulation employing 100 ns revealed noteworthy structural stability and compactness for all the top three candidates. The MM/GBSA binding free energy estimation corroborated the docking interactions where Yohimbine (-30.47 Kcal/mol) scored better than Diosgenin (-27.54 Kcal/mol) and Trifolin (-29.58 Kcal/mol). Target class prediction revealed enzymes in most of the cases and some FDA approved drugs were found as structurally similar analogs for Trifolin and Yohimbine. These findings could lead to the development of novel and effective anti-angiogenic agents.Communicated by Ramaswamy H. Sarma.

CircLIFR synergizes with MSH2 to attenuate chemoresistance via MutSα/ATM-p73 axis in bladder cancer

Background

Cisplatin (CDDP) has become a standard-of-care treatment for muscle-invasive bladder cancer (MIBC), while chemoresistance remains a major challenge. Accumulating evidence indicates that circular RNAs (circRNAs) are discrete functional entities. However, the regulatory functions as well as complexities of circRNAs in modulating CDDP-based chemotherapy in bladder cancer are yet to be well revealed.

Methods

Through analyzing the expression profile of circRNAs in bladder cancer tissues, RNA FISH, circRNA pull-down assay, mass spectrometry analysis and RIP, circLIFR was identified and its interaction with MSH2 was confirmed. The effects of circLIFR and MSH2 on CDDP-based chemotherapy were explored by flow cytometry and rescue experiments. Co-IP and Western blot were used to investigate the molecular mechanisms underlying the functions of circLIFR and MSH2. Biological implications of circLIFR and MSH2 in bladder cancer were implemented in tumor xenograft models and PDX models.

Results

CircLIFR was downregulated in bladder cancer and expression was positively correlated with favorable prognosis. Moreover, circLIFR synergizing with MSH2, which was a mediator of CDDP sensitivity in bladder cancer cells, positively modulated sensitivity to CDDP in vitro and in vivo. Mechanistically, circLIFR augmented the interaction between MutSα and ATM, ultimately contributing to stabilize p73, which triggered to apoptosis. Importantly, MIBC with high expression of circLIFR and MSH2 was more sensitive to CDDP-based chemotherapy in tumor xenograft models and PDX models.

Conclusions

CircLIFR could interact with MSH2 to positively modulate CDDP-sensitivity through MutSα/ATM-p73 axis in bladder cancer. CircLIFR and MSH2 might be act as promising therapeutic targets for CDDP-resistant bladder cancer.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12943-021-01360-4.

Reserpine inhibits DNA repair, cell proliferation, invasion and induces apoptosis in oral carcinogenesis via modulation of TGF-β signaling

Reserpine is a natural indole alkaloid isolated from Rauwolfia serpentina and has potent antioxidant, antimicrobial, and anti-mutagenic properties. Accordingly, this study aimed to investigate the effect of reserpine on DNA repair, cell proliferation, invasion and apoptosis in 7,12-dimethylbenz[a]anthracene(DMBA)-induced hamster buccal pouch (HBP) carcinogenesis. Transforming growth factor-β (TGF-β) was found to induce Smad2, 3 and 4 phosphorylation triggering Smad3/Snail mediated DNA repair proteins and Smad2/4 nuclear translocation. In contrast, reserpine inhibits TGF-β dependent Smad2/3/4 phosphorylation, thereby blockage Smad3/Snail activation and Smad2/4 nuclear translocation. Interruption of these oncogenic signaling pathways leads to downregulating ERCC1, XPF, Ku70, DNA-PKcs, PCNA, cyclin D1, HIF-1α, IL-6, Mcl-1 and stimulates Bax, cytochrome C, Apaf-1, caspase-9, caspase-3 and PARP protein expressions. This study provides therapeutic potential of reserpine in inhibiting DNA repair, cell proliferation, and invasion while simultaneously inducing apoptosis via modulation TGF-β signals.

All-atom molecular dynamics comparison of disease-associated zinc fingers

An important regulatory domain of NF-[Formula: see text]B Essential Modulator (NEMO) is a ubiquitin-binding zinc finger, with a tetrahedral CYS3HIS1 zinc-coordinating binding site. Two variations of NEMO's zinc finger are implicated in various disease states including ectodermal dysplasia and adult-onset glaucoma. To discern structural and dynamical differences between these disease states, we present results of 48-[Formula: see text]s of molecular dynamics simulations for three zinc finger systems each in two states, with and without zinc-bound and correspondingly appropriate cysteine thiol/thiolate configurations. The wild-type protein, often studied for its role in cancer, maintains the most rigid and conformationally stable zinc-bound configuration compared with the diseased counterparts. The glaucoma-related protein has persistent loss of secondary structure except within the dominant conformation. Conformational overlap between wild-type and glaucoma isoforms indicate a competitive binding mechanism may be substantial in the malfunctioning configuration, while the alpha-helical disruption of the ectodermal dysplasia suggests a loss of binding selectivity is responsible for aberrant function.

Allosteric binding sites in Rab11 for potential drug candidates

Rab11 is an important protein subfamily in the RabGTPase family. These proteins physiologically function as key regulators of intracellular membrane trafficking processes. Pathologically, Rab11 proteins are implicated in many diseases including cancers, neurodegenerative diseases and type 2 diabetes. Although they are medically important, no previous study has found Rab11 allosteric binding sites where potential drug candidates can bind to. In this study, by employing multiple clustering approaches integrating principal component analysis, independent component analysis and locally linear embedding, we performed structural analyses of Rab11 and identified eight representative structures. Using these representatives to perform binding site mapping and virtual screening, we identified two novel binding sites in Rab11 and small molecules that can preferentially bind to different conformations of these sites with high affinities. After identifying the binding sites and the residue interaction networks in the representatives, we computationally showed that these binding sites may allosterically regulate Rab11, as these sites communicate with switch 2 region that binds to GTP/GDP. These two allosteric binding sites in Rab11 are also similar to two allosteric pockets in Ras that we discovered previously.

Visualizing correlated motion with HDBSCAN clustering

Correlated motion analysis provides a method for understanding communication between and dynamic similarities of biopolymer residues and domains. The typical equal-time correlation matrices-frequently visualized with pseudo-colorings or heat maps-quickly convey large regions of highly correlated motion but hide more subtle similarities of motion. Here we propose a complementary method for visualizing correlations within proteins (or general biopolymers) that quickly conveys intuition about which residues have a similar dynamic behavior. For grouping residues, we use the recently developed non-parametric clustering algorithm HDBSCAN. Although the method we propose here can be used to group residues using correlation as a similarity matrix-the most straightforward and intuitive method-it can also be used to more generally determine groups of residues which have similar dynamic properties. We term these latter groups "Dynamic Domains", as they are based not on spatial closeness but rather closeness in the column space of a correlation matrix. We provide examples of this method across three human proteins of varying size and function-the Nf-Kappa-Beta essential modulator, the clotting promoter Thrombin and the mismatch repair protein (dimer) complex MutS-alpha. Although the examples presented here are from all-atom molecular dynamics simulations, this visualization technique can also be used on correlations matrices built from any ensembles of conformations from experiment or computation.

Molecular dynamics simulations of aptamer-binding reveal generalized allostery in thrombin

Thrombin is an attractive target for antithrombotic therapy due to its central role in thrombosis and hemostasis as well as its role in inducing tumor growth, metastasis, and tumor invasion. The thrombin-binding DNA aptamer (TBA), is under investigation for anticoagulant drugs. Although aptamer binding experiments have been revealed various effects on thrombin's enzymatic activities, the detailed picture of the thrombin's allostery from TBA binding is still unclear. To investigate thrombin's response to the aptamer-binding at the molecular level, we compare the mechanical properties and free energy landscapes of the free and aptamer-bound thrombin using microsecond-scale all-atom GPU-based molecular dynamics simulations. Our calculations on residue fluctuations and coupling illustrate the allosteric effects of aptamer-binding at the atomic level, highlighting the exosite II, 60s, γ and the sodium loops, and the alpha helix region in the light chains involved in the allosteric changes. This level of details clarifies the mechanisms of previous experimentally demonstrated phenomena, and provides a prediction of the reduced autolysis rate after aptamer-binding. The shifts in thrombin's ensemble of conformations and free energy surfaces after aptamer-binding demonstrate that the presence of bound-aptamer restricts the conformational freedom of thrombin suggesting that conformational selection, i.e. generalized allostery, is the dominant mechanism of thrombin-aptamer binding. The profound perturbation on thrombin's mechanical and thermodynamic properties due to the aptamer-binding, which was revealed comprehensively as a generalized allostery in this work, may be exploited in further drug discovery and development.

Mechanistic insights into thrombin's switch between "slow" and "fast" forms

Thrombin is a multifunctional enzyme that plays an important role in blood coagulation, cell growth, and metastasis. Depending upon the binding of sodium ions, thrombin presents significantly different enzymatic activities. In the environment with sodium ions, thrombin is highly active in cleaving the coagulated substrates and this is referred to as the "fast" form; in the environment without sodium ions, thrombin turns catalytically less active and is in the "slow" form. Although many experimental studies over the last two decades have attempted to reveal the structural and kinetic differences between these two forms, it remains vague and disputed how the functional switch between the "fast" and "slow" forms is mediated by Na+ cations. In this work, we employ microsecond-scale all-atom molecular dynamics simulations to investigate the differences in the structural ensembles in sodium-bound/unbound and potassium-bound/unbound thrombin. Our calculations indicate that the regulatory regions, including the 60s, γ loops, and exosite I and II, are primarily affected by both the bound and unbound cations. Conformational free energy surfaces, estimated from principal component analysis, further reveal the existence of multiple conformational states. The binding of a cation introduces changes in the distribution of these states. Through comparisons with potassium-binding, the binding of sodium ions appears to shift the population toward conformational states that might be catalytically favorable. Our study of thrombin in the presence of sodium/potassium ions suggests Na+-mediated generalized allostery is the mechanism of thrombin's functional switch between the "fast" and "slow" forms.

Capturing the photo-signaling state of a photoreceptor in a steady-state fashion by binding a transition metal complex

Binding a small molecule to proteins causes conformational changes, but often to a limited extent. Here, we demonstrate that the interaction of a CO-releasing molecule (CORM3) with a photoreceptor photoactive yellow protein (PYP) drives large structural changes in the latter. The interaction of CORM3 and a mutant of PYP, Met100Ala, not only trigger the isomerization of its chromophore, p-coumaric acid, from its anionic trans configuration to a protonated cis configuration, but also increases the content of β-sheet at the cost of α-helix and random coil in the secondary structure of the protein. The CORM3 derived Met100Ala is found to highly resemble the signaling state, which is one of the key photo-intermediates of this photoactive protein, in both protein local conformation and chromophore configuration. The organometallic reagents hold promise as protein engineering tools. This work highlights a novel approach to structurally accessing short lived intermediates of proteins in a steady-state fashion.

MutSα's Multi-Domain Allosteric Response to Three DNA Damage Types Revealed by Machine Learning

MutSα is a key component in the mismatch repair (MMR) pathway. This protein is responsible for initiating the signaling pathways for DNA repair or cell death. Herein we investigate this heterodimer's post-recognition, post-binding response to three types of DNA damage involving cytotoxic, anti-cancer agents-carboplatin, cisplatin, and FdU. Through a combination of supervised and unsupervised machine learning techniques along with more traditional structural and kinetic analysis applied to all-atom molecular dynamics (MD) calculations, we predict that MutSα has a distinct response to each of the three damage types. Via a binary classification tree (a supervised machine learning technique), we identify key hydrogen bond motifs unique to each type of damage and suggest residues for experimental mutation studies. Through a combination of a recently developed clustering (unsupervised learning) algorithm, RMSF calculations, PCA, and correlated motions we predict that each type of damage causes MutSα to explore a specific region of conformation space. Detailed analysis suggests a short range effect for carboplatin-primarily altering the structures and kinetics of residues within 10 angstroms of the damaged DNA-and distinct longer-range effects for cisplatin and FdU. In our simulations, we also observe that a key phenylalanine residue-known to stack with a mismatched or unmatched bases in MMR-stacks with the base complementary to the damaged base in 88.61% of MD frames containing carboplatinated DNA. Similarly, this Phe71 stacks with the base complementary to damage in 91.73% of frames with cisplatinated DNA. This residue, however, stacks with the damaged base itself in 62.18% of trajectory frames with FdU-substituted DNA and has no stacking interaction at all in 30.72% of these frames. Each drug investigated here induces a unique perturbation in the MutSα complex, indicating the possibility of a distinct signaling event and specific repair or death pathway (or set of pathways) for a given type of damage.

Uncovering Large-Scale Conformational Change in Molecular Dynamics without Prior Knowledge

As the length of molecular dynamics (MD) trajectories grows with increasing computational power, so does the importance of clustering methods for partitioning trajectories into conformational bins. Of the methods available, the vast majority require users to either have some a priori knowledge about the system to be clustered or to tune clustering parameters through trial and error. Here we present non-parametric uses of two modern clustering techniques suitable for first-pass investigation of an MD trajectory. Being non-parametric, these methods require neither prior knowledge nor parameter tuning. The first method, HDBSCAN, is fast-relative to other popular clustering methods-and is able to group unstructured or intrinsically disordered systems (such as intrinsically disordered proteins, or IDPs) into bins that represent global conformational shifts. HDBSCAN is also useful for determining the overall stability of a system-as it tends to group stable systems into one or two bins-and identifying transition events between metastable states. The second method, iMWK-Means, with explicit rescaling followed by K-Means, while slower than HDBSCAN, performs well with stable, structured systems such as folded proteins and is able to identify higher resolution details such as changes in relative position of secondary structural elements. Used in conjunction, these clustering methods allow a user to discern quickly and without prior knowledge the stability of a simulated system and identify both local and global conformational changes.

Cytotoxicity of the indole alkaloid reserpine from Rauwolfia serpentina against drug-resistant tumor cells

BACKGROUND

The antihypertensive reserpine is an indole alkaloid from Rauwolfia serpentina and exerts also profound activity against cancer cells in vitro and in vivo. The present investigation was undertaken to investigate possible modes of action to explain its activity toward drug-resistant tumor cells.

MATERIAL AND METHODS

Sensitive and drug-resistant tumor cell lines overexpressing P-glycoprotein (ABCB1/MDR1), breast cancer resistance protein (ABCG2/BCRP), mutation-activated epidermal growth factor receptor (EGFR), wild-type and p53-knockout cells as well as the NCI panel of cell lines from different tumor origin were analyzed. Reserpine's cytotoxicity was investigated by resazurin and sulforhodamine assays, flow cytometry, and COMPARE and hierarchical cluster analyses of transcriptome-wide microarray-based RNA expressions.

RESULTS

P-glycoprotein- or BCRP overexpressing tumor cells did not reveal cross-resistance to reserpine. EGFR-overexpressing cells were collateral sensitive and p53- Knockout cells cross-resistant to this drug compared to their wild-type parental cell lines. Reserpine increased the uptake of doxorubicin in P-glycoprotein-overexpressing cells, indicating that reserpine inhibited the efflux function of P-glycoprotein. Using molecular docking, we found that reserpine bound with even higher binding energy to P-glycoprotein and EGFR than the control drugs verapamil (P-glycoprotein inhibitor) and erlotinib (EGFR inhibitor). COMPARE and cluster analyses of microarray data showed that the mRNA expression of a panel of genes predicted the sensitivity or resistance of the NCI tumor cell line panel with statistical significance. The genes belonged to diverse pathways and biological functions, e.g. cell survival and apoptosis, EGFR activation, regulation of angiogenesis, cell mobility, cell adhesion, immunological functions, mTOR signaling, and Wnt signaling.

CONCLUSION

The lack of cross-resistance to most resistance mechanisms and the collateral sensitivity in EGFR-transfectants compared to wild-type cells speak for a promising role of reserpine in cancer chemotherapy. Reserpine deserves further consideration for cancer therapy in the clinical setting.

Destabilization of the MutSα's protein-protein interface due to binding to the DNA adduct induced by anticancer agent carboplatin via molecular dynamics simulations

DNA mismatch repair (MMR) proteins maintain genetic integrity in all organisms by recognizing and repairing DNA errors. Such alteration of hereditary information can lead to various diseases, including cancer. Besides their role in DNA repair, MMR proteins detect and initiate cellular responses to certain type of DNA damage. Its response to the damaged DNA has made the human MMR pathway a useful target for anticancer agents such as carboplatin. This study indicates that strong, specific interactions at the interface of MutSα in response to the mismatched DNA recognition are replaced by weak, non-specific interactions in response to the damaged DNA recognition. Data suggest a severe impairment of the dimerization of MutSα in response to the damaged DNA recognition. While the core of MutSα is preserved in response to the damaged DNA recognition, the loss of contact surface and the rearrangement of contacts at the protein interface suggest a different packing in response to the damaged DNA recognition. Coupled in response to the mismatched DNA recognition, interaction energies, hydrogen bonds, salt bridges, and solvent accessible surface areas at the interface of MutSα and within the subunits are uncoupled or asynchronously coupled in response to the damaged DNA recognition. These pieces of evidence suggest that the loss of a synchronous mode of response in the MutSα's surveillance for DNA errors would possibly be one of the mechanism(s) of signaling the MMR-dependent programed cell death much wanted in anticancer therapies. The analysis was drawn from dynamics simulations.

Non-specificity and synergy at the binding site of the carboplatin-induced DNA adduct via molecular dynamics simulations of the MutSα-DNA recognition complex

MutSα is the most abundant mismatch-binding factor of human DNA mismatch repair (MMR) proteins. MMR maintains genetic stability by recognizing and repairing DNA defects. Failure to accomplish their function may lead to cancer. In addition, MutSα recognizes at least some types of DNA damage making it a target for anticancer agents. Here, complementing scarce experimental data, we report unique hydrogen-bonding motifs associated with the recognition of the carboplatin induced DNA damage by MutSα. These data predict that carboplatin and cisplatin induced damaging DNA adducts are recognized by MutSα in a similar manner. Our simulations also indicate that loss of base pairing at the damage site results in (1) non-specific binding and (2) changes in the atomic flexibility at the lesion site and beyond. To further quantify alterations at MutSα-DNA interface in response to damage recognition, non-bonding interactions and salt bridges were investigated. These data indicate (1) possible different packing and (2) disruption of the salt bridges at the MutSα-DNA interface in the damaged complex. These findings (1) underscore the general observation of disruptions at the MutSα-DNA interface and (2) highlight the nature of the anticancer effect of the carboplatin agent. The analysis was carried out from atomistic simulations.

Transcription pattern of p53-targeted DNA repair genes in the hypoxia-tolerant subterranean mole rat Spalax

The tumor suppressor gene p53 induces growth arrest and/or apoptosis in response to DNA damage/hypoxia. Inactivation of p53 confers a selective advantage to tumor cells under a hypoxic microenvironment during tumor progression. The subterranean blind mole rat, Spalax, spends its life underground at low-oxygen tensions, hence developing a wide range of respiratory/molecular adaptations to hypoxic stress, including critical changes in p53 structure and signaling pathway. The highly conserved p53 Arg(R)-172 is substituted by lysine (K) in Spalax, identical with a tumor-associated mutation. Functionality assays revealed that Spalax p53 is unable to activate apoptotic target genes but is still capable of activating cell cycle arrest genes. Furthermore, we have shown that the transcription patterns of representative p53-induced genes (Apaf1 and Mdm2) in Spalax are influenced by hypoxia. Cell cycle arrest allows the cells to repair DNA damage via different DNA repair genes. We tested the transcription pattern of three p53-related DNA repair genes (p53R2, Mlh1, and Msh2) under normoxia and short-acute hypoxia in Spalax, C57BL/6 wild-type mice, and two strains of mutant C57BL/6 mice, each carrying a different mutation at the R172 position. Our results show that while wild-type/mutant mice exhibit strong hypoxia-induced reductions of repair gene transcript levels, no such inhibition is found in Spalax under hypoxia. Moreover, unlike mouse p53R2, Spalax p53R2 transcript levels are strongly elevated under hypoxia. These results suggest that critical repair functions, which are known to be inhibited under hypoxia in mice, remain active in Spalax, as part of its unique hypoxia tolerance mechanisms.

Computational and synthetic studies towards improving rescinnamine as an inducer of MSH2-dependent apoptosis in cancer treatment

We, and others, have previously shown that mismatch repair proteins, in addition to their repair function, contribute to cell death initiation. In response to some drugs, this cell death activity is independent of the repair function of the proteins. Rescinnamine, a derivative of the indole alkaloid reserpine, a drug used to treat hypertension several decades ago, was shown to target the cell death-initiating activity of mismatch repair proteins. When used in animals, the hypotensive action of this drug prevents applying appropriate concentrations for statistically significant tumor reduction. Using a combination of computational modeling, chemical synthesis and cell assays, we determine how rescinnamine can be structurally modified and what effect these modifications have on cell survival. These results inform further computational modeling to suggest new synthetic lead molecules to move toward further biological testing.

Molecular dynamics simulations of protein dynamics and their relevance to drug discovery

Molecular dynamics simulations have become increasingly useful in studying biological systems of biomedical interest, and not just in the study of model or toy systems. In this article, the methods and principles of all-atom molecular dynamics will be elucidated with several examples provided of their utility to investigators interested on drug discovery.

Parameters of Reserpine Analogs That Induce MSH2/MSH6-Dependent Cytotoxic Response

Mismatch repair proteins modulate the cytotoxicity of several chemotherapeutic agents. We have recently proposed a "death conformation" of the MutS homologous proteins that is distinguishable from their "repair conformation." This conformation can be induced by a small molecule, reserpine, leading to DNA-independent cell death. We investigated the parameters for a small reserpine-like molecule that are required to interact with MSH2/MSH6 to induce MSH2/MSH6-dependent cytotoxic response. A multidisciplinary approach involving structural modeling, chemical synthesis, and cell biology analyzed reserpine analogs and modifications. We demonstrate that the parameters controlling the induction of MSH2/MSH6-dependent cytotoxicity for reserpine-analogous molecules reside in the specific requirements for methoxy groups, the size of the molecule, and the orientation of molecules within the protein-binding pocket. Reserpine analog rescinnamine showed improved MSH2-dependent cytotoxicity. These results have important implications for the identification of compounds that require functional MMR proteins to exhibit their full cytotoxicity, which will avoid resistance in MMR-deficient cells.