Linear Control Theory for Gene Network Modeling

Systems biology is an interdisciplinary field that aims at understanding complex interactions in cells. Here we demonstrate that linear control theory can provide valuable insight and practical tools for the characterization of complex biological networks. We provide the foundation for such analyses through the study of several case studies including cascade and parallel forms, feedback and feedforward loops. We reproduce experimental results and provide rational analysis of the observed behavior. We demonstrate that methods such as the transfer function (frequency domain) and linear state-space (time domain) can be used to predict reliably the properties and transient behavior of complex network topologies and point to specific design strategies for synthetic networks.     

Towards a Rigorous Assessment of Systems Biology Models: The DREAM3 Challenges

Background

Systems biology has embraced computational modeling in response to the quantitative nature and increasing scale of contemporary data sets. The onslaught of data is accelerating as molecular profiling technology evolves. The Dialogue for Reverse Engineering Assessments and Methods (DREAM) is a community effort to catalyze discussion about the design, application, and assessment of systems biology models through annual reverse-engineering challenges.

Methodology and Principal Findings

We describe our assessments of the four challenges associated with the third DREAM conference which came to be known as the DREAM3 challenges: signaling cascade identification, signaling response prediction, gene expression prediction, and the DREAM3 in silico network challenge. The challenges, based on anonymized data sets, tested participants in network inference and prediction of measurements. Forty teams submitted 413 predicted networks and measurement test sets. Overall, a handful of best-performer teams were identified, while a majority of teams made predictions that were equivalent to random. Counterintuitively, combining the predictions of multiple teams (including the weaker teams) can in some cases improve predictive power beyond that of any single method.

Conclusions

DREAM provides valuable feedback to practitioners of systems biology modeling. Lessons learned from the predictions of the community provide much-needed context for interpreting claims of efficacy of algorithms described in the scientific literature.     

Role of the p38 mitogen-activated protein kinase pathway in the generation of the effects of imatinib mesylate (STI571) in BCR-ABL-expressing cells

Imatinib mesylate (STI571), a specific inhibitor of the BCR-ABL tyrosine kinase, exhibits potent antileukemic effects in vitro and in vivo. Despite the well established role of STI571 in the treatment of chronic myelogenous leukemia, the precise mechanisms by which inhibition of BCR-ABL tyrosine kinase activity results in generation of antileukemic responses remain unknown. In the present study we provide evidence that treatment of CML-derived BCR-ABL-expressing leukemia cells with STI571 results in activation of the p38 mitogen-activated protein (MAP) kinase signaling pathway. Our data indicate that STI571 induces phosphorylation of the p38 and activation of its kinase domain, in KT-1 cells and other BCR-ABL-expressing cell lines. We also identify the kinases MAP kinase-activated protein kinase-2 and Msk1 as two downstream effectors of p38, activated during inhibition of BCR-ABL activity by STI571. Importantly, pharmacological inhibition of p38 reverses the growth inhibitory effects of STI571 on primary leukemic colony-forming unit granulocyte/macrophage progenitors from patients with CML. Altogether, our data establish that activation of the p38 MAP kinase signaling cascade plays an important role in the generation of the effects of STI571 on BCR-ABL-expressing cells. They also suggest that, in addition to activation of mitogenic pathways, BCR-ABL promotes leukemogenesis by suppressing the function of growth inhibitory signaling cascades.     

Binding of human immunodeficiency virus type 1 to immature dendritic cells can occur independently of DC-SIGN and mannose binding C-type lectin receptors via a cholesterol-dependent pathway

Interactions of human immunodeficiency virus type 1 (HIV-1) with immature dendritic cells (DC) are believed to be multifactorial and involve binding to the CD4 antigen, DC-specific ICAM-3-grabbing nonintegrin (DC-SIGN), mannose binding C-type lectin receptors (MCLR), and heparan sulfate proteoglycans (HSPG). In this study we assessed the relative contributions of these previously defined virus attachment factors to HIV binding and accumulation in DC and the subsequent transfer of the bound virus particle to CD4(+) T cells. Using competitive inhibitors of HIV-1 attachment to DC, we have identified the existence of DC-SIGN-, MCLR-, and HSPG-independent mechanism(s) of HIV attachment and internalization. Furthermore, virus particles bound by DC independently of CD4, DC-SIGN, MCLR, and HSPG are efficiently transmitted to T cells. Treatment of virus particles with the protease subtilisin or treatment of immature DC with trypsin significantly reduced virus binding, thus demonstrating the role of HIV envelope glycoprotein interactions with unidentified DC-surface factor(s). Finally, this DC-mediated virus binding and internalization are dependent on lipid rafts. We propose that pathways to HIV-1 attachment and uptake in DC exhibit functional redundancy; that is, they are made up of multiple independent activities that can, at least in part, compensate for one another.     

