Pseudo-nullclines enable the analysis and prediction of signaling model dynamics

A powerful method to qualitatively analyze a 2D system is the use of nullclines, curves which separate regions of the plane where the sign of the time derivatives is constant, with their intersections corresponding to steady states. As a quick way to sketch the phase portrait of the system, they can be sufficient to understand the qualitative dynamics at play without integrating the differential equations. While it cannot be extended straightforwardly for dimensions higher than 2, sometimes the phase portrait can still be projected onto a 2-dimensional subspace, with some curves becoming pseudo-nullclines. In this work, we study cell signaling models of dimension higher than 2 with behaviors such as oscillations and bistability. Pseudo-nullclines are defined and used to qualitatively analyze the dynamics involved. Our method applies when a system can be decomposed into 2 modules, mutually coupled through 2 scalar variables. At the same time, it helps track bifurcations in a quick and efficient manner, key for understanding the different behaviors. Our results are both consistent with the expected dynamics, and also lead to new responses like excitability. Further work could test the method for other regions of parameter space and determine how to extend it to three-module systems.

Frequency-preference response in covalent modification cycles under substrate sequestration conditions

Covalent modification cycles (CMCs) are basic units of signaling systems and their properties are well understood. However, their behavior has been mostly characterized in situations where the substrate is in excess over the modifying enzymes. Experimental data on protein abundance suggest that the enzymes and their target proteins are present in comparable concentrations, leading to substrate sequestration by the enzymes. In this enzyme-in-excess regime, CMCs have been shown to exhibit signal termination, the ability of the product to return to a stationary value lower than its peak in response to constant stimulation, while this stimulation is still active, with possible implications for the ability of systems to adapt to environmental inputs. We characterize the conditions leading to signal termination in CMCs in the enzyme-in-excess regime. We also demonstrate that this behavior leads to a preferred frequency response (band-pass filters) when the cycle is subjected to periodic stimulation, whereas the literature reports that CMCs investigated so far behave as low-pass filters. We characterize the relationship between signal termination and the preferred frequency response to periodic inputs and we explore the dynamic mechanism underlying these phenomena. Finally, we describe how the behavior of CMCs is reflected in similar types of responses in the cascades of which they are part. Evidence of protein abundance in vivo shows that enzymes and substrates are present in comparable concentrations, thus suggesting that signal termination and frequency-preference response to periodic inputs are also important dynamic features of cell signaling systems, which have been overlooked.

Discriminating between negative cooperativity and ligand binding to independent sites using pre-equilibrium properties of binding curves

Negative cooperativity is a phenomenon in which the binding of a first ligand or substrate molecule decreases the rate of subsequent binding. This definition is not exclusive to ligand-receptor binding, it holds whenever two or more molecules undergo two successive binding events. Negative cooperativity turns the binding curve more graded and cannot be distinguished from two independent and different binding events based on equilibrium measurements only. The need of kinetic data for this purpose was already reported. Here, we study the binding response as a function of the amount of ligand, at different times, from very early times since ligand is added and until equilibrium is reached. Over those binding curves measured at different times, we compute the dynamic range: the fold change required in input to elicit a change from 10 to 90% of maximum output, finding that it evolves in time differently and controlled by different parameters in the two situations that are identical in equilibrium. Deciphering which is the microscopic model that leads to a given binding curve adds understanding on the molecular mechanisms at play, and thus, is a valuable tool. The methods developed in this article were tested both with simulated and experimental data, showing to be robust to noise and experimental constraints.

When two successive events occur, it may make sense to know if they affect somehow each other, particularly if the properties of the second event are modified by the occurrence of the first one. Two scenarios lead to the same overall outcome: first, the two events are identical but they interfere with each other, and second, the two events are independent but non identical. The interference caused in the first scenario produces the same result as having a second event with different properties. Now, let’s name these events as bindings, the interference as negative cooperativity, and the non-identical events as independent binding. In this work we focus on the dynamic process by which the two scenarios produce the same result. We selected a relevant but not characterized before property of the binding process, called its dynamic range, and found it behaves differently in these two scenarios and controlled by different parameters of interest. Based on this feature, we developed and algorithm to distinguish between negative cooperativity and independent binding based on the time evolution of the dynamic range. This tool allows to discover the microscopic model behind the data and may be useful in other similar problems in cell signaling.

Robustness in spatially driven bistability in signaling systems

Biological systems are spatially organized. This microscopic heterogeneity has been shown to produce emergent complex behaviors such as bistability. Even though the connection between spatiality and dynamic response is essential to understand biological output, its robustness and extent has not been sufficiently explored. This work focuses on a previously described system which is composed of two monostable modules acting on different cellular compartments and sharing species through linear shuttling reactions. One of the two main purposes of this paper is to quantify the frequency of occurrence of bistability throughout the parameter space and to identify which parameters and in which value ranges control the emergence and the properties of bistability. We found that a very small fraction of the sampled parameter space produced a bistable response. Most importantly, shuttling parameters were among the most influential ones to control this property. The other goal of this paper is to simplify the same system as much as possible without losing compartment-induced bistability. This procedure provided a simplified model that still connects two monostable systems by a reduced set of linear shuttling reactions that circulates all the species around the two compartments. Bistable systems are one of the main building blocks of more complex behaviors such as oscillations, memory, and digitalization. Therefore, we expect that the proposed minimal system provides insight into how these behaviors can arise from compartmentalization.

Cell cycle dynamics of mouse embryonic stem cells in the ground state and during transition to formative pluripotency

Mouse embryonic stem cells (mESCs) can be maintained as homogeneous populations in the ground state of pluripotency. Release from this state in minimal conditions allows to obtain cells that resemble those of the early post-implantation epiblast, providing an important developmental model to study cell identity transitions. However, the cell cycle dynamics of mESCs in the ground state and during its dissolution have not been extensively studied. By performing live imaging experiments of mESCs bearing cell cycle reporters, we show here that cells in the pluripotent ground state display a cell cycle structure comparable to the reported for mESCs in serum-based media. Upon release from self-renewal, the cell cycle is rapidly accelerated by a reduction in the length of the G1 phase and of the S/G2/M phases, causing an increased proliferation rate. Analysis of cell lineages indicates that cell cycle variables of sister cells are highly correlated, suggesting the existence of inherited cell cycle regulators from the parental cell. Together with a major morphological reconfiguration upon differentiation, our findings support a correlation between this in vitro model and early embryonic events.

Ultrasensitivity in signaling cascades revisited: Linking local and global ultrasensitivity estimations

