Recognition of protein complexation based on hydrophobicity distribution

The identification of the surface area able to generate the protein-protein complexation ligand and ion ligation is critical for the recognition of the biological function of particular proteins. The technique based on the analysis of the irregularity of hydrophobicity distribution is used as the criterion for the recognition of the interaction regions. Particularly, the exposure of hydrophobic residues on the surface of protein as well as the localization of the hydrophilic residues in the hydrophobic core is treated as potential area ready to interact with external molecules. The model based on the “fuzzy oil drop” approach treating the protein molecule as the drop of hydrophobicity concentrated in the central part of structure with the hydrophobicity close to zero on the surface according to 3-dimensional Gauss function. The comparison with the observed hydrophobicy in particular protein reveals some irregularities. These irregularities seem to represent the aim-oriented localization.     

The Hydrophobicity of the H3 Histone Fold differs from the Hydrophobicity of the other three Folds

The eukaryotic histone dimers, H3–H4 and H2A–H2B, are formed in the cytosol prior to being transported into the nucleus and assembled into the nucleosome. Residue side-chain distances from the interior of the histone dimers are obtained with an ellipsoidal spatial metric and structural information provided by X-ray analyses at atomic resolution of the nucleosome core particles. While the spatial hydrophobic moment profiles of the dimers are comparable with profiles obtained previously that characterize the hydrophobic core of single-chain, single-domain globular soluble proteins, correlation coefficients between the side-chain hydrophobicities and distances from the interior of the H3–H4 dimer and H2A–H2B dimer differ significantly. This difference is traced to the H3 histone fold, which segregates fewer hydrophobic residues within the protein interior than the three other folds. Examination of the correlation coefficient between residue hydrophobicity and side-chain distance from the dimer interior over local regions of the fold sequence shows that the region of reduced correlation is associated mainly with the residues at the carboxyl end of the H3 histone fold, the helical region of the fold involved in the H3–H3 binding of the (H3–H4)₂ tetramer of the nucleosome. Hydrophobic interactions apparently contribute to the binding of this fourfold helical bundle and this evolutionary requirement may trade off against the requirement for H3–H4 dimer stability. The present results provide a different view than previously proposed, albeit of similar origin, to account for the reduced stability of the H3–H4 dimer compared with the H2A–H2B dimer.Reviewing Editor: Dr. Martin Kreitman     

Better prediction of protein contact number using a support vector regression analysis of amino acid sequence

Background  

Protein tertiary structure can be partly characterized via each amino acid's contact number measuring how residues are spatially arranged. The contact number of a residue in a folded protein is a measure of its exposure to the local environment, and is defined as the number of C _β atoms in other residues within a sphere around the C _β atom of the residue of interest. Contact number is partly conserved between protein folds and thus is useful for protein fold and structure prediction. In turn, each residue's contact number can be partially predicted from primary amino acid sequence, assisting tertiary fold analysis from sequence data. In this study, we provide a more accurate contact number prediction method from protein primary sequence.     

Type II restriction endonuclease R.Eco29kI is a member of the GIY-YIG nuclease superfamily

Background

The majority of experimentally determined crystal structures of Type II restriction endonucleases (REases) exhibit a common PD-(D/E)XK fold. Crystal structures have been also determined for single representatives of two other folds: PLD (R.BfiI) and half-pipe (R.PabI), and bioinformatics analyses supported by mutagenesis suggested that some REases belong to the HNH fold. Our previous bioinformatic analysis suggested that REase R.Eco29kI shares sequence similarities with one more unrelated nuclease superfamily, GIY-YIG, however so far no experimental data were available to support this prediction. The determination of a crystal structure of the GIY-YIG domain of homing endonuclease I-TevI provided a template for modeling of R.Eco29kI and prompted us to validate the model experimentally.

Results

Using protein fold-recognition methods we generated a new alignment between R.Eco29kI and I-TevI, which suggested a reassignment of one of the putative catalytic residues. A theoretical model of R.Eco29kI was constructed to illustrate its predicted three-dimensional fold and organization of the active site, comprising amino acid residues Y49, Y76, R104, H108, E142, and N154. A series of mutants was constructed to generate amino acid substitutions of selected residues (Y49A, R104A, H108F, E142A and N154L) and the mutant proteins were examined for their ability to bind the DNA containing the Eco29kI site 5'-CCGCGG-3' and to catalyze the cleavage reaction. Experimental data reveal that residues Y49, R104, E142, H108, and N154 are important for the nuclease activity of R.Eco29kI, while H108 and N154 are also important for specific DNA binding by this enzyme.

