A systems model for immune cell interactions unravels the mechanism of inflammation in human skin.

Inflammation is characterized by altered cytokine levels produced by cell populations in a highly interdependent manner. To elucidate the mechanism of an inflammatory reaction, we have developed a mathematical model for immune cell interactions via the specific, dose-dependent cytokine production rates of cell populations. The model describes the criteria required for normal and pathological immune system responses and suggests that alterations in the cytokine production rates can lead to various stable levels which manifest themselves in different disease phenotypes. The model predicts that pairs of interacting immune cell populations can maintain homeostatic and elevated extracellular cytokine concentration levels, enabling them to operate as an immune system switch. The concept described here is developed in the context of psoriasis, an immune-mediated disease, but it can also offer mechanistic insights into other inflammatory pathologies as it explains how interactions between immune cell populations can lead to disease phenotypes.

Multiple calcium binding sites make calmodulin multifunctional.

Protein-protein or protein-ion interactions with multisite proteins are essential to the regulation of intracellular and extracellular events. There is, however, limited understanding of how ligand-multisite protein interactions selectively regulate the activities of multiple protein targets. In this paper, we focus on the important calcium (Ca(2+)) binding protein calmodulin (CaM), which has four Ca(2+) ion binding sites and regulates the activity of over 30 other proteins. Recent progress in structural studies has led to significant improvements in the understanding of Ca(2+)-CaM-dependent regulation mechanisms. However, no quantitative model is currently available that can fully explain how the structural diversity of protein interaction surfaces leads to selective activation of protein targets. In this paper, we analyze the multisite protein-ligand binding mechanism using mathematical modelling and experimental data for Ca(2+)-CaM-dependent protein targets. Our study suggests a potential mechanism for selective and differential activation of Ca(2+)-CaM targets by the same CaM molecules, which are involved in a variety of intracellular functions. The close agreement between model predictions and experimental dose-response curves for CaM targets available in the literature suggests that such activation is due to the selective activity of CaM conformations in complexes with variable numbers of Ca(2+) ions. Although the paper focuses on the Ca(2+)-CaM pair as a particularly data rich example, the proposed model predictions are quite general and can easily be extended to other multisite proteins. The results of the study may therefore be proposed as a general explanation for multifunctional target regulation by multisite proteins.

Structural and functional properties of the human notch-1 ligand binding region.

We present NMR structural and dynamics analysis of the putative ligand binding region of human Notch-1, comprising EGF-like domains 11-13. Functional integrity of an unglycosylated, recombinant fragment was confirmed by calcium-dependent binding of tetrameric complexes to ligand-expressing cells. EGF modules 11 and 12 adopt a well-defined, rod-like orientation rigidified by calcium. The interdomain tilt is similar to that found in previously studied calcium binding EGF pairs, but the angle of twist is significantly different. This leads to an extended double-stranded beta sheet structure, spanning the two EGF modules. Based on the conservation of residues involved in interdomain hydrophobic packing, we propose this arrangement to be prototypical of a distinct class of EGF linkages. On this premise, we have constructed a model of the 36 EGF modules of the Notch extracellular domain that enables predictions to be made about the general role of calcium binding to this region.

A systems model of phosphorylation for inflammatory signaling events.

Phosphorylation is a fundamental biochemical reaction that modulates protein activity in cells. While a single phosphorylation event is relatively easy to understand, multisite phosphorylation requires systems approaches for deeper elucidation of the underlying molecular mechanisms. In this paper we develop a mechanistic model for single- and multi-site phosphorylation. The proposed model is compared with previously reported studies. We compare the predictions of our model with experiments published in the literature in the context of inflammatory signaling events in order to provide a mechanistic description of the multisite phosphorylation-mediated regulation of Signal Transducer and Activator of Transcription 3 (STAT3) and Interferon Regulatory Factor 5 (IRF-5) proteins. The presented model makes crucial predictions for transcription factor phosphorylation events in the immune system. The model proposes potential mechanisms for T cell phenotype switching and production of cytokines. This study also provides a generic framework for the better understanding of a large number of multisite phosphorylation-regulated biochemical circuits.

