Virtual high-throughput ligand screening.

In Structural Genomics projects, virtual high-throughput ligand screening can be utilized to provide important functional details for newly determined protein structures. Using a variety of publicly available software tools, it is possible to computationally model, predict, and evaluate how different ligands interact with a given protein. At the Center for Structural Genomics of Infectious Diseases (CSGID) a series of protein analysis, docking and molecular dynamics software is scripted into a single hierarchical pipeline allowing for an exhaustive investigation of protein-ligand interactions. The ability to conduct accurate computational predictions of protein-ligand binding is a vital component in improving both the efficiency and economics of drug discovery. Computational simulations can minimize experimental efforts, the slowest and most cost prohibitive aspect of identifying new therapeutics.

Sensor domain of histidine kinase KinB of Pseudomonas: a helix-swapped dimer.

The overproduction of polysaccharide alginate is responsible for the formation of mucus in the lungs of cystic fibrosis patients. Histidine kinase KinB of the KinB-AlgB two-component system in Pseudomonas aeruginosa acts as a negative regulator of alginate biosynthesis. The modular architecture of KinB is similar to other histidine kinases. However, its periplasmic signal sensor domain is unique and is found only in the Pseudomonas genus. Here, we present the first crystal structures of the KinB sensor domain. The domain is a dimer in solution, and in the crystal it shows an atypical dimer of a helix-swapped four-helix bundle. A positively charged cavity is formed on the dimer interface and involves several strictly conserved residues, including Arg-60. A phosphate anion is bound asymmetrically in one of the structures. In silico docking identified several monophosphorylated sugars, including ß-D-fructose 6-phosphate and ß-D-mannose 6-phosphate, a precursor and an intermediate of alginate synthesis, respectively, as potential KinB ligands. Ligand binding was confirmed experimentally. Conformational transition from a symmetric to an asymmetric structure and decreasing dimer stability caused by ligand binding may be a part of the signal transduction mechanism of the KinB-AlgB two-component system.

Predicting HLA class I non-permissive amino acid residues substitutions.

Prediction of peptide binding to human leukocyte antigen (HLA) molecules is essential to a wide range of clinical entities from vaccine design to stem cell transplant compatibility. Here we present a new structure-based methodology that applies robust computational tools to model peptide-HLA (p-HLA) binding interactions. The method leverages the structural conservation observed in p-HLA complexes to significantly reduce the search space and calculate the system's binding free energy. This approach is benchmarked against existing p-HLA complexes and the prediction performance is measured against a library of experimentally validated peptides. The effect on binding activity across a large set of high-affinity peptides is used to investigate amino acid mismatches reported as high-risk factors in hematopoietic stem cell transplantation.

Bacillus anthracis inosine 5'-monophosphate dehydrogenase in action: the first bacterial series of structures of phosphate ion-, substrate-, and product-bound complexes.

Inosine 5'-monophosphate dehydrogenase (IMPDH) catalyzes the first unique step of the GMP branch of the purine nucleotide biosynthetic pathway. This enzyme is found in organisms of all three kingdoms. IMPDH inhibitors have broad clinical applications in cancer treatment, as antiviral drugs and as immunosuppressants, and have also displayed antibiotic activity. We have determined three crystal structures of Bacillus anthracis IMPDH, in a phosphate ion-bound (termed "apo") form and in complex with its substrate, inosine 5'-monophosphate (IMP), and product, xanthosine 5'-monophosphate (XMP). This is the first example of a bacterial IMPDH in more than one state from the same organism. Furthermore, for the first time for a prokaryotic enzyme, the entire active site flap, containing the conserved Arg-Tyr dyad, is clearly visible in the structure of the apoenzyme. Kinetic parameters for the enzymatic reaction were also determined, and the inhibitory effect of XMP and mycophenolic acid (MPA) has been studied. In addition, the inhibitory potential of two known Cryptosporidium parvum IMPDH inhibitors was examined for the B. anthracis enzyme and compared with those of three bacterial IMPDHs from Campylobacter jejuni, Clostridium perfringens, and Vibrio cholerae. The structures contribute to the characterization of the active site and design of inhibitors that specifically target B. anthracis and other microbial IMPDH enzymes.

Structure of transcription factor HetR required for heterocyst differentiation in cyanobacteria.

