Search

Tomasz Kostrzewa,1 Kamlesh K Sahu,2 Magdalena Gorska-Ponikowska,1 Jack A Tuszynski,3 Alicja Kuban-Jankowska1 1Department of Medical Chemistry, Medical University of Gdansk, Gdansk, Poland; 2Department of Medical Microbiology and Immunology, University of Alberta, Edmonton, AB, Canada; 3Department of Oncology, University of Alberta, Edmonton, AB, Canada Background: The protein tyrosine phosphatases PTP1B and SHP2 are promising drug targets in treatment design for breast cancer. Searching for specific inhibitors of their activity has recently become the challenge of many studies. Previous work has indicated that the promising PTP inhibitors may be small compounds that are able to bind and interact with amino residues from the binding site.Purpose: The main goal of our study was to synthesize and analyze the effect of selected small peptide inhibitors on oncogenic PTP1B and SHP2 enzymatic activity and viability of MCF7 breast cancer cells. We also performed computational analysis of peptides binding with allosteric sites of PTP1B and SHP2 phosphatases.Methods: We measured the inhibitory activity of compounds utilizing recombinant enzymes and MCF7 cell line. Computational analysis involved docking studies of binding conformation and interactions of inhibitors with allosteric sites of phosphatases.Results: The results showed that the tested compounds decrease the enzymatic activity of phosphatases PTP1B and SHP2 with IC50 values in micromolar ranges. We observed higher inhibitory activity of dipeptides than tripeptides. Phe-Asp was the most effective against SHP2 enzymatic activity, with IC50=5.2±0.4 µM. Micromolar concentrations of tested dipeptides also decreased the viability of MCF7 breast cancer cells, with higher inhibitory activity observed for the Phe-Asp peptide. Moreover, the peptides tested were able to bind and interact with allosteric sites of PTP1B and SHP2 phosphatases.Conclusion: Our research showed that small peptide compounds can be considered for the design of specific inhibitors of oncogenic protein tyrosine phosphatases. Keywords: breast cancer, protein tyrosine phosphatase PTP1B, SHP2, peptides, PTP inhibitors

Abstract Background A maximum entropy approach is proposed to predict the cytotoxic effects of a panel of colchicine derivatives in several human cancer cell lines. Data was obtained from cytotoxicity assays performed with 21 drug molecules from the same family of colchicine compounds and correlate these results with independent tubulin isoform expression measurements for several cancer cell lines. The maximum entropy method is then used in conjunction with computed relative binding energy values for each of the drug molecules against tubulin isotypes to which these compounds bind with different affinities. Results We have found by using our analysis that αβI and αβIII tubulin isoforms are the most important isoforms in establishing predictive response of cancer cell sensitivity to colchicine derivatives. However, since αβI tubulin is widely distributed in the human body, targeting it would lead to severe adverse side effects. Consequently, we have identified tubulin isotype αβIII as the most important molecular target for inhibition of microtubule polymerization and hence cancer cell cytotoxicity. Tubulin isotypes αβI and αβII are concluded to be secondary targets. Conclusions The benefit of being able to correlate expression levels of specific tubulin isotypes and the resultant cell death effect is that it will enable us to better understand the origin of drug resistance and hence design optimal structures for the elimination of cancer cells. The conclusion of the study described herein identifies tubulin isotype αβIII as a target for optimized chemotherapy drug design.

Abstract Background Since the transfer and application of modern sequencing technologies to the analysis of amplified fragment-length polymorphisms (AFLP), evolutionary biologists have included an increasing number of samples and markers in their studies. Although justified in this context, the use of automated scoring procedures may result in technical biases that weaken the power and reliability of further analyses. Results Using a new scoring algorithm, RawGeno, we show that scoring errors – in particular "bin oversplitting" (i.e. when variant sizes of the same AFLP marker are not considered as homologous) and "technical homoplasy" (i.e. when two AFLP markers that differ slightly in size are mistakenly considered as being homologous) – induce a loss of discriminatory power, decrease the robustness of results and, in extreme cases, introduce erroneous information in genetic structure analyses. In the present study, we evaluate several descriptive statistics that can be used to optimize the scoring of the AFLP analysis, and we describe a new statistic, the information content per bin (Ibin) that represents a valuable estimator during the optimization process. This statistic can be computed at any stage of the AFLP analysis without requiring the inclusion of replicated samples. Finally, we show that downstream analyses are not equally sensitive to scoring errors. Indeed, although a reasonable amount of flexibility is allowed during the optimization of the scoring procedure without causing considerable changes in the detection of genetic structure patterns, notable discrepancies are observed when estimating genetic diversities from differently scored datasets. Conclusion Our algorithm appears to perform as well as a commercial program in automating AFLP scoring, at least in the context of population genetics or phylogeographic studies. To our knowledge, RawGeno is the only freely available public-domain software for fully automated AFLP scoring, from electropherogram files to user-defined working binary matrices. RawGeno was implemented in an R CRAN package (with an user-friendly GUI) and can be found at http://sourceforge.net/projects/rawgeno.

The Davydov model was conjectured to describe how an amide I excitation created during ATP hydrolysis in myosin might be significant in providing energy to drive myosin's chemomechanical cycle. The free energy surfaces of the myosin relay helix peptide dissolved in 2,2,2-trifluoroethanol (TFE), determined by metadynamics simulations, demonstrate local minima differing in free energy by only ~2 kT, corresponding to broken and stabilized hydrogen bonds, respectively. Experimental pump-probe and 2D infrared spectroscopy were performed on the peptide dissolved in TFE. The relative heights of two peaks seen in the pump-probe data and the corresponding relative volumes of diagonal peaks seen in the 2D-IR spectra at time delays between 0.5 ps and 1 ps differ noticeably from what is seen at earlier or later time delays or in the linear spectrum, indicating that a vibrational excitation may influence the conformational state of this helix. Thus, it is possible that the presence of an amide I excitation may be a direct factor in the conformational state taken on by the myosin relay helix following ATP hydrolysis in myosin.

