Vesicular and non-vesicular transport feed distinct glycosylation pathways in the Golgi.

Newly synthesized proteins and lipids are transported across the Golgi complex via different mechanisms whose respective roles are not completely clear. We previously identified a non-vesicular intra-Golgi transport pathway for glucosylceramide (GlcCer)--the common precursor of the different series of glycosphingolipids-that is operated by the cytosolic GlcCer-transfer protein FAPP2 (also known as PLEKHA8) (ref. 1). However, the molecular determinants of the FAPP2-mediated transfer of GlcCer from the cis-Golgi to the trans-Golgi network, as well as the physiological relevance of maintaining two parallel transport pathways of GlcCer--vesicular and non-vesicular--through the Golgi, remain poorly defined. Here, using mouse and cell models, we clarify the molecular mechanisms underlying the intra-Golgi vectorial transfer of GlcCer by FAPP2 and show that GlcCer is channelled by vesicular and non-vesicular transport to two topologically distinct glycosylation tracks in the Golgi cisternae and the trans-Golgi network, respectively. Our results indicate that the transport modality across the Golgi complex is a key determinant for the glycosylation pattern of a cargo and establish a new paradigm for the branching of the glycosphingolipid synthetic pathway.

Threshold-controlled ubiquitination of the EGFR directs receptor fate.

How the cell converts graded signals into threshold-activated responses is a question of great biological relevance. Here, we uncover a nonlinear modality of epidermal growth factor receptor (EGFR)-activated signal transduction, by demonstrating that the ubiquitination of the EGFR at the PM is threshold controlled. The ubiquitination threshold is mechanistically determined by the cooperative recruitment of the E3 ligase Cbl, in complex with Grb2, to the EGFR. This, in turn, is dependent on the simultaneous presence of two phosphotyrosines, pY1045 and either one of pY1068 or pY1086, on the same EGFR moiety. The dose-response curve of EGFR ubiquitination correlate precisely with the non-clathrin endocytosis (NCE) mode of EGFR internalization. Finally, EGFR-NCE mechanistically depends on EGFR ubiquitination, as the two events can be simultaneously re-engineered on a phosphorylation/ubiquitination-incompetent EGFR backbone. Since NCE controls the degradation of the EGFR, our findings have implications for how the cell responds to increasing levels of EGFR signalling, by varying the balance of receptor signalling and degradation/attenuation.

Growth against entropy in bacterial metabolism: the phenotypic trade-off behind empirical growth rate distributions in E. coli.

The solution space of genome-scale models of cellular metabolism provides a map between physically viable flux configurations and cellular metabolic phenotypes described, at the most basic level, by the corresponding growth rates. By sampling the solution space of E. coli's metabolic network, we show that empirical growth rate distributions recently obtained in experiments at single-cell resolution can be explained in terms of a trade-off between the higher fitness of fast-growing phenotypes and the higher entropy of slow-growing ones. Based on this, we propose a minimal model for the evolution of a large bacterial population that captures this trade-off. The scaling relationships observed in experiments encode, in such frameworks, for the same distance from the maximum achievable growth rate, the same degree of growth rate maximization, and/or the same rate of phenotypic change. Being grounded on genome-scale metabolic network reconstructions, these results allow for multiple implications and extensions in spite of the underlying conceptual simplicity.

The Caenorhabditis elegans Elongator complex regulates neuronal alpha-tubulin acetylation.

Although acetylated alpha-tubulin is known to be a marker of stable microtubules in neurons, precise factors that regulate alpha-tubulin acetylation are, to date, largely unknown. Therefore, a genetic screen was employed in the nematode Caenorhabditis elegans that identified the Elongator complex as a possible regulator of alpha-tubulin acetylation. Detailed characterization of mutant animals revealed that the acetyltransferase activity of the Elongator is indeed required for correct acetylation of microtubules and for neuronal development. Moreover, the velocity of vesicles on microtubules was affected by mutations in Elongator. Elongator mutants also displayed defects in neurotransmitter levels. Furthermore, acetylation of alpha-tubulin was shown to act as a novel signal for the fine-tuning of microtubules dynamics by modulating alpha-tubulin turnover, which in turn affected neuronal shape. Given that mutations in the acetyltransferase subunit of the Elongator (Elp3) and in a scaffold subunit (Elp1) have previously been linked to human neurodegenerative diseases, namely Amyotrophic Lateral Sclerosis and Familial Dysautonomia respectively highlights the importance of this work and offers new insights to understand their etiology.

Low duplicability and network fragility of cancer genes.

We identified genomic and network properties of approximately 600 genes mutated in different cancer types. These genes tend not to duplicate but, unlike most human singletons, they encode central hubs of highly interconnected modules within the protein-protein interaction network (PIN). We find that cancer genes are fragile components of the human gene repertoire, sensitive to dosage modification. Furthermore, other nodes of the human PIN with similar properties are rare and probably enriched in candidate cancer genes.