Functional and topological properties in hepatocellular carcinoma transcriptome.

Hepatocellular carcinoma (HCC) is a leading cause of global cancer mortality. However, little is known about the precise molecular mechanisms involved in tumor formation and pathogenesis. The primary goal of this study was to elucidate genome-wide molecular networks involved in development of HCC with multiple etiologies by exploring high quality microarray data. We undertook a comparative network analysis across 264 human microarray profiles monitoring transcript changes in healthy liver, liver cirrhosis, and HCC with viral and alcoholic etiologies. Gene co-expression profiling was used to derive a consensus gene relevance network of HCC progression that consisted of 798 genes and 2,012 links. The HCC interactome was further confirmed to be phenotype-specific and non-random. Additionally, we confirmed that co-expressed genes are more likely to share biological function, but not sub-cellular localization. Analysis of individual HCC genes revealed that they are topologically central in a human protein-protein interaction network. We used quantitative RT-PCR in a cohort of normal liver tissue (n = 8), hepatitis C virus (HCV)-induced chronic liver disease (n = 9), and HCC (n = 7) to validate co-expressions of several well-connected genes, namely ASPM, CDKN3, NEK2, RACGAP1, and TOP2A. We show that HCC is a heterogeneous disorder, underpinned by complex cross talk between immune response, cell cycle, and mRNA translation pathways. Our work provides a systems-wide resource for deeper understanding of molecular mechanisms in HCC progression and may be used further to define novel targets for efficient treatment or diagnosis of this disease.

Aggregation pattern transitions by slightly varying the attractive/repulsive function.

Among collective behaviors of biological swarms and flocks, the attractive/repulsive (A/R) functional links between particles play an important role. By slightly changing the cutoff distance of the A/R function, a drastic transition between two distinct aggregation patterns is observed. More precisely, a large cutoff distance yields a liquid-like aggregation pattern where the particle density decreases monotonously from the inside to the outwards within each aggregated cluster. Conversely, a small cutoff distance produces a crystal-like aggregation pattern where the distance between each pair of neighboring particles remains constant. Significantly, there is an obvious spinodal in the variance curve of the inter-particle distances along the increasing cutoff distances, implying a legible transition pattern between the liquid-like and crystal-like aggregations. This work bridges the aggregation phenomena of physical particles and swarming of organisms in nature upon revealing some common mechanism behind them by slightly varying their inter-individual attractive/repulsive functions, and may find its potential engineering applications, for example, in the formation design of multi-robot systems and unmanned aerial vehicles (UAVs).

Computational modelling elucidates the mechanism of ciliary regulation in health and disease.

BACKGROUND: Ciliary dysfunction leads to a number of human pathologies, including primary ciliary dyskinesia, nephronophthisis, situs inversus pathology or infertility. The mechanism of cilia beating regulation is complex and despite extensive experimental characterization remains poorly understood. We develop a detailed systems model for calcium, membrane potential and cyclic nucleotide-dependent ciliary motility regulation. RESULTS: The model describes the intimate relationship between calcium and potassium ionic concentrations inside and outside of cilia with membrane voltage and, for the first time, describes a novel type of ciliary excitability which plays the major role in ciliary movement regulation. Our model describes a mechanism that allows ciliary excitation to be robust over a wide physiological range of extracellular ionic concentrations. The model predicts the existence of several dynamic modes of ciliary regulation, such as the generation of intraciliary Ca2+ spike with amplitude proportional to the degree of membrane depolarization, the ability to maintain stable oscillations, monostable multivibrator regimes, all of which are initiated by variability in ionic concentrations that translate into altered membrane voltage. CONCLUSIONS: Computational investigation of the model offers several new insights into the underlying molecular mechanisms of ciliary pathologies. According to our analysis, the reported dynamic regulatory modes can be a physiological reaction to alterations in the extracellular environment. However, modification of the dynamic modes, as a result of genetic mutations or environmental conditions, can cause a life threatening pathology.

