Evidence of Conformational Selection Driving the Formation of Ligand Binding Sites in Protein-Protein Interfaces

Many protein-protein interactions (PPIs) are compelling targets for drug discovery, and in a number of cases can be disrupted by small molecules. The main goal of this study is to examine the mechanism of binding site formation in the interface region of proteins that are PPI targets by comparing ligand-free and ligand-bound structures. To avoid any potential bias, we focus on ensembles of ligand-free protein conformations obtained by nuclear magnetic resonance (NMR) techniques and deposited in the Protein Data Bank, rather than on ensembles specifically generated for this study. The measures used for structure comparison are based on detecting binding hot spots, i.e., protein regions that are major contributors to the binding free energy. The main tool of the analysis is computational solvent mapping, which explores the surface of proteins by docking a large number of small “probe” molecules. Although we consider conformational ensembles obtained by NMR techniques, the analysis is independent of the method used for generating the structures. Finding the energetically most important regions, mapping can identify binding site residues using ligand-free models based on NMR data. In addition, the method selects conformations that are similar to some peptide-bound or ligand-bound structure in terms of the properties of the binding site. This agrees with the conformational selection model of molecular recognition, which assumes such pre-existing conformations. The analysis also shows the maximum level of similarity between unbound and bound states that is achieved without any influence from a ligand. Further shift toward the bound structure assumes protein-peptide or protein-ligand interactions, either selecting higher energy conformations that are not part of the NMR ensemble, or leading to induced fit. Thus, forming the sites in protein-protein interfaces that bind peptides and can be targeted by small ligands always includes conformational selection, although other recognition mechanisms may also be involved.     

Homology modeling and assigned functional annotation of an uncharacterized antitoxin protein from Streptomyces xinghaiensis

Streptomyces xinghaiensis is a Gram-positive, aerobic and non-motile bacterium. The bacterial genome is known. Therefore, it is of interest to study the uncharacterized proteins in the genome. An uncharacterized protein (gi|518540893|86 residues) in the genome was selected for a comprehensive computational sequence-structure-function analysis using available data and tools. Subcellular localization of the targeted protein with conserved residues and assigned secondary structures is documented. Sequence homology search against the protein data bank (PDB) and non-redundant GenBank proteins using BLASTp showed different homologous proteins with known antitoxin function. A homology model of the target protein was developed using a known template (PDB ID: 3CTO:A) with 62% sequence similarity in HHpred after assessment using programs PROCHECK and QMEAN6. The predicted active site using CASTp is analyzed for assigned anti-toxin function. This information finds specific utility in annotating the said uncharacterized protein in the bacterial genome.     

脊髓灰质炎病毒各代表株壳蛋白二级结构及亲水概率的电子计算机预测

本文用脊髓灰质炎病毒3个型6个强弱代表株壳蛋白一级结构,借助电子计算机预测和计算出病毒壳蛋白的二级结构和亲水概率。     

Antheprot: An interactive graphics software for analyzing protein structures from sequences

ANTHEPROT is a fully interactive program devoted to the analysis of protein structures using a graphics workstation. It presents four options: The first option can predict secondary structures using five methods, and hydrophobicity, solvent accessibility, flexibility and antigenicity profiles using eighteen scales. The user may introduce his own scales. The results displayed on the screen can be easily analyzed. The second option is for representing results concerning up to eight proteins by one method. To compare these proteins, it is possible to align the profiles or the predicted secondary structure according to various motifs. The secondary structure deduced from crystallographic data may also be introduced. The third option is designed to compare the primary structure of two proteins and to visualize on the screen regions that exhibit similarity. Six different comparison matrices may be used, but the user can also introduce his own matrices. The last option is for studying the proteolytic peptides resulting from a chemical or enzymatic digestion of a given protein. It is possible to analyze the protein cleavage using eleven chemical reagents or enzymes. The results are displayed on the screen as RP-HPLC chromatogram.     

In silico analysis of squalene synthase in Fabaceae family using bioinformatics tools

