Mechanistic regulation of planarian shape during growth and degrowth.

Adult planarians can grow when fed and degrow (shrink) when starved while maintaining their whole-body shape. It is unknown how the morphogens patterning the planarian axes are coordinated during feeding and starvation or how they modulate the necessary differential tissue growth or degrowth. Here, we investigate the dynamics of planarian shape together with a theoretical study of the mechanisms regulating whole-body proportions and shape. We found that the planarian body proportions scale isometrically following similar linear rates during growth and degrowth, but that fed worms are significantly wider than starved worms. By combining a descriptive model of planarian shape and size with a mechanistic model of anterior-posterior and medio-lateral signaling calibrated with a novel parameter optimization methodology, we theoretically demonstrate that the feedback loop between these positional information signals and the shape they control can regulate the planarian whole-body shape during growth. Furthermore, the computational model produced the correct shape and size dynamics during degrowth as a result of a predicted increase in apoptosis rate and pole signal during starvation. These results offer mechanistic insights into the dynamic regulation of whole-body morphologies.

Democratizing protein language models with parameter-efficient fine-tuning.

Proteomics has been revolutionized by large protein language models (PLMs), which learn unsupervised representations from large corpora of sequences. These models are typically fine-tuned in a supervised setting to adapt the model to specific downstream tasks. However, the computational and memory footprint of fine-tuning (FT) large PLMs presents a barrier for many research groups with limited computational resources. Natural language processing has seen a similar explosion in the size of models, where these challenges have been addressed by methods for parameter-efficient fine-tuning (PEFT). In this work, we introduce this paradigm to proteomics through leveraging the parameter-efficient method LoRA and training new models for two important tasks: predicting protein-protein interactions (PPIs) and predicting the symmetry of homooligomer quaternary structures. We show that these approaches are competitive with traditional FT while requiring reduced memory and substantially fewer parameters. We additionally show that for the PPI prediction task, training only the classification head also remains competitive with full FT, using five orders of magnitude fewer parameters, and that each of these methods outperform state-of-the-art PPI prediction methods with substantially reduced compute. We further perform a comprehensive evaluation of the hyperparameter space, demonstrate that PEFT of PLMs is robust to variations in these hyperparameters, and elucidate where best practices for PEFT in proteomics differ from those in natural language processing. All our model adaptation and evaluation code is available open-source at https://github.com/microsoft/peft_proteomics. Thus, we provide a blueprint to democratize the power of PLM adaptation to groups with limited computational resources.

Biologically meaningful regulatory logic enhances the convergence rate in Boolean networks and bushiness of their state transition graph.

Boolean models of gene regulatory networks (GRNs) have gained widespread traction as they can easily recapitulate cellular phenotypes via their attractor states. Their overall dynamics are embodied in a state transition graph (STG). Indeed, two Boolean networks (BNs) with the same network structure and attractors can have drastically different STGs depending on the type of Boolean functions (BFs) employed. Our objective here is to systematically delineate the effects of different classes of BFs on the structural features of the STG of reconstructed Boolean GRNs while keeping network structure and biological attractors fixed, and explore the characteristics of BFs that drive those features. Using $10$ reconstructed Boolean GRNs, we generate ensembles that differ in BFs and compute from their STGs the dynamics' rate of contraction or 'bushiness' and rate of 'convergence', quantified with measures inspired from cellular automata (CA) that are based on the garden-of-Eden (GoE) states. We find that biologically meaningful BFs lead to higher STG 'bushiness' and 'convergence' than random ones. Obtaining such 'global' measures gets computationally expensive with larger network sizes, stressing the need for feasible proxies. So we adapt Wuensche's $Z$-parameter in CA to BFs in BNs and provide four natural variants, which, along with the average sensitivity of BFs computed at the network level, comprise our descriptors of local dynamics and we find some of them to be good proxies for bushiness. Finally, we provide an excellent proxy for the 'convergence' based on computing transient lengths originating at random states rather than GoE states.

Neural parameter calibration for large-scale multiagent models.