Conclusion

Substitutions of residues Y49, R104, H108, E142 and N154 predicted by the model to be a part of the active site lead to mutant proteins with strong defects in the REase activity. These results are in very good agreement with the structural model presented in this work and with our prediction that R.Eco29kI belongs to the GIY-YIG superfamily of nucleases. Our study provides the first experimental evidence for a Type IIP REase that does not belong to the PD-(D/E)XK or HNH superfamilies of nucleases, and is instead a member of the unrelated GIY-YIG superfamily.     

Structural and mechanistic investigations on Salmonella typhimurium acetate kinase (AckA): identification of a putative ligand binding pocket at the dimeric interface

Background

A large number of studies have been carried out to obtain amino acid propensities for ??-helices and ??-sheets. The obtained propensities for ??-helices are consistent with each other, and the pair-wise correlation coefficient is frequently high. On the other hand, the ??-sheet propensities obtained by several studies differed significantly, indicating that the context significantly affects ??-sheet propensity.

Results

We calculated amino acid propensities for ??-helices and ??-sheets for 39 and 24 protein folds, respectively, and addressed whether they correlate with the fold. The propensities were also calculated for exposed and buried sites, respectively. Results showed that ??-helix propensities do not differ significantly by fold, but ??-sheet propensities are diverse and depend on the fold. The propensities calculated for exposed sites and buried sites are similar for ??-helix, but such is not the case for the ??-sheet propensities. We also found some fold dependence on amino acid frequency in ??-strands. Folds with a high Ser, Thr and Asn content at exposed sites in ??-strands tend to have a low Leu, Ile, Glu, Lys and Arg content (correlation coefficient = ?0.90) and to have flat ??-sheets. At buried sites in ??-strands, the content of Tyr, Trp, Gln and Ser correlates negatively with the content of Val, Ile and Leu (correlation coefficient = ?0.93). "All-??" proteins tend to have a higher content of Tyr, Trp, Gln and Ser, whereas "??/??" proteins tend to have a higher content of Val, Ile and Leu.

Conclusions

The ??-helix propensities are similar for all folds and for exposed and buried residues. However, ??-sheet propensities calculated for exposed residues differ from those for buried residues, indicating that the exposed-residue fraction is one of the major factors governing amino acid composition in ??-strands. Furthermore, the correlations we detected suggest that amino acid composition is related to folding properties such as the twist of a ??-strand or association between two ?? sheets.     

Critical Structural and Functional Roles for the N-Terminal Insertion Sequence in Surfactant Protein B Analogs

Background

Surfactant protein B (SP-B; 79 residues) belongs to the saposin protein superfamily, and plays functional roles in lung surfactant. The disulfide cross-linked, N- and C-terminal domains of SP-B have been theoretically predicted to fold as charged, amphipathic helices, suggesting their participation in surfactant activities. Earlier structural studies with Mini-B, a disulfide-linked construct based on the N- and C-terminal regions of SP-B (i.e., ∼residues 8–25 and 63–78), confirmed that these neighboring domains are helical; moreover, Mini-B retains critical in vitro and in vivo surfactant functions of the native protein. Here, we perform similar analyses on a Super Mini-B construct that has native SP-B residues (1–7) attached to the N-terminus of Mini-B, to test whether the N-terminal sequence is also involved in surfactant activity.

Methodology/Results

FTIR spectra of Mini-B and Super Mini-B in either lipids or lipid-mimics indicated that these peptides share similar conformations, with primary α-helix and secondary β-sheet and loop-turns. Gel electrophoresis demonstrated that Super Mini-B was dimeric in SDS detergent-polyacrylamide, while Mini-B was monomeric. Surface plasmon resonance (SPR), predictive aggregation algorithms, and molecular dynamics (MD) and docking simulations further suggested a preliminary model for dimeric Super Mini-B, in which monomers self-associate to form a dimer peptide with a “saposin-like” fold. Similar to native SP-B, both Mini-B and Super Mini-B exhibit in vitro activity with spread films showing near-zero minimum surface tension during cycling using captive bubble surfactometry. In vivo, Super Mini-B demonstrates oxygenation and dynamic compliance that are greater than Mini-B and compare favorably to full-length SP-B.

