Localization of ligand binding site in proteins identified in silico

Knowledge-based models for protein folding assume that the early-stage structural form of a polypeptide is determined by the backbone conformation, followed by hydrophobic collapse. Side chain–side chain interactions, mostly of hydrophobic character, lead to the formation of the hydrophobic core, which seems to stabilize the structure of the protein in its natural environment. The fuzzy-oil-drop model is employed to represent the idealized hydrophobicity distribution in the protein molecule. Comparing it with the one empirically observed in the protein molecule reveals that they are not in agreement. It is shown in this study that the irregularity of hydrophobic distributions is aim-oriented. The character and strength of these irregularities in the organization of the hydrophobic core point to the specificity of a particular protein’s structure/function. When the location of these irregularities is determined versus the idealized fuzzy-oil-drop, function-related areas in the protein molecule can be identified. The presented model can also be used to identify ways in which protein–protein complexes can possibly be created. Active sites can be predicted for any protein structure according to the presented model with the free prediction server at . The implication based on the model presented in this work suggests the necessity of active presence of ligand during the protein folding process simulation. Figure Fuzzy-oil-drop model applied to identify the ligation site in lysozyme complexed with N-acetylglucosamine (PDB ID:1LMQ) in form of hydrophobicity deficiency (ΔH) profile and three-dimensional distribution of on protein surface     

Hydrophobic peptide tags as tools in bioseparation

Hydrophobic interactions are highly selective, and differences in surface hydrophobicities between proteins can be used as an efficient handle to facilitate protein isolation. Aromatic amino acid residues are of particular importance for molecular recognition because they have a key role in several biological functions. The hydrophobicity of a protein can easily be altered with minor genetic modifications, such as site-directed mutagenesis or fusions of hydrophobic peptide tags. An important advantage of hydrophobic peptide tags over traditional affinity tags is the possibility of exploring simple and inexpensive bioseparation materials. Recent results demonstrate the potential of hydrophobic interaction chromatography and aqueous two-phase systems as tools to study relative hydrophobicities of recombinant proteins with only minor alterations. This review focuses on hydrophobic peptide tags as fusion partners, which can be used as important tools in bioseparation.     

Contribution to the Prediction of the Fold Code: Application to Immunoglobulin and Flavodoxin Cases

Background

Folding nucleus of globular proteins formation starts by the mutual interaction of a group of hydrophobic amino acids whose close contacts allow subsequent formation and stability of the 3D structure. These early steps can be predicted by simulation of the folding process through a Monte Carlo (MC) coarse grain model in a discrete space. We previously defined MIRs (Most Interacting Residues), as the set of residues presenting a large number of non-covalent neighbour interactions during such simulation. MIRs are good candidates to define the minimal number of residues giving rise to a given fold instead of another one, although their proportion is rather high, typically [15-20]% of the sequences. Having in mind experiments with two sequences of very high levels of sequence identity (up to 90%) but different folds, we combined the MIR method, which takes sequence as single input, with the “fuzzy oil drop” (FOD) model that requires a 3D structure, in order to estimate the residues coding for the fold. FOD assumes that a globular protein follows an idealised 3D Gaussian distribution of hydrophobicity density, with the maximum in the centre and minima at the surface of the “drop”. If the actual local density of hydrophobicity around a given amino acid is as high as the ideal one, then this amino acid is assigned to the core of the globular protein, and it is assumed to follow the FOD model. Therefore one obtains a distribution of the amino acids of a protein according to their agreement or rejection with the FOD model.

Results

We compared and combined MIR and FOD methods to define the minimal nucleus, or keystone, of two populated folds: immunoglobulin-like (Ig) and flavodoxins (Flav). The combination of these two approaches defines some positions both predicted as a MIR and assigned as accordant with the FOD model. It is shown here that for these two folds, the intersection of the predicted sets of residues significantly differs from random selection. It reduces the number of selected residues by each individual method and allows a reasonable agreement with experimentally determined key residues coding for the particular fold. In addition, the intersection of the two methods significantly increases the specificity of the prediction, providing a robust set of residues that constitute the folding nucleus.     