Actin filaments, as key components of the cytoskeleton, have aroused great interest due to their numerous functional roles in eukaryotic cells, including intracellular electrical signaling. The aim of this research is to characterize the alternating current (AC) conduction characteristics of both globular and polymerized actin and quantitatively compare their values to those theoretically predicted earlier. Actin filaments have been demonstrated to act as conducting bionanowires, forming a signaling network capable of transmitting ionic waves in cells. We performed conductivity measurements for different concentrations of actin, considering both unpolymerized and polymerized actin to identify potential differences in their electrical properties. These measurements revealed two relevant characteristics: first, the polymerized actin, arranged in filaments, has a lower impedance than its globular counterpart; second, an increase in the actin concentration leads to higher conductivities. Furthermore, from the data collected, we developed a quantitative model to represent the electrical properties of actin in a buffer solution. We hypothesize that actin filaments can be modeled as electrical resistor-inductor-capacitor (RLC) circuits, where the resistive contribution is due to the viscous ion flows along the filaments; the inductive contribution is due to the solenoidal flows along and around the helix-shaped filament and the capacitive contribution is due to the counterion layer formed around each negatively charged filament., Competing Interests: The authors declare no competing interests.

Combretastatins isolated from the Combretum caffrum tree belong to a group of closely related stilbenes. They are colchicine binding site inhibitors which disrupt the polymerization process of microtubules in tubulins, causing mitotic arrest. In vitro and in vivo studies have proven that some combretastatins exhibit antitumor properties, and among them, combretastatin A-4 is the most active mitotic inhibitor. In this study, a series of novel combretastatin A-4 analogs containing carboxylic acid, ester, and amide moieties were synthesized and their cytotoxic activity against six tumor cell lines was determined using sulforhodamine B assay. For the most cytotoxic compounds ( 8 and 20 ), further studies were performed. These compounds were shown to induce G0/G1 cell cycle arrest in MDA and A549 cells, in a concentration-dependent manner. Moreover, in vitro tubulin polymerization assays showed that both compounds are tubulin polymerization enhancers. Additionally, computational analysis of the binding modes and binding energies of the compounds with respect to the key human tubulin isotypes was performed. We have obtained a satisfactory correlation of the binding energies with the IC 50 values when weighted averages of the binding energies accounting for the abundance of tubulin isotypes in specific cancer cell lines were computed., Competing Interests: The authors declare no conflicts of interest.

Virtually all cells use energy-driven, ion-specific membrane pumps to maintain large transmembrane gradients of Na + , K + , Cl - , Mg ++ , and Ca ++ , but the corresponding evolutionary benefit remains unclear. We propose that these gradients enable a dynamic and versatile biological system that acquires, analyzes, and responds to environmental information. We hypothesize that environmental signals are transmitted into the cell by ion fluxes along pre-existing gradients through gated ion-specific membrane channels. The consequent changes in cytoplasmic ion concentration can generate a local response or orchestrate global/regional cellular dynamics through wire-like ion fluxes along pre-existing and self-assembling cytoskeleton to engage the endoplasmic reticulum, mitochondria, and nucleus., Competing Interests: D.N., J.A.T, and R.A.G. have no conflicts of interest. I.E.N is on the scientific advisory of Endectra, LLC. And cofounder of iRAY technologies LLC, act as deputy editor for the journal of Medical Physics, co-chief editor of BJR AI, and receives funding from NIH and DoD., (© 2024 The Authors.)

We investigate Quantum Electrodynamics corresponding to the holographic brain theory introduced by Pribram to describe memory in the human brain. First, we derive a super-radiance solution in Quantum Electrodynamics with non-relativistic charged bosons (a model of molecular conformational states of water) for coherent light sources of holograms. Next, we estimate memory capacity of a brain neocortex, and adopt binary holograms to manipulate optical information. Finally, we introduce a control theory to manipulate holograms involving biological water's molecular conformational states. We show how a desired waveform in holography is achieved in a hierarchical model using numerical simulations.

Protein-protein interaction occurs on surface patches with some degree of complementary geometric and chemical features. Building on this understanding, this study endeavors to characterize the spike protein of the SARS-CoV-2 virus at the morphological and geometrical levels in its Alpha, Delta, and Omicron variants. In particular, the affinity between different SARS-CoV-2 spike proteins and the ACE2 receptor present on the membrane of the human respiratory system cells is investigated. To achieve an adequate degree of geometrical accuracy, the 3D depth maps of the proteins in exam are filtered by developing an ad-hoc convolutional filter with a kernel implemented as a sphere of varying radius, simulating a ball rolling on the surface (similar to the 'rolling ball' filter). This ball ideally models a hypothetical molecule that could interface with the protein and is inspired by the geometric approach to macromolecule-ligand interactions proposed by Kuntz et al. in 1982. The aim is to mitigate the imperfections and to obtain a smoother surface that could be studied from a geometrical perspective for binding purposes. A set of geometric descriptors, borrowed from the 3D face analysis context is then mapped point-by-point onto protein depth maps. Following a feature extraction phase inspired by Histogram of Oriented Gradients and Local Binary Patterns, the final histogram features are used as input for a Support Vector Machine classifier to automatically classify the proteins according to their surface affinity, where a similarity in shape is observed between ACE2 and the spike protein of the SARS-CoV-2 Omicron variant. Finally, Root Mean Square Error analysis is used to quantify the geometrical affinity between the ACE2 receptor and the respective Receptor Binding Domains of the three SARS-CoV-2 variants, culminating in a geometrical explanation for the higher contagiousness of Omicron relative to the other variants under study., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 Elsevier Inc. All rights reserved.)

Volatile anesthetics (VAs) are medicinal chemistry compounds commonly used to enable surgical procedures for patients who undergo painful treatments and can be partially or fully sedated, remaining in an unconscious state during the operation. The specific molecular mechanism of anesthesia is still an open issue, but scientific evidence supports the hypothesis of the involvement of both putative hydrophobic cavities in membrane receptors as binding pockets and interactions between anesthetics and cytoplasmic proteins. Previous studies demonstrated the binding of VAs to tubulin. Since actin is the other major component of the cytoskeleton, this study involves an investigation of its interactions with four major anesthetics: halothane, isoflurane, sevoflurane, and desflurane. Molecular docking was implemented using the Molecular Operating Environment (MOE) software (version 2022.02) and applied to a G-actin monomer, extrapolating the relative binding affinities and root-mean-square deviation (RMSD) values. A comparison with the F-actin was also made to assess if the generally accepted idea about the enhanced F-to-G-actin transformation during anesthesia is warranted. Overall, our results confirm the solvent-like behavior of anesthetics, as evidenced by Van der Waals interactions as well as the relevant hydrogen bonds formed in the case of isoflurane and sevoflurane. Also, a comparison of the interactions of anesthetics with tubulin was made. Finally, the short- and long-term effects of anesthetics are discussed for their possible impact on the occurrence of mental disorders.

