Analysis on conservation of disulphide bonds and their structural features in homologous protein domain families

Background: Disulphide bridges are well known to play key roles in stability, folding and functions of proteins. Introduction or deletion of disulphides by site-directed mutagenesis have produced varying effects on stability and folding depending upon the protein and location of disulphide in the 3-D structure. Given the lack of complete understanding it is worthwhile to learn from an analysis of extent of conservation of disulphides in homologous proteins. We have also addressed the question of what structural interactions replaces a disulphide in a homologue in another homologue. Results: Using a dataset involving 34,752 pairwise comparisons of homologous protein domains corresponding to 300 protein domain families of known 3-D structures, we provide a comprehensive analysis of extent of conservation of disulphide bridges and their structural features. We report that only 54% of all the disulphide bonds compared between the homologous pairs are conserved, even if, a small fraction of the non-conserved disulphides do include cytoplasmic proteins. Also, only about one fourth of the distinct disulphides are conserved in all the members in protein families. We note that while conservation of disulphide is common in many families, disulphide bond mutations are quite prevalent. Interestingly, we note that there is no clear relationship between sequence identity between two homologous proteins and disulphide bond conservation. Our analysis on structural features at the sites where cysteines forming disulphide in one homologue are replaced by non-Cys residues show that the elimination of a disulphide in a homologue need not always result in stabilizing interactions between equivalent residues. Conclusion: We observe that in the homologous proteins, disulphide bonds are conserved only to a modest extent. Very interestingly, we note that extent of conservation of disulphide in homologous proteins is unrelated to the overall sequence identity between homologues. The non-conserved disulphides are often associated with variable structural features that were recruited to be associated with differentiation or specialisation of protein function.

Regolith mass balance inferred from combined mineralogical, geochemical and geophysical studies: Mule Hole gneissic watershed, South India

The aim of this study is to propose a method to assess the long-term chemical weathering mass balance for a regolith developed on a heterogeneous silicate substratum at the small experimental watershed scale by adopting a combined approach of geophysics, geochemistry and mineralogy. We initiated in 2003 a study of the steep climatic gradient and associated geomorphologic features of the edge of the rifted continental passive margin of the Karnataka Plateau, Peninsular India. In the transition sub-humid zone of this climatic gradient we have studied the pristine forested small watershed of Mule Hole (4.3 km(2)) mainly developed on gneissic substratum. Mineralogical, geochemical and geophysical investigations were carried out (i) in characteristic red soil profiles and (ii) in boreholes up to 60 m deep in order to take into account the effect of the weathering mantle roots. In addition, 12 Electrical Resistivity Tomography profiles (ERT), with an investigation depth of 30 m, were generated at the watershed scale to spatially characterize the information gathered in boreholes and soil profiles. The location of the ERT profiles is based on a previous electromagnetic survey, with an investigation depth of about 6 m. The soil cover thickness was inferred from the electromagnetic survey combined with a geological/pedological survey. Taking into account the parent rock heterogeneity, the degree of weathering of each of the regolith samples has been defined using both the mineralogical composition and the geochemical indices (Loss on Ignition, Weathering Index of Parker, Chemical Index of Alteration). Comparing these indices with electrical resistivity logs, it has been found that a value of 400 Ohm m delineates clearly the parent rocks and the weathered materials, Then the 12 inverted ERT profiles were constrained with this value after verifying the uncertainty due to the inversion procedure. Synthetic models based on the field data were used for this purpose. The estimated average regolith thickness at the watershed scale is 17.2 m, including 15.2 m of saprolite and 2 m of soil cover. Finally, using these estimations of the thicknesses, the long-term mass balance is calculated for the average gneiss-derived saprolite and red soil. In the saprolite, the open-system mass-transport function T indicates that all the major elements except Ca are depleted. The chlorite and biotite crystals, the chief sources for Mg (95%), Fe (84%), Mn (86%) and K (57%, biotite only), are the first to undergo weathering and the oligoclase crystals are relatively intact within the saprolite with a loss of only 18%. The Ca accumulation can be attributed to the precipitation of CaCO3 from the percolating solution due to the current and/or the paleoclimatic conditions. Overall, the most important losses occur for Si, Mg and Na with -286 x 10(6) mol/ha (62% of the total mass loss), -67 x 10(6) mol/ha (15% of the total mass loss) and -39 x 10(6) mol/ha (9% of the total mass loss), respectively. Al, Fe and K account for 7%, 4% and 3% of the total mass loss, respectively. In the red soil profiles, the open-system mass-transport functions point out that all major elements except Mn are depleted. Most of the oligoclase crystals have broken down with a loss of 90%. The most important losses occur for Si, Na and Mg with -55 x 10(6) mol/ha (47% of the total mass loss), -22 x 10(6) mol/ha (19% of the total mass loss) and -16 x 10(6) mol/ha (14% of the total mass loss), respectively. Ca, Al, K and Fe account for 8%, 6%, 4% and 2% of the total mass loss, respectively. Overall these findings confirm the immaturity of the saprolite at the watershed scale. The soil profiles are more evolved than saprolite but still contain primary minerals that can further undergo weathering and hence consume atmospheric CO2.

