Protein fold recognition score functions: Unusual construction strategies

We describe two ways of optimizing score functions for protein sequence to structure threading. The first method adjusts parameters to improve sequence to structure alignment. The second adjusts parameters so as to improve a score function's ability to rank alignments calculated in the first score function. Unlike those functions known as knowledge-based force fields, the resulting parameter sets do not rely on Boltzmann statistics, have no claim to representing free energies and are purely constructions for recognizing protein folds. The methods give a small improvement, but suggest that functions can be profitably optimized for very specific aspects of protein fold recognition, Proteins 1999;36:454-461. (C) 1999 Wiley-Liss, Inc.

Structure/function studies of S100A8/A9

The functional importance of members of the S100 Ca2+-binding protein family is recently emerging. A variety of activities, several of which are apparently opposing, are attributed to S100A8, a protein implicated in embryogenesis, growth, differentiation, and immune and inflammatory processes. Murine (m) S100A8 was initially described as a chemoattractant (CP-10) for myeloid cells. It is coordinately expressed with mS100A9 (MRP14) in neutrophils and the non-covalent heterodimer is presumed to be the functional intracellular species. The extracellular chemotactic activity of mS100A8, however, is not dependent on mS100A9 and occurs at concentrations (10(-13)-10(-11) M) at which the non-covalent heterodimer would probably dissociate. This review focuses on the structure and post-translational modifications of mS100A8/A9 and their effects on function, particularly chemotaxis.

Protein sequence threading, the alignment problem, and a two-step strategy

Conventionally, protein structure prediction via threading relies on some nonoptimal method to align a protein sequence to each member of a library of known structures. We show how a score function (force field) can be modified so as to allow the direct application of a dynamic programming algorithm to the problem. This involves an approximation whose damage can be minimized by an optimization process during score function parameter determination. The method is compared to sequence to structure alignments using a more conventional pair-wise score function and the frozen approximation. The new method produces results comparable to the frozen approximation, but is faster and has fewer adjustable parameters. It is also free of memory of the template's original amino acid sequence, and does not suffer from a problem of nonconvergence, which can be shown to occur with the frozen approximation. Alignments generated by the simplified score function can then be ranked using a second score function with the approximations removed. (C) 1999 John Wiley & Sons, Inc.

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.

Assessment of limb muscle and adipose tissue by dual-energy X-ray absorptiometry using magnetic resonance imaging for comparison

OBJECTIVE: To use magnetic resonance imaging (MRI) to validate estimates of muscle and adipose tissue (AT) in lower limb sections obtained by dual-energy X-ray absorptiometry (DXA) modelling. DESIGN: MRI measurements were used as reference for validating limb muscle and AT estimates obtained by DXA models that assume fat-free soft tissue (FFST) comprised mainly muscle: model A accounted for bone hydration only; model B also applied constants for FFST in bone and skin and fat in muscle and AT; model C was as model B but allowing for variable fat in muscle and AT. SUBJECTS: Healthy men (n = 8) and women (n = 8), ages 41 - 62 y; mean (s.d.) body mass indices (BMIs) of 28.6 (5.4) kg/m(2) and 25.1 (5.4) kg/m2, respectively. MEASUREMENTS: MRI scans of the legs and whole body DXA scans were analysed for muscle and AT content of thigh (20 cm) and lower leg (10 cm) sections; 24 h creatinine excretion was measured. RESULTS: Model A overestimated thigh muscle volume (MRI mean, 2.3 l) substantially (bias 0.36 l), whereas model B underestimated it by only 2% (bias 0.045 l). Lower leg muscle (MRI mean, 0.6 l) was better predicted using model A (bias 0.04 l, 7% overestimate) than model B (bias 0.1 l, 17% underestimate). The 95% limits of agreement were high for these models (thigh,+/- 20%; lower leg,+/- 47%). Model C predictions were more discrepant than those of model B. There was generally less agreement between MRI and all DXA models for AT. Measurement variability was generally less for DXA measurements of FFST (coefficient of variation 0.7 - 1.8%) and fat (0.8 - 3.3%) than model B estimates of muscle (0.5-2.6%) and AT (3.3 - 6.8%), respectively. Despite strong relationships between them, muscle mass was overestimated by creatinine excretion with highly variable predictability. CONCLUSION: This study has shown the value of DXA models for assessment of muscle and AT in leg sections, but suggests the need to re-evaluate some of the assumptions upon which they are based.

Visual function: How spiders find the right rock to crawl under

Animals that go on hunting expeditions face the problem of finding the way home at the end of the day. A group of hunting spiders has now been added to the list of animals that use the celestial pattern of polarized light as a compass for navigation. (C) 1999 Elsevier Science Ltd. All rights reserved.

What the evolution of the amyloid protein precursor supergene family tells us about its function

The Alzheimer's disease amyloid protein precursor (APP) gene is part of a multi-gene super-family from which sixteen homologous amyloid precursor-like proteins (APLP) and APP species homologues have been isolated and characterised. Comparison of exon structure (including the uncharacterised APL-1 gene), construction of phylogenetic trees, and analysis of the protein sequence alignment of known homologues of the APP super-family were performed to reconstruct the evolution of the family and to assess the functional significance of conserved protein sequences between homologues. This analysis supports an adhesion function for all members of the APP super family, with specificity determined by those sequences which are not conserved between APLP lineages, and provides evidence for an increasingly complex APP superfamily during evolution. The analysis also suggests that Drosophila APPL and Caenorhabdotids elegans APL-1 may be a fourth APLP lineage indicating that these proteins, while not functional homologues of human APP, are similarly likely to regulate cell adhesion. Furthermore, the beta A4 sequence is highly conserved only in APP orthologues, strongly suggesting this sequence is of significant functional importance in this lineage. (C) 2000 Elsevier Science Ltd. All rights reserved.

