Role of hydrophobic clusters and long-range contact networks in the folding of (alpha/beta)8 barrel proteins

Importance of Inter-residue Contacts for Understanding Protein Folding and Unfolding Rates, Remote Homology, and Drug Design

Inter-residue interactions in protein structures provide valuable insights into protein folding and stability. Understanding these interactions can be helpful in many crucial applications, including rational design of therapeutic small molecules and biologics, locating functional protein sites, and predicting protein-protein and protein-ligand interactions. The process of developing machine learning models incorporating inter-residue interactions has been improved recently. This review highlights the theoretical models incorporating inter-residue interactions in predicting folding and unfolding rates of proteins. Utilizing contact maps to depict inter-residue interactions aids researchers in developing computer models for detecting remote homologs and interface residues within protein-protein complexes which, in turn, enhances our knowledge of the relationship between sequence and structure of proteins. Further, the application of contact maps derived from inter-residue interactions is highlighted in the field of drug discovery. Overall, this review presents an extensive assessment of the significant models that use inter-residue interactions to investigate folding rates, unfolding rates, remote homology, and drug development, providing potential future advancements in constructing efficient computational models in structural biology.

De Novo Design of a Highly Stable Ovoid TIM Barrel: Unlocking Pocket Shape towards Functional Design

The ability to finely control the structure of protein folds is an important prerequisite to functional protein design. The TIM barrel fold is an important target for these efforts as it is highly enriched for diverse functions in nature. Although a TIM barrel protein has been designed de novo, the ability to finely alter the curvature of the central beta barrel and the overall architecture of the fold remains elusive, limiting its utility for functional design. Here, we report the de novo design of a TIM barrel with ovoid (twofold) symmetry, drawing inspiration from natural beta and TIM barrels with ovoid curvature. We use an autoregressive backbone sampling strategy to implement our hypothesis for elongated barrel curvature, followed by an iterative enrichment sequence design protocol to obtain sequences which yield a high proportion of successfully folding designs. Designed sequences are highly stable and fold to the designed barrel curvature as determined by a 2.1 Å resolution crystal structure. The designs show robustness to drastic mutations, retaining high melting temperatures even when multiple charged residues are buried in the hydrophobic core or when the hydrophobic core is ablated to alanine. As a scaffold with a greater capacity for hosting diverse hydrogen bonding networks and installation of binding pockets or active sites, the ovoid TIM barrel represents a major step towards the de novo design of functional TIM barrels.

The Stability Landscape of de novo TIM Barrels Explored by a Modular Design Approach

•

The TIM barrel is a versatile fold to understand structure-stability relationships.

•

A collection of de novo TIM barrels with improved hydrophobic cores was designed.

•

DeNovoTIMs are reversible in chemical and thermal unfolding, which is uncommon in TIM barrels.

•

Epistatic effects play a central role in DeNovoTIMs stabilization.

•

DeNovoTIMs navigate a previously uncharted region of the stability landscape.

The ability to design stable proteins with custom-made functions is a major goal in biochemistry with practical relevance for our environment and society. Understanding and manipulating protein stability provide crucial information on the molecular determinants that modulate structure and stability, and expand the applications of de novo proteins. Since the (β/⍺)₈-barrel or TIM-barrel fold is one of the most common functional scaffolds, in this work we designed a collection of stable de novo TIM barrels (DeNovoTIMs), using a computational fixed-backbone and modular approach based on improved hydrophobic packing of sTIM11, the first validated de novo TIM barrel, and subjected them to a thorough folding analysis. DeNovoTIMs navigate a region of the stability landscape previously uncharted by natural TIM barrels, with variations spanning 60 degrees in melting temperature and 22 kcal per mol in conformational stability throughout the designs. Significant non-additive or epistatic effects were observed when stabilizing mutations from different regions of the barrel were combined. The molecular basis of epistasis in DeNovoTIMs appears to be related to the extension of the hydrophobic cores. This study is an important step towards the fine-tuned modulation of protein stability by design.

Network Connectivity, Centrality and Fragmentation in the Greek-Key Protein Topology

Understanding and computationally predicting the protein folding process remains one of the most challenging scientific problems and has uniquely garnered the interdisciplinary efforts of researchers from both the biological, chemical, physical and computational disciplines. Previous studies have demonstrated the importance of long-range interactions in guiding the native structure. However, predicting how the native long-range interaction network forms to generate a specific topology from among all other conformations remains unresolved. The present research study conducts an exploratory study to identify amino acids and long-range interactions that have the potential to play a key role in building and maintaining the protein topology. Towards this end, the application of network science is utilized and developed to analyze the structures of a group of proteins that share a common Greek-key topology but differ in sequence, secondary structure and function. We investigate the idea that the residues with high betweeness centrality score are potentially significant in maintaining the protein network and in governing the Greek-key topology. This hypothesis is tested by two different computational methods: through a fragmentation test and by the analysis of diameter impacts. In summary, we find a subset of selected residues in similar geographical positions in all model proteins, which demonstrates the role of these specific residues and regions in governing the Greek-key topology from a network perspective.

Changes in hydrophobicity mainly promotes the aggregation tendency of ALS associated SOD1 mutants

