Protein contact prediction using patterns of correlation

We describe a new method for using neural networks to predict residue contact pairs in a protein. The main inputs to the neural network are a set of 25 measures of correlated mutation between all pairs of residues in two windows of size 5 centered on the residues of interest. While the individual pair-wise correlations are a relatively weak predictor of contact, by training the network on windows of correlation the accuracy of prediction is significantly improved. The neural network is trained on a set of 100 proteins and then tested on a disjoint set of 1033 proteins of known structure. An average predictive accuracy of 21.7% is obtained taking the best L/2 predictions for each protein, where L is the sequence length. Taking the best L/10 predictions gives an average accuracy of 30.7%. The predictor is also tested on a set of 59 proteins from the CASP5 experiment. The accuracy is found to be relatively consistent across different sequence lengths, but to vary widely according to the secondary structure. Predictive accuracy is also found to improve by using multiple sequence alignments containing many sequences to calculate the correlations. (C) 2004 Wiley-Liss, Inc.

Computational methods for the analysis of protein structure and function

The vast majority of known proteins have not yet been experimentally characterized and little is known about their function. The design and implementation of computational tools can provide insight into the function of proteins based on their sequence, their structure, their evolutionary history and their association with other proteins. Knowledge of the three-dimensional (3D) structure of a protein can lead to a deep understanding of its mode of action and interaction, but currently the structures of <1% of sequences have been experimentally solved. For this reason, it became urgent to develop new methods that are able to computationally extract relevant information from protein sequence and structure. The starting point of my work has been the study of the properties of contacts between protein residues, since they constrain protein folding and characterize different protein structures. Prediction of residue contacts in proteins is an interesting problem whose solution may be useful in protein folding recognition and de novo design. The prediction of these contacts requires the study of the protein inter-residue distances related to the specific type of amino acid pair that are encoded in the so-called contact map. An interesting new way of analyzing those structures came out when network studies were introduced, with pivotal papers demonstrating that protein contact networks also exhibit small-world behavior. In order to highlight constraints for the prediction of protein contact maps and for applications in the field of protein structure prediction and/or reconstruction from experimentally determined contact maps, I studied to which extent the characteristic path length and clustering coefficient of the protein contacts network are values that reveal characteristic features of protein contact maps. Provided that residue contacts are known for a protein sequence, the major features of its 3D structure could be deduced by combining this knowledge with correctly predicted motifs of secondary structure. In the second part of my work I focused on a particular protein structural motif, the coiled-coil, known to mediate a variety of fundamental biological interactions. Coiled-coils are found in a variety of structural forms and in a wide range of proteins including, for example, small units such as leucine zippers that drive the dimerization of many transcription factors or more complex structures such as the family of viral proteins responsible for virus-host membrane fusion. The coiled-coil structural motif is estimated to account for 5-10% of the protein sequences in the various genomes. Given their biological importance, in my work I introduced a Hidden Markov Model (HMM) that exploits the evolutionary information derived from multiple sequence alignments, to predict coiled-coil regions and to discriminate coiled-coil sequences. The results indicate that the new HMM outperforms all the existing programs and can be adopted for the coiled-coil prediction and for large-scale genome annotation. Genome annotation is a key issue in modern computational biology, being the starting point towards the understanding of the complex processes involved in biological networks. The rapid growth in the number of protein sequences and structures available poses new fundamental problems that still deserve an interpretation. Nevertheless, these data are at the basis of the design of new strategies for tackling problems such as the prediction of protein structure and function. Experimental determination of the functions of all these proteins would be a hugely time-consuming and costly task and, in most instances, has not been carried out. As an example, currently, approximately only 20% of annotated proteins in the Homo sapiens genome have been experimentally characterized. A commonly adopted procedure for annotating protein sequences relies on the "inheritance through homology" based on the notion that similar sequences share similar functions and structures. This procedure consists in the assignment of sequences to a specific group of functionally related sequences which had been grouped through clustering techniques. The clustering procedure is based on suitable similarity rules, since predicting protein structure and function from sequence largely depends on the value of sequence identity. However, additional levels of complexity are due to multi-domain proteins, to proteins that share common domains but that do not necessarily share the same function, to the finding that different combinations of shared domains can lead to different biological roles. In the last part of this study I developed and validate a system that contributes to sequence annotation by taking advantage of a validated transfer through inheritance procedure of the molecular functions and of the structural templates. After a cross-genome comparison with the BLAST program, clusters were built on the basis of two stringent constraints on sequence identity and coverage of the alignment. The adopted measure explicity answers to the problem of multi-domain proteins annotation and allows a fine grain division of the whole set of proteomes used, that ensures cluster homogeneity in terms of sequence length. A high level of coverage of structure templates on the length of protein sequences within clusters ensures that multi-domain proteins when present can be templates for sequences of similar length. This annotation procedure includes the possibility of reliably transferring statistically validated functions and structures to sequences considering information available in the present data bases of molecular functions and structures.

