The role of disorder in interaction networks: a structural analysis

Recent studies have emphasized the value of including structural information into the topological analysis of protein networks. Here, we utilized structural information to investigate the role of intrinsic disorder in these networks. Hub proteins tend to be more disordered than other proteins (i.e. the proteome average); however, we find this only true for those with one or two binding interfaces (‘single’‐interface hubs). In contrast, the distribution of disordered residues in multi‐interface hubs is indistinguishable from the overall proteome. Surprisingly, we find that the binding interfaces in single‐interface hubs are highly structured, as is the case for multi‐interface hubs. However, the binding partners of single‐interface hubs tend to have a higher level of disorder than the proteome average, suggesting that their binding promiscuity is related to the disorder of their binding partners. In turn, the higher level of disorder of single‐interface hubs can be partly explained by their tendency to bind to each other in a cascade. A good illustration of this trend can be found in signaling pathways and, more specifically, in kinase cascades. Finally, our findings have implications for the current controversy related to party and date‐hubs.     

Intrinsic disorder of the extracellular matrix

The extracellular matrix is very well organized at the supramolecular and tissue levels and little is known on the potential role of intrinsic disorder in promoting its organization. We predicted the amount of disorder and identified disordered regions in the human extracellular proteome with established computational tools. The extracellular proteome is significantly enriched in proteins comprising more than 50% of disorder compared to the complete human proteome. The enrichment is mostly due to long disordered regions containing at least 100 consecutive disordered residues. The amount of intrinsic disorder is heterogeneous in the extracellular protein families, with the most disordered being collagens and the small integrin-binding ligand N-linked glycoproteins. Although most domains found in extracellular proteins are structured, the fibronectin III domains contain a variable amount of disordered residues (up to 92%). Binding sites for heparin and integrins are found in disordered sequences of extracellular proteins. Intrinsic disorder is evenly distributed in hubs and ends in the interaction network of extracellular proteins with their extracellular partners. In contrast, extracellular hubs are significantly enriched in disorder in the network of extracellular proteins with their extracellular, membrane and intracellular partners. Disorder could thus provide the structural plasticity required for the hubs to interact with membrane and intracellular proteins. Organization and assembly of the extracellular matrix, development of mineralized tissues and cell-matrix adhesion are the biological processes overrepresented in the most disordered extracellular proteins. Extracellular disorder is associated with binding to growth factors, glycosaminoglycans and integrins at the molecular level.     

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.     

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.     

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.     

Adsorption equilibrium in immunoaffnity chromatography with antibodies to synthetic peptides

The effects of charged residues in peptide antigens on the binding characteristics of polyclonal antipeptide antibodies were studied using immunoadsorbents prepared by coupling the antibodies to CNBr-activated Sepharose 4B. Among the antipeptide antibodies, an antibody to the peptide without charged residues showed the most stable interaction with the peptide to the changes in pH. Conversely, the binding affinity of antibodies to the pep-tides with histidine residues having a unique pKa value of 6.0 decreased steeply with pH at around 6.0. The binding affinity of an antibody to the peptide with many charged residues decreased steeply with an increase in the ionic strength (adjusted by NaCl). Since circular dichroism (CD) spectrum measurements indicate that these peptides show disordered structures in the pH range of adsorption measurement, the dependence of peptide-antibody interaction on environmental conditions is attributed to the characteristics of side chains of the peptides. These results indicate that the dependence of the binding affinity of antipeptide antibodies on pH and the ionic strength is dominantly affected by the number and the pKa values of charged residues in the peptides.     

A dissection of the protein-protein interfaces in icosahedral virus capsids