Protein misfolding and aggregation due to mutations, are associated with fatal neurodegenerative disorders. The mutations in Cu/Zn superoxide dismutase (SOD1) causing its misfolding and aggregation are found linked to the motor neuron disorder, amyotrophic lateral sclerosis. Since the mutations are scattered throughout SOD1 structure, determining the exact molecular mechanism underlying the ALS pathology remains unresolved. In this study, we have investigated the major molecular factors that mainly contribute to SOD1 destabilization, intrinsic disorder, and misfolding using sequence and structural information. We have analysed 153 ALS causing SOD1 point mutants for aggregation tendency using four different aggregation prediction tools, viz., Aggrescan3D (A3D), CamSol, GAP and Zyggregator. Our results suggest that 74-79 mutants are susceptible to aggregation, due to distorted native interactions originated at the mutation site. Majority of the aggregation prone mutants are located in the buried regions of SOD1 molecule. Further, the mutations at the hydrophobic amino acids primarily promote the aggregation tendency of SOD1 protein through different destabilizing mechanisms including changes in hydrophobic free energy, loss of electrostatic interactions in the protein's surface and loss of hydrogen bonds that bridges the protein core and surface.

Pairwise contact energy statistical potentials can help to find probability of point mutations

To adopt a particular fold, a protein requires several interactions between its amino acid residues. The energetic contribution of these residue-residue interactions can be approximated by extracting statistical potentials from known high resolution structures. Several methods based on statistical potentials extracted from unrelated proteins are found to make a better prediction of probability of point mutations. We postulate that the statistical potentials extracted from known structures of similar folds with varying sequence identity can be a powerful tool to examine probability of point mutation. By keeping this in mind, we have derived pairwise residue and atomic contact energy potentials for the different functional families that adopt the (α/β)₈ TIM-Barrel fold. We carried out computational point mutations at various conserved residue positions in yeast Triose phosphate isomerase enzyme for which experimental results are already reported. We have also performed molecular dynamics simulations on a subset of point mutants to make a comparative study. The difference in pairwise residue and atomic contact energy of wildtype and various point mutations reveals probability of mutations at a particular position. Interestingly, we found that our computational prediction agrees with the experimental studies of Silverman et al. (Proc Natl Acad Sci 2001;98:3092-3097) and perform better prediction than i_Mutant and Cologne University Protein Stability Analysis Tool. The present work thus suggests deriving pairwise contact energy potentials and molecular dynamics simulations of functionally important folds could help us to predict probability of point mutations which may ultimately reduce the time and cost of mutation experiments. Proteins 2016; 85:54-64. © 2016 Wiley Periodicals, Inc.

Evaluation of thermo-stability of bluetongue virus recombinant VP7 antigen in indirect ELISA

This study shows the thermo-stability of lyophilized and purified recombinant VP7 bluetongue virus (BTV) protein in the presence of two sugar stabilizers (trehalose and mannitol) at different temperature. Truncated VP7 protein purified by nickel affinity column was lyophilized in the presence of trehalose and mannitol at 60 mM final concentration and then exposed to different temperature like 4, 25, 37 and 45 °C for various periods like 5 months, 7 weeks, 7 days and 48 h, respectively. After thermal treatment, the reactivity of the protein was evaluated in indirect ELISA. At 4 and 25 °C, the protein was stable up to 5 months and 7 weeks, respectively, irrespective of stabilizers used. At 37 °C, it was stable up to 3 days with both the stabilizers, after which it lost its stability and reactivity. At 45 °C, the protein was stable up to 30 and 24 h with trehalose and mannitol stabilizers, respectively. Both stabilizers found suitable for stability of the protein. However, trehalose appeared to have better stabilizing effect, particularly at higher temperatures than the mannitol. Trehalose could be used as stabilizer for freeze-drying the recombinant VP7 protein if an indirect ELISA kit based on the purified rVP7 protein is supplied to different laboratories of the country for detection of BTV antibody in sheep.

Specific non-local interactions are not necessary for recovering native protein dynamics

The elastic network model (ENM) is a widely used method to study native protein dynamics by normal mode analysis (NMA). In ENM we need information about all pairwise distances, and the distance between contacting atoms is restrained to the native value. Therefore ENM requires O(N²) information to realize its dynamics for a protein consisting of N amino acid residues. To see if (or to what extent) such a large amount of specific structural information is required to realize native protein dynamics, here we introduce a novel model based on only O(N) restraints. This model, named the ‘contact number diffusion’ model (CND), includes specific distance restraints for only local (along the amino acid sequence) atom pairs, and semi-specific non-local restraints imposed on each atom, rather than atom pairs. The semi-specific non-local restraints are defined in terms of the non-local contact numbers of atoms. The CND model exhibits the dynamic characteristics comparable to ENM and more correlated with the explicit-solvent molecular dynamics simulation than ENM. Moreover, unrealistic surface fluctuations often observed in ENM were suppressed in CND. On the other hand, in some ligand-bound structures CND showed larger fluctuations of buried protein atoms interacting with the ligand compared to ENM. In addition, fluctuations from CND and ENM show comparable correlations with the experimental B-factor. Although there are some indications of the importance of some specific non-local interactions, the semi-specific non-local interactions are mostly sufficient for reproducing the native protein dynamics.

Investigations on the role of π-π interactions and π-π networks in eNOS and nNOS proteins