Triterpenoid saponins are a diverse group of bioactive compounds, which are used for possessing of many biomedical and pharmaceutical products. Generally, squalene synthase (SQS) is defined as an emerging and essential branch point enzyme far from the major pathway of isoprenoids biosynthetic and a latent adjusting point, which manages carbon flux into triterpenes biosynthesis and sterols. The present study deals with the detailed characterization of SQS by bioinformatics approaches to evaluate physicochemical properties, structural characteristics including secondary and 3D structure prediction and functional analysis from eight plants related to Fabaceae family and Arabidopsis thaliana. Bioinformatics analysis revealed that SQS proteins have two transmembrane regions in the C-terminal. The predicted motifs were used to design universal degenerate primers for PCR analysis and other molecular applications. Phylogenetic analysis showed conserved regions at different stretches with maximum homology in amino acid residues within all SQSs. The secondary structure prediction results showed that the amino acid sequence of all squalene synthases had α helix and random coil as the main components. The reliability of the received model was confirmed using the ProSA and RAMPAGE programs. Determining of active site by CASTp proposes the possibility of using this protein as probable medication target. The findings of the present study may be useful for further assessments on characterization and cloning of squalene synthase.     

Three-Dimensional Models of Immune Cell Surface Proteins and Identification of Binding Sites

The extracellular regions of many cell surface proteins of the immune system contain distinct domains that may be linked in many different ways and are often only loosely tethered to the transmembrane segment. In efforts to identify regions critical for binding, molecular models of these domains are used to select residues for mutagenesis and to map binding sites. Many immune cell surface proteins belong to protein superfamilies and display only limited sequence identity compared to proteins of known three-dimensional (3D) structure, often 30% or less. Therefore, detailed 3D structures are difficult to predict, and structure-based sequence analysis and model assessment are particularly important components of the model building process. In some cases, experimentally determined structures have made it possible to assess the accuracy of predictions, which illustrates the opportunities and shortcomings of the approach. Herein the model-based identification of binding sites in cell surface proteins is described and representative examples are discussed.     

DNA-hybridization electron microscopy tertiary structure of 16 S rRNA

Seven regions of 16 S rRNA have been located on the surface of the 30 S ribosomal subunit by DNA-hybridization electron microscopy. This information has been incorporated into a model for the tertiary structure of 16 S rRNA, accounting for approximately 40% of the total 16 S rRNA. A structure labeled the platform ring is proposed for a region of rRNA within the central domain. This structure rings the edges of the platform and includes regions 655-751 and 769-810. Another region, the recognition complex, consists of nucleotides 500 to 545, and occupies a region on the exterior surface of the subunit near the elongation factor Tu binding site. Ribosomal proteins that have been mapped by immunoelectron microscopy are superimposed onto the model in order to examine possible regions of interaction. Good correlation between the model locations of ribosomal proteins, and regions of rRNA protected by ribosomal proteins provide independent support for this model.     

Predicted metal binding sites for phytoremediation

Metal ion binding domains are found in proteins that mediate transport, buffering or detoxification of metal ions. The objective of the study is to design and analyze metal binding motifs against the genes involved in phytoremediation. This is being done on the basis of certain pre-requisite amino-acid residues known to bind metal ions/metal complexes in medicinal and aromatic plants (MAP''s). Earlier work on MAP''s have shown that heavy metals accumulated by aromatic and medicinal plants do not appear in the essential oil and that some of these species are able to grow in metal contaminated sites. A pattern search against the UniProtKB/Swiss-Prot and UniProtKB/TrEMBL databases yielded true positives in each case showing the high specificity of the motifs designed for the ions of nickel, lead, molybdenum, manganese, cadmium, zinc, iron, cobalt and xenobiotic compounds. Motifs were also studied against PDB structures. Results of the study suggested the presence of binding sites on the surface of protein molecules involved. PDB structures of proteins were finally predicted for the binding sites functionality in their respective phytoremediation usage. This was further validated through CASTp server to study its physico-chemical properties. Bioinformatics implications would help in designing strategy for developing transgenic plants with increased metal binding capacity. These metal binding factors can be used to restrict metal update by plants. This helps in reducing the possibility of metal movement into the food chain.     