Ultrasensitive response motifs, capable of converting graded stimuli into binary responses, are well-conserved in signal transduction networks. Although it has been shown that a cascade arrangement of multiple ultrasensitive modules can enhance the system's ultrasensitivity, how a given combination of layers affects a cascade's ultrasensitivity remains an open question for the general case. Here, we introduce a methodology that allows us to determine the presence of sequestration effects and to quantify the relative contribution of each module to the overall cascade's ultrasensitivity. The proposed analysis framework provides a natural link between global and local ultrasensitivity descriptors and it is particularly well-suited to characterize and understand mathematical models used to study real biological systems. As a case study, we have considered three mathematical models introduced by O'Shaughnessy et al. to study a tunable synthetic MAPK cascade, and we show how our methodology can help modelers better understand alternative models.

Signaling cascades transmit information downstream and upstream but unlikely simultaneously

Background

Signal transduction is the process through which cells communicate with the external environment, interpret stimuli and respond to them. This mechanism is controlled by signaling cascades, which play the role of intracellular transmitter, being able to transmit biochemical information between cell membrane and nucleus. In theory as well as in practice, it has been shown that a perturbation can propagate upstream (and not only downstream) a cascade, by a mechanism known as retroactivity. This study aims to compare the conditions on biochemical parameters which favor one or the other direction of signaling in such a cascade.

Results

From a mathematical point of view, we show that the steady states of a cascade of arbitrary length n are described by an iterative map of second order, meaning that the cascade tiers are actually coupled three-by-three. We study the influence of the biochemical parameters in the control of the direction of transmission – upstream and/or downstream – along a signaling cascade. A numerical and statistical approach, based on the random scan of parameters describing a 3-tier signaling cascade, provides complementary findings to the analytical study. In particular, computing the likelihood of parameters with respect to various signaling regimes, we identify conditions on biochemical parameters which enhance a specific direction of propagation corresponding to forward or retro-signaling regimes. A compact graphical representation is designed to relay the gist of these conditions.

Conclusions

The values of biochemical parameters such as kinetic rates, Michaelis-Menten constants, total concentrations of kinases and of phosphatases, determine the propensity of a cascade to favor or impede downstream or upstream signal transmission. We found that generally there is an opposition between parameter sets favoring forward and retro-signaling regimes. Therefore, on one hand our study supports the idea that in most cases, retroactive effects can be neglected when a cascade which is efficient in forward signaling, is perturbed by an external ligand inhibiting the activation at some tier of the cascade. This result is relevant for therapeutic methodologies based on kinase inhibition. On the other hand, our study highlights a less-known part of the parameter space where, although the forward signaling is inefficient, the cascade can interestingly act as a retro-signaling device.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-016-0303-2) contains supplementary material, which is available to authorized users.

Utilization of extracellular information before ligand-receptor binding reaches equilibrium expands and shifts the input dynamic range

Cell signaling systems sense and respond to ligands that bind cell surface receptors. These systems often respond to changes in the concentration of extracellular ligand more rapidly than the ligand equilibrates with its receptor. We demonstrate, by modeling and experiment, a general "systems level" mechanism cells use to take advantage of the information present in the early signal, before receptor binding reaches a new steady state. This mechanism, pre-equilibrium sensing and signaling (PRESS), operates in signaling systems in which the kinetics of ligand-receptor binding are slower than the downstream signaling steps, and it typically involves transient activation of a downstream step. In the systems where it operates, PRESS expands and shifts the input dynamic range, allowing cells to make different responses to ligand concentrations so high as to be otherwise indistinguishable. Specifically, we show that PRESS applies to the yeast directional polarization in response to pheromone gradients. Consideration of preexisting kinetic data for ligand-receptor interactions suggests that PRESS operates in many cell signaling systems throughout biology. The same mechanism may also operate at other levels in signaling systems in which a slow activation step couples to a faster downstream step.

Characterization of the reconstituted UTase/UR-PII-NRII-NRI bicyclic signal transduction system that controls the transcription of nitrogen-regulated (Ntr) genes in Escherichia coli

A reconstituted UTase/UR-PII-NRII-NRI bicyclic cascade regulated PII uridylylation and NRI phosphorylation in response to glutamine. We examined the sensitivity and robustness of the responses of the individual cycles and of the bicyclic system. The sensitivity of the glutamine response of the upstream UTase/UR-PII monocycle depended upon the PII concentration, and we show that PII exerted substrate inhibition of the UTase activity of UTase/UR, potentially contributing to this dependence of sensitivity on PII. In the downstream NRII-NRI monocycle, PII controlled NRI phosphorylation state, and the response to PII was hyperbolic at both saturating and unsaturating NRI concentration. As expected from theory, the level of NRI∼P produced by the NRII-NRI monocycle was robust to changes in the NRII or NRI concentrations when NRI was in excess over NRII, as long as the NRII concentration was above a threshold value, an example of absolute concentration robustness (ACR). Because of the parameters of the system, at physiological protein levels and ratios of NRI to NRII, the level of NRI∼P depended upon both protein concentrations. In bicyclic UTase/UR-PII-NRII-NRI systems, the NRI phosphorylation state response to glutamine was always hyperbolic, regardless of the PII concentration or sensitivity of the upstream UTase/UR-PII cycle. In these bicyclic systems, NRI phosphorylation state was only robust to variation in the PII/NRII ratio within a narrow range; when PII was in excess NRI∼P was low, and when NRII was in excess NRI phosphorylation was elevated, throughout the physiological range of glutamine concentrations. Our results show that the bicyclic system produced a graded response of NRI phosphorylation to glutamine under a range of conditions, and that under most conditions the response of NRI phosphorylation state to glutamine levels depended on the concentrations of NRI, NRII, and PII.

Retroactive signaling in short signaling pathways

In biochemical signaling pathways without explicit feedback connections, the core signal transduction is usually described as a one-way communication, going from upstream to downstream in a feedforward chain or network of covalent modification cycles. In this paper we explore the possibility of a new type of signaling called retroactive signaling, offered by the recently demonstrated property of retroactivity in signaling cascades. The possibility of retroactive signaling is analysed in the simplest case of the stationary states of a bicyclic cascade of signaling cycles. In this case, we work out the conditions for which variables of the upstream cycle are affected by a change of the total amount of protein in the downstream cycle, or by a variation of the phosphatase deactivating the same protein. Particularly, we predict the characteristic ranges of the downstream protein, or of the downstream phosphatase, for which a retroactive effect can be observed on the upstream cycle variables. Next, we extend the possibility of retroactive signaling in short but nonlinear signaling pathways involving a few covalent modification cycles.

Abstract 4925: The natural enzyme sequestration in signaling cascades provides inherent opportunities for off-target effects induced by kinase inhibitors

