A novel parallel approach to the likelihood-based estimation of admixture in population genetics

Inferring population admixture from genetic data and quantifying it is a difficult but crucial task in evolutionary and conservation biology. Unfortunately state-of-the-art probabilistic approaches are computationally demanding. Effectively exploiting the computational power of modern multiprocessor systems can thus have a positive impact to Monte Carlo-based simulation of admixture modeling. A novel parallel approach is briefly described and promising results on its message passing interface (MPI)-based C implementation are reported.

Hepatitis C virus NS5A protein binds the SH3 domain of the Fyn tyrosine kinase with high affinity: mutagenic analysis of residues within the SH3 domain that contribute to the interaction

Background: The hepatitis C virus (HCV) non-structural 5A protein (NS5A) contains a highly conserved C-terminal polyproline motif with the consensus sequence Pro-X-X- Pro-X-Arg that is able to interact with the Src-homology 3 (SH3) domains of a variety of cellular proteins. Results: To understand this interaction in more detail we have expressed two N-terminally truncated forms of NS5A in E. coli and examined their interactions with the SH3 domain of the Src-family tyrosine kinase, Fyn. Surface plasmon resonance analysis revealed that NS5A binds to the Fyn SH3 domain with what can be considered a high affinity SH3 domain-ligand interaction (629 nM), and this binding did not require the presence of domain I of NS5A (amino acid residues 32-250). Mutagenic analysis of the Fyn SH3 domain demonstrated the requirement for an acidic cluster at the C-terminus of the RT-Src loop of the SH3 domain, as well as several highly conserved residues previously shown to participate in SH3 domain peptide binding. Conclusion: We conclude that the NS5A: Fyn SH3 domain interaction occurs via a canonical SH3 domain binding site and the high affinity of the interaction suggests that NS5A would be able to compete with cognate Fyn ligands within the infected cell.

Metalloestrogens: an emerging class of inorganic xenoestrogens with potential to add to the oestrogenic burden of the human breast

Many compounds in the environment have been shown capable of binding to cellular oestrogen receptors and then mimicking the actions of physiological oestrogens. The widespread origin and diversity in chemical structure of these environmental oestrogens is extensive but to date such compounds have been organic and in particular phenolic or carbon ring structures of varying structural complexity. Recent reports of the ability of certain metal ions to also bind to oestrogen receptors and to give rise to oestrogen agonist responses in vitro and in vivo has resulted in the realisation that environmental oestrogens can also be inorganic and such xenoestrogens have been termed metalloestrogens. This report highlights studies which show metalloestrogens to include aluminium, antimony, arsenite, barium, cadmium, chromium (Cr(II)), cobalt, copper, lead, mercury, nickel, selenite, tin and vanadate. The potential for these metal ions to add to the burden of aberrant oestrogen signalling within the human breast is discussed. Copyright (c) 2006 John Wiley & Sons, Ltd.

Argon predissociation spectroscopy of the OH-center dot H2O and Cl-center dot H2O complexes in the 1000-1900 cm(-1) region: Intramolecular bending transitions and the search for the shared-proton fundamental in the hydroxide monohydrate

We present argon predissociation vibrational spectra of the OH-.H2O and Cl-.H2O complexes in the 1000-1900 cm(-1) energy range, far below the OH stretching region reported in previous studies. This extension allows us to explore the fundamental transitions of the intramolecular bending vibrations associated with the water molecule, as well as that of the shared proton inferred from previous assignments of overtones in the higher energy region. Although the water bending fundamental in the Cl-.H2O spectrum is in very good agreement with expectations, the OH-.H2O spectrum is quite different than anticipated, being dominated by a strong feature at 1090 cm(-1). New full-diniensionality calculations of the OH-.H2O vibrational level structure using diffusion Monte Carlo and the VSCF/CI methods indicate this band arises from excitation of the shared proton.

Orientational ordering in the nematic phase of a polyethylene glycol-peptide conjugate in aqueous solution

