Predicting reliable regions in protein alignments from sequence profiles

For applications such as comparative modelling one major issue is the reliability of sequence alignments. Reliable regions in alignments can be predicted using sub-optimal alignments of the same pair of sequences. Here we show that reliable regions in alignments can also be predicted from multiple sequence profile information alone.Alignments were created for a set of remotely related pairs of proteins using five different test methods. Structural alignments were used to assess the quality of the alignments and the aligned positions were scored using information from the observed frequencies of amino acid residues in sequence profiles pre-generated for each template structure. High-scoring regions of these profile-derived alignment scores were a good predictor of reliably aligned regions.These profile-derived alignment scores are easy to obtain and are applicable to any alignment method. They can be used to detect those regions of alignments that are reliably aligned and to help predict the quality of an alignment. For those residues within secondary structure elements, the regions predicted as reliably aligned agreed with the structural alignments for between 92% and 97.4% of the residues. In loop regions just under 92% of the residues predicted to be reliable agreed with the structural alignments. The percentage of residues predicted as reliable ranged from 32.1% for helix residues to 52.8% for strand residues.This information could also be used to help predict conserved binding sites from sequence alignments. Residues in the template that were identified as binding sites, that aligned to an identical amino acid residue and where the sequence alignment agreed with the structural alignment were in highly conserved, high scoring regions over 80% of the time. This suggests that many binding sites that are present in both target and template sequences are in sequence-conserved regions and that there is the possibility of translating reliability to binding site prediction.     

Free energetics of rigid body association of ubiquitin binding domains: A biochemical model for binding mediated by hydrophobic interaction

Weak intermolecular interactions, such as hydrophobic associations, underlie numerous biomolecular recognition processes. Ubiquitin is a small protein that represents a biochemical model for exploring thermodynamic signatures of hydrophobic association as it is widely held that a major component of ubiquitin's binding to numerous partners is mediated by hydrophobic regions on both partners. Here, we use atomistic molecular dynamics simulations in conjunction with the Adaptive Biasing Force sampling method to compute potentials of mean force (the reversible work, or free energy, associated with the binding process) to investigate the thermodynamic signature of complexation in this well‐studied biochemical model of hydrophobic association. We observe that much like in the case of a purely hydrophobic solute (i.e., graphene, carbon nanotubes), association is favored by entropic contributions from release of water from the interprotein regions. Moreover, association is disfavored by loss of enthalpic interactions, but unlike in the case of purely hydrophobic solutes, in this case protein‐water interactions are lost and not compensated for by additional water‐water interactions generated upon release of interprotein and moreso, hydration, water. We further find that relative orientations of the proteins that mutually present hydrophobic regions of each protein to its partner are favored over those that do not. In fact, the free energy minimum as predicted by a force field based method recapitulates the experimental NMR solution structure of the complex. Proteins 2014; 82:1453–1468. © 2014 Wiley Periodicals, Inc.     

A comparative study of available software for high-accuracy homology modeling: from sequence alignments to structural models

An open question in protein homology modeling is, how well do current modeling packages satisfy the dual criteria of quality of results and practical ease of use? To address this question objectively, we examined homology‐built models of a variety of therapeutically relevant proteins. The sequence identities across these proteins range from 19% to 76%. A novel metric, the difference alignment index (DAI), is developed to aid in quantifying the quality of local sequence alignments. The DAI is also used to construct the relative sequence alignment (RSA), a new representation of global sequence alignment that facilitates comparison of sequence alignments from different methods. Comparisons of the sequence alignments in terms of the RSA and alignment methodologies are made to better understand the advantages and caveats of each method. All sequence alignments and corresponding 3D models are compared to their respective structure‐based alignments and crystal structures. A variety of protein modeling software was used. We find that at sequence identities >40%, all packages give similar (and satisfactory) results; at lower sequence identities (<25%), the sequence alignments generated by Profit and Prime, which incorporate structural information in their sequence alignment, stand out from the rest. Moreover, the model generated by Prime in this low sequence identity region is noted to be superior to the rest. Additionally, we note that DSModeler and MOE, which generate reasonable models for sequence identities >25%, are significantly more functional and easier to use when compared with the other structure‐building software.     