Due to the stimulation of neuronal membrane dipoles by action potentials, under suitable conditions coherent dipole oscillations can be formed. We argue that these dipole oscillations satisfy the weak Bose-Einstein condensate criteria of the Froehlich model of biological coherence. They can subsequently generate electromagnetic fields (EMFs) propagating in the inter-neuronal space. When neighboring neurons fire synchronously, EMFs can create interference patterns and hence form holographic images containing analog information about the sensory inputs that trigger neuronal activity. The mirror pattern projected by EMFs inside the neuron can encode information in the neuronal cytoskeleton. We outline an experimental verification of our hypothesis and its consequences for anesthesia, neurodegenerative diseases, and psychiatric states., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Cavaglià, Deriu and Tuszynski.)

Fear conditioning constitutes the best and most reproducible paradigm to study the neurobiological mechanisms underlying emotions. On the other hand, studies on the synaptic plasticity phenomena underlying fear conditioning present neural circuits enforcing this learning pattern related to post-traumatic stress disorder (PTSD). Notably, in both humans and the rodent model, fear conditioning and context rely on dependent neurocircuitry in the amygdala and prefrontal cortex, cingulate gyrus, and hippocampus. In this review, an overview of the role that classical neurotransmitters play in the contextual conditioning model of fear, and therefore in PTSD, was reported.

We describe a strategy for the development of a rational approach of neoplastic disease therapy based on the demonstration that scale-free networks are susceptible to specific attacks directed against its connective hubs. This strategy involves the (i) selection of up-regulated hubs of connectivity in the tumors interactome, (ii) drug repurposing of these hubs, (iii) RNA silencing of non-druggable hubs, (iv) in vitro hub validation, (v) tumor-on-a-chip, (vi) in vivo validation, and (vii) clinical trial. Hubs are protein targets that are assessed as targets for rational therapy of cancer in the context of personalized oncology. We confirmed the existence of a negative correlation between malignant cell aggressivity and the target number needed for specific drugs or RNA interference (RNAi) to maximize the benefit to the patient's overall survival. Interestingly, we found that some additional proteins not generally targeted by drug treatments might justify the addition of inhibitors designed against them in order to improve therapeutic outcomes. However, many proteins are not druggable, or the available pharmacopeia for these targets is limited, which justifies a therapy based on encapsulated RNAi.

Cancer is amenable to low-cost treatments, given that it has a significant metabolic component, which can be affected through diet and lifestyle change at minimal cost. The Warburg hypothesis states that cancer cells have an altered cell metabolism towards anaerobic glycolysis. Given this metabolic reprogramming in cancer cells, it is possible to target cancers metabolically by depriving them of glucose. In addition to dietary and lifestyle modifications which work on tumors metabolically, there are a panoply of nutritional supplements and repurposed drugs associated with cancer prevention and better treatment outcomes. These interventions and their evidentiary basis are covered in the latter half of this review to guide future cancer treatment.

One of the most important aspects of successful cancer therapy is the identification of a target protein for inhibition interaction. Conventionally, this consists of screening a panel of genes to assess which is mutated and then developing a small molecule to inhibit the interaction of two proteins or to simply inhibit a specific protein from all interactions. In previous work, we have proposed computational methods that analyze protein-protein networks using both topological approaches and thermodynamic quantification provided by Gibbs free energy. In order to make these approaches both easier to implement and free of arbitrary topological filtration criteria, in the present paper, we propose a modification of the topological-thermodynamic analysis, which focuses on the selection of the most thermodynamically stable proteins and their subnetwork interaction partners with the highest expression levels. We illustrate the implementation of the new approach with two specific cases, glioblastoma (glioma brain tumors) and chronic lymphatic leukoma (CLL), based on the publicly available patient-derived datasets. We also discuss how this can be used in clinical practice in connection with the availability of approved and investigational drugs.

Hypercoagulability and formation of extensive and difficult-to-lyse microclots are a hallmark of both acute COVID-19 and long COVID. Fibrinogen, when converted to fibrin, is responsible for clot formation, but abnormal structural and mechanical clot properties can lead to pathologic thrombosis. Recent experimental evidence suggests that the spike protein (SP) from severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) may directly bind to the blood coagulation factor fibrinogen and induce structurally abnormal blood clots with heightened proinflammatory activity. Accordingly, in this study, we used molecular docking and molecular dynamics simulations to explore the potential activity of the antiparasitic drug ivermectin (IVM) to prevent the binding of the SARS-CoV-2 SP to fibrinogen and reduce the occurrence of microclots. Our computational results indicate that IVM may bind with high affinity to multiple sites on the fibrinogen peptide, with binding more likely in the central, E region, and in the coiled-coil region, as opposed to the globular D region. Taken together, our in silico results suggest that IVM may interfere with SP-fibrinogen binding and, potentially, decrease the formation of fibrin clots resistant to degradation. Additional in vitro studies are warranted to validate whether IVM binding to fibrinogen is sufficiently stable to prevent interaction with the SP, and potentially reduce its thrombo-inflammatory effect in vivo.

Background: Patients with advanced hepatocellular carcinoma (HCC) and poor liver function lack effective systemic therapies. Low-energy electromagnetic fields (EMFs) can influence cell biological processes via non-thermal effects and may represent a new treatment option., Methods: This single-site feasibility trial enrolled patients with advanced HCC, Child-Pugh A and B, Eastern Cooperative Oncology Group 0-2. Patients underwent 90-min amplitude-modulated EMF exposure procedures every 2-4 weeks, using the AutEMdev (Autem Therapeutics). Patients could also receive standard care. The primary endpoints were safety and the identification of hemodynamic variability patterns. Exploratory endpoints included health-related quality of life (HRQoL), overall survival (OS). and objective response rate (ORR) using RECIST v1.1., Results: Sixty-six patients with advanced HCC received 539 AutEMdev procedures (median follow-up, 30 months). No serious adverse events occurred during procedures. Self-limiting grade 1 somnolence occurred in 78.7% of patients. Hemodynamic variability during EMF exposure was associated with specific amplitude-modulation frequencies. HRQoL was maintained or improved among patients remaining on treatment. Median OS was 11.3 months (95% confidence interval [CI]: 6.0, 16.6) overall (16.0 months [95% CI: 4.4, 27.6] and 12.0 months [6.4, 17.6] for combination therapy and monotherapy, respectively). ORR was 24.3% (32% and 17% for combination therapy and monotherapy, respectively)., Conclusion: AutEMdev EMF exposure has an excellent safety profile in patients with advanced HCC. Hemodynamic alterations at personalized frequencies may represent a surrogate of anti-tumor efficacy. NCT01686412., (© 2023 Autem Therapeutics. Cancer Medicine published by John Wiley & Sons Ltd.)