Modeling networks of coupled enzymatic reactions using the total quasi-steady state approximation.

In metabolic networks, metabolites are usually present in great excess over the enzymes that catalyze their interconversion, and describing the rates of these reactions by using the Michaelis-Menten rate law is perfectly valid. This rate law assumes that the concentration of enzyme-substrate complex (C) is much less than the free substrate concentration (S0). However, in protein interaction networks, the enzymes and substrates are all proteins in comparable concentrations, and neglecting C with respect to S0 is not valid. Borghans, DeBoer, and Segel developed an alternative description of enzyme kinetics that is valid when C is comparable to S0. We extend this description, which Borghans et al. call the total quasi-steady state approximation, to networks of coupled enzymatic reactions. First, we analyze an isolated Goldbeter-Koshland switch when enzymes and substrates are present in comparable concentrations. Then, on the basis of a real example of the molecular network governing cell cycle progression, we couple two and three Goldbeter-Koshland switches together to study the effects of feedback in networks of protein kinases and phosphatases. Our analysis shows that the total quasi-steady state approximation provides an excellent kinetic formalism for protein interaction networks, because (1) it unveils the modular structure of the enzymatic reactions, (2) it suggests a simple algorithm to formulate correct kinetic equations, and (3) contrary to classical Michaelis-Menten kinetics, it succeeds in faithfully reproducing the dynamics of the network both qualitatively and quantitatively.

Quantifying the entropic cost of cellular growth control.

Viewing the ways a living cell can organize its metabolism as the phase space of a physical system, regulation can be seen as the ability to reduce the entropy of that space by selecting specific cellular configurations that are, in some sense, optimal. Here we quantify the amount of regulation required to control a cell's growth rate by a maximum-entropy approach to the space of underlying metabolic phenotypes, where a configuration corresponds to a metabolic flux pattern as described by genome-scale models. We link the mean growth rate achieved by a population of cells to the minimal amount of metabolic regulation needed to achieve it through a phase diagram that highlights how growth suppression can be as costly (in regulatory terms) as growth enhancement. Moreover, we provide an interpretation of the inverse temperature ß controlling maximum-entropy distributions based on the underlying growth dynamics. Specifically, we show that the asymptotic value of ß for a cell population can be expected to depend on (i) the carrying capacity of the environment, (ii) the initial size of the colony, and (iii) the probability distribution from which the inoculum was sampled. Results obtained for E. coli and human cells are found to be remarkably consistent with empirical evidence.

Quantitative constraint-based computational model of tumor-to-stroma coupling via lactate shuttle.

Cancer cells utilize large amounts of ATP to sustain growth, relying primarily on non-oxidative, fermentative pathways for its production. In many types of cancers this leads, even in the presence of oxygen, to the secretion of carbon equivalents (usually in the form of lactate) in the cell's surroundings, a feature known as the Warburg effect. While the molecular basis of this phenomenon are still to be elucidated, it is clear that the spilling of energy resources contributes to creating a peculiar microenvironment for tumors, possibly characterized by a degree of toxicity. This suggests that mechanisms for recycling the fermentation products (e.g. a lactate shuttle) may be active, effectively inducing a mutually beneficial metabolic coupling between aberrant and non-aberrant cells. Here we analyze this scenario through a large-scale in silico metabolic model of interacting human cells. By going beyond the cell-autonomous description, we show that elementary physico-chemical constraints indeed favor the establishment of such a coupling under very broad conditions. The characterization we obtained by tuning the aberrant cell's demand for ATP, amino-acids and fatty acids and/or the imbalance in nutrient partitioning provides quantitative support to the idea that synergistic multi-cell effects play a central role in cancer sustainment.

Quantitative analysis reveals how EGFR activation and downregulation are coupled in normal but not in cancer cells.

Ubiquitination of the epidermal growth factor receptor (EGFR) that occurs when Cbl and Grb2 bind to three phosphotyrosine residues (pY1045, pY1068 and pY1086) on the receptor displays a sharp threshold effect as a function of EGF concentration. Here we use a simple modelling approach together with experiments to show that the establishment of the threshold requires both the multiplicity of binding sites and cooperative binding of Cbl and Grb2 to the EGFR. While the threshold is remarkably robust, a more sophisticated model predicted that it could be modulated as a function of EGFR levels on the cell surface. We confirmed experimentally that the system has evolved to perform optimally at physiological levels of EGFR. As a consequence, this system displays an intrinsic weakness that causes--at the supraphysiological levels of receptor and/or ligand associated with cancer--uncoupling of the mechanisms leading to signalling through phosphorylation and attenuation through ubiquitination.

Unravelling druggable signalling networks that control F508del-CFTR proteostasis.