德国小蠊泛素基因的克隆及序列分析

设计一对简并引物,从德国小蠊Blattellager manica细胞中克隆了泛素基因的编码区,GenBank登录号为AY501003。序列分析表明,该编码区的长度为228 bp,编码的多肽由76个氨基酸残基组成,相对分子质量为8.47 Kd ,其等电点为5.73。同源性比较发现,德国小蠊泛素基因与其他真核生物泛素基因在氨基酸水平上具有94%以上的相似性。     

Rat muscle acylphosphatase: Purification, amino sequence, and immunological characterization

Acylphosphatase was purified from rat skeletal muscle essentially by gel filtration and high-performance ion-exchange chromatography. The complete amino acid sequence was reconstructed by using the sequence data obtained from tryptic, peptic, andS. aureus V8 protease peptides. The protein consists of 96 amino acid residues and is acetylated at the NH₂-terminus. The immunological cross-reactivity of acylphosphatase from rat and horse skeletal muscle was examined by ELISA. The reaction with rabbit antiserum revealed the presence of at least five antigenic sites on rat enzyme, two of which are common to horse muscle enzyme. Anti-rat antibodies also recognize the peptide that corresponds to the initial part of the molecule, which varies greatly from equine enzyme. Two completely new antigenic sites are herein described: the first can be considered the main antigenic site and is located within positions 21–36, the second is in the COOH-terminal part of the molecule. A mixture of immunoreactive peptides gives strong antibody-antigen reaction inhibition (94%).     

Molecular dynamics of zinc-finger ubiquitin binding domains: a comparative study of histone deacetylase 6 and ubiquitin-specific protease 5

HDAC6 is a unique cytoplasmic histone deacetylase characterized by two deacetylase domains, and by a zinc-finger ubiquitin binding domain (ZnF-UBP) able to recognize ubiquitin (Ub). The latter has recently been demonstrated to be involved in the progression of neurodegenerative diseases and in mediating infection by the influenza A virus. Nowadays, understanding the dynamic and energetic features of HDAC6 ZnF-UBP-Ub recognition is considered as a crucial step for the conception of HDAC6 potential modulators. In this study, the atomic, solvent-related, and thermodynamic features behind HDAC6 ZnF-UBP-Ub recognition have been analyzed through molecular dynamics simulations. The behavior was then compared to the prototypical ZnF-UBP from ubiquitin-specific protease 5 (USP5) in order to spot relevant differences useful for selective drug design. Principal component analysis highlighted flapping motions of the L2A loop which were lowered down upon Ub binding in both systems. While polar and nonpolar interactions involving Ub G75 and G76 residues were also common features stabilizing both complexes, salt bridges showed a different pattern, more significant in HDAC6 ZnF-UBP-Ub, whose energetic contribution in USP5 ZnF-UBP-Ub was compensated by the presence of a more stable bridging water molecule. Whereas molecular mechanics/Poisson–Boltzmann surface area (MM-PBSA) free energies of binding were comparable for both systems, in agreement with experiments, computational alanine scanning and free energy decomposition data revealed that HDAC6 E1141 and D1178 are potential hotspots for the design of selective HDAC6 modulators.     

Structural and functional characterization of a ubiquitin variant engineered for tight and specific binding to an alpha‐helical ubiquitin interacting motif

Ubiquitin interacting motifs (UIMs) are short α‐helices found in a number of eukaryotic proteins. UIMs interact weakly but specifically with ubiquitin conjugated to other proteins, and in so doing, mediate specific cellular signals. Here we used phage display to generate ubiquitin variants (UbVs) targeting the N‐terminal UIM of the yeast Vps27 protein. Selections yielded UbV.v27.1, which recognized the cognate UIM with high specificity relative to other yeast UIMs and bound with an affinity more than two orders of magnitude higher than that of ubiquitin. Structural and mutational studies of the UbV.v27.1‐UIM complex revealed the molecular details for the enhanced affinity and specificity of UbV.v27.1, and underscored the importance of changes at the binding interface as well as at positions that do not contact the UIM. Our study highlights the power of the phage display approach for selecting UbVs with unprecedented affinity and high selectivity for particular α‐helical UIM domains within proteomes, and it establishes a general approach for the development of inhibitors targeting interactions of this type.     