Given its critical role in cell mitosis, the tubulin γ chain represents a viable chemotherapeutic target to solve the specificity issues associated with targeting α and β tubulin. Since γ tubulin is overexpressed in glioblastoma multiforme (GBM) and some breast lesions, the glaziovianin A derivative gatastatin, presented as a γ-tubulin-specific inhibitor, could yield a successful therapeutic strategy. The present work aims to identify the binding sites and modes of gatastatin and its derivatives through molecular-docking simulations. Computational binding free energy predictions were compared to experimental microscale thermophoresis assay results. The computational simulations did not reveal a strong preference toward γ tubulin, suggesting that further derivatization may be needed to increase its specificity.

Microtubules are cylindrical protein polymers assembled in the cytoplasm of all eukaryotic cells by polymerization of aβ tubulin dimers, which are involved in cell division, migration, signaling, and intracellular traffic. These functions make them essential in the proliferation of cancerous cells and metastases. Tubulin has been the molecular target of many anticancer drugs because of its crucial role in the cell proliferation process. By developing drug resistance, tumor cells severely limit the successful outcomes of cancer chemotherapy. Hence, overcoming drug resistance motivates the design of new anticancer therapeutics. Here, we retrieve short peptides obtained from the data repository of antimicrobial peptides (DRAMP) and report on the computational screening of their predicted tertiary structures for the ability to inhibit tubulin polymerization using multiple combinatorial docking programs, namely PATCHDOCK, FIREDOCK, and ClusPro. The interaction visualizations show that all the best peptides from the docking analysis bind to the interface residues of the tubulin isoforms αβl, αβll, αβlll, and αβlV, respectively. The docking studies were further confirmed by a molecular dynamics simulation, in which the computed root-mean-square deviation (RMSD), and root-mean-square fluctuation (RMSF), verified the stable nature of the peptide-tubulin complexes. Physiochemical toxicity and allergenicity studies were also performed. This present study suggests that these identified anticancer peptide molecules might destabilize the tubulin polymerization process and hence can be suitable candidates for novel drug development. It is concluded that wet-lab experiments are needed to validate these findings.

Over the past two decades, it was discovered that introducing synthetic small interfering RNAs (siRNAs) into the cytoplasm facilitates effective gene-targeted silencing. This compromises gene expression and regulation by repressing transcription or stimulating sequence-specific RNA degradation. Substantial investments in developing RNA therapeutics for disease prevention and treatment have been made. We discuss the application to proprotein convertase subtilisin/kexin type 9 (PCSK9), which binds to and degrades the low-density lipoprotein cholesterol (LDL-C) receptor, interrupting the process of LDL-C uptake into hepatocytes. PCSK9 loss-of-function modifications show significant clinical importance by causing dominant hypocholesterolemia and lessening the risk of cardiovascular disease (CVD). Monoclonal antibodies and small interfering RNA (siRNA) drugs targeting PCSK9 are a significant new option for managing lipid disorders and improving CVD outcomes. In general, monoclonal antibodies are restricted to binding with cell surface receptors or circulating proteins. Similarly, overcoming the intracellular and extracellular defenses that prevent exogenous RNA from entering cells must be achieved for the clinical application of siRNAs. N -acetylgalactosamine (GalNAc) conjugates are a simple solution to the siRNA delivery problem that is especially suitable for treating a broad spectrum of diseases involving liver-expressed genes. Inclisiran is a GalNAc-conjugated siRNA molecule that inhibits the translation of PCSK9. The administration is only required every 3 to 6 months, which is a significant improvement over monoclonal antibodies for PCSK9. This review provides an overview of siRNA therapeutics with a focus on detailed profiles of inclisiran, mainly its delivery strategies. We discuss the mechanisms of action, its status in clinical trials, and its prospects.

Exposure to Low-Energy Amplitude-Modulated Radiofrequency Electromagnetic Fields (LEAMRFEMF) represents a new treatment option for patients with advanced hepatocellular carcinoma (AHCC). We focus on two medical devices that modulate the amplitude of a 27.12 MHz carrier wave to generate envelope waves in the low Hz to kHz range. Each provides systemic exposure to LEAMRFEMF via an intrabuccal antenna. This technology differs from so-called Tumour Treating Fields because it uses different frequency ranges, uses electromagnetic rather than electric fields, and delivers energy systemically rather than locally. The AutemDev also deploys patient-specific frequencies. LEAMRFEMF devices use 100-fold less power than mobile phones and have no thermal effects on tissue. Tumour type-specific or patient-specific treatment frequencies can be derived by measuring haemodynamic changes induced by exposure to LEAMRFEMF. These specific frequencies inhibited growth of human cancer cell lines in vitro and in mouse xenograft models. In uncontrolled prospective clinical trials in patients with AHCC, minorities of patients experienced complete or partial tumour responses. Pooled comparisons showed enhanced overall survival in treated patients compared to historical controls. Mild transient somnolence was the only notable treatment-related adverse event. We hypothesize that intracellular oscillations of charged macromolecules and ion flows couple resonantly with LEAMRFEMF. This resonant coupling appears to disrupt cell division and subcellular trafficking of mitochondria. We provide an estimate of the contribution of the electromagnetic effects to the overall energy balance of an exposed cell by calculating the power delivered to the cell, and the energy dissipated through the cell due to EMF induction of ionic flows along microtubules. We then compare this with total cellular metabolic energy production and conclude that energy delivered by LEAMRFEMF may provide a beneficial shift in cancer cell metabolism away from aberrant glycolysis. Further clinical research may confirm that LEAMRFEMF has therapeutic value in AHCC., Competing Interests: JAT and FC were employed by Autem Therapeutics. JAT and FC have stock and honoraria from Autem Medical and patents relating to the AutEMdev device., (© 2022 Tuszynski and Costa.)

