Philip Winter

Coarse-Grained Molecular Simulation of Diffusion and Reaction Kinetics in a Crowded Virtual Cytoplasm

We present a general-purpose model for biomolecular simulations at the molecular level that incorporates stochasticity, spatial dependence, and volume exclusion, using diffusing and reacting particles with physical dimensions. To validate the model, we first established the formal relationship between the microscopic model parameters (timestep, move length, and reaction probabilities) and the macroscopic coefficients for diffusion and reaction rate. We then compared simulation results with Smoluchowski theory for diffusion-limited irreversible reactions and the best available approximation for diffusion-influenced reversible reactions. To simulate the volumetric effects of a crowded intracellular environment, we created a virtual cytoplasm composed of a heterogeneous population of particles diffusing at rates appropriate to their size. The particle-size distribution was estimated from the relative abundance, mass, and stoichiometries of protein complexes using an experimentally derived proteome catalog from Escherichia coli K12. Simulated diffusion constants exhibited anomalous behavior as a function of time and crowding. Although significant, the volumetric impact of crowding on diffusion cannot fully account for retarded protein mobility in vivo, suggesting that other biophysical factors are at play. The simulated effect of crowding on barnase-barstar dimerization, an experimentally characterized example of a bimolecular association reaction, reveals a biphasic time course, indicating that crowding exerts different effects over different timescales. These observations illustrate that quantitative realism in biosimulation will depend to some extent on mesoscale phenomena that are not currently well understood.

Comparison of SDS‐ and methanol‐assisted protein solubilization and digestion methods for Escherichia coli membrane proteome analysis by 2‐D LC‐MS/MS

Both organic solvent and surfactant have been used for dissolving membrane proteins for shotgun proteomics. In this work, two methods of protein solubilization, namely using 60% methanol or 1% SDS, to dissolve and analyze the inner membrane fraction of an Escherichia coli K12 cell lysate were compared. A total of 358 proteins (1417 unique peptides) from the methanol‐solubilized protein mixture and 299 proteins (892 peptides) from the SDS‐solubilized samplewere identified by using trypsin digestion and 2‐D LC‐ESI MS/MS. It was found that the methanol method detected more hydrophobic peptides, resulting in a greater number of proteins identified, than the SDS method. We found that 159 out of 358 proteins (44%) and 120 out of 299 proteins (40%) detected from the methanol‐ and SDS‐solubilized samples, respectively, are integral membrane proteins. Among the 190 integral membrane proteins 70 were identified exclusively in the methanol‐solubilized sample, 89 were identified by both methods, and only 31 proteins were exclusively identified by the SDS method. It is shown that the integral membrane proteins reflected the theoretical proteome for number of transmembrane helices, length, functional class, and topology, indicating there was no bias in the proteins identified.

Quantitative analysis of the effect of tubulin isotype expression on sensitivity of cancer cell lines to a set of novel colchicine derivatives

BackgroundA 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.ResultsWe 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.ConclusionsThe 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.

Modeling the Yew Tree Tubulin and a Comparison of its Interaction with Paclitaxel to Human Tubulin

ABSTRACTPurposeTo explore possible ways in which yew tree tubulin is naturally resistant to paclitaxel. While the yew produces a potent cytotoxin, paclitaxel, it is immune to paclitaxel’s cytotoxic action.MethodsTubulin sequence data for plant species were obtained from Alberta 1000 Plants Initiative. Sequences were assembled with Trinity de novo assembly program and tubulin identified. Homology modeling using MODELLER software was done to generate structures for yew tubulin. Molecular dynamics simulations and molecular mechanics Poisson–Boltzmann calculations were performed with the Amber package to determine binding affinity of paclitaxel to yew tubulin. ClustalW2 program and PHYLIP package were used to perform phylogenetic analysis on plant tubulin sequences.ResultsWe specifically analyzed several important regions in tubulin structure: the high-affinity paclitaxel binding site, as well as the intermediate binding site and microtubule nanopores. Our analysis indicates that the high-affinity binding site contains several substitutions compared to human tubulin, all of which reduce the binding energy of paclitaxel.ConclusionsThe yew has achieved a significant reduction of paclitaxel’s affinity for its tubulin by utilizing several specific residue changes in the binding pocket for paclitaxel.

Computational Design and Biological Testing of Highly Cytotoxic Colchicine Ring A Modifications

Microtubules are the primary target for many anti‐cancer drugs, the majority of which bind specifically to β‐tubulin. The existence of several β‐tubulin isotypes, coupled with their varied expression in normal and cancerous cells provides a platform upon which to construct selective chemotherapeutic agents. We have examined five prevalent human β‐tubulin isotypes and identified the colchicine‐binding site as the most promising for drug design based on specificity. Using this binding site as a template, we have designed several colchicine derivatives and computationally probed them for affinity to the β‐tubulin isotypes. These compounds were synthesized and subjected to cytotoxicity assays to determine their effectiveness against several cancerous cell lines. We observed a correlation between computational‐binding predictions and experimentally determined IC50 values, demonstrating the utility of computational screening in the design of more effective colchicine derivatives. The most promising derivative exhibited an IC50 approximately threefold lower than values previously reported for either colchicine or paclitaxel, demonstrating the utility of computational design and assessment of binding to tubulin.

Experimental and computational study of the interaction of novel colchicinoids with a recombinant human αI/βI-tubulin heterodimer.

