Structural and Functional Similarity between the Bacterial Type III Secretion System Needle Protein PrgI and the Eukaryotic Apoptosis Bcl-2 Proteins

Background: Functional similarity is challenging to identify when global sequence and structure similarity is low. Activesites or functionally relevant regions are evolutionarily more stable relative to the remainder of a protein structure and provide an alternative means to identify potential functional similarity between proteins. We recently developed the FASTNMR methodology to discover biochemical functions or functional hypotheses of proteins of unknown function by experimentally identifying ligand binding sites. FAST-NMR utilizes our CPASS software and database to assign a function based on a similarity in the structure and sequence of ligand binding sites between proteins of known and unknown function. Methodology/Principal Findings: The PrgI protein from Salmonella typhimurium forms the needle complex in the type III secretion system (T3SS). A FAST-NMR screen identified a similarity between the ligand binding sites of PrgI and the Bcl-2 apoptosis protein Bcl-xL. These ligand binding sites correlate with known protein-protein binding interfaces required for oligomerization. Both proteins form membrane pores through this oligomerization to release effector proteins to stimulate cell death. Structural analysis indicates an overlap between the PrgI structure and the pore forming motif of Bcl-xL. A sequence alignment indicates conservation between the PrgI and Bcl-xL ligand binding sites and pore formation regions. This active-site similarity was then used to verify that chelerythrine, a known Bcl-xL inhibitor, also binds PrgI. Conclusions/Significance: A structural and functional relationship between the bacterial T3SS and eukaryotic apoptosis was identified using our FAST-NMR ligand affinity screen in combination with a bioinformatic analysis based on our CPASS program. A similarity between PrgI and Bcl-xL is not readily apparent using traditional global sequence and structure analysis, but was only identified because of conservation in ligand binding sites. These results demonstrate the unique opportunity that ligand-binding sites provide for the identification of functional relationships when global sequence and structural information is limited

Structure and in Vitro Transcription of Selected Human Transfer RNA Genes

The purpose of this study was to investigate human tRNA gene structure, organization, and expression by isolating and analyzing several human transfer RNA genes

Estimating Protein-Ligand Binding Affinity using High- Throughput Screening by NMR

Many of today’s drug discovery programs utilize high-throughput screening methods that rely on quick evaluations of protein activity to rank potential chemical leads. By monitoring biologically relevant protein-ligand interactions, NMR can provide a means to validate these discovery leads and to optimize the drug discovery process. NMR-based screens typically use a change in chemical shift or linewidth to detect a protein-ligand interaction. However, the relatively low throughput of current NMR screens and their high demand on sample requirements generally makes it impractical to collect complete binding curves to measure the affinity for each compound in a large and diverse chemical library. As a result, NMR ligand screens are typically limited to identifying candidates that bind to a protein and do not give any estimate of the binding affinity. To address this issue, a methodology has been developed to rank binding affinities for ligands based on NMR-based screens that use 1D 1H NMR line-broadening experiments. This method was demonstrated by using it to estimate the dissociation equilibrium constants for twelve ligands with the protein human serum albumin (HSA). The results were found to give good agreement with previous affinities that have been reported for these same ligands with HSA

Bacterial Protein Structures Reveal Phylum Dependent Divergence

Protein sequence space is vast compared to protein fold space. This raises important questions about how structures adapt to evolutionary changes in protein sequences. A growing trend is to regard protein fold space as a continuum rather than a series of discrete structures. From this perspective, homologous protein structures within the same functional classification should reveal a constant rate of structural drift relative to sequence changes. The clusters of orthologous groups (COG) classification system was used to annotate homologous bacterial protein structures in the Protein Data Bank (PDB). The structures and sequences of proteins within each COG were compared against each other to establish their relatedness. As expected, the analysis demonstrates a sharp structural divergence between the bacterial phyla Firmicutes and Proteobacteria. Additionally, each COG had a distinct sequence/structure relationship, indicating that different evolutionary pressures affect the degree of structural divergence. However, our analysis also shows the relative drift rate between sequence identity and structure divergence remains constant

Estimating Protein-Ligand Binding Affinity using High- Throughput Screening by NMR

