Computational modeling to predict mechanical function of joints: application to the lower leg with simulation of two cadaver studies

Computational models of musculoskeletal joints and limbs can provide useful information about joint mechanics. Validated models can be used as predictive devices for understanding joint function and serve as clinical tools for predicting the outcome of surgical procedures. A new computational modeling approach was developed for simulating joint kinematics that are dictated by bone/joint anatomy, ligamentous constraints, and applied loading. Three-dimensional computational models of the lower leg were created to illustrate the application of this new approach. Model development began with generating three-dimensional surfaces of each bone from CT images and then importing into the three-dimensional solid modeling software SOLIDWORKS and motion simulation package COSMOSMOTION. Through SOLIDWORKS and COSMOSMOTION, each bone surface file was filled to create a solid object and positioned necessary components added, and simulations executed. Three-dimensional contacts were added to inhibit intersection of the bones during motion. Ligaments were represented as linear springs. Model predictions were then validated by comparison to two different cadaver studies, syndesmotic injury and repair and ankle inversion following ligament transection. The syndesmotic injury model was able to predict tibial rotation, fibular rotation, and anterior/posterior displacement. In the inversion simulation, calcaneofibular ligament extension and angles of inversion compared well. Some experimental data proved harder to simulate accurately, due to certain software limitations and lack of complete experimental data. Other parameters that could not be easily obtained experimentally can be predicted and analyzed by the computational simulations. In the syndesmotic injury study, the force generated in the tibionavicular and calcaneofibular ligaments reduced with the insertion of the staple, indicating how this repair technique changes joint function. After transection of the calcaneofibular ligament in the inversion stability study, a major increase in force was seen in several of the ligaments on the lateral aspect of the foot and ankle, indicating the recruitment of other structures to permit function after injury. Overall, the computational models were able to predict joint kinematics of the lower leg with particular focus on the ankle complex. This same approach can be taken to create models of other limb segments such as the elbow and wrist. Additional parameters can be calculated in the models that are not easily obtained experimentally such as ligament forces, force transmission across joints, and three-dimensional movement of all bones. Muscle activation can be incorporated in the model through the action of applied forces within the software for future studies.     

Analyzing multi-environment variety trials using randomization-derived mixed models

Of interest is the analysis of results of a series of experiments repeated at several environments with the same set of plant varieties. Suppose that the experiments, multi-environment variety trials, are all conducted in resolvable incomplete block (IB) designs. Following the randomization approach adopted in Caliński and Kageyama (2000, Lecture Notes in Statistics, 150), two models for analyzing such trial data can be considered. One is derived under a complete additivity assumption, the other takes into account possible different responses of the varieties to variable environmental conditions. The analysis under the first, the standard model, does not provide answers to questions related to the performance of the individual varieties at different environments. These can be considered when using the more general second model. The purpose of this article is to devise interesting parameter estimation and hypothesis testing procedures under that more realistic model. Its application is illustrated by a thorough analysis of a set of data from a winter wheat series of trials.     

Estimation of evolutionary parameters with phylogenetic trees

An important issue in the phylogenetic analysis of nucleotide sequence data using the maximum likelihood (ML) method is the underlying evolutionary model employed. We consider the problem of simultaneously estimating the tree topology and the parameters in the underlying substitution model and of obtaining estimates of the standard errors of these parameter estimates. Given a fixed tree topology and corresponding set of branch lengths, the ML estimates of standard evolutionary model parameters are asymptotically efficient, in the sense that their joint distribution is asymptotically normal with the variance–covariance matrix given by the inverse of the Fisher information matrix. We propose a new estimate of this conditional variance based on estimation of the expected information using a Monte Carlo sampling (MCS) method. Simulations are used to compare this conditional variance estimate to the standard technique of using the observed information under a variety of experimental conditions. In the case in which one wishes to estimate simultaneously the tree and parameters, we provide a bootstrapping approach that can be used in conjunction with the MCS method to estimate the unconditional standard error. The methods developed are applied to a real data set consisting of 30 papillomavirus sequences. This overall method is easily incorporated into standard bootstrapping procedures to allow for proper variance estimation.     