Replication of bovine papillomavirus type 1 (BPV-1) DNA in Saccharomyces cerevisiae following infection with BPV-1 virions

Saccharomyces cerevisiae protoplasts exposed to bovine papillomavirus type 1 (BPV-1) virions demonstrated uptake of virions on electron microscopy. S. cerevisiae cells looked larger after exposure to BPV-1 virions, and cell wall regeneration was delayed. Southern blot hybridization of Hirt DNA from cells exposed to BPV-1 virions demonstrated BPV-1 DNA, which could be detected over 80 days of culture and at least 13 rounds of division. Two-dimensional gel analysis of Hirt DNA showed replicative intermediates, confirming that the BPV-1 genome was replicating within S. cerevisiae. Nicked circle, linear, and supercoiled BPV-1 DNA species were observed in Hirt DNA preparations from S. cerevisiae cells infected for over 50 days, and restriction digestion showed fragments hybridizing to BPV-1 in accord with the predicted restriction map for circular BPV-1 episomes. These data suggest that BPV-1 can infect S. cerevisiae and that BPV-1 episomes can replicate in the infected S. cerevisiae cells.

Dynamic arterial elastance to predict arterial pressure response to volume loading in preload-dependent patients

INTRODUCTION Hemodynamic resuscitation should be aimed at achieving not only adequate cardiac output but also sufficient mean arterial pressure (MAP) to guarantee adequate tissue perfusion pressure. Since the arterial pressure response to volume expansion (VE) depends on arterial tone, knowing whether a patient is preload-dependent provides only a partial solution to the problem. The objective of this study was to assess the ability of a functional evaluation of arterial tone by dynamic arterial elastance (Ea(dyn)), defined as the pulse pressure variation (PPV) to stroke volume variation (SVV) ratio, to predict the hemodynamic response in MAP to fluid administration in hypotensive, preload-dependent patients with acute circulatory failure. METHODS We performed a prospective clinical study in an adult medical/surgical intensive care unit in a tertiary care teaching hospital, including 25 patients with controlled mechanical ventilation who were monitored with the Vigileo(®) monitor, for whom the decision to give fluids was made because of the presence of acute circulatory failure, including arterial hypotension (MAP ≤65 mmHg or systolic arterial pressure <90 mmHg) and preserved preload responsiveness condition, defined as a SVV value ≥10%. RESULTS Before fluid infusion, Ea(dyn) was significantly different between MAP responders (MAP increase ≥15% after VE) and MAP nonresponders. VE-induced increases in MAP were strongly correlated with baseline Ea(dyn) (r(2) = 0.83; P < 0.0001). The only predictor of MAP increase was Ea(dyn) (area under the curve, 0.986 ± 0.02; 95% confidence interval (CI), 0.84-1). A baseline Ea(dyn) value >0.89 predicted a MAP increase after fluid administration with a sensitivity of 93.75% (95% CI, 69.8%-99.8%) and a specificity of 100% (95% CI, 66.4%-100%). CONCLUSIONS Functional assessment of arterial tone by Ea(dyn), measured as the PVV to SVV ratio, predicted arterial pressure response after volume loading in hypotensive, preload-dependent patients under controlled mechanical ventilation.