Computational models have become a powerful tool in the quantitative sciences to understand the behavior of complex systems that evolve in time. However, they often contain a potentially large number of free parameters whose values cannot be obtained from theory but need to be inferred from data. This is especially the case for models in the social sciences, economics, or computational epidemiology. Yet, many current parameter estimation methods are mathematically involved and computationally slow to run. In this paper, we present a computationally simple and fast method to retrieve accurate probability densities for model parameters using neural differential equations. We present a pipeline comprising multiagent models acting as forward solvers for systems of ordinary or stochastic differential equations and a neural network to then extract parameters from the data generated by the model. The two combined create a powerful tool that can quickly estimate densities on model parameters, even for very large systems. We demonstrate the method on synthetic time series data of the SIR model of the spread of infection and perform an in-depth analysis of the Harris-Wilson model of economic activity on a network, representing a nonconvex problem. For the latter, we apply our method both to synthetic data and to data of economic activity across Greater London. We find that our method calibrates the model orders of magnitude more accurately than a previous study of the same dataset using classical techniques, while running between 195 and 390 times faster.

Neuronal morphology enhances robustness to perturbations of channel densities.

Biological neurons show significant cell-to-cell variability but have the striking ability to maintain their key firing properties in the face of unpredictable perturbations and stochastic noise. Using a population of multi-compartment models consisting of soma, neurites, and axon for the lateral pyloric neuron in the crab stomatogastric ganglion, we explore how rebound bursting is preserved when the 14 channel conductances in each model are all randomly varied. The coupling between the axon and other compartments is critical for the ability of the axon to spike during bursts and consequently determines the set of successful solutions. When the coupling deviates from a biologically realistic range, the neuronal tolerance of conductance variations is lessened. Thus, the gross morphological features of these neurons enhance their robustness to perturbations of channel densities and expand the space of individual variability that can maintain a desired output pattern.

A strategy combining denoising and cryo-EM single particle analysis.

In cryogenic electron microscopy (cryo-EM) single particle analysis (SPA), high-resolution three-dimensional structures of biological macromolecules are determined by iteratively aligning and averaging a large number of two-dimensional projections of molecules. Since the correlation measures are sensitive to the signal-to-noise ratio, various parameter estimation steps in SPA will be disturbed by the high-intensity noise in cryo-EM. However, denoising algorithms tend to damage high frequencies and suppress mid- and high-frequency contrast of micrographs, which exactly the precise parameter estimation relies on, therefore, limiting their application in SPA. In this study, we suggest combining a cryo-EM image processing pipeline with denoising and maximizing the signal's contribution in various parameter estimation steps. To solve the inherent flaws of denoising algorithms, we design an algorithm named MScale to correct the amplitude distortion caused by denoising and propose a new orientation determination strategy to compensate for the high-frequency loss. In the experiments on several real datasets, the denoised particles are successfully applied in the class assignment estimation and orientation determination tasks, ultimately enhancing the quality of biomacromolecule reconstruction. The case study on classification indicates that our strategy not only improves the resolution of difficult classes (up to 5 Å) but also resolves an additional class. In the case study on orientation determination, our strategy improves the resolution of the final reconstructed density map by 0.34 Å compared with conventional strategy. The code is available at https://github.com/zhanghui186/Mscale.

Sequential Bayesian Phylogenetic Inference.

The ideal approach to Bayesian phylogenetic inference is to estimate all parameters of interest jointly in a single hierarchical model. However, this is often not feasible in practice due to the high computational cost. Instead, phylogenetic pipelines generally consist of sequential analyses, whereby a single point estimate from a given analysis is used as input for the next analysis (e.g., a single multiple sequence alignment is used to estimate a gene tree). In this framework, uncertainty is not propagated from step to step, which can lead to inaccurate or spuriously confident results. Here, we formally develop and test a sequential inference approach for Bayesian phylogenetic inference, which uses importance sampling to generate observations for the next step of an analysis pipeline from the posterior distribution produced in the previous step. Our sequential inference approach presented here not only accounts for uncertainty between analysis steps, but also allows for greater flexibility in software choice (and hence model availability) and can be computationally more efficient than the traditional joint inference approach when multiple models are being tested. We show that our sequential inference approach is identical in practice to the joint inference approach only if sufficient information in the data is present (a narrow posterior distribution) and/or sufficiently many importance samples are used. Conversely, we show that the common practice of using a single point estimate can be biased, e.g., a single phylogeny estimate to transform an unrooted phylogeny into a time-calibrated phylogeny. We demonstrate the theory of sequential Bayesian inference using both a toy example and an empirical case study of divergence-time estimation in insects using a relaxed clock model from transcriptome data. In the empirical example, we estimate three posterior distributions of branch lengths from the same data (DNA character matrix with a GTR+Γ+I substitution model, an amino acid data matrix with empirical substitution models, and an amino acid data matrix with the PhyloBayes CAT-GTR model). Finally, we apply three different node-calibration strategies and show that divergence-time estimates are affected by both the data source and underlying substitution process to estimate branch lengths as well as the node-calibration strategies. Thus, our new sequential Bayesian phylogenetic inference provides the opportunity to efficiently test different approaches for divergence time estimation, including branch-length estimation from other software.

