Mitochondrial control region diversity of the houbara bustard Chlamydotis undulata complex and genetic structure along the Atlantic seaboard of North Africa

The houbara bustard, Chlamydotis undulata, is a declining cryptic desert bird whose range extends from North Africa to Central Asia. Three subspecies are currently recognized by geographical distribution and morphology: C.u.fuertaventurae, C.u.undulata and C.u.macqueenii. We have sequenced 854 bp of mitochondrial control region from 73 birds to describe their population genetic structure with a particular sampling focus on the connectivity between C.u.fuertaventurae and C.u.undulata along the Atlantic seaboard of North Africa. Nucleotide and haplotypic diversity varied among the subspecies being highest in C.u.undulata, lowest in C.u.fuertaventurae and intermediate in C.u.macqueenii. C.u.fuertaventurae and C.u.undulata are paraphyletic and an average nucleotide divergence of 2.08% splits the later from C.u.macqueenii. We estimate that C.u.fuertaventurae and C.u.undulata split from C.u.macqueenii approximately 430 000 years ago. C.u.fuertaventurae and C.u.undulata are weakly differentiated (F-ST = 0.27, N-m = 1.3), indicative of a recent shared history. Archaeological evidence indicates that houbara bustards have been present on the Canary Islands for 130-170 000 years. However, our genetic data point to a more recent separation of C.u.fuertaventurae and C.u.undulata at around 20-25 000 years. Concordant archaeological, climatic opportunities for colonization and genetic data point to a scenario of: (i) initial colonization of the Canary Islands about 130 000 years ago; (ii) a period of secondary contact 19-30 000 years ago homogenizing any pre-existing genetic structure followed by; (iii) a period of relative isolation that persists today.

Structural basis for control of secondary vessels in the long-finned eel, Anguilla reinhardtii

Histological sections of primary segmental arteries and associated interarterial anastomoses and secondary vessels from the long-finned eel Anguilla reinhardtii were examined by light and transmission electron microscopy. Interarterial anastomoses were found to originate from the primary vasculature as depressions through the tunica intima and media, from where they ran perpendicularly to the adventitial layer, before coiling extensively. From here the anastomoses travelled a relatively linear path in the outer margin of the adventitia to anastomose with a secondary vessel running in parallel with the primary counterpart. In contrast to findings from other species, secondary vessels had a structure quite similar to that of primary vessels; they were lined by endothelial cells on a continuous basement membrane, with a single layer of smooth muscle cells surrounding the vessel. Smooth muscle cells were also found in the vicinity of interarterial anastomoses in the adventitia, but these appeared more longitudinally orientated. The presence of smooth muscle cells on all aspects of the secondary circulation suggests that this vascular system is regulated in a similar manner as the primary vascular system. Because interarterial anastomoses are structurally integrated with the primary vessel from which they originate, it is anticipated that flow through secondary vessels to some extent is affected by the vascular tone of the primary vessel. Immunohistochemical studies showed that primary segmental arteries displayed moderate immunoreactivity to antibodies against 5-hydroxytryptamine and substance P, while interarterial anastomoses and secondary vessels showed dense immunoreactivity. No immunoreactivity was observed on primary or secondary arteries against neuropeptide Y or calcitonin gene-related peptide.

Accurate prediction of scorpion toxin functional properties from primary structures

Scorpion toxins are common experimental tools for studies of biochemical and pharmacological properties of ion channels. The number of functionally annotated scorpion toxins is steadily growing, but the number of identified toxin sequences is increasing at much faster pace. With an estimated 100,000 different variants, bioinformatic analysis of scorpion toxins is becoming a necessary tool for their systematic functional analysis. Here, we report a bioinformatics-driven system involving scorpion toxin structural classification, functional annotation, database technology, sequence comparison, nearest neighbour analysis, and decision rules which produces highly accurate predictions of scorpion toxin functional properties. (c) 2005 Elsevier Inc. All rights reserved.

Substrate specifity of protein kinases and computational prediction of substrates

