Development of a continuum mechanics model of passive skeletal muscle

Skeletal muscle force evaluation is difficult to implement in a clinical setting. Muscle force is typically assessed through either manual muscle testing, isokinetic/isometric dynamometry, or electromyography (EMG). Manual muscle testing is a subjective evaluation of a patient’s ability to move voluntarily against gravity and to resist force applied by an examiner. Muscle testing using dynamometers adds accuracy by quantifying functional mechanical output of a limb. However, like manual muscle testing, dynamometry only provides estimates of the joint moment. EMG quantifies neuromuscular activation signals of individual muscles, and is used to infer muscle function. Despite the abundance of work performed to determine the degree to which EMG signals and muscle forces are related, the basic problem remains that EMG cannot provide a quantitative measurement of muscle force. Intramuscular pressure (IMP), the pressure applied by muscle fibers on interstitial fluid, has been considered as a correlate for muscle force. Numerous studies have shown that an approximately linear relationship exists between IMP and muscle force. A microsensor has recently been developed that is accurate, biocompatible, and appropriately sized for clinical use. While muscle force and pressure have been shown to be correlates, IMP has been shown to be non-uniform within the muscle. As it would not be practicable to experimentally evaluate how IMP is distributed, computational modeling may provide the means to fully evaluate IMP generation in muscles of various shapes and operating conditions. The work presented in this dissertation focuses on the development and validation of computational models of passive skeletal muscle and the evaluation of their performance for prediction of IMP. A transversly isotropic, hyperelastic, and nearly incompressible model will be evaluated along with a poroelastic model.

Insulin-like growth factor I improves aspects of mycophenolate mofetil-impaired anastomotic healing in an experimental model

BACKGROUND: Patients taking immunosuppressants after transplantation may require intestinal surgery. Mycophenolate mofetil (MMF) has been found to impair the healing of colonic anastomoses in rats. This study examined whether insulin-like growth factor (IGF) I prevents MMF impairment of anastomotic healing. METHODS: Sixty-three rats were divided into three groups (MMF, MMF/IGF and control). Animals underwent a sigmoid colon anastomosis with a 6/0 suture, and were killed on days 2, 4 and 6 after surgery. Investigations included bursting pressure measurement, morphometric analysis, and assessment of mucosal proliferation by 5-bromo-2'-deoxyuridine and Ki67 immunohistochemistry of the anastomoses. RESULTS: The leak rate was three of 21, one of 20 and two of 20 in the MMF, MMF/IGF-I and control groups respectively. Anastomotic bursting pressures were significantly lower in the MMF group than in the control group on days 2 and 4, but there was no significant difference by day 6. Values in the MMF/IGF-I and control groups were similar. Colonic crypt depth was significantly reduced in MMF-treated animals on days 2 and 4, but this impairment was attenuated by IGF-I on day 4. Similarly, IGF-I reduced the negative impact of MMF on mucosal proliferation on days 2 and 6. CONCLUSION: Exogenous IGF-I improves some aspects of MMF-impaired anastomotic healing.

Assessments and Computational Simulations in Support of a Time-Varying Mass Flow Rate Measurement Technique for Pulsatile Gas Flow

This thesis covers the correction, and verification, development, and implementation of a computational fluid dynamics (CFD) model for an orifice plate meter. Past results were corrected and further expanded on with compressibility effects of acoustic waves being taken into account. One dynamic pressure difference transducer measures the time-varying differential pressure across the orifice meter. A dynamic absolute pressure measurement is also taken at the inlet of the orifice meter, along with a suitable temperature measurement of the mean flow gas. Together these three measurements allow for an incompressible CFD simulation (using a well-tested and robust model) for the cross-section independent time-varying mass flow rate through the orifice meter. The mean value of this incompressible mass flow rate is then corrected to match the mean of the measured flow rate( obtained from a Coriolis meter located up stream of the orifice meter). Even with the mean and compressibility corrections, significant differences in the measured mass flow rates at two orifice meters in a common flow stream were observed. This means that the compressibility effects associated with pulsatile gas flows is significant in the measurement of the time-varying mass flow rate. Future work (with the approach and initial runs covered here) will provide an indirect verification of the reported mass flow rate measurements.

Double-lung transplantation in the rat: an acute, syngeneic in situ model