Assessing parameter efficient methods for pre-trained language model in annotating scRNA-seq data.

Annotating cell types of single-cell RNA sequencing (scRNA-seq) data is crucial for studying cellular heterogeneity in the tumor microenvironment. Recently, large-scale pre-trained language models (PLMs) have achieved significant progress in cell-type annotation of scRNA-seq data. This approach effectively addresses previous methods' shortcomings in performance and generalization. However, fine-tuning PLMs for different downstream tasks demands considerable computational resources, rendering it impractical. Hence, a new research branch introduces parameter-efficient fine-tuning (PEFT). This involves optimizing a few parameters while leaving the majority unchanged, leading to substantial reductions in computational expenses. Here, we utilize scBERT, a large-scale pre-trained model, to explore the capabilities of three PEFT methods in scRNA-seq cell type annotation. Extensive benchmark studies across several datasets demonstrate the superior applicability of PEFT methods. Furthermore, downstream analysis using models obtained through PEFT showcases their utility in novel cell type discovery and model interpretability for potential marker genes. Our findings underscore the considerable potential of PEFT in PLM-based cell type annotation, presenting novel perspectives for the analysis of scRNA-seq data.

Modeling for COVID-19 college reopening decisions: Cornell, a case study.

We consider epidemiological modeling for the design of COVID-19 interventions in university populations, which have seen significant outbreaks during the pandemic. A central challenge is sensitivity of predictions to input parameters coupled with uncertainty about these parameters. Nearly 2 y into the pandemic, parameter uncertainty remains because of changes in vaccination efficacy, viral variants, and mask mandates, and because universities' unique characteristics hinder translation from the general population: a high fraction of young people, who have higher rates of asymptomatic infection and social contact, as well as an enhanced ability to implement behavioral and testing interventions. We describe an epidemiological model that formed the basis for Cornell University's decision to reopen for in-person instruction in fall 2020 and supported the design of an asymptomatic screening program instituted concurrently to prevent viral spread. We demonstrate how the structure of these decisions allowed risk to be minimized despite parameter uncertainty leading to an inability to make accurate point estimates and how this generalizes to other university settings. We find that once-per-week asymptomatic screening of vaccinated undergraduate students provides substantial value against the Delta variant, even if all students are vaccinated, and that more targeted testing of the most social vaccinated students provides further value.

Comparing methods for statistical inference with model uncertainty.

Probability models are used for many statistical tasks, notably parameter estimation, interval estimation, inference about model parameters, point prediction, and interval prediction. Thus, choosing a statistical model and accounting for uncertainty about this choice are important parts of the scientific process. Here we focus on one such choice, that of variables to include in a linear regression model. Many methods have been proposed, including Bayesian and penalized likelihood methods, and it is unclear which one to use. We compared 21 of the most popular methods by carrying out an extensive set of simulation studies based closely on real datasets that span a range of situations encountered in practical data analysis. Three adaptive Bayesian model averaging (BMA) methods performed best across all statistical tasks. These used adaptive versions of Zellner's g-prior for the parameters, where the prior variance parameter g is a function of sample size or is estimated from the data. We found that for BMA methods implemented with Markov chain Monte Carlo, 10,000 iterations were enough. Computationally, we found two of the three best methods (BMA with g=√n and empirical Bayes-local) to be competitive with the least absolute shrinkage and selection operator (LASSO), which is often preferred as a variable selection technique because of its computational efficiency. BMA performed better than Bayesian model selection (in which just one model is selected).

On-Chip Topological Photonic Crystal Nanobeam Filters.

We present an on-chip filter with a broad tailorable working wavelength and a single-mode operation. This is realized through the application of topological photonic crystal nanobeam filters employing synthesis parameter dimensions. By introducing the translation of air holes as a new synthetic parameter dimension, we obtained nanobeams with tunable Zak phases. Leveraging the bulk-edge correspondence, we identify the existence of topological cavity modes and establish a correlation between the cavity's interface morphology and working wavelength. Through experiments, we demonstrate filters with adjustable filtering wavelengths ranging from 1301 to 1570 nm. Our work illustrates the use of the synthetic translation dimension in the design of on-chip filters, and it holds potential for applications in other devices such as microcavities.

