Ten simple rules for failing successfully in academia.

Failure is an integral part of life and by extension academia. At the same time, failure is often ignored, with potentially negative consequences both for the science and the scientists involved. This article provides several strategies for learning from and dealing with failure instead of ignoring it. Hopefully, our recommendations are widely applicable, while still taking into account individual differences between academics. These simple rules allow academics to further develop their own strategies for failing successfully in academia.

A multi-state model of the CaMKII dodecamer suggests a role for calmodulin in maintenance of autophosphorylation.

Ca2+/calmodulin-dependent protein kinase II (CaMKII) accounts for up to 2 percent of all brain protein and is essential to memory function. CaMKII activity is known to regulate dynamic shifts in the size and signaling strength of neuronal connections, a process known as synaptic plasticity. Increasingly, computational models are used to explore synaptic plasticity and the mechanisms regulating CaMKII activity. Conventional modeling approaches may exclude biophysical detail due to the impractical number of state combinations that arise when explicitly monitoring the conformational changes, ligand binding, and phosphorylation events that occur on each of the CaMKII holoenzyme's subunits. To manage the combinatorial explosion without necessitating bias or loss in biological accuracy, we use a specialized syntax in the software MCell to create a rule-based model of a twelve-subunit CaMKII holoenzyme. Here we validate the rule-based model against previous experimental measures of CaMKII activity and investigate molecular mechanisms of CaMKII regulation. Specifically, we explore how Ca2+/CaM-binding may both stabilize CaMKII subunit activation and regulate maintenance of CaMKII autophosphorylation. Noting that Ca2+/CaM and protein phosphatases bind CaMKII at nearby or overlapping sites, we compare model scenarios in which Ca2+/CaM and protein phosphatase do or do not structurally exclude each other's binding to CaMKII. Our results suggest a functional mechanism for the so-called "CaM trapping" phenomenon, wherein Ca2+/CaM may structurally exclude phosphatase binding and thereby prolong CaMKII autophosphorylation. We conclude that structural protection of autophosphorylated CaMKII by Ca2+/CaM may be an important mechanism for regulation of synaptic plasticity.

Women are underrepresented in computational biology: An analysis of the scholarly literature in biology, computer science and computational biology.

While women are generally underrepresented in STEM fields, there are noticeable differences between fields. For instance, the gender ratio in biology is more balanced than in computer science. We were interested in how this difference is reflected in the interdisciplinary field of computational/quantitative biology. To this end, we examined the proportion of female authors in publications from the PubMed and arXiv databases. There are fewer female authors on research papers in computational biology, as compared to biology in general. This is true across authorship position, year, and journal impact factor. A comparison with arXiv shows that quantitative biology papers have a higher ratio of female authors than computer science papers, placing computational biology in between its two parent fields in terms of gender representation. Both in biology and in computational biology, a female last author increases the probability of other authors on the paper being female, pointing to a potential role of female PIs in influencing the gender balance.

Cooperativity: a competition of definitions.

Cooperativity as a concept is easy to grasp intuitively, but surprisingly hard to define. Two recent papers shed light on the issue and continue the debate on how best to define cooperative binding.

The systems biology format converter.

BACKGROUND: Interoperability between formats is a recurring problem in systems biology research. Many tools have been developed to convert computational models from one format to another. However, they have been developed independently, resulting in redundancy of efforts and lack of synergy. RESULTS: Here we present the System Biology Format Converter (SBFC), which provide a generic framework to potentially convert any format into another. The framework currently includes several converters translating between the following formats: SBML, BioPAX, SBGN-ML, Matlab, Octave, XPP, GPML, Dot, MDL and APM. This software is written in Java and can be used as a standalone executable or web service. CONCLUSIONS: The SBFC framework is an evolving software project. Existing converters can be used and improved, and new converters can be easily added, making SBFC useful to both modellers and developers. The source code and documentation of the framework are freely available from the project web site.

The quantitative methods boot camp: teaching quantitative thinking and computing skills to graduate students in the life sciences.