Aromatic cluster mutations produce focal modulations of β-sheet structure

Site-directed mutagenesis is a powerful tool for altering the structure and function of proteins in a focused manner. Here, we examined how a model β-sheet protein could be tuned by mutation of numerous surface-exposed residues to aromatic amino acids. We designed these aromatic side chain “clusters” at highly solvent-exposed positions in the flat, single-layer β-sheet of Borrelia outer surface protein A (OspA). This unusual β-sheet scaffold allows us to interrogate the effects of these mutations in the context of well-defined structure but in the absence of the strong scaffolding effects of globular protein architecture. We anticipated that the introduction of a cluster of aromatic amino acid residues on the β-sheet surface would result in large conformational changes and/or stabilization and thereby provide new means of controlling the properties of β-sheets. Surprisingly, X-ray crystal structures revealed that the introduction of aromatic clusters produced only subtle conformational changes in the OspA β-sheet. Additionally, despite burying a large degree of hydrophobic surface area, the aromatic cluster mutants were slightly less stable than the wild-type scaffold. These results thereby demonstrate that the introduction of aromatic cluster mutations can serve as a means for subtly modulating β-sheet conformation in protein design.     

In vitro reconstitution of substrate S-acylation by the zDHHC family of protein acyltransferases

Protein S-acylation, more commonly known as protein palmitoylation, is a biological process defined by the covalent attachment of long chain fatty acids onto cysteine residues of a protein, effectively altering the local hydrophobicity and influencing its stability, localization and overall function. Observed ubiquitously in all eukaryotes, this post translational modification is mediated by the 23-member family of zDHHC protein acyltransferases in mammals. There are thousands of proteins that are S-acylated and multiple zDHHC enzymes can potentially act on a single substrate. Since its discovery, numerous methods have been developed for the identification of zDHHC substrates and the individual members of the family that catalyse their acylation. Despite these recent advances in assay development, there is a persistent gap in knowledge relating to zDHHC substrate specificity and recognition, that can only be thoroughly addressed through in vitro reconstitution. Herein, we will review the various methods currently available for reconstitution of protein S-acylation for the purposes of identifying enzyme–substrate pairs with a particular emphasis on the advantages and disadvantages of each approach.     

Local hydrophobicity stabilizes secondary structures in proteins

The probability of occurrence of helix and β-sheet residues in 47 globular proteins was determined as a function of local hydrophobicity, which was defined by the sum of the Nozaki-Tanford transfer free energies at two nearest-neighbors on both sides of the amino acid sequence. In general, hydrophilic amino acids favor neither helix nor β-sheet formations when neighbor residues are also hydrophilic but favor helix formation at higher local hydrophobicity. On the other hand, some hydrophobic amino acids such as Met, Leu, and Ile favor helix formation when neighbor residues are hydrophilic. None of the hydrophobic amino acids favor β-sheet formation with hydrophilic neighbors, but most of them strongly favor β-sheet formation at high local hydrophobicity. When the average of 20 amino acids is taken, both helix and β-sheet residue probabilities are higher at higher local hydrophobicity, although the increase is steeper for β-sheets. Therefore, β-sheet formation is more influenced by local hydrophobicity than helix formation. Generally, helices are nearer the surface and tend to have hydrophilic and hydrophobic faces at opposite sides. The tendency of alternating regions of hydrophilic and hydrophobic residues in a helical sequence was revealed by calculating the correlation of the Nozaki-Tanford values. Such amphipathic helices may be important in protein–protein and protein–lipid interactions and in forming hydrophilic channels in the membrane. The choice of 30 nonhomologous proteins as the data set did not alter the above results.     