Off-target effects from targeted therapies are often attributed to cross-talk, which usually refers to inter-pathway molecular interactions that occur because of explicit regulatory feedback connections between two pathways. Recent experimental and theoretical studies have demonstrated, however, that covalently modified cascades naturally exhibit bidirectional signal propagation via a phenomenon termed retroactivity. This phenomenon arises due to enzyme sequestration where each cycle is coupled, not only to the next cycle, but also to the previous cycle. While retroactivity occurs naturally in covalently modified cascades, signaling pathways likely evolved to propagate information in a downstream manner. An important consequence of retroactivity, however, is that a downstream perturbation can induce an upstream response without the presence of regulatory feedback connections. We hypothesize that kinase inhibitors can produce off-target effects as a consequence of retroactivity alone via the following mechanism: a signal travels upstream from the site of a downstream perturbation through retroactivity and, upon reaching a shared upstream component, is delivered to an independent parallel pathway. To test the hypothesis we used a computational model to simulate the targeted inhibition of a specific kinase in a series signaling networks using physiologically and therapeutically relevant ranges for all parameters. Surprisingly, our results suggest that an off-target effect due to retroactive signaling is more likely when the first cycle in a non-inhibited cascade is “off” and not consuming large amounts of a shared up-stream activator. Our results also suggest that the kinetics governing covalently modified cycles in a cascade are more important for propagating an upstream off-target effect than the binding affinity of the drug to the targeted protein, which is a commonly optimized property in drug development. Finally, our results suggest that a single mutation has the capacity to produce a large spontaneous off-target effect without any direct regulatory connections between the targeted protein and the effected protein. Together, our results suggest that retroactivity may play an important role in the dysregulated signaling networks of cancer cells as well as the cellular response to targeted therapies. These findings have important implications for somatic evolution in cancer and the onset of therapeutic resistance, which has been widely reported for many targeted cancer therapeutics, including kinase inhibitors.

Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 103rd Annual Meeting of the American Association for Cancer Research; 2012 Mar 31-Apr 4; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2012;72(8 Suppl):Abstract nr 4925. doi:1538-7445.AM2012-4925

Response to Comment on “‘Load-Induced Modulation of Signal Transduction Networks’: Reconciling Ultrasensitivity with Bifunctionality?”

Complex biochemical and regulatory properties of a bifunctional enzyme mean that its activity cannot be modeled as a simple bifunctional system with distinct and reciprocally regulated states.

Load-induced modulation of signal transduction networks

Biological signal transduction networks are commonly viewed as circuits that pass along information--in the process amplifying signals, enhancing sensitivity, or performing other signal-processing tasks--to transcriptional and other components. Here, we report on a "reverse-causality" phenomenon, which we call load-induced modulation. Through a combination of analytical and experimental tools, we discovered that signaling was modulated, in a surprising way, by downstream targets that receive the signal and, in doing so, apply what in physics is called a load. Specifically, we found that non-intuitive changes in response dynamics occurred for a covalent modification cycle when load was present. Loading altered the response time of a system, depending on whether the activity of one of the enzymes was maximal and the other was operating at its minimal rate or whether both enzymes were operating at submaximal rates. These two conditions, which we call "limit regime" and "intermediate regime," were associated with increased or decreased response times, respectively. The bandwidth, the range of frequency in which the system can process information, decreased in the presence of load, suggesting that downstream targets participate in establishing a balance between noise-filtering capabilities and a circuit's ability to process high-frequency stimulation. Nodes in a signaling network are not independent relay devices, but rather are modulated by their downstream targets.

Kinase inhibitors can produce off-target effects and activate linked pathways by retroactivity

BACKGROUND

It has been shown in experimental and theoretical work that covalently modified signaling cascades naturally exhibit bidirectional signal propagation via a phenomenon known as retroactivity. An important consequence of retroactivity, which arises due to enzyme sequestration in covalently modified signaling cascades, is that a downstream perturbation can produce a response in a component upstream of the perturbation without the need for explicit feedback connections. Retroactivity may, therefore, play an important role in the cellular response to a targeted therapy. Kinase inhibitors are a class of targeted therapies designed to interfere with a specific kinase molecule in a dysregulated signaling pathway. While extremely promising as anti-cancer agents, kinase inhibitors may produce undesirable off-target effects by non-specific interactions or pathway cross-talk. We hypothesize that targeted therapies such as kinase inhibitors can produce off-target effects as a consequence of retroactivity alone.

RESULTS

We used a computational model and a series of simple signaling motifs to test the hypothesis. Our results indicate that within physiologically and therapeutically relevant ranges for all parameters, a targeted inhibitor can naturally induce an off-target effect via retroactivity. The kinetics governing covalent modification cycles in a signaling network were more important for propagating an upstream off-target effect in our models than the kinetics governing the targeted therapy itself. Our results also reveal the surprising and crucial result that kinase inhibitors have the capacity to turn "on" an otherwise "off" parallel cascade when two cascades share an upstream activator.

CONCLUSIONS

A proper and detailed characterization of a pathway's structure is important for identifying the optimal protein to target as well as what concentration of the targeted therapy is required to modulate the pathway in a safe and effective manner. We believe our results support the position that such characterizations should consider retroactivity as a robust potential source of off-target effects induced by kinase inhibitors and other targeted therapies.

p38γ promotes breast cancer cell motility and metastasis through regulation of RhoC GTPase, cytoskeletal architecture, and a novel leading edge behavior

Understanding the molecular alterations that confer cancer cells with motile, metastatic properties is needed to improve patient survival. Here, we report that p38γ motogen-activated protein kinase regulates breast cancer cell motility and metastasis, in part, by controlling expression of the metastasis-associated small GTPase RhoC. This p38γ-RhoC regulatory connection was mediated by a novel mechanism of modulating RhoC ubiquitination. This relationship persisted across multiple cell lines and in clinical breast cancer specimens. Using a computational mechanical model based on the finite element method, we showed that p38γ-mediated cytoskeletal changes are sufficient to control cell motility. This model predicted novel dynamics of leading edge actin protrusions, which were experimentally verified and established to be closely related to cell shape and cytoskeletal morphology. Clinical relevance was supported by evidence that elevated expression of p38γ is associated with lower overall survival of patients with breast cancer. Taken together, our results offer a detailed characterization of how p38γ contributes to breast cancer progression. Herein we present a new mechanics-based analysis of cell motility, and report on the discovery of a leading edge behavior in motile cells to accommodate modified cytoskeletal architecture. In summary, these findings not only identify a novel mechanism for regulating RhoC expression but also advance p38γ as a candidate therapeutic target.

Histone methyltransferase EZH2 induces Akt-dependent genomic instability and BRCA1 inhibition in breast cancer