We selected 49 icosahedral virus capsids whose crystal structures are reported in the Protein Data Bank. They belong to the T=1, T=3, pseudo T=3 and other lattice types. We identified in them 779 unique interfaces between pairs of subunits, all repeated by icosahedral symmetry. We analyzed the geometric and physical chemical properties of these interfaces and compared with interfaces in protein-protein complexes and homodimeric proteins, and with crystal packing contacts. The capsids contain one to 16 subunits implicated in three to 66 unique interfaces. Each subunit loses 40-60% of its accessible surface in contacts with an average of 8.5 neighbors. Many of the interfaces are very large with a buried surface area (BSA) that can exceed 10,000 A(2), yet 39% are small with a BSA<800 A(2) comparable to crystal packing contacts. Pairwise capsid interfaces overlap, so that one-third of the residues are part of more than one interface. Those with a BSA>800 A(2) resemble homodimer interfaces in their chemical composition. Relative to the protein surface, they are non-polar, enriched in aliphatic residues and depleted of charged residues, but not of neutral polar residues. They contain one H-bond per about 200 A(2) BSA. Small capsid interfaces (BSA<800 A(2)) are only slightly more polar. They have a similar amino acid composition, but they bury fewer atoms and contain fewer H-bonds for their size. Geometric parameters that estimate the quality of the atomic packing suggest that the small capsid interfaces are loosely packed like crystal packing contacts, whereas the larger interfaces are close-packed as in protein-protein complexes and homodimers. We discuss implications of these findings on the mechanism of capsid assembly, assuming that the larger interfaces form first to yield stable oligomeric species (capsomeres), and that medium-size interfaces allow the stepwise addition of capsomeres to build larger intermediates.     

BATMAS30: amino acid substitution matrix for alignment of bacterial transporters

Aligned amino acid sequences of three functionally independent samples of transmembrane (TM) transport proteins have been analyzed. The concept of TM-kernel is proposed as the most probable transmembrane region of a sequence. The average amino acid composition of TM-kernels differs from the published amino acid composition of transmembrane segments. TM-kernels contain more alanines, glycines, and less polar, charged, and aromatic residues in contrast to non-TM-proteins. There are also differences between TM-kernels of bacterial and eukaryotic proteins. We have constructed amino acid substitution matrices for bacterial TM-kernels, named the BATMAS (BActerial Transmembrane MAtrix of Substitutions) series. In TM-kernels, polar and charged residues, as well as proline and tyrosine, are highly conserved, whereas there are more substitutions within the group of hydrophobic residues, in contrast to non-TM-proteins that have fewer, relatively more conserved, hydrophobic residues. These results demonstrate that alignment of transmembrane proteins should be based on at least two amino acid substitution matrices, one for loops (e.g., the BLOSUM series) and one for TM-segments (the BATMAS series), and the choice of the TM-matrix should be different for eukaryotic and bacterial proteins.     

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.     

Phosphorylation Variation during the Cell Cycle Scales with Structural Propensities of Proteins

Phosphorylation at specific residues can activate a protein, lead to its localization to particular compartments, be a trigger for protein degradation and fulfill many other biological functions. Protein phosphorylation is increasingly being studied at a large scale and in a quantitative manner that includes a temporal dimension. By contrast, structural properties of identified phosphorylation sites have so far been investigated in a static, non-quantitative way. Here we combine for the first time dynamic properties of the phosphoproteome with protein structural features. At six time points of the cell division cycle we investigate how the variation of the amount of phosphorylation correlates with the protein structure in the vicinity of the modified site. We find two distinct phosphorylation site groups: intrinsically disordered regions tend to contain sites with dynamically varying levels, whereas regions with predominantly regular secondary structures retain more constant phosphorylation levels. The two groups show preferences for different amino acids in their kinase recognition motifs - proline and other disorder-associated residues are enriched in the former group and charged residues in the latter. Furthermore, these preferences scale with the degree of disorderedness, from regular to irregular and to disordered structures. Our results suggest that the structural organization of the region in which a phosphorylation site resides may serve as an additional control mechanism. They also imply that phosphorylation sites are associated with different time scales that serve different functional needs.     

Solution study of engineered quartz binding peptides using replica exchange molecular dynamics