Finiteness theorem for topological contact equivalence of map germs

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Exploring through cover - the integrated interpretation of high resolution aeromagnetic, airborne electromagnetic and ground gravity data from the Grant's Patch area, Eastern Goldfields Province, Archaean Yigarn Craton Part B; Gravity inversion as a

The cost and risk associated with mineral exploration in Australia increases significantly as companies move into deeper regolith-covered terrain. The ability to map the bedrock and the depth of weathering within an area has the potential to decrease this risk and increase the effectiveness of exploration programs. This paper is the second in a trilogy concerning the Grant's Patch area of the Eastern Goldfields. The recent development of the VPmg potential field inversion program in conjunction with the acquisition of high-resolution gravity data over an area with extensive drilling provided an opportunity to evaluate three-dimensional gravity inversion as a bedrock and regolith mapping tool. An apparent density model of the study area was constructed, with the ground represented as adjoining 200 m by 200 m vertical rectangular prisms. During inversion VPmg incrementally adjusted the density of each prism until the free-air gravity response of the model replicated the observed data. For the Grant's Patch study area, this image of the apparent density values proved easier to interpret than the Bouguer gravity image. A regolith layer was introduced into the model and realistic fresh-rock densities assigned to each basement prism according to its interpreted lithology. With the basement and regolith densities fixed, the VPmg inversion algorithm adjusted the depth to fresh basement until the misfit between the calculated and observed gravity response was minimised. The resulting geometry of the bedrock/regolith contact largely replicated the base of weathering indicated by drilling with predicted depth of weathering values from gravity inversion typically within 15% of those logged during RAB and RC drilling.

Film thickness in a ball-on-disc contact lubricated with greases, bleed oils and base oils

Three different lubricating greases and their bleed and base oils were compared in terms of film thickness in a ball-on-disc test rig through optical interferometry. The theoretical values calculated according to Hamrock's equation are in close agreement with the base oil film thickness measurements, which validates the selected experimental methodology. The grease and bleed oil film thickness under fully flooded lubrication conditions presented quite similar behaviour and levels. Therefore, the grease film thickness under full film conditions might be predicted using their bleed oil properties, namely the viscosity and pressure-viscosity coefficient. The base and bleed oil lubricant parameter LP are proportional to the measured film thickness. A relationship between grease and the corresponding bleed oil film thickness was evidenced.

The cover basement contact beneath the Rawil axial depression (Western Alps) - True amplitude seismic processing, petrophysical properties, and modeling

Since 1986, several near-vertical seismic reflection profiles have been recorded in Switzerland in order to map the deep geologic structure of the Alps. One objective of this endeavour has been to determine the geometries of the autochthonous basement and of the external crystalline massifs, important elements for understanding the geodynamics of the Alpine orogeny. The PNR-20 seismic line W1, located in the Rawil depression of the western Swiss Alps, provides important information on this subject. It extends northward from the `'Penninic front'' across the Helvetic nappes to the Prealps. The crystalline massifs do not outcrop along this profile. Thus, the interpretation of `'near-basement'' reflections has to be constrained by down-dip projections of surface geology, `'true amplitude'' processing, rock physical property studies and modelling. 3-D seismic modelling has been used to evaluate the seismic response of two alternative down-dip projection models. To constrain the interpretation in the southern part of the profile, `'true amplitude'' processing has provided information on the strength of the reflections. Density and velocity measurements on core samples collected up-dip from the region of the seismic line have been used to evaluate reflection coefficients of typical lithologic boundaries in the region. The cover-basement contact itself is not a source of strong reflections, but strong reflections arise from within the overlaying metasedimentary cover sequence, allowing the geometry of the top of the basement to be determined on the basis of `'near-basement'' reflections. The front of the external crystalline massifs is shown to extend beneath the Prealps, about 6 km north of the expected position. A 2-D model whose seismic response shows reflection patterns very similar to the observed is proposed.

Transcription factor map alignment of promoter regions

