Parameter estimation in systems biology models using spline approximation

Background  

Mathematical models for revealing the dynamics and interactions properties of biological systems play an important role in computational systems biology. The inference of model parameter values from time-course data can be considered as a "reverse engineering" process and is still one of the most challenging tasks. Many parameter estimation methods have been developed but none of these methods is effective for all cases and can overwhelm all other approaches. Instead, various methods have their advantages and disadvantages. It is worth to develop parameter estimation methods which are robust against noise, efficient in computation and flexible enough to meet different constraints.     

The KCNJ11 E23K polymorphism and progression of glycaemia in Southern Chinese: a long-term prospective study

Context

The KCNJ11 E23K variant is associated with type 2 diabetes mellitus (T2DM) in cross-sectional studies, but conflicting findings have been reported from prospective studies.

Objective

This study aimed to evaluate whether the E23K variant could predict glycaemic progression in a Southern Chinese population.

Methods/Principal Findings

We performed a long-term prospective study on 1912 subjects from the Hong Kong Cardiovascular Risk Factors Prevalence Study (CRISPS). The KCNJ11 E23K variant was associated with the progression to prediabetes after a median interval of 12 years on multinomial logistic regression analysis, even after adjustment for traditional risk factors (OR 1.29, P_{age, sex, BMI and fasting plasma glucose [FPG] adjusted}  = 0.02). Based on Cox proportional hazard regression analysis, the E23K variant also predicted incident prediabetes (HR 1.18, P_{age, sex, BMI and FPG adjusted}  = 0.021). However, E23K was not associated with the progression to T2DM in either multinomial or Cox regression analysis, and the association of E23K with glycaemic progression to either prediabetes or T2DM was significant only in unadjusted Cox regression analysis (P = 0.039). In a meta-analysis of eight prospective studies including our own, involving 15680 subjects, the E23K variant was associated with incident T2DM (fixed effect: OR 1.10, P = 4×10⁻³; random effect: OR 1.11, P = 0.035).

Conclusions

Our study has provided supporting evidence for the role of the E23K variant in glycaemic progression in Chinese, with its effect being more evident in the early stage of T2DM, as the subjects progressed from normal glucose tolerance to prediabetes.     

Microenvironmental geometry guides platelet adhesion and spreading: a quantitative analysis at the single cell level

To activate clot formation and maintain hemostasis, platelets adhere and spread onto sites of vascular injury. Although this process is well-characterized biochemically, how the physical and spatial cues in the microenvironment affect platelet adhesion and spreading remain unclear. In this study, we applied deep UV photolithography and protein micro/nanostamping to quantitatively investigate and characterize the spatial guidance of platelet spreading at the single cell level and with nanoscale resolution. Platelets adhered to and spread only onto micropatterned collagen or fibrinogen surfaces and followed the microenvironmental geometry with high fidelity and with single micron precision. Using micropatterned lines of different widths, we determined that platelets are able to conform to micropatterned stripes as thin as 0.6 μm and adopt a maximum aspect ratio of 19 on those protein patterns. Interestingly, platelets were also able to span and spread over non-patterned regions of up to 5 μm, a length consistent with that of maximally extended filopodia. This process appears to be mediated by platelet filopodia that are sensitive to spatial cues. Finally, we observed that microenvironmental geometry directly affects platelet biology, such as the spatial organization and distribution of the platelet actin cytoskeleton. Our data demonstrate that platelet spreading is a finely-tuned and spatially-guided process in which spatial cues directly influence the biological aspects of how clot formation is regulated.     

The antiangiogenic 16K prolactin impairs functional tumor neovascularization by inhibiting vessel maturation

Background

Angiogenesis, the formation of new blood vessels from existing vasculature, plays an essential role in tumor growth, invasion, and metastasis. 16K hPRL, the antiangiogenic 16-kDa N-terminal fragment of human prolactin was shown to prevent tumor growth and metastasis by modifying tumor vessel morphology.

Methodology/Principal Findings