π-π Interactions play an important role in the stability of protein structures. In the present study, we have analyzed the influence of π-π interactions in eNOS and nNOS proteins. The contribution of these π-π interacting residues in sequential separation, secondary structure involvement, solvent accessibility and stabilization centers has been evaluated. π-π interactions stabilize the core regions within eNOS and nNOS proteins. π-π interacting residues are evolutionary conserved. There is a significant number of π-π interactions in spite of the lesser natural occurrences of π-residues in eNOS and nNOS proteins. In addition to π-π interactions, π residues also form π-π networks in both eNOS and nNOS proteins which might play an important role in the structural stability of these protein structures.

Aromatic-aromatic interactions: analysis of π-π interactions in interleukins and TNF proteins

Aromatic-aromatic hydrogen bonds are important in many areas of chemistry, biology and materials science. In this study we have analyzed the roles played by the π-π interactions in interleukins (ILs) and tumor necrosis factor (TNF) proteins. Majority of π-π interacting residues are conserved in ILs and TNF proteins. The accessible surface area calculations in these proteins reveal that these interactions might be important in stabilizing the inner core regions of these proteins. In addition to π-π interactions, the aromatic residues also form π-networks in ILs and TNF proteins. The results obtained in the present study indicate that π-π interactions and π-π networks play important roles in the structural stability of ILs and TNF proteins.

Beta atomic contacts: identifying critical specific contacts in protein binding interfaces

Specific binding between proteins plays a crucial role in molecular functions and biological processes. Protein binding interfaces and their atomic contacts are typically defined by simple criteria, such as distance-based definitions that only use some threshold of spatial distance in previous studies. These definitions neglect the nearby atomic organization of contact atoms, and thus detect predominant contacts which are interrupted by other atoms. It is questionable whether such kinds of interrupted contacts are as important as other contacts in protein binding. To tackle this challenge, we propose a new definition called beta (β) atomic contacts. Our definition, founded on the β-skeletons in computational geometry, requires that there is no other atom in the contact spheres defined by two contact atoms; this sphere is similar to the van der Waals spheres of atoms. The statistical analysis on a large dataset shows that β contacts are only a small fraction of conventional distance-based contacts. To empirically quantify the importance of β contacts, we design βACV, an SVM classifier with β contacts as input, to classify homodimers from crystal packing. We found that our βACV is able to achieve the state-of-the-art classification performance superior to SVM classifiers with distance-based contacts as input. Our βACV also outperforms several existing methods when being evaluated on several datasets in previous works. The promising empirical performance suggests that β contacts can truly identify critical specific contacts in protein binding interfaces. β contacts thus provide a new model for more precise description of atomic organization in protein quaternary structures than distance-based contacts.

Protein fragments: functional and structural roles of their coevolution networks

Small protein fragments, and not just residues, can be used as basic building blocks to reconstruct networks of coevolved amino acids in proteins. Fragments often enter in physical contact one with the other and play a major biological role in the protein. The nature of these interactions might be multiple and spans beyond binding specificity, allosteric regulation and folding constraints. Indeed, coevolving fragments are indicators of important information explaining folding intermediates, peptide assembly, key mutations with known roles in genetic diseases, distinguished subfamily-dependent motifs and differentiated evolutionary pressures on protein regions. Coevolution analysis detects networks of fragments interaction and highlights a high order organization of fragments demonstrating the importance of studying at a deeper level this structure. We demonstrate that it can be applied to protein families that are highly conserved or represented by few sequences, enlarging in this manner, the class of proteins where coevolution analysis can be performed and making large-scale coevolution studies a feasible goal.

Arginine and Lysine interactions with π residues in metalloproteins

Metalloproteins have many different functions in cells such as enzymes; signal transduction, transport and storage proteins. About one third of all proteins require metals to carry out their functions. In the present study we have analyzed the roles played by Arg and Lys (cationic side chains) interactions with π (Phe, Tyr or Trp) residues and their role in the structural stability of metalloproteins. These interactions might play an important role in the global conformational stability in metalloproteins. In spite of its lower natural occurrence (1.76%) the number of Trp residues involved in energetically significant interactions is higher in metalloproteins.

Role of long- and short-range hydrophobic, hydrophilic and charged residues contact network in protein's structural organization

Background

The three-dimensional structure of a protein can be described as a graph where nodes represent residues and the strength of non-covalent interactions between them are edges. These protein contact networks can be separated into long and short-range interactions networks depending on the positions of amino acids in primary structure. Long-range interactions play a distinct role in determining the tertiary structure of a protein while short-range interactions could largely contribute to the secondary structure formations. In addition, physico chemical properties and the linear arrangement of amino acids of the primary structure of a protein determines its three dimensional structure. Here, we present an extensive analysis of protein contact subnetworks based on the London van der Waals interactions of amino acids at different length scales. We further subdivided those networks in hydrophobic, hydrophilic and charged residues networks and have tried to correlate their influence in the overall topology and organization of a protein.

Results

