Rapid model quality assessment for protein structure predictions using the comparison of multiple models without structural alignments

MOTIVATION: The accurate prediction of the quality of 3D models is a key component of successful protein tertiary structure prediction methods. Currently, clustering or consensus based Model Quality Assessment Programs (MQAPs) are the most accurate methods for predicting 3D model quality; however they are often CPU intensive as they carry out multiple structural alignments in order to compare numerous models. In this study, we describe ModFOLDclustQ - a novel MQAP that compares 3D models of proteins without the need for CPU intensive structural alignments by utilising the Q measure for model comparisons. The ModFOLDclustQ method is benchmarked against the top established methods in terms of both accuracy and speed. In addition, the ModFOLDclustQ scores are combined with those from our older ModFOLDclust method to form a new method, ModFOLDclust2, that aims to provide increased prediction accuracy with negligible computational overhead. RESULTS: The ModFOLDclustQ method is competitive with leading clustering based MQAPs for the prediction of global model quality, yet it is up to 150 times faster than the previous version of the ModFOLDclust method at comparing models of small proteins (<60 residues) and over 5 times faster at comparing models of large proteins (>800 residues). Furthermore, a significant improvement in accuracy can be gained over the previous clustering based MQAPs by combining the scores from ModFOLDclustQ and ModFOLDclust to form the new ModFOLDclust2 method, with little impact on the overall time taken for each prediction. AVAILABILITY: The ModFOLDclustQ and ModFOLDclust2 methods are available to download from: http://www.reading.ac.uk/bioinf/downloads/ CONTACT: l.j.mcguffin@reading.ac.uk.

An approximate Bayesian computation approach to overcome biases that arise when using amplified fragment length polymorphism markers to study population structure

There is great interest in using amplified fragment length polymorphism (AFLP) markers because they are inexpensive and easy to produce. It is, therefore, possible to generate a large number of markers that have a wide coverage of species genotnes. Several statistical methods have been proposed to study the genetic structure using AFLP's but they assume Hardy-Weinberg equilibrium and do not estimate the inbreeding coefficient, F-IS. A Bayesian method has been proposed by Holsinger and colleagues that relaxes these simplifying assumptions but we have identified two sources of bias that can influence estimates based on these markers: (i) the use of a uniform prior on ancestral allele frequencies and (ii) the ascertainment bias of AFLP markers. We present a new Bayesian method that avoids these biases by using an implementation based on the approximate Bayesian computation (ABC) algorithm. This new method estimates population-specific F-IS and F-ST values and offers users the possibility of taking into account the criteria for selecting the markers that are used in the analyses. The software is available at our web site (http://www-leca.uif-grenoble.fi-/logiciels.htm). Finally, we provide advice on how to avoid the effects of ascertainment bias.

Origin of avian genome size and structure in non-avian dinosaurs

Avian genomes are small and streamlined compared with those of other amniotes by virtue of having fewer repetitive elements and less non-coding DNA(1,2). This condition has been suggested to represent a key adaptation for flight in birds, by reducing the metabolic costs associated with having large genome and cell sizes(3,4). However, the evolution of genome architecture in birds, or any other lineage, is difficult to study because genomic information is often absent for long-extinct relatives. Here we use a novel bayesian comparative method to show that bone-cell size correlates well with genome size in extant vertebrates, and hence use this relationship to estimate the genome sizes of 31 species of extinct dinosaur, including several species of extinct birds. Our results indicate that the small genomes typically associated with avian flight evolved in the saurischian dinosaur lineage between 230 and 250 million years ago, long before this lineage gave rise to the first birds. By comparison, ornithischian dinosaurs are inferred to have had much larger genomes, which were probably typical for ancestral Dinosauria. Using comparative genomic data, we estimate that genome-wide interspersed mobile elements, a class of repetitive DNA, comprised 5 - 12% of the total genome size in the saurischian dinosaur lineage, but was 7 - 19% of total genome size in ornithischian dinosaurs, suggesting that repetitive elements became less active in the saurischian lineage. These genomic characteristics should be added to the list of attributes previously considered avian but now thought to have arisen in non-avian dinosaurs, such as feathers(5), pulmonary innovations 6, and parental care and nesting

