Sausage: protein threading with flexible force fields

Sausage is a protein sequence threading program, but with remarkable run-time flexibility. Using different scripts, it can calculate protein sequence-structure alignments, search structure libraries, swap force fields, create models form alignments, convert file formats and analyse results. There are several different force fields which might be classed as knowledge-based, although they do not rely on Boltzmann statistics. Different force fields are used for alignment calculations and subsequent ranking of calculated models.

Heat flow for Yang-Mills-Higgs fields, part I

The Yang-Mills-Higgs field generalizes the Yang-Mills field. The authors establish the local existence and uniqueness of the weak solution to the heat flow for the Yang-Mills-Higgs field in a vector bundle over a compact Riemannian 4-manifold, and show that the weak solution is gauge-equivalent to a smooth solution and there are at most finite singularities at the maximum existing time.

Saturation of a two-level atom in polychromatic fields

The response of a two-level atom in a strong polychromatic field composed of a large number of equidistant frequency components is investigated. We calculate numerically, as well as analytically,:the stationary population inversion and show that the saturation of the atomic transition strongly depends on whether or not there is a central (resonant) frequency component in the driving field. We find that, in the presence of the central component, the atom can remain in the ground state even for a strong Rabi frequency of the driving field. In addition, we find that the inversion is sensitive to the relative phase between the frequency components. When the central component is suppressed, the atomic transition saturates with the Rabi frequency independent of the relative phase.

Content and context of monocular regions determine perceived depth in random dot, unpaired background and phantom stereograms

Perceived depth was measured for three-types of stereograms with the colour/texture of half-occluded (monocular) regions either similar to or dissimilar to that of binocular regions or background. In a two-panel random dot stereogram the monocular region was filled with texture either similar or different to the far panel or left blank. In unpaired background stereograms the monocular region either matched the background or was different in colour or texture and in phantom stereograms the monocular region matched the partially occluded object or was a different colour or texture. In all three cases depth was considerably impaired when the monocular texture did not match either the background or the more distant surface. The content and context of monocular regions as well as their position are important in determining their role as occlusion cues and thus in three-dimensional layout. We compare coincidence and accidental view accounts of these effects. (C) 2002 Elsevier Science Ltd. All rights reserved.

Variable expression of mental retardation, autism, seizures, and dystonic hand movements in two families with an identical ARX gene mutation

Two families, originally diagnosed as having nonsyndromic X-linked mental retardation (NSXLMR), were reviewed when it was shown that they had a 24-bp duplication (428-45 1dup(24bp)) in the ARX gene [Stromme et al., 2002: Nat Genet 30:441-445]. This same duplication had also been found in three other families: one with X-linked infantile spasms and hypsarrhythmia (X-linked West syndrome, MIM 308350) and two with XLMR and dystonic movements of the hands (Partington syndrome, MIM 309510). On review, manifestations of both West and Partington syndromes were found in some individuals from both families. In addition, it was found that one individual had autism and two had autistic behavior, one of whom had epilepsy. The degree of mental retardation ranged from mild to severe. A GCG trinucleotide expansion (GCG)10+7 and a deletion of 1,517 by in the ARX gene have also been found in association with the West syndrome, and a missense mutation (1058C >T) in a family with a newly recognized form of myoclonic epilepsy, severe mental retardation, and spastic paraplegia [Scheffer et al., 2002: Neurology, in press]. Evidently all these disorders are expressions of mutations in the same gene. It remains to be seen what proportions of patients with infantile spasms, focal dystonia, autism, epilepsy, and nonsyndromic mental retardation are accounted for by mutations in the ARX gene. (C) 2002 Wiley-Liss, Inc.

Modelling the Distribution of Ischaemic Stroke-Specific Survival Time Using an EM-based Mixture Approach with Random Effects Adjustment

A two-component survival mixture model is proposed to analyse a set of ischaemic stroke-specific mortality data. The survival experience of stroke patients after index stroke may be described by a subpopulation of patients in the acute condition and another subpopulation of patients in the chronic phase. To adjust for the inherent correlation of observations due to random hospital effects, a mixture model of two survival functions with random effects is formulated. Assuming a Weibull hazard in both components, an EM algorithm is developed for the estimation of fixed effect parameters and variance components. A simulation study is conducted to assess the performance of the two-component survival mixture model estimators. Simulation results confirm the applicability of the proposed model in a small sample setting. Copyright (C) 2004 John Wiley Sons, Ltd.

