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.

Amino acid pairing preferences in parallel beta-sheets in proteins

Statistical approaches have been applied to examine amino acid pairing preferences within parallel beta-sheets. The main chain hydrogen bonding pattern in parallel beta-sheets means that, for each residue pair, only one of the residues is involved in main chain hydrogen bonding with the strand containing the partner residue. We call this the hydrogen bonded (HB) residue and the partner residue the non-hydrogen bonded (nHB) residue, and differentiate between the favorability of a pair and that of its reverse pair, e.g. Asn(HB)-Thr(nHB)versus Thr(HB)-Asn(nHB). Significantly (p < or = 0.000001) favoured pairings were rationalised using stereochemical arguments. For instance, Asn(HB)-Thr(nHB) and Arg(HB)-Thr(nHB) were favoured pairs, where the residues adopted favoured chi1 rotamer positions that allowed side-chain interactions to occur. In contrast, Thr(HB)-Asn(nHB) and Thr(HB)-Arg(nHB) were not significantly favoured, and could only form side-chain interactions if the residues involved adopted less favourable chi1 conformations. The favourability of hydrophobic pairs e.g. Ile(HB)-Ile(nHB), Val(HB)-Val(nHB) and Leu(HB)-Ile(nHB) was explained by the residues adopting their most preferred chi1 and chi2 conformations, which enabled them to form nested arrangements. Cysteine-cysteine pairs are significantly favoured, although these do not form intrasheet disulphide bridges. Interactions between positively and negatively charged residues were asymmetrically preferred: those with the negatively charged residue at the HB position were more favoured. This trend was accounted for by the presence of general electrostatic interactions, which, based on analysis of distances between charged atoms, were likely to be stronger when the negatively charged residue is the HB partner. The Arg(HB)-Asp(nHB) interaction was an exception to this trend and its favorability was rationalised by the formation of specific side-chain interactions. This research provides rules that could be applied to protein structure prediction, comparative modelling and protein engineering and design. The methods used to analyse the pairing preferences are automated and detailed results are available (http://www.rubic.rdg.ac.uk/betapairprefsparallel/).

Amino acid pairing preferences in parallel beta-sheets in proteins

Statistical approaches have been applied to examine amino acid pairing preferences within parallel beta-sheets. The main chain hydrogen bonding pattern in parallel beta-sheets means that, for each residue pair, only one of the residues is involved in main chain hydrogen bonding with the strand containing the partner residue. We call this the hydrogen bonded (HB) residue and the partner residue the non-hydrogen bonded (nHB) residue, and differentiate between the favourability of a pair and that of its reverse pair, e.g. Asn(HB)-Thr(nHB) versus Thr(HB)-Asn(nHB). Significantly (p <= 0.000001) favoured pairings were rationalised using stereochemical arguments. For instance, Asn(HB)-Thr(nHB) and Arg(HB)-Thr(nHB) were favoured pairs, where the residues adopted favoured chi(1) rotamer positions that allowed side-chain interactions to occur. In contrast, Thr(HB)-Asn(nHB) and Thr(HB)-Arg(nHB) were not significantly favoured, and could only form side-chain interactions if the residues involved adopted less favourable chi(1) conformations. The favourability of hydrophobic pairs e.g. Ile(HB)-Ile(nHB), Val(HB)-Val(nHB) and Leu(HB)-Ile(nHB) was explained by the residues adopting their most preferred chi(1) and chi(2) conformations, which enabled them to form nested arrangements. Cysteine-cysteine pairs are significantly favoured, although these do not form intrasheet disulphide bridges. Interactions between positively and negatively charged residues were asymmetrically preferred: those with the negatively charged residue at the HB position were more favoured. This trend was accounted for by the presence of general electrostatic interactions, which, based on analysis of distances between charged atoms, were likely to be stronger when the negatively charged residue is the HB partner. The Arg(HB)-Asp(nHB) interaction was an exception to this trend and its favourability was rationalised by the formation of specific side-chain interactions. This research provides rules that could be applied to protein structure prediction, comparative modelling and protein engineering and design. The methods used to analyse the pairing preferences are automated and detailed results are available (http:// www.rubic.rdg.ac.uk/betapairprefsparallel/). (c) 2005 Elsevier Ltd. All rights reserved.

A suite of parallel vectors for baculovirus expression

The expression of proteins using recombinant baculoviruses is a mature and widely used technology. However, some aspects of the technology continue to detract from high throughput use and the basis of the final observed expression level is poorly understood. Here, we describe the design and use of a set of vectors developed around a unified cloning strategy that allow parallel expression of target proteins in the baculovirus system as N-terminal or C-terminal fusions. Using several protein kinases as tests we found that amino-terminal fusion to maltose binding protein rescued expression of the poorly expressed human kinase Cot but had only a marginal effect on expression of a well-expressed kinase IKK-2. In addition, MBP fusion proteins were found to be secreted from the expressing cell. Use of a carboxyl-terminal GFP tagging vector showed that fluorescence measurement paralleled expression level and was a convenient readout in the context of insect cell expression, an observation that was further supported with additional non-kinase targets. The expression of the target proteins using the same vectors in vitro showed that differences in expression level were wholly dependent on the environment of the expressing cell and an investigation of the time course of expression showed it could affect substantially the observed expression level for poorly but not well-expressed proteins. Our vector suite approach shows that rapid expression survey can be achieved within the baculovirus system and in addition, goes some way to identifying the underlying basis of the expression level obtained. (c) 2006 Elsevier Inc. All rights reserved.

Rapid parallel expression in E-coli and insect cells: Analysis of five lef gene products of the Autographa californica multiple nuclear polyhedrosis virus (AcMNPV)

A number of strategies are emerging for the high throughput (HTP) expression of recombinant proteins to enable structural and functional study. Here we describe a workable HTP strategy based on parallel protein expression in E. coli and insect cells. Using this system we provide comparative expression data for five proteins derived from the Autographa californica polyhedrosis virus genome that vary in amino acid composition and in molecular weight. Although the proteins are part of a set of factors known to be required for viral late gene expression, the precise function of three of the five, late expression factors (lefs) 6, 7 and 10, is unknown. Rapid expression and characterisation has allowed the determination of their ability to bind DNA and shown a cellular location consistent with their properties. Our data point to the utility of a parallel expression strategy to rapidly obtain workable protein expression levels from many open reading frames (ORFs).

Copper complexes of new benzodioxotetraaza macrocycles with potential applications in nuclear medicine

Two novel benzodioxotetraaza macrocycles [2,9-dioxo-1,4,7,10-tetraazabicyclo[10.4.0]1,11-hexadeca-1(11),13,15-triene (H(2)L1) and 2,10-dioxo-1,4,8,11-tetraazabicyclo[11.4.0]1,12-heptadeca-1(12),14,16-triene (H(2)L2)] were synthesized by a [1 + 1] crablike cyclization. The protonation constants of both ligands were determined by H-1 NMR titration and by potentiometry at 25.0 degrees C in 0.10 M ionic strength in KNO3. The latter method was also used to ascertain the stability constants of their copper(II) complexes. These studies showed that the CuL1 complex has a much lower thermodynamic stability than the CuL2, and the H(2)L2 displays an excellent affinity for copper(II), due to the good fit of copper(II) into its cavity. The copper complexes of the novel ligands were characterized by electronic spectroscopy in solution and by crystal X-ray diffraction. These studies indicated that the copper center in the CuL1 complex adopts a square-pyramidal geometry with the four nitrogen atoms of the macrocycle forming the equatorial plane and a water molecule at axial position, and the copper in the CuL2 complex is square-planar. Several labeling conditions were tested, and only H(2)L2 could be labeled with Cu-67 efficiently (> 98%) in mild conditions (39 degrees C, 15 min) to provide a slightly hydrophilic radioligand (log D = -0.19 +/- 0.03 at pH 7.4). The in vitro stability was studied in the presence of different buffers or with an excess of diethylenetriamine-pentaethanoic acid. Very high stability was shown under these conditions for over 5 days. The incubation of the radiocopper complex in human serum showed 6% protein binding.

Supramolecular parallel beta-sheet and amyloid-like fibril forming peptides using delta-aminovaleric acid residue

Four terminally blocked tripeptides containing delta-aminovaleric acid residue self-assemble to form supramolecular beta-sheet structures as are revealed from their FT-IR data. Single crystal X-ray diffraction studies of two representative peptides also show that they form parallel beta-sheet structures. Self-aggregation of these beta-sheet forming peptides leads to the formation of fibrillar structures, as is evident from scanning electron microscopic (SEM) and transmission electron microscopic (TEM) images. These peptide fibrils bind to a physiological dye, Congo red and exhibit a typical green-gold birefringence under polarized light, showing close resemblance to neurodegenerative disease causing amyloid fibrils. (c) 2005 Elsevier Ltd. All rights reserved.

Studies on the parallel synthesis and evaluation of new heterocyclic extractants for the partitioning of minor actinides

Multiple parallel synthesis and evaluation have been combined in order to identify new nitrogen heterocycles for the partitioning of minor actinides(III) such as americium(III) from lanthanides such as europium(Ill). An array of triazine-containing molecules was made using multiple parallel syntheses from diketones and amide hydrazides. An excess of each of the resulting purified reagents was dissolved in 1,1,2,2-tetrachloroethane containing 2-bromodecanoic acid, and equilibrated with an aqueous solution containing the radiotracers Eu-152 and Am-241 in nitric acid ([Eu] + [Am] < 400 nanomol dm(-3)). Gamma counting of the organic and aqueous phases led to the identification of several new reagents for the selective extraction of americium(III). In particular, 6-(2-pyridyl)-2-(5,6-dialkyl-1,2,4-triazaphenyl)pyridines were found to be effective reagents for the separation of americium(III) from europium(III), (SFAm/Eu was ca. 30 in [HNO3] = 0.013 mol/L).

The synthesis and characterisation of two new iron thioantimonates: [Fe(en)(3)](2)Sb2S5 center dot 0.55H(2)O and [Fe(en)(3)](2)Sb4S8

Two new iron thioantimonates, [Fe(en)(3)](2)Sb2S5 (.) 0.55H(2)O (1) and [Fe(en)(3)](2)Sb4S8 (2). were synthesised under solvothermal conditions from the reactions of Sb2S3, FeCl2 and S in the presence of ethylenediamine at 413 and 438 K, respectively. The products were characterised by single-crystal X-ray diffraction, elemental analysis and SQUID magnetometry. Compound 1 is unusual in containing isolated Sb2S54- anions formed from two corner-sharing SbS33- trigonal pyramids. These units are arranged in rippled layers, 4 A apart, parallel to the bc-plane. Octahedrally coordinated [Fe(en)(3)](2+) cations lie in depressions within these anionic layers. In compound (2), repeated corner linking of SbS33- trigonal pyramids generates SbS2- chains, which may be considered as a polymerised form of the Sb2S54- anions in 1. The SbS2- chains are separated by [Fe(en)(3)](2+) cations. In both compounds, there is an extensive network of hydrogen bonds between the nitrogen atoms of the ethylenediamine ligands and the sulfur atoms of the anions and, in the case of 1, the uncoordinated water molecule. (c) 2005 Elsevier Ltd. All rights reserved.

Methylammonium antimony sulfide

Bis(methylammonium) octaantimony(III) dodecasulfide persulfide, (CH3NH3)(2)[Sb8S12(S-2)], contains pairs of [Sb4S7](2-) chains joined through an unusual persulfide bond to create infinite double [Sb8S14](2-) chains. The double chains are interlocked by longer Sb center dot center dot center dot S interactions to form sheets approximately parallel to the (101) crystallographic plane. Methylammonium cations, formed by decomposition of 2-methylpropane-1,2-diamine during the synthesis, are located in large (Sb8S10) hetero-ring apertures created within the double chains.

New dioxadiaza- and trioxadiaza-macrocycles containing one dibenzofuran unit with two amino pendant arms: synthesis, protonation and complexation studies

New dioxadiaza- and trioxadiaza-macrocycles containing one rigid dibenzofuran unit (DBF) and N-(2-aminoethyl) pendant arms were synthesized, N,N'-bis(2-aminoethyl)-[17]( DBF) N2O2 (L-1) and N,N'-bis(2-aminoethyl)-[22](DBF)N2O3 (L-2), respectively. The binding properties of both macrocycles to metal ions and structural studies of their metal complexes were carried out. The protonation constants of both compounds and the stability constants of their complexes with Co2+, Ni2+, Cu2+, Zn2+, Cd2+, and Pb2+ were determined at 298.2 K, in aqueous solutions, and at ionic strength 0.10 mol dm(-3) in KNO3. Mononuclear complexes with both ligands were formed, and dinuclear complexes were only found for L-2. The thermodynamic binding affinities of the metal complexes of L-2 are lower than those of L-1 as expected, but the Pb2+ complexes of both macrocycles exhibit close stability constant values. On the other hand, the binding affinities of Cd2+ and Pb2+ for L-1 are very high, when compared to those of Co2+, Ni2+ and Zn2+. These interesting properties were explained by the presence of the rigid DBF moiety in the backbone of the macrocycle and to the special match between the macrocyclic cavity size and the studied larger metal ions. To elucidate the adopted structures of complexes in solution, the nickel(II) and copper( II) complexes with both ligands were further studied by UV-vis-MR spectroscopy in DMSO-H2O 1 : 1 (v/v) solution. The copper(II) complexes were also studied by EPR spectroscopy in the same mixture of solvents. The crystal structure of the copper complex of L-1 was also determined. The copper(II) displays an octahedral geometry, the four nitrogen atoms forming the equatorial plane and two oxygen atoms, one from the DBF unit and the other one from the ether oxygen, in axial positions. One of the ether oxygens of the macrocycle is out of the coordination sphere. Our results led us to suggest that this geometry is also adopted by the Co2+ to Zn2+ complexes, and only the larger Cd2+ and Pb2+ manage to form complexes with the involvement of all the oxygen atoms of the macrocyclic backbone.

2D parallel interpenetration of [M-2(bpp)(4)X-4] [M, Fe(II)/Co(II); bpp, 4,4 '-trimethylenedipyridine; X, SCN-SeCN- and N-3(-)] complexes: pseudohalide-dependent conformation of bpp

Three coordination complexes of Co(II)/Fe(II) with 4,4'-trimethylenedipyridine (bpp) and pseudohalides (SCN-, SeCN- and N-3(-)) have been synthesized. The complexes have been characterized by X-ray single crystal structure determination. They are isomorphous having 2D layers in which two independent wavy nets display parallel interwoven structures. Pseudohalide binds metal centers through N terminal and occupies the trans axial positions of the octahedral metal coordination environment. Pseudohalide remains pendant on both sides of the polymeric layer and help the stacking through hydrogen bonding. The conformation of bpp in the interpenetrated nets is observed to be dependent on the choice of pseudohalide. (C) 2008 Elsevier Inc. All rights reserved.

The parallel programming of voluntary and reflexive saccades

A novel two-step paradigm was used to investigate the parallel programming of consecutive, stimulus-elicited ('reflexive') and endogenous ('voluntary') saccades. The mean latency of voluntary saccades, made following the first reflexive saccades in two-step conditions, was significantly reduced compared to that of voluntary saccades made in the single-step control trials. The latency of the first reflexive saccades was modulated by the requirement to make a second saccade: first saccade latency increased when a second voluntary saccade was required in the opposite direction to the first saccade, and decreased when a second saccade was required in the same direction as the first reflexive saccade. A second experiment confirmed the basic effect and also showed that a second reflexive saccade may be programmed in parallel with a first voluntary saccade. The results support the view that voluntary and reflexive saccades can be programmed in parallel on a common motor map. (c) 2006 Elsevier Ltd. All rights reserved.

Disparity with respect to a local reference plane as a dominant cue for stereoscopic depth relief

Earlier studies showed that the disparity with respect to other visible points could not explain stereoacuity performance, nor could various spatial derivatives of disparity [Glennerster, A., McKee, S. P., & Birch, M. D. (2002). Evidence of surface-based processing of binocular disparity. Current Biology, 12:825-828; Petrov, Y., & Glennerster, A. (2004). The role of the local reference in stereoscopic detection of depth relief. Vision Research, 44:367-376.] Two possible cues remain: (i) local changes in disparity gradient or (ii) disparity with respect to an interpolated line drawn through the reference points. Here, we aimed to distinguish between these two cues. Subjects judged.. in a two AFC paradigm, whether a target dot was in front of a plane defined by three reference dots or, in other experiments, in front of a line defined by two reference dots. We tested different slants of the reference line or plane and different locations of the target relative to the reference points. For slanted reference lines or plane, stereoacuity changed little as the target position was varied. For judgments relative to a frontoparallel reference line, stereoacuity did vary with target position, but less than would be predicted by disparity gradient change. This provides evidence that disparity with respect to the reference plane is an important cue. We discuss the potential advantages of this measure in generating a representation of surface relief that is invariant to viewpoint transformations. (c) 2006 Elsevier Ltd. All rights reserved.

Parallel and distributed computing with Java

The Java language first came to public attention in 1995. Within a year, it was being speculated that Java may be a good language for parallel and distributed computing. Its core features, including being objected oriented and platform independence, as well as having built-in network support and threads, has encouraged this view. Today, Java is being used in almost every type of computer-based system, ranging from sensor networks to high performance computing platforms, and from enterprise applications through to complex research-based.simulations. In this paper the key features that make Java a good language for parallel and distributed computing are first discussed. Two Java-based middleware systems, namely MPJ Express, an MPI-like Java messaging system, and Tycho, a wide-area asynchronous messaging framework with an integrated virtual registry are then discussed. The paper concludes by highlighting the advantages of using Java as middleware to support distributed applications.