Functional and receptor binding characterization of recombinant murine macrophage inflammatory protein 2: sequence analysis and mutagenesis identify receptor binding epitopes.

Murine macrophage inflammatory protein-2 (MIP-2), a member of the alpha-chemokine family, is one of several proteins secreted by cells in response to lipopolysaccharide. Many of the alpha-chemokines, such as interleukin-8, gro-alpha/MGSA, and neutrophil activating peptide-2 (NAP-2), are associated with neutrophil activation and chemotaxis. We describe the expression, purification, and characterization of murine MIP-2 from Pichia pastoris. Circular dichroism spectroscopy reveals that MIP-2 exhibits a highly ordered secondary structure consistent with the alpha/beta structures of other chemokines. Recombinant MIP-2 is chemotactic for human and murine neutrophils and up-regulates cell surface expression of Mac-1. MIP-2 binds to human and murine neutrophils with dissociation constants of 6.4 nM and 2.9 nM, respectively. We further characterize the binding of MIP-2 to the human types A and B IL-8 receptors and the murine homologue of the IL-8 receptor. MIP-2 displays low-affinity binding to the type A IL-8 receptor (Kd > 120 nM) and high-affinity binding to the type B IL-8 receptor (Kd 5.7 nM) and the murine receptor (Kd 6.8 nM). The three-dimensional structure of IL-8 and sequence analysis of six chemokines (IL-8, gro-alpha, NAP-2, ENA-78, KC, and MIP-2) that display high-affinity binding to the IL-8 type B receptor are used to identify an extended N-terminal surface that interacts with this receptor. Two mutants of MIP-2 establish that this region is also involved in binding and activating the murine homologue of the IL-8 receptor. Differences in the sequence between IL-8 and related chemokines identify a unique hydrophobic/aromatic region surrounded by charged residues that is likely to impart specificity to IL-8 for binding to the type A receptor.     

Differential ubiquitin binding of the UBA domains from human c-Cbl and Cbl-b: NMR structural and biochemical insights

The Cbl proteins, RING-type E3 ubiquitin ligases, are responsible for ubiquitinating the activated tyrosine kinases and targeting them for degradation. Both c-Cbl and Cbl-b have a UBA (ubiquitin-associated) domain at their C-terminal ends, and these two UBA domains share a high sequence similarity (75%). However, only the UBA from Cbl-b, but not from c-Cbl, can bind ubiquitin (Ub). To understand the mechanism by which the UBA domains specifically interact with Ub with different affinities, we determined the solution NMR structures of these two UBA domains, cUBA from human c-Cbl and UBAb from Cbl-b. Their structures show that these two UBA domains share the same fold, a compact three-helix bundle, highly resembling the typical UBA fold. Chemical shift perturbation experiments reveal that the helix-1 and loop-1 of UBAb form a predominately hydrophobic surface for Ub binding. By comparing the Ub-interacting surface on UBAb and its counterpart on cUBA, we find that the hydrophobic patch on cUBA is interrupted by a negatively charged residue Glu12. Fluorescence titration data show that the Ala12Glu mutant of UBAb completely loses the ability to bind Ub, whereas the mutation disrupting the dimerization has no significant effect on Ub binding. This study provides structural and biochemical insights into the Ub binding specificities of the Cbl UBA domains, in which the hydrophobic surface distribution on the first helix plays crucial roles in their differential affinities for Ub binding. That is, the amino acid residue diversity in the helix-1 region, but not the dimerization, determines the abilities of various UBA domains binding with Ub.     

NEMO specifically recognizes K63‐linked poly‐ubiquitin chains through a new bipartite ubiquitin‐binding domain