Here we investigated the effect of 16K hPRL on tumor vessel maturation and on the related signaling pathways. We show that 16K hPRL treatment leads, in a murine B16-F10 tumor model, to a dysfunctional tumor vasculature with reduced pericyte coverage, and disruption of the PDGF-B/PDGFR-B, Ang/Tie2, and Delta/Notch pathways. In an aortic ring assay, 16K hPRL impairs endothelial cell and pericyte outgrowth from the vascular ring. In addition, 16K hPRL prevents pericyte migration to endothelial cells. This event was independent of a direct inhibitory effect of 16K hPRL on pericyte viability, proliferation, or migration. In endothelial cell-pericyte cocultures, we found 16K hPRL to disturb Notch signaling.

Conclusions/Significance

Taken together, our data show that 16K hPRL impairs functional tumor neovascularization by inhibiting vessel maturation and for the first time that an endogenous antiangiogenic agent disturbs Notch signaling. These findings provide new insights into the mechanisms of 16K hPRL action and highlight its potential for use in anticancer therapy.     

RET mutational spectrum in Hirschsprung disease: evaluation of 601 Chinese patients

Rare (RVs) and common variants of the RET gene contribute to Hirschsprung disease (HSCR; congenital aganglionosis). While RET common variants are strongly associated with the commonest manifestation of the disease (males; short-segment aganglionosis; sporadic), rare coding sequence (CDS) variants are more frequently found in the lesser common and more severe forms of the disease (females; long/total colonic aganglionosis; familial).Here we present the screening for RVs in the RET CDS and intron/exon boundaries of 601 Chinese HSCR patients, the largest number of patients ever reported. We identified 61 different heterozygous RVs (50 novel) distributed among 100 patients (16.64%). Those include 14 silent, 29 missense, 5 nonsense, 4 frame-shifts, and one in-frame amino-acid deletion in the CDS, two splice-site deletions, 4 nucleotide substitutions and a 22-bp deletion in the intron/exon boundaries and 1 single-nucleotide substitution in the 5' untranslated region. Exonic variants were mainly clustered in RET the extracellular domain. RET RVs were more frequent among patients with the most severe phenotype (24% vs. 15% in short-HSCR). Phasing RVs with the RET HSCR-associated haplotype suggests that RVs do not underlie the undisputable association of RET common variants with HSCR. None of the variants were found in 250 Chinese controls.     

Apoptosis,cytokine and chemokine induction by non-structural 1 (NS1) proteins encoded by different influenza subtypes

Background

Influenza pandemic remains a serious threat to human health. Viruses of avian origin, H5N1, H7N7 and H9N2, have repeatedly crossed the species barrier to infect humans. Recently, a novel strain originated from swine has evolved to a pandemic. This study aims at improving our understanding on the pathogenic mechanism of influenza viruses, in particular the role of non-structural (NS1) protein in inducing pro-inflammatory and apoptotic responses.

Methods

Human lung epithelial cells (NCI-H292) was used as an in-vitro model to study cytokine/chemokine production and apoptosis induced by transfection of NS1 mRNA encoded by seven infleunza subtypes (seasonal and pandemic H1, H2, H3, H5, H7, and H9), respectively.

Results

The results showed that CXCL-10/IP10 was most prominently induced (> 1000 folds) and IL-6 was slightly induced (< 10 folds) by all subtypes. A subtype-dependent pattern was observed for CCL-2/MCP-1, CCL3/MIP-1α, CCL-5/RANTES and CXCL-9/MIG; where induction by H5N1 was much higher than all other subtypes examined. All subtypes induced a similar temporal profile of apoptosis following transfection. The level of apoptosis induced by H5N1 was remarkably higher than all others. The cytokine/chemokine and apoptosis inducing ability of the 2009 pandemic H1N1 was similar to previous seasonal strains.

Conclusions

In conclusion, the NS1 protein encoded by H5N1 carries a remarkably different property as compared to other avian and human subtypes, and is one of the keys to its high pathogenicity. NCI-H292 cells system proves to be a good in-vitro model to delineate the property of NS1 proteins.

     

A simple, versatile and efficient method to genetically modify human monocyte-derived dendritic cells with HIV-1-derived lentiviral vectors