The binding free energies on human tubulin of selected colchicine and thiocolchicine compounds were determined. Two methods were used for the determination of binding free energies: one is based on theoretical prediction simulating the dissociation of the compound from tubulin using a series of molecular dynamics simulations, and the other method involves a series of experiments that measured the affinity of the compound on a synthetically expressed and purified tubulin protein using a spectrofluorometric technique.

Novel mutations involving βI-, βIIA-, or βIVB-tubulin isotypes with functional resemblance to βIII-tubulin in breast cancer

Tubulin is the target for very widely used anti-tumor drugs, including Vinca alkaloids, taxanes, and epothilones, which are an important component of chemotherapy in breast cancer and other malignancies. Paclitaxel and other tubulin-targeting drugs bind to the β subunit of tubulin, which is a heterodimer of α and β subunits. β-Tubulin exists in the form of multiple isotypes, which are differentially expressed in normal and neoplastic cells and differ in their ability to bind to drugs. Among them, the βIII isotype is overexpressed in many aggressive and metastatic cancers and may serve as a prognostic marker in certain types of cancer. The underpinning mechanisms accounting for the overexpression of this isotype in cancer cells are unclear. To better understand the role of β-tubulin isotypes in cancer, we analyzed over 1000 clones from 90 breast cancer patients, sequencing their β-tubulin isotypes, in search of novel mutations. We have elucidated two putative emerging molecular subgroups of invasive breast cancer, each of which involve mutations in the βI-, βIIA-, or βIVB isotypes of tubulin that increase their structural, and possibly functional, resemblance to the βIII isotype. A unifying feature of the first of the two subgroups is the mutation of the highly reactive C239 residue of βI- or βIVB-tubulin to L239, R239, Y239, or P239, culminating in probable conversion of these isotypes from ROS-sensitive to ROS-resistant species. In the second subgroup, βI, βIIA, and βIVB have up to seven mutations to the corresponding residues in βIII-tubulin. Given that βIII-tubulin has emerged as a pro-survival factor, overexpression of this isotype may confer survival advantages to certain cancer cell types. In this mini-review, we bring attention to a novel mechanism by which cancer cells may undergo adaptive mutational changes involving alternate β-tubulin isotypes to make them acquire some of the pro-survival properties of βIII-tubulin. These “hybrid” tubulins, combining the sequences and/or properties of two wild-type tubulins (βIII and either βI, βIIA, or βIVB), are novel isotypes expressed solely in cancer cells and may contribute to the molecular understanding and stratification of invasive breast cancer and provide novel molecular targets for rational drug development.

Activation of Hydrogen Peroxide to Peroxytetradecanoic Acid Is Responsible for Potent Inhibition of Protein Tyrosine Phosphatase CD45

Hydrogen peroxide induces oxidation and consequently inactivation of many protein tyrosine phosphatases. It was found that hydrogen peroxide, in the presence of carboxylic acids, was efficiently activated to form even more potent oxidant - peroxy acid. We have found that peroxytetradecanoic acid decreases the enzymatic activity of CD45 phosphatase significantly more than hydrogen peroxide. Our molecular docking computational analysis suggests that peroxytetradecanoic acid has a higher binding affinity to the catalytic center of CD45 than hydrogen peroxide.

Novel Colchicine Derivatives and their Anti-Cancer Activity

In this paper we provide an overview of the status of various colchicine derivatives in preclinical development with special focus on their anti-cancer activity. We discuss several groups of compounds that have been designed to differentially bind with specific affinities for tubulin β isotypes, especially in regard to βIII, which is commonly over-expressed in cancer. Computational prediction, protein-based and cell-based assays are summarized as well as some animal tests conducted on these compounds. It is concluded that an untapped potential exists for exploiting the colchicine scaffold as a pharmacophore with the possibility of increasing its affinity for tubulin isotypes overexpressed in cancer and decreasing it for normal cells thereby widening the therapeutic window.

Inactivation of Protein Tyrosine Phosphatases by Peracids Correlates with the Hydrocarbon Chain Length

Background/Aims: Protein tyrosine phosphatases are crucial enzymes controlling numerous physiological and pathophysiological events and can be regulated by oxidation of the catalytic domain cysteine residue. Peracids are highly oxidizing compounds, and thus may induce inactivation of PTPs. The aim of the present study was to evaluate the inhibitory effect of peracids with different length of hydrocarbon chain on the activity of selected PTPs. Methods: The enzymatic activity of human CD45, PTP1B, LAR, bacterial YopH was assayed under the cell-free conditions, and activity of cellular CD45 in human Jurkat cell lysates. The molecular docking and molecular dynamics were performed to evaluate the peracids binding to the CD45 active site. Results: Here we demonstrate that peracids reduce enzymatic activity of recombinant CD45, PTP1B, LAR, YopH and cellular CD45. Our studies indicate that peracids are more potent inhibitors of CD45 than hydrogen peroxide (with an IC50 value equal to 25 nM for peroctanoic acid and 8 µM for hydrogen peroxide). The experimental data show that the inactivation caused by peracids is dependent on hydrocarbon chain length of peracids with maximum inhibitory effect of medium-chain peracids (C8-C12 acyl chain), which correlates with calculated binding affinities to the CD45 active site. Conclusion: Peracids are potent inhibitors of PTPs with the strongest inhibitory effect observed for medium-chain peracids.