Heterocyclic compounds are a class of compounds of natural origin with favorable properties and hence have major pharmaceutical significance. They have an exceptional adroitness favoring their use as diverse smart biomimetics, in addition to possessing an active pharmacophore in a complex structure. This has made them an indispensable motif in the drug discovery field. Heterocyclic compounds are usually classified according to the ring size, type, and the number of heteroatoms present in the ring. Among different heterocyclic ring systems, nitrogen heterocyclic compounds are more abundant in nature. They also have considerable pharmacological significance. This review highlights recent pioneering studies in the biological assessment of nitrogen-containing compounds, namely: triazoles, tetrazoles, imidazole/benzimidazoles, pyrimidines, and quinolines. It explores publications between April 2020 and February 2022 and will benefit researchers in medicinal chemistry and pharmacology. The present work is organized based on the size of the heterocyclic ring.

Scoulerine is a natural compound that is known to bind to tubulin and has anti-mitotic properties demonstrated in various cancer cells. Its molecular mode of action has not been precisely known. In this work, we perform computational prediction and experimental validation of the mode of action of scoulerine. Based on the existing data in the Protein Data Bank (PDB) and using homology modeling, we create human tubulin structures corresponding to both free tubulin dimers and tubulin in a microtubule. We then perform docking of the optimized structure of scoulerine and find the highest affinity binding sites located in both the free tubulin and in a microtubule. We conclude that binding in the vicinity of the colchicine binding site and near the laulimalide binding site are the most likely locations for scoulerine interacting with tubulin. Thermophoresis assays using scoulerine and tubulin in both free and polymerized form confirm these computational predictions. We conclude that scoulerine exhibits a unique property of a dual mode of action with both microtubule stabilization and tubulin polymerization inhibition, both of which have similar affinity values.

The SARS-CoV-2 virus invades and replicates within host cells by "hijacking" biomolecular machinery, gaining control of the microtubule cytoskeleton. After attaching to membrane receptors and entering cells, the SARS-CoV-2 virus co-opts the dynamic intra-cellular cytoskeletal network of microtubules, actin, and the microtubule-organizing center, enabling three factors that lead to clinical pathology: (1) viral load due to intra-cellular trafficking, (2) cell-to-cell spread by filopodia, and (3) immune dysfunction, ranging from hyper-inflammatory cytokine storm to ineffective or absent response. These factors all depend directly on microtubules and the microtubule-organizing center, as do cell functions such as mitosis and immune cell movement. Here we consider how the SARS-CoV-2 virus may "hijack" cytoskeletal functions by docking inside the microtubule-organizing center's centriole "barrels", enabling certain interactions between the virus's positively charged spike ("S") proteins and negatively charged C-termini of the microtubules that the centriole comprises, somewhat like fingers on a keyboard. This points to the potential benefit of therapies aimed not directly at the virus but at the microtubules and microtubule-organizing center of the host cell on which the virus depends. These therapies could range from anti-microtubule drugs to low-intensity ultrasound (megahertz mechanical vibrations) externally applied to the vagus nerve at the neck and/or to the spleen (since both are involved in mediating inflammatory response). Given that ultrasound imaging machines suitable for vagal/splenic ultrasound are available for clinical trials in every hospital, we recommend an alternative therapeutic approach for COVID-19 based on addressing and normalizing the host cell microtubules and microtubule-organizing centers co-opted by the SARS-CoV-2 virus.

Pyrazoles are five-membered heterocyclic compounds that contain nitrogen. They are an important class of compounds for drug development; thus, they have attracted much attention. In the meantime, pyrazole derivatives have been synthesized as target structures and have demonstrated numerous biological activities such as antituberculosis, antimicrobial, antifungal, and anti-inflammatory. This review summarizes the results of published research on pyrazole derivatives synthesis and biological activities. The published research works on pyrazole derivatives synthesis and biological activities between January 2018 and December 2021 were retrieved from the Scopus database and reviewed accordingly.

The existence of an exclusion zone in which particles of a colloidal suspension in water are repelled from hydrophilic surfaces has been experimentally demonstrated in numerous studies, especially in the case of Nafion surfaces. Various explanations have been proposed for the origin of this phenomenon, which is not completely understood yet. In particular, the existence of a fourth phase of water has been proposed by G. Pollack and if this theory is proven correct, its implications on our understanding of the properties of water, especially in biological systems, would be profound and could give rise to new medical therapies. Here, a simple approach based on the linearized Poisson-Boltzmann equation is developed in order to study the repulsive forces mediated by ordered water and involving the following interacting biomolecules: 1) microtubule and a tubulin dimer, 2) two tubulin dimers and 3) a tubulin sheet and a tubulin dimer. The choice of microtubules in this study is motivated because they could be a good candidate for the generation of an exclusion zone in the cell and these models could be a starting point for detailed experimental investigations of this phenomenon., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Chierici, Dogariu and Tuszynski.)

Microtubules are cylindrical protein polymers formed from α β -tubulin heterodimers in the cytoplasm of eukaryotic cells. Microtubule disturbance may cause cell cycle arrest in the G2/M phase, and anomalous mitotic spindles will form. Microtubules are an important target for cancer drug action because of their critical role in mitosis. Several microtubule-targeting agents with vast therapeutic advantages have been developed, but they often lead to multidrug resistance and adverse side effects. Thus, single-target therapy has drawbacks in the effective control of tubulin polymerization. Molecular hybridization, based on the amalgamation of two or more pharmacophores of bioactive conjugates to engender a single molecular structure with enhanced pharmacokinetics and biological activity, compared to their parent molecules, has recently become a promising approach in drug development. The practical application of combined active scaffolds targeting tubulin polymerization inhibitors has been corroborated in the past few years. Meanwhile, different designs and syntheses of novel anti-tubulin hybrids have been broadly studied, illustrated, and detailed in the literature. This review describes various molecular hybrids with their reported structural-activity relationships (SARs) where it is possible in an effort to generate efficacious tubulin polymerization inhibitors. The aim is to create a platform on which new active scaffolds can be modeled for improved tubulin polymerization inhibitory potency and hence, the development of new therapeutic agents against cancer.