HetR is an essential regulator of heterocyst development in cyanobacteria. HetR binds to a DNA palindrome upstream of the hetP gene. We report the crystal structure of HetR from Fischerella at 3.0 Å. The protein is a dimer comprised of a central DNA-binding unit containing the N-terminal regions of the two subunits organized with two helix-turn-helix motifs; two globular flaps extending in opposite directions; and a hood over the central core formed from the C-terminal subdomains. The flaps and hood have no structural precedent in the protein database, therefore representing new folds. The structural assignments are supported by site-directed mutagenesis and DNA-binding studies. We suggest that HetR serves as a scaffold for assembly of transcription components critical for heterocyst development.

Structure of surface layer homology (SLH) domains from Bacillus anthracis surface array protein.

Surface (S)-layers, para-crystalline arrays of protein, are deposited in the envelope of most bacterial species. These surface organelles are retained in the bacterial envelope through the non-covalent association of proteins with cell wall carbohydrates. Bacillus anthracis, a Gram-positive pathogen, produces S-layers of the protein Sap, which uses three consecutive repeats of the surface-layer homology (SLH) domain to engage secondary cell wall polysaccharides (SCWP). Using x-ray crystallography, we reveal here the structure of these SLH domains, which assume the shape of a three-prong spindle. Each SLH domain contributes to a three-helical bundle at the spindle base, whereas another α-helix and its connecting loops generate the three prongs. The inter-prong grooves contain conserved cationic and anionic residues, which are necessary for SLH domains to bind the B. anthracis SCWP. Modeling experiments suggest that the SLH domains of other S-layer proteins also fold into three-prong spindles and capture bacterial envelope carbohydrates by a similar mechanism.

Assisted assignment of ligands corresponding to unknown electron density.

A semi-automated computational procedure to assist in the identification of bound ligands from unknown electron density has been developed. The atomic surface surrounding the density blob is compared to a library of three-dimensional ligand binding surfaces extracted from the Protein Data Bank (PDB). Ligands corresponding to surfaces which share physicochemical texture and geometric shape similarities are considered for assignment. The method is benchmarked against a set of well represented ligands from the PDB, in which we show that we can identify the correct ligand based on the corresponding binding surface. Finally, we apply the method during model building and refinement stages from structural genomics targets in which unknown density blobs were discovered. A semi-automated computational method is described which aims to assist crystallographers with assigning the identity of a ligand corresponding to unknown electron density. Using shape and physicochemical similarity assessments between the protein surface surrounding the density and a database of known ligand binding surfaces, a plausible list of candidate ligands are identified for consideration. The method is validated against highly observed ligands from the Protein Data Bank and results are shown from its use in a high-throughput structural genomics pipeline.

Protein functional surfaces: global shape matching and local spatial alignments of ligand binding sites.

BACKGROUND: Protein surfaces comprise only a fraction of the total residues but are the most conserved functional features of proteins. Surfaces performing identical functions are found in proteins absent of any sequence or fold similarity. While biochemical activity can be attributed to a few key residues, the broader surrounding environment plays an equally important role. RESULTS: We describe a methodology that attempts to optimize two components, global shape and local physicochemical texture, for evaluating the similarity between a pair of surfaces. Surface shape similarity is assessed using a three-dimensional object recognition algorithm and physicochemical texture similarity is assessed through a spatial alignment of conserved residues between the surfaces. The comparisons are used in tandem to efficiently search the Global Protein Surface Survey (GPSS), a library of annotated surfaces derived from structures in the PDB, for studying evolutionary relationships and uncovering novel similarities between proteins. CONCLUSION: We provide an assessment of our method using library retrieval experiments for identifying functionally homologous surfaces binding different ligands, functionally diverse surfaces binding the same ligand, and binding surfaces of ubiquitous and conformationally flexible ligands. Results using surface similarity to predict function for proteins of unknown function are reported. Additionally, an automated analysis of the ATP binding surface landscape is presented to provide insight into the correlation between surface similarity and function for structures in the PDB and for the subset of protein kinases.

The 1.38 A crystal structure of DmsD protein from Salmonella typhimurium, a proofreading chaperone on the Tat pathway.

The DmsD protein is necessary for the biogenesis of dimethyl sulphoxide (DMSO) reductase in many prokaryotes. It performs a critical chaperone function initiated through its binding to the twin-arginine signal peptide of DmsA, the catalytic subunit of DMSO reductase. Upon binding to DmsD, DmsA is translocated to the periplasm via the so-called twin-arginine translocation (Tat) pathway. Here we report the 1.38 A crystal structure of the protein DmsD from Salmonella typhimurium and compare it with a close functional homolog, TorD. DmsD has an all-alpha fold structure with a notable helical extension located at its N-terminus with two solvent exposed hydrophobic residues. A major difference between DmsD and TorD is that TorD structure is a domain-swapped dimer, while DmsD exists as a monomer. Nevertheless, these two proteins have a number of common features suggesting they function by using similar mechanisms. A possible signal peptide-binding site is proposed based on structural similarities. Computational analysis was used to identify a potential GTP binding pocket on similar surfaces of DmsD and TorD structures.

