Prediction of MHC class II-binding peptides using an evolutionary algorithm and artificial neural network

Motivation: Prediction methods for identifying binding peptides could minimize the number of peptides required to be synthesized and assayed, and thereby facilitate the identification of potential T-cell epitopes. We developed a bioinformatic method for the prediction of peptide binding to MHC class II molecules. Results: Experimental binding data and expert knowledge of anchor positions and binding motifs were combined with an evolutionary algorithm (EA) and an artificial neural network (ANN): binding data extraction --> peptide alignment --> ANN training and classification. This method, termed PERUN, was implemented for the prediction of peptides that bind to HLA-DR4(B1*0401). The respective positive predictive values of PERUN predictions of high-, moderate-, low- and zero-affinity binder-a were assessed as 0.8, 0.7, 0.5 and 0.8 by cross-validation, and 1.0, 0.8, 0.3 and 0.7 by experimental binding. This illustrates the synergy between experimentation and computer modeling, and its application to the identification of potential immunotheraaeutic peptides.

Hepatic structure-pharmacokinetic relationships: The hepatic disposition and metabolite kinetics of a homologous series of O-acyl derivatives of salicylic acid

1 The hepatic disposition and metabolite kinetics of a homologous series of O-acyl (acetyl, propionyl, butanoyl, pentanoyl, hexanoyl and octanoyl) esters of salicylic acid (C2SA, C3SA, C4SA, C5SA, C6SA and C8SA, respectively) was determined using a single-pass, in-sills rat liver preparation. 2 The hepatic venous outflow profiles for the parent esters and the generated metabolite, salicylic acid (SA) were analysed by HPLC. Non-parametric moments analysis was used to determine the area under the curve (AUC'), mean transit time (MTT) and normalized variance (CV2) for the parent esters and generated SA. 3 Pregenerated SA ([C-14]-salicylic acid) was injected into each liver with the parent ester to determine its distribution characteristics. 4 The overall recovery of ester plus metabolite was 89% of the ester dose injected and independent of the ester carbon number, suggesting that ester extraction was due to hepatic metabolism to salicylic acid. 5 The metabolite AUC' value increased directly with the lipophilicity of the parent ester (from 0.12 for C2SA to 0.95 for C8SA). By contrast, the parent AUC' decreased with the lipophilicity (from 0.85 for C2SA to zero for C8SA). The metabolite MTT value also showed a trend to increase with the lipophilicity of the parent ester (from 15.72 s for C3SA to 61.97 s for C8SA). However, the parent MTT value shows no significant change across the series. 6 The two-compartment dispersion model was used to derive the kinetic parameters for parent ester, pregenerated SA and generated SA. Consequently, these parameters were used to estimate the values of AUG', MITT and CV2 for the parent ester and metabolite. The moments values obtained using the two-compartment dispersion model show similar trends to the corresponding moments values obtained from the outflow profiles using a non-parametric approach. 7 The more lipophilic aspirin analogues are more confined to the portal circulation after oral administration than aspirin due to their more extensive hepatic elimination avoiding systemic prostacyclin inhibition. Given that aspirin's selectivity as an anti-thrombotic agent has been postulated to be due to selective anti-platelet effects in the portal circulation, the more lipophilic and highly extracted analogues are potentially more selective anti-thrombotic agents than aspirin.

Spectroscopic identification of a dinuclear metal centre in manganese(II)-activated aminopeptidase P from Escherichia coli: implications for human prolidase

Electron paramagnetic resonance (EPR) spectra and X-ray absorption (EXAFS and XANES) data have been recorded for the manganese enzyme aminopeptidase P (AMPP, PepP protein) from Escherichia coli. The biological function of the protein, a tetramer of 50-kDa subunits, is the hydrolysis of N-terminal Xaa-Pro peptide bonds. Activity assays confirm that the enzyme is activated by treatment with Mn2+. The EPR spectrum of Mn2+-activated AMPP at liquid-He temperature is characteristic of an exchange-coupled dinuclear Mn(II) site, the Mn-Mn separation calculated from the zero-field splitting D of the quintet state being 3.5 (+/- 0.1) Angstrom. In the X-ray absorption spectrum of Mn2+-activated AMPP at the Mn K edge, the near-edge features are consistent with octahedrally coordinated Mn atoms in oxidation state +2. EXAFS data, limited to k less than or equal to 12 Angstrom(-1) by traces of Fe in the protein, are consistent with a single coordination shell occupied predominantly by O donor atoms at an average Mn-ligand distance of 2.15 Angstrom, but the possibility of a mixture of O and N donor atoms is not excluded. The Mn-Mn interaction at 3.5 Angstrom, is not detected in the EXAFS, probably due to destructive interference from light outer-shell atoms. The biological function, amino acid sequence and metal-ion dependence of E. coli AMPP are closely related to those of human prolidase, an enzyme that specifically cleaves Xaa-Pro dipeptides. Mutations that lead to human prolidase deficiency and clinical symptoms have been identified. Several known inhibitors of prolidase also inhibit AMPP. When these inhibitors are added to Mn2+-activated AMPP, the EPR spectrum and EXAFS remain unchanged. It can be inferred that the inhibitors either do not bind directly to the Mn centres, or substitute for existing Mn ligands without a significant change in donor atoms or coordination geometry. The conclusions from the spectroscopic measurements on AMPP have been verified by, and complement, a recent crystal structure analysis.

Success of surgery for primary aldosteronism judged by residual autonomous aldosterone production

Since February 1996 we have prospectively assessed residual adrenal autonomy by the fludrocortisone suppression test (FST) in 23 patients 3 months after unilateral adrenalectomy for Conn syndrome and in 45 patients after a longer interval. In regard to blood pressure, 36 (53%) patients were cured of hypertension and the remaining 32 (47%) patients had improved hypertension control at the time of their latest postoperative clinical assessment. In regard to the outcome of surgery, patients who achieved normal suppressibility of aldosterone were regarded as cured, and those who had greater suppressibility after surgery were considered improved. Time since surgery for the whole group averaged 26 months. By these biochemical criteria, 42 patients (62%) were cured by surgery, and the rest improved; 16 (76%) of 21 women were cured, and 26 (55%) of 47 men. The women (mean +/- SD age 47 +/- 11 years) were significantly (p < 0.05) younger than the men (52 +/- 9 Sears). Preoperative aldosterone levels before and after FST were similar in the cured and improved groups and fell significantly (p < 0.01) in both groups following surgery. After surgical reduction of autonomous aldosterone production, mean plasma renin activity levels increased sixfold in the cured group and threefold in the improved group. Surgical mortality in this group of 68 patients with Conn syndrome was zero.

Eigenvalues of Casimir operators for gl(m/infinity)

A full set of Casimir operators for the Lie superalgebra gl(m/infinity) is constructed and shown to be well defined in the category O-FS generated by the highest-weight irreducible representations with only a finite number of non-zero weight components. The eigenvalues of these Casimir operators are determined explicitly in terms of the highest weight. Characteristic identities satisfied by certain (infinite) matrices with entries from gl(m/infinity) are also determined.

Fluctuations in heroin purity and the incidence of fatal heroin overdose

In order to determine the role played by heroin purity in fatal heroin overdoses, time series analyses were conducted on the purity of street heroin seizures in south western Sydney and overdose fatalities in that region. A total of 322 heroin samples were analysed in fortnightly periods between February 1993 to January 1995. A total of 61 overdose deaths occurred in the region in the study period. Cross correlation plots revealed a significant correlation of 0.57 at time lag zero between mean purity of heroin samples per fortnight and number of overdose fatalities. Similarly, there was a significant correlation of 0.50 at time lag zero between the highest heroin purity per fortnight and number of overdose fatalities. The correlation between range of heroin purity and number of deaths per fortnight was 0.40. A simultaneous multiple regression on scores adjusted for first order correlation indicated both the mean level of heroin purity and the range of heroin purity were independent predictors of the number of deaths per fortnight. The results indicate that the occurrence of overdose fatalities was moderately associated with both the average heroin purity and the range of heroin purity over the study period. (C) 1999 Elsevier Science Ireland Ltd. All rights reserved.

Bounds on integrals of the Wigner function

The integral of the Wigner function over a subregion of the phase space of a quantum system may be less than zero or greater than one. It is shown that for systems with 1 degree of freedom, the problem of determining the best possible upper and lower bounds on such an integral, over an possible states, reduces to the problem of finding the greatest and least eigenvalues of a Hermitian operator corresponding to the subregion. The problem is solved exactly in the case of an arbitrary elliptical region. These bounds provide checks on experimentally measured quasiprobability distributions.

Experimental examination of solute transport by surface runoff from low-angle slopes

The removal of chemicals in solution, by overland flow from agricultural land has the potential to be a significant source of chemical loss from zero-till and surface mulched farming systems. The objective of this study was to determine the magnitude of solute loss by surface runoff from agricultural systems. Previous experiments have enhanced the understanding of the exchange process, but the initial soil conditions together with the tracer application method in these experiments have meant that in some cases the results have limited applicability to field situations. In this study, two different sets of experiments were carried out to determine the magnitude of solute loss by surface runoff. These experiments entailed the surface application of bromide to (1) field scale plots 18 m long by 2 m wide and (2) repacked soil cores 236 mm in diameter; followed by the application of simulated rainfall in both cases. The most substantial finding of the field experiments was that the quantities of solute in surface runoff varied greatly with soil type and structure (0.07-14.9% of the applied bromide). Also, on some soils, large quantities of tracer were measured in the surface runoff even after several hours of infiltration. The experiments on soil cores showed that soil structure plays an important role in the quantity of chemical that may be transported in the surface runoff. These field results showed that, in certain systems, solute movement by overland flow is an important transport mechanism, which should be considered when budgeting for chemical loss. (C) 2000 Elsevier Science B.V. All rights reserved.

Early diagnosis of lymphedema in postsurgery breast cancer patients

Lymphedema is an accumulation of lymph fluid in the limb resulting from an insufficiency of the lymphatic system. It is commonly associated with surgical or radiotherapy treatment for breast cancer. As with many progressively debilitating disorders, the effectiveness of treatment is significantly improved by earlier intervention. Multiple frequency bioelectrical impedance analysis (MFBIA) previously was shown to provide accurate relative measures of lymphedema in the upper limb in patients after treatment for breast cancer, This presentation reports progress to date on a three-year prospective study to evaluate the efficacy of MFBIA to predict the early onset of lymphedema in breast cancer patients following treatment. Bioelectrical impedance measurements of each upper limb were recorded in a group of healthy control subjects (n = 50) to determine the ratio of extracellular limb-fluid volumes. From this population, the expected normal range of asymmetry (99.7% confidence) between the limbs was determined, Patients undergoing surgery to treat breast cancer were recruited into the study, and MFBIA measurements were recorded presurgery, at one month and three months after surgery, and then at two-month intervals for up to 24 months postsurgery, When patients had an MFBIA measure outside the 99.7% range of the control group, they were referred to their physician for clinical assessment. Results to date: Over 100 patients were recruited into the study over the past two years; at present, 19 have developed lymphedema and, of these, 12 are receiving treatment. In each of these 19 cases, MFBIA predicted the onset of the condition up to four months before it could be clinically diagnosed. The false-negative rate currently is zero, The study will continue to monitor patients over the remaining year to accurately ascertain estimates of specificity and sensitivity of the procedure.

A numerical study of flexural buckling of foliated rock slopes

The occurrence of foliated rock masses is common in mining environment. Methods employing continuum approximation in describing the deformation of such rock masses possess a clear advantage over methods where each rock layer and each inter-layer interface (joint) is explicitly modelled. In devising such a continuum model it is imperative that moment (couple) stresses and internal rotations associated with the bending of the rock layers be properly incorporated in the model formulation. Such an approach will lead to a Cosserat-type theory. In the present model, the behaviour of the intact rock layer is assumed to be linearly elastic and the joints are assumed to be elastic-perfectly plastic. Condition of slip at the interfaces are determined by a Mohr-Coulomb criterion with tension cut off at zero normal stress. The theory is valid for large deformations. The model is incorporated into the finite element program AFENA and validated against an analytical solution of elementary buckling problems of a layered medium under gravity loading. A design chart suitable for assessing the stability of slopes in foliated rock masses against flexural buckling failure has been developed. The design chart is easy to use and provides a quick estimate of critical loading factors for slopes in foliated rock masses. It is shown that the model based on Euler's buckling theory as proposed by Cavers (Rock Mechanics and Rock Engineering 1981; 14:87-104) substantially overestimates the critical heights for a vertical slope and underestimates the same for sub-vertical slopes. Copyright (C) 2001 John Wiley & Sons, Ltd.

Analytical solution to the zero-inertia problem for surge flow phenomena in nonprismatic channels

Surge flow phenomena. e.g.. as a consequence of a dam failure or a flash flood, represent free boundary problems. ne extending computational domain together with the discontinuities involved renders their numerical solution a cumbersome procedure. This contribution proposes an analytical solution to the problem, It is based on the slightly modified zero-inertia (ZI) differential equations for nonprismatic channels and uses exclusively physical parameters. Employing the concept of a momentum-representative cross section of the moving water body together with a specific relationship for describing the cross sectional geometry leads, after considerable mathematical calculus. to the analytical solution. The hydrodynamic analytical model is free of numerical troubles, easy to run, computationally efficient. and fully satisfies the law of volume conservation. In a first test series, the hydrodynamic analytical ZI model compares very favorably with a full hydrodynamic numerical model in respect to published results of surge flow simulations in different types of prismatic channels. In order to extend these considerations to natural rivers, the accuracy of the analytical model in describing an irregular cross section is investigated and tested successfully. A sensitivity and error analysis reveals the important impact of the hydraulic radius on the velocity of the surge, and this underlines the importance of an adequate description of the topography, The new approach is finally applied to simulate a surge propagating down the irregularly shaped Isar Valley in the Bavarian Alps after a hypothetical dam failure. The straightforward and fully stable computation of the flood hydrograph along the Isar Valley clearly reflects the impact of the strongly varying topographic characteristics on the How phenomenon. Apart from treating surge flow phenomena as a whole, the analytical solution also offers a rigorous alternative to both (a) the approximate Whitham solution, for generating initial values, and (b) the rough volume balance techniques used to model the wave tip in numerical surge flow computations.

Convective instability of 3-D fluid-saturated geological fault zones heated from below

We conduct a theoretical analysis to investigate the convective instability of 3-D fluid-saturated geological fault zones when they are heated uniformly from below. In particular, we have derived exact analytical solutions for the critical Rayleigh numbers of different convective flow structures. Using these critical Rayleigh numbers, three interesting convective flow structures have been identified in a geological fault zone system. It has been recognized that the critical Rayleigh numbers of the system have a minimum value only for the fault zone of infinite length, in which the corresponding convective flow structure is a 2-D slender-circle flow. However, if the length of the fault zone is finite, the convective flow in the system must be 3-D. Even if the length of the fault zone is infinite, since the minimum critical Rayleigh number for the 2-D slender-circle flow structure is so close to that for the 3-D convective flow structure, the system may have almost the same chance to pick up the 3-D convective flow structures. Also, because the convection modes are so close for the 3-D convective flow structures, the convective flow may evolve into the 3-D finger-like structures, especially for the case of the fault thickness to height ratio approaching zero. This understanding demonstrates the beautiful aspects of the present analytical solution for the convective instability of 3-D geological fault zones, because the present analytical solution is valid for any value of the ratio of the fault height to thickness. Using the present analytical solution, the conditions, under which different convective flow structures may take place, can be easily determined.

Theoretical and numerical analyses of convective instability in porous media with temperature-dependent viscosity

Exact analytical solutions of the critical Rayleigh numbers have been obtained for a hydrothermal system consisting of a horizontal porous layer with temperature-dependent viscosity. The boundary conditions considered are constant temperature and zero vertical Darcy velocity at both the top and bottom of the layer. Not only can the derived analytical solutions be readily used to examine the effect of the temperature-dependent viscosity on the temperature-gradient driven convective flow, but also they can be used to validate the numerical methods such as the finite-element method and finite-difference method for dealing with the same kind of problem. The related analytical and numerical results demonstrated that the temperature-dependent viscosity destabilizes the temperature-gradient driven convective flow and therefore, may affect the ore body formation and mineralization in the upper crust of the Earth. Copyright (C) 2003 John Wiley Sons, Ltd.

Reliability of multiple frequency bioelectrical impedance analysis: An intermachine comparison

The technical reliability (i.e., interinstrument and interoperator reliability) of three SEAC-swept frequency bioimpedance monitors was assessed for both errors of measurement and associated analyses. In addition, intraoperator and intrainstrument variability was evaluated for repeat measures over a 4-hour period. The measured impedance values from a range of resistance-capacitance circuits were accurate to within 3% of theoretical values over a range of 50-800 ohms. Similarly, phase was measured over the range 1 degrees-19 degrees with a maximum deviation of 1.3 degrees from the theoretical value. The extrapolated impedance at zero frequency was equally well determined (+/-3%). However, the accuracy of the extrapolated value at infinite frequency was decreased, particularly at impedances below 50 ohms (approaching the lower limit of the measurement range of the instrument). The interinstrument/operator variation for whole body measurements were recorded on human volunteers with biases of less than +/-1% for measured impedance values and less than 3% for phase. The variation in the extrapolated values of impedance at zero and infinite frequencies included variations due to operator choice of the analysis parameters but was still less than +/-0.5%. (C) 1997 Wiley-Liss, Inc.

Tenoxicam controls pain without altering orthodontic movement of maxillary canines

Arantes GM, Arantes VMN, Ashmawi HA, Posso IP To study the efficacy of tenoxicam for pain control, its potential for preemptive analgesia, and its influence on the orthodontic movement of upper canine teeth. This was a randomized controlled double-blind cross-over study. The patients were divided into three groups. Two groups received tenoxicam in daily doses of 20 mg orally for 3 days. Group A received the first dose of the drug before orthodontic activation and group B, just afterwards. Group C (control) received a placebo for 3 days. All groups had access to 750 mg of paracetamol up to four times a day. Three orthodontic activations were performed at 30-day intervals. Each patient belonged to two different groups. Pain intensity was assessed using a descriptive Pain Scale and a Visual Analog Scale. Private clinic; 36 patients undergoing bilateral canine tooth retraction. The statistical analysis did not show any difference in movement between the active groups and the control at any time. There was no statistical difference between the groups that received tenoxicam. Pain intensity in these groups was lower than in the placebo group. The difference in pain intensity between the active groups and the control was greatest at the assessment made 12 h after activation and it tended to zero, 72 h after activation. Tenoxicam did not influence orthodontic movement of the upper canines. It was effective for pain control and did not present any preemptive analgesic effect.