The orientational ordering of the nematic phase of a polyethylene glycol (PEG)-peptide block copolymer in aqueous solution is probed by small-angle neutron scattering (SANS), with the sample subjected to steady shear in a Couette cell. The PEG-peptide conjugate forms fibrils that behave as semiflexible rodlike chains. The orientational order parameters (P) over bar (2) and (P) over bar (4) are obtained by modeling the data using a series expansion approach to the form factor of uniform cylinders. The method used is independent of assumptions on the form of the singlet orientational distribution function. Good agreement with the anisotropic two-dimensional SANS patterns is obtained. The results show shear alignment starting at very low shear rates, and the orientational order parameters reach a plateau at higher shear rates with a pseudologarithmic dependence on shear rate. The most probable distribution functions correspond to fibrils parallel to the flow direction under shear, but a sample at rest shows a bimodal distribution with some of the rodlike peptide fibrils oriented perpendicular to the flow direction.

Structural studies on acridine derivatives binding to telomeric DNA

Acridine derivatives can inhibit a variety of nuclear enzymes by binding or intercalating to DNA. This class of compounds is of great interest in the development of novel anticancer agents. Despite the availability of crystallographic data for some of the compounds complexed with DNA, uncertainties remain about the mechanisms of action, binding preferences and biological targets. To investigate the intercalation of several acridine derivatives, a variety of techniques are being employed. Single-crystal X-ray diffraction is being used to determine the high resolution three-dimensional structure of short sequences of quadruplex telomeric DNA with bound drug. This will be compared to the effect of drug binding to long segments of double-stranded DNA using fibre diffraction, with neutron diffraction studies planned to analyse the hydrogen bonding patterns of the DNA-drug complexes. Small-angle neutron scattering (SANS) will also be applied to study drug binding to both short and long sequences of quadruplex and double-stranded DNA in solution. Initial SANS measurements of the telomeric repeat d(TGGGGT) imply that this hexamer is present as a quadruplex. (c) 2006 Elsevier B.V. All rights reserved.

Conformational and hydration effects of site-selective sodium, calcium and strontium ion binding to the DNA Holliday junction structure d(TCGGTACCGA)(4)

The role of metal ions in determining the solution conformation of the Holliday junction is well established, but to date the picture of metal ion binding from structural studies of the four-way DNA junction is very incomplete. Here we present two refined structures of the Holliday junction formed by the sequence d(TCGGTACCGA) in the presence of Na+ and Ca2+, and separately with Sr2+ to resolutions of 1.85 Angstrom and 1.65 Angstrom, respectively. This sequence includes the ACC core found to promote spontaneous junction formation, but its structure has not previously been reported. Almost complete hydration spheres can be defined for each metal cation. The Na+ sites, the most convincing observation of such sites in junctions to date, are one on either face of the junction crossover region, and stabilise the ordered hydration inside the junction arms. The four Ca2+ sites in the same structure are at the CG/CG steps in the minor groove. The Sr2+ ions occupy the TC/AG, GG/CC, and TA/TA sites in the minor groove, giving ten positions forming two spines of ions, spiralling through the minor grooves within each arm of the stacked-X structure. The two structures were solved in the two different C2 lattices previously observed, with the Sr2+ derivative crystallising in the more highly symmetrical form with two-fold symmetry at its centre. Both structures show an opening of the minor groove face of the junction of 8.4degrees in the Ca2+ and Na+ containing structure, and 13.4degrees in the Sr2+ containing structure. The crossover angles at the junction are 39.3degrees and 43.3degrees, respectively. In addition to this, a relative shift in the base pair stack alignment of the arms of 2.3 Angstrom is observed for the Sr2+ containing structure only. Overall these results provide an insight into the so-far elusive stabilising ion structure for the DNA Holliday junction. (C) 2003 Elsevier Science Ltd. All rights reserved.

Horses for courses: Relating innovation diffusion concepts to the stage of the diffusion process

The research uses a sociological perspective to build an improved, context specific understanding of innovation diffusion within the UK construction industry. It is argued there is an iterative interplay between actors and the social system they occupy that directly influences the diffusion process as well as the methodology adopted. The research builds upon previous findings that argued a level of best fit for the three innovation diffusion concepts of cohesion, structural equivalence and thresholds. That level of best fit is analysed here using empirical data from the UK construction industry. This analysis allows an understanding of how the relative importance of these concepts' actually varies within the stages of the innovation diffusion process. The conclusion that the level of relevance fluctuates in relation to the stages of the diffusion process is a new development in the field.