Familial partial epilepsy with variable foci: Clinical features and linkage to chromosome 22q12

Background: Familial partial epilepsy with variable foci (FPEVF) is an autosomal dominant syndrome characterized by partial seizures originating from different brain regions in different family members in the absence of detectable structural abnormalities. A gene for FPEVF was mapped to chromosome 22q12 in two distantly related French-Canadian families. Methods: We describe the clinical features and performed a linkage analysis in a Spanish kindred and in a third French-Canadian family distantly related to the original pedigrees. Results: Onset of seizures was typically in middle childhood, and attacks were usually easy to control. Seizure semiology varied among family members but was constant for each individual. In some, a pattern of nocturnal frontal lobe seizures led to consideration of the diagnosis of autosomal dominant nocturnal frontal lobe epilepsy (ADNFLE). The Spanish family was mapped to chromosome 22q (multipoint lod score, 3.4), and the new French-Canadian family had a multipoint lod score of 2.97 and shared the haplotype of the original French-Canadian families. Conclusions: Identification of the various forms of familial partial epilepsy is challenging, particularly in small families, in which insufficient individuals exist to identify a specific pattern. We provide clinical guidelines for this task, which will eventually be supplanted by specific molecular diagnosis. We confirmed linkage of FPEVF to chromosome 22q 12 and redefined the region to a 5.2-Mb segment of DNA.

Variable effects of larval size on post-metamorphic performance in the field

Larval quality may be capable of explaining much of the variation in the recruitment and subsequent population dynamics of benthic marine invertebrates. Whilst the effects of larval nutritional condition on adult performance have received the most attention, recent work has shown that larval size may also be an important and ubiquitous source of variation in larval quality. We examined the effects of variation in larval size on the post-metamorphic survival and growth of Watersipora subtorquata in 2 very different habitats - experimental substrata and pier pilings. We found strong effects of larval size on colony performance, although these varied among experiments. For colonies on experimental substrata, larval size positively affected adult survival and, initially, growth. However, after 3 wk in the field, there was no relationship between larval size and colony size, possibly because colonies were completely surrounded by newly settled organisms. Larval size also positively affected post-metamorphic growth of colonies on pier pilings, but, surprisingly, colonies that came from larger larvae had lower survival than colonies from smaller larvae. Overall, variation in larval size will strongly affect the recruitment and subsequent performance of adults in this species, although this may vary among different habitats. This study highlights the importance of examining the effects of larval quality on adult performance in as realistic conditions as possible, because of the strong interaction between larval size effects and the environment.

Characterization of optically driven fluid stress fields with optical tweezers

We present a controlled stress microviscometer with applications to complex fluids. It generates and measures microscopic fluid velocity fields, based on dual beam optical tweezers. This allows an investigation of bulk viscous properties and local inhomogeneities at the probe particle surface. The accuracy of the method is demonstrated in water. In a complex fluid model (hyaluronic acid), we observe a strong deviation of the flow field from classical behavior. Knowledge of the deviation together with an optical torque measurement is used to determine the bulk viscosity. Furthermore, we model the observed deviation and derive microscopic parameters.

Three-dimensional model for multi-component reactive transport with variable density groundwater flow

PHWAT is a new model that couples a geochemical reaction model (PHREEQC-2) with a density-dependent groundwater flow and solute transport model (SEAWAT) using the split-operator approach. PHWAT was developed to simulate multi-component reactive transport in variable density groundwater flow. Fluid density in PHWAT depends not on only the concentration of a single species as in SEAWAT, but also the concentrations of other dissolved chemicals that can be subject to reactive processes. Simulation results of PHWAT and PHREEQC-2 were compared in their predictions of effluent concentration from a column experiment. Both models produced identical results, showing that PHWAT has correctly coupled the sub-packages. PHWAT was then applied to the simulation of a tank experiment in which seawater intrusion was accompanied by cation exchange. The density dependence of the intrusion and the snow-plough effect in the breakthrough curves were reflected in the model simulations, which were in good agreement with the measured breakthrough data. Comparison simulations that, in turn, excluded density effects and reactions allowed us to quantify the marked effect of ignoring these processes. Next, we explored numerical issues involved in the practical application of PHWAT using the example of a dense plume flowing into a tank containing fresh water. It was shown that PHWAT could model physically unstable flow and that numerical instabilities were suppressed. Physical instability developed in the model in accordance with the increase of the modified Rayleigh number for density-dependent flow, in agreement with previous research. (c) 2004 Elsevier Ltd. All rights reserved.