The largest connected component (LCC) of long (LRN)-, short (SRN)- and all-range (ARN) networks within proteins exhibit a transition behaviour when plotted against different interaction strengths of edges among amino acid nodes. While short-range networks having chain like structures exhibit highly cooperative transition; long- and all-range networks, which are more similar to each other, have non-chain like structures and show less cooperativity. Further, the hydrophobic residues subnetworks in long- and all-range networks have similar transition behaviours with all residues all-range networks, but the hydrophilic and charged residues networks don’t. While the nature of transitions of LCC’s sizes is same in SRNs for thermophiles and mesophiles, there exists a clear difference in LRNs. The presence of larger size of interconnected long-range interactions in thermophiles than mesophiles, even at higher interaction strength between amino acids, give extra stability to the tertiary structure of the thermophiles. All the subnetworks at different length scales (ARNs, LRNs and SRNs) show assortativity mixing property of their participating amino acids. While there exists a significant higher percentage of hydrophobic subclusters over others in ARNs and LRNs; we do not find the assortative mixing behaviour of any the subclusters in SRNs. The clustering coefficient of hydrophobic subclusters in long-range network is the highest among types of subnetworks. There exist highly cliquish hydrophobic nodes followed by charged nodes in LRNs and ARNs; on the other hand, we observe the highest dominance of charged residues cliques in short-range networks. Studies on the perimeter of the cliques also show higher occurrences of hydrophobic and charged residues’ cliques.

Conclusions

The simple framework of protein contact networks and their subnetworks based on London van der Waals force is able to capture several known properties of protein structure as well as can unravel several new features. The thermophiles do not only have the higher number of long-range interactions; they also have larger cluster of connected residues at higher interaction strengths among amino acids, than their mesophilic counterparts. It can reestablish the significant role of long-range hydrophobic clusters in protein folding and stabilization; at the same time, it shed light on the higher communication ability of hydrophobic subnetworks over the others. The results give an indication of the controlling role of hydrophobic subclusters in determining protein’s folding rate. The occurrences of higher perimeters of hydrophobic and charged cliques imply the role of charged residues as well as hydrophobic residues in stabilizing the distant part of primary structure of a protein through London van der Waals interaction.

EXPLORING THE ROLE OF C–H … π INTERACTIONS ON THE STRUCTURAL STABILITY OF ANTIMICROBIAL PEPTIDES

A computational analysis on the C – H … π interactions in a group of 53 antimicrobial peptides was investigated. A total of 162 C – H … π interactions were observed. Side-chain to side-chain C – H … π interactions are the predominant type of interactions in antimicrobial peptides data set. There was an average of one significant C – H … π interaction for every 7 residues in the antimicrobial peptides investigated. Long-range C – H … π interactions are the predominant type of interactions. The secondary structure preference, solvent accessibility and stabilization centers of these of C – H … π interacting residues were estimated. It is likely that the C – H … π interactions contribute significantly to the overall stability of antimicrobial peptides. These interactions were observed after a molecular dynamics study on these set of antimicrobial peptides using CHARMM force field.

Structural stability studies in adhesion molecules--role of cation-π interactions

Cell adhesion molecules are important for their various roles in many cellular events and responses. In the present study, we have analyzed the roles played by cation-π interactions in the structural stability of adhesion molecules. These interactions are mainly formed by long-range contacts. The occurrence of arginine is higher than lysine to form cation-π interactions. The secondary structure preferences of interacting residues are independent of amino acid class. Cation-π interactions might stabilize the interface between the terminus and core in this class of proteins. The results obtained in the present study will be useful in understanding the contribution of cation-π interactions to the overall stability of adhesion proteins.

Fractal symmetry of protein interior: what have we learned?

The application of fractal dimension-based constructs to probe the protein interior dates back to the development of the concept of fractal dimension itself. Numerous approaches have been tried and tested over a course of (almost) 30 years with the aim of elucidating the various facets of symmetry of self-similarity prevalent in the protein interior. In the last 5 years especially, there has been a startling upsurge of research that innovatively stretches the limits of fractal-based studies to present an array of unexpected results on the biophysical properties of protein interior. In this article, we introduce readers to the fundamentals of fractals, reviewing the commonality (and the lack of it) between these approaches before exploring the patterns in the results that they produced. Clustering the approaches in major schools of protein self-similarity studies, we describe the evolution of fractal dimension-based methodologies. The genealogy of approaches (and results) presented here portrays a clear picture of the contemporary state of fractal-based studies in the context of the protein interior. To underline the utility of fractal dimension-based measures further, we have performed a correlation dimension analysis on all of the available non-redundant protein structures, both at the level of an individual protein and at the level of structural domains. In this investigation, we were able to separately quantify the self-similar symmetries in spatial correlation patterns amongst peptide-dipole units, charged amino acids, residues with the π-electron cloud and hydrophobic amino acids. The results revealed that electrostatic environments in the interiors of proteins belonging to 'α/α toroid' (all-α class) and 'PLP-dependent transferase-like' domains (α/β class) are highly conducive. In contrast, the interiors of 'zinc finger design' ('designed proteins') and 'knottins' ('small proteins') were identified as folds with the least conducive electrostatic environments. The fold 'conotoxins' (peptides) could be unambiguously identified as one type with the least stability. The same analyses revealed that peptide-dipoles in the α/β class of proteins, in general, are more correlated to each other than are the peptide-dipoles in proteins belonging to the all-α class. Highly favorable electrostatic milieu in the interiors of TIM-barrel, α/β-hydrolase structures could explain their remarkably conserved (evolutionary) stability from a new light. Finally, we point out certain inherent limitations of fractal constructs before attempting to identify the areas and problems where the implementation of fractal dimension-based constructs can be of paramount help to unearth latent information on protein structural properties.