Protein‐spanning water networks and implications for prediction of protein–protein interactions mediated through hydrophobic effects

Hydrophobic effects, often conflated with hydrophobic forces, are implicated as major determinants in biological association and self‐assembly processes. Protein–protein interactions involved in signaling pathways in living systems are a prime example where hydrophobic effects have profound implications. In the context of protein–protein interactions, a priori knowledge of relevant binding interfaces (i.e., clusters of residues involved directly with binding interactions) is difficult. In the case of hydrophobically mediated interactions, use of hydropathy‐based methods relying on single residue hydrophobicity properties are routinely and widely used to predict propensities for such residues to be present in hydrophobic interfaces. However, recent studies suggest that consideration of hydrophobicity for single residues on a protein surface require accounting of the local environment dictated by neighboring residues and local water. In this study, we use a method derived from percolation theory to evaluate spanning water networks in the first hydration shells of a series of small proteins. We use residue‐based water density and single‐linkage clustering methods to predict hydrophobic regions of proteins; these regions are putatively involved in binding interactions. We find that this simple method is able to predict with sufficient accuracy and coverage the binding interface residues of a series of proteins. The approach is competitive with automated servers. The results of this study highlight the importance of accounting of local environment in determining the hydrophobic nature of individual residues on protein surfaces. Proteins 2014; 82:3312–3326. © 2014 Wiley Periodicals, Inc.     

The Structure of the Cytomegalovirus-Encoded m04 Glycoprotein,a Prototypical Member of the m02 Family of Immunoevasins

The ability of CMVs to evade the immune system of the host is dependent on the expression of a wide array of glycoproteins, many of which interfere with natural killer cell function. In murine CMV, two large protein families mediate this immune-evasive function. Although it is established that the m145 family members mimic the structure of MHC-I molecules, the structure of the m02 family remains unknown. The most extensively studied m02 family member is m04, a glycoprotein that escorts newly assembled MHC-I molecules to the cell surface, presumably to avoid “missing self” recognition. Here we report the crystal structure of the m04 ectodomain, thereby providing insight into this large immunoevasin family. m04 adopted a β-sandwich immunoglobulin variable (Ig-V)-like fold, despite sharing very little sequence identity with the Ig-V superfamily. In addition to the Ig-V core, m04 possesses several unique structural features that included an unusual β-strand topology, a number of extended loops and a prominent α-helix. The m04 interior was packed by a myriad of hydrophobic residues that form distinct clusters around two conserved tryptophan residues. This hydrophobic core was well conserved throughout the m02 family, thereby indicating that murine CMV encodes a number of Ig-V-like molecules. We show that m04 binds a range of MHC-I molecules with low affinity in a peptide-independent manner. Accordingly, the structure of m04, which represents the first example of an murine CMV encoded Ig-V fold, provides a basis for understanding the structure and function of this enigmatic and large family of immunoevasins.     

NMR structure of the N-terminal coiled coil domain of the Andes hantavirus nucleocapsid protein

The hantaviruses are emerging infectious viruses that in humans can cause a cardiopulmonary syndrome or a hemorrhagic fever with renal syndrome. The nucleocapsid (N) is the most abundant viral protein, and during viral assembly, the N protein forms trimers and packages the viral RNA genome. Here, we report the NMR structure of the N-terminal domain (residues 1-74, called N1-74) of the Andes hantavirus N protein. N1-74 forms two long helices (alpha1 and alpha2) that intertwine into a coiled coil domain. The conserved hydrophobic residues at the helix alpha1-alpha2 interface stabilize the coiled coil; however, there are many conserved surface residues whose function is not known. Site-directed mutagenesis, CD spectroscopy, and immunocytochemistry reveal that a point mutation in the conserved basic surface formed by Arg22 or Lys26 lead to antibody recognition based on the subcellular localization of the N protein. Thus, Arg22 and Lys26 are likely involved in a conformational change or molecular recognition when the N protein is trafficked from the cytoplasm to the Golgi, the site of viral assembly and maturation.     