Predicting gene ontology functions from protein's regional surface structures

Background  

Annotation of protein functions is an important task in the post-genomic era. Most early approaches for this task exploit only the sequence or global structure information. However, protein surfaces are believed to be crucial to protein functions because they are the main interfaces to facilitate biological interactions. Recently, several databases related to structural surfaces, such as pockets and cavities, have been constructed with a comprehensive library of identified surface structures. For example, CASTp provides identification and measurements of surface accessible pockets as well as interior inaccessible cavities.     

PROFALIGN algorithm identifies the regions containing folding determinants by scoring pairs of hydrophobic profiles of remotely related proteins.

Profile comparison methods have been shown to be very powerful in creating accurate alignments of protein sequences, especially in the case of remotely related proteins (RRP). These methods take advantage of the observation that hydrophobic profiles are more conserved than the corresponding amino acid sequences. Here, we present the PROFALIGN algorithm, which allows one to perform a detailed comparative analysis, at both local and global levels of two protein sequence profiles. The user can either choose among four different hydrophobic scales (Miyazawa-Jernigan, Eisenberg, Engelman-Steiz, and Kyte-Doolittle) or can add a personal scale. The interface is designed for a wide range of users, including those who are not involved in protein research. It allows one to vary the alignment parameters (such as gap penalties, embedding, and profile smoothness). Secondary structure propensity is added as an optional alignment filter. Similar segments of two proteins are singled out on the basis of score. We have tested the algorithm with different Src homology 3 (SH3) domain fragments sharing low sequence homology but very similar three-dimensional (3D) structures. By using the Miyazawa-Jernigan hydrophobic scale, PROFALIGN was able to detect the strong correlation between the regions that are known to be crucial for SH3 transition state topology. PROFALIGN seems able to identify most of the mutual alignment of structures on the basis of their hydrophobic profiles, delimiting the regions containing the key determinants of folding. Therefore, the present methodology may be useful for the detection of the most structurally relevant positions inside remote related proteins.     

Structure of C-terminal fragment of merozoite surface protein-1 from Plasmodium vivax determined by homology modeling and molecular dynamics refinement

One current vaccine candidate against Plasmodium vivax targeting asexual blood stage is the major merozoite surface protein-1 of P. vivax (PvMSP-1). Vaccine trials with PvMSP-1₁₉ and PvMSP-1₃₃ have succeeded in protecting monkeys and a large proportion of individuals, naturally exposed to P. vivax transmission, develop specific antibodies to PvMSP-1₁₉. This study presents a model for the three-dimensional structure of the C-terminal 19 kDa fragment of P. vivax MSP-1 determined by means of homology modeling and molecular dynamics refinement. The structure proved to be consistent with MSP-1₁₉ of known crystal or solution structures. The presence of a main binding pocket, well suited for protein–protein interactions, was determined by CASTp. Corrections reported to the sequence of PvMSP-1₁₉ Belem strain were also inspected. Our model is currently used as a basis to understand antibody interactions with PvMSP-1₁₉.     

Development of rules for vaccine engineering based on the variability of peptide fragments in closely related proteins

Based on the protein sequence data bank (PIR), the "variable fragment" bank, comprising pairs of closely related proteins, containing one or more strongly differing sites of primary structures was formed. The bank includes 465 "variable fragments" of 383 protein pairs. Amino acid residues composition of "variable fragments" was examined and indexes of potential amino acid residues variability was formed. An analysis of amino acid fragments replaceability was carried out by substituting the N-, C-terminal, or middle part of a chain), the fragments length differences and physico-chemical properties of residues, such as volume, hydrophobicity, polarity, isoelectric point, etc. Some general empirical rules of peptide insertions in carrier-proteins were created based on these analyses. The rules are directed for performing modifications maintaining the common structure and function of the carrier-protein molecule. The selection scheme for determining the regions suitable for modification and the criteria for defining the width of acceptable modifications in this regions were suggested. The use of potential variability profile for detecting regions suitable for peptide insertion was considered on the model of hepatitis B surface protein.     

The chromosome periphery during mitosis