An important property of NEMO, the core element of the IKK complex involved in NF‐κB activation, resides in its ability to specifically recognize poly‐ubiquitin chains. A small domain called NOA/UBAN has been suggested to be responsible for this property. We recently demonstrated that the C‐terminal Zinc Finger (ZF) of NEMO is also able to bind ubiquitin. We show here by ZF swapping and mutagenesis that this represents its only function. While neither NOA nor ZF shows any preference for K63‐linked chains, we demonstrate that together they form a bipartite high‐affinity K63‐specific ubiquitin‐binding domain. A similar domain can be found in two other proteins, Optineurin and ABIN2, and can be freely exchanged with that of NEMO without interfering with its activity. This suggests that the main function of the C‐terminal half of NEMO is to specifically bind K63‐linked poly‐ubiquitin chains. We also demonstrate that the recently described binding of NEMO to linear poly‐ubiquitin chains is dependent on the NOA alone and does not require the presence of the ZF.     

Insights into the binding mode of sulphamates and sulphamides to hCA II: crystallographic studies and binding free energy calculations

Sulphamate and sulphamide derivatives have been largely investigated as carbonic anhydrase inhibitors (CAIs) by means of different experimental techniques. However, the structural determinants responsible for their different binding mode to the enzyme active site were not clearly defined so far. In this paper, we report the X-ray crystal structure of hCA II in complex with a sulphamate inhibitor incorporating a nitroimidazole moiety. The comparison with the structure of hCA II in complex with its sulphamide analogue revealed that the two inhibitors adopt a completely different binding mode within the hCA II active site. Starting from these results, we performed a theoretical study on sulphamate and sulphamide derivatives, demonstrating that electrostatic interactions with residues within the enzyme active site play a key role in determining their binding conformation. These findings open new perspectives in the design of effective CAIs using the sulphamate and sulphamide zinc binding groups as lead compounds.     

RosettaDDGPrediction for high-throughput mutational scans: From stability to binding

Reliable prediction of free energy changes upon amino acid substitutions (ΔΔGs) is crucial to investigate their impact on protein stability and protein–protein interaction. Advances in experimental mutational scans allow high-throughput studies thanks to multiplex techniques. On the other hand, genomics initiatives provide a large amount of data on disease-related variants that can benefit from analyses with structure-based methods. Therefore, the computational field should keep the same pace and provide new tools for fast and accurate high-throughput ΔΔG calculations. In this context, the Rosetta modeling suite implements effective approaches to predict folding/unfolding ΔΔGs in a protein monomer upon amino acid substitutions and calculate the changes in binding free energy in protein complexes. However, their application can be challenging to users without extensive experience with Rosetta. Furthermore, Rosetta protocols for ΔΔG prediction are designed considering one variant at a time, making the setup of high-throughput screenings cumbersome. For these reasons, we devised RosettaDDGPrediction, a customizable Python wrapper designed to run free energy calculations on a set of amino acid substitutions using Rosetta protocols with little intervention from the user. Moreover, RosettaDDGPrediction assists with checking completed runs and aggregates raw data for multiple variants, as well as generates publication-ready graphics. We showed the potential of the tool in four case studies, including variants of uncertain significance in childhood cancer, proteins with known experimental unfolding ΔΔGs values, interactions between target proteins and disordered motifs, and phosphomimetics. RosettaDDGPrediction is available, free of charge and under GNU General Public License v3.0, at https://github.com/ELELAB/RosettaDDGPrediction .     

CHDL: a cadherin-like domain in Proteobacteria and Cyanobacteria

We identified a cadherin-like domain (CHDL) using computational analysis. The CHDL domain is mostly distributed in Proteobacteria and Cyanobacteria, although it is also found in some eukaryotic proteins. Prediction of three-dimensional protein folding indicated that the CHDL domain has an immunoglobulin beta-sandwich fold and belongs to the cadherin superfamily. The CHDL domain does not have LDRE and DxNDN motifs, which are conserved in the cadherin domain, but has three other motifs: PxAxxD, DxDxD and YT-V/I-S/T-D, which might contribute to forming a calcium-binding site. The identification of this cadherin-like domain indicates that the cadherin superfamily may exhibit wider sequence and structural diversity than previously appreciated. Domain architecture analysis revealed that the CHDL domain is also associated with other adhesion domains as well as enzyme domains. Based on computational analysis and previous experimental data, we predict that the CHDL domain has calcium-binding and also carbohydrate-binding activity.     