Conclusion

Super Mini-B shows enhanced surfactant activity, probably due to the self-assembly of monomer peptide into dimer Super Mini-B that mimics the functions and putative structure of native SP-B.     

Influence of medium- and long-range interactions in different folding types of globular proteins

Recognition of protein fold from amino acid sequence is a challenging task. The structure and stability of proteins from different fold are mainly dictated by inter-residue interactions. In our earlier work, we have successfully used the medium- and long-range contacts for predicting the protein folding rates, discriminating globular and membrane proteins and for distinguishing protein structural classes. In this work, we analyze the role of inter-residue interactions in commonly occurring folds of globular proteins in order to understand their folding mechanisms. In the medium-range contacts, the globin fold and four-helical bundle proteins have more contacts than that of DNA-RNA fold although they all belong to all-alpha class. In long-range contacts, only the ribonuclease fold prefers 4-10 range and the other folding types prefer the range 21-30 in alpha/beta class proteins. Further, the preferred residues and residue pairs influenced by these different folds are discussed. The information about the preference of medium- and long-range contacts exhibited by the 20 amino acid residues can be effectively used to predict the folding type of each protein.     

Variation of amino acid properties in all-beta globular and outer membrane protein structures

Discriminating outer membrane (OM) proteins from globular proteins is an important task. The structural analysis of β-strands dominating globular (all-β) proteins and OM proteins provides useful insight to distinguish between them. In this work, we analyze the characteristic features of the 20 amino acid residues in all-β and OM proteins. We set up numerical indices for several properties of amino acid residues, such as, conformational parameters, surrounding hydrophobicity, accessible surface area and reduction in accessibility, and inter-residue contacts. We found that all the aromatic residues prefer to be in β-strands of both globular and OM proteins. The surrounding hydrophobicity of aromatic and non-polar amino acid residues in globular proteins is significantly higher than that of OM proteins. The residues Trp, Arg, Phe and Gln show a remarkable difference of reduction in accessibility between all-β globular (βG) and OM proteins. The positively charged residues, Lys and Arg in the membrane part of OM proteins have more number of contacts than globular proteins. Further, the behavior of the 20 amino acid residues in β-strand segments of globular and OM proteins have been discussed. The parameters developed in this work can be used for identifying transmembrane β-strands in OM proteins and for discriminating βG proteins from OM proteins.     

Structural and Mutational Analysis of Escherichia coli AlkB Provides Insight into Substrate Specificity and DNA Damage Searching

Background

In Escherichia coli, cytotoxic DNA methyl lesions on the N1 position of purines and N3 position of pyrimidines are primarily repaired by the 2-oxoglutarate (2-OG) iron(II) dependent dioxygenase, AlkB. AlkB repairs 1-methyladenine (1-meA) and 3-methylcytosine (3-meC) lesions, but it also repairs 1-methylguanine (1-meG) and 3-methylthymine (3-meT) at a much less efficient rate. How the AlkB enzyme is able to locate and identify methylated bases in ssDNA has remained an open question.

Methodology/Principal Findings

We determined the crystal structures of the E. coli AlkB protein holoenzyme and the AlkB-ssDNA complex containing a 1-meG lesion. We coupled this to site-directed mutagenesis of amino acids in and around the active site, and tested the effects of these mutations on the ability of the protein to bind both damaged and undamaged DNA, as well as catalyze repair of a methylated substrate.

Conclusions/Significance

A comparison of our substrate-bound AlkB-ssDNA complex with our unliganded holoenzyme reveals conformational changes of residues within the active site that are important for binding damaged bases. Site-directed mutagenesis of these residues reveals novel insight into their roles in DNA damage recognition and repair. Our data support a model that the AlkB protein utilizes at least two distinct conformations in searching and binding methylated bases within DNA: a “searching” mode and “repair” mode. Moreover, we are able to functionally separate these modes through mutagenesis of residues that affect one or the other binding state. Finally, our mutagenesis experiments show that amino acid D135 of AlkB participates in both substrate specificity and catalysis.     

Prediction of polypeptide fragments exposed to the solvent

