by Giovanni Settanni, David Serquera, Piotr E. Marszalek, Emanuele Paci, Laura S. Itzhaki

Ankyrin repeat proteins are elastic materials that unfold and refold sequentially, repeat by repeat, under force. Herein we use atomistic molecular dynamics to compare the mechanical properties of the 7-ankyrin-repeat oncoprotein Gankyrin in isolation and in complex with its binding partner S6-C. We show that the bound S6-C greatly increases the resistance of Gankyrin to mechanical stress. The effect is specific to those repeats of Gankyrin directly in contact with S6-C, and the mechanical ‘hot spots’ of the interaction map to the same repeats as the thermodynamic hot spots. A consequence of stepwise nature of unfolding and the localized nature of ligand binding is that it impacts on all aspects of the protein's mechanical behavior, including the order of repeat unfolding, the diversity of unfolding pathways accessed, the nature of partially unfolded intermediates, the forces required and the work transferred to the system to unfold the whole protein and its parts. Stepwise unfolding thus provides the means to buffer repeat proteins and their binding partners from mechanical stress in the cell. Our results illustrate how ligand binding can control the mechanical response of proteins. The data also point to a cellular mechano-switching mechanism whereby binding between two partner macromolecules is regulated by mechanical stress.

by Sandra V. Bennun, Kevin J. Yarema, Michael J. Betenbaugh, Frederick J. Krambeck

Abnormalities in glycan biosynthesis have been conclusively linked to many diseases but the complexity of glycosylation has hindered the analysis of glycan data in order to identify glycoforms contributing to disease. To overcome this limitation, we developed a quantitative N-glycosylation model that interprets and integrates mass spectral and transcriptomic data by incorporating key glycosylation enzyme activities. Using the cancer progression model of androgen-dependent to androgen-independent Lymph Node Carcinoma of the Prostate (LNCaP) cells, the N-glycosylation model identified and quantified glycan structural details not typically derived from single-stage mass spectral or gene expression data. Differences between the cell types uncovered include increases in H(II) and Ley epitopes, corresponding to greater activity of α2-Fuc-transferase (FUT1) in the androgen-independent cells. The model further elucidated limitations in the two analytical platforms including a defect in the microarray for detecting the GnTV (MGAT5) enzyme. Our results demonstrate the potential of systems glycobiology tools for elucidating key glycan biomarkers and potential therapeutic targets. The integration of multiple data sets represents an important application of systems biology for understanding complex cellular processes.

by Mikael Sunnåker, Alberto Giovanni Busetto, Elina Numminen, Jukka Corander, Matthieu Foll, Christophe Dessimoz

Approximate Bayesian computation (ABC) constitutes a class of computational methods rooted in Bayesian statistics. In all model-based statistical inference, the likelihood function is of central importance, since it expresses the probability of the observed data under a particular statistical model, and thus quantifies the support data lend to particular values of parameters and to choices among different models. For simple models, an analytical formula for the likelihood function can typically be derived. However, for more complex models, an analytical formula might be elusive or the likelihood function might be computationally very costly to evaluate. ABC methods bypass the evaluation of the likelihood function. In this way, ABC methods widen the realm of models for which statistical inference can be considered. ABC methods are mathematically well-founded, but they inevitably make assumptions and approximations whose impact needs to be carefully assessed. Furthermore, the wider application domain of ABC exacerbates the challenges of parameter estimation and model selection. ABC has rapidly gained popularity over the last years and in particular for the analysis of complex problems arising in biological sciences (e.g., in population genetics, ecology, epidemiology, and systems biology).

by Stefka Tyanova, Jürgen Cox, Jesper Olsen, Matthias Mann, Dmitrij Frishman