To ensure signalling fidelity, kinases must act only on a defined subset of cellular targets. Appreciating the basis for this substrate specificity is essential for understanding the role of an individual protein kinase in a particular cellular process. The specificity in the cell is determined by a combination of peptide specificity of the kinase (the molecular recognition of the sequence surrounding the phosphorylation site), substrate recruitment and phosphatase activity. Peptide specificity plays a crucial role and depends on the complementarity between the kinase and the substrate and therefore on their three-dimensional structures. Methods for experimental identification of kinase substrates and characterization of specificity are expensive and laborious, therefore, computational approaches are being developed to reduce the amount of experimental work required in substrate identification. We discuss the structural basis of substrate specificity of protein kinases and review the experimental and computational methods used to obtain specificity information. (c) 2005 Elsevier B.V. All rights reserved.

Multiagent realization of prediction-based diagnosis and loss prevention

A multiagent diagnostic system implemented in a Protege-JADE-JESS environment interfaced with a dynamic simulator and database services is described in this paper. The proposed system architecture enables the use of a combination of diagnostic methods from heterogeneous knowledge sources. The process ontology and the process agents are designed based on the structure of the process system, while the diagnostic agents implement the applied diagnostic methods. A specific completeness coordinator agent is implemented to coordinate the diagnostic agents based on different methods. The system is demonstrated on a case study for diagnosis of faults in a granulation process based on HAZOP and FMEA analysis.

Prediction of HLA-DQ3.2 beta ligands: evidence of multiple registers in class II binding peptides

Motivation: While processing of MHC class II antigens for presentation to helper T-cells is essential for normal immune response, it is also implicated in the pathogenesis of autoimmune disorders and hypersensitivity reactions. Sequence-based computational techniques for predicting HLA-DQ binding peptides have encountered limited success, with few prediction techniques developed using three-dimensional models. Methods: We describe a structure-based prediction model for modeling peptide-DQ3.2 beta complexes. We have developed a rapid and accurate protocol for docking candidate peptides into the DQ3.2 beta receptor and a scoring function to discriminate binders from the background. The scoring function was rigorously trained, tested and validated using experimentally verified DQ3.2 beta binding and non-binding peptides obtained from biochemical and functional studies. Results: Our model predicts DQ3.2 beta binding peptides with high accuracy [area under the receiver operating characteristic (ROC) curve A(ROC) > 0.90], compared with experimental data. We investigated the binding patterns of DQ3.2 beta peptides and illustrate that several registers exist within a candidate binding peptide. Further analysis reveals that peptides with multiple registers occur predominantly for high-affinity binders.

SCORPION2: A database for structure-function analysis of scorpion toxins

Scorpion toxins are important experimental tools for characterization of vast array of ion channels and serve as scaffolds for drug design. General public database entries contain limited annotation whereby rich structure-function information from mutation studies is typically not available. SCORPION2 contains more than 800 records of native and mutant toxin sequences enriched with binding affinity and toxicity information, 624 three-dimensional structures and some 500 references. SCORPION2 has a set of search and prediction tools that allow users to extract and perform specific queries: text searches of scorpion toxin records, sequence similarity search, extraction of sequences, visualization of scorpion toxin structures, analysis of toxic activity, and functional annotation of previously uncharacterized scorpion toxins. The SCORPION2 database is available at http://sdmc.i2r.a-star.edu.sg/scorpion/. (c) 2006 Elsevier Ltd. All rights reserved.

Solution structure and characterization of the LGR8 receptor binding surface of insulin-like peptide 3

Insulin-like peptide 3 (INSL3), a member of the relaxin peptide family, is produced in testicular Leydig cells and ovarian thecal cells. Gene knock-out experiments have identified a key biological role in initiating testes descent during fetal development. Additionally, INSL3 has an important function in mediating male and female germ cell function. These actions are elicited via its recently identified receptor, LGR8, a member of the leucine-rich repeat-containing G-protein- coupled receptor family. To identify the structural features that are responsible for the interaction of INSL3 with its receptor, its solution structure was determined by NMR spectroscopy together with in vitro assays of a series of B-chain alanine-substituted analogs. Synthetic human INSL3 was found to adopt a characteristic relaxin/ insulin-like fold in solution but is a highly dynamic molecule. The four termini of this two-chain peptide are disordered, and additional conformational exchange is evident in the molecular core. Alanine-substituted analogs were used to identify the key residues of INSL3 that are responsible for the interaction with the ectodomain of LGR8. These include Arg(B16) and Val(B19), with His(B12) and Arg(B20) playing a secondary role, as evident from the synergistic effect on the activity in double and triple mutants involving these residues. Together, these amino acids combine with the previously identified critical residue, Trp(B27), to form the receptor binding surface. The current results provide clear direction for the design of novel specific agonists and antagonists of this receptor.