A method is presented to predict those polypeptide segments within a globular protein that are more likely to be exposed to the solvent. The protein amino acidic sequence is the only information needed by this new algorithm. It uses a consensus hydrophobicity scale, derived from 28 known scales, and it is based on the comparison between the average hydrophobicity of a polypeptide fragment and the average hydrophobicity expected for a segment containing the same number of residues. The latter values are pre-computed from a non-redundant set of single chain protein structural domains. The comparison between the two average values results in a t value that readily provides the prediction with a statistical significance. A jack-knife validation analysis indicates that the protein segment predicted to be the most solvent exposed is actually solvent exposed and amongst the fragments that are most exposed. The source of a stand-alone program, written in C language, that allows the prediction of the most likely solvent exposed segment in a globular protein is available from the author.     

Increased sequence hydrophobicity reduces conformational specificity: A mutational case study of the Arc repressor protein

The amino-acid sequences of soluble, globular proteins must have hydrophobic residues to form a stable core, but excess sequence hydrophobicity can lead to loss of native state conformational specificity and aggregation. Previous studies of polar-to-hydrophobic mutations in the β-sheet of the Arc repressor dimer showed that a single substitution at position 11 (N11L) leads to population of an alternate dimeric fold in which the β-sheet is replaced by helix. Two additional hydrophobic mutations at positions 9 and 13 (Q9V and R13V) lead to population of a differently folded octamer along with both dimeric folds. Here we conduct a comprehensive study of the sequence determinants of this progressive loss of fold specificity. We find that the alternate dimer-fold specifically results from the N11L substitution and is not promoted by other hydrophobic substitutions in the β-sheet. We also find that three highly hydrophobic substitutions at positions 9, 11, and 13 are necessary and sufficient for oligomer formation, but the oligomer size depends on the identity of the hydrophobic residue in question. The hydrophobic substitutions increase thermal stability, illustrating how increased hydrophobicity can increase folding stability even as it degrades conformational specificity. The oligomeric variants are predicted to be aggregation-prone but may be hindered from doing so by proline residues that flank the β-sheet region. Loss of conformational specificity due to increased hydrophobicity can manifest itself at any level of structure, depending upon the specific mutations and the context in which they occur.     

The Refined Crystal Structure of an Eel Pout Type III Antifreeze Protein RD1 at 0.62-Å Resolution Reveals Structural Microheterogeneity of Protein and Solvation

RD1 is a 7-kDa globular protein from the Antarctic eel pout Lycodichthys dearborni. It belongs to type III of the four types of antifreeze proteins (AFPs) found in marine fishes living at subzero temperatures. For type III AFP, a potential ice-binding flat surface has been identified and is imbedded with side chains capable of making hydrogen bonds with a specific lattice plane on ice. So far, all crystallographic studies on type III AFPs were carried out using the Atlantic ocean pout Macrozoarces americanus as the source organism. Here we present the crystal structure of a type III AFP from a different zoarcid fish, and at an ultra-high resolution of 0.62 Å. The protein fold of RD1 comprises a compact globular domain with two internal tandem motifs arranged about a pseudo-dyad symmetry. Each motif of the “pretzel fold” includes four short β-strands and a 3₁₀ helix. There is a novel internal cavity of 45 ų surrounded by eight conserved nonpolar residues. The model contains several residues with alternate conformations, and a number of split water molecules, probably caused by alternate interactions with the protein molecule. After extensive refinement that includes hydrogen atoms, significant residual electron densities associated with the electrons of peptides and many other bonds could be visualized.     

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

Background

The members of cupin superfamily exhibit large variations in their sequences, functions, organization of domains, quaternary associations and the nature of bound metal ion, despite having a conserved β-barrel structural scaffold. Here, an attempt has been made to understand structure-function relationships among the members of this diverse superfamily and identify the principles governing functional diversity. The cupin superfamily also contains proteins for which the structures are available through world-wide structural genomics initiatives but characterized as “hypothetical”. We have explored the feasibility of obtaining clues to functions of such proteins by means of comparative analysis with cupins of known structure and function.

Methodology/Principal Findings

A 3-D structure-based phylogenetic approach was undertaken. Interestingly, a dendrogram generated solely on the basis of structural dissimilarity measure at the level of domain folds was found to cluster functionally similar members. This clustering also reflects an independent evolution of the two domains in bicupins. Close examination of structural superposition of members across various functional clusters reveals structural variations in regions that not only form the active site pocket but are also involved in interaction with another domain in the same polypeptide or in the oligomer.

Conclusions/Significance