Structure and magnetotransport properties of the layered manganites REl.2Srl.8Mn,07(R E = La, Pr, Nd)

We have examined the magnetotransport properties and the structure, by Rietveld refinement of powder X-ray data, of the phases RE(1.2)Sr(1.8)Mn(2)O(7) (RE = La, Pr, Nd). We find that on cooling, La1.2Sr1.8Mn2O7 undergoes a transition to a nearly perfect ferromagnet with 90% magnetization at 1.45 T, as reported by earlier workers, but the Pr and Nd phases show only a small magnetization that grows gradually as the temperature is decreased. There seems to be significant correlation between electrical transport and the Jahn-Teller elongation of the apical Mn-O bonds in these systems. The elongation of the apical Mn-O bonds forces the nine-coordinate rock-salt site to be occupied preferentially by the smaller rare-earth-metal cations. This preferential occupation is reliably obtained from the X-ray refinement. All three title phases show a magnetoresistance ratio of about 4(corresponding to a magnetoresistance, [R(0)-R(H)]/R(0), of about 75%) at a field of 7 T and temperatures around 100 K.

Contribution of geophysical surveys to groundwater modelling of a porous aquifer in semiarid Niger: An overview

Subsurface geophysical surveys were carried out using a large range of methods in an unconfined sandstone aquifer in semiarid south-western Niger for improving both the conceptual model of water flow through the unsaturated zone and the parameterization of numerical a groundwater model of the aquifer. Methods included: electromagnetic mapping, electrical resistivity tomography (ERT), resistivity logging, time domain electromagnetic sounding (TDEM), and magnetic resonance sounding (MRS). Analyses of electrical conductivities, complemented by geochemical measurements, allowed us to identify preferential pathways for infiltration and drainage beneath gullies and alluvial fans. The mean water content estimated by MRS (13%) was used for computing the regional groundwater recharge from long-term change in the water table. The ranges in permeability and water content obtained with MRS allowed a reduction of the degree of freedom of aquifer parameters used in groundwater modelling.

Formation and preservation of pedogenic carbonates in South India, links with paleo-monsoon and pedological conditions: Clues from Sr isotopes, U-Th series and REEs

The influence of the pedogenic and climatic contexts on the formation and preservation of pedogenic carbonates in a climosequence in the Western Ghats (Karnataka Plateau, South West India) has been studied. Along the climosequence, the current mean annual rainfall (MAR) varies within a 80 km transect from 6000 mm at the edge of the Plateau to 500 mm inland. Pedogenic carbonates occur in the MAR range of 500-1200 mm. In the semi-arid zone (MAR: 500-900 mm), carbonates occur (i) as rhick hardpan calcretes on pediment slopes and (ii) as nodular horizons in polygenic black soils (i.e. vertisols). In the sub-humid zone (MAR: 900-1500 mm), pedogenic carbonates are disseminated in the black soil matrices either as loose, irregular and friable nodules of millimetric size or as indurated botryoidal nodules of centimetric to pluricentimetric size. They also occur at the top layers of the saprolite either as disseminated pluricentimetric indurated nodules or carbonate-cemented lumps of centimetric to decimetric size. Chemical and isotopic (Sr-87/Sr-86) compositions of the carbonate fraction were determined after leaching with 0.25 N HCl. The corresponding residual fractions containing both primary minerals and authigenic clays were digested separately and analyzed. The trend defined by the Sr-87/Sr-86 signatures of both labile carbonate fractions and corresponding residual fractions indicates that a part of the labile carbonate fraction is genetically linked to the local soil composition. Considering the residual fraction of each sample as the most likely lithogenic source of Ca in carbonates, it is estimated that from 24% to 82% (55% on average) of Ca is derived from local bedrock weathering, leading to a consumption of an equivalent proportion of atmospheric CO2. These values indicate that climatic conditions were humid enough to allow silicate weathering: MAR at the time of carbonate formation likely ranged from 400 to 700 mm, which is 2- to 3-fold less than the current MAR at these locations. The Sr, U and Mg contents and the (U-234/U-238) activity ratio in the labile carbonate fraction help to understand the conditions of carbonate formation. The relatively high concentrations of Sr, U and Mg in black soil carbonates may indicate fast growth and accumulation compared to carbonates in saprolite, possibly due to a better confinement of the pore waters which is supported by their high (U-234/U-238) signatures, and/or to higher content of dissolved carbonates in the pore waters. The occurrence of Ce, Mn and Fe oxides in the cracks of carbonate reflects the existence of relatively humid periods after carbonate formation. The carbonate ages determined by the U-Th method range from 1.33 +/- 0.84 kyr to 7.5 +/- 2.7 kyr and to a cluster of five ages around 20 kyr, i.e. the Last Glacial Maximum period. The young occurrences are only located in the black soils, which therefore constitute sensitive environments for trapping and retaining atmospheric CO2 even on short time scales. The maximum age of carbonates depends on their location in the climatic gradient: from about 20 kyr for centimetric nodules at Mule Hole (MAR = 1100 mm/yr) to 200 kyr for the calcrete at Gundlupet (MAR = 700 mm/yr, Durand et al., 2007). The intensity of rainfall during wet periods would indeed control the lifetime of pedogenic carbonates and thus the duration of inorganic carbon storage in soils. (C) 2010 Elsevier Ltd. All rights reserved.

