De novo determination of peptide structure with solid-state magic-angle spinning NMR spectroscopy

The three-dimensional structure of the chemotactic peptide N-formyl-l-Met-l-Leu-l-Phe-OH was determined by using solid-state NMR (SSNMR). The set of SSNMR data consisted of 16 (13)C-(15)N distances and 18 torsion angle constraints (on 10 angles), recorded from uniformly (13)C,(15)N- and (15)N-labeled samples. The peptide's structure was calculated by means of simulated annealing and a newly developed protocol that ensures that all of conformational space, consistent with the structural constraints, is searched completely. The result is a high-quality structure of a molecule that has thus far not been amenable to single-crystal diffraction studies. The extensions of the SSNMR techniques and computational methods to larger systems appear promising.

Local rule-based theory of virus shell assembly

A local rule-based theory is developed which shows that the self-assembly of icosahedral virus shells may depend on only the lower-level interactions of a protein subunit with its neighbors--i.e., on local rules rather than on larger structural building blocks. The local rule theory provides a framework for understanding the assembly of icosahedral viruses. These include both viruses that fall in the quasiequivalence theory of Caspar and Klug and the polyoma virus structure, which violates quasi-equivalence and has puzzled researchers since it was first observed. Local rules are essentially templates for energetically favorable arrangements. The tolerance margins for these rules are investigated through computer simulations. When these tolerance margins are exceeded in a particular way, the result is a "spiraling" malformation that has been observed in nature.

Membrane tension controls adhesion positioning at the leading edge of cells

Pontes et al. show that plasma membrane mechanics exerts an upstream control during cell motility. Variations in plasma membrane tension orchestrate the behavior of the cell leading edge, with increase–decrease cycles in tension promoting adhesion row positioning.

Cell migration is dependent on adhesion dynamics and actin cytoskeleton remodeling at the leading edge. These events may be physically constrained by the plasma membrane. Here, we show that the mechanical signal produced by an increase in plasma membrane tension triggers the positioning of new rows of adhesions at the leading edge. During protrusion, as membrane tension increases, velocity slows, and the lamellipodium buckles upward in a myosin II–independent manner. The buckling occurs between the front of the lamellipodium, where nascent adhesions are positioned in rows, and the base of the lamellipodium, where a vinculin-dependent clutch couples actin to previously positioned adhesions. As membrane tension decreases, protrusion resumes and buckling disappears, until the next cycle. We propose that the mechanical signal of membrane tension exerts upstream control in mechanotransduction by periodically compressing and relaxing the lamellipodium, leading to the positioning of adhesions at the leading edge of cells.

Engrailed (Gln50-->Lys) homeodomain-DNA complex at 1.9 A resolution: structural basis for enhanced affinity and altered specificity

BACKGROUND

The homeodomain is one of the key DNA-binding motifs used in eukaryotic gene regulation, and homeodomain proteins play critical roles in development. The residue at position 50 of many homeodomains appears to determine the differential DNA-binding specificity, helping to distinguish among binding sites of the form TAATNN. However, the precise role(s) of residue 50 in the differential recognition of alternative sites has not been clear. None of the previously determined structures of homeodomain-DNA complexes has shown evidence for a stable hydrogen bond between residue 50 and a base, and there has been much discussion, based in part on NMR studies, about the potential importance of water-mediated contacts. This study was initiated to help clarify some of these issues.

RESULTS

The crystal structure of a complex containing the engrailed Gln50-->Lys variant (QK50) with its optimal binding site TAATCC (versus TAATTA for the wild-type protein) has been determined at 1.9 A resolution. The overall structure of the QK50 variant is very similar to that of the wild-type complex, but the sidechain of Lys50 projects directly into the major groove and makes several hydrogen bonds to the O6 and N7 atoms of the guanines at base pairs 5 and 6. Lys50 also makes an additional water-mediated contact with the guanine at base pair 5 and has an alternative conformation that allows a hydrogen bond with the O4 of the thymine at base pair 4.

CONCLUSIONS