Phosphorylation at specific residues can activate a protein, lead to its localization to particular compartments, be a trigger for protein degradation and fulfill many other biological functions. Protein phosphorylation is increasingly being studied at a large scale and in a quantitative manner that includes a temporal dimension. By contrast, structural properties of identified phosphorylation sites have so far been investigated in a static, non-quantitative way. Here we combine for the first time dynamic properties of the phosphoproteome with protein structural features. At six time points of the cell division cycle we investigate how the variation of the amount of phosphorylation correlates with the protein structure in the vicinity of the modified site. We find two distinct phosphorylation site groups: intrinsically disordered regions tend to contain sites with dynamically varying levels, whereas regions with predominantly regular secondary structures retain more constant phosphorylation levels. The two groups show preferences for different amino acids in their kinase recognition motifs - proline and other disorder-associated residues are enriched in the former group and charged residues in the latter. Furthermore, these preferences scale with the degree of disorderedness, from regular to irregular and to disordered structures. Our results suggest that the structural organization of the region in which a phosphorylation site resides may serve as an additional control mechanism. They also imply that phosphorylation sites are associated with different time scales that serve different functional needs.

by Cibele Vieira Falkenberg, Leslie M. Loew

The Rho family of GTPases control actin organization during diverse cellular responses (migration, cytokinesis and endocytosis). Although the primary members of this family (RhoA, Rac and Cdc42) have different downstream effects on actin remodeling, the basic mechanism involves targeting to the plasma membrane and activation by GTP binding. Our hypothesis is that the details of GTPase cycling between membrane and cytosol are key to the differential upstream regulation of these biochemical switches. Accordingly, we developed a modeling framework to analyze experimental data for these systems. This analysis can reveal details of GDI-mediated cycling and help distinguish between GDI-dependent and -independent mechanisms, including vesicle trafficking and direct association-dissociation of GTPase with membrane molecules. Analysis of experimental data for Rac membrane cycling reveals that the lower apparent affinity of GDI for RacGTP compared to RacGDP can be fully explained by the faster dissociation of the latter from the membrane. Non-dimensional steady-state solutions for membrane fraction of GTPase are presented in multidimensional charts. This methodology is then used to analyze glucose stimulated Rac cycling in pancreatic β-cells. The charts are used to illustrate the effects of GEFs/GAPs and regulated affinities between GTPases and membrane and/or GDI on the amount of membrane bound GTPase. In a similar fashion, the charts can be used as a guide in assessing how targeted modifications may compensate for altered GTPase-GDI balance in disease scenarios.

by Chiara Pastrello, David Otasek, Kristen Fortney, Giuseppe Agapito, Mario Cannataro, Elize Shirdel, Igor Jurisica

High-throughput technologies produce massive amounts of data. However, individual methods yield data specific to the technique used and biological setup. The integration of such diverse data is necessary for the qualitative analysis of information relevant to hypotheses or discoveries. It is often useful to integrate these datasets using pathways and protein interaction networks to get a broader view of the experiment. The resulting network needs to be able to focus on either the large-scale picture or on the more detailed small-scale subsets, depending on the research question and goals. In this tutorial, we illustrate a workflow useful to integrate, analyze, and visualize data from different sources, and highlight important features of tools to support such analyses.

by Nol Chindapol, Jaap A. Kaandorp, Carolina Cronemberger, Tali Mass, Amatzia Genin

The growth of scleractinian corals is strongly influenced by the effect of water motion. Corals are known to have a high level of phenotypic variation and exhibit a diverse range of growth forms, which often contain a high level of geometric complexity. Due to their complex shape, simulation models represent an important option to complement experimental studies of growth and flow. In this work, we analyzed the impact of flow on coral's morphology by an accretive growth model coupled with advection-diffusion equations. We performed simulations under no-flow and uni-directional flow setup with the Reynolds number constant. The relevant importance of diffusion to advection was investigated by varying the diffusion coefficient, rather than the flow speed in P??clet number. The flow and transport equations were coupled and solved using COMSOL Multiphysics. We then compared the simulated morphologies with a series of Computed Tomography (CT) scans of scleractinian corals Pocillopora verrucosa exposed to various flow conditions in the in situ controlled flume setup. As a result, we found a similar trend associated with the increasing P??clet for both simulated forms and in situ corals; that is uni-directional current tends to facilitate asymmetrical growth response resulting in colonies with branches predominantly developed in the upstream direction. A closer look at the morphological traits yielded an interesting property about colony symmetry and plasticity induced by uni-directional flow. Both simulated and in situ corals exhibit a tendency where the degree of symmetry decreases and compactification increases in conjunction with the augmented P??clet thus indicates the significant importance of hydrodynamics.