A complex structure, visible by electron microscopy, surrounds each chromosome during mitosis. The organization of this structure is distinct from that of the chromosomes and the cytoplasm. It forms a perichromosomal layer that can be isolated together with the chromosomes. This layer covers the chromosomes except in centromeric regions. The perichromosomal layer includes nuclear and nucleolar proteins as well as ribonucleoproteins (RNPs). The list of proteins and RNAs identified includes nuclear matrix proteins (perichromin, peripherin), nucleolar proteins (perichro-monucleolin, Ki-67 antigen, B23 protein, fibrillarin, p103, p52), ribosomal proteins (S1) and snRNAs (U3 RNAs). Only limited information is available about how and when the perichromosomal layer is formed. During early prophase, the proteins extend from the nucleoli towards the periphery of the nucleus. Thin cordon-like structures reach the nuclear envelope delimiting areas in which chromosomes condense. At telophase, the proteins are associated with the part of the chromosomes remaining condensed and accumulate in newly formed nucleoli in regions where chromatin is already decondensed. The perichromosomal layer contains several different classes of proteins and RNPs and it has been attributed various roles: (1) in chromosome organization, (2) as a barrier around the chromosomes, (3) involvement in compartmentation of the cells in prophase and telophase and (4) a binding site for chromosomal passenger proteins necessary to the early process of nuclear assembly.     

Intrinsic disorder in the Protein Data Bank

The Protein Data Bank (PDB) is the preeminent source of protein structural information. PDB contains over 32,500 experimentally determined 3-D structures solved using X-ray crystallography or nuclear magnetic resonance spectroscopy. Intrinsically disordered regions fail to form a fixed 3-D structure under physiological conditions. In this study, we compare the amino-acid sequences of proteins whose structures are determined by X-ray crystallography with the corresponding sequences from the Swiss-Prot database. The analyzed dataset includes 16,370 structures, which represent 18,101 PDB chains and 5,434 different proteins from 910 different organisms (2,793 eukaryotic, 2,109 bacterial, 288 viral, and 244 archaeal). In this dataset, on average, each Swiss-Prot protein is represented by 7 PDB chains with 76% of the crystallized regions being represented by more than one structure. Intriguingly, the complete sequences of only approximately 7% of proteins are observed in the corresponding PDB structures, and only approximately 25% of the total dataset have >95% of their lengths observed in the corresponding PDB structures. This suggests that the vast majority of PDB proteins is shorter than their corresponding Swiss-Prot sequences and/or contain numerous residues, which are not observed in maps of electron density. To determine the prevalence of disordered regions in PDB, the residues in the Swiss-Prot sequences were grouped into four general categories, "Observed" (which correspond to structured regions), "Not observed" (regions with missing electron density, potentially disordered), "Uncharacterized," and "Ambiguous," depending on their appearance in the corresponding PDB entries. This non-redundant set of residues can be viewed as a 'fragment' or empirical domain database that contains a set of experimentally determined structured regions or domains and a set of experimentally verified disordered regions or domains. We studied the propensities and properties of residues in these four categories and analyzed their relations to the predictions of disorder using several algorithms. "Non-observed," "Ambiguous," and "Uncharacterized" regions were shown to possess the amino acid compositional biases typical of intrinsically disordered proteins. The application of four different disorder predictors (PONDR(R) VL-XT, VL3-BA, VSL1P, and IUPred) revealed that the vast majority of residues in the "Observed" dataset are ordered, and that the "Not observed" regions are mostly disordered. The "Uncharacterized" regions possess some tendency toward order, whereas the predictions for the short "Ambiguous" regions are really ambiguous. Long "Ambiguous" regions (>70 amino acid residues) are mostly predicted to be ordered, suggesting that they are likely to be "wobbly" domains. Overall, we showed that completely ordered proteins are not highly abundant in PDB and many PDB sequences have disordered regions. In fact, in the analyzed dataset approximately 10% of the PDB proteins contain regions of consecutive missing or ambiguous residues longer than 30 amino-acids and approximately 40% of the proteins possess short regions (> or =10 and < 30 amino-acid long) of missing and ambiguous residues.     