Structured modeling of a microbial system: I. A theoretical study of lactic acid fermentation

Most fermentation models presented in the literature are unstructured, i.e., the biomass composition is assumed constant during all operating conditions. These models are unable to simulate experiments carried out at widely different operating conditions. It is therefore interesting to examine simple structured models where knowledge of the cell physiology is taken into account in the modeling phase. In this article, a simple structured model is presented. The model is based on experimental work with the lactic acid bacteria Streptococcus cremoris, but due to the similarities in basic metabolism for many microorganisms it is applicable also for other fermentation system. The basic assumption in the model is that the biomass can be divided into two parts (compartments)-an active part and a mainly inactive structural part. The size of the active part has a pivotal role in the model.     

A Biophysical Model for Integration of Electrical, Osmotic, and pH Regulation in the Human Bronchial Epithelium

A dynamical biophysical model for the functioning of an epithelium is presented. This model integrates the electrical and osmotic behaviors of the epithelium, taking into account intracellular conditions. The specific tissue modeled is the human bronchial epithelium, which is of particular interest, as it is the location of the most common lethal symptoms of cystic fibrosis. The model is implemented in a modular form to facilitate future application of the code to other epithelial tissue by inputting different transporters, channels, and geometric parameters. The model includes pH regulation as an integral component of overall regulation of epithelial function, through the interdependence of pH, bicarbonate concentration, and current. The procedures for specification, the validation of the model, and parametric studies are presented using available experimental data of cultured human bronchial epithelium. Parametric studies are performed to elucidate a), the contribution of basolateral chloride channels to the short-circuit current functional form, and b), the role that regulation of basolateral potassium conductance plays in epithelial function.     

Parameterization for In-Silico Modeling of Ion Channel Interactions with Drugs

Since the first Hodgkin and Huxley ion channel model was described in the 1950s, there has been an explosion in mathematical models to describe ion channel function. As experimental data has become richer, models have concomitantly been improved to better represent ion channel kinetic processes, although these improvements have generally resulted in more model complexity and an increase in the number of parameters necessary to populate the models. Models have also been developed to explicitly model drug interactions with ion channels. Recent models of drug-channel interactions account for the discrete kinetics of drug interaction with distinct ion channel state conformations, as it has become clear that such interactions underlie complex emergent kinetics such as use-dependent block. Here, we describe an approach for developing a model for ion channel drug interactions. The method describes the process of extracting rate constants from experimental electrophysiological function data to use as initial conditions for the model parameters. We then describe implementation of a parameter optimization method to refine the model rate constants describing ion channel drug kinetics. The algorithm takes advantage of readily available parallel computing tools to speed up the optimization. Finally, we describe some potential applications of the platform including the potential for gaining fundamental mechanistic insights into ion channel function and applications to in silico drug screening and development.     

A joint frailty model for survival and gap times between recurrent events

Therapy for patients with a recurrent disease focuses on delaying disease recurrence and prolonging survival. A common analysis approach for such data is to estimate the distribution of disease-free survival, that is, the time to the first disease recurrence or death, whichever happens first. However, treating death similarly as disease recurrence may give misleading results. Also considering only the first recurrence and ignoring subsequent ones can result in loss of statistical power. We use a joint frailty model to simultaneously analyze disease recurrences and survival. Separate parameters for disease recurrence and survival are used in the joint model to distinguish treatment effects on these two types of events. The correlation between disease recurrences and survival is taken into account by a shared frailty. The effect of disease recurrence on survival can also be estimated by this model. The EM algorithm is used to fit the model, with Markov chain Monte Carlo simulations in the E-steps. The method is evaluated by simulation studies and illustrated through a study of patients with heart failure. Sensitivity analysis for the parametric assumption of the frailty distribution is assessed by simulations.     

Effect of joint laxity on polyethylene wear in total knee replacement