Search for identical octapeptides in unrelated proteins: Structural plasticity revisited

Since proteins are dynamic in nature, they can alter their local structure in response to changes in their environment factors such as temperature, pH, phosphorylation, and binding of other small molecules. These conformational changes are extremely important for the correct folding and functioning of proteins. There are also a number of diseases associated with protein conformational change such as amyloid diseases. To stimulate research into the above factors which specify one conformation over another, different theoretical models have been proposed and tested against sequence similar distant structure protein fragments. In order to simplify the computational complexity of identifying conformational changes in proteins, various local sequence search algorithms were employed and the structural plasticity in unrelated proteins was examined by various research groups. In the present work, we revisit the mechanism of structural plasticity in unrelated proteins with increased number of structures in Protein Data Bank by comparing identical octapeptides in unrelated proteins with dictionary of protein secondary structure extracted from existing experimental data. Our goal is to bring out the influence of hydrophobic residues, hydrophilic residues, flanking residues, difference in secondary structural propensities of surrounding residues, difference in phi-psi angles and local and nonlocal interactions in identical octapeptides adopting different conformations. Also we have used surrounding hydrophobicity, environment dependent interaction energy, atomic mean force potential, structural unit contacts and difference profiles models to explore the factors which cause structural plasticity. The results discussed here may provide insights into protein folding, design and function.

Characterizing the regularity of tetrahedral packing motifs in protein tertiary structure

MOTIVATION

While protein secondary structure is well understood, representing the repetitive nature of tertiary packing in proteins remains difficult. We have developed a construct called the relative packing group (RPG) that applies the clique concept from graph theory as a natural basis for defining the packing motifs in proteins. An RPG is defined as a clique of residues, where every member contacts all others as determined by the Delaunay tessellation. Geometrically similar RPGs define a regular element of tertiary structure or tertiary motif (TerMo). This intuitive construct provides a simple approach to characterize general repetitive elements of tertiary structure.

RESULTS

A dataset of over 4 million tetrahedral RPGs was clustered using different criteria to characterize the various aspects of regular tertiary structure in TerMos. Grouping this data within the SCOP classification levels of Family, Superfamily, Fold, Class and PDB showed that similar packing is shared across different folds. Classification of RPGs based on residue sequence locality reveals topological preferences according to protein sizes and secondary structure. We find that larger proteins favor RPGs with three local residues packed against a non-local residue. Classifying by secondary structure, helices prefer mostly local residues, sheets favor at least two local residues, while turns and coil populate with more local residues. To depict these TerMos, we have developed 2 complementary and intuitive representations: (i) Dirichlet process mixture density estimation of the torsion angle distributions and (ii) kernel density estimation of the Cartesian coordinate distribution. The TerMo library and representations software are available upon request.

Sequence and structural analysis of two designed proteins with 88% identity adopting different folds

Protein folding is a natural phenomenon by which a sequence of amino acids folds into a unique functional three-dimensional structure. Although the sequence code that governs folding remains a mystery, one can identify key inter-residue contacts responsible for a given topology. In nature, there are many pairs of proteins of a given length that share little or no sequence identity. Similarly, there are many proteins that share a common topology but lack significant evidence of homology. In order to tackle this problem, protein engineering studies have been used to determine the minimal number of amino acid residues that codes for a particular fold. In recent years, the coupling of theoretical models and experiments in the study of protein folding has resulted in providing some fruitful clues. He et al. have designed two proteins with 88% sequence identity, which adopt different folds and functions. In this work, we have systematically analysed these two proteins by performing pentapeptide search, secondary structure predictions, variation in inter-residue interactions and residue-residue pair preferences, surrounding hydrophobicity computations, conformational switching and energy computations. We conclude that the local secondary structural preference of the two designed proteins at the Nand C-terminal ends to adopt either coil or strand conformation may be a crucial factor in adopting the different folds. Early on during the process of folding, both proteins may choose different energetically favourable pathways to attain the different folds.

Influence of cation-pi interactions to the structural stability of prokaryotic and eukaryotic translation elongation factors

We have investigated the role of cation-pi interactions on translation elongation factors. In our investigation, an average of four significant cation-pi interactions were found, that is, an average of one cation-pi interaction per 44 residues in the ten elongation factors were observed. The analysis on the influence of short (< + or - 4), medium (> + or - 4 to < + or - 20) and long (>20) range contacts showed that cation-pi interactions are mainly formed by medium and long-range contacts. Arg-Tyr pair was found largest in number but energetic contribution of Arg-Trp pair was found most. Preferred secondary structural conformation analysis of the residues involved in cation-pi interaction indicates that the cationic Arg prefers to be in helix and Lys having equal probability for helix and strand, whereas the aromatic Phe and Trp were found mostly in helix while Tyr in strand regions. The cation-pi interaction residues involved in these proteins were found highly conserved with 48.86% residues having conservation score of > or = 6. Analysis of secondary structure preference of the energetically significant cation-pi residues in different solvent accessible range indicates that most of the pi residues are found buried or partially buried whereas cationic residues were found mostly at the protein surface. The results presented in this study will be useful for structural stability studies in translation elongation factors.

Prediction of inter-residue contact clusters from hydrophobic cores