Epitope Tags beside the N-Terminal Cytoplasmic Tail of Human BST-2 Alter Its Intracellular Trafficking and HIV-1 Restriction

BST-2 blocks the particle release of various enveloped viruses including HIV-1, and this antiviral activity is dependent on the topological arrangement of its four structural domains. Several functions of the cytoplasmic tail (CT) of BST-2 have been previously discussed, but the exact role of this domain remains to be clearly defined. In this study, we investigated the impact of truncation and commonly-used tags addition into the CT region of human BST-2 on its intracellular trafficking and signaling as well as its anti-HIV-1 function. The CT-truncated BST-2 exhibited potent inhibition on Vpu-defective HIV-1 and even wild-type HIV-1. However, the N-terminal HA-tagged CT-truncated BST-2 retained little antiviral activity and dramatically differed from its original protein in the cell surface level and intracellular localization. Further, we showed that the replacement of the CT domain with a hydrophobic tag altered BST-2 function possibly by preventing its normal vesicular trafficking. Notably, we demonstrated that a positive charged motif “KRXK” in the conjunctive region between the cytotail and the transmembrane domain which is conserved in primate BST-2 is important for the protein trafficking and the antiviral function. These results suggest that although the CT of BST-2 is not essential for its antiviral activity, the composition of residues in this region may play important roles in its normal trafficking which subsequently affected its function. These observations provide additional implications for the structure-function model of BST-2.     

Replica Exchange Molecular Dynamics Simulations Provide Insight into Substrate Recognition by Small Heat Shock Proteins

The small heat shock proteins (sHSPs) are a virtually ubiquitous and diverse group of molecular chaperones that can bind and protect unfolding proteins from irreversible aggregation. It has been suggested that intrinsic disorder of the N-terminal arm (NTA) of sHSPs is important for substrate recognition. To investigate conformations of the NTA that could recognize substrates we performed replica exchange molecular dynamics simulations. Behavior at normal and stress temperatures of the dimeric building blocks of dodecameric HSPs from wheat (Ta16.9) and pea (Ps18.1) were compared because they display high sequence similarity, but Ps18.1 is more efficient in binding specific substrates. In our simulations, the NTAs of the dimer are flexible and dynamic; however, rather than exhibiting highly extended conformations they retain considerable α-helical character and contacts with the conserved α-crystallin domain (ACD). Network analysis and clustering methods reveal that there are two major conformational forms designated either “open” or “closed” based on the relative position of the two NTAs and their hydrophobic solvent accessible surface area. The equilibrium constant for the closed to open transition is significantly different for Ta16.9 and Ps18.1, with the latter showing more open conformations at elevated temperature correlated with its more effective chaperone activity. In addition, the Ps18.1 NTAs have more hydrophobic solvent accessible surface than those of Ta16.9. NTA hydrophobic patches are comparable in size to the area buried in many protein-protein interactions, which would enable sHSPs to bind early unfolding intermediates. Reduced interactions of the Ps18.1 NTAs with each other and with the ACD contribute to the differences in dynamics and hydrophobic surface area of the two sHSPs. These data support a major role for the conformational equilibrium of the NTA in substrate binding and indicate features of the NTA that contribute to sHSP chaperone efficiency.     

New approaches for predicting protein retention time in hydrophobic interaction chromatography

Hydrophobic interaction chromatography (HIC) is an important technique for the purification of proteins. In this paper, we review three different approaches for predicting protein retention time in HIC, based either on a protein's structure or on its amino-acidic composition, and we have extended one of these approaches. The first approach correlates the protein retention time in HIC with the protein average surface hydrophobicity. This methodology is based on the protein three-dimensional structure data and considers the hydrophobic contribution of the exposed amino acid residues as a weighted average. The second approach, which we have extended, is based on the high correlation level between the average surface hydrophobicity of a protein's hydrophobic interacting zone and its retention time in HIC. Finally, a third approach carries out a prediction of the average surface hydrophobicity of a protein, using only its amino-acidic composition, without knowing its three-dimensional structure. These models would make it possible to test different operating conditions for the purification of a target protein by computer simulations, and thus make it easier to select the optimal conditions, contributing to the rational design and optimization of the process.     