DoSA: Database of Structural Alignments

Protein structure alignment is a crucial step in protein structure-function analysis. Despite the advances in protein structure alignment algorithms, some of the local conformationally similar regions are mislabeled as structurally variable regions (SVRs). These regions are not well superimposed because of differences in their spatial orientations. The Database of Structural Alignments (DoSA) addresses this gap in identification of local structural similarities obscured in global protein structural alignments by realigning SVRs using an algorithm based on protein blocks. A set of protein blocks is a structural alphabet that abstracts protein structures into 16 unique local structural motifs. DoSA provides unique information about 159 780 conformationally similar and 56 140 conformationally dissimilar SVRs in 74 705 pairwise structural alignments of homologous proteins. The information provided on conformationally similar and dissimilar SVRs can be helpful to model loop regions. It is also conceivable that conformationally similar SVRs with conserved residues could potentially contribute toward functional integrity of homologues, and hence identifying such SVRs could be helpful in understanding the structural basis of protein function.

The Escherichia coli Phosphotyrosine Proteome Relates to Core Pathways and Virulence

While phosphotyrosine modification is an established regulatory mechanism in eukaryotes, it is less well characterized in bacteria due to low prevalence. To gain insight into the extent and biological importance of tyrosine phosphorylation in Escherichia coli, we used immunoaffinity-based phosphotyrosine peptide enrichment combined with high resolution mass spectrometry analysis to comprehensively identify tyrosine phosphorylated proteins and accurately map phosphotyrosine sites. We identified a total of 512 unique phosphotyrosine sites on 342 proteins in E. coli K12 and the human pathogen enterohemorrhagic E. coli (EHEC) O157:H7, representing the largest phosphotyrosine proteome reported to date in bacteria. This large number of tyrosine phosphorylation sites allowed us to define five phosphotyrosine site motifs. Tyrosine phosphorylated proteins belong to various functional classes such as metabolism, gene expression and virulence. We demonstrate for the first time that proteins of a type III secretion system (T3SS), required for the attaching and effacing (A/E) lesion phenotype characteristic for intestinal colonization by certain EHEC strains, are tyrosine phosphorylated by bacterial kinases. Yet, A/E lesion and metabolic phenotypes were unaffected by the mutation of the two currently known tyrosine kinases, Etk and Wzc. Substantial residual tyrosine phosphorylation present in an etk wzc double mutant strongly indicated the presence of hitherto unknown tyrosine kinases in E. coli. We assess the functional importance of tyrosine phosphorylation and demonstrate that the phosphorylated tyrosine residue of the regulator SspA positively affects expression and secretion of T3SS proteins and formation of A/E lesions. Altogether, our study reveals that tyrosine phosphorylation in bacteria is more prevalent than previously recognized, and suggests the involvement of phosphotyrosine-mediated signaling in a broad range of cellular functions and virulence.

Correlation between local structural dynamics of proteins inferred from NMR ensembles and evolutionary dynamics of homologues of known structure

Conformational changes in proteins are extremely important for their biochemical functions. Correlation between inherent conformational variations in a protein and conformational differences in its homologues of known structure is still unclear. In this study, we have used a structural alphabet called Protein Blocks (PBs). PBs are used to perform abstraction of protein 3-D structures into a 1-D strings of 16 alphabets (a-p) based on dihedral angles of overlapping pentapeptides. We have analyzed the variations in local conformations in terms of PBs represented in the ensembles of 801 protein structures determined using NMR spectroscopy. In the analysis of concatenated data over all the residues in all the NMR ensembles, we observe that the overall nature of inherent local structural variations in NMR ensembles is similar to the nature of local structural differences in homologous proteins with a high correlation coefficient of .94. High correlation at the alignment positions corresponding to helical and beta-sheet regions is only expected. However, the correlation coefficient by considering only the loop regions is also quite high (.91). Surprisingly, segregated position-wise analysis shows that this high correlation does not hold true to loop regions at the structurally equivalent positions in NMR ensembles and their homologues of known structure. This suggests that the general nature of local structural changes is unique; however most of the local structural variations in loop regions of NMR ensembles do not correlate to their local structural differences at structurally equivalent positions in homologues.