Heterogenisation of Os species on MCM-41 structure for epoxidation of trans-stilbene

Metal organic chemical vapour deposition technique (MOCVD) has been used to immobilise Os species onto the internal porous structure of MCM-41. Evidence suggests that volatile Os-3(CO)(12) cluster reacts with surface silanol groups of the MCM-41 via an oxidative addition reaction to yield a trinuclear HOs3(CO)(10)(OSi-) surface species. After heat treatment in air or at their very low surface coverage, these triangular sites break up to partially oxidised mononuclear surface species. In the presence of tert-butyl hydroperoxide (TBHP) as an oxidant, we demonstrate that the mononuclear species form extremely active species that catalyse the oxidation of trans-stilbene selectively to the corresponding epoxide. By carefully controlling the parameters of the MOCVD method (loading and calcination temperature), we report a new class of optimised MCM-41 porous heterogeneous catalysts carrying isolated but active Os sites for the selective oxidation of trans-stilbene in liquid phase. The reaction selectivity of the solid supported Os is apparently higher than the soluble homogeneous Os-3(CO)(12) cluster. It is envisaged that our solid supported catalysts not only facilitate separation from products but also offer an excellent utilisation of Os for catalysis. (C) 2003 Elsevier Science B.V. All rights reserved.

Macrocyclic amines as structure-directing agents for the synthesis of three-dimensional antimony-sulfide frameworks

A new family of antimony sulfides, incorporating the macrocyclic tetramine 1,4,8,11-tetraazacyclotetradecane ( cyclam), has been prepared by a hydrothermal method. [C10N4H26][Sb4S7] (1), [Ni(C10N4H24)][Sb4S7] (2), and [Co(C10N4H24)](x)[C10N4H26](1-x)[Sb4S7] (0.08 <= x <= 0.74) (3) have been characterized by single-crystal X-ray diffraction, elemental analysis, thermogravimetry, and analytical electron microscopy. All three materials possess the same novel three-dimensional Sb4S72- framework, constructed from layers of parallel arrays of Sb4S84- chains stacked at 90 to one another. In 1, doubly protonated macrocyclic cations reside in the channel structure of the antimonysulfide framework. In 2 and 3, the cyclam acts as a ligand, chelating the divalent transition- metal cation. Analytical and X-ray diffraction data indicate that the level of metal incorporation in 2 is effectively complete, whereas in 3, both metalated and nonmetalated forms of the macrocycle coexist within the structure.

Simple general method of generating non-oxo, non-amavadine model octacoordinated vanadium(IV) complexes of some tetradentate ONNO chelating ligands from various oxovanadium(IV/V) compounds and structural characterization of one of them

A simple general route of obtaining very stable octacoordinated non-oxovanadium( IV) complexes of the general formula VL2 (where H2L is a tetradentate ONNO donor) is presented. Six such complexes (1-6) are adequately characterized by elemental analysis, mass spectrometry, and various spectroscopic techniques. One of these compounds (1) has been structurally characterized. The molecule has crystallographic 4 symmetry and has a dodecahedral structure existing in a tetragonal space group P4n2. The non-oxo character and VL2 stoichiometry for all of the complexes are established from analytical and mass spectrometric data. In addition, the non-oxo character is clearly indicated by the complete absence of the strong nu(v=o) band in the 925-1025 cm(-1) region, which is a signature of all oxovanadium species. The complexes are quite stable in open air in the solid state and in solution, a phenomenon rarely observed in non-oxovanadium(IV) or bare vanadium(IV) complexes.

Influence of disaccharide structure on prebiotic selectivity in vitro