The internal cranial anatomy of the Plesiosauria (Reptilia, Sauropterygia): evidence for a functional secondary palate

In the late 19th Century, the choanae (or internal nares) of the Plesiosauria were identified as a pair of palatal openings located rostral to the external nares, implying a rostrally directed respiratory duct and air path inside the rostrum. Despite obvious functional shortcomings, this idea was firmly established in the scientific literature by the first decade of the 20th Century. The functional consequences of this morphology were only re-examined by the end of the 20th Century, leading to the conclusion that the choanae were not involved in respiration but instead in underwater olfaction, the animals supposedly breathing with the mouth agape. Re-evaluation of the palatal and internal cranial anatomy of the Plesiosauria reveals that the traditional identification of the choanae as a pair of fenestrae situated rostral to the external nares appears erroneous. These openings more likely represent the bony apertures of ducts that lead to internal salt glands situated inside the maxillary rostrum. The 'real' functional choanae (or caudal interpterygoid vacuities), are situated at the caudal end of the bony palate between the sub-temporal fossae, as was suggested in the mid-19th Century. The existence of a functional secondary palate in the Plesiosauria is therefore strongly supported, and the anatomical, physiological, and evolutionary implications of such a structure are discussed.

Vibrational dynamics and structure of natural feldspars

The present study gives a contribution to the knowledge on the Na-feldspar and plagioclases, extending the database of the Raman spectra of plagioclases with different chemical compositions and structural orders. This information may be used for the future planetary explorations by “rovers”, for the investigation of ceramics nanocrystal materials and for the mineralogical phase identification in sediments. Na-feldspar and plagioclase solid solution have been investigated by Raman spectroscopy in order to determine the relationships between the vibrational changes and the plagioclase crystal chemistry and structure. We focused on the Raman micro-spectroscopy technique, being a non-destructive method, suited for contactless analysis with high spatial resolution. Chemical and structural analyses have been performed on natural samples to test the usefulness of Raman spectroscopy as a tool in the study of the pressure-induced structural deformations, the disordering processes due to change in the Al-Si distribution in the tetrahedral sites and, finally, in the determination of the anorthitic content (Anx) in plagioclase minerals. All the predicted 39 Ag Raman active modes have been identified and assigned to specific patterns of atomic vibrational motion. A detailed comparison between experimental and computed Raman spectra has been performed and previous assignments have been revised, solving some discrepancies reported in recent literature. The ab initio calculation at the hybrid HF/DFT level with the WC1LYP Hamiltonian has proven to give excellent agreement between calculated and experimentally measured Raman wavenumbers and intensities in triclinic minerals. A short digression on the 36 infrared active modes of Na-feldspar has been done too. The identification of all 39 computed Raman modes in the experimentally measured spectra of the fully ordered Na-feldspar, known as low albite, along with the detailed description of each vibrational mode, has been essential to extend the comparative analysis to the high pressure and high temperature structural forms of albite, which reflect the physical–chemical conditions of the hosting rocks. The understanding of feldspar structure response to pressure and temperature is crucial in order to constrain crustal behaviour. The compressional behaviour of the Na-feldspar has been investigated for the first time by Raman spectroscopy. The absence of phase transitions and the occurrence of two secondary compression mechanisms acting at different pressures have been confirmed. Moreover, Raman data suggest that the internal structural changes are confined to a small pressure interval, localized around 6 GPa, not spread out from 4 to 8 GPa as suggested by previous X-rays studies on elasticity. The dominant compression mechanisms act via tetrahedral tilting, while the T-O bond lengths remain nearly constant at moderate compressional regimes. At the spectroscopic level, this leads to the strong pressure dependencies of T-O-T bending modes, as found for the four modes at 478, 508, 578 and 815 cm-1. The Al-Si distribution in the tetrahedral sites affects also the Raman spectrum of Na-feldspar. In particular, peak broadening is more sensitive than peak position to changes in the degree of order. Raman spectroscopy is found to be a good probe for local ordering, in particular being sensitive to the first annealing steps, when the macroscopic order parameter is still high. Even though Raman data are scattered and there are outliers in the estimated values of the degree of order, the average peak linewidths of the Na-feldspar characteristic doublet band, labelled here as υa and υb, as a function of the order parameter Qod show interesting trends: both peak linewidths linearly increase until saturation. From Qod values lower than 0.6, peak broadening is no more affected by the Al-Si distribution. Moreover, the disordering process is found to be heterogeneous. SC-XRD and Raman data have suggested an inter-crystalline inhomogeneity of the samples, i.e., the presence of regions with different defect density on the micrometric scale. Finally, the influence of Ca-Na substitution in the plagioclase Raman spectra has been investigated. Raman spectra have been collected on a series of well characterized natural, low structural plagioclases. The variations of the Raman modes as a function of the chemical composition and the structural order have been determined. The number of the observed Raman bands at each composition gives information about the unit-cell symmetry: moving away from the C1 structures, the number of the Raman bands enhances, as the number of formula units in the unit cell increases. The modification from an “albite-like” Raman spectrum to a more “anorthite-like” spectrum occurs from sample An78 onwards, which coincides with the appearance of c reflections in the diffraction patterns of the samples. The evolution of the Raman bands υa and υb displays two changes in slope at ~An45 and ~An75: the first one occurs between e2 and e1 plagioclases, the latter separates e1 and I1 plagioclases with only b reflections in their diffraction patterns from I1 and P1 samples having b and c reflections too. The first variation represents exactly the e2→e1 phase transitions, whereas the second one corresponds in good approximation to the C1→I1 transition, which has been determined at ~An70 by previous works. The I1→P1 phase transition in the anorthite-rich side of the solid solution is not highlighted in the collected Raman spectra. Variations in peak broadening provide insights into the behaviour of the order parameter on a local scale, suggesting an increase in the structural disorder within the solid solution, as the structures have to incorporate more Al atoms to balance the change from monovalent to divalent cations. All the information acquired on these natural plagioclases has been used to produce a protocol able to give a preliminary estimation of the chemical composition of an unknown plagioclase from its Raman spectrum. Two calibration curves, one for albite-rich plagioclases and the other one for the anorthite-rich plagioclases, have been proposed by relating the peak linewidth of the most intense Raman band υa and the An content. It has been pointed out that the dependence of the composition from the linewidth can be obtained only for low structural plagioclases with a degree of order not far away from the references. The proposed tool has been tested on three mineralogical samples, two of meteoric origin and one of volcanic origin. Chemical compositions by Raman spectroscopy compare well, within an error of about 10%, with those obtained by elemental techniques. Further analyses on plagioclases with unknown composition will be necessary to validate the suggested method and introduce it as routine tool for the determination of the chemical composition from Raman data in planetary missions.