The structural context provided by the folding and docking of the engrailed homeodomain allows Lys50 to make remarkably favorable contacts with the guanines at base pairs 5 and 6 of the binding site. Although many different residues occur at position 50 in different homeodomains, and although numerous position 50 variants have been constructed, the most striking examples of altered specificity usually involve introducing or removing a lysine sidechain from position 50. This high-resolution structure also confirms the critical role of Asn51 in homeodomain-DNA recognition and further clarifies the roles of water molecules near residues 50 and 51.

A prediction tool for hospital-acquired pressure ulcers among surgical patients: Surgical pressure ulcer risk score

Surgical patients are prone to developing hospital-acquired pressure ulcers (HAPU). Therefore, a better prediction tool is needed to predict risk using preoperative data. This study aimed to determine, from previously published HAPU risk factors, which factors are significant among our surgical population and to develop a prediction tool that identifies pressure ulcer risk before the operation. A literature review was first performed to elicit all the published HAPU risk factors before conducting a retrospective case-control study using medical records. The known HAPU risks were compared between patients with HAPU and without HAPU who underwent operations during the same period (July 2015-December 2016). A total of 80 HAPU cases and 189 controls were analysed. Multivariate logistic regression analyses identified eight significant risk factors: age ≥ 75 years, female gender, American Society of Anaesthesiologists ≥ 3, body mass index < 23, preoperative Braden score ≤ 14, anaemia, respiratory disease, and hypertension. The model had bootstrap-corrected c-statistic 0.78 indicating good discrimination. A cut-off score of ≥6 is strongly predictive, with a positive predictive value of 73.2% (confidence interval [CI]: 59.7%-84.2%) and a negative predictive value of 80.7% (CI: 74.3%-86.1%). SPURS contributes to the preoperative identification of pressure ulcer risk that could help nurses implement preventive measures earlier.

Combination Therapy and the Evolution of Resistance: The Theoretical Merits of Synergism and Antagonism in Cancer

The search for effective combination therapies for cancer has focused heavily on synergistic combinations because they exhibit enhanced therapeutic efficacy at lower doses. Although synergism is intuitively attractive, therapeutic success often depends on whether drug resistance develops. The impact of synergistic combinations (vs. antagonistic or additive combinations) on the process of drug-resistance evolution has not been investigated. In this study, we use a simplified computational model of cancer cell numbers in a population of drug-sensitive, singly-resistant, and fully-resistant cells to simulate the dynamics of resistance evolution in the presence of two-drug combinations. When we compared combination therapies administered at the same combination of effective doses, simulations showed synergistic combinations most effective at delaying onset of resistance. Paradoxically, when the therapies were compared using dose combinations with equal initial efficacy, antagonistic combinations were most successful at suppressing expansion of resistant subclones. These findings suggest that, although synergistic combinations could suppress resistance through early decimation of cell numbers (making them "proefficacy" strategies), they are inherently fragile toward the development of single resistance. In contrast, antagonistic combinations suppressed the clonal expansion of singly-resistant cells, making them "antiresistance" strategies. The distinction between synergism and antagonism was intrinsically connected to the distinction between offensive and defensive strategies, where offensive strategies inflicted early casualties and defensive strategies established protection against anticipated future threats. Our findings question the exclusive focus on synergistic combinations and motivate further consideration of nonsynergistic combinations for cancer therapy.Significance: Computational simulations show that if different combination therapies have similar initial efficacy in cancers, then nonsynergistic drug combinations are more likely than synergistic drug combinations to provide a long-term defense against the evolution of therapeutic resistance. Cancer Res; 78(9); 2419-31. ©2018 AACR.

Wnt proteins synergize to activate β-catenin signaling

Wnt ligands are involved in diverse signaling pathways that are active during development, maintenance of tissue homeostasis and in various disease states. While signaling regulated by individual Wnts has been extensively studied, Wnts are rarely expressed alone, and the consequences of Wnt gene co-expression are not well understood. Here, we studied the effect of co-expression of Wnts on the β-catenin signaling pathway. While some Wnts are deemed 'non-canonical' due to their limited ability to activate β-catenin when expressed alone, unexpectedly, we find that multiple Wnt combinations can synergistically activate β-catenin signaling in multiple cell types. WNT1- and WNT7B-mediated synergistic Wnt signaling requires FZD5, FZD8 and LRP6, as well as the WNT7B co-receptors GPR124 (also known as ADGRA2) and RECK. Unexpectedly, this synergistic signaling occurs downstream of β-catenin stabilization, and is correlated with increased lysine acetylation of β-catenin. Wnt synergy provides a general mechanism to confer increased combinatorial control over this important regulatory pathway.