To obtain structure-function information of a range of carbohydrates, which are available only in very small quantities, an in vitro fermentation method using 7 mg of carbohydrate, 0.7 mL of basal medium, and 1% (w/v) of fecal bacteria was validated against a pH-controlled batch culture with 150 mL of basal medium and 1.5g of test carbohydrate. This method was used to determine the influence of different glycosidic linkages and monosaccharide compositions of disaccharides on the selectivity of microbial fermentation. A prebiotic index (PI) was calculated for each disaccharide. Generally, disaccharides with linkages of 1-2, 1-4, and 1-6 generated a high PI score, with kojibiose and sophorose showing the greatest values (21.62 and 18.63, respectively). Apart from 6 alpha-mannobiose, mannose-containing disaccharicles gave a low PI due to low numbers of bifidobacteria and lactobacilli and an increase in bacteroides. The structure-function information obtained in this study may lead to a predictive understanding of how specific structures are fermented by the human gut microflora.

Training conditional and cumulative risk judgements: the role of frequencies, problem-structure and Einstellung

Two experiments implement and evaluate a training scheme for learning to apply frequency formats to probability judgements couched in terms of percentages. Results indicate that both conditional and cumulative probability judgements can be improved in this manner, however the scheme is insufficient to promote any deeper understanding of the problem structure. In both experiments, training on one problem type only (either conditional or cumulative risk judgements) resulted in an inappropriate transfer of a learned method at test. The obstacles facing a frequency-based training programme for teaching appropriate use of probability data are discussed. Copyright (c) 2006 John Wiley & Sons, Ltd.

The factor structure of the Edinburgh postnatal depression scale in a South African peri-urban settlement