Lentiviral vectors derived from the human immunodeficiency type 1 virus (HIV-1 LV) are among the finest tools available today for the genetic modification of human monocyte-derived dendritic cells (MDDCs). However, this process is largely inefficient because MDDCs show a strong resistance to HIV-1 transduction. Here we describe a step-by-step protocol from the production of LVs to cell transduction that allows the efficient genetic modification of MDDCs. This protocol can be completed in 23 d from the initial phase of LV production to the final analysis of the results of MDDC transduction. The method relies on the simultaneous addition of HIV-1 LVs along with noninfectious virion-like particles carrying Vpx, a nonstructural protein encoded by the simian immunodeficiency virus (Vpx-VLPs). When thus provided in target cells, Vpx exerts a strong positive effect on incoming LVs by counteracting the restriction present in MDDCs; accordingly, 100% of cells can be transduced with low viral inputs. Vpx-VLPs will improve the efficiency of LV-mediated transduction of MDDCs with vectors for both ectopic gene expression and depletion studies.     

A retrospective analysis of the haemodynamic and metabolic effects of fluid resuscitation in Vietnamese adults with severe falciparum malaria

Background

Optimising the fluid resuscitation of patients with severe malaria is a simple and potentially cost-effective intervention. Current WHO guidelines recommend central venous pressure (CVP) guided, crystalloid based, resuscitation in adults.

Methods

Prospectively collected haemodynamic data from intervention trials in Vietnamese adults with severe malaria were analysed retrospectively to assess the responses to fluid resuscitation.

Results

43 patients were studied of whom 24 received a fluid load. The fluid load resulted in an increase in cardiac index (mean increase: 0.75 L/min/m² (95% Confidence interval (CI): 0.41 to 1.1)), but no significant change in acid-base status post resuscitation (mean increase base deficit 0.6 mmol/L (95% CI: −0.1 to 1.3). The CVP and PAoP (pulmonary artery occlusion pressure) were highly inter-correlated (r_s = 0.7, p<0.0001), but neither were correlated with acid-base status (arterial pH, serum bicarbonate, base deficit) or respiratory status (PaO₂/FiO₂ ratio). There was no correlation between the oxygen delivery (DO₂) and base deficit at the 63 time-points where they were assessed simultaneously (r_s = −0.09, p = 0.46).

Conclusions

In adults with severe falciparum malaria there was no observed improvement in patient outcomes or acid-base status with fluid loading. Neither CVP nor PAoP correlated with markers of end-organ perfusion or respiratory status, suggesting these measures are poor predictors of their fluid resuscitation needs.     

A splice isoform of DNedd4, DNedd4-long, negatively regulates neuromuscular synaptogenesis and viability in Drosophila

Background

Neuromuscular (NM) synaptogenesis is a tightly regulated process. We previously showed that in flies, Drosophila Nedd4 (dNedd4/dNedd4S) is required for proper NM synaptogenesis by promoting endocytosis of commissureless from the muscle surface, a pre-requisite step for muscle innervation. DNedd4 is an E3 ubiquitin ligase comprised of a C2-WW(x3)-Hect domain architecture, which includes several splice isoforms, the most prominent ones are dNedd4-short (dNedd4S) and dNedd4-long (dNedd4Lo).

Methodology/Principal Findings

We show here that while dNedd4S is essential for NM synaptogenesis, the dNedd4Lo isoform inhibits this process and causes lethality. Our results reveal that unlike dNedd4S, dNedd4Lo cannot rescue the lethality of dNedd4 null (DNedd4^T121FS) flies. Moreover, overexpression of UAS-dNedd4Lo specifically in wildtype muscles leads to NM synaptogenesis defects, impaired locomotion and larval lethality. These negative effects of dNedd4Lo are ameliorated by deletion of two regions (N-terminus and Middle region) unique to this isoform, and by inactivating the catalytic activity of dNedd4Lo, suggesting that these unique regions, as well as catalytic activity, are responsible for the inhibitory effects of dNedd4Lo on synaptogenesis. In accord with these findings, we demonstrate by sqRT-PCR an increase in dNedd4S expression relative to the expression of dNedd4Lo during embryonic stages when synaptogenesis takes place.

Conclusion/Significance

Our studies demonstrate that splice isoforms of the same dNedd4 gene can lead to opposite effects on NM synaptogenesis.