The past decade has seen a rapid increase in the ability of biologists to collect large amounts of data. It is therefore vital that research biologists acquire the necessary skills during their training to visualize, analyze, and interpret such data. To begin to meet this need, we have developed a "boot camp" in quantitative methods for biology graduate students at Harvard Medical School. The goal of this short, intensive course is to enable students to use computational tools to visualize and analyze data, to strengthen their computational thinking skills, and to simulate and thus extend their intuition about the behavior of complex biological systems. The boot camp teaches basic programming using biological examples from statistics, image processing, and data analysis. This integrative approach to teaching programming and quantitative reasoning motivates students' engagement by demonstrating the relevance of these skills to their work in life science laboratories. Students also have the opportunity to analyze their own data or explore a topic of interest in more detail. The class is taught with a mixture of short lectures, Socratic discussion, and in-class exercises. Students spend approximately 40% of their class time working through both short and long problems. A high instructor-to-student ratio allows students to get assistance or additional challenges when needed, thus enhancing the experience for students at all levels of mastery. Data collected from end-of-course surveys from the last five offerings of the course (between 2012 and 2014) show that students report high learning gains and feel that the course prepares them for solving quantitative and computational problems they will encounter in their research. We outline our course here which, together with the course materials freely available online under a Creative Commons License, should help to facilitate similar efforts by others.

Multi-state modeling of biomolecules.

Multi-state modeling of biomolecules refers to a series of techniques used to represent and compute the behavior of biological molecules or complexes that can adopt a large number of possible functional states. Biological signaling systems often rely on complexes of biological macromolecules that can undergo several functionally significant modifications that are mutually compatible. Thus, they can exist in a very large number of functionally different states. Modeling such multi-state systems poses two problems: the problem of how to describe and specify a multi-state system (the "specification problem") and the problem of how to use a computer to simulate the progress of the system over time (the "computation problem"). To address the specification problem, modelers have in recent years moved away from explicit specification of all possible states and towards rule-based formalisms that allow for implicit model specification, including the κ-calculus, BioNetGen, the Allosteric Network Compiler, and others. To tackle the computation problem, they have turned to particle-based methods that have in many cases proved more computationally efficient than population-based methods based on ordinary differential equations, partial differential equations, or the Gillespie stochastic simulation algorithm. Given current computing technology, particle-based methods are sometimes the only possible option. Particle-based simulators fall into two further categories: nonspatial simulators, such as StochSim, DYNSTOC, RuleMonkey, and the Network-Free Stochastic Simulator (NFSim), and spatial simulators, including Meredys, SRSim, and MCell. Modelers can thus choose from a variety of tools, the best choice depending on the particular problem. Development of faster and more powerful methods is ongoing, promising the ability to simulate ever more complex signaling processes in the future.

Cooperative binding.

Molecular binding is an interaction between molecules that results in a stable association between those molecules. Cooperative binding occurs if the number of binding sites of a macromolecule that are occupied by a specific type of ligand is a nonlinear function of this ligand's concentration. This can be due, for instance, to an affinity for the ligand that depends on the amount of ligand bound. Cooperativity can be positive (supralinear) or negative (infralinear). Cooperative binding is most often observed in proteins, but nucleic acids can also exhibit cooperative binding, for instance of transcription factors. Cooperative binding has been shown to be the mechanism underlying a large range of biochemical and physiological processes.

Validation Through Collaboration: Encouraging Team Efforts to Ensure Internal and External Validity of Computational Models of Biochemical Pathways.

Computational modelling of biochemical reaction pathways is an increasingly important part of neuroscience research. In order to be useful, computational models need to be valid in two senses: First, they need to be consistent with experimental data and able to make testable predictions (external validity). Second, they need to be internally consistent and independently reproducible (internal validity). Here, we discuss both types of validity and provide a brief overview of tools and technologies used to ensure they are met. We also suggest the introduction of new collaborative technologies to ensure model validity: an incentivised experimental database for external validity and reproducibility audits for internal validity. Both rely on FAIR principles and on collaborative science practices.

An allosteric model of calmodulin explains differential activation of PP2B and CaMKII.

Calmodulin plays a vital role in mediating bidirectional synaptic plasticity by activating either calcium/calmodulin-dependent protein kinase II (CaMKII) or protein phosphatase 2B (PP2B) at different calcium concentrations. We propose an allosteric model for calmodulin activation, in which binding to calcium facilitates the transition between a low-affinity [tense (T)] and a high-affinity [relaxed (R)] state. The four calcium-binding sites are assumed to be nonidentical. The model is consistent with previously reported experimental data for calcium binding to calmodulin. It also accounts for known properties of calmodulin that have been difficult to model so far, including the activity of nonsaturated forms of calmodulin (we predict the existence of open conformations in the absence of calcium), an increase in calcium affinity once calmodulin is bound to a target, and the differential activation of CaMKII and PP2B depending on calcium concentration.