Chapter 4. Predicting and characterizing protein functions through matching geometric and evolutionary patterns of binding surfaces.

Predicting protein functions from structures is an important and challenging task. Although proteins are often thought to be packed as tightly as solids, closer examination based on geometric computation reveals that they contain numerous voids and pockets. Most of them are of random nature, but some are binding sites providing surfaces to interact with other molecules. A promising approach for function prediction is to infer functions through discovery of similarity in local binding pockets, as proteins binding to similar substrates/ligands and carrying out similar functions have similar physical constraints for binding and reactions. In this chapter, we describe computational methods to distinguish those surface pockets that are likely to be involved in important biological functions, and methods to identify key residues in these pockets. We further describe how to predict protein functions at large scale from structures by detecting binding surfaces similar in residue make-ups, shape, and orientation. We also describe a Bayesian Monte Carlo method that can separate selection pressure due to biological function from pressure due to protein folding. We show how this method can be used to reconstruct the evolutionary history of binding surfaces for detecting similar binding surfaces. In addition, we briefly discuss how the negative image of a binding pocket can be casted, and how such information can be used to facilitate drug discovery.

Protein surface analysis for function annotation in high-throughput structural genomics pipeline.

Structural genomics (SG) initiatives are expanding the universe of protein fold space by rapidly determining structures of proteins that were intentionally selected on the basis of low sequence similarity to proteins of known structure. Often these proteins have no associated biochemical or cellular functions. The SG success has resulted in an accelerated deposition of novel structures. In some cases the structural bioinformatics analysis applied to these novel structures has provided specific functional assignment. However, this approach has also uncovered limitations in the functional analysis of uncharacterized proteins using traditional sequence and backbone structure methodologies. A novel method, named pvSOAR (pocket and void Surface of Amino Acid Residues), of comparing the protein surfaces of geometrically defined pockets and voids was developed. pvSOAR was able to detect previously unrecognized and novel functional relationships between surface features of proteins. In this study, pvSOAR is applied to several structural genomics proteins. We examined the surfaces of YecM, BioH, and RpiB from Escherichia coli as well as the CBS domains from inosine-5'-monosphate dehydrogenase from Streptococcus pyogenes, conserved hypothetical protein Ta549 from Thermoplasm acidophilum, and CBS domain protein mt1622 from Methanobacterium thermoautotrophicum with the goal to infer information about their biochemical function.

pvSOAR: detecting similar surface patterns of pocket and void surfaces of amino acid residues on proteins.

Detecting similar protein surfaces provides an important route for discovering unrecognized or novel functional relationship between proteins. The web server pvSOAR (pocket and void Surfaces Of Amino acid Residues) provides an online resource to identify similar protein surface regions. pvSOAR can take a structure either uploaded by a user or obtained from the Protein Data Bank, and identifies similar surface patterns based on geometrically defined pockets and voids. It provides several search modes to compare protein surfaces by similarity in local sequence, local shape and local orientation. Statistically significant search results are reported for visualization and interactive exploration. pvSOAR can be used to predict biological functions of proteins with known three-dimensional structures but unknown biological roles. It can also be used to study functional relationship between proteins and for exploration of the evolutionary origins of structural elements important for protein function. We present an example using pvSOAR to explore the biological roles of a protein whose structure was solved by the structural genomics project. The pvSOAR web server is available at http://pvsoar.bioengr.uic.edu/.

topoSNP: a topographic database of non-synonymous single nucleotide polymorphisms with and without known disease association.

The database of topographic mapping of Single Nucleotide Polymorphism (topoSNP) provides an online resource for analyzing non-synonymous SNPs (nsSNPs) that can be mapped onto known 3D structures of proteins. These include disease- associated nsSNPs derived from the Online Mendelian Inheritance in Man (OMIM) database and other nsSNPs derived from dbSNP, a resource at the National Center for Biotechnology Information that catalogs SNPs. TopoSNP further classifies each nsSNP site into three categories based on their geometric location: those located in a surface pocket or an interior void of the protein, those on a convex region or a shallow depressed region, and those that are completely buried in the interior of the protein structure. These unique geometric descriptions provide more detailed mapping of nsSNPs to protein structures. The current release also includes relative entropy of SNPs calculated from multiple sequence alignment as obtained from the Pfam database (a database of protein families and conserved protein motifs) as well as manually adjusted multiple alignments obtained from ClustalW. These structural and conservational data can be useful for studying whether nsSNPs in coding regions are likely to lead to phenotypic changes. TopoSNP includes an interactive structural visualization web interface, as well as downloadable batch data. The database will be updated at regular intervals and can be accessed at: http://gila.bioengr.uic.edu/snp/toposnp.