A substitution matrix for structural alphabet based on structural alignment of homologous proteins and its applications

Analysis of protein structures based on backbone structural patterns known as structural alphabets have been shown to be very useful. Among them, a set of 16 pentapeptide structural motifs known as protein blocks (PBs) has been identified and upon which backbone model of most protein structures can be built. PBs allows simplification of 3D space onto 1D space in the form of sequence of PBs. Here, for the first time, substitution probabilities of PBs in a large number of aligned homologous protein structures have been studied and are expressed as a simplified 16 x 16 substitution matrix. The matrix was validated by benchmarking how well it can align sequences of PBs rather like amino acid alignment to identify structurally equivalent regions in closely or distantly related proteins using dynamic programming approach. The alignment results obtained are very comparable to well established structure comparison methods like DALI and STAMP. Other interesting applications of the matrix have been investigated. We first show that, in variable regions between two superimposed homologous proteins, one can distinguish between local conformational differences and rigid-body displacement of a conserved motif by comparing the PBs and their substitution scores. Second, we demonstrate, with the example of aspartic proteinases, that PBs can be efficiently used to detect the lobe/domain flexibility in the multidomain proteins. Lastly, using protein kinase as an example, we identify regions of conformational variations and rigid body movements in the enzyme as it is changed to the active state from an inactive state.     

Intrinsic Disorder in the Protein Data Bank

Abstract

The Protein Data Bank (PDB) is the preeminent source of protein structural information. PDB contains over 32,500 experimentally determined 3-D structures solved using X-ray crystallography or nuclear magnetic resonance spectroscopy. Intrinsically disordered regions fail to form a fixed 3-D structure under physiological conditions. In this study, we compare the amino-acid sequences of proteins whose structures are determined by X-ray crystallography with the corresponding sequences from the Swiss-Prot database. The analyzed dataset includes 16,370 structures, which represent 18,101 PDB chains and 5,434 different proteins from 910 different organisms (2,793 eukaryotic, 2,109 bacterial, 288 viral, and 244 archaeal). In this dataset, on average, each Swiss-Prot protein is represented by 7 PDB chains with 76% of the crystallized regions being represented by more than one structure. Intriguingly, the complete sequences of only ~7% of proteins are observed in the corresponding PDB structures, and only ~25% of the total dataset have >95% of their lengths observed in the corresponding PDB structures. This suggests that the vast majority of PDB proteins is shorter than their corresponding Swiss-Prot sequences and/or contain numerous residues, which are not observed in maps of electron density. To determine the prevalence of disordered regions in PDB, the residues in the Swiss-Prot sequences were grouped into four general categories, “Observed” (which correspond to structured regions), “Not observed” (regions with missing electron density, potentially disordered), “Uncharacterized,” and “Ambiguous,” depending on their appearance in the corresponding PDB entries. This non-redundant set of residues can be viewed as a ‘fragment’ or empirical domain database that contains a set of experimentally determined structured regions or domains and a set of experimentally verified disordered regions or domains. We studied the propensities and properties of residues in these four categories and analyzed their relations to the predictions of disorder using several algorithms. “Non-observed,” “Ambiguous,” and “Uncharacterized” regions were shown to possess the amino acid compositional biases typical of intrinsically disordered proteins. The application of four different disorder predictors (PONDR® VL-XT, VL3-BA, VSL1P, and IUPred) revealed that the vast majority of residues in the “Observed” dataset are ordered, and that the “Not observed” regions are mostly disordered. The “Uncharacterized” regions possess some tendency toward order, whereas the predictions for the short “Ambiguous” regions are really ambiguous. Long “Ambiguous” regions (>70 amino acid residues) are mostly predicted to be ordered, suggesting that they are likely to be “wobbly” domains.

Overall, we showed that completely ordered proteins are not highly abundant in PDB and many PDB sequences have disordered regions. In fact, in the analyzed dataset ~10% of the PDB proteins contain regions of consecutive missing or ambiguous residues longer than 30 amino-acids and ~40% of the proteins possess short regions (≥10 and <30 amino-acid long) of missing and ambiguous residues.     

Homology model of surface antigen OmpC from Salmonella typhi and its functional implications

