Proceedings of Computational Methods for Coupled Problems in Science and Engineering

Computational Methods for Coupled Problems in Science and Engineering

Methods for prediction of peptide binding to MHC molecules: A comparative study

Background: A variety of methods for prediction of peptide binding to major histocompatibility complex (MHC) have been proposed. These methods are based on binding motifs, binding matrices, hidden Markov models (HMM), or artificial neural networks (ANN). There has been little prior work on the comparative analysis of these methods. Materials and Methods: We performed a comparison of the performance of six methods applied to the prediction of two human MHC class I molecules, including binding matrices and motifs, ANNs, and HMMs. Results: The selection of the optimal prediction method depends on the amount of available data (the number of peptides of known binding affinity to the MHC molecule of interest), the biases in the data set and the intended purpose of the prediction (screening of a single protein versus mass screening). When little or no peptide data are available, binding motifs are the most useful alternative to random guessing or use of a complete overlapping set of peptides for selection of candidate binders. As the number of known peptide binders increases, binding matrices and HMM become more useful predictors. ANN and HMM are the predictive methods of choice for MHC alleles with more than 100 known binding peptides. Conclusion: The ability of bioinformatic methods to reliably predict MHC binding peptides, and thereby potential T-cell epitopes, has major implications for clinical immunology, particularly in the area of vaccine design.

Computational methods for prediction of T-cell epitopes - a framework for modelling, testing, and applications

Computational models complement laboratory experimentation for efficient identification of MHC-binding peptides and T-cell epitopes. Methods for prediction of MHC-binding peptides include binding motifs, quantitative matrices, artificial neural networks, hidden Markov models, and molecular modelling. Models derived by these methods have been successfully used for prediction of T-cell epitopes in cancer, autoimmunity, infectious disease, and allergy. For maximum benefit, the use of computer models must be treated as experiments analogous to standard laboratory procedures and performed according to strict standards. This requires careful selection of data for model building, and adequate testing and validation. A range of web-based databases and MHC-binding prediction programs are available. Although some available prediction programs for particular MHC alleles have reasonable accuracy, there is no guarantee that all models produce good quality predictions. In this article, we present and discuss a framework for modelling, testing, and applications of computational methods used in predictions of T-cell epitopes. (C) 2004 Elsevier Inc. All rights reserved.

Timing of rotator cuff activation during shoulder external rotation in throwers with and without symptoms of pain

Study Design: Fine-wire EMG rotator cuff onset time analysis in 2 matched groups of throwers with and without pain. Objective: To identify if there is a difference in the activation patterns of the rotator cuff muscles during a rapid shoulder external rotation task between throwers with and without pain. Background: The coordinated action of the rotator cuff is recognized as essential for glenohumeral joint control in the throwing athlete. Identification of abnormalities occurring in muscle activation patterns for injured athletes is relevant when prescribing rehabilitative exercises. Methods and Measures: Twelve throwers with shoulder pain were compared to a matched group of 11 asymptomatic throwers. Participants were matched for age, height, body mass, and habitual activity. Fine-wire EMG electrodes were inserted into the subscapularis, supraspinatus, and infraspinatus. EMG activity was measured during a reaction time task of rapid shoulder external rotation in a seated position. The timing of onset of EMG activity was analyzed in relation to visualization of a light (reaction time) and to the onset of infraspinatus activity (relative latency). Results: In the group with shoulder pain, the onset of subscapularis activity was found to be significantly delayed (reaction time, P = .0018; relative latency, P = .0005) from the onset of infraspinatus activity when compared to the control group. Conclusions: The presence of shoulder pain in these athletes was associated with a difference in the onset of subscapularis EMG activity during a rapid shoulder external rotation movement. This was an initial step in the understanding of the joint protection mechanisms of the glenohumeral joint and the problems that occur in throwers. This information may assist in providing future guidelines for more effective rehabilitation and prevention strategies for this condition.

Birth of Koalas Phascolarctos cinereus at Lone Pine Koala Sanctuary following artificial insemination

