Estimating and evaluating the statistics of gapped local-alignment scores

We present a novel maximum-likelihood-based algorithm for estimating the distribution of alignment scores from the scores of unrelated sequences in a database search. Using a new method for measuring the accuracy of p-values, we show that our maximum-likelihood-based algorithm is more accurate than existing regression-based and lookup table methods. We explore a more sophisticated way of modeling and estimating the score distributions (using a two-component mixture model and expectation maximization), but conclude that this does not improve significantly over simply ignoring scores with small E-values during estimation. Finally, we measure the classification accuracy of p-values estimated in different ways and observe that inaccurate p-values can, somewhat paradoxically, lead to higher classification accuracy. We explain this paradox and argue that statistical accuracy, not classification accuracy, should be the primary criterion in comparisons of similarity search methods that return p-values that adjust for target sequence length.

Determinants of tissue Doppler measures of regional diastolic function during dobutamine stress echocardiography

Background Diastolic dysfunction induced by ischemia may alter transmitral blood flow, but this reflects global ventricular function, and pseudonormalization may occur with increased preload. Tissue Doppler may assess regional diastolic function and is relatively load-independent, but limited data exist regarding its application to stress testing. We sought to examine the stress response of regional diastolic parameters to dobutomine echocardiography (DbE). Methods Sixty-three patients underwent study with DbE: 20 with low probability of coronary artery disease (CAD) and 43 with CAD who underwent angiography. A standard DbE protocol was used, and segments were categorized as ischemic, scar, or normal. Color tissue Doppler was acquired at baseline and peak stress, and waveforms in the basal and mid segments were used to measure early filling (Em), late filling (Am), and E deceleration time. Significant CAD was defined by stenoses >50% vessel diameter. Results Diastolic parameters had limited feasibility because of merging of Em and Am waves at high heart rates and limited reproducibility. Nonetheless, compared with normal segments, segments subtended with significant stenoses showed a lower Em velocity at rest (6.2 +/- 2.6 cm/s vs 4.8 +/- 2.2 cm/s, P < .0001) and peak (7.5 +/- 4.2 cm/s vs 5.1 +/- 3.6 cm/s, P < .0001), Abnormal segments also showed a shorter E deceleration time (51 +/- 27 ms vs 41 +/- 27 ms, P = .0001) at base and peak. No changes were documented in Am. The same pattern was seen with segments identified as ischemic with wall motion score. However, in the absence of ischemia, segments of patients with left ventricular hypertrophy showed a lower Em velocity, with blunted Em responses to stress. Conclusion Regional diastolic function is sensitive to ischemia. However, a number of practical limitations limit the applicability of diastolic parameters for the quantification of stress echocardiography.

Development of a fully quantitative approach to the interpretation of stress echocardiography using radial and longitudinal myocardial velocities

Background: Tissue Doppler may be used to quantify regional left ventricular function but is limited by segmental variation of longitudinal velocity from base to apex and free to septal walls. We sought to overcome this by developing a composite of longitudinal and radial velocities. Methods and Results. We examined 82 unselected patients undergoing a standard dobutamine echocardiogram. Longitudinal velocity was obtained in the basal and mid segments of each wall using tissue Doppler in the apical views. Radial velocities were derived in the same segments using an automated border detection system and centerline method with regional chords grouped according to segment location and temporally averaged. In 25 patients at low probability of coronary disease, the pattern of regional variation in longitudinal velocity (higher in the septum) was the opposite of radial velocity (higher in the free wall) and the combination was homogenous. In 57 patients undergoing angiography, velocity in abnormal segments was less than normal segments using longitudinal (6.0 +/- 3.6 vs 9.0 +/- 2.2 cm/s, P = .01) and radial velocity (6.0 +/- 4.0 vs 8.0 +/- 3.9 cm/s, P = .02). However, the composite velocity permitted better separation of abnormal and normal segments (13.3 +/- 5.6 vs 17.5 +/- 4.2 cm/s, P = .001). There was no significant difference between the accuracy of this quantitative approach and expert visual wall motion analysis (81% vs 84%, P = .56). Conclusion: Regional variation of uni-dimensional myocardial velocities necessitates site-specific normal ranges, probably because of different fiber directions. Combined analysis of longitudinal and radial velocities allows the derivation of a composite velocity, which is homogenous in all segments and may allow better separation of normal and abnormal myocardium.