Cell-delivery therapeutics for liver regeneration

For acute, chronic, or hereditary diseases of the liver, cell transplantation therapies can stimulate liver regeneration or serve as a bridge until liver transplantation can be performed. Recently, fetal hepatocytes, stem cells, liver progenitor cells, or other primitive and proliferative cell types have been employed for cell transplantation therapies, in an effort to improve the survival, proliferation, and engraftment of the transplanted cells. Reviewing earlier studies, which achieved success by transplanting mature hepatocytes, we propose that there is a switch-like regulation of liver regeneration that changes state according to a stimulus threshold of extracellular influences such as cytokines, matrices and neighboring cells. Important determinants of a successful clinical outcome include sufficient quantities and functional levels of the transplanted cells (even for short periods to alter the environment), rather than just engraftment levels or survival durations of the exogenously transplanted cells. The relative importance of these determining factors will impact future choices of cell sources, delivery vehicles, and sites of cell transplantation to stimulate liver regeneration for patients with severe liver diseases.

The CD47-binding peptide of thrombospondin-1 induces defenestration of liver sinusoidal endothelial cells

BACKGROUND & AIMS

A fenestrated phenotype is characteristic of liver sinusoidal endothelial cells (LSECs), but liver sinusoids become defenestrated during fibrosis and other liver diseases. Thrombospondin-1 (TSP1) is a matrix glycoprotein with pro-fibrotic effects, and the CD47-binding fragment of TSP1 also has anti-angiogenic effects in endothelial cells. We hypothesized that the CD47-binding fragment of TSP1 could induce defenestration in LSECs through the Rho-Rho kinase (ROCK)-myosin pathway.

METHODS

Freshly isolated rat LSECs were treated with TSP1 or CD47-binding peptides of TSP1. LSEC fenestration was assessed with scanning electron microscopy, and myosin phosphorylation was assessed with immuno-fluorescence.

RESULTS

Treating LSECs with TSP1 caused a dose-dependent loss of fenestrae, and this effect could not be blocked by SB-431542, the TGF-β1 receptor inhibitor. A CD47-binding fragment of TSP1, p4N1, was able to induce defenestration, and a CD47-blocking antibody, B6H12, was able to suppress p4N1-induced defenestration. The p4N1 fragment also caused contraction of fenestra size, correlated with an increase in myosin activation. Pretreatment with Y-237642 (a ROCK inhibitor) prevented p4N1-induced myosin activation and fenestrae decrease. Simvastatin has also been shown to antagonize Rho-ROCK signalling, and we found that simvastatin pretreatment protected LSECs from p4N1-induced myosin activation and defenestration.

CONCLUSIONS

We conclude that CD47 signals through the Rho-ROCK-myosin pathway to induce defenestration in LSECs. In addition, our results show that simvastatin and Y-237642 have a beneficial impact on fenestration in vitro, providing an additional explanation for the efficacy of these compounds for regression of liver fibrosis.

Delineation of lipopolysaccharide (LPS)-binding sites on hemoglobin: from in silico predictions to biophysical characterization

Hemoglobin (Hb) functions as a frontline defense molecule during infection by hemolytic microbes. Binding to LPS induces structural changes in cell-free Hb, which activates the redox activity of the protein for the generation of microbicidal free radicals. Although the interaction between Hb and LPS has implications for innate immune defense, the precise LPS-interaction sites on Hb remain unknown. Using surface plasmon resonance, we found that both the Hb α and β subunits possess high affinity LPS-binding sites, with K(D) in the nanomolar range. In silico analysis of Hb including phospho-group binding site prediction, structure-based sequence comparison, and docking to model the protein-ligand interactions showed that Hb possesses evolutionarily conserved surface cationic patches that could function as potential LPS-binding sites. Synthetic Hb peptides harboring predicted LPS-binding sites served to validate the computational predictions. Surface plasmon resonance analysis differentiated LPS-binding peptides from non-binders. Binding of the peptides to lipid A was further substantiated by a fluorescent probe displacement assay. The LPS-binding peptides effectively neutralized the endotoxicity of LPS in vitro. Additionally, peptide B59 spanning residues 59-95 of Hbβ attached to the surface of Gram-negative bacteria as shown by flow cytometry and visualized by immunogold-labeled scanning electron microscopy. Site-directed mutagenesis of the Hb subunits further confirmed the function of the predicted residues in binding to LPS. In summary, the integration of computational predictions and biophysical characterization has enabled delineation of multiple LPS-binding hot spots on the Hb molecule.

