Structure-Based Phylogeny as a Diagnostic for Functional Characterization of Proteins with a Cupin Fold

Background: The members of cupin superfamily exhibit large variations in their sequences, functions, organization of domains, quaternary associations and the nature of bound metal ion, despite having a conserved beta-barrel structural scaffold. Here, an attempt has been made to understand structure-function relationships among the members of this diverse superfamily and identify the principles governing functional diversity. The cupin superfamily also contains proteins for which the structures are available through world-wide structural genomics initiatives but characterized as ``hypothetical''. We have explored the feasibility of obtaining clues to functions of such proteins by means of comparative analysis with cupins of known structure and function. Methodology/Principal Findings: A 3-D structure-based phylogenetic approach was undertaken. Interestingly, a dendrogram generated solely on the basis of structural dissimilarity measure at the level of domain folds was found to cluster functionally similar members. This clustering also reflects an independent evolution of the two domains in bicupins. Close examination of structural superposition of members across various functional clusters reveals structural variations in regions that not only form the active site pocket but are also involved in interaction with another domain in the same polypeptide or in the oligomer. Conclusions/Significance: Structure-based phylogeny of cupins can influence identification of functions of proteins of yet unknown function with cupin fold. This approach can be extended to other proteins with a common fold that show high evolutionary divergence. This approach is expected to have an influence on the function annotation in structural genomics initiatives.

A Search for Energy Minimized Sequences of Proteins

In this paper, we present numerical evidence that supports the notion of minimization in the sequence space of proteins for a target conformation. We use the conformations of the real proteins in the Protein Data Bank (PDB) and present computationally efficient methods to identify the sequences with minimum energy. We use edge-weighted connectivity graph for ranking the residue sites with reduced amino acid alphabet and then use continuous optimization to obtain the energy-minimizing sequences. Our methods enable the computation of a lower bound as well as a tight upper bound for the energy of a given conformation. We validate our results by using three different inter-residue energy matrices for five proteins from protein data bank (PDB), and by comparing our energy-minimizing sequences with 80 million diverse sequences that are generated based on different considerations in each case. When we submitted some of our chosen energy-minimizing sequences to Basic Local Alignment Search Tool (BLAST), we obtained some sequences from non-redundant protein sequence database that are similar to ours with an E-value of the order of 10(-7). In summary, we conclude that proteins show a trend towards minimizing energy in the sequence space but do not seem to adopt the global energy-minimizing sequence. The reason for this could be either that the existing energy matrices are not able to accurately represent the inter-residue interactions in the context of the protein environment or that Nature does not push the optimization in the sequence space, once it is able to perform the function.

Inhibition of Tyrosine Phosphorylation of Sperm Flagellar Proteins, Outer Dense Fiber Protein-2 and Tektin-2, Is Associated With Impaired Motility During Capacitation of Hamster Spermatozoa

In mammals, acquisition of fertilization competence of spermatozoa is dependent on the phenomenon of sperm capacitation. One of the critical molecular events of sperm capacitation is protein tyrosine phosphorylation. In a previous study, we demonstrated that a specific epidermal growth factor receptor (EGFR)-tyrosine kinase inhibitor, tyrphostin-A47, inhibited hamster sperm capacitation, accompanied by a reduced sperm protein tyrosine phosphorylation. Interestingly, a high percentage of tyrphostin-A47-treated spermatozoa exhibited circular motility, which was associated with a distinct hypo-tyrosine phosphorylation of flagellar proteins, predominantly of Mr 45,000-60,000. In this study, we provide evidence on the localization of capacitation-associated tyrosine-phosphorylated proteins to the nonmembranous, structural components of the sperm flagellum. Consistent with this, we show their ultrastructural localization in the outer dense fiber, axoneme, and fibrous sheath of spermatozoa. Among hypo-tyrosine phosphorylated major proteins of tyrphostin-A47-treated spermatozoa, we identified the 45 kDa protein as outer dense fiber protein-2 and the 51 kDa protein as tektin-2, components of the sperm outer dense fiber and axoneme, respectively. This study shows functional association of hypo-tyrosine-phosphorylation status of outer dense fiber protein-2 and tektin-2 with impaired flagellar bending of spermatozoa, following inhibition of EGFR-tyrosine kinase, thereby showing the critical importance of flagellar protein tyrosine phosphorylation during capacitation and hyperactivation of hamster spermatozoa.

Molecular characterisation of ABC transporter type FtsE and FtsX proteins of Mycobacterium tuberculosis.