We use replica-exchange molecular dynamics (REMD) to interrogate molecular structures and properties of four engineered dodecapeptides (in solution, in the absence of a surface) that have been shown to bind to quartz with different propensities. We find that all of the strong-binding peptides feature some polyproline type II secondary structure, have less conformational freedom, and feature fewer intrapeptide hydrogen bonds compared with the weak binder. The regions of contiguous proline content in a given sequence appear to play a role in fostering some of these properties of the strong binders. For preliminary insights into quartz binding, we perform lattice-matching studies between a grid corresponding with the quartz (100) surface and the strong-binding peptide REMD structures. Our findings indicate a commonality among the putative contact residues, even for peptide structures with very different backbone conformations. Furthermore, interpeptide interactions in solution are studied. Our preliminary findings indicate that the strong-binder interpeptide contacts are dominated by weak, nonspecific hydrophobic interactions, while the weak-binding peptide shows more variable behavior due to the distribution of charged residues. In summary, the solution structures of peptides appear to be significant. We propose that these differences in their intra- and interpeptide interactions can influence their propensity to bind onto a solid substrate.     

Non-covalent interactions across subunit interfaces in Sm proteins

The distinguishing property of Sm protein associations is their high stability. In order to understand this property, we analyzed the interface non-covalent interactions and compared the properties of the Sm protein interfaces with those of a test set, Binding Interface Database (BID). The comparison revealed that the main differences between interfaces of Sm proteins and those of the BID set are the content of charged residues, hydrogen bonds, salt bridges, and conservation scores of interface residues. In Sm proteins, the interfaces have more hydrophobic and fewer charged residues than the surface, which is also the case for the BID test set and other proteins. However, in the interfaces, the content of charged residues in Sm proteins (26%) is substantially larger than that in the BID set (22%). Both interfaces of Sm proteins and of test set have a similar number of hydrophobic interactions per 100 Å². The interfaces of Sm proteins have substantially more hydrogen bonds than the interfaces in test set. The results show clearly that the interfaces of Sm proteins form more salt bridges compared with test set. On average, there are about 16 salt bridges per interface. The high conservation score of amino acids that are involved in non-covalent interactions in protein interfaces is an additional strong argument for their importance. The overriding conclusion from this study is that the non-covalent interactions in Sm protein interfaces considerably contribute to stability of higher order structures.     

Interior and surface of monomeric proteins

The solvent-accessible surface area (As) of 46 monomeric proteins is calculated using atomic co-ordinates from high-resolution and well-refined crystal structures. The As of these proteins can be determined to within 1 to 2% and that of their individual residues to within 10 to 20%. The As values of proteins are correlated with their molecular weight (Mr) in the range 4000 to 35,000: the power law As = 6.3 M0.73 predicts protein As values to within 4% on average. The average water-accessible surface is found to be 57% non-polar, 24% polar and 19% charged, with 5% root-mean-square variations. The molecular surface buried inside the protein is 58% non-polar, 39% polar and 4% charged. The buried surface contains more uncharged polar groups (mostly peptides) than the surface that remains accessible, but many fewer charged groups. On average, 15% of residues in small proteins and 32% in larger ones may be classed as "buried residues", having less than 5% of their surface accessible to the solvent. The accessibilities of most other residues are evenly distributed in the range 5 to 50%. Although the fraction of buried residues increases with molecular weight, the amino acid compositions of the protein interior and surface show no systematic variation with molecular weight, except for small proteins that are often very rich in buried cysteines. From amino acid compositions of protein surfaces and interiors we calculate an effective coefficient of partition for each type of residue, and derive an implied set of transfer free energy values. This is compared with other sets of partition coefficients derived directly from experimental data. The extent to which groups of residues (charged, polar and non-polar) are buried within proteins correlates well with their hydrophobicity derived from amino acid transfer experiments. Within these three groups, the correlation is low.     

Structure of three Humanin peptides with different activities upon interaction with liposome

We have recently shown that a 24 amino acid Humanin (HN) adopts an anti-parallel β-sheet structure in the presence of a negatively charged 1,2-dioleoyl-sn-glycero-3-phosphoglycerol (DOPG) and suggested a possibility that it interacts with lipid membranes and thereby exerts neuroprotective effects through the target cell surface receptors or the intracellular signaling molecules following membrane interaction events. The structures of two HN analogs, having either a S7A mutation or a S14G mutation, were examined under the identical conditions, as the S7A analog is inactive and the S14G analog is 1000-fold more active than the wild type HN. These analogs showed a secondary structure indistinguishable from the structure of HN in the presence of DOPG liposome, while unrelated peptides were disordered with and without DOPG. It thus appeared that HN and the analogs, regardless of the biological activities, have an ability to interact with DOPG liposome and form an anti-parallel β-sheet structure. While the wild type HN and the S7A and S14G analogs were largely disordered in buffer, the S14G analog showed greater stability as a disordered structure in the buffer at a physiological temperature, suggesting that it maintains the disordered structure presumably required for the interaction with the DOPG liposome and thereby greater neuroprotective activity.     