Increased levels of EZH2, a critical regulator of cellular memory, signal the presence of metastasis and poor outcome in breast cancer patients. High levels of EZH2 are associated with nuclear pleomorphism, lack of estrogen receptor expression, and decreased nuclear levels of BRCA1 tumor suppressor protein in invasive breast carcinomas. The mechanism by which EZH2 overexpression promotes the growth of poorly differentiated invasive carcinomas remains to be defined. Here, we show that EZH2 controls the intracellular localization of BRCA1 protein. Conditional doxycycline-induced upregulation of EZH2 in benign mammary epithelial cells results in nuclear export of BRCA1 protein, aberrant mitoses with extra centrosomes, and genomic instability. EZH2 inhibition in CAL51 breast cancer cells induces BRCA1 nuclear localization and rescues defects in ploidy and mitosis. Mechanistically, EZH2 overexpression is sufficient for activation of the phosphoinositide 3-kinase/Akt (PI3K/Akt) pathway specifically through activation of Akt isoform 1. EZH2-induced BRCA1 nuclear export, aneuploidy, and mitotic defects were prevented by treatment with the PI3K inhibitors LY294002 or wortmannin. Targeted inhibition of Akt-1, Akt-2, and Akt-3 isoforms revealed that the EZH2-induced phenotype requires specific activation of Akt-1. The relevance of our studies to human breast cancer is highlighted by the finding that high EZH2 protein levels are associated with upregulated expression of phospho-Akt-1 (Ser473) and decreased nuclear expression of phospho-BRCA1 (Ser1423) in 39% of invasive breast carcinomas. These results enable us to pinpoint one mechanism by which EZH2 regulates BRCA1 expression and genomic stability mediated by the PI3K/Akt-1 pathway.

Abstract 1485: Identification of p38γ as a metastatic oncogene: Mechanism of action and clinical implications

Understanding the molecular alterations that confer metastatic properties to otherwise benign cells is essential to controlling breast cancer and significantly improving patient survival.

p38γ is a member of the p38 MAPK family and is normally expressed predominantly in muscle tissue. Developmentally, the primary role of p38γ is to promote myoblast differentiation into myotubes. Viewing the function of p38γ through the perspective of its developmental role led us to hypothesize that p38γ is involved in enabling mesenchymal-like behavior in breast cancer cells by controlling their motility properties.

To elucidate the role of p38γ in breast cancer progression we used RNAi to knock down (KD) p38γ in three aggressive breast cancer cell lines. We then subjected the KD cells to a panel of in vitro assays targeting the hallmarks of cancer progression. We investigated in detail the relationship between cytoskeletal defects and motility in p38γ KD cells using an innovative fusion of mathematical modeling and experimental cell biology. To evaluate the in vivo relevance of our in vitro observations, stable p38γ KD cells were orthotopically xenografted into athymic nude mice. Finally, we analyzed p38γ and RhoC expression in over 300 clinical breast cancer samples.

We found that p38γ phosphorylation is elevated in all three breast cancer cell lines compared to two control cell lines. p38γ KD also significantly impairs cell motility and invasion while increasing proliferation, and also drastically alters actin cytoskeletal architecture. By developing an in silico mechanical model of cell motility we show that the modified cytoskeletal architecture of p38γ KD cells is a driving force behind the observed altered motility. Our analysis also uncovered a novel leading edge behavior used by both p38γ KD and control cells to accommodate varying cytoskeletal architectures. Biochemically, p38γ concurrently alters RhoC expression, and rescuing RhoC expression restores aggressiveness to p38γ KD cells. When xenografted, p38γ KD cells are significantly less metastatic than scrambled control cells despite forming larger primary tumors. p38γ and RhoC expression correlate across clinical breast cancer specimens and may be predictive of disease severity.

We demonstrate for the first time that p38γ is a metastatic oncogene that specifically acts as a switch between primary tumor proliferation and metastatic cell motility and thus presumably between ductal carcinoma in situ and invasive breast cancer. p38γ elicits its effects, at least in part, through cytoskeletal remodeling and by modulating expression of another known metastatic oncogene, RhoC. These findings have important therapeutic implications as p38γ and RhoC expression were found to be concurrently altered in patient tissue.

Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 102nd Annual Meeting of the American Association for Cancer Research; 2011 Apr 2-6; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2011;71(8 Suppl):Abstract nr 1485. doi:10.1158/1538-7445.AM2011-1485

Abstract 4907: Elucidating the role of retroactive signaling and kinase inhibitors on off-target drug effects

The primary objective of targeted cancer therapies is to modulate cancer progression by perturbing specific molecules involved in aberrant proliferation and invasion. Kinase inhibitors are targeted therapies which are designed to interfere with a specific kinase molecule in a dysregulated oncogenic signaling cascade. While extremely promising as anti-cancer agents, such inhibitors may have undesirable off-target effects, whether by non-specific interactions or by effects from pathway cross-talk. We have shown in published experimental and theoretical work that covalently modified signaling cascades naturally exhibit bidirectional signal propagation. This phenomenon is termed retroactivity and challenges the widespread notion that information in cascades only flows from the cell surface to the nucleus. Previous work has demonstrated that increasing the concentration of a phosphatase in the terminal cycle of a covalently modified cascade may result in a measurable decrease in the concentration of the previous cycle's activated kinase. Thus, a downstream perturbation in a signaling cascade can produce a reverse (or retroactive) response without the need for direct negative feedback connections. This led us to hypothesize that the use of an inhibitory drug in a signaling network may cause an upstream off-target effect simply by inhibiting the activation or deactivation of a downstream kinase. To test the hypothesis that retroactivity contributes to off-target effects, we extended our previous work to a computational model that tested a series of signaling networks. The objective of our approach was two-fold: (1) to probe the effect of retroactivity on a kinase inhibitor in a signaling network and (2) to test whether retroactivity is likely to produce a measurable off-target effect under physiologically realistic conditions. Specifically, our model simulates the targeted inhibition of an activated kinase in a series of multi-cycle networks. The results of our work indicate that at physiologically and therapeutically relevant concentrations, a targeted inhibitor may induce a measurable off-target effect via retroactivity. We also performed local sensitivity analyses to predict the kinetic parameters that most affect the off-target response. Surprisingly, the drug disassociation constant is predicted to have very little effect while parameters such as the enzyme saturation and maximum velocity of some cycles are predicted to be very important. A proper characterization of a pathway's structure is important for identifying which protein in the pathway represents the optimal drug target as well as what concentration of the targeted therapy is likely to modulate the pathway in the manner desired. We believe our results support the position that such characterizations should consider the role of retroactivity as a source of a potential off-target effects by kinase inhibitors.

Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 102nd Annual Meeting of the American Association for Cancer Research; 2011 Apr 2-6; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2011;71(8 Suppl):Abstract nr 4907. doi:10.1158/1538-7445.AM2011-4907

Multisite phosphorylation provides an effective and flexible mechanism for switch-like protein degradation