Elicitation of drug resistance and various survival strategies inside host macrophages have been the hallmarks of Mycobacterium tuberculosis as a successful pathogen. ATP Binding Cassette (ABC) transporter type proteins are known to be involved in the efflux of drugs in bacterial and mammalian systems. FtsE, an ABC transporter type protein, in association with the integral membrane protein FtsX, is involved in the assembly of potassium ion transport proteins and probably of cell division proteins as well, both of which being relevant to tubercle bacillus. In this study, we cloned ftsE gene of M. tuberculosis, overexpressed and purified. The recombinant MtFtsE-6xHis protein and the native MtFtsE protein were found localized on the membrane of E. coli and M. tuberculosis cells, respectively. MtFtsE-6xHis protein showed ATP binding in vitro, for which the K42 residue in the Walker A motif was found essential. While MtFtsE-6xHis protein could partially complement growth defect of E. coli ftsE temperature-sensitive strain MFT1181, co-expression of MtFtsE and MtFtsX efficiently complemented the growth defect, indicating that the MtFtsE and MtFtsX proteins might be performing an associated function. MtFtsE and MtFtsX-6xHis proteins were found to exist as a complex on the membrane of E. coli cells co-expressing the two proteins.

Characterization of single-stranded DNA-binding proteins from Mycobacteria. The carboxyl-terminal of domain of SSB is essential for stable association with its cognate RecA protein

Single-stranded DNA-binding proteins (SSB) play an important role in most aspects of DNA metabolism including DNA replication, repair, and recombination. We report here the identification and characterization of SSB proteins of Mycobacterium smegmatis and Mycobacterium tuberculosis. Sequence comparison of M. smegmatis SSB revealed that it is homologous to M. tuberculosis SSB, except for a small spacer connecting the larger amino-terminal domain with the extreme carboxyl-terminal tail. The purified SSB proteins of mycobacteria bound single-stranded DNA with high affinity, and the association and dissociation constants were similar to that of the prototype SSB. The proteolytic signatures of free and bound forms of SSB proteins disclosed that DNA binding was associated with structural changes at the carboxyl-terminal domain. Significantly, SSB proteins from mycobacteria displayed high affinity for cognate RecA, whereas Escherichia coli SSB did not under comparable experimental conditions. Accordingly, SSB and RecA were coimmunoprecipitated from cell lysates, further supporting an interaction between these proteins in vivo. The carboxyl-terminal domain of M. smegmatis SSB, which is not essential for interaction with ssDNA, is the site of binding of its cognate RecA. These studies provide the first evidence for stable association of eubacterial SSB proteins with their cognate RecA, suggesting that these two proteins might function together during DNA repair and/or recombination.

High lithium storage in micrometre sized mesoporous spherical self-assembly of anatase titania nanospheres and carbon

Composite of anatase titania (TiO2) nanospheres and carbon grown and self-assembled into micron-sized mesoporous spheres via a solvothermal synthesis route are discussed here in the context of rechargeable lithium-ion battery. The morphology and carbon content and hence the electrochemical performance are observed to be significantly influenced by the synthesis parameters. Synthesis conditions resulting in a mesoporous arrangement of an optimized amount carbon and TiO2 exhibited the best lithium battery performance. The first discharge cycle capacity of carbon-titania mesoporous spheres (solvothermal reaction at 150 degrees C at 6 h, calcination at 500 degrees C under air, BET surface area 80 m(2)g(-1)) was 334 mAhg(-1) (approximately 1 Li) at current rate of 0.066 Ag-1. High storage capacity and good cyclability is attributed to the nanostructuring of TiO2 (mesoporosity) as well as due to formation of a percolation network of carbon around the TiO2 nanoparticles. The micron-sized mesoporous spheres of carbon-titania composite nanoparticles also show good rate cyclability in the range (0.066-6.67) Ag-1.

A Method to Find Sequentially Separated Motifs in Biological Sequences (SSMBS)

Sequence motifs occurring in a particular order in proteins or DNA have been proved to be of biological interest. In this paper, a new method to locate the occurrences of up to five user-defined motifs in a specified order in large proteins and in nucleotide sequence databases is proposed. It has been designed using the concept of quantifiers in regular expressions and linked lists for data storage. The application of this method includes the extraction of relevant consensus regions from biological sequences. This might be useful in clustering of protein families as well as to study the correlation between positions of motifs and their functional sites in DNA sequences.

Water balance modelling in a tropical watershed under deciduous forest (Mule Hole, India): Regolith matric storage buffers the groundwater recharge process