Systems biology of cell behavior.

Systems Biology approaches to drug discovery largely focus on the increasing understanding of intracellular and cellular circuits, by computational representation of a molecular system followed by parameter validation against experimental data. This chapter outlines a universal approach to systems biology that allows the linking of intracellular molecular machinery and cellular activity. This procedure is achieved by applying mathematical modeling to molecular modules of a cell in the light of systems biology techniques.

Crosstalk between G-protein and Ca2+ pathways switches intracellular cAMP levels.

Cyclic adenosine monophosphate and cyclic guanosine monophosphate are universal intracellular messengers whose concentrations are regulated by molecular networks comprised of different isoforms of the synthases adenylate cyclase or guanylate cyclase and the phosphodiesterases which degrade these compounds. In this paper, we employ a systems biology approach to develop mathematical models of these networks that, for the first time, take into account the different biochemical properties of the isoforms involved. To investigate the mechanisms underlying the joint regulation of cAMP and cGMP, we apply our models to analyse the regulation of cilia beat frequency in Paramecium by Ca(2+). Based on our analysis of these models, we propose that the diversity of isoform combinations that occurs in living cells provides an explanation for the huge variety of intracellular processes that are dependent on these networks. The inclusion of both G-protein receptor and Ca(2+)-dependent regulation of AC in our models allows us to propose a new explanation for the switching properties of G-protein subunits involved in nucleotide regulation. Analysis of the models suggests that, depending on whether the G-protein subunit is bound to AC, Ca(2+) can either activate or inhibit AC in a concentration-dependent manner. The resulting analysis provides an explanation for previous experimental results that showed that alterations in Ca(2+) concentrations can either increase or decrease cilia beat frequency over particular Ca(2+) concentration ranges.

Computational modelling suggests dynamic interactions between Ca2+, IP3 and G protein-coupled modules are key to robust Dictyostelium aggregation.

Under conditions of starvation, Dictyostelium cells begin a programme of development during which they aggregate to form a multicellular structure by chemotaxis, guided by propagating waves of cyclic AMP that are relayed robustly from cell to cell. In this paper, we develop and analyse a new model for the intracellular and extracellular cAMP dependent processes that regulate Dictyostelium migration. The model allows, for the first time, a quantitative analysis of the dynamic interactions between calcium, IP(3) and G protein-dependent modules that are shown to be key to the generation of robust cAMP oscillations in Dictyostelium cells. The model provides a mechanistic explanation for the transient increase in cytosolic free Ca(2+) concentration seen in recent experiments with the application of the calmodulin inhibitor calmidazolium (R24571) to Dictyostelium cells, and also allows elucidation of the effects of varying both the conductivity of stretch-activated channels and the concentration of external phosphodiesterase on the oscillatory regime of an individual cell. A rigorous analysis of the robustness of the new model shows that interactions between the different modules significantly reduce the sensitivity of the resulting cAMP oscillations to variations in the kinetics of different Dictyostelium cells, an essential requirement for the generation of the spatially and temporally synchronised chemoattractant cAMP waves that guide Dictyostelium aggregation.

Electrostatic and functional analysis of the seven-bladed WD beta-propellers.

beta-propeller domains composed of WD repeats are highly ubiquitous and typically used as multi-site docking platforms to coordinate and integrate the activities of groups of proteins. Here, we have used extensive homology modelling of the WD40-repeat family of seven-bladed beta-propellers coupled with subsequent structural classification and clustering of these models to define subfamilies of beta-propellers with common structural, and probable, functional characteristics. We show that it is possible to assign seven-bladed WD beta-propeller proteins into functionally different groups based on the information gained from homology modelling. We examine general structural diversity within the WD40-repeat family of seven-bladed beta-propellers and demonstrate that seven-bladed beta-propellers composed of WD-repeats are structurally distinct from other seven-bladed beta-propellers. We further provide some insights into the multifunctional diversity of the seven-bladed WD beta-propeller surfaces. This report once again reinforces the importance of structural data and the usefulness of homology models in functional classification.