Structure-based phylogeny of cupins can influence identification of functions of proteins of yet unknown function with cupin fold. This approach can be extended to other proteins with a common fold that show high evolutionary divergence. This approach is expected to have an influence on the function annotation in structural genomics initiatives.     

Nuclear magnetic resonance structure of the N-terminal domain of nonstructural protein 3 from the severe acute respiratory syndrome coronavirus

This paper describes the structure determination of nsp3a, the N-terminal domain of the severe acute respiratory syndrome coronavirus (SARS-CoV) nonstructural protein 3. nsp3a exhibits a ubiquitin-like globular fold of residues 1 to 112 and a flexibly extended glutamic acid-rich domain of residues 113 to 183. In addition to the four beta-strands and two alpha-helices that are common to ubiquitin-like folds, the globular domain of nsp3a contains two short helices representing a feature that has not previously been observed in these proteins. Nuclear magnetic resonance chemical shift perturbations showed that these unique structural elements are involved in interactions with single-stranded RNA. Structural similarities with proteins involved in various cell-signaling pathways indicate possible roles of nsp3a in viral infection and persistence.     

Structure of Membrane-active Toxin from Crab Spider Heriaeus melloteei Suggests Parallel Evolution of Sodium Channel Gating Modifiers in Araneomorphae and Mygalomorphae

We present a structural and functional study of a sodium channel activation inhibitor from crab spider venom. Hm-3 is an insecticidal peptide toxin consisting of 35 amino acid residues from the spider Heriaeus melloteei (Thomisidae). We produced Hm-3 recombinantly in Escherichia coli and determined its structure by NMR spectroscopy. Typical for spider toxins, Hm-3 was found to adopt the so-called “inhibitor cystine knot” or “knottin” fold stabilized by three disulfide bonds. Its molecule is amphiphilic with a hydrophobic ridge on the surface enriched in aromatic residues and surrounded by positive charges. Correspondingly, Hm-3 binds to both neutral and negatively charged lipid vesicles. Electrophysiological studies showed that at a concentration of 1 μm Hm-3 effectively inhibited a number of mammalian and insect sodium channels. Importantly, Hm-3 shifted the dependence of channel activation to more positive voltages. Moreover, the inhibition was voltage-dependent, and strong depolarizing prepulses attenuated Hm-3 activity. The toxin is therefore concluded to represent the first sodium channel gating modifier from an araneomorph spider and features a “membrane access” mechanism of action. Its amino acid sequence and position of the hydrophobic cluster are notably different from other known gating modifiers from spider venom, all of which are described from mygalomorph species. We hypothesize parallel evolution of inhibitor cystine knot toxins from Araneomorphae and Mygalomorphae suborders.     

Identification of mannose interacting residues using local composition

Background

Mannose binding proteins (MBPs) play a vital role in several biological functions such as defense mechanisms. These proteins bind to mannose on the surface of a wide range of pathogens and help in eliminating these pathogens from our body. Thus, it is important to identify mannose interacting residues (MIRs) in order to understand mechanism of recognition of pathogens by MBPs.

Results

This paper describes modules developed for predicting MIRs in a protein. Support vector machine (SVM) based models have been developed on 120 mannose binding protein chains, where no two chains have more than 25% sequence similarity. SVM models were developed on two types of datasets: 1) main dataset consists of 1029 mannose interacting and 1029 non-interacting residues, 2) realistic dataset consists of 1029 mannose interacting and 10320 non-interacting residues. In this study, firstly, we developed standard modules using binary and PSSM profile of patterns and got maximum MCC around 0.32. Secondly, we developed SVM modules using composition profile of patterns and achieved maximum MCC around 0.74 with accuracy 86.64% on main dataset. Thirdly, we developed a model on a realistic dataset and achieved maximum MCC of 0.62 with accuracy 93.08%. Based on this study, a standalone program and web server have been developed for predicting mannose interacting residues in proteins (http://www.imtech.res.in/raghava/premier/).

Conclusions

Compositional analysis of mannose interacting and non-interacting residues shows that certain types of residues are preferred in mannose interaction. It was also observed that residues around mannose interacting residues have a preference for certain types of residues. Composition of patterns/peptide/segment has been used for predicting MIRs and achieved reasonable high accuracy. It is possible that this novel strategy may be effective to predict other types of interacting residues. This study will be useful in annotating the function of protein as well as in understanding the role of mannose in the immune system.