Alterations in the mechanical properties of the human chondrocyte pericellular matrix with osteoarthritis

In articular cartilage, chondrocytes are surrounded by a pericellular matrix (PCM), which together with the chondrocyte have been termed the "chondron." While the precise function of the PCM is not know there has been considerable speculation that it plays a role in regulating the biomechanical environment of the chondrocyte. In this study, we measured the Young's modulus of the PCM from normal and osteoarthritic cartilage using the micropipette aspiration technique, coupled with a newly developed axisymmetric elastic layered half-space model of the experimental configuration. Viable, intact chondrons were extracted from human articular cartilage using a new microaspiration-based isolation technique. In normal cartilage, the Young's modulus of the PCM was similar in chondrons isolated from the surface zone (68.9 +/- 18.9 kPa) as compared to the middle and deep layers (62.0 +/- 30.5 kPa). However, the mean Young's modulus of the PCM (pooled for the two zones) was significantly decreased in osteoarthritic cartilage (66.5 +/- 23.3 kPa versus 41.3 +/- 21.1 kPa, p < 0.001). In combination with previous theoretical models of cell-matrix interactions in cartilage, these findings suggest that the PCM has an important influence on the stress-strain environment of the chondrocyte that potentially varies with depth from the cartilage surface. Furthermore, the significant loss of PCM stiffness that was observed in osteoarthritic cartilage may affect the magnitude and distribution of biomechanical signals perceived by the chondrocytes.     

The binding of beta- and gamma-cyclodextrins to glycogen phosphorylase b: kinetic and crystallographic studies

A number of regulatory binding sites of glycogen phosphorylase (GP), such as the catalytic, the inhibitor, and the new allosteric sites are currently under investigation as targets for inhibition of hepatic glycogenolysis under high glucose concentrations; in some cases specific inhibitors are under evaluation in human clinical trials for therapeutic intervention in type 2 diabetes. In an attempt to investigate whether the storage site can be exploited as target for modulating hepatic glucose production, alpha-, beta-, and gamma-cyclodextrins were identified as moderate mixed-type competitive inhibitors of GPb (with respect to glycogen) with K(i) values of 47.1, 14.1, and 7.4 mM, respectively. To elucidate the structural basis of inhibition, we determined the structure of GPb complexed with beta- and gamma-cyclodextrins at 1.94 A and 2.3 A resolution, respectively. The structures of the two complexes reveal that the inhibitors can be accommodated in the glycogen storage site of T-state GPb with very little change of the tertiary structure and provide a basis for understanding their potency and subsite specificity. Structural comparisons of the two complexes with GPb in complex with either maltopentaose (G5) or maltoheptaose (G7) show that beta- and gamma-cyclodextrins bind in a mode analogous to the G5 and G7 binding with only some differences imposed by their cyclic conformations. It appears that the binding energy for stabilization of enzyme complexes derives from hydrogen bonding and van der Waals contacts to protein residues. The binding of alpha-cyclodextrin and octakis (2,3,6-tri-O-methyl)-gamma-cyclodextrin was also investigated, but none of them was bound in the crystal; moreover, the latter did not inhibit the phosphorylase reaction.     

Crystal structure of rabbit muscle glycogen phosphorylase a in complex with a potential hypoglycaemic drug at 2.0 A resolution

CP320626 has been identified as a potent inhibitor, synergistic with glucose, of human liver glycogen phosphorylase a (LGPa), a possible target for type 2 diabetes therapy. CP320626 is also a potent inhibitor of human muscle GPa. In order to elucidate the structural basis of the mechanism of CP320626 inhibition, the structures of T state rabbit muscle GPa (MGPa) in complex with glucose and in complex with both glucose and CP320626 were determined at 2.0 A resolution, and refined to crystallographic R values of 0.179 (R(free)=0.218) and 0.207 (R(free)=0.235), respectively. CP320626 binds at the new allosteric site, some 33 A from the catalytic site, where glucose binds. The binding of CP320626 to MGPa does not promote extensive conformational changes except for small shifts of the side chain atoms of residues R60, V64, and K191. Both CP320626 and glucose promote the less active T state, while structural comparisons of MGPa-glucose-CP320626 complex with LGPa complexed with a related compound (CP403700) and a glucose analogue inhibitor indicate that the residues of the new allosteric site, conserved in the two isozymes, show no significant differences in their positions.