Molecular dissection,tissue localization and Ca2+ binding of the ryanodine receptor of Caenorhabditis elegans

We investigated the possible role of residues at the Ccap position in an alpha-helix on protein stability. A set of 431 protein alpha-helices containing a C'-Gly from the Protein Data Bank (PDB) was analyzed, and the normalized frequencies for finding particular residues at the Ccap position, the average fraction of buried surface area, and the hydrogen bonding patterns of the Ccap residue side-chain were calculated. We found that on average the Ccap position is 70% buried and noted a significant correlation (R=0.8) between the relative burial of this residue and its hydrophobicity as defined by the Gibbs energy of transfer from octanol or cyclohexane to water. Ccap residues with polar side-chains are commonly involved in hydrogen bonding. The hydrogen bonding pattern is such that, the longer side-chains of Glu, Gln, Arg, Lys, His form hydrogen bonds with residues distal (>+/-4) in sequence, while the shorter side-chains of Asp, Asn, Ser, Thr exhibit hydrogen bonds with residues close in sequence (<+/-4), mainly involving backbone atoms. Experimentally we determined the thermodynamic propensities of residues at the Ccap position using the protein ubiquitin as a model system. We observed a large variation in the stability of the ubiquitin variants depending on the nature of the Ccap residue. Furthermore, the measured changes in stability of the ubiquitin variants correlate with the hydrophobicity of the Ccap residue. The experimental results, together with the statistical analysis of protein structures from the PDB, indicate that the key hydrophobic capping interactions between a helical residue (C3 or C4) and a residue outside the helix (C", C3' or C4') are frequently enhanced by the hydrophobic interactions with Ccap residues.     

Gauss-function-Based model of hydrophobicity density in proteins

The model adopting the three-dimensional Gauss function to express the hydrophobicity distribution in proteins is presented in this paper. The tendency to create the hydrophobic center during protein folding is expressed in form of an external force field of the form of three-dimensional Gauss function which directs the folding polypeptide to locate the hydrophobic residues in a central part of the molecule and hydrophilic ones exposed toward the molecular surface. The decrease of the differences between hydrophobicity distribution as it appears at each step of the folding simulation and the expected hydrophobicity distribution (three-dimensional Gauss function) is the convergence criterion together with traditional non-bonding interaction optimization. The model was applied to fold the hypothetical membrane protein (target protein in CASP6) TA0354_69_121 from Thermoplasma acidophilum.     

The hydrophobic moment and its use in the classification of amphiphilic structures (review)

Amphiphilic alpha-helices play a major role in membrane dependent processes and are manifested in the primary structure of a protein by the periodic appearance of hydrophobic residues. Based on these periodic sequences, the hydrophobic moment was introduced, , which essentially treats the hydrophobicity of amino acid residues as a two-dimensional vector sum and provides a measure of amphiphilicity within regular repeat structures. To identify putative amphiphilic alpha-helix forming sequences, hydrophobic moment analysis assumes an amino acid residue periodicity of 100 and scans protein primary structures to find the 11-residue window with maximal . Taken with the window's mean hydrophobicity, , hydrophobic moment plot analysis uses the coordinate pair, [, ] to classify alpha-helices as either surface active, globular or transmembrane. More recently, this latter analysis has been extended to recognize candidate oblique orientated alpha-helices. Here, the hydrophobic moment is reviewed and data to query the logic of using a fixed window length and a fixed residue angular periodicity in hydrophobic moment analysis are provided. In addition, problems associated with the use of such analysis to predict alpha-helix structure/function relationships are considered.     