We address the problem of comparing and characterizing the promoter regions of genes with similar expression patterns. This remains a challenging problem in sequence analysis, because often the promoter regions of co-expressed genes do not show discernible sequence conservation. In our approach, thus, we have not directly compared the nucleotide sequence of promoters. Instead, we have obtained predictions of transcription factor binding sites, annotated the predicted sites with the labels of the corresponding binding factors, and aligned the resulting sequences of labels—to which we refer here as transcription factor maps (TF-maps). To obtain the global pairwise alignment of two TF-maps, we have adapted an algorithm initially developed to align restriction enzyme maps. We have optimized the parameters of the algorithm in a small, but well-curated, collection of human–mouse orthologous gene pairs. Results in this dataset, as well as in an independent much larger dataset from the CISRED database, indicate that TF-map alignments are able to uncover conserved regulatory elements, which cannot be detected by the typical sequence alignments.

Effect of tibial tray inclination on the femoro-tibial contact pattern of a mobile total knee arthroplasty

Introduction: The posterior inclination of the tibial component is an important factor that can affect the success of total knee arthroplasty. It can reduce the posterior impingement and thus increase the range of flexion, but it may also induce instability in flexion, anterior impingement between the polyethylene of postero-stabilizing knee prosthesis, and anterior conflict with the cortical bone and the stem. Although the problem is identified, there is still a debate on the ideal inclination angle and the surgical technique to avoid an excessive posterior inclination. The aim of this study was to predict the effect of a posterior inclination of the tibial component on the contact pattern on the tibial insert, using a numerical musculoskeletal model of the knee joint. Methods: A 3D finite element model of the knee joint was developed to simulate an active and loaded squat movement after total knee arthroplasty. Flexion was actively controlled by the quadriceps muscle and muscle activations were estimated from EMG data and were synchronized by a feedback algorithm. Two inclinations of the tibial tray were considered: a posterior inclination of 0° or 10°. During the entire range of flexion, the following quantities were calculated: the tibiofemoral and patello-femoral contact force, and the contact pattern on polyethylene insert. The antero-posterior displacement of the contact pattern was also measured. Abaqus 6.7 was used for all analyses. Results: The tibio-femoral and patello-femoral contact forces increased during flexion and reached respectively 4 and 7 BW (bodyweight) at 90° of flexion. They were slightly affected by the inclination of the tibial tray. Without posterior inclination, the contact pattern on the tibial insert remained centered. The contact pressure was lower than 5 MPa below 60° of flexion, but exceeded 20 MPa at 90° of flexion. The posterior inclination displaced the contact point posteriorly by 2 to 4 mm. Conclusion: The inclination of the tibial tray displaced the contactpattern towards the posterior border of the tibial insert. However, even for 10° of inclination, the contact center remained far from the posterior border (12 mm). There was no instability predicted for this movement.

CFD Modelling of Direct Contact Condensation in Suppression Pools by Applying Condensation Models of Separated Flow

