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.     

Natively unfolded proteins

It is now clear that a significant fraction of eukaryotic genomes encode proteins with substantial regions of disordered structure. In spite of the lack of structure, these proteins nevertheless are functional; many are involved in critical steps of the cell cycle and regulatory processes. In general, intrinsically disordered proteins interact with a target ligand (often DNA) and undergo a structural transition to a folded form when bound. Several features of intrinsically disordered proteins make them well suited to interacting with multiple targets and to cell regulation. New algorithms have been developed to identify disordered regions of proteins and have demonstrated their presence in cancer-associated proteins and proteins regulated by phosphorylation.     

Structural characterization of an intrinsically unfolded mini-HBX protein from hepatitis B virus

The hepatitis B virus x protein (HBX) is expressed in HBVinfected liver cells and can interact with a wide range of cellular proteins. In order to understand such promiscuous behavior of HBX we expressed a truncated mini-HBX protein (named Tr-HBX) (residues 18-142) with 5 Cys → Ser mutations and characterized its structural features using circular dichroism (CD) spectropolarimetry, NMR spectroscopy as well as bioinformatics tools for predicting disorder in intrinsically unstructured proteins (IUPs). The secondary structural content of Tr-HBX from CD data suggests that Tr-HBX is only partially folded. The protein disorder prediction by IUPred reveals that the unstructured region encompasses its N-terminal ～30 residues of Tr-HBX. A two-dimensional (1)H-(15)N HSQC NMR spectrum exhibits fewer number of resonances than expected, suggesting that Tr-HBX is a hybrid type IUP where its folded C-terminal half coexists with a disordered N-terminal region. Many IUPs are known to be capable of having promiscuous interactions with a multitude of target proteins. Therefore the intrinsically disordered nature of Tr-HBX revealed in this study provides a partial structural basis for the promiscuous structure-function behavior of HBX.     

Short peptides as predictors for the structure of polyarginine sequences in disordered proteins

Intrinsically disordered proteins and intrinsically disordered regions are frequently enriched in charged amino acids. Intrinsically disordered regions are regularly involved in important biological processes in which one or more charged residues is the driving force behind a protein-biomolecule interaction. Several lines of experimental and computational evidence suggest that polypeptides and proteins that carry high net charges have a high preference for extended conformations with average end-to-end distances exceeding expectations for self-avoiding random coils. Here, we show that charged arginine residues even in short glycine-capped model peptides (GRRG and GRRRG) significantly affect the conformational propensities of each other when compared with the intrinsic propensities of a mostly unperturbed arginine in the tripeptide GRG. A conformational analysis based on experimentally determined J-coupling constants from heteronuclear NMR spectroscopy and amide I′ band profiles from vibrational spectroscopy reveals that nearest-neighbor interactions stabilize extended β-strand conformations at the expense of polyproline II and turn conformations. The results from molecular dynamics simulations with a CHARMM36m force field and TIP3P water reproduce our results only to a limited extent. The use of the Ramachandran distribution of the central residue of GRRRG in a calculation of end-to-end distances of polyarginines of different length yielded the expected power law behavior. The scaling coefficient of 0.66 suggests that such peptides would be more extended than predicted by a self-avoiding random walk. Our findings thus support in principle theoretical predictions.     

Influence of electrostatic forces on the association kinetics and conformational ensemble of an intrinsically disordered protein

Recent work has revealed that the association of a disordered region of a protein with a folded binding partner can occur as rapidly as association between two folded proteins. This is the case for the phosphatase calcineurin (CaN) and its association with its activator calmodulin. Calmodulin binds to the intrinsically disordered regulatory domain of CaN. Previous studies have shown that electrostatic steering can accelerate the binding of folded proteins with disordered ligands. Given that electrostatic forces are strong determinants of disordered protein ensembles, the relationship between electrostatics, conformational ensembles, and quaternary interactions is unclear. Here, we employ experimental approaches to explore the impact of electrostatic interactions on the association of calmodulin with the disordered regulatory region of CaN. We find that estimated association rate constants of calmodulin with our chosen calmodulin-substrates are within the diffusion-limited regime. The association rates are dependent on the ionic strength, indicating that favorable electrostatic forces increase the rate of association. Further, we show that charged amino acids outside the calmodulin-binding site modulate the binding rate. Conformational ensembles obtained from computer simulations suggest that electrostatic interactions within the regulatory domain might bias the conformational ensemble such that the calmodulin binding region is readily accessible. Given the prevalence of charged residues in disordered protein chains, our findings are likely relevant to many protein-protein interactions.     