Improving Proteomic Identification Using Narrow Isolation Windows with Zeno SWATH Data-Independent Acquisition.

Data-independent acquisition (DIA) techniques such as sequential window acquisition of all theoretical mass spectra (SWATH) acquisition have emerged as the preferred strategies for proteomic analyses. Our study optimized the SWATH-DIA method using a narrow isolation window placement approach, improving its proteomic performance. We optimized the acquisition parameter combinations of narrow isolation windows with different widths (1.9 and 2.9 Da) on a ZenoTOF 7600 (Sciex); the acquired data were analyzed using DIA-NN (version 1.8.1). Narrow SWATH (nSWATH) identified 5916 and 7719 protein groups on the digested peptides, corresponding to 400 ng of protein from mouse liver and HEK293T cells, respectively, improving identification by 7.52 and 4.99%, respectively, compared to conventional SWATH. The median coefficient of variation of the quantified values was less than 6%. We further analyzed 200 ng of benchmark samples comprising peptides from known ratios ofEscherichia coli, yeast, and human peptides using nSWATH. Consequently, it achieved accuracy and precision comparable to those of conventional SWATH, identifying an average of 95,456 precursors and 9342 protein groups across three benchmark samples, representing 12.6 and 9.63% improved identification compared to conventional SWATH. The nSWATH method improved identification at various loading amounts of benchmark samples, identifying 40.7% more protein groups at 25 ng. These results demonstrate the improved performance of nSWATH, contributing to the acquisition of deeper proteomic data from complex biological samples.

Phase plane dynamics of ERK phosphorylation.

The extracellular signal-regulated kinase (ERK) controls multiple critical processes in the cell and is deregulated in human cancers, congenital abnormalities, immune diseases, and neurodevelopmental syndromes. Catalytic activity of ERK requires dual phosphorylation by an upstream kinase, in a mechanism that can be described by two sequential Michaelis-Menten steps. The estimation of individual reaction rate constants from kinetic data in the full mechanism has proved challenging. Here, we present an analytically tractable approach to parameter estimation that is based on the phase plane representation of ERK activation and yields two combinations of six reaction rate constants in the detailed mechanism. These combinations correspond to the ratio of the specificities of two consecutive phosphorylations and the probability that monophosphorylated substrate does not dissociate from the enzyme before the second phosphorylation. The presented approach offers a language for comparing the effects of mutations that disrupt ERK activation and function in vivo. As an illustration, we use phase plane representation to analyze dual phosphorylation under heterozygous conditions, when two enzyme variants compete for the same substrate.

Addressing current limitations of household transmission studies by collecting contact data.

Modeling studies of household transmission data have helped characterize the role of children in influenza and COVID-19 epidemics. However, estimates from these studies may be biased since they do not account for the heterogeneous nature of household contacts. Here, we quantified the impact of contact heterogeneity between household members on the estimation of child relative susceptibility and infectivity. We simulated epidemics of SARS-CoV-2-like and influenza-like infections in a synthetic population of 1,000 households assuming heterogeneous contact levels. Relative contact frequencies were derived from a household contact study according to which contacts are more frequent in the father-mother pair, followed by the child-mother, child-child, and finally child-father pairs. Child susceptibility and infectivity were then estimated while accounting for heterogeneous contacts or not. When ignoring contact heterogeneity, child relative susceptibility was underestimated by approximately 20% in the two disease scenarios. Child relative infectivity was underestimated by 20% when children and adults had different infectivity levels. These results are sensitive to our assumptions of European-style household contact patterns; but they highlight that household studies collecting both disease and contact data are needed to assess the role of complex household contact behavior on disease transmission and improve estimation of key biological parameters.

Prediction of order parameters based on protein NMR structure ensemble and machine learning.

The fast motions of proteins at the picosecond to nanosecond timescale, known as fast dynamics, are closely related to protein conformational entropy and rearrangement, which in turn affect catalysis, ligand binding and protein allosteric effects. The most used NMR approach to study fast protein dynamics is the model free method, which uses order parameter S2 to describe the amplitude of the internal motion of local group. However, to obtain order parameter through NMR experiments is quite complex and lengthy. In this paper, we present a machine learning approach for predicting backbone 1H-15N order parameters based on protein NMR structure ensemble. A random forest model is used to learn the relationship between order parameters and structural features. Our method achieves high accuracy in predicting backbone 1H-15N order parameters for a test dataset of 10 proteins, with a Pearson correlation coefficient of 0.817 and a root-mean-square error of 0.131.