Accurate estimations of water balance are needed in semi-arid and sub-humid tropical regions, where water resources are scarce compared to water demand. Evapotranspiration plays a major role in this context, and the difficulty to quantify it precisely leads to major uncertainties in the groundwater recharge assessment, especially in forested catchments. In this paper, we propose to assess the importance of deep unsaturated regolith and water uptake by deep tree roots on the groundwater recharge process by using a lumped conceptual model (COMFORT). The model is calibrated using a 5 year hydrological monitoring of an experimental watershed under dry deciduous forest in South India (Mule Hole watershed). The model was able to simulate the stream discharge as well as the contrasted behaviour of groundwater table along the hillslope. Water balance simulated for a 32 year climatic time series displayed a large year-to-year variability, with alternance of dry and wet phases with a time period of approximately 14 years. On an average, input by the rainfall was 1090 mm year(-1) and the evapotranspiration was about 900 mm year(-1) out of which 100 mm year(-1) was uptake from the deep saprolite horizons. The stream flow was 100 mm year(-1) while the groundwater underflow was 80 mm year(-1). The simulation results suggest that (i) deciduous trees can uptake a significant amount of water from the deep regolith, (ii) this uptake, combined with the spatial variability of regolith depth, can account for the variable lag time between drainage events and groundwater rise observed for the different piezometers and (iii) water table response to recharge is buffered due to the long vertical travel time through the deep vadose zone, which constitutes a major water reservoir. This study stresses the importance of long term observations for the understanding of hydrological processes in tropical forested ecosystems. (C) 2009 Elsevier B.V. All rights reserved.

CsrA interacting small RNAs in Haemophilus spp genomes: a theoretical analysis

The csrA is a carbon storage regulator gene that encodes a protein with multiple RNA interaction sites. Bacterial non-coding small RNAs like csrB, csrC and their counterparts in diverse bacterial genus are identified to control the regulatory activities of CsrA and its orthologs. An attempt has been made in this study to identify 'novel' non-coding small RNAs that are involved in the regulatory activities of csrA gene. All CsrA-interacting small RNAs are computationally fingerprinted to have multiple occurrence of 7-nucleotide CsrA interacting repeats [CAGGA(U/A/C)G] along with a 18-nucleotide upstream binding site. However, in several of the genomes like Haemophilus spp, the upstream binding site is not identified. The current methodology overcomes this difficulty by identifying small RNA-specific orphan transcriptional units within the intergenic regions of the genome. The results could identify all known CsrA-interacting small RNAs in E. coli, Vibrio cholerae and Pseudomonas aeruginosa genomes and additionally has picked six new possible CsrA-interacting small RNA regions in E. coli. Our computational analysis indicates that known rygD and rprA sRNAs in E. coli could possibly interact with CsrA proteins. The study is extended to three of the Haemophilus genomes that could identify seven new possible CsrA interacting small RNAs.

Evolutionary Significance of Self-acylation Property in Acyl Carrier Proteins

Acyl carrier protein is an integral component of many cellular metabolic processes. A number of studies have reported self-acylation behavior in acyl carrier proteins. Although AM exhibit high levels of similarity in their primary and tertiary structures, self-acylation behavior is restricted to only some ACPs that can be classified into two major families based on their function. The first family of ACPs is involved in polyketide biosynthesis, whereas the second family participates in fatty acid synthesis. Facilitated by the growing number of genome sequences available for analyses, large-scale phylogenetic studies were used in these studies to uncover as to how self-acylation behavior of acyl carrier proteins is linked with the evolution of metabolic pathways in organisms. These studies show that self-acylation behavior in acyl carrier proteins was lost during the course of evolution, with certain organisms and organelles viz. plastids, retaining it for specified functions. (C) 2009 IUBMB IUBMB Life, 61(8): 853-859, 2009

Inter-helical Interactions in Membrane Proteins: Analysis Based on the Local Backbone Geometry and the Side Chain Interactions

The availability of a significant number of the Structures of helical membrane proteins has prompted us to investigate the mode of helix-helix packing. In the present study, we have considered a dataset of alpha-helical membrane proteins representing Structures solved from all the known superfamilies. We have described the geometry of all the helical residues in terms of local coordinate axis at the backbone level. Significant inter-helical interactions have been considered as contacts by weighing the number of atom-atom contacts, including all the side-chain atoms. Such a definition of local axis and the contact criterion has allowed us to investigate the inter-helical interaction in a systematic and quantitative manner. We show that a single parameter (designated as alpha), which is derived from the parameters representing the Mutual orientation of local axes, is able to accurately Capture the details of helix-helix interaction. The analysis has been carried Out by dividing the dataset into parallel, anti-parallel, and perpendicular orientation of helices. The study indicates that a specific range of alpha value is preferred for interactions among the anti-parallel helices. Such a preference is also seen among interacting residues of parallel helices, however to a lesser extent. No such preference is seen in the case of perpendicular helices, the contacts that arise mainly due to the interaction Of Surface helices with the end of the trans-membrane helices. The Study Supports the prevailing view that the anti-parallel helices are well packed. However, the interactions between helices of parallel orientation are non-trivial. The packing in alpha-helical membrane proteins, which is systematically and rigorously investigated in this study, may prove to be useful in modeling of helical membrane proteins.