Order independent structural alignment of circularly permuted proteins.

Circular permutation connects the N and C termini of a protein and concurrently cleaves elsewhere in the chain, providing an important mechanism for generating novel protein fold and functions. However, their in genomes is unknown because current detection methods can miss many occurrences, mistaking random repeats as circular permutation. Here we develop a method for detecting circularly permuted proteins from structural comparison. Sequence order independent alignment of protein structures can be regarded as a special case of the maximum-weight independent set problem, which is known to be computationally hard. We develop an efficient approximation algorithm by repeatedly solving relaxations of an appropriate intermediate integer programming formulation, we show that the approximation ratio is much better than the theoretical worst case ratio of r=1/4. Circularly permuted proteins reported in literature can be identified rapidly with our method, while they escape the detection by publicly available servers for structural alignment.

Inferring functional relationships of proteins from local sequence and spatial surface patterns.

We describe a novel approach for inferring functional relationship of proteins by detecting sequence and spatial patterns of protein surfaces. Well-formed concave surface regions in the form of pockets and voids are examined to identify similarity relationship that might be directly related to protein function. We first exhaustively identify and measure analytically all 910,379 surface pockets and interior voids on 12,177 protein structures from the Protein Data Bank. The similarity of patterns of residues forming pockets and voids are then assessed in sequence, in spatial arrangement, and in orientational arrangement. Statistical significance in the form of E and p-values is then estimated for each of the three types of similarity measurements. Our method is fully automated without human intervention and can be used without input of query patterns. It does not assume any prior knowledge of functional residues of a protein, and can detect similarity based on surface patterns small and large. It also tolerates, to some extent, conformational flexibility of functional sites. We show with examples that this method can detect functional relationship with specificity for members of the same protein family and superfamily, as well as remotely related functional surfaces from proteins of different fold structures. We envision that this method can be used for discovering novel functional relationship of protein surfaces, for functional annotation of protein structures with unknown biological roles, and for further inquiries on evolutionary origins of structural elements important for protein function.

CASTp: Computed Atlas of Surface Topography of proteins.

Computed Atlas of Surface Topography of proteins (CASTp) provides an online resource for locating, delineating and measuring concave surface regions on three-dimensional structures of proteins. These include pockets located on protein surfaces and voids buried in the interior of proteins. The measurement includes the area and volume of pocket or void by solvent accessible surface model (Richards' surface) and by molecular surface model (Connolly's surface), all calculated analytically. CASTp can be used to study surface features and functional regions of proteins. CASTp includes a graphical user interface, flexible interactive visualization, as well as on-the-fly calculation for user uploaded structures. CASTp is updated daily and can be accessed at http://cast.engr.uic.edu.

Higher-order interhelical spatial interactions in membrane proteins.

Higher-order interactions are important for protein folding and assembly. We introduce the concept of interhelical three-body interactions as derived from Delaunay triangulation and alpha shapes of protein structures. In addition to glycophorin A, where triplets are strongly correlated with protein stability, we found that tight interhelical triplet interactions exist extensively in other membrane proteins, where many types of triplets occur far more frequently than in soluble proteins. We developed a probabilistic model for estimating the value of membrane helical interaction triplet (MHIT) propensity. Because the number of known structures of membrane proteins is limited, we developed a bootstrap method for determining the 95% confidence intervals of estimated MHIT values. We identified triplets that have high propensity for interhelical interactions and are unique to membrane proteins, e.g. AGF, AGG, GLL, GFF and others. A significant fraction (32%) of triplet types contains triplets that may be involved in interhelical hydrogen bond interactions, suggesting the prevalent and important roles of H-bond in the assembly of TM helices. There are several well-defined spatial conformations for triplet interactions on helices with similar parallel or antiparallel orientations and with similar right-handed or left-handed crossing angles. Often, they contain small residues and correspond to the regions of the closest contact between helices. Sequence motifs such as GG4 and AG4 can be part of the three-body interactions that have similar conformations, which in turn can be part of a higher-order cooperative four residue spatial motif observed in helical pairs from different proteins. In many cases, spatial motifs such as serine zipper and polar clamp are part of triplet interactions. On the basis of the analysis of the archaeal rhodopsin family of proteins, tightly packed triplet interactions can be achieved with several different choices of amino acid residues.