The factor structure of the Edinburgh Postnatal Depression scale (EPDS) and similar instruments have received little attention in the literature. The researchers set out to investigate the construct validity and reliability of the EPDS amongst impoverished South African women. The EPDS was translated into isiXhosa (using Brislin's back translation method) and administered by trained interviewers to 147 women in Khayelitsha, South Africa. Responses were subjected to maximum likelihood confirmatory factor analysis. A single factor structure was found, consistent with the theory on which the EPDS was based. Internal consistency was satisfactory (a = 0.89).

A new method for evaluating regional air quality forecasts

A Kriging interpolation method is combined with an object-based evaluation measure to assess the ability of the UK Met Office's dispersion and weather prediction models to predict the evolution of a plume of tracer as it was transported across Europe. The object-based evaluation method, SAL, considers aspects of the Structure, Amplitude and Location of the pollutant field. The SAL method is able to quantify errors in the predicted size and shape of the pollutant plume, through the structure component, the over- or under-prediction of the pollutant concentrations, through the amplitude component, and the position of the pollutant plume, through the location component. The quantitative results of the SAL evaluation are similar for both models and close to a subjective visual inspection of the predictions. A negative structure component for both models, throughout the entire 60 hour plume dispersion simulation, indicates that the modelled plumes are too small and/or too peaked compared to the observed plume at all times. The amplitude component for both models is strongly positive at the start of the simulation, indicating that surface concentrations are over-predicted by both models for the first 24 hours, but modelled concentrations are within a factor of 2 of the observations at later times. Finally, for both models, the location component is small for the first 48 hours after the start of the tracer release, indicating that the modelled plumes are situated close to the observed plume early on in the simulation, but this plume location error grows at later times. The SAL methodology has also been used to identify differences in the transport of pollution in the dispersion and weather prediction models. The convection scheme in the weather prediction model is found to transport more pollution vertically out of the boundary layer into the free troposphere than the dispersion model convection scheme resulting in lower pollutant concentrations near the surface and hence a better forecast for this case study.

Gas-phase rate coefficients for the reactions of O(3P), S(3P), Se(3P), and Te(3P) with alkenes: application of perturbation frontier molecular orbital theory, correlations, and structure-activity relations (SARs)

The kinetics of the reactions of the atoms O(P-3), S(P-3), Se(P-3), and Te((3)p) with a series of alkenes are examined for correlations relating the logarithms of the rate coefficients to the energies of the highest occupied molecular orbitals (HOMOs) of the alkenes. These correlations may be employed to predict rate coefficients from the calculated HOMO energy of any other alkene of interest. The rate coefficients obtained from the correlations were used to formulate structure-activity relations (SARs) for reactions of O((3)p), S(P-3), Se (P-3), and Te((3)p) with alkenes. A comparison of the values predicted by both the correlations and the SARs with experimental data where they exist allowed us to assess the reliability of our method. We demonstrate the applicability of perturbation frontier molecular orbital theory to gas-phase reactions of these atoms with alkenes. The correlations are apparently not applicable to reactions of C(P-3), Si(P-3), N(S-4), and Al(P-2) atoms with alkenes, a conclusion that could be explained in terms of a different mechanism for reaction of these atoms.

Algorithmic statistical analysis of electrophysiological data for the investigation of structure-activity relationship in single neurons

We are developing computational tools supporting the detailed analysis of the dependence of neural electrophysiological response on dendritic morphology. We approach this problem by combining simulations of faithful models of neurons (experimental real life morphological data with known models of channel kinetics) with algorithmic extraction of morphological and physiological parameters and statistical analysis. In this paper, we present the novel method for an automatic recognition of spike trains in voltage traces, which eliminates the need for human intervention. This enables classification of waveforms with consistent criteria across all the analyzed traces and so it amounts to reduction of the noise in the data. This method allows for an automatic extraction of relevant physiological parameters necessary for further statistical analysis. In order to illustrate the usefulness of this procedure to analyze voltage traces, we characterized the influence of the somatic current injection level on several electrophysiological parameters in a set of modeled neurons. This application suggests that such an algorithmic processing of physiological data extracts parameters in a suitable form for further investigation of structure-activity relationship in single neurons.

Spectroscopic properties of alkenylferrocenes; crystal and molecular structure of (Z)-(1,2-diphenylethenyl)ferrocene

Studies of the 1H n.m.r. and electronic spectra of a series of alkenylferrocenes including (E) and (Z) stereoisomers of various styrylferrocenes, have provided methods of structure elucidation. Crystals of the title compound are monoclinic, space group P21/c with Z= 4 in a unit cell of dimensions a= 17.603(2), b= 10.218(2), c= 10.072 Å, β= 103.27(2)°. The structure has been determined by the heavy-atom method from diffractometer data and refind by full-matrix least-squares techniques to R= 0.043 for 2 219 unique reflections.

Unsaturated σ-hydrocarbyl transition-metal complexes. Part 4. Crystal and molecular structure of trans-chlorobis(diethylphenylphosphine)-(phenylethynyl)platinum(II) and comments on the relative trans influence of various carbon ligands

The molecular structure of trans-[PtCl(CCPh)(PEt2Ph)2] has been determined by X-ray diffraction methods. The crystals are monoclinic, space group P21, with a= 12.359(3), b= 13.015(3), c= 9.031(2)Å, β= 101.65(2)°, and Z= 2. The structure has been solved by the heavy-atom method and refined by full-matrix least squares to R 0.046 for 1 877 diffractometric intensity data. The crystals contain discrete molecules in which the platinum coordination is square planar. The phenylethynyl group is non-linear, with a Pt–CC angle of 163(2)°. Selected bond lengths are Pt–Cl 2.407(5) and Pt–C 1.98(2)Å. The structural trans influences of CCPh, CHCH2, and CH2SiMe3 ligands in platinum(II) complexes are compared; there is only a small dependence on hybridization at the ligating carbon atom.

Unsaturated σ-hydrocarbyl transition-metal complexes. Part 3. Crystal and molecular structure of trans-chlorobis(diethylphenylphosphine)(vinyl)platinum(II)

The molecular structure of trans-[PtCl(CHCH2)(PEt2Ph)2] has been determined by X-ray diffraction methods. The crystals are orthorhombic, space group Pbcn, with a= 10.686(2), b= 13.832(4), c= 16.129(4)Å, and Z= 4. The structure has been solved by the heavy-atom method and refined by full-matrix least squares to R 0.044 for 1 420 diffractometric intensity data. The crystals contain discrete molecules in which the platinum co-ordination is square planar. The Pt–Cl bond vector coincides with a crystallographic diad axis about which the atoms of the vinyl group are disordered. Selected bond lengths (Å) are Pt–Cl 2.398(4), Pt–P 2.295(3), and Pt–C 2.03(2). The Pt–CC angle is 127(2)°. From a survey of the available structural data it is concluded that there is little, if any, back donation from platinum to carbon in platinum–alkenyl linkages.