Quantum dynamics of two coupled qubits

We investigate the difference between classical and quantum dynamics of coupled magnetic dipoles. We prove that in general the dynamics of the classical interaction Hamiltonian differs from the corresponding quantum model, regardless of the initial state. The difference appears as nonpositive-definite diffusion terms in the quantum evolution equation of an appropriate positive phase-space probability density. Thus, it is not possible to express the dynamics in terms of a convolution of a positive transition probability function and the initial condition as can be done in the classical case. It is this feature that enables the quantum system to evolve to an entangled state. We conclude that the dynamics are a quantum element of nuclear magnetic resonance quantum-information processing. There are two limits where our quantum evolution coincides with the classical one: the short-time limit before spin-spin interaction sets in and the long-time limit when phase diffusion is incorporated.

Effects of paternal drinking, conduct disorder and childhood home environment on the development of alcohol use disorders in a Thai population

Aims To identify influences on the development of alcohol use disorders in a Thai population, particularly parental drinking and childhood environment. Design Case-control study. Setting A university hospital, a regional hospital and a community hospital in southern Thailand. Participants Ninety-one alcohol-dependents and 177 hazardous/harmful drinkers were recruited as cases and 144 non-or infrequent drinkers as controls. Measurements Data on parental drinking, family demographic characteristics, family activities, parental disciplinary practice, early religious life and conduct disorder were obtained using a structured interview questionnaire. The main outcome measure was the subject's classification as alcohol-dependent, hazardous/harmful drinker or non-/infrequent drinker. Findings A significant relationship was found between having a drinking father and the occurrence of hazardous/harmful drinking or alcohol dependence in the subjects. Childhood factors (conduct disorder and having been a temple boy, relative probability ratios, RPRs and 95% CI: 6.39, 2.81-14.55 and 2.21, 1.19-4.08, respectively) also significantly predicted alcohol dependence, while perceived poverty and ethnic alienation was reported less frequently by hazardous/harmful drinkers and alcohol-dependents (RPRS and 95% CIs = 0.34, 0.19-0.62 and 0.59, 0.38-0.93, respectively) than the controls. The relative probability ratio for the effect of the father's infrequent drinking on the son's alcohol dependence was 2.92 (95% CI = 1.42-6.02) and for the father's heavy or dependent drinking 2.84 (95% CI=1.31-6.15). Conclusions Being exposed to a light-drinking, father increases the risk of a son's alcohol use disorders exhibited either as hazardous-harmful or dependent drinking. However, exposure to a heavy- or dependent-drinking father is associated more uniquely with an increased risk of his son being alcohol-dependent. The extent to which this is seen in other cultures is worthy of exploration.

Discovery and structures of the cyclotides: novel macrocyclic peptides from plants

Circular disulfide-rich polypeptides were unknown a decade ago but over recent years a large family of such molecules has been discovered, which we now refer to as the cyclotides. They are typically about 30 amino acids in size, contain an N- to C-cyclised backbone and incorporate three disulfide bonds arranged in a cystine knot motif. In this motif, an embedded ring in the structure formed by two disulfide bonds and their connecting backbone segments is penetrated by the third disulfide bond. The combination of this knotted and strongly braced structure with a circular backbone renders the cyclotides impervious to enzymatic breakdown and makes them exceptionally stable. This article describes the discovery of the cyclotides in plants from the Rubiaceae and Violaceae families, their chemical synthesis, folding, structural characterisation, and biosynthetic origin. The cyclotides have a diverse range of biological applications, ranging from uterotonic action, to anti-HIV and neurotensin antagonism. Certain plants from which they are derived have a history of uses in native medicine, with activity being observed after oral ingestion of a tea made from the plants. This suggests the possibility that the cyclotides may be orally bioavailable. They therefore have a range of potential applications as a stable peptide framework.

An atomistic simulation of solid state sintering using Monte Carlo methods