A contact map is a key factor representing a specific protein structure. To simplify the protein contact map prediction, we predict the inter-residue contact clusters centered at the groups of their surrounding inter-residue contacts. In this paper, we adopt a Support Vector Machine (SVM)-based approach to predict the inter-residue contact cluster centers. The input of the SVM predictor includes sequence profile, evolutionary rate and predicted secondary structure. The SVM predictor is based on hydrophobic cores that may be considered as locations of the inter-residue contact clusters. About 35% of clustering centers of inter-residue contacts can be predicted accurately.

A role for confined water in chaperonin function

Chaperonins engulf other proteins and accelerate their folding by an unknown mechanism. Here, we combine all-atom molecular dynamics simulations with data from experimental assays of the activity of the bacterial chaperonin GroEL to demonstrate that a chaperonin’s ability to facilitate folding is correlated with the affinity of its interior surface for water. Our results suggest a novel view of the behavior of confined water for models of in vivo protein folding scenarios.

Investigations on unconventional hydrogen bonds in RNA binding proteins: The role of CH⋯OC interactions

We have investigated the roles played by C-H...O=C interactions in RNA binding proteins. There was an average of 78 CH...O=C interactions per protein and also there was an average of one significant CH...O=C interactions for every 6 residues in the 59 RNA binding proteins studied. Main chain-Main chain (MM) CH...O=C interactions are the predominant type of interactions in RNA binding proteins. The donor atom contribution to CH...O=C interactions was mainly from aliphatic residues. The acceptor atom contribution for MM CH...O=C interactions was mainly from Val, Phe, Leu, Ile, Arg and Ala. The secondary structure preference analysis of CH...O=C interacting residues showed that, Arg, Gln, Glu and Tyr preferred to be in helix, while Ala, Asp, Cys, Gly, Ile, Leu, Lys, Met, Phe, Trp and Val preferred to be in strand conformation. Most of the CH...O=C interacting polar amino acid residues were solvent exposed while, majority of the CH...O=C interacting non polar residues were excluded from the solvent. Long and medium-range CH...O=C interactions are the predominant type of interactions in RNA binding proteins. More than 50% of CH...O=C interacting residues had a higher conservation score. Significant percentage of CH...O=C interacting residues had one or more stabilization centers. Sixty-six percent of the theoretically predicted stabilizing residues were also involved in CH...O=C interactions and hence these residues may also contribute additional stability to RNA binding proteins.

Investigations on C–H⋯π interactions in RNA binding proteins

We have investigated the roles played by C-Hcdots, three dots, centeredpi interactions in RNA binding proteins. There was an average of 55 C-Hcdots, three dots, centeredpi interactions per protein and also there was an average of one significant C-Hcdots, three dots, centeredpi interaction for every nine residues in the 59 RNA binding proteins studied. Main-chain to side-chain C-Hcdots, three dots, centeredpi interactions is the predominant type of interactions in RNA binding proteins. The donor atom contribution to C-Hcdots, three dots, centeredpi interactions was mainly from Phe, Tyr, Trp, Pro, Gly, Lys, His and Ala residues. The acceptor atom contribution to main-chain to side-chain C-Hcdots, three dots, centeredpi and side-chain to side-chain C-Hcdots, three dots, centeredpi interactions was mainly from Phe and Tyr residues. On the contrary, the acceptor atoms of Trp residues contributed to all the four types of C-Hcdots, three dots, centeredpi interactions. Also, Trp contributed both donor and acceptor atoms in main-chain to side-chain, main-chain to side-chain five-member aromatic ring and side-chain to side-chain C-Hcdots, three dots, centeredpi interactions. The secondary structure preference analysis of C-Hcdots, three dots, centeredpi interacting residues showed that, Arg, Gln, Glu, His, Ile, Leu, Lys, Met, Phe and Tyr preferred to be in helix, while Ala, Asp, Cys, Gly, Trp and Val preferred to be in strand conformation. Long-range C-Hcdots, three dots, centeredpi interactions are the predominant type of interactions in RNA binding proteins. More than 50% of C-Hcdots, three dots, centeredpi interacting residues had a higher conservation score. Significant percentage of C-Hcdots, three dots, centeredpi interacting residues had one or more stabilization centers. Seven percent of the theoretically predicted stabilizing residues were also involved in C-Hcdots, three dots, centeredpi interactions and hence these residues may also contribute additional stability to RNA binding proteins.

Influence of cation-π interactions on RNA-binding proteins

The energy contribution due to cation-pi interactions has been computed for 37 RNA binding proteins. The contribution of these cation-pi interacting residues in sequential separation, secondary structure involvement, solvent accessibility, and stabilization centers has been evaluated. Sequential separation of the cation-pi involving residues show that, long range contacts predominates in all the proteins studied. Lys and Arg prefers to be in helical structures. Of the cation-pi interacting residues, Arg and Lys were in the exposed regions and the aromatic residues (Phe, Tyr and Trp) were in the buried and partially buried regions in the protein structures. Stabilization centers for these proteins showed that all the five residues found in cation-pi interactions are important in locating one or more of such centers. On the whole, the results presented in this work will be very useful for further investigations on the specificity and selectivity of RNA binding proteins and also for their structural studies.

Periodic distributions of hydrophobic amino acids allows the definition of fundamental building blocks to align distantly related proteins