Exome Sequencing Reveals Germline SMAD9 Mutation That Reduces Phosphatase and Tensin Homolog Expression and Is Associated With Hamartomatous Polyposis and Gastrointestinal Ganglioneuromas

Hamartomatous polyposis syndromes (HPS) account for a small but appreciable proportion of inherited gastrointestinal cancer predisposition syndromes; patients with HPS have an increased risk for colon and extracolonic malignancies. We present a unique case of familial juvenile polyposis syndrome associated with gastrointestinal ganglioneuromas of unknown etiology. The patient was tested for HPS-associated genes, but no mutation was detected. Exome sequencing identified a germline heterozygous mutation in SMAD9 (SMAD9(V90M)). This mutation was predicted to be an activating mutation. HEK cells transfected to express SMAD9(V90M) had reduced expression of phosphatase and tensin homolog; this reduction was also observed in a polyp from the patient. We have therefore identified a new susceptibility locus for HPS. Patients with hamartomatous polyposis in the colon associated with ganglioneuromatosis should be referred for genetic assessments.

Hepatic stellate cell-targeted delivery of hepatocyte growth factor transgene via bile duct infusion enhances its expression at fibrotic foci to regress dimethylnitrosamine-induced liver fibrosis

Liver fibrosis generates fibrotic foci with abundant activated hepatic stellate cells and excessive collagen deposition juxtaposed with healthy regions. Targeted delivery of antifibrotic therapeutics to hepatic stellate cells (HSCs) might improve treatment outcomes and reduce adverse effects on healthy tissue. We delivered the hepatocyte growth factor (HGF) gene specifically to activated hepatic stellate cells in fibrotic liver using vitamin A-coupled liposomes by retrograde intrabiliary infusion to bypass capillarized hepatic sinusoids. The antifibrotic effects of DsRed2-HGF vector encapsulated within vitamin A-coupled liposomes were validated by decreases in fibrotic markers in vitro. Fibrotic cultures transfected with the targeted transgene showed a significant decrease in fibrotic markers such as transforming growth factor-β1. In rats, dimethylnitrosamine-induced liver fibrosis is manifested by an increase in collagen deposition and severe defenestration of sinusoidal endothelial cells. The HSC-targeted transgene, administered via retrograde intrabiliary infusion in fibrotic rats, successfully reduced liver fibrosis markers alpha-smooth muscle actin and collagen, accompanied by an increase in the expression of DsRed2-HGF near the fibrotic foci. Thus, targeted delivery of HGF gene to hepatic stellate cells increased the transgene expression at the fibrotic foci and strongly enhanced its antifibrotic effects.

An expectation-maximization algorithm enables accurate ecological modeling using longitudinal microbiome sequencing data

BACKGROUND

The dynamics of microbial communities is driven by a range of interactions from symbiosis to predator-prey relationships, the majority of which are poorly understood. With the increasing availability of high-throughput microbiome taxonomic profiling data, it is now conceivable to directly learn the ecological models that explicitly define microbial interactions and explain community dynamics. The applicability of these approaches is severely limited by the lack of accurate absolute cell density measurements (biomass).

METHODS

We present a new computational approach that resolves this key limitation in the inference of generalized Lotka-Volterra models (gLVMs) by coupling biomass estimation and model inference with an expectation-maximization algorithm (BEEM).

RESULTS