Phosphorylation-triggered degradation is a common strategy for elimination of regulatory proteins in many important cell signaling processes. Interesting examples include cyclin-dependent kinase inhibitors such as p27 in human and Sic1 in yeast, which play crucial roles during the G1/S transition in the cell cycle. In this work, we have modeled and analyzed the dynamics of multisite-phosphorylation-triggered protein degradation systematically. Inspired by experimental observations on the Sic1 protein and a previous intriguing theoretical conjecture, we develop a model to examine in detail the degradation dynamics of a protein featuring multiple phosphorylation sites and a threshold site number for elimination in response to a kinase signal. Our model explains the role of multiple phosphorylation sites, compared to a single site, in the regulation of protein degradation. A single-site protein cannot convert a graded input of kinase increase to much sharper output, whereas multisite phosphorylation is capable of generating a highly switch-like temporal profile of the substrate protein with two characteristics: a temporal threshold and rapid decrease beyond the threshold. We introduce a measure termed temporal response coefficient to quantify the extent to which a response in the time domain is switch-like and further investigate how this property is determined by various factors including the kinase input, the total number of sites, the threshold site number for elimination, the order of phosphorylation, the kinetic parameters, and site preference. Some interesting and experimentally verifiable predictions include that the non-degradable fraction of the substrate protein exhibits a more switch-like temporal profile; a sequential system is more switch-like, while a random system has the advantage of increased robustness; all the parameters, including the total number of sites, the threshold site number for elimination and the kinetic parameters synergistically determine the exact extent to which the degradation profile is switch-like. Our results suggest design principles for protein degradation switches which might be a widespread mechanism for precise regulation of cellular processes such as cell cycle progression.

From in vitro to in silico and back again: using biological and mathematical synergy to decipher breast cancer cell motility

The complexity of biological systems is often prohibitive in testing specific hypotheses from first physical principles. To circumvent these limitations we used biological data to inform a mathematical model of breast cancer cell motility. Using this in silico model we were able to accurately assess the influence of actin cytoskeletal architecture on the motility of a genetically modified breast cancer cell line. Furthermore, using the in silico model revealed a biological phenomenon that has not been previously described in live cell movement. Fusing biology and mathematics as presented here represents a new direction for biomedical research in which advances in each field synergistically drive discoveries in the other.

Signaling properties of a covalent modification cycle are altered by a downstream target

We used a model system of purified components to explore the effects of a downstream target on the signaling properties of a covalent modification cycle, an example of retroactivity. In the experimental system used, a bifunctional enzyme catalyzed the modification and demodification of its substrate protein, with both activities regulated by a small molecule stimulus. Here we examined how a downstream target for one or both forms of the substrate of the covalent modification cycle affected the steady-state output of the system, the sensitivity of the response to the stimulus, and the concentration of the stimulus required to provide the half-maximal response (S(50)). When both the modified and unmodified forms of the substrate protein were sequestered by the downstream target, the sensitivity of the response was dramatically decreased, but the S(50) was only modestly affected. Conversely, when the downstream target only sequestered the unmodified form of the substrate protein, significant effects were observed on both system sensitivity and S(50). Behaviors of the experimental systems were well approximated both by simple models allowing analytical solutions and by a detailed model based on the known interactions and enzymatic activities. Modeling and experimentation indicated that retroactivity may result in subsensitive responses, even if the covalent modification cycle displays significant ultrasensitivity in the absence of retroactivity. Thus, we provide examples of how a downstream target can alter the signaling properties of an upstream signal transduction covalent modification cycle.

Abstract 3174: Visualization of the activation/deactivation cycle of RhoC in inflammatory breast cancer

A key event in the development of Inflammatory Breast Cancer, a particularly aggressive, metastatic form of breast cancer, is the over-expression of the GTPase protein RhoC. RhoC plays a role in regulating cell shape, attachment and motility and over-expression is likely associated with tumor proliferation. RhoC is activated into a GTP-bound state by the regulatory proteins GEFs (guanine nucleotide exchange factor proteins) and is deactivated to a GDP-bound state by GAPs (GTPase activating proteins). Stimulation of IBC cells with LPA (lysophophatidic acid), an activator of the cycle, showed a transient RhoC-GTP increase, peaking at 2 minutes and then diminishing until 20 minutes. It has been suggested that LPA can activate both GEF and GAP proteins, and that a delay in the GAP protein activation could explain this RhoC-GTP behavior. Because RhoC activity occurs at the plasma membrane, the translocation of GEF and GAP proteins was studied as an indication of RhoC interaction. Immunofluorescent stains for the three key proteins_RhoC, GAP, and GEF_were performed on an IBC cell line before and after stimulation with LPA and observed using confocal microscopy. Preliminary images revealed a co-localization of three proteins RhoC, p190B (a GAP protein), and PDZ (a GEF protein) in the cytosol, forming a gradient of high to low protein concentration from the nuclear membrane to the plasma membrane. Further analysis with ImageJ software showed slight differences in protein concentration at the two membranes upon stimulation with LPA. RhoC protein levels increased at both membranes, PDZ quantities increased slightly at the plasma membrane, and the p190B levels decreased at both membranes as a result of stimulation. This data suggests that after 5 minutes, LPA has a positive regulatory effect on the GEF protein and a negative one on the GAP protein at the plasma membrane. In addition, co-localization of the proteins RhoC and p190B was analyzed at different time points. The data indicates a decrease in co-localization after LPA stimulation that is slightly recovered 20 minutes after treatment. A possible explanation for these results is that 5 minutes after stimulation, RhoC undergoes translocation to the nuclear and plasma membranes, and is not followed by the GAP protein until later times. As of now, these conclusions are preliminary. These results will be verified using a new RhoC antibody, and additional co-localization analysis.

Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 101st Annual Meeting of the American Association for Cancer Research; 2010 Apr 17-21; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2010;70(8 Suppl):Abstract nr 3174.

Abstract 5193: p38 gamma is a novel driver of breast cancer metastasis by modulating cellular motility

Breast cancer is the second leading cause of cancer-related deaths in women in the US each year. The vast majority of these fatalities are not caused by primary tumor burden but rather by metastases to vital organs. The clinical shift from localized to metastatic breast cancer entails a requirement that cancer cells activate an invasive program and be able to adapt to changing extracellular stimuli. The p38 mitogen activated protein kinase (MAPK) pathway represents a potential signaling switch for the transition from primary to metastatic cancer. p38 is a member of the MAPK family of stress and mitogen-responsive protein kinases and consists of four closely related isoforms: alpha, beta, gamma, and delta. p38 serves as a major signaling hub in the cell, integrating signals from a variety of signaling pathways and channeling these stimuli into cellular responses through an array of effector proteins. The four isoforms show unique expression patterns in normal tissue: alpha, beta, and delta are the most ubiquitously expressed, while gamma is restricted primarily to skeletal muscle. Surprisingly, we found high levels of phosphorylated p38 gamma in a myoepithelial-derived aggressive breast cancer cell line, MDA-MB-231, compared to the non-tumorigenic mammary epithelial cell line MCF-10A, despite both cell lines having similar levels of total p38 gamma. This led us to hypothesize that, based on its normal function in skeletal muscle, p38 gamma may play a role in mediating metastatic behavior in basally-derived breast cancers. To test this hypothesis, we used isoform-specific shRNA to knock down p38 gamma expression in MDA-MB-231. Inhibition of p38 gamma (referred to as shGamma) resulted in decreased cell motility and invasion, two hallmarks of cancer metastasis. Analysis of individual cells revealed a change in the mode of motility from mesenchymal-like in control cells to lamellipodia-driven motility in the shGamma cells. This altered motility appears to be the result of a dramatically modified actin cytoskeleton, characterized by loss of lengthwise stress fibers and bundled, rather than branched, actin in the lamellipodia. When we investigated gene expression changes responsible for this phenotype we discovered that expression of the metastatic oncogene RhoC GTPase is significantly reduced at the protein, but not mRNA, level in shGamma cells due to decreased RhoC protein stability. These findings have great potential clinical significance, as early evidence indicates that p38 gamma is overexpressed in cases of invasive breast cancer, and appears to correlate with RhoC expression in these samples. Our current and future work includes determining which breast cancer subtypes utilize p38 gamma for invasion and metastasis, and characterizing in more detail the dynamic actin cytoskeletal changes affected by p38 gamma.

Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 101st Annual Meeting of the American Association for Cancer Research; 2010 Apr 17-21; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2010;70(8 Suppl):Abstract nr 5193.