Competitive binding of UBPY and ubiquitin to the STAM2 SH3 domain revealed by NMR

To date, the signal transducing adaptor molecule 2 (STAM2) was shown to harbour two ubiquitin binding domains (UBDs) known as the VHS and UIM domains, while the SH3 domain of STAM2 was reported to interact with deubiquitinating enzymes (DUBs) like UBPY and AMSH. In the present study, NMR evidences the interaction of the STAM2 SH3 domain with ubiquitin, demonstrating that SH3 constitutes the third UBD of STAM2. Furthermore, we show that a UBPY-derived peptide can outcompete ubiquitin for SH3 binding and vice versa. These results suggest that the SH3 domain of STAM2 plays versatile roles in the context of ubiquitin mediated receptor sorting.     

棉卷叶野螟泛素基因的克隆、序列分析及原核表达

本研究用RT-PCR方法,克隆了棉卷叶野螟Haritalodes derogata (Fabricius)泛素基因编码区,GenBank登录号为EU580145。序列分析表明,该编码区长228 bp,编码76个氨基酸,推测的编码蛋白的相对分子质量和等电点分别为8.53 kD和5.83。同源性比较发现,棉卷叶野螟泛素基因与其他10种昆虫泛素基因在氨基酸水平上具有93%以上的相似性。系统发育树显示棉卷叶野螟与斜纹夜蛾Spodoptera litura (Fabricius)遗传距离较近,通过同源建模获得了该棉卷叶野螟基因编码蛋白的理论三维结构。将棉卷叶野螟泛素基因与pET-32a(+)连接,构建原核表达载体pET-32a-ub,经IPTG诱导,棉卷叶野螟泛素基因在大肠杆菌BL21(DE3) 中高效表达。本研究成功克隆了棉卷叶野螟泛素基因的编码区,并经Western blotting分析证明实现了该基因的原核表达,为进一步研究其在该昆虫体内的作用机理奠定了基础。     

Solution structure of the N‐terminal domain of DC‐UbP/UBTD2 and its interaction with ubiquitin

DC‐UbP/UBTD2 is a ubiquitin (Ub) domain‐containing protein first identified from dendritic cells, and is implicated in ubiquitination pathway. The solution structure and backbone dynamics of the C‐terminal Ub‐like (UbL) domain were elucidated in our previous work. To further understand the biological function of DC‐UbP, we then solved the solution structure of the N‐terminal domain of DC‐UbP (DC‐UbP_N) and studied its Ub binding properties by NMR techniques. The results show that DC‐UbP_N holds a novel structural fold and acts as a Ub‐binding domain (UBD) but with low affinity. This implies that the DC‐UbP protein, composing of a combination of both UbL and UBD domains, might play an important role in regulating protein ubiquitination and delivery of ubiquitinated substrates in eukaryotic cells.     

Identification and evaluation of universal epitopes in Plasmodium vivax Duffy binding protein

Selected PvDBP-derived synthetic peptides were tested in competition assays with HLA molecules in order to identify and evaluate their binding to a wide range of MHC class II molecules. Binding was evaluated as the peptide’s ability to displace the biotinylated control peptide (HA_306-318) and was detected by a conventional ELISA. Thus, one epitope for the HLA-DR1 molecule, two epitopes for the HLA-DR4 molecule, six epitopes for the HLA-DR7 molecule and three epitopes for the HLA-DR11 molecule displaying a high binding percentage (above 50%) were experimentally obtained. The in vitro results were compared with the epitope prediction results. Two peptides behaved as universal epitopes since they bound to a larger number of HLA-DR molecules. Given that these peptides are located in the conserved PvDBP region II, they could be considered good candidates to be included in the design of a synthetic vaccine against Plasmodium vivax malaria.     