by Anatoly Rinberg, Adam L. Taylor, Eve Marder

The crab Cancer borealis undergoes large daily fluctuations in environmental temperature (8–24°C) and must maintain appropriate neural function in the face of this perturbation. In the pyloric circuit of the crab stomatogastric ganglion, we pharmacologically isolated the pacemaker kernel (the AB and PD neurons) and characterized its behavior in response to temperature ramps from 7°C to 31°C. For moderate temperatures, the pacemaker displayed a frequency-temperature curve statistically indistinguishable from that of the intact circuit, and like the intact circuit maintained a constant duty cycle. At high temperatures (above 23°C), a variety of different behaviors were seen: in some preparations the pacemaker increased in frequency, in some it slowed, and in many preparations the pacemaker stopped oscillating (“crashed”). Furthermore, these crashes seemed to fall into two qualitatively different classes. Additionally, the animal-to-animal variability in frequency increased at high temperatures. We used a series of Morris-Lecar mathematical models to gain insight into these phenomena. The biophysical components of the final model have temperature sensitivities similar to those found in nature, and can crash via two qualitatively different mechanisms that resemble those observed experimentally. The crash type is determined by the precise parameters of the model at the reference temperature, 11°C, which could explain why some preparations seem to crash in one way and some in another. Furthermore, even models with very similar behavior at the reference temperature diverge greatly at high temperatures, resembling the experimental observations.

by Satomi Matsuoka, Tatsuo Shibata, Masahiro Ueda

The molecular mechanisms that underlie asymmetric PTEN distribution at the posterior of polarized motile cells and regulate anterior pseudopod formation were addressed by novel single-molecule tracking analysis. Heterogeneity in the lateral mobility of PTEN on a membrane indicated the existence of three membrane-binding states with different diffusion coefficients and membrane-binding lifetimes. The stochastic state transition kinetics of PTEN among these three states were suggested to be regulated spatially along the cell polarity such that only the stable binding state is selectively suppressed at the anterior membrane to cause local PTEN depletion. By incorporating experimentally observed kinetic parameters into a simple mathematical model, the asymmetric PTEN distribution can be explained quantitatively to illustrate the regulatory mechanisms for cellular asymmetry based on an essential causal link between individual stochastic reactions and stable localizations of the ensemble.

by Omry Koren, Dan Knights, Antonio Gonzalez, Levi Waldron, Nicola Segata, Rob Knight, Curtis Huttenhower, Ruth E. Ley

Recent analyses of human-associated bacterial diversity have categorized individuals into ‘enterotypes’ or clusters based on the abundances of key bacterial genera in the gut microbiota. There is a lack of consensus, however, on the analytical basis for enterotypes and on the interpretation of these results. We tested how the following factors influenced the detection of enterotypes: clustering methodology, distance metrics, OTU-picking approaches, sequencing depth, data type (whole genome shotgun (WGS) vs.16S rRNA gene sequence data), and 16S rRNA region. We included 16S rRNA gene sequences from the Human Microbiome Project (HMP) and from 16 additional studies and WGS sequences from the HMP and MetaHIT. In most body sites, we observed smooth abundance gradients of key genera without discrete clustering of samples. Some body habitats displayed bimodal (e.g., gut) or multimodal (e.g., vagina) distributions of sample abundances, but not all clustering methods and workflows accurately highlight such clusters. Because identifying enterotypes in datasets depends not only on the structure of the data but is also sensitive to the methods applied to identifying clustering strength, we recommend that multiple approaches be used and compared when testing for enterotypes.