Several studies on large and small families of proteins proved in a general manner that hydrophobic amino acids are globally conserved even if they are subjected to high rate substitution. Statistical analysis of amino acids evolution within blocks of hydrophobic amino acids detected in sequences suggests their usage as a basic structural pattern to align pairs of proteins of less than 25% sequence identity, with no need of knowing their 3D structure. The authors present a new global alignment method and an automatic tool for Proteins with HYdrophobic Blocks ALignment (PHYBAL) based on the combinatorics of overlapping hydrophobic blocks. Two substitution matrices modeling a different selective pressure inside and outside hydrophobic blocks are constructed, the Inside Hydrophobic Blocks Matrix and the Outside Hydrophobic Blocks Matrix, and a 4D space of gap values is explored. PHYBAL performance is evaluated against Needleman and Wunsch algorithm run with Blosum 30, Blosum 45, Blosum 62, Gonnet, HSDM, PAM250, Johnson and Remote Homo matrices. PHYBAL behavior is analyzed on eight randomly selected pairs of proteins of >30% sequence identity that cover a large spectrum of structural properties. It is also validated on two large datasets, the 127 pairs of the Domingues dataset with >30% sequence identity, and 181 pairs issued from BAliBASE 2.0 and ranked by percentage of identity from 7 to 25%. Results confirm the importance of considering substitution matrices modeling hydrophobic contexts and a 4D space of gap values in aligning distantly related proteins. Two new notions of local and global stability are defined to assess the robustness of an alignment algorithm and the accuracy of PHYBAL. A new notion, the SAD-coefficient, to assess the difficulty of structural alignment is also introduced. PHYBAL has been compared with Hydrophobic Cluster Analysis and HMMSUM methods.

Computation of non-covalent interactions in TNF proteins and interleukins

The roles played by the non-covalent interactions have been investigated for a set of six TNF proteins and nine Interleukins. The stabilizing residues have been identified by a consensus approach using the concepts of available surface area, medium and long-range interactions and conservation of amino acid residues. The cation-pi interactions have been computed based on a geometric approach such as distance and energy criteria. We identified an average of 1 energetically significant cation-pi interactions in every 94 residues in TNF proteins and 1 in every 62 residues in Interleukins. In TNF proteins, the cationic groups Lys preferred to be in helix while Arg preferred to be in strand regions while in Interleukins the Arg residues preferred to be in helix and Lys preferred to be in strand regions. From the available surface area calculations, we found that, almost all the cation and pi residues in TNF proteins and Interleukins were either in buried or partially buried regions and none of them in the exposed regions. Medium and long-range interactions were predominant in both TNF proteins and Interleukins. It was observed that the percentage of stabilizing centers were more in TNF proteins as compared to the Interleukins, while the percentage of conserved residues were more in Interleukins than in TNF proteins. In the stabilizing residues Lys was observed to be a stabilizing residue in both TNF proteins and Interleukins. Among the aromatic group, Phe was seen to be a stabilizing residue in both TNF and Interleukins. We suggest that this study on the computation of cation-pi interactions in TNF proteins and Interleukins would be very helpful in further understanding the structure, stability and functional similarity of these proteins.

Structural analysis of cation-pi interactions in DNA binding proteins

Cation-pi interactions play an important role in the stability of protein structures. In this work, we have analyzed the influence of cation-pi interactions in DNA binding proteins. We observed cation-pi interactions in 45 out of 62 DNA binding proteins and there is no significant correlation between the number of amino acid residues and number of cation-pi interactions. These interactions are mainly formed by long-range contacts, and the role of short and medium-range contacts is minimal. The preference of Arg is higher than Lys to form cation-pi interactions. The pair-wise cation-pi interaction energy between aromatic and positively charged residues shows that Arg-Tyr energy is the strongest among the possible six pairs. The structural analysis of cation-pi interaction forming residues shows that Lys, Trp, and Tyr prefer to be in the binding site of protein-DNA complexes. Further, the accessible surface areas of cation-pi interaction forming cationic residues are significantly less than that of other residues. The preference of cation-pi interaction forming residues in different secondary structures shows that Lys prefers to be in strand and Phe prefers to be in turn regions. The results obtained in the present study will be useful in understanding the contribution of cation-pi interactions to the stability and specificity of protein-DNA complexes.

Folding Mechanism of All‐β Globular Proteins

Elucidating the mechanism for the fast folding of proteins is a challenging task. In our earlier work, we introduced the concept of "long-range order" and related it successfully to protein folding rates. In this article, we propose a new hypothesis for the folding of two-state all-beta proteins. The mechanism is based on the formation of a hydrophobic core, propagation of beta-strands, and the establishment of hydrogen bonds. Our hypothesis has been strengthened by the observation of a folding nucleus in beta-strands and the hydrogen-bonding network between residues in beta-strands. Our insights on protein folding show an excellent agreement with experimental observations.

Inter-residue interactions in protein folding and stability

During the process of protein folding, the amino acid residues along the polypeptide chain interact with each other in a cooperative manner to form the stable native structure. The knowledge about inter-residue interactions in protein structures is very helpful to understand the mechanism of protein folding and stability. In this review, we introduce the classification of inter-residue interactions into short, medium and long range based on a simple geometric approach. The features of these interactions in different structural classes of globular and membrane proteins, and in various folds have been delineated. The development of contact potentials and the application of inter-residue contacts for predicting the structural class and secondary structures of globular proteins, solvent accessibility, fold recognition and ab initio tertiary structure prediction have been evaluated. Further, the relationship between inter-residue contacts and protein-folding rates has been highlighted. Moreover, the importance of inter-residue interactions in protein-folding kinetics and for understanding the stability of proteins has been discussed. In essence, the information gained from the studies on inter-residue interactions provides valuable insights for understanding protein folding and de novo protein design.