Quantifying calcium fluxes underlying calcium puffs in Xenopus laevis oocytes

We determine the calcium fluxes through inositol 1,4,5-trisphosphate receptor/channels underlying calcium puffs of Xenopus laevis oocytes using a simplified version of the algorithm of Ventura et al. An analysis of 130 puffs obtained with Fluo-4 indicates that Ca2+ release comes from a region of width approximately 450 nm, that the release duration is peaked around 18 s and that the underlying Ca2+ currents range between 0.12 and 0.95 pA. All these parameters are independent of IP(3) concentration. We explore what distributions of channels that open during a puff, N(p), and what relations between current and number of open channels, I(N(p)), are compatible with our findings and with the distribution of puff-to-trigger amplitude ratio reported in Rose et al. To this end, we use simple "mean field" models in which all channels open and close simultaneously. We find that the variability among clusters plays an important role in shaping the observed puff amplitude distribution and that a model for which I(N(p)) approximately N(p) for small N(p) and I(N(p)) approximately N(p)(1/alpha) (alpha > 1) for large N(p), provides the best agreement. Simulations of more detailed models in which channels open and close stochastically show that this nonlinear behavior can be attributed to the limited time resolution of the observations and to the averaging procedure that is implicit in the mean-field models. These conclusions are also compatible with observations of approximately 400 puffs obtained using the dye Oregon green.

Preclinical development of a bifunctional cancer cell homing, PKCepsilon inhibitory peptide for the treatment of head and neck cancer

Head and neck squamous cell carcinoma (HNSCC) is the sixth most frequent cancer worldwide, comprising approximately 50% of all malignancies in some developing nations. Our recent work identified protein kinase Cepsilon (PKCepsilon) as a critical and causative player in establishing an aggressive phenotype in HNSCC. In this study, we investigated the specificity and efficacy of HN1-PKCepsilon, a novel bifunctional cancer cell homing, PKCepsilon inhibitory peptide, as a treatment for HNSCC. HN1-PKCepsilon peptide was designed by merging two separate technologies and synthesized as a capped peptide with two functional modules, HN1 (cancer cell homing) and PKCepsilon (specific PKCepsilon inhibitory), connected by a novel linker module. HN1-PKCepsilon preferentially internalized into UMSCC1 and UMSCC36 cells, two HNSCC cell lines, in comparison with oral epithelial cells: 82.1% positive for UMSCC1 and 86.5% positive for UMSCC36 compared with 1.2% positive for oral epithelial cells. In addition, HN1-PKCepsilon penetrated HNSCC cells in a dose- and time-dependent manner. Consistent with these in vitro observations, systemic injection of HN1-PKCepsilon resulted in selective delivery of HN1-PKCepsilon into UMSCC1 xenografts in nude mice. HN1-PKCepsilon blocked the translocation of active PKCepsilon in UMSCC1 cells, confirming HN1-PKCepsilon as a PKCepsilon inhibitor. HN1-PKCepsilon inhibited cell invasion by 72 +/- 2% (P < 0.001, n = 12) and cell motility by 56 +/- 2% (P < 0.001, n = 5) in UMSCC1 cells. Moreover, in vivo bioluminescence imaging showed that HN1-PKCepsilon significantly (83 +/- 1% inhibition; P < 0.02) retards the growth of UMSCC1 xenografts in nude mice. Our work indicates that the bifunctional HN1-PKCepsilon inhibitory peptide represents a promising novel therapeutic strategy for HNSCC.

On the role of cell signaling models in cancer research

The main objective of this review is to emphasize the role and importance of the careful mathematical/computational modeling of signaling networks for the understanding of aberrant signaling in cancer and for the development of targeted therapies.

A hidden feedback in signaling cascades is revealed

Cycles involving covalent modification of proteins are key components of the intracellular signaling machinery. Each cycle is comprised of two interconvertable forms of a particular protein. A classic signaling pathway is structured by a chain or cascade of basic cycle units in such a way that the activated protein in one cycle promotes the activation of the next protein in the chain, and so on. Starting from a mechanistic kinetic description and using a careful perturbation analysis, we have derived, to our knowledge for the first time, a consistent approximation of the chain with one variable per cycle. The model we derive is distinct from the one that has been in use in the literature for several years, which is a phenomenological extension of the Goldbeter-Koshland biochemical switch. Even though much has been done regarding the mathematical modeling of these systems, our contribution fills a gap between existing models and, in doing so, we have unveiled critical new properties of this type of signaling cascades. A key feature of our new model is that a negative feedback emerges naturally, exerted between each cycle and its predecessor. Due to this negative feedback, the system displays damped temporal oscillations under constant stimulation and, most important, propagates perturbations both forwards and backwards. This last attribute challenges the widespread notion of unidirectionality in signaling cascades. Concrete examples of applications to MAPK cascades are discussed. All these properties are shared by the complete mechanistic description and our simplified model, but not by previously derived phenomenological models of signaling cascades.

Genetic Determinants of Aggressive Breast Cancer

The development and spread of breast and other human cancers are caused by the overexpression, mutation, and/or deletion of specific genes that drive these events. Finding genetic and molecular differences between cancerous and healthy cells can reveal the genetic determinants of cancer. This knowledge results in a better understanding of the carcinogenic process and improved predictive power, with implications for identifying new drug targets, designing novel therapies, and improving preclinical and clinical studies. We review the concepts of biomarker, genetic marker, and genetic determinant in cancer, with particular focus on the most aggressive and lethal form of breast cancer, termed inflammatory breast cancer (IBC). Using IBC as an example, we describe in detail the approaches to identify the genes that are responsible for-and not merely associated with-this disease.