This paper presents results on the simulation of the solid state sintering of copper wires using Monte Carlo techniques based on elements of lattice theory and cellular automata. The initial structure is superimposed onto a triangular, two-dimensional lattice, where each lattice site corresponds to either an atom or vacancy. The number of vacancies varies with the simulation temperature, while a cluster of vacancies is a pore. To simulate sintering, lattice sites are picked at random and reoriented in terms of an atomistic model governing mass transport. The probability that an atom has sufficient energy to jump to a vacant lattice site is related to the jump frequency, and hence the diffusion coefficient, while the probability that an atomic jump will be accepted is related to the change in energy of the system as a result of the jump, as determined by the change in the number of nearest neighbours. The jump frequency is also used to relate model time, measured in Monte Carlo Steps, to the actual sintering time. The model incorporates bulk, grain boundary and surface diffusion terms and includes vacancy annihilation on the grain boundaries. The predictions of the model were found to be consistent with experimental data, both in terms of the microstructural evolution and in terms of the sintering time. (C) 2002 Elsevier Science B.V. All rights reserved.

Levene tests of homogeneity of variance for general block and treatment designs

This article develops a weighted least squares version of Levene's test of homogeneity of variance for a general design, available both for univariate and multivariate situations. When the design is balanced, the univariate and two common multivariate test statistics turn out to be proportional to the corresponding ordinary least squares test statistics obtained from an analysis of variance of the absolute values of the standardized mean-based residuals from the original analysis of the data. The constant of proportionality is simply a design-dependent multiplier (which does not necessarily tend to unity). Explicit results are presented for randomized block and Latin square designs and are illustrated for factorial treatment designs and split-plot experiments. The distribution of the univariate test statistic is close to a standard F-distribution, although it can be slightly underdispersed. For a complex design, the test assesses homogeneity of variance across blocks, treatments, or treatment factors and offers an objective interpretation of residual plot.

Behavior of argon gas release from manganese oxide minerals as revealed by Ar-40/Ar-39 laser incremental heating analysis

Manganese oxides in association with paleo-weathering may provide significant insights into the multiple factors affecting the formation and evolution of weathering profiles, such as temperature, precipitation, and biodiversity. Laser probe step-heating analysis of supergene hollandite and cryptomelane samples collected from central Queensland, Australia, yield well-defined plateaus and consistent isochron ages, confirming the feasibility, dating very-fined supergene manganese oxides by Ar-40/(39) Ar technique. Two distinct structural sites hosting Ar isotopes can be identified in light of their degassing behaviors obtained by incremental heating analyses. The first site, releasing its gas fraction at the laser power 0.2-0.4 W, yields primarily Ar-40(atm), Ar-38(atm), and Ar-36(atm), (atmospheric Ar isotopes). The second sites yield predominantly Ar-40* (radiogenic Ar-40), Ar-39(K), and Ar-38(K) (nucleogenic components), at similar to0.5-1.0 W. There is no significant Ar gas released at the laser power higher than 1.0 W, indicating the breakdown of the tunnel sites hosting the radiogenic and nucleogenic components. The excellent match between the degassing behaviors of Ar-40*, Ar-39(K), and Ar-38(K) suggests that these isotopes occupy the same crystallographic sites and that Ar-39(K) loss from the tunnel site by recoil during neutron irradiation and/or bake-out procedure preceding isotopic analysis does not occur. Present investigation supports that neither the overwhelming atmospheric Ar-40 nor the very-fined nature of the supergene manganese oxides poses problems in extracting meaningful weathering geo-chronological information by analyzing supergene manganese oxides minerals.

Risk and protective factors for depressive symptomatology among a community sample of adolescents and young adults