BACKGROUND. The high rate of reperfusion injury in clinical lung transplantation mandates significant improvements in lung preservation. Innovations should be validated using standardized and low-cost experimental models. METHODS. The model introduced here is analyzed by comparing global lung function after varying ischemic times (2, 4, 8, 16, and 24 hours). A rat double-lung block is flush-perfused, and the main pulmonary artery and left atrium are connected to the left pulmonary artery and vein of a syngeneic recipient using a T-shaped stent. With pressure side ports and incorporated flow crystals, measurement of vascular resistance and graft oxygenation can be performed. The transplant is ventilated separately, and compliance and resistance are determined. RESULTS. The increase in the ischemic interval from 2 to 24 hours caused an increase in the alveolar arterial oxygen difference from 220 +/- 20 to 600 +/- 34 mm Hg, pulmonary vascular resistance from 198 +/- 76 to 638 +/- 212 mm Hg.mL-1.min-1, and resistance to airflow from 274 +/- 50 to 712 +/- 30 cm H2O/L H2O, and a decrease in pulmonary compliance from 0.4 +/- 0.05 to 0.12 +/- 0.06 mL/cm H2O. CONCLUSIONS. This in situ, syngeneic rat lung transplantation model offers an alternative to large animal models for verification of lung preservation solutions and for modification of donor or recipient treatment regimens.

Multivariate modelling of responses to conditional items: New possibilities for latent class analysis

Questionnaire data may contain missing values because certain questions do not apply to all respondents. For instance, questions addressing particular attributes of a symptom, such as frequency, triggers or seasonality, are only applicable to those who have experienced the symptom, while for those who have not, responses to these items will be missing. This missing information does not fall into the category 'missing by design', rather the features of interest do not exist and cannot be measured regardless of survey design. Analysis of responses to such conditional items is therefore typically restricted to the subpopulation in which they apply. This article is concerned with joint multivariate modelling of responses to both unconditional and conditional items without restricting the analysis to this subpopulation. Such an approach is of interest when the distributions of both types of responses are thought to be determined by common parameters affecting the whole population. By integrating the conditional item structure into the model, inference can be based both on unconditional data from the entire population and on conditional data from subjects for whom they exist. This approach opens new possibilities for multivariate analysis of such data. We apply this approach to latent class modelling and provide an example using data on respiratory symptoms (wheeze and cough) in children. Conditional data structures such as that considered here are common in medical research settings and, although our focus is on latent class models, the approach can be applied to other multivariate models.

Earnings Management and Measurement Error

In this paper I consider the impact of a noisy indicator regarding a manager’s manipulative behavior on optimal effort incentives and the extent of earnings management. The analysis in this paper extends a twotask, single performance measure LEN model by including a binary random variable. I show that contracting on the noisy indicator variable is not always useful. More specifically, the principal uses the indicator variable to prevent earnings management only under conditions where manipulative behavior is not excessive. Thus, under conditions of excessive earnings management, accounting adjustments that yield a more congruent overall performance measure can be more effective than an appraisal of the existence of earnings management to mitigate the earnings management problem.

Meta-model for calculating the mean energy demand of automated storage and retrieval systems

Energy efficiency has become an important research topic in intralogistics. Especially in this field the focus is placed on automated storage and retrieval systems (AS/RS) utilizing stacker cranes as these systems are widespread and consume a significant portion of the total energy demand of intralogistical systems. Numerical simulation models were developed to calculate the energy demand rather precisely for discrete single and dual command cycles. Unfortunately these simulation models are not suitable to perform fast calculations to determine a mean energy demand value of a complete storage aisle. For this purpose analytical approaches would be more convenient but until now analytical approaches only deliver results for certain configurations. In particular, for commonly used stacker cranes equipped with an intermediate circuit connection within their drive configuration there is no analytical approach available to calculate the mean energy demand. This article should address this research gap and present a calculation approach which enables planners to quickly calculate the energy demand of these systems.

Accounting for baseline differences and measurement error in the analysis of change over time

If change over time is compared in several groups, it is important to take into account baseline values so that the comparison is carried out under the same preconditions. As the observed baseline measurements are distorted by measurement error, it may not be sufficient to include them as covariate. By fitting a longitudinal mixed-effects model to all data including the baseline observations and subsequently calculating the expected change conditional on the underlying baseline value, a solution to this problem has been provided recently so that groups with the same baseline characteristics can be compared. In this article, we present an extended approach where a broader set of models can be used. Specifically, it is possible to include any desired set of interactions between the time variable and the other covariates, and also, time-dependent covariates can be included. Additionally, we extend the method to adjust for baseline measurement error of other time-varying covariates. We apply the methodology to data from the Swiss HIV Cohort Study to address the question if a joint infection with HIV-1 and hepatitis C virus leads to a slower increase of CD4 lymphocyte counts over time after the start of antiretroviral therapy.

ANOVA kernels and RKHS of zero mean functions for model-based sensitivity analysis