Inhibition of DNA repair enzymes is an attractive target for increasing the efficacy of DNA damaging chemotherapies. The ERCC1-XPF heterodimer is a key endonuclease in numerous single and double strand break repair processes, and inhibition of the heterodimerization has previously been shown to sensitize cancer cells to DNA damage. In this work, the previously reported ERCC1-XPF inhibitor 4 was used as the starting point for an in silico study of further modifications of the piperazine side-chain. A selection of the best scoring hits from the in silico screen were synthesized using a late stage functionalization strategy which should allow for further iterations of this class of inhibitors to be readily synthesized. Of the synthesized compounds, compound 6 performed the best in the in vitro fluorescence based endonuclease assay. The success of compound 6 in inhibiting ERCC1-XPF endonuclease activity in vitro translated well to cell-based assays investigating the inhibition of nucleotide excision repair and disruption of heterodimerization. Subsequently compound 6 was shown to sensitize HCT-116 cancer cells to treatment with UVC, cyclophosphamide, and ionizing radiation. This work serves as an important step towards the synergistic use of DNA repair inhibitors with chemotherapeutic drugs., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Weilbeer, Jay, Donnelly, Gentile, Karimi-Busheri, Yang, Mani, Yu, Elmenoufy, Barakat, Tuszynski, Weinfeld and West.)

Very few drugs in clinical practice feature the chemical diversity, narrow therapeutic window, unique route of administration, and reversible cognitive effects of volatile anesthetics. The correlation between their hydrophobicity and their potency and the increasing amount of evidence suggesting that anesthetics exert their action on transmembrane proteins, justifies the investigation of their effects on phospholipid bilayers at the molecular level, given the strong functional and structural link between transmembrane proteins and the surrounding lipid matrix. Molecular dynamics simulations of a model lipid bilayer in the presence of ethylene, desflurane, methoxyflurane, and the nonimmobilizer 1,2-dichlorohexafluorocyclobutane (also called F6 or 2N) at different concentrations highlight the structural consequences of VA partitioning in the lipid phase, with a decrease of lipid order and bilayer thickness, an increase in overall lipid lateral mobility and area-per-lipid, and a marked reduction in the mechanical stiffness of the membrane, that strongly correlates with the compounds' hydrophobicity., Competing Interests: The authors declare no competing interests., (© 2022 The Author(s).)

The assembly/disassembly of biological macromolecules plays an important role in their biological functionalities. Although the dynamics of tubulin polymers and their super-assembly into microtubule structures is critical for many cellular processes, details of their cyclical polymerization/depolymerization are not fully understood. Here, we use a specially designed light scattering technique to continuously examine the effects of temperature cycling on the process of microtubule assembly/disassembly. We observe a thermal hysteresis loop during tubulin assembly/disassembly, consistently with earlier reports on the coexistence of tubulin and microtubules as a phase transition. In a cyclical process, the structural hysteresis has a kinetic component that depends on the rate of temperature change but also an intrinsic thermodynamic component that depends on the protein topology, possibly related to irreversible processes. Analyzing the evolution of such thermal hysteresis loops over successive cycles, we found that the assembly/disassembly ceases after some time, which is indicative of protein aging leading to its inability to self-assemble after a finite number of temperature cycles. The emergence of assembly-incompetent tubulin could have major consequences for human pathologies related to microtubules, including aging, neurodegenerative diseases and cancer., Competing Interests: The authors declare no conflict of interest., (© 2021 Published by Elsevier B.V.)

The underlying mechanism determining the size of a particular cell is one of the fundamental unknowns in cell biology. Here, using a new approach that could be used for most of unicellular species, we show that the protein synthesis and cell size are interconnected biophysically and that protein synthesis may be the chief mechanism in establishing size limitations of unicellular organisms. This result is obtained based on the free energy balance equation of protein synthesis and the second law of thermodynamics. Our calculations show that protein synthesis involves a considerable amount of entropy reduction due to polymerization of amino acids depending on the cytoplasmic volume of the cell. The amount of entropy reduction will increase with cell growth and eventually makes the free energy variations of the protein synthesis positive (that is, forbidden thermodynamically). Within the limits of the second law of thermodynamics we propose a framework to estimate the optimal cell size at division.

Competing Interests: The authors declare no conflict of interest. MC, MD and JAT are serving as Guest Editors of this journal. We declare that they had no involvement in the peer review of this article and has no access to information regarding its peer review. Full responsibility for the editorial process for this article was delegated to RRP.

Background: The emergence and rapid spread of SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) in thelate 2019 has caused a devastating global pandemic of the severe pneumonia-like disease coronavirus disease 2019 (COVID-19). Although vaccines have been and are being developed, they are not accessible to everyone and not everyone can receive these vaccines. Also, it typically takes more than 10 years until a new therapeutic agent is approved for usage. Therefore, repurposing of known drugs can lend itself well as a key approach for significantly expediting the development of new therapies for COVID-19., Methods: We have incorporated machine learning-based computational tools and in silico models into the drug discovery process to predict Adsorption, Distribution, Metabolism, Excretion, and Toxicity (ADMET) profiles of 90 potential drugs for COVID-19 treatment identified from two independent studies mainly with the purpose of mitigating late-phase failures because of inferior pharmacokinetics and toxicity., Results: Here, we summarize the cardiotoxicity and general toxicity profiles of 90 potential drugs for COVID-19 treatment and outline the risks of repurposing and propose a stratification of patients accordingly. We shortlist a total of five compounds based on their non-toxic properties., Conclusion: In summary, this manuscript aims to provide a potentially useful source of essential knowledge on toxicity assessment of 90 compounds for healthcare practitioners and researchers to find off-label alternatives for the treatment for COVID-19. The majority of the molecules discussed in this manuscript have already moved into clinical trials and thus their known pharmacological and human safety profiles are expected to facilitate a fast track preclinical and clinical assessment for treating COVID-19., (© 2021. The Author(s).)