fingerprint: visual depiction of variation in multiple sequence alignments

There is a lack of programs available that focus on providing an overview of an aligned set of sequences such that the comparison of homologous sites becomes comprehensible and intuitive. Being able to identify similarities, differences, and patterns within a multiple sequence alignment is biologically valuable because it permits visualization of the distribution of a particular feature and inferences about the structure, function, and evolution of the sequences in question. We have therefore created a web server, fingerprint, which combines the characteristics of existing programs that represent identity, variability, charge, hydrophobicity, solvent accessibility, and structure along with new visualizations based on composition, heterogeneity, heterozygosity, d_N/d_S and nucleotide diversity. fingerprint is easy to use and globally accessible through any computer using any major browser. fingerprint is available at http://evol.mcmaster.ca/fingerprint/ .     

Empirical analysis of protein insertions and deletions determining parameters for the correct placement of gaps in protein sequence alignments

To understand how protein segments are inserted and deleted during divergent evolution, a set of pairwise alignments contained exactly one gap, and therefore arising from the first insertion-deletion (indel) event in the time separating the homologs, was examined. The alignments showed that "structure breaking" amino acids (PGDNS) were preferred within and flanking gapped regions, as are two residues with hydrophilic side-chains (QE) that frequently occur at the surface of protein folds. Conversely, hydrophobic residues (FMILYVW) occur infrequently within and flanking the gapped region. These preferences are modestly different in protein pairs separated by an episode of adaptive evolution, than in pairs diverging under strong functional constraints. Surprisingly, regions near an indel have not evolved more rapidly than the sequence pair overall, showing no evidence that an indel event must be compensated by local amino acid replacement. The gap-lengths are best approximated by a Zipfian distribution, with the probability of a gap of length L decreasing as a function of L(-1.8). These features are largely independent of the length of the gap and the extent of divergence (measured by both silent and non-silent sequence changes) separating the two proteins. Surprisingly, amino acid repeats were discovered in more than a third of the polypeptide segments in and around the gap. These correspond to repeats in the DNA sequence. This suggests that a signature of the mechanism by which indels occur in the DNA sequence remains in the encoded protein sequences. These data suggest specific tools to score gap placement in an alignment. They also suggest tools that distinguish true indels from gaps created by mistaken gene finding, including under-predicted and over-predicted introns. By providing mechanisms to identify errors, the tools will enhance the value of genome sequence databases in support of integrated paleogenomics strategies used to extract functional information in a post-genomic environment.     

RosettaDDGPrediction for high‐throughput mutational scans: From stability to binding

Reliable prediction of free energy changes upon amino acid substitutions (ΔΔGs) is crucial to investigate their impact on protein stability and protein–protein interaction. Advances in experimental mutational scans allow high‐throughput studies thanks to multiplex techniques. On the other hand, genomics initiatives provide a large amount of data on disease‐related variants that can benefit from analyses with structure‐based methods. Therefore, the computational field should keep the same pace and provide new tools for fast and accurate high‐throughput ΔΔG calculations. In this context, the Rosetta modeling suite implements effective approaches to predict folding/unfolding ΔΔGs in a protein monomer upon amino acid substitutions and calculate the changes in binding free energy in protein complexes. However, their application can be challenging to users without extensive experience with Rosetta. Furthermore, Rosetta protocols for ΔΔG prediction are designed considering one variant at a time, making the setup of high‐throughput screenings cumbersome. For these reasons, we devised RosettaDDGPrediction, a customizable Python wrapper designed to run free energy calculations on a set of amino acid substitutions using Rosetta protocols with little intervention from the user. Moreover, RosettaDDGPrediction assists with checking completed runs and aggregates raw data for multiple variants, as well as generates publication‐ready graphics. We showed the potential of the tool in four case studies, including variants of uncertain significance in childhood cancer, proteins with known experimental unfolding ΔΔGs values, interactions between target proteins and disordered motifs, and phosphomimetics. RosettaDDGPrediction is available, free of charge and under GNU General Public License v3.0, at https://github.com/ELELAB/RosettaDDGPrediction.