Many of today’s drug discovery programs utilize high-throughput screening methods that rely on quick evaluations of protein activity to rank potential chemical leads. By monitoring biologically relevant protein-ligand interactions, NMR can provide a means to validate these discovery leads and to optimize the drug discovery process. NMR-based screens typically use a change in chemical shift or linewidth to detect a protein-ligand interaction. However, the relatively low throughput of current NMR screens and their high demand on sample requirements generally makes it impractical to collect complete binding curves to measure the affinity for each compound in a large and diverse chemical library. As a result, NMR ligand screens are typically limited to identifying candidates that bind to a protein and do not give any estimate of the binding affinity. To address this issue, a methodology has been developed to rank binding affinities for ligands based on NMR-based screens that use 1D 1H NMR line-broadening experiments. This method was demonstrated by using it to estimate the dissociation equilibrium constants for twelve ligands with the protein human serum albumin (HSA). The results were found to give good agreement with previous affinities that have been reported for these same ligands with HSA

A Correlation between Protein Function and Ligand Binding Profiles

We report that proteins with the same function bind the same set of small molecules from a standardized chemical library. This observation led to a quantifiable and rapidly adaptable method for protein functional analysis using experimentally-derived ligand binding profiles. Ligand binding is measured using a high-throughput NMR ligand affinity screen with a structurally diverse chemical library. The method was demonstrated using a set of 19 proteins with a range of functions. A statistically significant similarity in ligand binding profiles was only observed between the two functionally identical albumins and between the five functionally similar amylases. This new approach is independent of sequence, structure or evolutionary information, and therefore, extends our ability to analyze and functionally annotate novel genes

Analysis of Metabolomic PCA Data using Tree Diagrams

Large amounts of data from high throughput metabolomic experiments are commonly visualized using a principal component analysis (PCA) 2D scores plot. The question of the similarity or difference between multiple metabolic states then becomes a question of the degree of overlap between their respective data point clusters in PC scores space. A qualitative visual inspection of the clustering pattern in PCA score plots is a common protocol. This report describes the application of tree diagrams and bootstrapping techniques for an improved quantitative analysis of metabolic PCA data clustering. Our PCAtoTree program creates a distance matrix with 100 bootstrap steps that describes the separation of all clusters in a metabolic dataset. Using accepted phylogenetic software, the distance matrix resulting from the various metabolic states is organized into a phylogenetic-like tree format, where bootstrap values ≥ 50 indicate a statistically relevant branch separation. PCAtoTree analysis of two previously published data sets demonstrates the improved resolution of metabolic state differences using tree diagrams. In addition, for metabolomic studies of large numbers of different metabolic states, the tree format provides a better description of similarities and differences between each metabolic state. The approach is also tolerant of sample size variations between different metabolic states

Analysis of Metabolomic PCA Data using Tree Diagrams

Large amounts of data from high throughput metabolomic experiments are commonly visualized using a principal component analysis (PCA) 2D scores plot. The question of the similarity or difference between multiple metabolic states then becomes a question of the degree of overlap between their respective data point clusters in PC scores space. A qualitative visual inspection of the clustering pattern in PCA score plots is a common protocol. This report describes the application of tree diagrams and bootstrapping techniques for an improved quantitative analysis of metabolic PCA data clustering. Our PCAtoTree program creates a distance matrix with 100 bootstrap steps that describes the separation of all clusters in a metabolic dataset. Using accepted phylogenetic software, the distance matrix resulting from the various metabolic states is organized into a phylogenetic-like tree format, where bootstrap values ≥ 50 indicate a statistically relevant branch separation. PCAtoTree analysis of two previously published data sets demonstrates the improved resolution of metabolic state differences using tree diagrams. In addition, for metabolomic studies of large numbers of different metabolic states, the tree format provides a better description of similarities and differences between each metabolic state. The approach is also tolerant of sample size variations between different metabolic states

PROFESS: a PROtein Function, Evolution, Structure and Sequence database

The proliferation of biological databases and the easy access enabled by the Internet is having a beneficial impact on biological sciences and transforming the way research is conducted. There are ∼1100 molecular biology databases dispersed throughout the Internet. To assist in the functional, structural and evolutionary analysis of the abundant number of novel proteins continually identified from whole-genome sequencing, we introduce the PROFESS (PROtein Function, Evolution, Structure and Sequence) database. Our database is designed to be versatile and expandable and will not confine analysis to a pre-existing set of data relationships. A fundamental component of this approach is the development of an intuitive query system that incorporates a variety of similarity functions capable of generating data relationships not conceived during the creation of the database. The utility of PROFESS is demonstrated by the analysis of the structural drift of homologous proteins and the identification of potential pancreatic cancer therapeutic targets based on the observation of protein–protein interaction networks