BEEM outperforms the state-of-the-art methods for inferring gLVMs, while simultaneously eliminating the need for additional experimental biomass data as input. BEEM's application to previously inaccessible public datasets (due to the lack of biomass data) allowed us to construct ecological models of microbial communities in the human gut on a per-individual basis, revealing personalized dynamics and keystone species.

CONCLUSIONS

BEEM addresses a key bottleneck in "systems analysis" of microbiomes by enabling accurate inference of ecological models from high throughput sequencing data without the need for experimental biomass measurements.

Plasmin triggers a switch-like decrease in thrombospondin-dependent activation of TGF-β1

Transforming growth factor-β1 (TGF-β1) is a potent regulator of extracellular matrix production, wound healing, differentiation, and immune response, and is implicated in the progression of fibrotic diseases and cancer. Extracellular activation of TGF-β1 from its latent form provides spatiotemporal control over TGF-β1 signaling, but the current understanding of TGF-β1 activation does not emphasize cross talk between activators. Plasmin (PLS) and thrombospondin-1 (TSP1) have been studied individually as activators of TGF-β1, and in this work we used a systems-level approach with mathematical modeling and in vitro experiments to study the interplay between PLS and TSP1 in TGF-β1 activation. Simulations and steady-state analysis predicted a switch-like bistable transition between two levels of active TGF-β1, with an inverse correlation between PLS and TSP1. In particular, the model predicted that increasing PLS breaks a TSP1-TGF-β1 positive feedback loop and causes an unexpected net decrease in TGF-β1 activation. To test these predictions in vitro, we treated rat hepatocytes and hepatic stellate cells with PLS, which caused proteolytic cleavage of TSP1 and decreased activation of TGF-β1. The TGF-β1 activation levels showed a cooperative dose response, and a test of hysteresis in the cocultured cells validated that TGF-β1 activation is bistable. We conclude that switch-like behavior arises from natural competition between two distinct modes of TGF-β1 activation: a TSP1-mediated mode of high activation and a PLS-mediated mode of low activation. This switch suggests an explanation for the unexpected effects of the plasminogen activation system on TGF-β1 in fibrotic diseases in vivo, as well as novel prognostic and therapeutic approaches for diseases with TGF-β dysregulation.

Superpixel-based segmentation of muscle fibers in multi-channel microscopy

BACKGROUND

Confetti fluorescence and other multi-color genetic labelling strategies are useful for observing stem cell regeneration and for other problems of cell lineage tracing. One difficulty of such strategies is segmenting the cell boundaries, which is a very different problem from segmenting color images from the real world. This paper addresses the difficulties and presents a superpixel-based framework for segmentation of regenerated muscle fibers in mice.

RESULTS

We propose to integrate an edge detector into a superpixel algorithm and customize the method for multi-channel images. The enhanced superpixel method outperforms the original and another advanced superpixel algorithm in terms of both boundary recall and under-segmentation error. Our framework was applied to cross-section and lateral section images of regenerated muscle fibers from confetti-fluorescent mice. Compared with "ground-truth" segmentations, our framework yielded median Dice similarity coefficients of 0.92 and higher.

CONCLUSION

Our segmentation framework is flexible and provides very good segmentations of multi-color muscle fibers. We anticipate our methods will be useful for segmenting a variety of tissues in confetti fluorecent mice and in mice with similar multi-color labels.

Extracellular haemoglobin upregulates and binds to tissue factor on macrophages: Implications for coagulation and oxidative stress

The mechanisms of crosstalk between haemolysis, coagulation and innate immunity are evolutionarily conserved from the invertebrate haemocyanin to the vertebrate haemoglobin (Hb). In vertebrates, extracellular Hb resulting from haemolytic infections binds bacterial lipopolysaccharide (LPS) to unleash the antimicrobial redox activity of Hb. Because bacterial invasion also upregulates tissue factor (TF), the vertebrate coagulation initiator, we asked whether there may be functional interplay between the redox activity of Hb and the procoagulant activity of TF. Using real-time PCR, TF-specific ELISA, flow cytometry and TF activity assay, we found that Hb upregulated the expression of functional TF in macrophages. ELISA, flow cytometry and immunofluorescence microscopy showed binding between Hb and TF, in isolation and in situ. Bioinformatic analysis of Hb and TF protein sequences showed co-evolution across species, suggesting that Hbβ binds TF. Empirically, TF suppressed the LPS-induced activation of Hb redox activity. Furthermore, Hb desensitised TF to the effects of antioxidants like glutathione or serum. This bi-directional regulation between Hb and TF constitutes a novel link between coagulation and innate immunity. In addition, induction of TF by Hb is a potentially central mechanism for haemolysis to trigger coagulation.