THE KNICKKOPF DOMON DOMAIN IS ESSENTIAL FOR CUTICLE DIFFERENTIATION IN Drosophila melanogaster

The dopamine monoxygenase N-terminal (DOMON) domain is found in extracellular proteins across several eukaryotic and prokaryotic taxa. It has been proposed that this domain binds to heme or sugar moieties. Here, we have analyzed the role of four highly conserved amino acids in the DOMON domain of the Drosophila melanogaster Knickkopf protein that is inserted into the apical plasma membrane and assists extracellular chitin organization. In principal, we generated Knickkopf versions with exchanged residues tryptophan(299,) methionine(333), arginine(401), or histidine(437), and scored for the ability of the respective engineered protein to normalize the knickkopf mutant phenotype. Our results confirm the absolute necessity of tryptophan(299,) methionine(333), and histidine(437) for Knickkopf function and stability, the latter two being predicted to be critical for heme binding. In contrast, arginine(401) is required for full efficiency of Knickkopf activity. Taken together, our genetic data support the prediction of these residues to mediate the function of Knickkopf during cuticle differentiation in insects. Hence, the DOMON domain is apparently an essential factor contributing to the construction of polysaccharide-based extracellular matrices.

Use of a structural alphabet to find compatible folds for amino acid sequences

The structural annotation of proteins with no detectable homologs of known 3D structure identified using sequence-search methods is a major challenge today. We propose an original method that computes the conditional probabilities for the amino-acid sequence of a protein to fit to known protein 3D structures using a structural alphabet, known as Protein Blocks (PBs). PBs constitute a library of 16 local structural prototypes that approximate every part of protein backbone structures. It is used to encode 3D protein structures into 1D PB sequences and to capture sequence to structure relationships. Our method relies on amino acid occurrence matrices, one for each PB, to score global and local threading of query amino acid sequences to protein folds encoded into PB sequences. It does not use any information from residue contacts or sequence-search methods or explicit incorporation of hydrophobic effect. The performance of the method was assessed with independent test datasets derived from SCOP 1.75A. With a Z-score cutoff that achieved 95% specificity (i.e., less than 5% false positives), global and local threading showed sensitivity of 64.1% and 34.2%, respectively. We further tested its performance on 57 difficult CASP10 targets that had no known homologs in PDB: 38 compatible templates were identified by our approach and 66% of these hits yielded correctly predicted structures. This method scales-up well and offers promising perspectives for structural annotations at genomic level. It has been implemented in the form of a web-server that is freely available at http://www.bo-protscience.fr/forsa.

Key challenges for the creation and maintenance of specialist protein resources

As the volume of data relating to proteins increases, researchers rely more and more on the analysis of published data, thus increasing the importance of good access to these data that vary from the supplemental material of individual articles, all the way to major reference databases with professional staff and long-term funding. Specialist protein resources fill an important middle ground, providing interactive web interfaces to their databases for a focused topic or family of proteins, using specialized approaches that are not feasible in the major reference databases. Many are labors of love, run by a single lab with little or no dedicated funding and there are many challenges to building and maintaining them. This perspective arose from a meeting of several specialist protein resources and major reference databases held at the Wellcome Trust Genome Campus (Cambridge, UK) on August 11 and 12, 2014. During this meeting some common key challenges involved in creating and maintaining such resources were discussed, along with various approaches to address them. In laying out these challenges, we aim to inform users about how these issues impact our resources and illustrate ways in which our working together could enhance their accuracy, currency, and overall value. Proteins 2015; 83:1005-1013. (c) 2015 The Authors. Proteins: Structure, Function, and Bioinformatics Published by Wiley Periodicals, Inc.

Creating a specialist protein resource network: a meeting report for the protein bioinformatics and community resources retreat

During 11-12 August 2014, a Protein Bioinformatics and Community Resources Retreat was held at the Wellcome Trust Genome Campus in Hinxton, UK. This meeting brought together the principal investigators of several specialized protein resources (such as CAZy, TCDB and MEROPS) as well as those from protein databases from the large Bioinformatics centres (including UniProt and RefSeq). The retreat was divided into five sessions: (1) key challenges, (2) the databases represented, (3) best practices for maintenance and curation, (4) information flow to and from large data centers and (5) communication and funding. An important outcome of this meeting was the creation of a Specialist Protein Resource Network that we believe will improve coordination of the activities of its member resources. We invite further protein database resources to join the network and continue the dialogue.