Elucidating the mechanisms of cooperative calcium-calmodulin interactions: a structural systems biology approach.

BACKGROUND: Calmodulin is an important multifunctional molecule that regulates the activities of a large number of proteins in the cell. Calcium binding induces conformational transitions in calmodulin that make it specifically active to particular target proteins. The precise mechanisms underlying calcium binding to calmodulin are still, however, quite poorly understood. RESULTS: In this study, we adopt a structural systems biology approach and develop a mathematical model to investigate various types of cooperative calcium-calmodulin interactions. We compare the predictions of our analysis with physiological dose-response curves taken from the literature, in order to provide a quantitative comparison of the effects of different mechanisms of cooperativity on calcium-calmodulin interactions. The results of our analysis reduce the gap between current understanding of intracellular calmodulin function at the structural level and physiological calcium-dependent calmodulin target activation experiments. CONCLUSION: Our model predicts that the specificity and selectivity of CaM target regulation is likely to be due to the following factors: variations in the target-specific Ca2+ dissociation and cooperatively effected dissociation constants, and variations in the number of Ca2+ ions required to bind CaM for target activation.

A calcium dependent de-adhesion mechanism regulates the direction and rate of cell migration: a mathematical model.

Cell migration has long been studied by a variety of techniques and many proteins have been implicated in its regulation. Integrins, key proteins that link the cell to the extracellular matrix, are central to adhesion complexes whose turnover defines the rate of cell locomotion. The formation and disassembly of these adhesions is regulated by both intracellular and extracellular factors. In this study we have focused on the Ca2+-dependent protein network (module) that disassembles the adhesion complexes. We have developed a mathematical model that includes the Ca2+-dependent enzymes micro-calpain and phospholipase C (PLC) as well as IP3 receptors and stretch activated Ca2+ channels, all of which have been reported to regulate migration. The model also considers the spatial effects of Ca2+ propagation into lamella. Our model predicts differential activation of calpain at the leading and trailing edges of the cell. Since disassembly of integrin adhesive contacts is proportional to the degree of calpain activation, this leads to cell migration in a preferred direction. We show how the dynamics of Ca2+ spiking affects calpain activation and thus changes the disassembly rate of adhesions. The spiking is controlled by PLC activity and currents through stretch-activated Ca2+ channels. Our model thus combines the effects of various molecular factors and leads to a consistent explanation of the regulation of the rate and direction of cell migration.

Solution structure and dynamics of a prototypical chordin-like cysteine-rich repeat (von Willebrand Factor type C module) from collagen IIA.

Chordin-like cysteine-rich (CR) repeats (also referred to as von Willebrand factor type C (VWC) modules) have been identified in approximately 200 extracellular matrix proteins. These repeats, named on the basis of amino acid conservation of 10 cysteine residues, have been shown to bind members of the transforming growth factor-beta (TGF-beta) superfamily and are proposed to regulate growth factor signaling. Here we describe the intramolecular disulfide bonding, solution structure, and dynamics of a prototypical chordin-like CR repeat from procollagen IIA (CR(ColIIA)), which has been previously shown to bind TGF-beta1 and bone morphogenetic protein-2. The CR(ColIIA) structure manifests a two sub-domain architecture tethered by a flexible linkage. Initial structures were calculated using RosettaNMR, a de novo prediction method, and final structure calculations were performed using CANDID within CYANA. The N-terminal region contains mainly beta-sheet and the C-terminal region is more irregular with the fold constrained by disulfide bonds. Mobility between the N- and C-terminal sub-domains on a fast timescale was confirmed using NMR relaxation measurements. We speculate that the mobility between the two sub-domains may decrease upon ligand binding. Structure and sequence comparisons have revealed an evolutionary relationship between the N-terminal sub-domain of the CR module and the fibronectin type 1 domain, suggesting that these domains share a common ancestry. Based on the previously reported mapping of fibronectin binding sites for vascular endothelial growth factor to regions containing fibronectin type 1 domains, we discuss the possibility that this structural homology might also have functional relevance.