Galectin-3 is a carbohydrate-binding protein and the most studied member of the galectin family. It regulates several functions throughout the body, among which are inflammation and post-injury remodelling. Recent studies have highlighted the similarity between Galectin-3's carbohydrate recognition domain and the so-called "galectin fold" present on the N-terminal domain of the S1 sub-unit of the SARS-CoV-2 spike protein. Sialic acids binding to the N-terminal domain of the Spike protein are known to be crucial for viral entry into humans, and the role of Galectin-3 as a mediator of lung fibrosis has long been the object of study since its levels have been found to be abnormally high in alveolar macrophages following lung injury. In this context, the discovery of a double inhibitor may both prevent viral entry and reduce post-infection pulmonary fibrosis. In this study, we use a database of 56 compounds, among which 37 have known experimental affinity with Galectin-3. We carry out virtual screening of this database with respect to Galectin-3 and Spike protein. Several ligands are found to exhibit promising binding affinity and interaction with the Spike protein's N-terminal domain as well as with Galectin-3. This finding strongly suggests that existing Galectin-3 inhibitors possess dual-binding capabilities to disrupt Spike-ACE2 interactions. Herein we identify the most promising inhibitors of Galectin-3 and Spike proteins, of which five emerge as potential dual effective inhibitors. Our preliminary results warrant further in vitro and in vivo testing of these putative inhibitors against SARS-CoV-2 with the hope of being able to halt the spread of the virus in the future.

The biological effects of terahertz (THz) radiation have been observed across multiple levels of biological organization, however the sub-cellular mechanisms underlying the phenotypic changes remain to be elucidated. Filamentous protein complexes such as microtubules are essential cytoskeletal structures that regulate diverse biological functions, and these may be an important target for THz interactions underlying THz-induced effects observed at the cellular or tissue level. Here, we show disassembly of microtubules within minutes of exposure to extended trains of intense, picosecond-duration THz pulses. Further, the rate of disassembly depends on THz intensity and spectral content. As inhibition of microtubule dynamics is a mechanism of clinically-utilized anti-cancer agents, disruption of microtubule networks may indicate a potential therapeutic mechanism of intense THz pulses., Competing Interests: The authors declare no conflicts of interest., (© 2021 Optical Society of America under the terms of the OSA Open Access Publishing Agreement.)

The Janus Kinase signalling pathway is implicated in the pathogenesis of immune-related diseases. The potency of small-molecule Janus Kinase inhibitors in the treatment of inflammatory diseases demonstrates that this pathway can be successfully targeted for therapeutic purposes. The outstanding relevant questions concerning drugs' efficacy and toxicity challenge the research to enhance the selectivity of these drugs. The promising results of computational techniques, such as Molecular Dynamics and Molecular Docking, coupled with experimental studies, can improve the understanding of the molecular mechanism of Janus Kinase pathway and thus enable the rational design of new more selective inhibitor molecules., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2020. Published by Elsevier Inc.)

Microtubules are highly negatively charged proteins which have been shown to behave as bio-nanowires capable of conducting ionic currents. The electrical characteristics of microtubules are highly complicated and have been the subject of previous work; however, the impact of the ionic concentration of the buffer solution on microtubule electrical properties has often been overlooked. In this work we use the non-linear Poisson Boltzmann equation, modified to account for a variable permittivity and a Stern Layer, to calculate counterion concentration profiles as a function of the ionic concentration of the buffer. We find that for low-concentration buffers ([KCl] from 10 μ M to 10 mM ) the counterion concentration is largely independent of the buffer's ionic concentration, but for physiological-concentration buffers ([KCl] from 100 to 500 mM ) the counterion concentration varies dramatically with changes in the buffer's ionic concentration. We then calculate the conductivity of microtubule-counterion complexes, which are found to be more conductive than the buffer when the buffer's ionic concentrations is less than ≈100 mM and less conductive otherwise. These results demonstrate the importance of accounting for the ionic concentration of the buffer when analyzing microtubule electrical properties both under laboratory and physiological conditions. We conclude by calculating the basic electrical parameters of microtubules over a range of ionic buffer concentrations applicable to nanodevice and medical applications., Competing Interests: The authors declare that this study received funding from Novocure Inc. The funder was not involved in the study design, collection, analysis, interpretation of data, the writing of this article or the decision to submit it for publication., (Copyright © 2021 Eakins, Patel, Kalra, Rezania, Shankar and Tuszynski.)

The coronavirus causing the COVID-19 pandemic, SARS-CoV-2, uses -1 programmed ribosomal frameshifting (-1 PRF) to control the relative expression of viral proteins. As modulating -1 PRF can inhibit viral replication, the RNA pseudoknot stimulating -1 PRF may be a fruitful target for therapeutics treating COVID-19. We modeled the unusual 3-stem structure of the stimulatory pseudoknot of SARS-CoV-2 computationally, using multiple blind structural prediction tools followed by μs-long molecular dynamics simulations. The results were compared for consistency with nuclease-protection assays and single-molecule force spectroscopy measurements of the SARS-CoV-1 pseudoknot, to determine the most likely conformations. We found several possible conformations for the SARS-CoV-2 pseudoknot, all having an extended stem 3 but with different packing of stems 1 and 2. Several conformations featured rarely-seen threading of a single strand through junctions formed between two helices. These structural models may help interpret future experiments and support efforts to discover ligands inhibiting -1 PRF in SARS-CoV-2., Competing Interests: The authors have declared that no competing interests exist.

This paper presents a historical overview of quantum physics methodology's development and application to various science fields beyond physics, especially biology and consciousness. Following a successful interpretation of several early 20t⁢h century experiments, quantum physics gradually provided a conceptual framework for molecular bonds via quantum chemistry. In recent years individual biological phenomena such as photosynthesis and bird navigation have been experimentally and theoretically analyzed using quantum methods, building conceptual foundations for quantum physics' entry into biology. Quantum concepts have also been recently employed to explain physiology's allometric scaling laws by introducing quantum metabolism theory. In the second part of this work, we discuss how quantum physics may also be pivotal to our understanding of consciousness, which has been touted by some researchers as the last frontier of modern science. Others believe that consciousness does not belong within the realm of science at all. Several hypotheses, especially the Orch OR theory, have been suggested over the past two decades to introduce a scientific basis to consciousness theory. We discuss the merits and potential extensions of these approaches., Competing Interests: The author declares no conflict of interest., (© 2020 Tuszynski Published by IMR Press.)