Evolutionarily conserved regions and hydrophobic contacts at the superfamily level: The case of the fold-type I, pyridoxal-5'-phosphate-dependent enzymes

The wealth of biological information provided by structural and genomic projects opens new prospects of understanding life and evolution at the molecular level. In this work, it is shown how computational approaches can be exploited to pinpoint protein structural features that remain invariant upon long evolutionary periods in the fold-type I, PLP-dependent enzymes. A nonredundant set of 23 superposed crystallographic structures belonging to this superfamily was built. Members of this family typically display high-structural conservation despite low-sequence identity. For each structure, a multiple-sequence alignment of orthologous sequences was obtained, and the 23 alignments were merged using the structural information to obtain a comprehensive multiple alignment of 921 sequences of fold-type I enzymes. The structurally conserved regions (SCRs), the evolutionarily conserved residues, and the conserved hydrophobic contacts (CHCs) were extracted from this data set, using both sequence and structural information. The results of this study identified a structural pattern of hydrophobic contacts shared by all of the superfamily members of fold-type I enzymes and involved in native interactions. This profile highlights the presence of a nucleus for this fold, in which residues participating in the most conserved native interactions exhibit preferential evolutionary conservation, that correlates significantly (r = 0.70) with the extent of mean hydrophobic contact value of their apolar fraction.

Locating the stabilizing residues in (alpha/beta)8 barrel proteins based on hydrophobicity, long-range interactions, and sequence conservation

In nature, 1 out of every 10 proteins has an (alpha/beta)(8) (TIM)-barrel fold, and in most cases, pairwise comparisons show no sequence similarity between them. Hence, delineating the key residues that induce very different sequences to share a common fold is important for understanding the folding and stability of TIM-barrel domains. In this work, we propose a new consensus approach for locating these stabilizing residues based on long-range interactions, hydrophobicity, and conservation of amino acid residues. We have identified 957 stabilizing residues in 63 proteins from a nonredundant set of 71 TIM-barrel domains. Most of these residues are located in the 8-stranded beta-sheet, with nearly one half of them oriented toward the interior of the barrel and the other half oriented toward the surrounding alpha-helices. Several stabilizing residues are found in the N- and C-terminal loops, whereas very few appear in the alpha-helices that surround the internal beta-sheet. Further, these 957 residues are placed in 434 stabilizing segments of various sizes, and each domain contains 1-10 of these segments. We found that 8 segments per domain is the most abundant one, and two thirds of the proteins have 7-9 stabilizing segments. Finally, we verified the identified residues with experimental temperature factors and found that these residues are among the ones with less mobility in the considered proteins. We suggest that our new protocol serves as a powerful tool to identify the stabilizing residues in TIM-barrel domains, which can be used as potential candidates for studying protein folding and stability by means of protein engineering experiments.

Importance of hydrophobic cluster formation through long-range contacts in the folding transition state of two-state proteins

Understanding the folding pathways of proteins is a challenging task. The Phi value approach provides a detailed understanding of transition-state structures of folded proteins. In this work, we have computed the hydrophobicity associated with each residue in the folded state of 16 two-state proteins and compared the Phi values of each mutant residue. We found that most of the residues with high Phi value coincide with local maximum in surrounding hydrophobicity, or have nearby residues that show such maximum in hydrophobicity, indicating the importance of hydrophobic interactions in the transition state. We have tested our approach to different structural classes of proteins, such as alpha-helical, SH3 domains of all-beta proteins, beta-sandwich, and alpha/beta proteins, and we observed a good agreement with experimental results. Further, we have proposed a hydrophobic contact network pattern to relate the Phi values with long-range contacts, which will be helpful to understand the transition-state structures of folded proteins. The present approach could be used to identify potential hydrophobic clusters that may form through long-range contacts during the transition state.

Five Hierarchical Levels of Sequence-Structure Correlation in Proteins

This article reviews recent work towards modelling protein folding pathways using a bioinformatics approach. Statistical models have been developed for sequence-structure correlations in proteins at five levels of structural complexity: (i) short motifs; (ii) extended motifs; (iii) nonlocal pairs of motifs; (iv) 3-dimensional arrangements of multiple motifs; and (v) global structural homology. We review statistical models, including sequence profiles, hidden Markov models (HMMs) and interaction potentials, for the first four levels of structural detail. The I-sites (folding Initiation sites) Library models short local structure motifs. Each succeeding level has a statistical model, as follows: HMMSTR (HMM for STRucture) is an HMM for extended motifs; HMMSTR-CM (Contact Maps) is a model for pairwise interactions between motifs; and SCALI-HMM (HMMs for Structural Core ALIgnments) is a set of HMMs for the spatial arrangements of motifs. The parallels between the statistical models and theoretical models for folding pathways are discussed in this article; however, global sequence models are not discussed because they have been extensively reviewed elsewhere. The data used and algorithms presented in this article are available at http://www.bioinfo.rpi.edu/~bystrc/ (click on "servers" or "downloads") or by request to bystrc@rpi.edu .