FLIP: a flop for execution signals

Resistance to apoptosis is one of the established hallmarks of cancer cells. This is a function of an imbalance between the proteins that facilitate death execution and those that inhibit apoptosis or promote cell proliferation. The anti-apoptotic protein, FLICE inhibitory protein (FLIP), first identified as a viral protein, is over-expressed in a variety of human pathologies. Initial observations linked FLIP expression to inhibition of death receptor induced apoptosis, due to its structural homology to the cysteine protease, caspase-8. FLIP impedes full processing of pro-caspase-8 to its active form and its release to the cytosol, and by doing so blocks apoptotic signaling downstream of the membrane death initiating signaling complex (DISC). Recent observations have highlighted the complex regulation of this protein and its cross talk with diverse signaling networks and metabolic processes. As FLIP expression is directly associated with chemotherapy resistance, a better understanding of its genomic organization, gene transcription, as well as post-transcriptional regulation could yield novel targets with potential therapeutic implications against drug refractory cancers. In this short review, we provide a brief overview of the structural and functional biology of this somewhat complex protein with direct relevance to carcinogenesis.

Immune cell phenotypes associated with disease severity and long-term neutralizing antibody titers after natural dengue virus infection

Prior immunological exposure to dengue virus can be both protective and disease-enhancing during subsequent infections with different dengue virus serotypes. We provide here a systematic, longitudinal analysis of B cell, T cell, and antibody responses in the same patients. Antibody responses as well as T and B cell activation differentiate primary from secondary responses. Hospitalization is associated with lower frequencies of activated, terminally differentiated T cells and higher percentages of effector memory CD4 T cells. Patients with more severe disease tend to have higher percentages of plasmablasts. This does not translate into long-term antibody titers, since neutralizing titers after 6 months correlate with percentages of specific memory B cells, but not with acute plasmablast activation. Overall, our unbiased analysis reveals associations between cellular profiles and disease severity, opening opportunities to study immunopathology in dengue disease and the potential predictive value of these parameters.

Development of a classification model for non-alcoholic steatohepatitis (NASH) using confocal Raman micro-spectroscopy

Non-alcoholic fatty liver disease (NAFLD) is the most common liver disorder in developed countries [1]. A subset of individuals with NAFLD progress to non-alcoholic steatohepatitis (NASH), an advanced form of NAFLD which predisposes individuals to cirrhosis, liver failure and hepatocellular carcinoma. The current gold standard for NASH diagnosis and staging is based on histological evaluation, which is largely semi-quantitative and subjective. To address the need for an automated and objective approach to NASH detection, we combined Raman micro-spectroscopy and machine learning techniques to develop a classification model based on a well-established NASH mouse model, using spectrum pre-processing, biochemical component analysis (BCA) and logistic regression. By employing a selected pool of biochemical components, we identified biochemical changes specific to NASH and show that the classification model is capable of accurately detecting NASH (AUC=0.85-0.87) in mice. The unique biochemical fingerprint generated in this study may serve as a useful criterion to be leveraged for further validation in clinical samples.

Functional coding haplotypes and machine-learning feature elimination identifies predictors of Methotrexate Response in Rheumatoid Arthritis patients

BACKGROUND

Major challenges in large scale genetic association studies include not only the identification of causative single nucleotide polymorphisms (SNPs), but also accounting for SNP-SNP interactions. This study thus proposes a novel feature engineering approach integrating potentially functional coding haplotypes (pfcHap) with machine-learning (ML) feature selection to identify biologically meaningful, possibly causative genetic factors, that take into consideration potential SNP-SNP interactions within the pfcHap, to best predict for methotrexate (MTX) response in rheumatoid arthritis (RA) patients.

METHODS

