A Step-by-Step Guide to Using BioNetFit.

BioNetFit is a software tool designed for solving parameter identification problems that arise in the development of rule-based models. It solves these problems through curve fitting (i.e., nonlinear regression). BioNetFit is compatible with deterministic and stochastic simulators that accept BioNetGen language (BNGL)-formatted files as inputs, such as those available within the BioNetGen framework. BioNetFit can be used on a laptop or stand-alone multicore workstation as well as on many Linux clusters, such as those that use the Slurm Workload Manager to schedule jobs. BioNetFit implements a metaheuristic population-based global optimization procedure, an evolutionary algorithm (EA), to minimize a user-defined objective function, such as a residual sum of squares (RSS) function. BioNetFit also implements a bootstrapping procedure for determining confidence intervals for parameter estimates. Here, we provide step-by-step instructions for using BioNetFit to estimate the values of parameters of a BNGL-encoded model and to define bootstrap confidence intervals. The process entails the use of several plain-text files, which are processed by BioNetFit and BioNetGen. In general, these files include (1) one or more EXP files, which each contains (experimental) data to be used in parameter identification/bootstrapping; (2) a BNGL file containing a model section, which defines a (rule-based) model, and an actions section, which defines simulation protocols that generate GDAT and/or SCAN files with model predictions corresponding to the data in the EXP file(s); and (3) a CONF file that configures the fitting/bootstrapping job and that defines algorithmic parameter settings.

Further along the Road Less Traveled: AMBER ff15ipq, an Original Protein Force Field Built on a Self-Consistent Physical Model.

We present the AMBER ff15ipq force field for proteins, the second-generation force field developed using the Implicitly Polarized Q (IPolQ) scheme for deriving implicitly polarized atomic charges in the presence of explicit solvent. The ff15ipq force field is a complete rederivation including more than 300 unique atomic charges, 900 unique torsion terms, 60 new angle parameters, and new atomic radii for polar hydrogens. The atomic charges were derived in the context of the SPC/Eb water model, which yields more-accurate rotational diffusion of proteins and enables direct calculation of nuclear magnetic resonance (NMR) relaxation parameters from molecular dynamics simulations. The atomic radii improve the accuracy of modeling salt bridge interactions relative to contemporary fixed-charge force fields, rectifying a limitation of ff14ipq that resulted from its use of pair-specific Lennard-Jones radii. In addition, ff15ipq reproduces penta-alanine J-coupling constants exceptionally well, gives reasonable agreement with NMR relaxation rates, and maintains the expected conformational propensities of structured proteins/peptides, as well as disordered peptides-all on the microsecond (µs) time scale, which is a critical regime for drug design applications. These encouraging results demonstrate the power and robustness of our automated methods for deriving new force fields. All parameters described here and the mdgx program used to fit them are included in the AmberTools16 distribution.

The receptor tyrosine kinase Pvr promotes tissue closure by coordinating corpse removal and epidermal zippering.

A leading cause of human birth defects is the incomplete fusion of tissues, often manifested in the palate, heart or neural tube. To investigate the molecular control of tissue fusion, embryonic dorsal closure and pupal thorax closure in Drosophila are useful experimental models. We find that Pvr mutants have defects in dorsal midline closure with incomplete amnioserosa internalization and epidermal zippering, as well as cardia bifida. These defects are relatively mild in comparison to those seen with other signaling mutants, such as in the JNK pathway, and we demonstrate that JNK signaling is not perturbed by altering Pvr receptor tyrosine kinase activity. Rather, modulation of Pvr levels in the ectoderm has an impact on PIP3 membrane accumulation, consistent with a link to PI3K signal transduction. Polarized PI3K activity influences protrusive activity from the epidermal leading edge and the protrusion area changes in accord with Pvr signaling intensity, providing a possible mechanism to explain Pvr mutant phenotypes. Tissue-specific rescue experiments indicate a partial requirement in epithelial tissue, but confirm the essential role of Pvr in hemocytes for embryonic survival. Taken together, we argue that inefficient removal of the internalizing amnioserosa tissue by mutant hemocytes coupled with impaired midline zippering of mutant epithelium creates a situation in some embryos whereby dorsal midline closure is incomplete. Based on these observations, we suggest that efferocytosis (corpse clearance) could contribute to proper tissue closure and thus might underlie some congenital birth defects.