JD419, a Staphylococcus aureus Phage With a Unique Morphology and Broad Host Range.

Phage therapy represents a possible treatment option to cure infections caused by multidrug-resistant bacteria, including methicillin and vancomycin-resistant Staphylococcus aureus, to which most antibiotics have become ineffective. In the present study, we report the isolation and complete characterization of a novel phage named JD219 exhibiting a broad host range able to infect 61 of 138 clinical strains of S. aureus tested, which included MRSA strains as well. The phage JD419 exhibits a unique morphology with an elongated capsid and a flexible tail. To evaluate the potential of JD419 to be used as a therapeutic phage, we tested the ability of the phage particles to remain infectious after treatment exceeding physiological pH or temperature. The activity was retained at pH values of 6.0-8.0 and below 50°C. As phages can contain virulence genes, JD419's complete genome was sequenced. The 45509 bp genome is predicted to contain 65 ORFs, none of which show homology to any known virulence or antibiotic resistance genes. Genome analysis indicates that JD419 is a temperate phage, despite observing rapid replication and lysis of host strains. Following the recent advances in synthetic biology, JD419 can be modified by gene engineering to remove prophage-related genes, preventing potential lysogeny, in order to be deployed as a therapeutic phage.

High neural activity accelerates the decline of cognitive plasticity with age in Caenorhabditis elegans.

The ability to learn progressively declines with age. Neural hyperactivity has been implicated in impairing cognitive plasticity with age, but the molecular mechanisms remain elusive. Here, we show that chronic excitation of the Caenorhabditis elegans O2-sensing neurons during ageing causes a rapid decline of experience-dependent plasticity in response to environmental O2 concentration, whereas sustaining lower activity of O2-sensing neurons retains plasticity with age. We demonstrate that neural activity alters the ageing trajectory in the transcriptome of O2-sensing neurons, and our data suggest that high-activity neurons redirect resources from maintaining plasticity to sustaining continuous firing. Sustaining plasticity with age requires the K+-dependent Na+/Ca2+ (NCKX) exchanger, whereas the decline of plasticity with age in high-activity neurons acts through calmodulin and the scaffold protein Kidins220. Our findings demonstrate directly that the activity of neurons alters neuronal homeostasis to govern the age-related decline of neural plasticity and throw light on the mechanisms involved.

The first 10 years of the international coordination network for standards in systems and synthetic biology (COMBINE).

This paper presents a report on outcomes of the 10th Computational Modeling in Biology Network (COMBINE) meeting that was held in Heidelberg, Germany, in July of 2019. The annual event brings together researchers, biocurators and software engineers to present recent results and discuss future work in the area of standards for systems and synthetic biology. The COMBINE initiative coordinates the development of various community standards and formats for computational models in the life sciences. Over the past 10 years, COMBINE has brought together standard communities that have further developed and harmonized their standards for better interoperability of models and data. COMBINE 2019 was co-located with a stakeholder workshop of the European EU-STANDS4PM initiative that aims at harmonized data and model standardization for in silico models in the field of personalized medicine, as well as with the FAIRDOM PALs meeting to discuss findable, accessible, interoperable and reusable (FAIR) data sharing. This report briefly describes the work discussed in invited and contributed talks as well as during breakout sessions. It also highlights recent advancements in data, model, and annotation standardization efforts. Finally, this report concludes with some challenges and opportunities that this community will face during the next 10 years.

FindSim: A Framework for Integrating Neuronal Data and Signaling Models.

Current experiments touch only small but overlapping parts of very complex subcellular signaling networks in neurons. Even with modern optical reporters and pharmacological manipulations, a given experiment can only monitor and control a very small subset of the diverse, multiscale processes of neuronal signaling. We have developed FindSim (Framework for Integrating Neuronal Data and SIgnaling Models) to anchor models to structured experimental datasets. FindSim is a framework for integrating many individual electrophysiological and biochemical experiments with large, multiscale models so as to systematically refine and validate the model. We use a structured format for encoding the conditions of many standard physiological and pharmacological experiments, specifying which parts of the model are involved, and comparing experiment outcomes with model output. A database of such experiments is run against successive generations of composite cellular models to iteratively improve the model against each experiment, while retaining global model validity. We suggest that this toolchain provides a principled and scalable way to tackle model complexity and diversity of data sources.