Adapting the International System of Units to the twenty-first century

We review the proposal of the International Committee for Weights and Measures (Comité International des Poids et Mesures, CIPM), currently being considered by the General Conference on Weights and Measures (Conférences Générales des Poids et Mesures, CGPM), to revise the International System of Units (Le Système International d’Unitès, SI). The proposal includes new definitions for four of the seven base units of the SI, and a new form of words to present the definitions of all the units. The objective of the proposed changes is to adopt definitions referenced to constants of nature, taken in the widest sense, so that the definitions may be based on what are believed to be true invariants. In particular, whereas in the current SI the kilogram, ampere, kelvin and mole are linked to exact numerical values of the mass of the international prototype of the kilogram, the magnetic constant (permeability of vacuum), the triple-point temperature of water and the molar mass of carbon-12, respectively, in the new SI these units are linked to exact numerical values of the Planck constant, the elementary charge, the Boltzmann constant and the Avogadro constant, respectively. The new wording used expresses the definitions in a simple and unambiguous manner without the need for the distinction between base and derived units. The importance of relations among the fundamental constants to the definitions, and the importance of establishing a mise en pratique for the realization of each definition, are also discussed.

The key role of the western boundary in linking the AMOC strength to the north-south pressure gradient

A key idea in the study of the Atlantic meridional overturning circulation (AMOC) is that its strength is proportional to the meridional density gradient, or more precisely, to the strength of the meridional pressure gradient. A physical basis that would tell us how to estimate the relevant meridional pressure gradient locally from the density distribution in numerical ocean models to test such an idea, has been lacking however. Recently, studies of ocean energetics have suggested that the AMOC is driven by the release of available potential energy (APE) into kinetic energy (KE), and that such a conversion takes place primarily in the deep western boundary currents. In this paper, we develop an analytical description linking the western boundary current circulation below the interface separating the North Atlantic Deep Water (NADW) and Antarctic Intermediate Water (AAIW) to the shape of this interface. The simple analytical model also shows how available potential energy is converted into kinetic energy at each location, and that the strength of the transport within the western boundary current is proportional to the local meridional pressure gradient at low latitudes. The present results suggest, therefore, that the conversion rate of potential energy may provide the necessary physical basis for linking the strength of the AMOC to the meridional pressure gradient, and that this could be achieved by a detailed study of the APE to KE conversion in the western boundary current.

The primary somatosensory cortex largely contributes to the early part of the cortical response elicited by nociceptive stimuli.

Research on the cortical sources of nociceptive laser-evoked brain potentials (LEPs) began almost two decades ago (Tarkka and Treede, 1993). Whereas there is a large consensus on the sources of the late part of the LEP waveform (N2 and P2 waves), the relative contribution of the primary somatosensory cortex (S1) to the early part of the LEP waveform (N1 wave) is still debated. To address this issue we recorded LEPs elicited by the stimulation of four limbs in a large population (n=35). Early LEP generators were estimated both at single-subject and group level, using three different approaches: distributed source analysis, dipolar source modeling, and probabilistic independent component analysis (ICA). We show that the scalp distribution of the earliest LEP response to hand stimulation was maximal over the central-parietal electrodes contralateral to the stimulated side, while that of the earliest LEP response to foot stimulation was maximal over the central-parietal midline electrodes. Crucially, all three approaches indicated hand and foot S1 areas as generators of the earliest LEP response. Altogether, these findings indicate that the earliest part of the scalp response elicited by a selective nociceptive stimulus is largely explained by activity in the contralateral S1, with negligible contribution from the secondary somatosensory cortex (S2).

Stratospheric water vapor and climate: sensitivity to the representation in radiation codes