The condensation rate has to be high in the safety pressure suppression pool systems of Boiling Water Reactors (BWR) in order to fulfill their safety function. The phenomena due to such a high direct contact condensation (DCC) rate turn out to be very challenging to be analysed either with experiments or numerical simulations. In this thesis, the suppression pool experiments carried out in the POOLEX facility of Lappeenranta University of Technology were simulated. Two different condensation modes were modelled by using the 2-phase CFD codes NEPTUNE CFD and TransAT. The DCC models applied were the typical ones to be used for separated flows in channels, and their applicability to the rapidly condensing flow in the condensation pool context had not been tested earlier. A low Reynolds number case was the first to be simulated. The POOLEX experiment STB-31 was operated near the conditions between the ’quasi-steady oscillatory interface condensation’ mode and the ’condensation within the blowdown pipe’ mode. The condensation models of Lakehal et al. and Coste & Lavi´eville predicted the condensation rate quite accurately, while the other tested ones overestimated it. It was possible to get the direct phase change solution to settle near to the measured values, but a very high resolution of calculation grid was needed. Secondly, a high Reynolds number case corresponding to the ’chugging’ mode was simulated. The POOLEX experiment STB-28 was chosen, because various standard and highspeed video samples of bubbles were recorded during it. In order to extract numerical information from the video material, a pattern recognition procedure was programmed. The bubble size distributions and the frequencies of chugging were calculated with this procedure. With the statistical data of the bubble sizes and temporal data of the bubble/jet appearance, it was possible to compare the condensation rates between the experiment and the CFD simulations. In the chugging simulations, a spherically curvilinear calculation grid at the blowdown pipe exit improved the convergence and decreased the required cell count. The compressible flow solver with complete steam-tables was beneficial for the numerical success of the simulations. The Hughes-Duffey model and, to some extent, the Coste & Lavi´eville model produced realistic chugging behavior. The initial level of the steam/water interface was an important factor to determine the initiation of the chugging. If the interface was initialized with a water level high enough inside the blowdown pipe, the vigorous penetration of a water plug into the pool created a turbulent wake which invoked the chugging that was self-sustaining. A 3D simulation with a suitable DCC model produced qualitatively very realistic shapes of the chugging bubbles and jets. The comparative FFT analysis of the bubble size data and the pool bottom pressure data gave useful information to distinguish the eigenmodes of chugging, bubbling, and pool structure oscillations.

Mathematical model to estimate of the deterioration of wooden poles in contact with soil used in rural areas

In São Paulo State, mainly in rural areas, the utilization of wooden poles is observed for different purposes. In this context, wood in contact with the ground presents faster deterioration, which is generally associated to environmental factors and, especially to the presence of fungi and insects. With the use of mathematical models, the useful life of wooden structures can be predicted by obtaining "climatic indexes" to indicate, comparatively among the areas studied, which have more or less tendency to fungi and insects attacks. In this work, by using climatological data of several cities at São Paulo State, a simplified mathematical model was obtained to measure the aggressiveness of the wood in contact with the soil.

Trans-Ancestral Studies Fine Map the SLE-Susceptibility Locus TNFSF4

We previously established an 80 kb haplotype upstream of TNFSF4 as a susceptibility locus in the autoimmune disease SLE. SLE-associated alleles at this locus are associated with inflammatory disorders, including atherosclerosis and ischaemic stroke. In Europeans, the TNFSF4 causal variants have remained elusive due to strong linkage disequilibrium exhibited by alleles spanning the region. Using a trans-ancestral approach to fine-map the locus, utilising 17,900 SLE and control subjects including Amerindian/Hispanics (1348 cases, 717 controls), African-Americans (AA) (1529, 2048) and better powered cohorts of Europeans and East Asians, we find strong association of risk alleles in all ethnicities; the AA association replicates in African-American Gullah (152,122). The best evidence of association comes from two adjacent markers: rs2205960-T (P = 1.71×10-34, OR = 1.43[1.26-1.60]) and rs1234317-T (P = 1.16×10-28, OR = 1.38[1.24-1.54]). Inference of fine-scale recombination rates for all populations tested finds the 80 kb risk and non-risk haplotypes in all except African-Americans. In this population the decay of recombination equates to an 11 kb risk haplotype, anchored in the 5′ region proximal to TNFSF4 and tagged by rs2205960-T after 1000 Genomes phase 1 (v3) imputation. Conditional regression analyses delineate the 5′ risk signal to rs2205960-T and the independent non-risk signal to rs1234314-C. Our case-only and SLE-control cohorts demonstrate robust association of rs2205960-T with autoantibody production. The rs2205960-T is predicted to form part of a decameric motif which binds NF-κBp65 with increased affinity compared to rs2205960-G. ChIP-seq data also indicate NF-κB interaction with the DNA sequence at this position in LCL cells. Our research suggests association of rs2205960-T with SLE across multiple groups and an independent non-risk signal at rs1234314-C. rs2205960-T is associated with autoantibody production and lymphopenia. Our data confirm a global signal at TNFSF4 and a role for the expressed product at multiple stages of lymphocyte dysregulation during SLE pathogenesis. We confirm the validity of trans-ancestral mapping in a complex trait. © 2013 Manku et al.