Intrinsically disordered protein from a pathogenic mesophile Mycobacterium tuberculosis adopts structured conformation at high temperature

Compared to eukaryotes, the occurrence of "intrinsically disordered" or "natively unfolded" proteins in prokaryotes has not been explored extensively. Here, we report the occurrence of an intrinsically disordered protein from the mesophilic human pathogen Mycobacterium tuberculosis. The Histidine-tagged recombinant Rv3221c biotin-binding protein is intrinsically disordered at ambient and physiological growth temperatures as revealed by circular dichroism and Fourier transform infrared (FTIR) spectroscopic studies. However, an increase in temperature induces a transition from disordered to structured state with a folding temperature of approximately 53 degrees C. Addition of a structure inducing solvent trifluoroethanol (TFE) causes the protein to fold at lower temperatures suggesting that TFE fosters hydrophobic interactions, which drives protein folding. Differential Scanning Calorimetry studies revealed that folding is endothermic and the transition from a disordered to structured state is continuous (higher-order), implying existence of intermediates during folding process. Secondary structure analysis revealed that the protein has propensity to form beta-sheets. This is in conformity with FTIR spectrum that showed an absorption peak at wave number of 1636 cm(-1), indicative of disordered beta-sheet conformation in the native state. These data suggest that although Rv3221c may be disordered under ambient or optimal growth temperature conditions, it has the potential to fold into ordered structure at high temperature driven by increased hydrophobic interactions. In contrast to the generally known behavior of other intrinsically disordered proteins folding at high temperature, Rv3221c does not appear to oligomerize or aggregate as revealed through numerous experiments including Congo red binding, Thioflavin T-binding, turbidity measurements, and examining molar ellipticity as a function of protein concentration. The amino acid composition of Rv3221c reveals that it has 24% charged and 54.9% hydrophobic amino acid residues. In this respect, this protein, although belonging to the class of intrinsically disordered proteins, has distinct features. The intrinsically disordered state and the biotin-binding feature of this protein suggest that it may participate in many biochemical processes requiring biotin as a cofactor and adopt suitable conformations upon binding other folded targets.     

Structured proteins and proteins with intrinsic disorder

Until recently, the point of view that the unique tertiary structure is necessary for protein function has prevailed. However, recent data have demonstrated that many cell proteins do not possess such structure in isolation, although displaying a distinct function under physiological conditions. These proteins were named the naturally, or intrinsically, disordered proteins. The fraction of intrinsically disordered regions in such proteins may vary from several amino acid residues to a completely unordered sequence of several tens or even several hundreds of residues. The main distinction of these proteins from structured (globular) proteins is that they have no unique tertiary structure in isolation and acquire it only upon interaction with their partners. The conformation of these proteins in a complex is determined not only by their own amino acid sequence (as is typical of structured, or globular, proteins) but also by the interacting partner. This review discusses the structure-function relationships in structured and intrinsically disordered proteins. The intricateness of this problem and the possible ways to solve it are illustrated by the example of the EF1A elongation factor family.     

Identification of transient hub proteins and the possible structural basis for their multiple interactions

Proteins that can interact with multiple partners play central roles in the network of protein-protein interactions. They are called hub proteins, and recently it was suggested that an abundance of intrinsically disordered regions on their surfaces facilitates their binding to multiple partners. However, in those studies, the hub proteins were identified as proteins with multiple partners, regardless of whether the interactions were transient or permanent. As a result, a certain number of hub proteins are subunits of stable multi-subunit proteins, such as supramolecules. It is well known that stable complexes and transient complexes have different structural features, and thus the statistics based on the current definition of hub proteins will hide the true nature of hub proteins. Therefore, in this paper, we first describe a new approach to identify proteins with multiple partners dynamically, using the Protein Data Bank, and then we performed statistical analyses of the structural features of these proteins. We refer to the proteins as transient hub proteins or sociable proteins, to clarify the difference with hub proteins. As a result, we found that the main difference between sociable and nonsociable proteins is not the abundance of disordered regions, in contrast to the previous studies, but rather the structural flexibility of the entire protein. We also found greater predominance of charged and polar residues in sociable proteins than previously reported.     