Coalescence of secondary dispersions

A study has been made of the coalescence of secondary dispersions in beds of woven meshes. The variables investigated were superficial velocity, bed depth, mesh geometry and fibre material; the effects of presoaking the bed in the dispersed phase before operation were also considered. Equipment was design~d to generate a 0.1% phase ratio toluene in water dispersion whose mean drop size was determined using a Coulter Counter. The coalesced drops were sized by photography and a novel holographic technique was developed to evaluate the mean diameter of the effluent secondary drops. Previous models describing single phase flow in porous media are reviewed and it was found that the experimental data obtained in this study is best represented by Keller's equation which is based on a physical model similar to the internal structure of the meshes. Statistical analysis of two phase data produced a correlation, for each mesh tested, relating the pressure drop to superficial velocity and bed depth. The flow parameter evaluated from the single phase model is incorporated into a theoretical comparison of drop capture mechanisms which indicated that direct and indirect interception are predominant. The resulting equation for drop capture efficiericy is used to predict the initial, local drop capture rate in a coalescer. A mathematical description of the saturation profiles was formulated and verified by average saturation data. Based 6n the Blake-Kozeny equation, an expression is derived analytically to predict the two phase pressure drop using the parameters which characterise the saturation profiles. By specifying the local saturation at the inlet face for a given velocity, good agreement between experimental pressure drop data and the model predictions was obtained.