The roles of contact residue disorder and domain composition in characterizing protein-ligand binding specificity and promiscuity

Most protein chains interact with only one ligand but a small number of protein chains can bind several ligands, and many examples are available in the protein-ligand complex database of PDB. Among these proteins, some show preferences for the ligands or types of ligands they bind; however, so far we have only poor understanding of what determines protein-ligand binding and its specificity. Here we investigate the structural and functional properties of proteins in protein-ligand complexes. Analysis of the protein-ligand complex dataset from the PDB structure database reveals that proteins with more interactions have more disordered contact residues. Those proteins containing few disordered contact residues that bind multiple ligands have a tendency to consist of several domains. Analysis of physicochemical properties of hub contact residues binding multiple ligands indicates that they are enriched for hydrophilic, charged, polar and His-Asp catalytic triad residues. Finally, in order to differentiate proteins binding different classes of ligands, we mapped the three most prominent classes of ligands onto different superfamily domains. Our results demonstrate that contact residue disorder and ordered multiple domains are complementary factors that play a crucial role in determining ligand binding specificity and promiscuity.     

Systematic Mutational Analysis of the Yeast Act1 Gene

We report the isolation and characterization of a synoptic set of site-directed mutations distributed throughout the single actin gene of Saccharomyces cerevisiae. Mutations were systematically targeted to the surface of the protein by identifying clusters of 2 or more charged residues in the primary sequence; every charged residue in a cluster was replaced with alanine. Mutations were recovered in high yield (34 of 36 constructed) as heterozygous diploids. Mutant phenotypes were examined in haploid segregants: 11 were recessive lethal, 16 conditional-lethal (including temperature-sensitive and salt-sensitive) and 7 had no discernible phenotype. Genetic analysis suggested that the two mutations constructed but not recovered in yeast may have a dominant defective phenotype. Location of the mutant residues on the three-dimensional structure of the rabbit muscle actin monomer confirmed that most (81%) of the charged residues we altered lie at or near the surface of the protein, confirming a key assumption of the method. Many of the new act1 alleles have properties readily interpreted in light of the actin structure and should prove useful in both genetic and biochemical studies of actin function.     

X-ray structures of the maltose-maltodextrin-binding protein of the thermoacidophilic bacterium Alicyclobacillus acidocaldarius provide insight into acid stability of proteins

Maltose-binding proteins act as primary receptors in bacterial transport and chemotaxis systems. We report here crystal structures of the thermoacidostable maltose-binding protein from Alicyclobacillus acidocaldarius, and explore its modes of binding to maltose and maltotriose. Further, comparison with the structures of related proteins from Escherichia coli (a mesophile), and two hyperthermophiles (Pyrococcus furiosus and Thermococcus litoralis) allows an investigation of the basis of thermo- and acidostability in this family of proteins.The thermoacidophilic protein has fewer charged residues than the other three structures, which is compensated by an increase in the number of polar residues. Although the content of acidic and basic residues is approximately equal, more basic residues are exposed on its surface whereas most acidic residues are buried in the interior. As a consequence, this protein has a highly positive surface charge. Fewer salt bridges are buried than in the other MBP structures, but the number exposed on its surface does not appear to be unusual. These features appear to be correlated with the acidostability of the A. acidocaldarius protein rather than its thermostability.An analysis of cavities within the proteins shows that the extremophile proteins are more closely packed than the mesophilic one. Proline content is slightly higher in the hyperthermophiles and thermoacidophiles than in mesophiles, and this amino acid is more common at the second position of beta-turns, properties that are also probably related to thermostability. Secondary structural content does not vary greatly in the different structures, and so is not a contributing factor.