Objectives: The study explores the risk and protective factors for current depressive symptomatology in a large community sample of 15-to-24-year-olds. Methods: The study was designed as a cross-sectional household survey, which used telephone recruitment followed by an anonymous self-report postal questionnaire. The final sample included 3,082 adolescents and young adults from Queensland, Australia. Results: The vast majority of measured risk and protective factors were associated with current depressive symptomatology. Key risk factors included high levels of neuroticism, perceived problems with parents, sexual abuse, relationship breakups, educational failure and sexual identity conflict. A different profile of protective factors was evident for each of these high-risk groups. Of particular note was the importance of well-developed intrapersonal skills as protective for both males and females. The significance of social connectedness as a protective factor for the males overall and across a range of high-risk groups was a central finding. Conclusions and implications: The implications of these findings in relation to a range of mental health promotion and mental illness prevention and early intervention initiatives are discussed. Supported initiatives include parenting programs that consider the realities of modern families, increasing community awareness of the impact on young people of the breakdown of their intimate relationships, initiatives in educational settings and workplaces to increase tolerance of gay/lesbian and bisexual lifestyles and the enhancement of social connectedness.

A probabilistic analysis of Clostridium perfringens growth during food service operations

The purpose of this study was threefold: first, the study was designed to illustrate the use of data and information collected in food safety surveys in a quantitative risk assessment. In this case, the focus was on the food service industry; however, similar data from other parts of the food chain could be similarly incorporated. The second objective was to quantitatively describe and better understand the role that the food service industry plays in the safety of food. The third objective was to illustrate the additional decision-making information that is available when uncertainty and variability are incorporated into the modelling of systems. (C) 2002 Elsevier Science B.V. All rights reserved.

An on-line investigation of lexical ambiguity processing in schizophrenia

The processing of lexical ambiguity in context was investigated in eight individuals with schizophrenia and a matched control group. Participants made speeded lexical decisions on the third word in auditory word triplets representing concordant (coin-bank-money), discordant (river-bank-money). neutral (day-bank-money), and unrelated (river-day-money) conditions. When the interstimulus interval (ISI) between the words was 100 ms. individuals with schizophrenia demonstrated priming consistent with selective. context-based lexical activation. At 1250 ms ISI a pattern of nonselective meaning facilitation was obtained. These results suggest an attentional breakdown in the sustained inhibition of meanings on the basis of lexical context. (C) 2002 Elsevier Science (USA).

Progressive Age-Related Changes Similar to Age-Related Macular Degeneration in a Transgenic Mouse Model

Age-related macular degeneration (AMD) is the major cause of blindness in the developed world. its pathomechanism is unknown and its late onset, complex genetics and strong environmental components have all hampered investigations. Here we demonstrate the development of an animal model for AMD that reproduces features associated with geographic atrophy, a transgenic mouse line (mcd/mcd) expressing a mutated form of cathepsin D that is enzymatically inactive thus impairing processing of phagocytosed photoreceptor outer segments in the retinal pigment epithelial (RPE) cells. Pigmentary changes indicating RPE cell atrophy and a decreased response to flash electroretinograms were observed in 11- to 12-month-old mcd/mcd mice. Histological studies showed RPE cell proliferation, photoreceptor degeneration, shortening of photoreceptor outer segments, and accumulation of immunoreactive photoreceptor breakdown products in the RPE cells. An accelerated photoreceptor cell death was detected in 12-month-old mcd/mcd mice. Transmission electron microscopy demonstrated presence of basal laminar and linear deposits that are considered to be the hallmarks of AMD. Small hard drusen associated with human age-related maculopathy were absent in the mcd/mcd mouse model at the ages analyzed. in summary, this model presents several features of AMD, thus providing a valuable tool for investigating the underlying biological processes and pathomechanism of AMD.

Prediction of Golgi Type II membrane proteins based on their transmembrane domains

Motivation: A major issue in cell biology today is how distinct intracellular regions of the cell, like the Golgi Apparatus, maintain their unique composition of proteins and lipids. The cell differentially separates Golgi resident proteins from proteins that move through the organelle to other subcellular destinations. We set out to determine if we could distinguish these two types of transmembrane proteins using computational approaches. Results: A new method has been developed to predict Golgi membrane proteins based on their transmembrane domains. To establish the prediction procedure, we took the hydrophobicity values and frequencies of different residues within the transmembrane domains into consideration. A simple linear discriminant function was developed with a small number of parameters derived from a dataset of Type II transmembrane proteins of known localization. This can discriminate between proteins destined for Golgi apparatus or other locations (post-Golgi) with a success rate of 89.3% or 85.2%, respectively on our redundancy-reduced data sets.