Predicted soil organic carbon stocks and changes in Jordan between 2000 and 2030 made using the GEFSOC modelling system

Estimates of soil organic carbon (SOC) stocks and changes under different land use systems can help determine vulnerability to land degradation. Such information is important for countries in and areas with high susceptibility to desertification. SOC stocks, and predicted changes between 2000 and 2030, were determined at the national scale for Jordan using The Global Environment Facility Soil Organic Carbon (GEFSOC) Modelling System. For the purpose of this study, Jordan was divided into three natural regions (The Jordan Valley, the Uplands and the Badia) and three developmental regions (North, Middle and South). Based on this division, Jordan was divided into five zones (based on the dominant land use): the Jordan Valley, the North Uplands, the Middle Uplands, the South Uplands and the Badia. This information was merged using GIS, along with a map of rainfall isohyets, to produce a map with 498 polygons. Each of these was given a unique ID, a land management unit identifier and was characterized in terms of its dominant soil type. Historical land use data, current land use and future land use change scenarios were also assembled, forming major inputs of the modelling system. The GEFSOC Modelling System was then run to produce C stocks in Jordan for the years 1990, 2000 and 2030. The results were compared with conventional methods of estimating carbon stocks, such as the mapping based SOTER method. The results of these comparisons showed that the model runs are acceptable, taking into consideration the limited availability of long-term experimental soil data that can be used to validate them. The main findings of this research show that between 2000 and 2030, SOC may increase in heavily used areas under irrigation and will likely decrease in grazed rangelands that cover most of Jordan giving an overall decrease in total SOC over time if the land is indeed used under the estimated forms of land use. (C) 2007 Elsevier B.V. All rights reserved.

Predicted soil organic carbon stocks and changes in the Brazilian Amazon between 2000 and 2030

Currently we have little understanding of the impacts of land use change on soil C stocks in the Brazilian Amazon. Such information is needed to determine impacts'6n the global C cycle and the sustainability of agricultural systems that are replacing native forest. The aim of this study was to predict soil carbon stocks and changes in the Brazilian Amazon during the period between 2000 and 2030, using the GEFSOC soil carbon (C) modelling system. In order to do so, we devised current and future land use scenarios for the Brazilian Amazon, taking into account: (i) deforestation, rates from the past three decades, (ii) census data on land use from 1940 to 2000, including the expansion and intensification of agriculture in the region, (iii) available information on management practices, primarily related to well managed pasture versus degraded pasture and conventional systems versus no-tillage systems for soybean (Glycine max) and (iv) FAO predictions on agricultural land use and land use changes for the years 2015 and 2030. The land use scenarios were integrated with spatially explicit soils data (SOTER database), climate, potential natural vegetation and land management units using the recently developed GEFSOC soil C modelling system. Results are presented in map, table and graph form for the entire Brazilian Amazon for the current situation (1990 and 2000) and the future (2015 and 2030). Results include soil organic C (SOC) stocks and SOC stock change rates estimated by three methods: (i) the Century ecosystem model, (ii) the Rothamsted C model and (iii) the intergovernmental panel on climate change (IPCC) method for assessing soil C at regional scale. In addition, we show estimated values of above and belowground biomass for native vegetation, pasture and soybean. The results on regional SOC stocks compare reasonably well with those based on mapping approaches. The GEFSOC system provided a means of efficiently handling complex interactions among biotic-edapho-climatic conditions (> 363,000 combinations) in a very large area (similar to 500 Mha) such as the Brazilian Amazon. All of the methods used showed a decline in SOC stock for the period studied; Century and RothC simulated values for 2030 being about 7% lower than those in 1990. Values from Century and RothC (30,430 and 25,000 Tg for the 0-20 cm layer for the Brazilian Amazon region were higher than those obtained from the IPCC system (23,400 Tg in the 0-30 cm layer). Finally; our results can help understand the major biogeochemical cycles that influence soil fertility and help devise management strategies that enhance the sustainability of these areas and thus slow further deforestation. (C) 2007 Elsevier B.V. All rights reserved.