Vibration-Controlled Transient Elastography-Based Parameters Predict Clinical Outcomes in Liver Transplant Recipients.

BACKGROUND AND AIMS: Vibration-controlled transient elastography (VCTE) is used in clinical practice to risk-stratify liver transplant (LT) recipients; however, there are currently little data demonstrating the relationship between VCTE and clinical outcomes. METHODS: A total of 362 adult LT recipients with successful VCTE examination between 2015 and 2022 were included. Presence of advanced fibrosis was defined as liver stiffness measurement (LSM) ≥10.5 kPa and hepatic steatosis as controlled attenuation parameter (CAP) ≥270 dB/m. The outcomes of interest included all-cause mortality, myocardial infarction (MI), and graft cirrhosis using cumulative incidence analysis that accounted for the competing risks of these outcomes. RESULTS: The LSM was elevated in 64 (18%) and CAP in 163 (45%) LT recipients. The baseline LSM values were similar in patients with elevated vs normal CAP values. After a median follow-up of 65 (interquartile range, 20-140) months from LT to baseline VCTE, 66 (18%) patients died, 12 (3%) developed graft cirrhosis, and 18 (5%) experienced an MI. Baseline high LSM was independently associated with all-cause mortality (hazard ratio [HR], 1.97; 95% confidence interval [CI], 1.11-3.50; P = .02) and new onset cirrhosis (HR, 6.74; 95% CI, 2.08-21.79; P < .01). A higher CAP value was significantly and independently associated with increased risk of experiencing a MI over study follow-up (HR, 4.14; 95% CI, 1.29-13.27; P = .017). CONCLUSIONS: The VCTE-based parameters are associated with clinical outcomes and offer the potential to be incorporated into clinical risk-stratification strategies to improve outcomes among LT recipients.

Mathematical modeling of the interaction between yolk utilization and fish growth in zebrafish, Danio rerio.

Optimal embryonic development plays a major role in the health of an individual beyond the developmental stage. Nutritional perturbation during development is associated with cardiovascular and metabolic disease later in life. With both nutritional uptake and overall growth being risk factors for eventual health, it is necessary to understand not only the behavior of the processes during development but also their interactions. In this study, we used differential equations, image analyses, curve fittings, parameter estimation and laboratory experiments to quantify the rate of yolk absorption and its effect on early development of a vertebrate model (Danio rerio). Findings from this study establish a nonlinear functional relationship between nutrient absorption and early fish growth. We found that the rate of change in fish length and yolk utilization is logistic, that is the yolk decays rapidly for a period of time before leveling out. An interesting finding from this study is that yolk utilization reaches its maximum at 84 h post-fertilization. We validated our mathematical models against experimental observations, making them powerful tools for replication and future simulations.

Identifying dryland-resilient chickpea genotypes for autumn sowing, with a focus on multi-trait stability parameters and biochemical enzyme activity.