An electrostatic basis for the stability of thermophilic proteins

Two factors provide key contributions to the stability of thermophilic proteins relative to their mesophilic homologues: electrostatic interactions of charged residues in the folded state and the dielectric response of the folded protein. The dielectric response for proteins in a "thermophilic series" globally modulates the thermal stability of its members, with the calculated dielectric constant for the protein increasing from mesophiles to hyperthermophiles. This variability results from differences in the distribution of charged residues on the surface of the protein, in agreement with structural and genetic observations. Furthermore, the contribution of electrostatic interactions to the stability of the folded state is more favorable for thermophilic proteins than for their mesophilic homologues. This leads to the conclusion that electrostatic interactions play an important role in determining the stability of proteins at high temperatures. The interplay between electrostatic interactions and dielectric response also provides further rationalization for the enhanced stability of thermophilic proteins with respect to cold-denaturation. Taken together, the distribution of charged residues and their fluctuations have been shown to be factors in modulating protein stability over the entire range of biologically relevant temperatures.     

Hydrophobicity scale for proteins based on inverse temperature transitions.

In general, proteins fold with hydrophobic residues buried, away from water. Reversible protein folding due to hydrophobic interactions results from inverse temperature transitions where folding occurs on raising the temperature. Because homoiothermic animals constitute an infinite heat reservoir, it is the transition temperature, Tt, not the endothermic heat of the transition, that determines the hydrophobically folded state of polypeptides at body temperature. Reported here is a new hydrophobicity scale based on the values of Tt for each amino acid residue as a guest in a natural repeating peptide sequence, the high polymers of which exhibit reversible inverse temperature transitions. Significantly, a number of ways have been demonstrated for changing Tt such that reversibly lowering Tt from above to below physiological temperature becomes a means of isothermally and reversibly driving hydrophobic folding. Accordingly, controlling Tt becomes a mechanism whereby proteins can be induced to carry out isothermal free energy transduction.     

Introducing intrinsic disorder reduces electrostatic steering in protein-protein interactions

Protein-protein interactions underlie many critical biology functions, such as cellular signaling and gene expression, in which electrostatic interactions can play a critical role in mediating the specificity and stability of protein complexes. A substantial portion of proteins are intrinsically disordered, and the influences of structural disorder on binding kinetics and thermodynamics have been widely investigated. However, whether the effect of electrostatic steering depends on structural disorder remains unexplored. In this work, we addressed the consequence of introducing intrinsic disorder in the electrostatic steering of the E3/Im3 complex using molecular dynamics simulation. Our results recapitulated the experimental observations that the responses of stability and kinetics to salt concentration for the ordered E3/Im3 complex were larger than those for the disordered E3/Im3 complex. Mechanistic analysis revealed that the native contact interactions involved in the encounter state and the transition state were essentially identical for both ordered and disordered E3. Therefore, the observed difference in electrostatic steering between ordered E3 and disordered E3 may result from their difference in conformation rather than their difference in binding mechanism. Because charged residues are frequently involved in protein-protein interactions, our results suggest that increasing structural disorder is expected to generally modulate the effect of electrostatic steering.     

Sequence characteristics responsible for protein-protein interactions in the intrinsically disordered regions of caseins,amelogenins, and small heat-shock proteins

Milk caseins and dental amelogenins are intrinsically disordered proteins (IDPs) that associate with themselves and others. Paradoxically, they are also described as hydrophobic proteins, which is difficult to reconcile with a solvent-exposed conformation. We attempt to resolve this paradox. We show that caseins and amelogenins are not hydrophobic proteins but they are more hydrophobic than most IDPs. Remarkably, uncharged residues from different regions of these mature proteins have a nearly constant average hydropathy but these regions exhibit different charged residue frequencies. A novel sequence analysis method was developed to identify hydrophobic and order-promoting regions that would favor conformational collapse. We found that such regions were uncommon; most hydrophobic and order-promoting residues were adjacent to hydrophilic or disorder-promoting residues. A further reason why caseins and amelogenins do not collapse is their high proportion of disorder-promoting proline residues. We conclude that in these proteins the hydrophobic effect is not large enough to cause conformational collapse but it can contribute, along with polar interactions, to protein-protein interactions. This behaviour is similar to the interaction of the disordered N-terminal region of small heat-shock proteins with either themselves during oligomer formation or other, unfolding, proteins during chaperone action.     