Identification and thermodynamic characterization of molten globule states of periplasmic binding proteins

Molten globule-like intermediates have been shown to occur during protein folding and are thought to be involved in protein translocation and membrane insertion. However, the determinants of molten globule stability and the extent of specific packing in molten globules is currently unclear. Using far- and near-UV CD and intrinsic and ANS fluorescence, we show that four periplasmic binding proteins (LBP, LIVBP, MBP, and RBP) form molten globules at acidic pH values ranging from 3.0 to 3.4. Only two of these (LBP and LIVBP) have similar sequences, but all four proteins adopt similar three-dimensional structures. We found that each of the four molten globules binds to its corresponding ligand without conversion to the native state. Ligand binding affinity measured by isothermal titration calorimetry for the molten globule state of LIVBP was found to be comparable to that of the corresponding native state, whereas for LBP, MBP, and RBP, the molten globules bound ligand with approximately 5-30-fold lower affinity than the corresponding native states. All four molten globule states exhibited cooperative thermal unfolding assayed by DSC. Estimated values of Delta C-p of unfolding show that these molten globule states contain 28-67% of buried surface area relative to the native states. The data suggest that molten globules of these periplasmic binding proteins retain a considerable degree of long range order. The ability of these sequentially unrelated proteins to form highly ordered molten globules may be related to their large size as well as an intrinsic property of periplasmic binding protein folds.

Amino acid interaction preferences in proteins

Understanding the key factors that influence the interaction preferences of amino acids in the folding of proteins have remained a challenge. Here we present a knowledge-based approach for determining the effective interactions between amino acids based on amino acid type, their secondary structure, and the contact based environment that they find themselves in the native state structure as measured by their number of neighbors. We find that the optimal information is approximately encoded in a 60 x 60 matrix describing the 20 types of amino acids in three distinct secondary structures (helix, beta strand, and loop). We carry out a clustering scheme to understand the similarity between these interactions and to elucidate a nonredundant set. We demonstrate that the inferred energy parameters can be used for assessing the fit of a given sequence into a putative native state structure.

How Effective Is The Data On Co-Occurrence Of Domains In Multi-Domain Proteins In Prediction Of Protein-Protein Interactions?

We explore the fuse of information on co-occurrence of domains in multi-domain proteins in predicting protein-protein interactions. The basic premise of our work is the assumption that domains co-occurring in a polypeptide chain undergo either structural or functional interactions among themselves. In this study we use a template dataset of domains in multidomain proteins and predict protein-protein interactions in a target organism. We note that maximum number of correct predictions of interacting protein domain families (158) is made in S. cerevisiae when the dataset of closely related organisms is used as the template followed by the more diverse dataset of bacterial proteins (48) and a dataset of randomly chosen proteins (23). We conclude that use of multi-domain information from organisms closely-related to the target can aid prediction of interacting protein families.

Rotavirus nonstructural proteins: a structural perspective

Rotavirus is a major cause of acute infantile diarrhoea worldwide. The virus genome consists of 11 segments of double-stranded RNA that codesfor six structural proteins (VP1-6) and six non-structural proteins(NSP1-6). NSPs are proteins expressed from the virus genome in the infected cell, but are not incorporated into the mature virus article. NSPs play an essential role in virus replication, morphogenesis and pathogenesis, and most of them exhibit multifunctional properties. Structure-function analysis of the NSPs is essential for understanding the molecular mechanisms by which the virus circumvents host innate immune responses, inhibits cellular protein synthesis, hijacks the protein synthetic machinery for its own propagation and manifests the disease process. Because of their essential roles in virus biology, NSPs represent potential targets for the development of antiviral agents. Determination of the three-dimensional structure of NSPs has been hindered due to low-level expression and aggregation. To date, the complete three-dimensional structure of only NSP2 has been determined. The structures of the N- and C-terminal domains of NSP3 and the diarrhoea-inducing domain of NSP4 have also been determined. This review primarily covers the structural and biological functions of the NSPs whose three-dimensional structural aspects have been fully or partially understood, but provides a brief account of other NSPs and the structural features of the mature virion as determined by electron cryomicroscopy.