This paper documents the successful development of an artificial insemination (AI) programme for the Koala Phascolurctos cinereus. The protocols for trials involving two methods to induce ovulation and two insemination techniques are described. In Trial 1, interrupted coitus using a 'teaser'♂ successfully induced ovulation in nine Koalas. Five ♀♀ were inseminated while conscious using a modified 'foley catheter' (Cook insemination catheter) resulting in the births of two offspring. The other four ♀♀ were anaesthetized and inseminated using a technique which allowed visualization of the most cranial portion of the urogenital sinus, where semen was deposited using a 3.5 Fr. 'Tom-cat catheter' (urogen-itoscopic insemination). Three of the four ♀♀ inseminated by this technique produced pouch young. Microsatellite analysis of DNA from the pouch young excluded the teaser ♀♀ as possible sires, confirming that all offspring were sired by donor sperm. In Trial 2, eight ♀♀ were induced to ovulate by injecting them with 250 International Units of human chorionic gonadotrophin (hCG). A luteal phase was confirmed in all eight ♀♀ but only one gave birth following urogenitoscopic insemination. The Koala pouch young in this study are the first of any marsupial to be conceived and born following A1 procedures. Details of the A1 procedures used are presented and the significance of A1 to the conservation biology of P. cinereus discussed.

Numerical methods for strong solutions of stochastic differential equations: an overview

This paper gives a review of recent progress in the design of numerical methods for computing the trajectories (sample paths) of solutions to stochastic differential equations. We give a brief survey of the area focusing on a number of application areas where approximations to strong solutions are important, with a particular focus on computational biology applications, and give the necessary analytical tools for understanding some of the important concepts associated with stochastic processes. We present the stochastic Taylor series expansion as the fundamental mechanism for constructing effective numerical methods, give general results that relate local and global order of convergence and mention the Magnus expansion as a mechanism for designing methods that preserve the underlying structure of the problem. We also present various classes of explicit and implicit methods for strong solutions, based on the underlying structure of the problem. Finally, we discuss implementation issues relating to maintaining the Brownian path, efficient simulation of stochastic integrals and variable-step-size implementations based on various types of control.

Folate: Methods of analysis

Folates and its derivatives occur as polyglutamates in nature. The multiplicity of forms and the generally low levels in foods makes quantitative analysis of folate a difficult task. The assay of folates from foods generally involves three steps: liberation of folates from the cellular matrix; deconjugation from the polyglutamate to the mono and di-glutamate forms; and the detection of the biological activity or chemical concentration of the resulting folates. The detection methods used are the microbiological assay relying on the turbidimetric bacterial growth of Lactobacillus rhamnosus which is by far the most commonly used method; the HPLC and LC/MS techniques and bio-specific procedures. This review attempts to describe the methods along with the merits and demerits of using each of these methods.

Do different surrogate methods detect lateral genetic transfer events of different relative ages?

Non-tree-based ('surrogate') methods have been used to identify instances of lateral genetic transfer in microbial genomes but agreement among predictions of different methods can be poor. It has been proposed that this disagreement arises because different surrogate methods are biased towards the detection of certain types of transfer events. This conjecture is supported by a rigorous phylogenetic analysis of 3776 proteins in Escherichia coli K12 MG1655 to map the ages of transfer events relative to one another.

Quantitative analysis of DNA-protein interactions using double-labeled native gel electrophoresis and fluorescence-based imaging

We have developed a sensitive, non-radioactive method to assess the interaction of transcription factors/DNA-binding proteins with DNA. We have modified the traditional radiolabeled DNA gel mobility shift assay to incorporate a DNA probe end-labeled with a Texas-red fluorophore and a DNA-binding protein tagged with the green fluorescent protein to monitor precisely DNA-protein complexation by native gel electrophoresis. We have applied this method to the DNA-binding proteins telomere release factor-1 and the sex-determining region-Y, demonstrating that the method is sensitive (able to detect 100 fmol of fluorescently labeled DNA), permits direct visualization of both the DNA probe and the DNA-binding protein, and enables quantitative analysis of DNA and protein complexation, and thereby an estimation of the stoichiometry of protein-DNA binding.

Analysis of trigonometric implicit Runge-Kutta methods

Using generalized collocation techniques based on fitting functions that are trigonometric (rather than algebraic as in classical integrators), we develop a new class of multistage, one-step, variable stepsize, and variable coefficients implicit Runge-Kutta methods to solve oscillatory ODE problems. The coefficients of the methods are functions of the frequency and the stepsize. We refer to this class as trigonometric implicit Runge-Kutta (TIRK) methods. They integrate an equation exactly if its solution is a trigonometric polynomial with a known frequency. We characterize the order and A-stability of the methods and establish results similar to that of classical algebraic collocation RK methods. (c) 2006 Elsevier B.V. All rights reserved.