Given a reproducing kernel Hilbert space (H,〈.,.〉)(H,〈.,.〉) of real-valued functions and a suitable measure μμ over the source space D⊂RD⊂R, we decompose HH as the sum of a subspace of centered functions for μμ and its orthogonal in HH. This decomposition leads to a special case of ANOVA kernels, for which the functional ANOVA representation of the best predictor can be elegantly derived, either in an interpolation or regularization framework. The proposed kernels appear to be particularly convenient for analyzing the effect of each (group of) variable(s) and computing sensitivity indices without recursivity.

Measurement of structural anisotropy in femoral trabecular bone using clinical-resolution CT images

Discrepancies in finite-element model predictions of bone strength may be attributed to the simplified modeling of bone as an isotropic structure due to the resolution limitations of clinical-level Computed Tomography (CT) data. The aim of this study is to calculate the preferential orientations of bone (the principal directions) and the extent to which bone is deposited more in one direction compared to another (degree of anisotropy). Using 100 femoral trabecular samples, the principal directions and degree of anisotropy were calculated with a Gradient Structure Tensor (GST) and a Sobel Structure Tensor (SST) using clinical-level CT. The results were compared against those calculated with the gold standard Mean-Intercept-Length (MIL) fabric tensor using micro-CT. There was no significant difference between the GST and SST in the calculation of the main principal direction (median error=28°), and the error was inversely correlated to the degree of transverse isotropy (r=−0.34, p<0.01). The degree of anisotropy measured using the structure tensors was weakly correlated with the MIL-based measurements (r=0.2, p<0.001). Combining the principal directions with the degree of anisotropy resulted in a significant increase in the correlation of the tensor distributions (r=0.79, p<0.001). Both structure tensors were robust against simulated noise, kernel sizes, and bone volume fraction. We recommend the use of the GST because of its computational efficiency and ease of implementation. This methodology has the promise to predict the structural anisotropy of bone in areas with a high degree of anisotropy, and may improve the in vivo characterization of bone.

AMPA RECEPTOR ALLOSTERISM: MEASUREMENT OF THE CONFORMATIONAL CHANGES IN THE LIGAND BINDING DOMAIN OF A FUNCTIONAL RECEPTOR

Ionotropic glutamate receptors are important excitatory neurotransmitter receptors in the mammalian central nervous system that have been implicated in a number of neuropathologies such as epilepsy, ischemia, and amyotrophic lateral sclerosis. Glutamate binding to an extracellular ligand binding domain initiates a series of structural changes that leads to the formation of a cation selective transmembrane channel, which consequently closes due to desensitization of the receptor. The crystal structures of the AMPA subtype of the glutamate receptor have been particularly useful in providing initial insight into the conformational changes in the ligand binding domain; however, these structures are limited by crystallographic constraint. To gain a clear picture of how agonist binding is coupled to channel activation and desensitization, it is essential to study changes in the ligand binding domain in a dynamic, physiological state. In this dissertation, a technique called Luminescence Resonance Energy Transfer was used to determine the conformational changes associated with activation and desensitization in a functional AMPA receptor (ÄN*-AMPA) that contains the ligand binding domain and transmembrane segments; ÄN*-AMPA has been modified such that fluorophores can be introduced at specific sites to serve as a readout of cleft closure or to establish intersubunit distances. Previous structural studies of cleft closure of the isolated ligand binding domain in conjunction with functional studies of the full receptor suggest that extent of cleft closure correlates with extent of activation. Here, LRET has been used to show that a similar relationship between cleft closure and activation is observed in the “full length” receptor showing that the isolated ligand binding domain is a good model of the domain in the full length receptor for changes within a subunit. Similar LRET investigations were used to study intersubunit distances specifically to probe conformational changes between subunits within a dimer in the tetrameric receptor. These studies show that the dimer interface is coupled in the open state, and decoupled in the desensitized state, similar to the isolated ligand binding domain crystal structure studies. However, we show that the apo state dimer interface is not pre-formed as in the crystal structure, hence suggesting a mechanism for functional transitions within the receptor based on LRET distances obtained.

PREDICTIVE VARIABLES FOR WEIGHT CHANGE 2 YEARS FOLLOWING PARTICIPATION IN A WEIGHT LOSS INTERVENTION PROGRAM