HLA-DP2 binding prediction by molecular dynamics simulations

Major histocompatibility complex (MHC) II proteins bind peptide fragments derived from pathogen antigens and present them at the cell surface for recognition by T cells. MHC proteins are divided into Class I and Class II. Human MHC Class II alleles are grouped into three loci: HLA-DP, HLA-DQ, and HLA-DR. They are involved in many autoimmune diseases. In contrast to HLA-DR and HLA-DQ proteins, the X-ray structure of the HLA-DP2 protein has been solved quite recently. In this study, we have used structure-based molecular dynamics simulation to derive a tool for rapid and accurate virtual screening for the prediction of HLA-DP2-peptide binding. A combinatorial library of 247 peptides was built using the "single amino acid substitution" approach and docked into the HLA-DP2 binding site. The complexes were simulated for 1 ns and the short range interaction energies (Lennard-Jones and Coulumb) were used as binding scores after normalization. The normalized values were collected into quantitative matrices (QMs) and their predictive abilities were validated on a large external test set. The validation shows that the best performing QM consisted of Lennard-Jones energies normalized over all positions for anchor residues only plus cross terms between anchor-residues.

AIDS: structure, replication and spread of the AIDS virus

AIDS (Acquired Immune Deficiency Syndrome)was first, described as a new disease of humans in 1981. The origins of the disease are controversial. AIDS is caused by a retrovirus, a type of virus which rarely attacks human cells. The first virus of this type recorded in humans is reponsible for a type of leukaemia and was identified in 1978. AIDS is thus the third type of human retrovirus to be discovered and hence, is referred to as T-lymphotrophic virus III (HTLV-III). For viruses to replicate, they have to invade a host cell which in this case is a T4-lymphocyte, a type of white blood cell that regulates the immune system. The problems of the disease result directly from the death of these cells. As a consequence, the immune system is compromised leading to a number of opportunistic secondary infections and other disorders.

Peptide binding prediction for the human class II MHC allele HLA-DP2:a molecular docking approach

MHC class II proteins bind oligopeptide fragments derived from proteolysis of pathogen antigens, presenting them at the cell surface for recognition by CD4+ T cells. Human MHC class II alleles are grouped into three loci: HLA-DP, HLA-DQ and HLA-DR. In contrast to HLA-DR and HLA-DQ, HLA-DP proteins have not been studied extensively, as they have been viewed as less important in immune responses than DRs and DQs. However, it is now known that HLA-DP alleles are associated with many autoimmune diseases. Quite recently, the X-ray structure of the HLA-DP2 molecule (DPA*0103, DPB1*0201) in complex with a self-peptide derived from the HLA-DR a-chain has been determined. In the present study, we applied a validated molecular docking protocol to a library of 247 modelled peptide-DP2 complexes, seeking to assess the contribution made by each of the 20 naturally occurred amino acids at each of the nine binding core peptide positions and the four flanking residues (two on both sides).

Using secondary knowledge to support decision tree classification of retrospective clinical data

Retrospective clinical data presents many challenges for data mining and machine learning. The transcription of patient records from paper charts and subsequent manipulation of data often results in high volumes of noise as well as a loss of other important information. In addition, such datasets often fail to represent expert medical knowledge and reasoning in any explicit manner. In this research we describe applying data mining methods to retrospective clinical data to build a prediction model for asthma exacerbation severity for pediatric patients in the emergency department. Difficulties in building such a model forced us to investigate alternative strategies for analyzing and processing retrospective data. This paper describes this process together with an approach to mining retrospective clinical data by incorporating formalized external expert knowledge (secondary knowledge sources) into the classification task. This knowledge is used to partition the data into a number of coherent sets, where each set is explicitly described in terms of the secondary knowledge source. Instances from each set are then classified in a manner appropriate for the characteristics of the particular set. We present our methodology and outline a set of experiential results that demonstrate some advantages and some limitations of our approach. © 2008 Springer-Verlag Berlin Heidelberg.