Experimental simulator studies are frequently performed to evaluate wear behavior in total knee replacement. It is vital that the simulation conditions match the physiological situation as closely as possible. To date, few experimental wear studies have examined the effects of joint laxity on wear and joint kinematics and the absence of the anterior cruciate ligament has not been sufficiently taken into account in simulator wear studies.The aim of this study was to investigate different ligament and soft tissue models with respect to wear and kinematics.A virtual soft tissue control system was used to simulate different motion restraints in a force-controlled knee wear simulator.The application of more realistic and sophisticated ligament models that considered the absence of anterior cruciate ligament lead to a significant increase in polyethylene wear (p=0.02) and joint kinematics (p<0.01). We recommend the use of more complex ligament models to appropriately simulate the function of the human knee joint and to evaluate the wear behavior of total knee replacements. A feasible simulation model is presented.     

A random variance model for detection of differential gene expression in small microarray experiments

MOTIVATION: Microarray techniques provide a valuable way of characterizing the molecular nature of disease. Unfortunately expense and limited specimen availability often lead to studies with small sample sizes. This makes accurate estimation of variability difficult, since variance estimates made on a gene by gene basis will have few degrees of freedom, and the assumption that all genes share equal variance is unlikely to be true. RESULTS: We propose a model by which the within gene variances are drawn from an inverse gamma distribution, whose parameters are estimated across all genes. This results in a test statistic that is a minor variation of those used in standard linear models. We demonstrate that the model assumptions are valid on experimental data, and that the model has more power than standard tests to pick up large changes in expression, while not increasing the rate of false positives. AVAILABILITY: This method is incorporated into BRB-ArrayTools version 3.0 (http://linus.nci.nih.gov/BRB-ArrayTools.html). SUPPLEMENTARY MATERIAL: ftp://linus.nci.nih.gov/pub/techreport/RVM_supplement.pdf     

Epithelial transport parameters: an analysis of experimental strategies

A set of experiments was simulated on a computer version of the Koefoed-Johnsen & Ussing model for high-resistance epithelia. The results obtained were analysed according to procedures commonly applied to the analyses of experimental data and interpreted in terms of the model parameters. Although the computer model encodes a stoichiometry of 3:2 for Na-K exchange through the Na pump, the simulation of published experimental procedures yields different figures in almost every case. We show that ENa as originally defined by Ussing & Zerahn (Acta physiol. scand. 23, 110-127 (1951)) and as obtained from flux-ratio experiments has different values under different experimental conditions with unchanged system parameters and that it is distinct from ENa measured by other methods. We also show that unless the pump is saturated with internal Na an increase in the rate of pumping cannot cause a substantial increase in the rate of transepithelial Na transport.     

A comparison between a new model and current models for estimating trunk segment inertial parameters

Modeling of the body segments to estimate segment inertial parameters is required in the kinetic analysis of human motion. A new geometric model for the trunk has been developed that uses various cross-sectional shapes to estimate segment volume and adopts a non-uniform density function that is gender-specific. The goal of this study was to test the accuracy of the new model for estimating the trunk's inertial parameters by comparing it to the more current models used in biomechanical research. Trunk inertial parameters estimated from dual X-ray absorptiometry (DXA) were used as the standard. Twenty-five female and 24 male college-aged participants were recruited for the study. Comparisons of the new model to the accepted models were accomplished by determining the error between the models’ trunk inertial estimates and that from DXA. Results showed that the new model was more accurate across all inertial estimates than the other models. The new model had errors within 6.0% for both genders, whereas the other models had higher average errors ranging from 10% to over 50% and were much more inconsistent between the genders. In addition, there was little consistency in the level of accuracy for the other models when estimating the different inertial parameters. These results suggest that the new model provides more accurate and consistent trunk inertial estimates than the other models for both female and male college-aged individuals. However, similar studies need to be performed using other populations, such as elderly or individuals from a distinct morphology (e.g. obese). In addition, the effect of using different models on the outcome of kinetic parameters, such as joint moments and forces needs to be assessed.     

Towards Predictive Computational Models of Oncolytic Virus Therapy: Basis for Experimental Validation and Model Selection

Oncolytic viruses are viruses that specifically infect cancer cells and kill them, while leaving healthy cells largely intact. Their ability to spread through the tumor makes them an attractive therapy approach. While promising results have been observed in clinical trials, solid success remains elusive since we lack understanding of the basic principles that govern the dynamical interactions between the virus and the cancer. In this respect, computational models can help experimental research at optimizing treatment regimes. Although preliminary mathematical work has been performed, this suffers from the fact that individual models are largely arbitrary and based on biologically uncertain assumptions. Here, we present a general framework to study the dynamics of oncolytic viruses that is independent of uncertain and arbitrary mathematical formulations. We find two categories of dynamics, depending on the assumptions about spatial constraints that govern that spread of the virus from cell to cell. If infected cells are mixed among uninfected cells, there exists a viral replication rate threshold beyond which tumor control is the only outcome. On the other hand, if infected cells are clustered together (e.g. in a solid tumor), then we observe more complicated dynamics in which the outcome of therapy might go either way, depending on the initial number of cells and viruses. We fit our models to previously published experimental data and discuss aspects of model validation, selection, and experimental design. This framework can be used as a basis for model selection and validation in the context of future, more detailed experimental studies. It can further serve as the basis for future, more complex models that take into account other clinically relevant factors such as immune responses.     

An Algorithm for Calculation of n-Dimensional Nonlinear Dependencies – Behaviour Tests on Rats

The evaluation of the data obtained during the behaviour tests always leads to the problem of multiple correlation, very often with non-linear dependencies on the target. All mathematical and statistical procedures that have been used so far are based on the assumption of an equation for the desired correlation for which parameters and related statistical equivalents are determined eventually. The MODAK system applied here (MODAK = algorithms of modelling for the calculation of multi-dimensional non-linear mathematical models) breaks down a complex correlation into individual dependencies in a mathematical and statistical way and selects suitable equations for each of them independently and determines the corresponding parameters. The numerical example evaluates data of behaviour tests on rats. First results obtained on the correlations of various behaviour tests indicate both the possibility of selecting suitable tests independent of each other and a better interpretation of the observed patterns of behaviour taking into account the interrelations between the tests. In addition, MODAK is a method which can be applied as a matter of course in a general way to all cases which call for the reduction and analysis of data occurring in process and system analysis and in the evaluation of test results requiring statistical modelling. So far, MODAK applications range from engineering sciences to medicine.     

A parametric copula model for analysis of familial binary data

Modeling the joint distribution of a binary trait (disease) within families is a tedious challenge, owing to the lack of a general statistical model with desirable properties such as the multivariate Gaussian model for a quantitative trait. Models have been proposed that either assume the existence of an underlying liability variable, the reality of which cannot be checked, or provide estimates of aggregation parameters that are dependent on the ordering of family members and on family size. We describe how a class of copula models for the analysis of exchangeable categorical data can be incorporated into a familial framework. In this class of models, the joint distribution of binary outcomes is characterized by a function of the given marginals. This function, referred to as a "copula," depends on an aggregation parameter that is weakly dependent on the marginal distributions. We propose to decompose a nuclear family into two sets of equicorrelated data (parents and offspring), each of which is characterized by an aggregation parameter (alphaFM and alphaSS, respectively). The marginal probabilities are modeled through a logistic representation. The advantage of this model is that it provides estimates of the aggregation parameters that are independent of family size and does not require any arbitrary ordering of sibs. It can be incorporated easily into segregation or combined segregation-linkage analysis and does not require extensive computer time. As an illustration, we applied this model to a combined segregation-linkage analysis of levels of plasma angiotensin I-converting enzyme (ACE) dichotomized into two classes according to the median. The conclusions of this analysis were very similar to those we had reported in an earlier familial analysis of quantitative ACE levels.     

A task-specific validation of homogeneous non-linear optimisation approaches

In biomechanics, musculoskeletal models are typically redundant. This situation is referred to as the distribution problem. Often, static, non-linear optimisation methods of the form “min: φ(f) subject to mechanical and muscular constraints” have been used to extract a unique set of muscle forces. Here, we present a method for validating this class of non-linear optimisation approaches where the homogeneous cost function, φ(f), is used to solve the distribution problem. We show that the predicted muscle forces for different loading conditions are scaled versions of each other if the joint loading conditions are just scaled versions. Therefore, we can calculate the theoretical muscle forces for different experimental conditions based on the measured muscle forces and joint loadings taken from one experimental condition and assuming that all input into the optimisation (e.g., moment arms, muscle attachment sites, size, fibre type distribution) and the optimisation approach are perfectly correct. Thus predictions of muscle force for other experimental conditions are accurate if the optimisation approach is appropriate, independent of the musculoskeletal geometry and other input required for the optimisation procedure. By comparing the muscle forces predicted in this way to the actual muscle forces obtained experimentally, we conclude that convex homogeneous non-linear optimisation approaches cannot predict individual muscle forces properly, as force-sharing among synergistic muscles obtained experimentally are not just scaled versions of joint loading, not even in a first approximation.     

Multi-locus Analysis of Genomic Time Series Data from Experimental Evolution

Genomic time series data generated by evolve-and-resequence (E&R) experiments offer a powerful window into the mechanisms that drive evolution. However, standard population genetic inference procedures do not account for sampling serially over time, and new methods are needed to make full use of modern experimental evolution data. To address this problem, we develop a Gaussian process approximation to the multi-locus Wright-Fisher process with selection over a time course of tens of generations. The mean and covariance structure of the Gaussian process are obtained by computing the corresponding moments in discrete-time Wright-Fisher models conditioned on the presence of a linked selected site. This enables our method to account for the effects of linkage and selection, both along the genome and across sampled time points, in an approximate but principled manner. We first use simulated data to demonstrate the power of our method to correctly detect, locate and estimate the fitness of a selected allele from among several linked sites. We study how this power changes for different values of selection strength, initial haplotypic diversity, population size, sampling frequency, experimental duration, number of replicates, and sequencing coverage depth. In addition to providing quantitative estimates of selection parameters from experimental evolution data, our model can be used by practitioners to design E&R experiments with requisite power. We also explore how our likelihood-based approach can be used to infer other model parameters, including effective population size and recombination rate. Then, we apply our method to analyze genome-wide data from a real E&R experiment designed to study the adaptation of D. melanogaster to a new laboratory environment with alternating cold and hot temperatures.     

One-dimensional and three-dimensional models of cerebrovascular flow

The Circle of Willis is a ring-like structure of blood vessels found beneath the hypothalamus at the base of the brain. Its main function is to distribute oxygen-rich arterial blood to the cerebral mass. One-dimensional (1D) and three-dimensional (3D) computational fluid dynamics (CFD) models of the Circle of Willis have been created to provide a simulation tool which can potentially be used to identify at-risk cerebral arterial geometries and conditions and replicate clinical scenarios, such as occlusions in afferent arteries and absent circulus vessels. Both models capture cerebral haemodynamic autoregulation using a proportional-integral (PI) controller to modify efferent artery resistances to maintain optimal efferent flow rates for a given circle geometry and afferent blood pressure. The models can be used to identify at-risk cerebral arterial geometries and conditions prior to surgery or other clinical procedures. The 1D model is particularly relevant in this instance, with its fast solution time suitable for real-time clinical decisions. Results show the excellent correlation between models for the transient efferent flux profile. The assumption of strictly Poiseuille flow in the 1D model allows more flow through the geometrically extreme communicating arteries than the 3D model. This discrepancy was overcome by increasing the resistance to flow in the anterior communicating artery in the 1D model to better match the resistance seen in the 3D results.