Cytochrome P450 monooxygenases (CYPs/P450s) are well known for their role in organisms' primary and secondary metabolism. Among 20 P450s of the tuberculosis-causing Mycobacterium tuberculosis H37Rv, CYP128A1 is particularly important owing to its involvement in synthesizing electron transport molecules such as menaquinone-9 (MK9). This study employs different in silico approaches to understand CYP128 P450 family's distribution and structural aspects. Genome data-mining of 4250 mycobacterial species has revealed the presence of 2674 CYP128 P450s in 2646 mycobacterial species belonging to six different categories. Contrast features were observed in the CYP128 gene distribution, subfamily patterns, and characteristics of the secondary metabolite biosynthetic gene cluster (BGCs) between M. tuberculosis complex (MTBC) and other mycobacterial category species. In all MTBC species (except one) CYP128 P450s belong to subfamily A, whereas subfamily B is predominant in another four mycobacterial category species. Of CYP128 P450s, 78% was a part of BGCs with CYP124A1 , or together with CYP124A1 and CYP121A1 . The CYP128 family ranked fifth in the conservation ranking. Unique amino acid patterns are present at the EXXR and CXG motifs. Molecular dynamic simulation studies indicate that the CYP128A1 bind to MK9 with the highest affinity compared to the azole drugs analyzed. This study provides comprehensive comparative analysis and structural insights of CYP128A1 in M. tuberculosis .

Colchicine belongs to a large group of microtubule polymerization inhibitors. Although the anti-cancer activity of colchicine and its derivatives has been established, none of them has found commercial application in cancer treatment due to side effects. Therefore, we designed and synthesized a series of six triple-modified 4-chlorothiocolchicine analogues with amide moieties and one urea derivative. These novel derivatives were tested against several different cancer cell lines (A549, MCF-7, LoVo, LoVo/DX) and primary acute lymphoblastic leukemia (ALL) cells and they showed activity in the nanomolar range. The obtained IC 50 values for novel derivatives were lower than those obtained for unmodified colchicine and common anticancer drugs such as doxorubicin and cisplatin. Further studies of colchicine and selected analogues were undertaken to indicate that they induced apoptotic cell death in ALL-5 cells. We also performed in silico studies to predict binding modes of the 4-chlorothiocolchicine derivatives to different β tubulin isotypes. The results indicate that select triple-modified 4-chlorothiocolchicine derivatives represent highly promising novel cancer chemotherapeutics., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2020 The Authors. Published by Elsevier Inc. All rights reserved.)

We analyzed protein expression data for Lupus patients, which have been obtained from publicly available databases. A combination of systems biology and statistical thermodynamics approaches was used to extract topological properties of the associated protein-protein interaction networks for each of the 291 patients whose samples were used to provide the molecular data. We have concluded that among the many proteins that appear to play critical roles in this pathology, most of them are either ribosomal proteins, ubiquitination pathway proteins or heat shock proteins. We propose some of the proteins identified in this study to be considered for drug targeting., Competing Interests: The authors have declared that no competing interests exist.

Membrane-active agents (MAAs), such as antimicrobial peptides (AMPs) and chemotherapy drugs (CDs), induce ion pores/channels inside lipid bilayer membrane, as confirmed by standard electrophysiology experiments. A novel experimental method is described which detects agents directly at the membrane as confirmed for MAAs, CDs and aptamers. MAAs exhibit characteristic 'charge based' interactions with lipids. Electrostatic (ES) and van der Waals (vdW) contributions to the interaction energies have been estimated using molecular dynamics (MD) simulations. These results are consistent with the screened Coulomb interaction predictions recently developed for lipid bilayer binding of integral AMP channels. Energy- and distance-dependence of MAA-lipid interactions from MD simulations are represented by universal probability functions. A generalized model of MAA-lipid interactions is developed based on the charge and geometrical profiles of the participating lipids and AMPs. The corresponding driving force correlates directly with the stability of MAA-lipid structures as observed in electrophysiology experiments. We conclude that MAAs and similar agents that target lipid membranes exhibit physiological effects mainly due to ES and vdW interactions determined by their charge profiles., (Copyright © 2019 Elsevier Inc. All rights reserved.)

The cardinal role of microtubules in cell mitosis makes them interesting drug targets for many pharmacological treatments, including those against cancer. Moreover, different expression patterns between cell types for several tubulin isotypes represent a great opportunity to improve the selectivity and specificity of the employed drugs and to design novel compounds with higher activity only on cells of interest. In this context, tubulin isotype βIII represents an excellent target for anti-tumoral therapies since it is overexpressed in most cancer cells and correlated with drug resistance. Colchicine is a well-known antimitotic agent, which is able to bind the tubulin dimer and to halt the mitotic process. However, it shows high toxicity also on normal cells and it is not specific for isotype βIII. In this context, the search for colchicine derivatives is a matter of great importance in cancer research. In this study, homology modeling techniques, molecular docking, and molecular dynamics simulations have been employed to characterize the interaction between 55 new promising colchicine derivatives and tubulin isotype βIII. These compounds were screened and ranked based on their binding affinity and conformational stability in the colchicine binding site of tubulin βIII. Results from this study point the attention on an amide of 4-chlorine thiocolchicine. This colchicine-derivative is characterized by a unique mode of interaction with tubulin, compared to all other compounds considered, which is primarily characterized by the involvement of the α-T5 loop, a key player in the colchicine binding site. Information provided by the present study may be particularly important in the rational design of colchicine-derivatives targeting drug resistant cancer phenotypes., (Copyright © 2020 Pallante, Rocca, Klejborowska, Huczynski, Grasso, Tuszynski and Deriu.)

Cancer chemotherapy agents are assessed for their therapeutic utility primarily by their ability to cause apoptosis of cancer cells and their potency is given by an IC50 value. Chemotherapy uses both target-specific and systemic-action drugs and drug combinations to treat cancer. It is important to judiciously choose a drug type, its dosage and schedule for optimized drug selection and administration. Consequently, the precise mathematical formulation of cancer cells' response to chemotherapy may assist in the selection process. In this paper, we propose a mathematical description of the cancer cell response to chemotherapeutic agent exposure based on a time-tested physical model of two-state multiple-component systems near criticality. We describe the Ising model methodology and apply it to a diverse panel of cytotoxic drugs administered against numerous cancer cell lines in a dose-response manner. The analysed dataset was generated by the Netherlands Translational Research Center B.V. (Oncolines). This approach allows for an accurate and consistent analysis of cytotoxic agents' effects on cancer cell lines and reveals the presence or absence of the bystander effect through the interaction constant. By calculating the susceptibility function, we see the value of IC50 coinciding with the peak of this measure of the system's sensitivity to external perturbations., Competing Interests: We declare we have no competing interests., (© 2020 The Authors.)