More than just tails: intrinsic disorder in histone proteins

Many biologically active proteins are disordered as a whole, or contain long disordered regions. These intrinsically disordered proteins/regions are very common in nature, abundantly found in all organisms, where they carry out important biological functions. The functions of these proteins complement the functional repertoire of "normal" ordered proteins, and many protein functional classes are heavily dependent on intrinsic disorder. Among these disorder-centric functions are interactions with nucleic acids and protein complex assembly. In this study, we present the results of comprehensive bioinformatics analyses of the abundance and roles of intrinsic disorder in 2007 histones from 746 species. We show that all the members of the histone family are intrinsically disordered proteins. Furthermore, intrinsic disorder is not only abundant in histones, but is absolutely necessary for various histone functions, starting from heterodimerization to formation of higher order oligomers, to interactions with DNA and other proteins, and to posttranslational modifications.     

The origin and impact of bound water around intrinsically disordered proteins

Proteins and water couple dynamically over a wide range of time scales. Motivated by their central role in protein function, protein-water dynamics and thermodynamics have been extensively studied for structured proteins, where correspondence to structural features has been made. However, properties controlling intrinsically disordered protein (IDP)-water dynamics are not yet known. We report results of megahertz-to-terahertz dielectric spectroscopy and molecular dynamics simulations of a group of IDPs with varying charge content along with structured proteins of similar size. Hydration water around IDPs is found to exhibit more heterogeneous rotational and translational dynamics compared with water around structured proteins of similar size, yielding on average more restricted dynamics around individual residues of IDPs, charged or neutral, compared with structured proteins. The on-average slower water dynamics is found to arise from excess tightly bound water in the first hydration layer, which is related to greater exposure to charged groups. The more tightly bound water to IDPs correlates with the smaller hydration shell found experimentally, and affects entropy associated with protein-water interactions, the contribution of which we estimate based on the dielectric measurements and simulations. Water-IDP dynamic coupling at terahertz frequencies is characterized by the dielectric measurements and simulations.     

The pairwise energy content estimated from amino acid composition discriminates between folded and intrinsically unstructured proteins

The structural stability of a protein requires a large number of interresidue interactions. The energetic contribution of these can be approximated by low-resolution force fields extracted from known structures, based on observed amino acid pairing frequencies. The summation of such energies, however, cannot be carried out for proteins whose structure is not known or for intrinsically unstructured proteins. To overcome these limitations, we present a novel method for estimating the total pairwise interaction energy, based on a quadratic form in the amino acid composition of the protein. This approach is validated by the good correlation of the estimated and actual energies of proteins of known structure and by a clear separation of folded and disordered proteins in the energy space it defines. As the novel algorithm has not been trained on unstructured proteins, it substantiates the concept of protein disorder, i.e. that the inability to form a well-defined 3D structure is an intrinsic property of many proteins and protein domains. This property is encoded in their sequence, because their biased amino acid composition does not allow sufficient stabilizing interactions to form. By limiting the calculation to a predefined sequential neighborhood, the algorithm was turned into a position-specific scoring scheme that characterizes the tendency of a given amino acid to fall into an ordered or disordered region. This application we term IUPred and compare its performance with three generally accepted predictors, PONDR VL3H, DISOPRED2 and GlobPlot on a database of disordered proteins.     

Nucleoporins’ exclusive amino acid sequence features regulate their transient interaction with and selectivity of cargo complexes in the nuclear pore