Simple data-driven models of intracellular calcium dynamics with predictive power

Biology is complex. However, it is not clear how much of this complexity must necessarily translate into complicated mathematical models of biological processes. Simple models can be appealing to physicists but are usually deceiving for biologists. Complicated models, on the other hand, depend on too many parameters whose values are frequently unknown. Therefore, complicated models, although in principle more realistic, can lead to erroneous results if they are sensitive to these unknown parameter values. Intracellular calcium signals provide an example of utmost biological importance in which the issue of "simple vs complex" can be explored. In this paper we show that simple models describing the dynamics of intracellular calcium can be directly inferred from experimental data, without no a priori information on unknown parameters. A similar approach can be followed to study other reaction-diffusion systems. In spite of their simplicity, these models can provide quantitative information on some of the processes that shape calcium signals, such as the calcium current that underlies an experimental observation. This shows that simple models of biological systems are not limited to qualitative descriptions.

Calcium oscillations and waves generated by multiple release mechanisms in pancreatic acinar cells

We explore the dynamic behavior of a model of calcium oscillations and wave propagation in the basal region of pancreatic acinar cells [Sneyd, J., et al., Biophys. J. 85: 1392-1405, 2003]. Since it is known that two principal calcium release pathways are involved, inositol trisphosphate receptors (IPR) and ryanodine receptors (RyR), we study how the model behavior depends on the density of each receptor type. Calcium oscillations can be mediated either by IPR or RyR. Continuous increases in either RyR or IPR density can lead to the appearance and disappearance of oscillations multiple times, and the two receptor types interact via their common effect on cytoplasmic calcium concentration and the subsequent effect on the total amount of calcium inside the cell. Increases in agonist concentration can stimulate oscillations via the RyR by increasing calcium influx. Using a two time-scale approach, we explain these complex behaviors by treating the total amount of cellular calcium as a slow parameter. Oscillations are controlled by the shape of the slow manifold and where it intersects the nullcline of the slow variable. When calcium diffusion is included, the existence of traveling waves in the model equation is strongly dependent on the interplay between the total amount of calcium in the cell and membrane transport, a feature that can be experimentally tested. Our results help us understand the behavior of a model that includes both receptors in comparison to the properties of each receptor type in isolation.

A model-independent algorithm to derive Ca2+ fluxes underlying local cytosolic Ca2+ transients

Local intracellular Ca(2+) signals result from Ca(2+) flux into the cytosol through individual channels or clusters of channels. To gain a mechanistic understanding of these events we need to know the magnitude and spatial distribution of the underlying Ca(2+) flux. However, this is difficult to infer from fluorescence Ca(2+) images because the distribution of Ca(2+)-bound dye is affected by poorly characterized processes including diffusion of Ca(2+) ions, their binding to mobile and immobile buffers, and sequestration by Ca(2+) pumps. Several methods have previously been proposed to derive Ca(2+) flux from fluorescence images, but all require explicit knowledge or assumptions regarding these processes. We now present a novel algorithm that requires few assumptions and is largely model-independent. By testing the algorithm with both numerically generated image data and experimental images of sparklets resulting from Ca(2+) flux through individual voltage-gated channels, we show that it satisfactorily reconstructs the magnitude and time course of the underlying Ca(2+) currents.

Saltatory and continuous calcium waves and the rapid buffering approximation

Calcium waves propagate inside cells due to a regenerative mechanism known as calcium-induced calcium release. Buffer-mediated calcium diffusion in the cytosol plays a crucial role in the process. However, most models of calcium waves either treat buffers phenomenologically or assume that they are in equilibrium with calcium (the rapid buffering approximation). In this article we address the issue of whether this approximation provides a good description of wave propagation. We first compare the timescales present in the problem, and determine the situations in which the equilibrium hypothesis fails. We then present a series of numerical studies based on the simple fire-diffuse-fire model of wave propagation. We find that the differences between the full and reduced descriptions may lead to errors that are above experimental resolution even for relatively fast buffers in the case of saltatory waves. Conversely, in the case of continuous waves, the approximation may give accurate results even for relatively slow buffers.

Generic two-variable model of excitability

We present a simple model that displays all classes of two-dimensional excitable regimes. One of the variables of the model displays the usual spikes observed in excitable systems. Since the model is written in terms of a "standard" vector field, it is always possible to fit it to experimental data displaying spikes in an algorithmic way. In fact, we use it to fit a series of membrane potential recordings obtained in the medicinal leech and time series generated with the FitzHugh-Nagumo equations and the excitability model of Eguía et al. [Phys. Rev. E 58, 2636 (1998)]. In each case, we determine the excitability class of the corresponding system.

Pentoxifylline influences the autocrine function of organ cultured donor corneas and enhances endothelial cell survival

BACKGROUND/AIMS

Scientific interest in pentoxifylline has been reawakened owing to the recognised effects of this drug on immune functions, particularly its influence on cytokine production. In a previous study, the authors demonstrated that spiking of organ culture media with endotoxin elicited a marked enhancement in the release of IL-6 and IL-8 from corneal tissue and that these events coincided with degenerative changes in endothelial cells and a higher incidence of actual loss among this population. Since traces of donor derived endotoxin can be detected in up to 50% of corneal organ cultures, this substance may have a direct influence on graft viability or trigger inflammatory responses in the host. They, therefore, wished to ascertain whether supplementation of media with pentoxifylline improved endothelial cell survival in organ cultured donor corneas.

METHODS

12 fellow pairs of donor corneas were cultured for 20 days, with a change of medium on day 10: One of each pair was incubated in the absence, and the other in the presence, of pentoxifylline (25 microg/ml). Samples of medium were withdrawn at regular intervals during the course of incubation and screened for cytokines IL-6, IL-8, and prostaglandin E2 by ELISA. Endothelial cell morphology and numerical density were assessed on days 0, 10 and 20.

RESULTS

Addition of pentoxifylline to organ culture media led to a significant improvement in endothelial cell survival. This drug also elicited a significant increase in the level of IL-6 and marginally suppressed that of IL-8 during the initial 10 day phase of incubation. During the second 10-20 day phase, the level of both IL-6 and IL-8 decreased significantly in the presence of pentoxifylline, the relation between these two cytokines being the inverse of that observed in the absence of the drug. No significant changes in the level of prostaglandin E2 were apparent.

CONCLUSION

The addition of pentoxifylline to organ culture media leads, ultimately, to a suppression of IL-6 and IL-8 secretion by corneal tissue. The potentially damaging effects of these cytokines are thereby quelled, as evidenced by the improvement in endothelial cell survival.

Central corneal thickness measurements in patients with normal tension glaucoma, primary open angle glaucoma, pseudoexfoliation glaucoma, or ocular hypertension

BACKGROUND/AIMS