Two-dimensional nonlinear dynamics of evanescent-wave guided atoms in a hollow fiber

We describe the classical and quantum two-dimensional nonlinear dynamics of large blue-detuned evanescent-wave guiding cold atoms in hollow fiber. We show that chaotic dynamics exists for classic dynamics, when the intensity of the beam is periodically modulated. The two-dimensional distributions of atoms in (x,y) plane are simulated. We show that the atoms will accumulate on several annular regions when the system enters a regime of global chaos. Our simulation shows that, when the atomic flux is very small, a similar distribution will be obtained if we detect the atomic distribution once each the modulation period and integrate the signals. For quantum dynamics, quantum collapses, and revivals appear. For periodically modulated optical potential, the variance of atomic position will be suppressed compared to the no modulation case. The atomic angular momentum will influence the evolution of wave function in two-dimensional quantum system of hollow fiber.

Quantum chaos in atom optics

I shall discuss the quantum and classical dynamics of a class of nonlinear Hamiltonian systems. The discussion will be restricted to systems with one degree of freedom. Such systems cannot exhibit chaos, unless the Hamiltonians are time dependent. Thus we shall consider systems with a potential function that has a higher than quadratic dependence on the position and, furthermore, we shall allow the potential function to be a periodic function of time. This is the simplest class of Hamiltonian system that can exhibit chaotic dynamics. I shall show how such systems can be realized in atom optics, where very cord atoms interact with optical dipole potentials of a far-off resonance laser. Such systems are ideal for quantum chaos studies as (i) the energy of the atom is small and action scales are of the order of Planck's constant, (ii) the systems are almost perfectly isolated from the decohering effects of the environment and (iii) optical methods enable exquisite time dependent control of the mechanical potentials seen by the atoms.

Unifying principles in terrestrial locomotion: Do hopping Australian marsupials fit in?

Mammalian terrestrial locomotion has many unifying principles. However, the Macropodoidea are a particularly interesting group that exhibit a number of significant deviations from the principles that seem to apply to other mammals. While the properties of materials that comprise the musculoskeletal system of mammals are similar, evidence suggests that tendon properties in macropodoid marsupials may be size or function dependent, in contrast to the situation in placental mammals. Postural differences related to hopping versus running have a dramatic effect on the scaling of the pelvic limb musculoskeletal system. Ratios of muscle fibre to tendon cross-sectional areas for ankle extensors and digital flexors scale with positive allometry in all mammals, but exponents are significantly higher in macropods. Tendon safety factors decline with increasing body mass in mammals, with eutherians at risk of ankle extensor tendon rupture at a body mass of about 150 kg, whereas kangaroos encounter similar problems at a body mass of approximately 35 kg. Tendon strength appears to limit locomotor performance in these animals. Elastic strain energy storage in tendons is mass dependent in all mammals, but exponents are significantly larger in macropodid. Tibial stresses may scale with positive allometry in kangaroos, which result in lower bone safety factors in macropods compared to eutherian mammals.

Accelerating seismic energy release and evolution of event time and size statistics: Results from two heterogeneous cellular automaton models

The evolution of event time and size statistics in two heterogeneous cellular automaton models of earthquake behavior are studied and compared to the evolution of these quantities during observed periods of accelerating seismic energy release Drier to large earthquakes. The two automata have different nearest neighbor laws, one of which produces self-organized critical (SOC) behavior (PSD model) and the other which produces quasi-periodic large events (crack model). In the PSD model periods of accelerating energy release before large events are rare. In the crack model, many large events are preceded by periods of accelerating energy release. When compared to randomized event catalogs, accelerating energy release before large events occurs more often than random in the crack model but less often than random in the PSD model; it is easier to tell the crack and PSD model results apart from each other than to tell either model apart from a random catalog. The evolution of event sizes during the accelerating energy release sequences in all models is compared to that of observed sequences. The accelerating energy release sequences in the crack model consist of an increase in the rate of events of all sizes, consistent with observations from a small number of natural cases, however inconsistent with a larger number of cases in which there is an increase in the rate of only moderate-sized events. On average, no increase in the rate of events of any size is seen before large events in the PSD model.

Free energy approximations in simple lattice proteins

This work addresses the question of whether it is possible to define simple pairwise interaction terms to approximate free energies of proteins or polymers. Rather than ask how reliable a potential of mean force is, one can ask how reliable it could possibly be. In a two-dimensional, infinite lattice model system one can calculate exact free energies by exhaustive enumeration. A series of approximations were fitted to exact results to assess the feasibility and utility of pairwise free energy terms. Approximating the true free energy with pairwise interactions gives a poor fit with little transferability between systems of different size. Adding extra artificial terms to the approximation yields better fits, but does not improve the ability to generalize from one system size to another. Furthermore, one cannot distinguish folding from nonfolding sequences via the approximated free energies. Most usefully, the methodology shows how one can assess the utility of various terms in lattice protein/polymer models. (C) 2001 American Institute of Physics.