Exome sequencing from 349 RA patients were analysed, of which they were split into training and unseen test set. Inferred pfcHaps were combined with 30 non-genetic features to undergo ML recursive feature elimination with cross-validation using the training set. Predictive capacity and robustness of the selected features were assessed using six popular machine learning models through a train set cross-validation and evaluated in an unseen test set.

FINDINGS

Significantly, 100 features (95 pfcHaps, 5 non-genetic factors) were identified to have good predictive performance (AUC: 0.776-0.828; Sensitivity: 0.656-0.813; Specificity: 0.684-0.868) across all six ML models in an unseen test dataset for the prediction of MTX response in RA patients.

INTERPRETATION

Majority of the predictive pfcHap SNPs were predicted to be potentially functional and some of the genes in which the pfcHap resides in were identified to be associated with previously reported MTX/RA pathways.

FUNDING

Singapore Ministry of Health's National Medical Research Council (NMRC) [NMRC/CBRG/0095/2015; CG12Aug17; CGAug16M012; NMRC/CG/017/2013]; National Cancer Center Research Fund and block funding Duke-NUS Medical School.; Singapore Ministry of Education Academic Research Fund Tier 2 grant MOE2019-T2-1-138.

Systematic parameter estimation in data-rich environments for cell signalling dynamics

Motivation: Computational models of biological signalling networks, based on ordinary differential equations (ODEs), have generated many insights into cellular dynamics, but the model-building process typically requires estimating rate parameters based on experimentally observed concentrations. New proteomic methods can measure concentrations for all molecular species in a pathway; this creates a new opportunity to decompose the optimization of rate parameters.

Results: In contrast with conventional parameter estimation methods that minimize the disagreement between simulated and observed concentrations, the SPEDRE method fits spline curves through observed concentration points, estimates derivatives and then matches the derivatives to the production and consumption of each species. This reformulation of the problem permits an extreme decomposition of the high-dimensional optimization into a product of low-dimensional factors, each factor enforcing the equality of one ODE at one time slice. Coarsely discretized solutions to the factors can be computed systematically. Then the discrete solutions are combined using loopy belief propagation, and refined using local optimization. SPEDRE has unique asymptotic behaviour with runtime polynomial in the number of molecules and timepoints, but exponential in the degree of the biochemical network. SPEDRE performance is comparatively evaluated on a novel model of Akt activation dynamics including redox-mediated inactivation of PTEN (phosphatase and tensin homologue).

Availability and implementation: Web service, software and supplementary information are available at www.LtkLab.org/SPEDRE

Supplementary information:Supplementary data are available at Bioinformatics online.

Contact: LisaTK@nus.edu.sg

Superoxide induced inhibition of death receptor signaling is mediated via induced expression of apoptosis inhibitory protein cFLIP

The death inhibitory proteins, cFLIP and Bcl-2, canonically act at different steps to regulate receptor-mediated apoptosis in cancer cells. Here we report that pharmacological or genetic means to effect an increase in intracellular superoxide result in cFLIP upregulation. Interestingly, Bcl-2 overexpression is associated with a concomitant increase in cFLIP, and reducing superoxide sensitizes Bcl-2 overexpressing cancer cells to receptor-mediated apoptosis via downregulation of cFLIP. Moreover, inhibiting glycolytic flux overcomes apoptosis resistance by superoxide-dependent downregulation of cFLIP. Superoxide-induced upregulation of cFLIP is a function of enhanced transcription, as evidenced by increases in cFLIP promoter activity and mRNA abundance. The positive effect of superoxide on cFLIP is mediated through its reaction with nitric oxide to generate peroxynitrite. Corroborating these findings in cell lines, subjecting primary cells derived from lymphoma patients to glucose deprivation ex vivo, as a means to decrease superoxide, not only reduced cFLIP expression but also significantly enhanced death receptor sensitivity. Based on this novel mechanistic insight into the redox regulation of cancer cell fate, modulation of intracellular superoxide could have potential therapeutic implications in cancers in which these two death inhibitory proteins present a therapeutic challenge.

Steady states and dynamics of urokinase-mediated plasmin activation in silico and in vitro