Cystic fibrosis (CF) is caused by mutations in CF transmembrane conductance regulator (CFTR). The most frequent mutation (F508del-CFTR) results in altered proteostasis, that is, in the misfolding and intracellular degradation of the protein. The F508del-CFTR proteostasis machinery and its homeostatic regulation are well studied, while the question whether 'classical' signalling pathways and phosphorylation cascades might control proteostasis remains barely explored. Here, we have unravelled signalling cascades acting selectively on the F508del-CFTR folding-trafficking defects by analysing the mechanisms of action of F508del-CFTR proteostasis regulator drugs through an approach based on transcriptional profiling followed by deconvolution of their gene signatures. Targeting multiple components of these signalling pathways resulted in potent and specific correction of F508del-CFTR proteostasis and in synergy with pharmacochaperones. These results provide new insights into the physiology of cellular proteostasis and a rational basis for developing effective pharmacological correctors of the F508del-CFTR defect.

Velocity fluctuations and dispersion in a simple porous medium.

We model a fluid-filled disordered porous medium by a lattice-Boltzmann system with randomly broken links. The broken links exert a friction on the fluid without excluding volume. Such a model closely mimics the idealized picture of a porous medium, which is often used in the theoretical analysis of hydrodynamic dispersion. We find that the Brinkman equation describes both the mean flow characteristics and the spatial decay of velocity fluctuations in the system. However, the temporal decay of the velocity correlations (that a particle experiences as it moves with the fluid), cannot be simply related to the spatial decay. It is this temporal decay that determines the dispersivity. Thus, hydrodynamic dispersion is generally greater than theories based on spatial correlations would imply. This is particularly true at high densities, where such theories considerably underestimate both the magnitude and transient time scale for dispersion. Nonetheless, temporal velocity correlations are still ultimately screened and the hydrodynamic dispersion coefficient converges exponentially. The long-lived transients reported for more realistic systems must therefore be due explicitly to the presence of excluded volume.

Counting and correcting thermodynamically infeasible flux cycles in genome-scale metabolic networks.

Thermodynamics constrains the flow of matter in a reaction network to occur through routes along which the Gibbs energy decreases, implying that viable steady-state flux patterns should be void of closed reaction cycles. Identifying and removing cycles in large reaction networks can unfortunately be a highly challenging task from a computational viewpoint. We propose here a method that accomplishes it by combining a relaxation algorithm and a Monte Carlo procedure to detect loops, with ad hoc rules (discussed in detail) to eliminate them. As test cases, we tackle (a) the problem of identifying infeasible cycles in the E. coli metabolic network and (b) the problem of correcting thermodynamic infeasibilities in the Flux-Balance-Analysis solutions for 15 human cell-type-specific metabolic networks. Results for (a) are compared with previous analyses of the same issue, while results for (b) are weighed against alternative methods to retrieve thermodynamically viable flux patterns based on minimizing specific global quantities. Our method, on the one hand, outperforms previous techniques and, on the other, corrects loopy solutions to Flux Balance Analysis. As a byproduct, it also turns out to be able to reveal possible inconsistencies in model reconstructions.

Lattice-Boltzmann Simulations of Ionic Current Modulation by DNA Translocation.

We present a numerical study of the effect of DNA translocation on the ionic current through a nanopore. We use a coarse-grained model to solve the electrokinetic equations at the Poisson-Boltzmann level for the microions, coupled to a lattice-Boltzmann equation for the solvent hydrodynamics. In most cases, translocation leads to a reduction in the ionic current. However, at low salt concentrations (large screening lengths) we find ionic current enhancement due to translocation. In an unstructured pore, translocation of the helical charge distribution of the DNA has no effect on the ionic current. However, if a localized charge probe is placed on the wall of the nanopore, we observe ionic current modulations that, though weak, should be experimentally observable.

Lattice-Boltzmann simulation of the sedimentation of charged disks.

We report a series of lattice-Boltzmann simulations of the sedimentation velocity of charged disks. In these simulations, we explicitly account for the hydrodynamic and electrostatic forces on disks and on their electrical double layer. By comparing our results with those for spheres with equal surface and charge, we can clarify the effect of the particle shape on the sedimentation process. We find that disks and spheres exhibit a different dependence of the sedimentation velocity on the Debye screening length. An analysis of the behavior of highly charged disks (beyond the scope of the linearized Poisson-Boltzmann equation) shows that, in that regime, the charge dependence of the sedimentation velocity of disks and spheres is similar. This suggests that, at high charge, the effective hydrodynamic shape of the disks becomes more spherical.

Discrete solution of the electrokinetic equations.

We present a robust scheme for solving the electrokinetic equations. This goal is achieved by combining the lattice-Boltzmann method with a discrete solution of the convection-diffusion equation for the different charged and neutral species that compose the fluid. The method is based on identifying the elementary fluxes between nodes, which ensures the absence of spurious fluxes in equilibrium. We show how the model is suitable to study electro-osmotic flows. As an illustration, we show that, by introducing appropriate dynamic rules in the presence of solid interfaces, we can compute the sedimentation velocity (and hence the sedimentation potential) of a charged sphere. Our approach does not assume linearization of the Poisson-Boltzmann equation and allows us for a wide variation of the Peclet number.