Cooperative binding mitigates the high-dose hook effect.

BACKGROUND: The high-dose hook effect (also called prozone effect) refers to the observation that if a multivalent protein acts as a linker between two parts of a protein complex, then increasing the amount of linker protein in the mixture does not always increase the amount of fully formed complex. On the contrary, at a high enough concentration range the amount of fully formed complex actually decreases. It has been observed that allosterically regulated proteins seem less susceptible to this effect. The aim of this study was two-fold: First, to investigate the mathematical basis of how allostery mitigates the prozone effect. And second, to explore the consequences of allostery and the high-dose hook effect using the example of calmodulin, a calcium-sensing protein that regulates the switch between long-term potentiation and long-term depression in neurons. RESULTS: We use a combinatorial model of a "perfect linker protein" (with infinite binding affinity) to mathematically describe the hook effect and its behaviour under allosteric conditions. We show that allosteric regulation does indeed mitigate the high-dose hook effect. We then turn to calmodulin as a real-life example of an allosteric protein. Using kinetic simulations, we show that calmodulin is indeed subject to a hook effect. We also show that this effect is stronger in the presence of the allosteric activator Ca 2+/calmodulin-dependent kinase II (CaMKII), because it reduces the overall cooperativity of the calcium-calmodulin system. It follows that, surprisingly, there are conditions where increased amounts of allosteric activator actually decrease the activity of a protein. CONCLUSIONS: We show that cooperative binding can indeed act as a protective mechanism against the hook effect. This will have implications in vivo where the extent of cooperativity of a protein can be modulated, for instance, by allosteric activators or inhibitors. This can result in counterintuitive effects of decreased activity with increased concentrations of both the allosteric protein itself and its allosteric activators.

Biophysical properties of presynaptic short-term plasticity in hippocampal neurons: insights from electrophysiology, imaging and mechanistic models.

Hippocampal neurons show different types of short-term plasticity (STP). Some exhibit facilitation of their synaptic responses and others depression. In this review we discuss presynaptic biophysical properties behind heterogeneity in STP in hippocampal neurons such as alterations in vesicle priming and docking, fusion, neurotransmitter filling and vesicle replenishment. We look into what types of information electrophysiology, imaging and mechanistic models have given about the time scales and relative impact of the different properties on STP with an emphasis on the use of mechanistic models as complementary tools to experimental procedures. Taken together this tells us that it is possible for a multitude of different mechanisms to underlie the same STP pattern, even though some are more important in specific cases, and that mechanistic models can be used to integrate the biophysical properties to see which mechanisms are more important in specific cases of STP.

Structural analysis and stochastic modelling suggest a mechanism for calmodulin trapping by CaMKII.

Activation of CaMKII by calmodulin and the subsequent maintenance of constitutive activity through autophosphorylation at threonine residue 286 (Thr286) are thought to play a major role in synaptic plasticity. One of the effects of autophosphorylation at Thr286 is to increase the apparent affinity of CaMKII for calmodulin, a phenomenon known as "calmodulin trapping". It has previously been suggested that two binding sites for calmodulin exist on CaMKII, with high and low affinities, respectively. We built structural models of calmodulin bound to both of these sites. Molecular dynamics simulation showed that while binding of calmodulin to the supposed low-affinity binding site on CaMKII is compatible with closing (and hence, inactivation) of the kinase, and could even favour it, binding to the high-affinity site is not. Stochastic simulations of a biochemical model showed that the existence of two such binding sites, one of them accessible only in the active, open conformation, would be sufficient to explain calmodulin trapping by CaMKII. We can explain the effect of CaMKII autophosphorylation at Thr286 on calmodulin trapping: It stabilises the active state and therefore makes the high-affinity binding site accessible. Crucially, a model with only one binding site where calmodulin binding and CaMKII inactivation are strictly mutually exclusive cannot reproduce calmodulin trapping. One of the predictions of our study is that calmodulin binding in itself is not sufficient for CaMKII activation, although high-affinity binding of calmodulin is.