Recent studies have revealed patients with ocular hypertension to have thicker than normal central corneas and those with normal tension glaucoma to have thinner than normal ones, as determined by ultrasonic pachymetry. Since corneal thickness measurements and applanation tonometric estimates of intraocular pressure (IOP) correlate positively, monitoring of the former parameter have served as the basis for adjusting readings pertaining to the latter, with the consequence that many patients have had to be reclassified. With a view to validating these pachymetric studies, the central corneal thickness was determined in patients with normal tension glaucoma, primary open angle glaucoma, pseudoexfoliation glaucoma, or ocular hypertension, as well as that of normal subjects, using optical low coherence reflectometry, which is a new and more precise method than ultrasonic pachymetry.

METHODS

34 patients with normal tension glaucoma, 20 with primary open angle glaucoma, 13 with pseudoexfoliation glaucoma, and 12 with ocular hypertension, together with 21 control subjects, were included in this observational, concurrent case-control study. One eye per individual was randomly selected for investigation. IOP was measured by Goldmann applanation tonometry and central corneal thickness by optical low coherence reflectometry.

RESULTS

Central corneal thickness was significantly higher (p < or =0.001) in patients with ocular hypertension than in normal individuals or in subjects with either normal tension glaucoma, primary open angle glaucoma, or pseudoexfoliation glaucoma, there being no significant differences between the latter four groups. Patients with ocular hypertension were also significantly younger (p < or =0.003) than those within any of the three glaucomatous groups.

CONCLUSION

This study confirms that a significant number of patients with ocular hypertension have normal IOPs after the appropriate adjustments have been made for deviations from normal in their central corneal thickness. The accurate measurement of this latter parameter is important not only for individual patient care, in permitting more precise estimations of IOP, but also for clinical studies, in assuring a more reliable classification of subjects.

Corneal thickness and endothelial density before and after cataract surgery

BACKGROUND/AIMS

Deturgescence of the corneal stroma is controlled by the pumping action of the endothelial layer and can be monitored by measurement of central corneal thickness (pachymetry). Loss or damage of endothelial cells leads to an increase in corneal thickness, which may ultimately induce corneal decompensation and loss of vision. Little is known about the effect of moderate reductions in endothelial cell number on the thickness of the corneal stroma. This study aimed to investigate this matter further using patients who had incurred moderate decreases in their endothelial cell counts as a result of cataract surgery.

METHODS

Central corneal thickness was measured 1 day before surgery, 1 day after surgery, and again at 3 months or 1 year. Endothelial cell counts were also performed 1 day before surgery and thereafter at 3 months or 1 year after surgery. The relationship between these two parameters was assessed statistically. Precise measurements of central corneal thickness were made by optical low coherence reflectometry. For comparative purposes, this parameter was also determined by ultrasonic pachymetry. Central corneal endothelial cell numerical density was estimated on photomicrographs taken with a specular microscope.

RESULTS

All patients had significant postoperative corneal swelling on the day after surgery; preoperative values were restored by 3 and 12 months, even though significant endothelial cell losses had occurred. No correlation existed between central corneal thickness and central corneal endothelial cell numerical density. Measurements estimated by ultrasonic pachymetry were more variable and significantly higher than those determined by optical low coherence reflectometry.

CONCLUSION

As long as the numerical density of the corneal endothelial cells does not fall below the physiological threshold, a moderate decrease in this parameter does not compromise the pumping activity of the layer as a whole.

Toxicity of pentoxifylline on monolayers of highly proliferative cells of epithelial origin

Interest in pentoxifylline has been recently reawakened owing to its suppressive effect on cell cytokine production. In this capacity, it may be of value as a routine supplement for culture media containing donor corneas. The purpose of the present study was to evaluate the toxic effects of pentoxifylline on two standardized cell lines of epithelial origin. Vero and Chang cells were incubated with various concentrations of pentoxifylline. Acute toxicity (4 hr) was assessed by monitoring the permeability of cells to propidium iodide; chronic toxicity (7 days) was determined by monitoring the effect of pentoxifylline on esterase activity and cell proliferation. The viability of cells was also assessed by microscopic inspection. Signs of acute toxicity became manifest at a pentoxifylline concentration of 100 mg/l in both Chang and Vero cells. Indications of chronic toxicity were observed at a drug concentration of 10 mg/l in Chang cells but at 1 mg/l in Vero ones. Proliferation was suppressed at pentoxifylline concentrations of 100 mg/l and 10 mg/l in Chang and Vero cells, respectively. Degenerative morphological changes were observed at a drug concentration of 100 mg/l in both cell types. At a concentration of 0.1 mg/l, pentoxifylline elicited no signs of acute or chronic toxicity in either Chang or Vero cells. At this dose, the drug is therefore unlikely to have deleterious effects on cultured donor corneas.

Detection of endotoxin in media from sterile corneal organ cultures

Lipopolysaccharides (LPS) are elements of the cell wall of gram-negative bacteria. They are also called endotoxins and are known to stimulate various inflammatory reactions by interaction with cytokines and macrophages. Additionally, they may have a direct effect on other cells. As contamination of sterile organ-culture media with bacterial substances may influence the donor-tissue prognosis, we investigated a series of culture media drawn from organ culture for the presence of endotoxin. A total of 341 samples of sterile organ-culture media from 3 different cornea banks using either organ culture or the refrigerated-storage technique were tested for endotoxin. The assay was performed using the Limulus amebocyte-lysate test. A level of endotoxin above the background threshold was found in 99/341 (29%) culture media. The incidence of endotoxin ranged from 14% up to 50%, depending on the cornea bank and culture system used. Endotoxin detected in sterile corneal organ cultures probably derives from nonreplicating bacterial postmortem donor-tissue contamination. The presence of endotoxin-positive cultures varied between the eye banks but was not related to a given storage method. As endotoxin may directly influence graft viability or trigger inflammatory host responses, these findings may have significance for the clinical results of corneal grafting.

[HLA typing of donor corneas with extended post mortem time]

In case of high-risk patients cornea transplantation should be carried out using HLA-matched donor corneas to minimize the risk of rejection. HLA typing using blood lymphocytes of the donor is impossible because of the long post-mortem times. Alternatively, HLA typing can be performed using retinal pigment epithelium (RPE). Nevertheless, this method is limited by increasing post-mortem times. The aim of the present study was the optimization of culture conditions for RPE cells isolated from donor eyes with long post-mortem times. The HLA type should be evaluated during an acceptable period of organ preservation of the corresponding cornea. In different steps the method for isolation of the cells was optimized and a growth medium for RPE cells was established. Various supplements, including uvea-conditioned medium, were assessed using growth assays. The optimization of culture conditions led to an increase in the estimation of the complete HLA I and HLA II antigens from 36% to 74%. The time needed for the typing-procedure could be reduced from 38 to 17 days (average). At present 30% of the donor tissue with long post-mortem times can be typed in less than 14 days.