A historical prospective study was designed to assess the man weight status of subjects who participated in a behavioral weight reduction program in 1983 and to determine whether there was an association between the dependent variable weight change and any of 31 independent variables after a 2 year follow-up period. Data was obtained by abstracting the subjects records and from a follow-up questionnaire administered 2 years following program participation. Five hundred nine subjects (386 females and 123 males) of 1460 subjects who participated in the program, completed and returned the questionnaire. Results showed that mean weight was significantly different (p < 0.001) between the measurement at baseline and after a 2 year follow-up period. The mean weight loss of the group was 5.8 pounds, 10.7 pounds for males and 4.2 pounds for females after a 2 year follow-up period. A total of 63.9% of the group, 69.9% of males and 61.9% of females were still below their initial weight after the 2 year follow-up period. Sixteen of the 31 variables assessed utilizing bivariate analyses were found to be significantly (p (LESSTHEQ) 0.05) associated with weight change after a 2 year follow-up period. These variables were then entered into a multivariate linear regression model. A total of 37.9% of the variance of the dependent variable, weight change, was accounted for by all 16 variables. Eight of these variables were found to be significantly (p (LESSTHEQ) 0.05) predictive of weight change in the stepwise multivariate process accounting for 37.1% of the variance. These variables included: Two baseline variables (percent over ideal body weight at enrollment and occupation) and six follow-up variables (feeling in control of eating habits, percent of body weight lost during treatment, frequency of weight measurement, physical activity, eating in response to emotions, and number of pounds of weight gain needed to resume a diet). It was concluded that a greater amount of emphasis should be placed on the six follow-up variables by clinicians involved in the treatment of obesity, and by the subjects themselves to enhance their chances of success at long-term weight loss. ^

Reactive oxygen intermediates contribute to necrotic and apoptotic neuronal injury in an infant rat model of bacterial meningitis due to group B streptococci.

Reactive oxygen intermediates (ROI) contribute to neuronal injury in cerebral ischemia and trauma. In this study we explored the role of ROI in bacterial meningitis. Meningitis caused by group B streptococci in infant rats led to two distinct forms of neuronal injury, areas of necrosis in the cortex and neuronal loss in the dentate gyrus of the hippocampus, the latter showing evidence for apoptosis. Staining of brain sections with diaminobenzidine after perfusion with manganese buffer and measurement of lipid peroxidation products in brain homogenates both provided evidence that meningitis led to the generation of ROI. Treatment with the radical scavenger alpha-phenyl-tert-butyl nitrone (PBN) (100 mg/kg q8h i.p.) beginning at the time of infection completely abolished ROI detection and the increase in lipidperoxidation. Cerebral cortical perfusion was reduced in animals with meningitis to 37.5+/-21.0% of uninfected controls (P < 0.05), and PBN restored cortical perfusion to 72.0+/-8.1% of controls (P < 0.05 vs meningitis). PBN also completely prevented neuronal injury in the cortex and hippocampus, when started at the time of infection (P < 0.02), and significantly reduced both forms of injury, when started 18 h after infection together with antibiotics (P < 0.004 for cortex and P < 0.001 for hippocampus). These data indicate that the generation of ROI is a major contributor to cerebral ischemia and necrotic and apoptotic neuronal injury in this model of neonatal meningitis.

Measurement of Top Quark Polarization in Top-Antitop Events from Proton-Proton Collisions at s√=7 TeV Using the ATLAS Detector

his Letter presents measurements of the polarization of the top quark in top-antitop quark pair events, using 4.7  fb−1 of proton-proton collision data recorded with the ATLAS detector at the Large Hadron Collider at s√=7  TeV. Final states containing one or two isolated leptons (electrons or muons) and jets are considered. Two measurements of αℓP, the product of the leptonic spin-analyzing power and the top quark polarization, are performed assuming that the polarization is introduced by either a CP conserving or a maximally CP violating production process. The measurements obtained, αℓPCPC=−0.035±0.014(stat)±0.037(syst) and αℓPCPV=0.020±0.016(stat)+0.013−0.017(syst), are in good agreement with the standard model prediction of negligible top quark polarization.

Measurement of the distributions of event-by-event flow harmonics in lead-lead collisions at = 2.76 TeV with the ATLAS detector at the LHC

The distributions of event-by-event harmonic flow coefficients v_n for n=2-4 are measured in sqrt(s_NN)=2.76 TeV Pb+Pb collisions using the ATLAS detector at the LHC. The measurements are performed using charged particles with transverse momentum pT> 0.5 GeV and in the pseudorapidity range |eta|<2.5 in a dataset of approximately 7 ub^-1 recorded in 2010. The shapes of the v_n distributions are described by a two-dimensional Gaussian function for the underlying flow vector in central collisions for v_2 and over most of the measured centrality range for v_3 and v_4. Significant deviations from this function are observed for v_2 in mid-central and peripheral collisions, and a small deviation is observed for v_3 in mid-central collisions. It is shown that the commonly used multi-particle cumulants are insensitive to the deviations for v_2. The v_n distributions are also measured independently for charged particles with 0.51 GeV. When these distributions are rescaled to the same mean values, the adjusted shapes are found to be nearly the same for these two pT ranges. The v_n distributions are compared with the eccentricity distributions from two models for the initial collision geometry: a Glauber model and a model that includes corrections to the initial geometry due to gluon saturation effects. Both models fail to describe the experimental data consistently over most of the measured centrality range.