Plasmin (PLS) and urokinase-type plasminogen activator (UPA) are ubiquitous proteases that regulate the extracellular environment. Although they are secreted in inactive forms, they can activate each other through proteolytic cleavage. This mutual interplay creates the potential for complex dynamics, which we investigated using mathematical modeling and in vitro experiments. We constructed ordinary differential equations to model the conversion of precursor plasminogen into active PLS, and precursor urokinase (scUPA) into active urokinase (tcUPA). Although neither PLS nor UPA exhibits allosteric cooperativity, modeling showed that cooperativity occurred at the system level because of substrate competition. Computational simulations and bifurcation analysis predicted that the system would be bistable over a range of parameters for cooperativity and positive feedback. Cell-free experiments with recombinant proteins tested key predictions of the model. PLS activation in response to scUPA stimulus was found to be cooperative in vitro. Finally, bistability was demonstrated in vitro by the presence of two significantly different steady-state levels of PLS activation for the same levels of stimulus. We conclude that ultrasensitive, bistable activation of UPA-PLS is possible in the presence of substrate competition. An ultrasensitive threshold for activation of PLS and UPA would have ramifications for normal and disease processes, including angiogenesis, metastasis, wound healing, and fibrosis.

Non-canonical Activation of Akt in Serum-Stimulated Fibroblasts, Revealed by Comparative Modeling of Pathway Dynamics

The dynamic behaviors of signaling pathways can provide clues to pathway mechanisms. In cancer cells, excessive phosphorylation and activation of the Akt pathway is responsible for cell survival advantages. In normal cells, serum stimulation causes brief peaks of extremely high Akt phosphorylation before reaching a moderate steady-state. Previous modeling assumed this peak and decline behavior (i.e., “overshoot”) was due to receptor internalization. In this work, we modeled the dynamics of the overshoot as a tool for gaining insight into Akt pathway function. We built an ordinary differential equation (ODE) model describing pathway activation immediately upstream of Akt phosphorylation at Thr³⁰⁸ (Aktp³⁰⁸). The model was fit to experimental measurements of Aktp³⁰⁸, total Akt, and phosphatidylinositol (3,4,5)-trisphosphate (PIP3), from mouse embryonic fibroblasts with serum stimulation. The canonical Akt activation model (the null hypothesis) was unable to recapitulate the observed delay between the peak of PIP3 (at 2 minutes), and the peak of Aktp³⁰⁸ (at 30–60 minutes). From this we conclude that the peak and decline behavior of Aktp³⁰⁸ is not caused by PIP3 dynamics. Models for alternative hypotheses were constructed by allowing an arbitrary dynamic curve to perturb each of 5 steps of the pathway. All 5 of the alternative models could reproduce the observed delay. To distinguish among the alternatives, simulations suggested which species and timepoints would show strong differences. Time-series experiments with membrane fractionation and PI3K inhibition were performed, and incompatible hypotheses were excluded. We conclude that the peak and decline behavior of Aktp³⁰⁸ is caused by a non-canonical effect that retains Akt at the membrane, and not by receptor internalization. Furthermore, we provide a novel spline-based method for simulating the network implications of an unknown effect, and we demonstrate a process of hypothesis management for guiding efficient experiments.

Influential pathways of cell signalling (such as PI3K/Akt) are routinely communicated using simple textbook-like diagrams that show only the most widely-accepted steps of the pathway. At the same time, there are countless other molecular influences relevant to each pathway, documented in the published literature, and more are being published every week. It should perhaps come as little surprise that during a routine observation of the Akt activation pathway, a simulation of the canonical model was mathematically incompatible with our observed dynamics. To progress beyond the standard, simplified model without testing an unreasonable number of molecular candidates individually, we employed computational modeling to analyze the dynamics of pathway activation. We asked when and where a non-canonical deviation could occur, relative to the canonical pathway. We used the timing of downstream activation to solve for the possible times of upstream initiation. By categorizing unknown effects by their dynamics, we were able to prune away implausible hypotheses using an efficient number of in vitro experiments. At the end we had a single plausible explanation for non-canonical Akt activation in our cells, and we confirmed experimentally that Akt is retained at the membrane after PIP3 is no longer present.