BACKGROUND: Chickpea is a key pulse crop grown in the spring in dryland regions. The cold resistance potential of chickpeas allows for the development of genotypes with varying sowing dates to take advantage of autumn and winter rainfall, particularly in dryland regions. In this study, we assessed grain yield, plant height, 100-seed weight, days to maturity, and days to flowering of 17 chickpea genotypes in five autumn-sown dryland regions from 2019 to 2021. Additionally, the response of selected chickpea genotypes to cold stress was examined at temperatures of -4 °C, 4 °C, and 22 °C by analyzing biochemical enzymes. RESULTS: Mixed linear model of ANOVA revealed a significant genotype × environment interaction for all traits measured, indicating varying reactions of genotypes across test environments. This study reported low estimates of broad-sense heritability for days to flowering (0.34), days to maturity (0.13), and grain yield (0.08). Plant height and seed weight exhibited the highest heritability, with genotypic selection accuracies of 0.73 and 0.92, respectively. Moreover, partial least square regression highlighted the impactful role of rainfall during all months except of October, November, and February on grain yield and its interaction with environments in autumn-planted chickpeas. Among the genotypes studied, G9, G10, and G17 emerged as superior based on stability parameters and grain yield. In particular, genotype G9 stood out as a promising genotype for dryland regions, considering both MTSI and genotype by yield*trait aproaches. The cold assay indicated that - 4 °C is crucial for distinguishing between susceptible and resistant genotypes. The results showed the important role of the enzymes CAT and GPX in contributing to the cold tolerance of genotype G9 in autumn-sown chickpeas. CONCLUSIONS: Significant G×E for agro-morphological traits of chickpea shows prerequisite for multi-trial analysis. Chickpea`s direct root system cause that monthly rainfall during plant establishment has no critical role in its yield interaction with dryland environment. Considering the importance of agro-morphological traits and their direct and indirect effects on grain yield, the utilization of multiple-trait stability approches is propose. Evaluation of chickpea germplasm reaction against cold stress is necessary for autumn-sowing. Finally, autumn sowing of genotype FLIP 10-128 C in dryland conditions can led to significant crop performance.

Transcription factor TaNF-YB2 interacts with partners TaNF-YA7/YC7 and transcriptionally activates distinct stress-defensive genes to modulate drought tolerance in T. Aestivum.

BACKGROUND: Drought stress limits significantly the crop productivity. However, plants have evolved various strategies to cope with the drought conditions by adopting complex molecular, biochemical, and physiological mechanisms. Members of the nuclear factor Y (NF-Y) transcription factor (TF) family constitute one of the largest TF classes and are involved in plant responses to abiotic stresses. RESULTS: TaNF-YB2, a NY-YB subfamily gene in T. aestivum, was characterized in this study focusing on its role in mediating plant adaptation to drought stress. Yeast two-hybrid (Y-2 H), biomolecular fluoresence complementation (BiFC), and Co-immunoprecipitation (Co-IP) assays indicated that TaNF-YB2 interacts with the NF-YA member TaNF-YA7 and NF-YC family member TaNF-YC7, which constitutes a heterotrimer TaNF-YB2/TaNF-YA7/TaNF-YC7. The TaNF-YB2 transcripts are induced in roots and aerial tissues upon drought signaling; GUS histochemical staining analysis demonstrated the roles of cis-regulatory elements ABRE and MYB situated in TaNF-YB2 promoter to contribute to target gene response to drought. Transgene analysis on TaNF-YB2 confirmed its functions in regulating drought adaptation via modulating stomata movement, osmolyte biosynthesis, and reactive oxygen species (ROS) homeostasis. TaNF-YB2 possessed the abilities in transcriptionally activating TaP5CS2, the P5CS family gene involving proline biosynthesis and TaSOD1, TaCAT5, and TaPOD5, the genes encoding antioxidant enzymes. Positive correlations were found between yield and the TaNF-YB2 transcripts in a core panel constituting 45 wheat cultivars under drought condition, in which two types of major haplotypes including TaNF-YB2-Hap1 and -Hap2 were included, with the former conferring more TaNF-YB2 transcripts and stronger plant drought tolerance. CONCLUSIONS: TaNF-YB2 is transcriptional response to drought stress. It is an essential regulator in mediating plant drought adaptation by modulating the physiological processes associated with stomatal movement, osmolyte biosynthesis, and reactive oxygen species (ROS) homeostasis, depending on its role in transcriptionally regulating stress response genes. Our research deepens the understanding of plant drought stress underlying NF-Y TF family and provides gene resource in efforts for molecular breeding the drought-tolerant cultivars in T. aestivum.

A protocol for dynamic model calibration.

Ordinary differential equation models are nowadays widely used for the mechanistic description of biological processes and their temporal evolution. These models typically have many unknown and nonmeasurable parameters, which have to be determined by fitting the model to experimental data. In order to perform this task, known as parameter estimation or model calibration, the modeller faces challenges such as poor parameter identifiability, lack of sufficiently informative experimental data and the existence of local minima in the objective function landscape. These issues tend to worsen with larger model sizes, increasing the computational complexity and the number of unknown parameters. An incorrectly calibrated model is problematic because it may result in inaccurate predictions and misleading conclusions. For nonexpert users, there are a large number of potential pitfalls. Here, we provide a protocol that guides the user through all the steps involved in the calibration of dynamic models. We illustrate the methodology with two models and provide all the code required to reproduce the results and perform the same analysis on new models. Our protocol provides practitioners and researchers in biological modelling with a one-stop guide that is at the same time compact and sufficiently comprehensive to cover all aspects of the problem.