Nucleocytoplasmic traffic of nucleic acids and proteins across the nuclear envelop via the nuclear pore complexes (NPCs) is vital for eukaryotic cells. NPCs screen transported macromolecules based on their morphology and surface chemistry. This selective nature of the NPC-mediated traffic is essential for regulating the fundamental functions of the nucleus, such as gene regulation, protein synthesis, and mechanotransduction. Despite the fundamental role of the NPC in cell and nuclear biology, the detailed mechanisms underlying how the NPC works have remained largely unknown. The critical components of NPCs enabling their selective barrier function are the natively unfolded phenylalanine- and glycine-rich proteins called “FG-nucleoporins” (FG Nups). These intrinsically disordered proteins are tethered to the inner wall of the NPC, and together form a highly dynamic polymeric meshwork whose physicochemical conformation has been the subject of intense debate. We observed that specific sequence features (called largest positive like-charge regions, or lpLCRs), characterized by extended subsequences that only possess positively charged amino acids, significantly affect the conformation of FG Nups inside the NPC. Here we investigate how the presence of lpLCRs affects the interactions between FG Nups and their interactions with the cargo complex. We combine coarse-grained molecular dynamics simulations with time-resolved force distribution analysis to disordered proteins to explore the behavior of the system. Our results suggest that the number of charged residues in the lpLCR domain directly governs the average distance between Phe residues and the intensity of interaction between them. As a result, the number of charged residues within lpLCR determines the balance between the hydrophobic interaction and the electrostatic repulsion and governs how dense and disordered the hydrophobic network formed by FG Nups is. Moreover, changing the number of charged residues in an lpLCR domain can interfere with ultrafast and transient interactions between FG Nups and the cargo complex.     

Network representation of protein interactions: Theory of graph description and analysis

A methodological framework is presented for the graph theoretical interpretation of NMR data of protein interactions. The proposed analysis generalizes the idea of network representations of protein structures by expanding it to protein interactions. This approach is based on regularization of residue‐resolved NMR relaxation times and chemical shift data and subsequent construction of an adjacency matrix that represents the underlying protein interaction as a graph or network. The network nodes represent protein residues. Two nodes are connected if two residues are functionally correlated during the protein interaction event. The analysis of the resulting network enables the quantification of the importance of each amino acid of a protein for its interactions. Furthermore, the determination of the pattern of correlations between residues yields insights into the functional architecture of an interaction. This is of special interest for intrinsically disordered proteins, since the structural (three‐dimensional) architecture of these proteins and their complexes is difficult to determine. The power of the proposed methodology is demonstrated at the example of the interaction between the intrinsically disordered protein osteopontin and its natural ligand heparin.     

Purification and reconstitution of the connexin43 carboxyl terminus attached to the 4th transmembrane domain in detergent micelles

In recent years, reports have identified that many eukaryotic proteins contain disordered regions spanning greater than 30 consecutive residues in length. In particular, a number of these intrinsically disordered regions occur in the cytoplasmic segments of plasma membrane proteins. These intrinsically disordered regions play important roles in cell signaling events, as they are sites for protein–protein interactions and phosphorylation. Unfortunately, in many crystallographic studies of membrane proteins, these domains are removed because they hinder the crystallization process. Therefore, a purification procedure was developed to enable the biophysical and structural characterization of these intrinsically disordered regions while still associated with the lipid environment. The carboxyl terminal domain from the gap junction protein connexin43 attached to the 4th transmembrane domain (TM4-Cx43CT) was used as a model system (residues G178-I382). The purification was optimized for structural analysis by nuclear magnetic resonance (NMR) because this method is well suited for small membrane proteins and proteins that lack a well-structured three-dimensional fold. The TM4-Cx43CT was purified to homogeneity with a yield of 6 mg/L from C41(DE3) bacterial cells, reconstituted in the anionic detergent 1-palmitoyl-2-hydroxy-sn-glycero-3-[phospho-RAC-(1-glycerol)], and analyzed by circular dichroism and NMR to demonstrate that the TM4-Cx43CT was properly folded into a functional conformation by its ability to form α-helical structure and associate with a known binding partner, the c-Src SH3 domain, respectively.     

Conformational propensities and residual structures in unfolded peptides and proteins

Ample evidence gathered over the last ten years indicates that unfolded and naturally disordered proteins and peptides can show local order in that short segments can adopt turn or polyproline II-like conformations. These findings show that unfolded states cannot be described by the so-called random coil model which assumes that individual amino acid residues sample the entire sterically accessible parts of the Ramachandran with very similar probabilities. This article reviews the experimental evidence for the notion that amino acid residues have different propensities for polyproline II, β-strand, helical and turn conformations in water. These propensities are changed by interactions with nearest neighbours. We show that for a substantial number of residues the conformational propensities in the unfolded state correlate with values for helix propagation and the Chou-Fasman propensities for β-strands. Based on the presented results we hypothesize that the conformational distributions of a representative set of short peptides could be used for predicting structural distributions of disordered peptides and proteins in the future.