There has been considerable interest in the climate impact of trends in stratospheric water vapor (SWV). However, the representation of the radiative properties of water vapor under stratospheric conditions remains poorly constrained across different radiation codes. This study examines the sensitivity of a detailed line-by-line (LBL) code, a Malkmus narrow-band model and two broadband GCM radiation codes to a uniform perturbation in SWV in the longwave spectral region. The choice of sampling rate in wave number space (Δν) in the LBL code is shown to be important for calculations of the instantaneous change in heating rate (ΔQ) and the instantaneous longwave radiative forcing (ΔFtrop). ΔQ varies by up to 50% for values of Δν spanning 5 orders of magnitude, and ΔFtrop varies by up to 10%. In the three less detailed codes, ΔQ differs by up to 45% at 100 hPa and 50% at 1 hPa compared to a LBL calculation. This causes differences of up to 70% in the equilibrium fixed dynamical heating temperature change due to the SWV perturbation. The stratosphere-adjusted radiative forcing differs by up to 96% across the less detailed codes. The results highlight an important source of uncertainty in quantifying and modeling the links between SWV trends and climate.

Factors contributing to the development of extreme North Atlantic cyclones and their relationship with the NAO

The occurrence of extreme cyclones is analysed in terms of their relationship to the NAO phase and the dominating environmental variables controlling their intensification. These are latent energy (equivalent potential temperature 850 hPa is used as an indicator), upper-air baroclinicity, horizontal divergence and jet stream strength. Cyclones over the North Atlantic are identified and tracked using a numerical algorithm, permitting a detailed analysis of their life cycles. Extreme cyclones are selected as the 10% most severe in terms of intensity. Investigations focus on the main strengthening phase of each cyclone. The environmental factors are related to the NAO, which affects the location and orientation of the cyclone tracks, thus explaining why extreme cyclones occur more (less) frequently during strong positive (negative) NAO phases. The enhanced number of extreme cyclones in positive NAO phases can be explained by the larger area with suitable growth conditions, which is better aligned with the cyclone tracks and is associated with increased cyclone life time and intensity. Moreover, strong intensification of cyclones is frequently linked to the occurrence of extreme values of growth factors in the immediate vicinity of the cyclone centre. Similar results are found for ECHAM5/OM1 for present day conditions, demonstrating that relationships between the environment factors and cyclones are also valid in the GCM. For future climate conditions (following the SRES A1B scenario), the results are similar, but a small increase of the frequency of extreme values is detected near the cyclone cores. On the other hand, total cyclone numbers decrease by 10% over the North Atlantic. An exception is the region near the British Isles, which features increased track density and intensity of extreme cyclones irrespective of the NAO phase. These changes are associated with an intensified jet stream close to Europe. Moreover, an enhanced frequency of explosive developments over the British Isles is found, leading to more frequent windstorms affecting Europe.

Satellite observations and modeling of transport in the upper troposphere through the lower mesosphere during the 2006 major stratospheric sudden warming

An unusually strong and prolonged stratospheric sudden warming (SSW) in January 2006 was the first major SSW for which globally distributed long-lived trace gas data are available covering the upper troposphere through the lower mesosphere. We use Aura Microwave Limb Sounder (MLS), Atmospheric Chemistry Experiment-Fourier Transform Spectrometer (ACE-FTS) data, the SLIMCAT Chemistry Transport Model (CTM), and assimilated meteorological analyses to provide a comprehensive picture of transport during this event. The upper tropospheric ridge that triggered the SSW was associated with an elevated tropopause and layering in trace gas profiles in conjunction with stratospheric and tropospheric intrusions. Anomalous poleward transport (with corresponding quasi-isentropic troposphere-to-stratosphere exchange at the lowest levels studied) in the region over the ridge extended well into the lower stratosphere. In the middle and upper stratosphere, the breakdown of the polar vortex transport barrier was seen in a signature of rapid, widespread mixing in trace gases, including CO, H2O, CH4 and N2O. The vortex broke down slightly later and more slowly in the lower than in the middle stratosphere. In the middle and lower stratosphere, small remnants with trace gas values characteristic of the pre-SSW vortex lingered through the weak and slow recovery of the vortex. The upper stratospheric vortex quickly reformed, and, as enhanced diabatic descent set in, CO descended into this strong vortex, echoing the fall vortex development. Trace gas evolution in the SLIMCAT CTM agrees well with that in the satellite trace gas data from the upper troposphere through the middle stratosphere. In the upper stratosphere and lower mesosphere, the SLIMCAT simulation does not capture the strong descent of mesospheric CO and H2O values into the reformed vortex; this poor CTM performance in the upper stratosphere and lower mesosphere results primarily from biases in the diabatic descent in assimilated analyses.