Najl V. Valeyev

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.

Transcriptome classification reveals molecular subtypes in psoriasis

BackgroundPsoriasis is an immune-mediated disease characterised by chronically elevated pro-inflammatory cytokine levels, leading to aberrant keratinocyte proliferation and differentiation. Although certain clinical phenotypes, such as plaque psoriasis, are well defined, it is currently unclear whether there are molecular subtypes that might impact on prognosis or treatment outcomes.ResultsWe present a pipeline for patient stratification through a comprehensive analysis of gene expression in paired lesional and non-lesional psoriatic tissue samples, compared with controls, to establish differences in RNA expression patterns across all tissue types. Ensembles of decision tree predictors were employed to cluster psoriatic samples on the basis of gene expression patterns and reveal gene expression signatures that best discriminate molecular disease subtypes. This multi-stage procedure was applied to several published psoriasis studies and a comparison of gene expression patterns across datasets was performed.ConclusionOverall, classification of psoriasis gene expression patterns revealed distinct molecular sub-groups within the clinical phenotype of plaque psoriasis. Enrichment for TGFb and ErbB signaling pathways, noted in one of the two psoriasis subgroups, suggested that this group may be more amenable to therapies targeting these pathways. Our study highlights the potential biological relevance of using ensemble decision tree predictors to determine molecular disease subtypes, in what may initially appear to be a homogenous clinical group. The R code used in this paper is available upon request.

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

BackgroundCalmodulin 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.ResultsIn 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.ConclusionOur 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.

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).

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.

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.

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 β-Propellers

β-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 β-propellers coupled with subsequent structural classification and clustering of these models to define subfamilies of β-propellers with common structural, and probable, functional characteristics. We show that it is possible to assign seven-bladed WD β-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 β-propellers and demonstrate that seven-bladed β-propellers composed of WD-repeats are structurally distinct from other seven-bladed β-propellers. We further provide some insights into the multifunctional diversity of the seven-bladed WD β-propeller surfaces. This report once again reinforces the importance of structural data and the usefulness of homology models in functional classification.

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

BackgroundCiliary 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.ResultsThe 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.ConclusionsComputational 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.