Reversible and irreversible modifications ofβ-lactoglobulin upon exposure to heat

Modifications in the exposure to the solvent of hydrophobic residues, changes in their organization into surface hydrophobic patches, and alterations in the dimerization equilibrium ofβ-lactoglobulin upon thermal treatment at neutralpH were studied. Exposure of tryptophan residues was temperature dependent and was essentially completed on the time scale of seconds. Reorganization of generic hydrophobic protein patches on the protein surface was monitored through binding of 1,8-anilinonaphthalenesulfonate, and was much slower than changes in tryptophan exposure. Different phases in surface hydrophobicity changes were related to the swelling and the subsequent collapse of the protein, which formed a metastable swollen intermediate. Heat treatment ofβ-lactoglobulin also resulted in the formation of soluble oligomeric aggregates. The aggregation process was studied as a function of temperature, demonstrating that (i) dimer dissociation was a necessary step in a sequential polymerization mechanism and (ii) cohesion of hydrophobic patches was the major driving force for aggregation.     

Mutation matrices and physical-chemical properties: Correlations and implications

To investigate how the properties of individual amino acids result in proteins with particular structures and functions, we have examined the correlations between previously derived structure-dependent mutation rates and changes in various physical-chemical properties of the amino acids such as volume, charge, α-helical and β-sheet propensity, and hydrophobicity. In most cases we found the ΔG of transfer from octanol to water to be the best model for evolutionary constraints, in contrast to the much weaker correlation with the ΔG of transfer from cyclohexane to water, a property found to be highly correlated to changes in stability in site-directed mutagenesis studies. This suggests that natural evolution may follow different rules than those suggested by results obtained in the laboratory. A high degree of conservation of a surface residue's relative hydrophobicity was also observed, a fact that cannot be explained by constraints on protein stability but that may reflect the consequences of the reverse-hydrophobic effect. Local propensity, especially α-helical propensity, is rather poorly conserved during evolution, indicating that non-local interactions dominate protein structure formation. We found that changes in volume were important in specific cases, most significantly in transitions among the hydrophobic residues in buried locations. To demonstrate how these techniques could be used to understand particular protein families, we derived and analyzed mutation matrices for the hypervariable and framework regions of antibody light chain V regions. We found a surprisingly high conservation of hydrophobicity in the hypervariable region, possibly indicating an important role for hydrophobicity in antigen recognition. Proteins 27:336–344, 1997. © 1997 Wiley-Liss, Inc.     

Effect of site-directed point mutations on protein misfolding: A simulation study

A Monte Carlo simulation based sequence design method is proposed to investigate the role of site-directed point mutations in protein misfolding. Site-directed point mutations are incorporated in the designed sequences of selected proteins. While most mutated sequences correctly fold to their native conformation, some of them stabilize in other nonnative conformations and thus misfold/unfold. The results suggest that a critical number of hydrophobic amino acid residues must be present in the core of the correctly folded proteins, whereas proteins misfold/unfold if this number of hydrophobic residues falls below the critical limit. A protein can accommodate only a particular number of hydrophobic residues at the surface, provided a large number of hydrophilic residues are present at the surface and critical hydrophobicity of the core is preserved. Some surface sites are observed to be equally sensitive toward site-directed point mutations as the core sites. Point mutations with highly polar and charged amino acids increases the misfold/unfold propensity of proteins. Substitution of natural amino acids at sites with different number of nonbonded contacts suggests that both amino acid identity and its respective site-specificity determine the stability of a protein. A clash-match method is developed to calculate the number of matching and clashing interactions in the mutated protein sequences. While misfolded/unfolded sequences have a higher number of clashing and a lower number of matching interactions, the correctly folded sequences have a lower number of clashing and a higher number of matching interactions. These results are valid for different SCOP classes of proteins.