Resonance fluorescence and Autler-Townes spectra of a two-level atom driven by two fields of equal frequencies

We study the effects of driving a two-level atom by two intense field modes that have equal frequencies but are otherwise distinguishable; the intensity of one mode is also assumed to be greater than that of the other. We calculate first the dressed states of the system, and then its resonance fluorescence and Autler-Townes absorption spectra. We find that the energy spectrum of the doubly dressed atom consists of a ladder of doublet continua. These continua manifest themselves in the fluorescence spectrum, where they produce continua at the positions of the Mellow sideband frequencies omega(L)+/-2 Omega of the strong field, and in the Autler-Townes absorption spectrum, which becomes a two-continuum doublet.

Calculating current densities and fields produced by shielded magnetic resonance imaging probes

A method is presented for computing the fields produced by radio frequency probes of the type used in magnetic resonance imaging. The effects of surrounding the probe with a shielding coil, intended to eliminate stray fields produced outside the probe, are included. An essential feature of these devices is the fact that the conducting rungs of the probe are of finite width relative to the coil radius, and it is therefore necessary to find the distribution of current within the conductors as part of the solution process. This is done here using a numerical method based on the inverse finite Hilbert transform, applied iteratively to the entire structure including its shielding coils. It is observed that the fields are influenced substantially by the width of the conducting rungs of the probe, since induced eddy currents within the rungs become more pronounced as their width is increased. The shield is also shown to have a significant effect on both the primary current density and the resultant fields. Quality factors are computed for these probes and compared with values measured experimentally.

Conformational dynamics of thyroid hormones by variable temperature nuclear magnetic resonance: The role of side chain rotations and cisoid/transoid interconversions

H-1 NMR spectra of the thyroid hormone thyroxine recorded at low temperature and high field show splitting into two peaks of the resonance due to the H2,6 protons of the inner (tyrosyl) ring. A single resonance is observed in 600 MHz spectra at temperatures above 185 K. An analysis of the line shape as a function of temperature shows that the coalescence phenomenon is due to an exchange process with a barrier of 37 kJ mol(-1). This is identical to the barrier for coalescence of the H2',6' protons of the outer (phenolic) ring reported previously for the thyroid hormones and their analogues. It is proposed that the separate peaks at low temperature are due to resonances for H2,6 in cisoid and transoid conformers which are populated in approximately equal populations. These two peaks are averaged resonances for the individual H2 and H6 protons. Conversion of cisoid to transoid forms can occur via rotation of either the alanyl side chain or the outer ring, from one face of the inner ring to the other. It is proposed that the latter process is the one responsible for the observed coalescence phenomenon. The barrier to rotation of the alanyl side chain is greater than or equal to 37 kJ mol(-1), which is significantly larger than has previously been reported for Csp(2)-Csp(3) bonds in other Ph-CH2-X systems. The recent crystal structure of a hormone agonist bound to the ligand-binding domain of the rat thyroid hormone receptor (Wagner et al. Nature 1995, 378, 690-697) shows the transoid form to be the bound conformation. The significant energy barrier to cisoid/transoid interconversion determined in the current study combined with the tight fit of the hormone to its receptor suggests that interconversion between the forms cannot occur at the receptor site but that selection for the preferred bound form occurs from the 50% population of the transoid form in solution.

Rapid computation of static fields produced by thick circular solenoids

A straightforward method is proposed for computing the magnetic field produced by a circular coil that contains a large number of turns wound onto a solenoid of rectangular cross section. The coil is thus approximated by a circular ring containing a continuous constant current density, which is very close to the real situation when sire of rectangular cross section is used. All that is required is to evaluate two functions, which are defined as integrals of periodic quantities; this is done accurately and efficiently using trapezoidal-rule quadrature. The solution can be obtained so rapidly that this procedure is ideally suited for use in stochastic optimization, An example is given, in which this approach is combined with a simulated annealing routine to optimize shielded profile coils for NMR.