Homology based 3D structural model of the immunodominant major surface antigen OmpC from Salmonella typhi, an obligatory human pathogen, was built to understand the possible unique conformational features of its antigenic loops with respect to other immunologically cross reacting porins. The homology model was built based on the known crystal structures of the E. coli porins OmpF and PhoE. Structure based sequence alignment helped to define the structurally conserved regions (SCRs). The SCR regions of OmpC were modelled using the coordinates of corresponding regions from reference proteins. Surface exposed variable regions were modelled based on the sequence similarity and loop search in PDB. Structural refinement based on symmetry restrained energy minimization resulted in an agreeable model for the trimer of OmpC. The resulting model was compared with other porin structures, having b-barrel fold with 16 transmembrane beta-strands, and found that the variable regions are unique in terms of sequence and structure. A ranking of the loops taking into account the antigenic index, the sequence variability, the surface accessibility in the context of the trimer, and the structural variability suggests that loop 4 (151-172), loop 5 (194-218) and loop 6 (237-264) are the best ranked B-cell epitopes. The model provides possible explanations for the functional and unique immunological properties associated with the surface exposed regions and outlines the implications for structure based experimental design.     

Demonstration of the involvement of outer surface protein E coiled coil structural domains and higher order structural elements in the binding of infection-induced antibody and the complement-regulatory protein, factor H

Factor H (fH) is an important regulator of the alternative complement cascade. Several human pathogens have been shown to bind fH to their surface, a process that facilitates immune evasion or cell to cell interaction. Among the pathogens that bind fH are some Borrelia species associated with Lyme disease and relapsing fever. The fH-binding proteins of the Lyme spirochetes form two classes (I and II). In Borrelia burgdorferi B31MI, class I includes the outer surface protein E (OspE) paralogs, L39, N38, and P38, whereas the class II group includes A68 and additional proteins that have not yet been identified. To identify the OspE determinants involved in fH and OspE-targeting infection-induced Ab (iAb) binding, deletion, random, and site-directed mutagenesis of L39 were performed. Mutations in several different regions of L39 abolished fH and or iAb binding, indicating that separable domains and residues of OspE are required for ligand binding. Some of the mutants that lost the ability to bind fH, iAb, or both had only a single amino acid change. Site-directed mutagenesis of three putative coiled coil motifs of OspE revealed that these higher order structures are required for fH binding but not for iAb binding. The data presented within demonstrate that the binding of fH and iAb to the OspE protein is mediated by higher order structures and protein conformation. These studies advance our understanding of fH binding as a virulence mechanism and facilitate ongoing efforts to use fH-binding proteins in the development of microbial vaccines.     

A quantitative analysis of interfacial amino acid conservation in protein-protein hetero complexes

A long-standing question in molecular biology is whether interfaces of protein-protein complexes are more conserved than the rest of the protein surfaces. Although it has been reported that conservation can be used as an indicator for predicting interaction sites on proteins, there are recent reports stating that the interface regions are only slightly more conserved than the rest of the protein surfaces, with conservation signals not being statistically significant enough for predicting protein-protein binding sites. In order to properly address these controversial reports we have studied a set of 28 well resolved hetero complex structures of proteins that consists of transient and non-transient complexes. The surface positions were classified into four conservation classes and the conservation index of the surface positions was quantitatively analyzed. The results indicate that the surface density of highly conserved positions is significantly higher in the protein-protein interface regions compared with the other regions of the protein surface. However, the average conservation index of the patches in the interface region is not significantly higher compared with other surface regions of the protein structures. This finding demonstrates that the number of conserved residue positions is a more appropriate indicator for predicting protein-protein binding sites than the average conservation index in the interacting region. We have further validated our findings on a set of 59 benchmark complex structures. Furthermore, an analysis of 19 complexes of antigen-antibody interactions shows that there is no conservation of amino acid positions in the interacting regions of these complexes, as expected, with the variable region of the immunoglobulins interacting mostly with the antigens. Interestingly, antigen interacting regions also have a higher number of non-conserved residue positions in the interacting region than the rest of the protein surface.