Robust protein protein interactions in crowded cellular environments

Modulation of α-synuclein aggregation amid diverse environmental perturbation

Intrinsically disordered protein α-synuclein (αS) is implicated in Parkinson's disease due to its aberrant aggregation propensity. In a bid to identify the traits of its aggregation, here we computationally simulate the multi-chain association process of αS in aqueous as well as under diverse environmental perturbations. In particular, the aggregation of αS in aqueous and varied environmental condition led to marked concentration differences within protein aggregates, resembling liquid-liquid phase separation (LLPS). Both saline and crowded settings enhanced the LLPS propensity. However, the surface tension of αS droplet responds differently to crowders (entropy-driven) and salt (enthalpy-driven). Conformational analysis reveals that the IDP chains would adopt extended conformations within aggregates and would maintain mutually perpendicular orientations to minimize inter-chain electrostatic repulsions. The droplet stability is found to stem from a diminished intra-chain interactions in the C-terminal regions of αS, fostering inter-chain residue-residue interactions. Intriguingly, a graph theory analysis identifies small-world-like networks within droplets across environmental conditions, suggesting the prevalence of a consensus interaction patterns among the chains. Together these findings suggest a delicate balance between molecular grammar and environment-dependent nuanced aggregation behavior of αS.

Diffusive intracellular interactions: On the role of protein net charge and functional adaptation

A striking feature of nucleic acids and lipid membranes is that they all carry net negative charge and so is true for the majority of intracellular proteins. It is suggested that the role of this negative charge is to assure a basal intermolecular repulsion that keeps the cytosolic content suitably 'fluid' for function. We focus in this review on the experimental, theoretical and genetic findings which serve to underpin this idea and the new questions they raise. Unlike the situation in test tubes, any functional protein-protein interaction in the cytosol is subject to competition from the densely crowded background, i.e. surrounding stickiness. At the nonspecific limit of this stickiness is the 'random' protein-protein association, maintaining profuse populations of transient and constantly interconverting complexes at physiological protein concentrations. The phenomenon is readily quantified in studies of the protein rotational diffusion, showing that the more net negatively charged a protein is the less it is retarded by clustering. It is further evident that this dynamic protein-protein interplay is under evolutionary control and finely tuned across organisms to maintain optimal physicochemical conditions for the cellular processes. The emerging picture is then that specific cellular function relies on close competition between numerous weak and strong interactions, and where all parts of the protein surfaces are involved. The outstanding challenge is now to decipher the very basics of this many-body system: how the detailed patterns of charged, polar and hydrophobic side chains not only control protein-protein interactions at close- and long-range but also the collective properties of the cellular interior as a whole.

Key aspects of the past 30 years of protein design

Proteins are the workhorse of life. They are the building infrastructure of living systems; they are the most efficient molecular machines known, and their enzymatic activity is still unmatched in versatility by any artificial system. Perhaps proteins' most remarkable feature is their modularity. The large amount of information required to specify each protein's function is analogically encoded with an alphabet of just ∼20 letters. The protein folding problem is how to encode all such information in a sequence of 20 letters. In this review, we go through the last 30 years of research to summarize the state of the art and highlight some applications related to fundamental problems of protein evolution.

Fever as an evolutionary agent to select immune complexes interfaces

We herein analyzed all available protein–protein interfaces of the immune complexes from the Protein Data Bank whose antigens belong to pathogens or cancers that are modulated by fever in mammalian hosts. We also included, for comparison, protein interfaces from immune complexes that are not significantly modulated by the fever response. We highlight the distribution of amino acids at these viral, bacterial, protozoan and cancer epitopes, and at their corresponding paratopes that belong strictly to monoclonal antibodies. We identify the “hotspots”, i.e. residues that are highly connected at such interfaces, and assess the structural, kinetic and thermodynamic parameters responsible for complex formation. We argue for an evolutionary pressure for the types of residues at these protein interfaces that may explain the role of fever as a selective force for optimizing antibody binding to antigens.

Role of non-specific interactions in the phase-separation and maturation of macromolecules

Phase separation of biomolecules could be mediated by both specific and non-specific interactions. How the interplay between non-specific and specific interactions along with polymer entropy influences phase separation is an open question. We address this question by simulating self-associating molecules as polymer chains with a short core stretch that forms the specifically interacting functional interface and longer non-core regions that participate in non-specific/promiscuous interactions. Our results show that the interplay of specific (strength, ϵsp) and non-specific interactions (strength, ϵns) could result in phase separation of polymers and its transition to solid-like aggregates (mature state). In the absence of ϵns, the polymer chains do not dwell long enough in the vicinity of each other to undergo phase separation and transition into a mature state. On the other hand, in the limit of strong ϵns, the assemblies cannot transition into the mature state and form a non-specific assembly, suggesting an optimal range of interactions favoring mature multimers. In the scenario where only a fraction (Nfrac) of the non-core regions participate in attractive interactions, we find that slight modifications to either ϵns or Nfrac can result in dramatically altered self-assembled states. Using a combination of heterogeneous and homogeneous mix of polymers, we establish how this interplay between interaction energies dictates the propensity of biomolecules to find the correct binding partner at dilute concentrations in crowded environments.

Integrating protein copy numbers with interaction networks to quantify stoichiometry in clathrin-mediated endocytosis

Proteins that drive processes like clathrin-mediated endocytosis (CME) are expressed at copy numbers within a cell and across cell types varying from hundreds (e.g. auxilin) to millions (e.g. clathrin). These variations contain important information about function, but without integration with the interaction network, they cannot capture how supply and demand for each protein depends on binding to shared and distinct partners. Here we construct the interface-resolved network of 82 proteins involved in CME and establish a metric, a stoichiometric balance ratio (SBR), that quantifies whether each protein in the network has an abundance that is sub- or super-stoichiometric dependent on the global competition for binding. We find that highly abundant proteins (like clathrin) are super-stoichiometric, but that not all super-stoichiometric proteins are highly abundant, across three cell populations (HeLa, fibroblast, and neuronal synaptosomes). Most strikingly, within all cells there is significant competition to bind shared sites on clathrin and the central AP-2 adaptor by other adaptor proteins, resulting in most being in excess supply. Our network and systematic analysis, including response to perturbations of network components, show how competition for shared binding sites results in functionally similar proteins having widely varying stoichiometries, due to variations in both abundance and their unique network of binding partners.

Prediction of Protein-Protein Binding Affinities from Unbound Protein Structures

Proteins are the workhorses of cells to carry out sophisticated and complex cellular processes. Such processes require a coordinated and regulated interactions between proteins that are both time and location specific. The strength, or binding affinity, of protein-protein interactions ranges between the micro- and the nanomolar association constant, often dictating the molecular mechanisms underlying the interaction and the longevity of the complex, i.e., transient or permanent. In consequence, there is a need to quantify the strength of protein-protein interactions for biological, biomedical, and biotechnological applications. While experimental methods are labor intensive and costly, computational ones are useful tools to predict the affinity of protein-protein interactions. In this chapter, we review the methods developed by us to address this question. We briefly present two methods to comprehend the structure of the protein complex derived by either comparative modeling or docking. Then we introduce BADOCK, a method to predict the binding energy without requiring the structure of the protein complex, thus overcoming one of the major limitations of structure-based methods for the prediction of binding affinity. BADOCK utilizes the structure of unbound proteins and the protein docking sampling space to predict protein-protein binding affinities. We present step-by-step protocols to utilize these methods, describing the inputs and potential pitfalls as well as their respective strengths and limitations.

Abundance Imparts Evolutionary Constraints of Similar Magnitude on the Buried, Surface, and Disordered Regions of Proteins

An understanding of the forces shaping protein conservation is key, both for the fundamental knowledge it represents and to allow for optimal use of evolutionary information in practical applications. Sequence conservation is typically examined at one of two levels. The first is a residue-level, where intra-protein differences are analyzed and the second is a protein-level, where inter-protein differences are studied. At a residue level, we know that solvent-accessibility is a prime determinant of conservation. By inverting this logic, we inferred that disordered regions are slightly more solvent-accessible on average than the most exposed surface residues in domains. By integrating abundance information with evolutionary data within and across proteins, we confirmed a previously reported strong surface-core association in the evolution of structured regions, but we found a comparatively weak association between disordered and structured regions. The facts that disordered and structured regions experience different structural constraints and evolve independently provide a unique setup to examine an outstanding question: why is a protein’s abundance the main determinant of its sequence conservation? Indeed, any structural or biophysical property linked to the abundance-conservation relationship should increase the relative conservation of regions concerned with that property (e.g., disordered residues with mis-interactions, domain residues with misfolding). Surprisingly, however, we found the conservation of disordered and structured regions to increase in equal proportion with abundance. This observation implies that either abundance-related constraints are structure-independent, or multiple constraints apply to different regions and perfectly balance each other.

Cellular Sticking Can Strongly Reduce Complex Binding by Speeding Dissociation

While extensive studies have been carried out to determine protein-RNA binding affinities, mechanisms, and dynamics in vitro, such studies do not take into consideration the effect of the many weak nonspecific interactions in a cell filled with potential binding partners. Here we experimentally tested the role of the cellular environment on affinity and binding dynamics between a protein and RNA in living U-2 OS cells. Our model system is the spliceosomal protein U1A and its binding partner SL2 of the U1 snRNA. The binding equilibrium was perturbed by a laser-induced temperature jump and monitored by Förster resonance energy transfer. The apparent binding affinity in live cells was reduced by up to 2 orders of magnitude compared to in vitro. The measured in-cell dissociation rate coefficients were up to 2 orders of magnitude larger, whereas no change in the measured association rate coefficient was observed. The latter is not what would be anticipated due to macromolecular crowding or nonspecific sticking of the uncomplexed U1A and SL2 in the cell. A quantitative model fits our experimental results, with the major cellular effect being that U1A and SL2 sticking to cellular components are capable of binding, just not as strongly as the free complex. This observation suggests that high binding affinities measured or designed in vitro are necessary for proper binding in vivo, where competition with many nonspecific interactions exists, especially for strongly interacting species with high charge or large hydrophobic surface areas.

The intracellular environment affects protein-protein interactions

Protein-protein interactions are essential for life but rarely thermodynamically quantified in living cells. In vitro efforts show that protein complex stability is modulated by high concentrations of cosolutes, including synthetic polymers, proteins, and cell lysates via a combination of hard-core repulsions and chemical interactions. We quantified the stability of a model protein complex, the A34F GB1 homodimer, in buffer, Escherichia coli cells and Xenopus laevis oocytes. The complex is more stable in cells than in buffer and more stable in oocytes than E. coli Studies of several variants show that increasing the negative charge on the homodimer surface increases stability in cells. These data, taken together with the fact that oocytes are less crowded than E. coli cells, lead to the conclusion that chemical interactions are more important than hard-core repulsions under physiological conditions, a conclusion also gleaned from studies of protein stability in cells. Our studies have implications for understanding how promiscuous-and specific-interactions coherently evolve for a protein to properly function in the crowded cellular environment.

Measuring the effects of macromolecular crowding on antibody function with biolayer interferometry

Biotherapeutic proteins are commonly dosed at high concentrations into the blood, which is an inherently complex, crowded solution with substantial protein content. The effects of macromolecular crowding may lead to an appreciable level of non-specific hetero-association in this physiological environment. Therefore, developing a method to characterize the diverse consequences of non-specific interactions between proteins under such non-ideal, crowded conditions, which deviate substantially from those commonly employed for in vitro characterization, is vital to achieving a more complete picture of antibody function in a biological context. In this study, we investigated non-specific interactions between human serum albumin (HSA) and two monoclonal antibodies (mAbs) by static light scattering and determined these interactions are both ionic strength-dependent and mAb-dependent. Using biolayer interferometry (BLI), we assessed the effect of HSA on antigen binding by mAbs, demonstrating that these non-specific interactions have a functional impact on mAb:antigen interactions, particularly at low ionic strength. While this effect is mitigated at physiological ionic strength, our in vitro data support the notion that HSA in the blood may lead to non-specific interactions with mAbs in vivo, with a potential impact on their interactions with antigen. Furthermore, the BLI method offers a high-throughput advantage compared to orthogonal techniques such as analytical ultracentrifugation and is amenable to a greater variety of solution conditions compared to nuclear magnetic resonance spectroscopy. Our study demonstrates that BLI is a viable technology for examining the impact of non-specific interactions on specific biologically relevant interactions, providing a direct method to assess binding events in crowded conditions.

Computational Assessment of Bacterial Protein Structures Indicates a Selection Against Aggregation

The aggregation of proteins compromises cell fitness, either because it titrates functional proteins into non-productive inclusions or because it results in the formation of toxic assemblies. Accordingly, computational proteome-wide analyses suggest that prevention of aggregation upon misfolding plays a key role in sequence evolution. Most proteins spend their lifetimes in a folded state; therefore, it is conceivable that, in addition to sequences, protein structures would have also evolved to minimize the risk of aggregation in their natural environments. By exploiting the AGGRESCAN3D structure-based approach to predict the aggregation propensity of >600 Escherichia coli proteins, we show that the structural aggregation propensity of globular proteins is connected with their abundance, length, essentiality, subcellular location and quaternary structure. These data suggest that the avoidance of protein aggregation has contributed to shape the structural properties of proteins in bacterial cells.

Self-assembly of a dimer system

In the self-assembly process which drives the formation of cellular membranes, micelles, and capsids, a collection of separated subunits spontaneously binds together to form functional and more ordered structures. In this work, we study the statistical physics of self-assembly in a simpler scenario: the formation of dimers from a system of monomers. The properties of the model allow us to frame the microstate counting as a combinatorial problem whose solution leads to an exact partition function. From the associated equilibrium conditions, we find that such dimer systems come in two types: "search-limited" and "combinatorics-limited," only the former of which has states where partial assembly can be dominated by correct contacts. Using estimates of biophysical quantities in systems of single-stranded DNA dimerization, transcription factor and DNA interactions, and protein-protein interactions, we find that all of these systems appear to be of the search-limited type, i.e., their fully correct dimerization regimes are more limited by the ability of monomers to find one another in the constituent volume than by the combinatorial disadvantage of correct dimers. We derive the parameter requirements for fully correct dimerization and find that rather than the ratio of particle number and volume (i.e., number density) being the relevant quantity, it is the product of particle diversity and volume that is constrained. Ultimately, this work contributes to an understanding of self-assembly by using the simple case of a system of dimers to analytically study the combinatorics of assembly.

Match fitness landscapes for macromolecular interaction networks: Selection for translational accuracy and rate can displace tRNA-binding interfaces of non-cognate aminoacyl-tRNA synthetases

Advances in structural biology of aminoacyl-tRNA synthetases (aaRSs) have revealed incredible diversity in how aaRSs bind their tRNA substrates. The causes of this diversity remain mysterious. We developed a new class of highly rugged fitness landscape models called match landscapes, through which genes encode the assortative interactions of their gene products through the complementarity and identifiability of their structural features. We used results from coding theory to prove bounds and equalities on fitness in match landscapes assuming additive interaction energies, macroscopic aminoacylation kinetics including proofreading, site-specific modifiers of interaction, and selection for translational accuracy in multiple, perfectly encoded site-types. Using genotypes based on extended Hamming codes we show that over a wide array of interface sizes and numbers of encoded cognate pairs, selection for translational accuracy alone is insufficient to displace the tRNA-binding interfaces of aaRSs. Yet, under combined selection for translational accuracy and rate, site-specific modifiers are selected to adaptively displace the tRNA-binding interfaces of non-cognate aaRS-tRNA pairs. We describe a remarkable correspondence between the lengths of perfect RNA (quaternary) codes and the modal sizes of small non-coding RNA families.

Shaping the topology of folding pathways in mechanical systems

Disordered mechanical systems, when strongly deformed, have complex configuration spaces with multiple stable states and pathways connecting them. The topology of such pathways determines which states are smoothly accessible from any part of configuration space. Controlling this topology would allow us to limit access to undesired states and select desired behaviors in metamaterials. Here, we show that the topology of such pathways, as captured by bifurcation diagrams, can be tuned using imperfections such as stiff hinges in elastic networks and creased thin sheets. We derive Linear Programming-like equations for designing desirable pathway topologies. These ideas are applied to eliminate the exponentially many ways of misfolding self-folding sheets by making some creases stiffer than others. Our approach allows robust folding by entire classes of external folding forces. Finally, we find that the bifurcation diagram makes pathways accessible only at specific folding speeds, enabling speed-dependent selection of different folded states.

Self-folding origami have applications for mechanical metamaterials but one of their pitfalls is that many undesirable folding modes exist. Here the authors propose an algorithm to determine which folding joints to make stiffer in order to ensure that the sheet is folded into the chosen state.

Effect of non-specific interactions on formation and stability of specific complexes

We introduce a simple model to describe the interplay between specific and non-specific interactions. We study the influence of various physical factors on the static and dynamic properties of the specific interactions of our model and show that contrary to intuitive expectations, non-specific interactions can assist in the formation of specific complexes and increase their stability. We then discuss the relevance of these results for biological systems.

Phase Transitions in Biological Systems with Many Components

Biological mixtures such as the cytosol may consist of thousands of distinct components. There is now a substantial body of evidence showing that, under physiological conditions, intracellular mixtures can phase separate into spatially distinct regions with differing compositions. In this article we present numerical evidence indicating that such spontaneous compartmentalization exploits general features of the phase diagram of a multicomponent biomolecular mixture. In particular, we show that demixed domains are likely to segregate when the variance in the intermolecular interaction strengths exceeds a well-defined threshold. Multiple distinct phases are likely to become stable under very similar conditions, which can then be tuned to achieve multiphase coexistence. As a result, only minor adjustments to the composition of the cytosol or the strengths of the intermolecular interactions are needed to regulate the formation of different domains with specific compositions, implying that phase separation is a robust mechanism for creating spatial organization. We further predict that this functionality is only weakly affected by increasing the number of components in the system. Our model therefore suggests that, for purely physico-chemical reasons, biological mixtures are naturally poised to undergo a small number of demixing phase transitions.

Sucrose modulates insulin amyloid-like fibril formation: effect on the aggregation mechanism and fibril morphology

Sucrose modifies the human insulin fibrillation pathways, affecting the fibril morphology.

Transient protein-protein interactions perturb E. coli metabolome and cause gene dosage toxicity

Gene dosage toxicity (GDT) is an important factor that determines optimal levels of protein abundances, yet its molecular underpinnings remain unknown. Here, we demonstrate that overexpression of DHFR in E. coli causes a toxic metabolic imbalance triggered by interactions with several functionally related enzymes. Though deleterious in the overexpression regime, surprisingly, these interactions are beneficial at physiological concentrations, implying their functional significance in vivo. Moreover, we found that overexpression of orthologous DHFR proteins had minimal effect on all levels of cellular organization - molecular, systems, and phenotypic, in sharp contrast to E. coli DHFR. Dramatic difference of GDT between 'E. coli's self' and 'foreign' proteins suggests the crucial role of evolutionary selection in shaping protein-protein interaction (PPI) networks at the whole proteome level. This study shows how protein overexpression perturbs a dynamic metabolon of weak yet potentially functional PPI, with consequences for the metabolic state of cells and their fitness.

Protein social behavior makes a stronger signal for partner identification than surface geometry

Cells are interactive living systems where proteins movements, interactions and regulation are substantially free from centralized management. How protein physico‐chemical and geometrical properties determine who interact with whom remains far from fully understood. We show that characterizing how a protein behaves with many potential interactors in a complete cross‐docking study leads to a sharp identification of its cellular/true/native partner(s). We define a sociability index, or S‐index, reflecting whether a protein likes or not to pair with other proteins. Formally, we propose a suitable normalization function that accounts for protein sociability and we combine it with a simple interface‐based (ranking) score to discriminate partners from non‐interactors. We show that sociability is an important factor and that the normalization permits to reach a much higher discriminative power than shape complementarity docking scores. The social effect is also observed with more sophisticated docking algorithms. Docking conformations are evaluated using experimental binding sites. These latter approximate in the best possible way binding sites predictions, which have reached high accuracy in recent years. This makes our analysis helpful for a global understanding of partner identification and for suggesting discriminating strategies. These results contradict previous findings claiming the partner identification problem being solvable solely with geometrical docking. Proteins 2016; 85:137–154. © 2016 Wiley Periodicals, Inc.

Great interactions: How binding incorrect partners can teach us about protein recognition and function

Protein–protein interactions play a key part in most biological processes and understanding their mechanism is a fundamental problem leading to numerous practical applications. The prediction of protein binding sites in particular is of paramount importance since proteins now represent a major class of therapeutic targets. Amongst others methods, docking simulations between two proteins known to interact can be a useful tool for the prediction of likely binding patches on a protein surface. From the analysis of the protein interfaces generated by a massive cross‐docking experiment using the 168 proteins of the Docking Benchmark 2.0, where all possible protein pairs, and not only experimental ones, have been docked together, we show that it is also possible to predict a protein's binding residues without having any prior knowledge regarding its potential interaction partners. Evaluating the performance of cross‐docking predictions using the area under the specificity‐sensitivity ROC curve (AUC) leads to an AUC value of 0.77 for the complete benchmark (compared to the 0.5 AUC value obtained for random predictions). Furthermore, a new clustering analysis performed on the binding patches that are scattered on the protein surface show that their distribution and growth will depend on the protein's functional group. Finally, in several cases, the binding‐site predictions resulting from the cross‐docking simulations will lead to the identification of an alternate interface, which corresponds to the interaction with a biomolecular partner that is not included in the original benchmark. Proteins 2016; 84:1408–1421. © 2016 The Authors Proteins: Structure, Function, and Bioinformatics Published by Wiley Periodicals, Inc.

A Simple Model of Protein Domain Swapping in Crowded Cellular Environments

Domain swapping in proteins is an important mechanism of functional and structural innovation. However, despite its ubiquity and importance, the physical mechanisms that lead to domain swapping are poorly understood. Here, we present a simple two-dimensional coarse-grained model of protein domain swapping in the cytoplasm. In our model, two-domain proteins partially unfold and diffuse in continuous space. Monte Carlo multiprotein simulations of the model reveal that domain swapping occurs at intermediate temperatures, whereas folded dimers and folded monomers prevail at low temperatures, and partially unfolded monomers predominate at high temperatures. We use a simplified amino acid alphabet consisting of four residue types, and find that the oligomeric state at a given temperature depends on the sequence of the protein. We also show that hinge strain between domains can promote domain swapping, consistent with experimental observations for real proteins. Domain swapping depends nonmonotonically on the protein concentration, with domain-swapped dimers occurring at intermediate concentrations and nonspecific interactions between partially unfolded proteins occurring at high concentrations. For folded proteins, we recover the result obtained in three-dimensional lattice simulations, i.e., that functional dimerization is most prevalent at intermediate temperatures and nonspecific interactions increase at low temperatures.

Relative Cosolute Size Influences the Kinetics of Protein-Protein Interactions

Protein signaling occurs in crowded intracellular environments, and while high concentrations of macromolecules are postulated to modulate protein-protein interactions, analysis of their impact at each step of the reaction pathway has not been systematically addressed. Potential cosolute-induced alterations in target association are particularly important for a signaling molecule like calmodulin (CaM), where competition among >300 targets governs which pathways are selectively activated. To explore how high concentrations of cosolutes influence CaM-target affinity and kinetics, we methodically investigated each step of the CaM-target binding mechanism under crowded or osmolyte-rich environments mimicked by ficoll-70, dextran-10, and sucrose. All cosolutes stabilized compact conformers of CaM and modulated association kinetics by affecting diffusion and rates of conformational change; however, the results showed that differently sized molecules had variable effects to enhance or impede unique steps of the association pathway. On- and off-rates were modulated by all cosolutes in a compensatory fashion, producing little change in steady-state affinity. From this work insights were gained on how high concentrations of inert crowding agents and osmolytes fit into a kinetic framework to describe protein-protein interactions relevant for cellular signaling.

Interaction between IGFBP7 and insulin: a theoretical and experimental study

Insulin-like growth factor binding protein 7 (IGFBP7) can bind to insulin with high affinity which inhibits the early steps of insulin action. Lack of recognition mechanism impairs our understanding of insulin regulation before it binds to insulin receptor. Here we combine computational simulations with experimental methods to investigate the interaction between IGFBP7 and insulin. Molecular dynamics simulations indicated that His200 and Arg198 in IGFBP7 were key residues. Verified by experimental data, the interaction remained strong in single mutation systems R198E and H200F but became weak in double mutation system R198E-H200F relative to that in wild-type IGFBP7. The results and methods in present study could be adopted in future research of discovery of drugs by disrupting protein-protein interactions in insulin signaling. Nevertheless, the accuracy, reproducibility, and costs of free-energy calculation are still problems that need to be addressed before computational methods can become standard binding prediction tools in discovery pipelines.

An inter-species protein-protein interaction network across vast evolutionary distance

In cellular systems, biophysical interactions between macromolecules underlie a complex web of functional interactions. How biophysical and functional networks are coordinated, whether all biophysical interactions correspond to functional interactions, and how such biophysical-versus-functional network coordination is shaped by evolutionary forces are all largely unanswered questions. Here, we investigate these questions using an "inter-interactome" approach. We systematically probed the yeast and human proteomes for interactions between proteins from these two species and functionally characterized the resulting inter-interactome network. After a billion years of evolutionary divergence, the yeast and human proteomes are still capable of forming a biophysical network with properties that resemble those of intra-species networks. Although substantially reduced relative to intra-species networks, the levels of functional overlap in the yeast-human inter-interactome network uncover significant remnants of co-functionality widely preserved in the two proteomes beyond human-yeast homologs. Our data support evolutionary selection against biophysical interactions between proteins with little or no co-functionality. Such non-functional interactions, however, represent a reservoir from which nascent functional interactions may arise.

Evolution of specificity in protein-protein interactions

Hub proteins are proteins that maintain promiscuous molecular recognition. Because they are reported to play essential roles in cellular control, there has been a special interest in the study of their structural and functional properties, yet the mechanisms by which they evolve to maintain functional interactions are poorly understood. By combining biophysical simulations of coarse-grained proteins and analysis of proteins-complex crystallographic structures, we seek to elucidate those mechanisms. We focus on two types of hub proteins: Multi hubs, which interact with their partners through different interfaces, and Singlish hubs, which do so through a single interface. We show that loss of structural stability is required for the evolution of protein-protein-interaction (PPI) networks, and it is more profound in Singlish hub systems. In addition, different ratios of hydrophobic to electrostatic interfacial amino acids are shown to support distinct network topologies (i.e., Singlish and Multi systems), and therefore underlie a fundamental design principle of PPI in a crowded environment. We argue that the physical nature of hydrophobic and electrostatic interactions, in particular, their favoring of either same-type interactions (hydrophobic-hydrophobic), or opposite-type interactions (negatively-positively charged) plays a key role in maintaining the network topology while allowing the protein amino acid sequence to evolve.

Is catalytic activity of chaperones a selectable trait for the emergence of heat shock response?

Although heat shock response is ubiquitous in bacterial cells, the underlying physical chemistry behind heat shock response remains poorly understood. To study the response of cell populations to heat shock we employ a physics-based ab initio model of living cells where protein biophysics (i.e., folding and protein-protein interactions in crowded cellular environments) and important aspects of proteins homeostasis are coupled with realistic population dynamics simulations. By postulating a genotype-phenotype relationship we define a cell division rate in terms of functional concentrations of proteins and protein complexes, whose Boltzmann stabilities of folding and strengths of their functional interactions are exactly evaluated from their sequence information. We compare and contrast evolutionary dynamics for two models of chaperon action. In the active model, foldase chaperones function as nonequilibrium machines to accelerate the rate of protein folding. In the passive model, holdase chaperones form reversible complexes with proteins in their misfolded conformations to maintain their solubility. We find that only cells expressing foldase chaperones are capable of genuine heat shock response to the increase in the amount of unfolded proteins at elevated temperatures. In response to heat shock, cells' limited resources are redistributed differently for active and passive models. For the active model, foldase chaperones are overexpressed at the expense of downregulation of high abundance proteins, whereas for the passive model; cells react to heat shock by downregulating their high abundance proteins, as their low abundance proteins are upregulated.

Mutation rates and evolution of multiple coding in RNA-based protocells

RNA has a myriad of biological roles in contemporary life. We use the RNA paradigm for genotype-phenotype mappings to study the evolution of multiple coding in dependence to mutation rates. We study three different one-to-many genotype-phenotype mappings which have the potential to encode the information for multiple functions on a single sequence. These three different maps are (i) cofolding, where two sequences can bind and “cofold,” (ii) suboptimal folding, where the alternative foldings within a certain range of the native state of sequences are considered, and (iii) adapter-based folding, in which protocells can evolve adapter-mediated alternative foldings. We study how protocells with a set of sequences can code for a set of predefined functional structures, while avoiding all other structures, which are considered to be misfoldings. Note that such misfolded structures are far more prevalent than functional ones. Our results highlight the flexibility of the RNA sequence to secondary structure mapping and the power of evolution to shape the genotype-phenotype mapping. We show that high fitness can be achieved even at high mutation rates. Mutation rates affect genome size, but differently depending on which folding method is used. We observe that cofolding limits the possibility to avoid misfolded structures and that adapters are always beneficial for fitness, but even more beneficial at low mutation rates. In all cases, the evolution procedure selects for molecules that can form additional structures. Our results indicate that inherent properties of RNA molecules and their interactions allow the evolution of complexity even at high mutation rates.

Molecular organization of cytochrome c2 near the binding domain of cytochrome bc1 studied by electron spin-lattice relaxation enhancement

Measurements of specific interactions between proteins are challenging. In redox systems, interactions involve surfaces near the attachment sites of cofactors engaged in interprotein electron transfer (ET). Here we analyzed binding of cytochrome c₂ to cytochrome bc₁ by measuring paramagnetic relaxation enhancement (PRE) of spin label (SL) attached to cytochrome c₂. PRE was exclusively induced by the iron atom of heme c₁ of cytochrome bc₁, which guaranteed that only the configurations with SL to heme c₁ distances up to ∼30 Å were detected. Changes in PRE were used to qualitatively and quantitatively characterize the binding. Our data suggest that at low ionic strength and under an excess of cytochrome c₂ over cytochrome bc₁, several cytochrome c₂ molecules gather near the binding domain forming a “cloud” of molecules. When the cytochrome bc₁ concentration increases, the cloud disperses to populate additional available binding domains. An increase in ionic strength weakens the attractive forces and the average distance between cytochrome c₂ and cytochrome bc₁ increases. The spatial arrangement of the protein complex at various ionic strengths is different. Above 150 mM NaCl the lifetime of the complexes becomes so short that they are undetectable. All together the results indicate that cytochrome c₂ molecules, over the range of salt concentration encompassing physiological ionic strength, do not form stable, long-lived complexes but rather constantly collide with the surface of cytochrome bc₁ and ET takes place coincidentally with one of these collisions.

Experimental evolution of protein-protein interaction networks

The modern synthesis of evolutionary theory and genetics has enabled us to discover underlying molecular mechanisms of organismal evolution. We know that in order to maximize an organism's fitness in a particular environment, individual interactions among components of protein and nucleic acid networks need to be optimized by natural selection, or sometimes through random processes, as the organism responds to changes and/or challenges in the environment. Despite the significant role of molecular networks in determining an organism's adaptation to its environment, we still do not know how such inter- and intra-molecular interactions within networks change over time and contribute to an organism's evolvability while maintaining overall network functions. One way to address this challenge is to identify connections between molecular networks and their host organisms, to manipulate these connections, and then attempt to understand how such perturbations influence molecular dynamics of the network and thus influence evolutionary paths and organismal fitness. In the present review, we discuss how integrating evolutionary history with experimental systems that combine tools drawn from molecular evolution, synthetic biology and biochemistry allow us to identify the underlying mechanisms of organismal evolution, particularly from the perspective of protein interaction networks.

Rigid-body motions of interacting proteins dominate multispecific binding of ubiquitin in a shape-dependent manner

To understand the dynamic aspects of multispecificity of ubiquitin, we studied nine ubiquitin-ligand (partner protein) complexes by normal mode analysis based on an elastic network model. The coupling between ubiquitin and ligand motions was analyzed by decomposing it into rigid-body (external) and vibrational (internal) motions of each subunit. We observed that in total the external motions in one of the subunits largely dominated the coupling. The combination of external motions of ubiquitin and the ligands showed general trends of rotations and translations. Moreover, we observed that the rotational motions of ubiquitin were correlated to the ligand orientations. We also identified ubiquitin atomic vibrations that differentiated the orientation of the ligand molecule. We observed that the extents of coupling were correlated to the shapes of the ligands, and this trend was more pronounced when the coupling involved vibrational motions of the ligand. In conclusion, an intricate interplay between internal and external motions of ubiquitin and the ligands help understand the dynamics of multispecificity, which is mostly guided by the shapes of the ligands and the complex.

Evolutionary pressure on the topology of protein interface interaction networks

The densely connected structure of protein-protein interaction (PPI) networks reflects the functional need of proteins to cooperate in cellular processes. However, PPI networks do not adequately capture the competition in protein binding. By contrast, the interface interaction network (IIN) studied here resolves the modular character of protein-protein binding and distinguishes between simultaneous and exclusive interactions that underlie both cooperation and competition. We show that the topology of the IIN is under evolutionary pressure, and we connect topological features of the IIN to specific biological functions. To reveal the forces shaping the network topology, we use a sequence-based computational model of interface binding along with network analysis. We find that the more fragmented structure of IINs, in contrast to the dense PPI networks, arises in large part from the competition between specific and nonspecific binding. The need to minimize nonspecific binding favors specific network motifs, including a minimal number of cliques (i.e., fully connected subgraphs) and many disconnected fragments. Validating the model, we find that these network characteristics are closely mirrored in the IIN of clathrin-mediated endocytosis. Features unexpected on the basis of our motif analysis are found to indicate either exceptional binding selectivity or important regulatory functions.

A one-shot germinal center model under protein structural stability constraints

The germinal center reaction is the process by which low-affinity B cells evolve into potent, immunoglobulin-secreting plasma and memory B cells. Since the recycling hypothesis was created, experimental studies have both tracked movement of a small population of B cells from the light zone into the dark zone, supporting the recycling model, and parallel to the light zone-dark zone interface, indicating a one-way trajectory. We present a novel, sequence-based ab initio model of protein stability and protein interactions. Our model contains a dark zone region of clonal expansion and somatic hypermutation and a light zone site of antigenic selection. We show not only that a one-shot model is sufficient to achieve biologically-realistic rates of affinity growth, population dynamics, and silent:non-silent mutation ratios in the complementary determining region and framework region of antibodies, but also that a stochastic recycling program with or without realistic constraints on the structural stabilities of GC antibodies cannot produce biologically-observed affinity growth, population dynamics or silent:non-silent mutation profiles. The effect of recycling erases affinity gains made by potent antibodies cycling back from the light zone and causes B cells to pool in the dark zone under high replication rates.

Analysis of genetic interaction networks shows that alternatively spliced genes are highly versatile

Alternative splicing has the potential to increase the diversity of the transcriptome and proteome. Where more than one transcript arises from a gene they are often so different that they are quite unlikely to have the same function. However, it remains unclear if alternative splicing generally leads to a gene being involved in multiple biological processes or whether it alters the function within a single process. Knowing that genetic interactions occur between functionally related genes, we have used them as a proxy for functional versatility, and have analysed the sets of genes of two well-characterised model organisms: Caenorhabditis elegans and Drosophila melanogaster. Using network analyses we find that few genes are functionally homogenous (only involved in a few functionally-related biological processes). Moreover, there are differences between alternatively spliced genes and genes with a single transcript; specifically, genes with alternatively splicing are, on average, involved in more biological processes. Finally, we suggest that factors other than specific functional classes determine whether a gene is alternatively spliced.

Counterbalance of ligand- and self-coupled motions characterizes multispecificity of ubiquitin

Date hub proteins are a type of proteins that show multispecificity in a time-dependent manner. To understand dynamic aspects of such multispecificity we studied Ubiquitin as a typical example of a date hub protein. Here we analyzed 9 biologically relevant Ubiquitin-protein (ligand) heterodimer structures by using normal mode analysis based on an elastic network model. Our result showed that the self-coupled motion of Ubiquitin in the complex, rather than its ligand-coupled motion, is similar to the motion of Ubiquitin in the unbound condition. The ligand-coupled motions are correlated to the conformational change between the unbound and bound conditions of Ubiquitin. Moreover, ligand-coupled motions favor the formation of the bound states, due to its in-phase movements of the contacting atoms at the interface. The self-coupled motions at the interface indicated loss of conformational entropy due to binding. Therefore, such motions disfavor the formation of the bound state. We observed that the ligand-coupled motions are embedded in the motions of unbound Ubiquitin. In conclusion, multispecificity of Ubiquitin can be characterized by an intricate balance of the ligand- and self-coupled motions, both of which are embedded in the motions of the unbound form.

RNA catalysis through compartmentalization

RNA performs important cellular functions in contemporary life forms. Its ability to act both as a catalyst and a storage mechanism for genetic information is also an important part of the RNA world hypothesis. Compartmentalization within modern cells allows the local concentration of RNA to be controlled and it has been suggested that this was also important in early life forms. Here, we mimic intracellular compartmentalization and macromolecular crowding by partitioning RNA in an aqueous two-phase system (ATPS). We show that the concentration of RNA is enriched by up to 3,000-fold in the dextran-rich phase of a polyethylene glycol/dextran ATPS and demonstrate that this can lead to approximately 70-fold increase in the rate of ribozyme cleavage. This rate enhancement can be tuned by the relative volumes of the two phases in the ATPS. Our observations support the importance of compartmentalization in the attainment of function in an RNA World as well as in modern biology.

Protein stickiness, rather than number of functional protein-protein interactions, predicts expression noise and plasticity in yeast

Background

A hub protein is one that interacts with many functional partners. The annotation of hub proteins, or more generally the protein-protein interaction “degree” of each gene, requires quality genome-wide data. Data obtained using yeast two-hybrid methods contain many false positive interactions between proteins that rarely encounter each other in living cells, and such data have fallen out of favor.

Results

We find that protein “stickiness”, measured as network degree in ostensibly low quality yeast two-hybrid data, is a more predictive genomic metric than the number of functional protein-protein interactions, as assessed by supposedly higher quality high throughput affinity capture mass spectrometry data. In the yeast Saccharomyces cerevisiae, a protein’s high stickiness, but not its high number of functional interactions, predicts low stochastic noise in gene expression, low plasticity of gene expression across different environments, and high probability of forming a homo-oligomer. Our results are robust to a multiple regression analysis correcting for other known predictors including protein abundance, presence of a TATA box and whether a gene is essential. Once the higher stickiness of homo-oligomers is controlled for, we find that homo-oligomers have noisier and more plastic gene expression than other proteins, consistent with a role for homo-oligomerization in mediating robustness.

Conclusions

Our work validates use of the number of yeast two-hybrid interactions as a metric for protein stickiness. Sticky proteins exhibit low stochastic noise in gene expression, and low plasticity in expression across different environments.

The interface of protein structure, protein biophysics, and molecular evolution

Abstract The interface of protein structural biology, protein biophysics, molecular evolution, and molecular population genetics forms the foundations for a mechanistic understanding of many aspects of protein biochemistry. Current efforts in interdisciplinary protein modeling are in their infancy and the state-of-the art of such models is described. Beyond the relationship between amino acid substitution and static protein structure, protein function, and corresponding organismal fitness, other considerations are also discussed. More complex mutational processes such as insertion and deletion and domain rearrangements and even circular permutations should be evaluated. The role of intrinsically disordered proteins is still controversial, but may be increasingly important to consider. Protein geometry and protein dynamics as a deviation from static considerations of protein structure are also important. Protein expression level is known to be a major determinant of evolutionary rate and several considerations including selection at the mRNA level and the role of interaction specificity are discussed. Lastly, the relationship between modeling and needed high-throughput experimental data as well as experimental examination of protein evolution using ancestral sequence resurrection and in vitro biochemistry are presented, towards an aim of ultimately generating better models for biological inference and prediction.

Analytic markovian rates for generalized protein structure evolution

A general understanding of the complex phenomenon of protein evolution requires the accurate description of the constraints that define the sub-space of proteins with mutations that do not appreciably reduce the fitness of the organism. Such constraints can have multiple origins, in this work we present a model for constrained evolutionary trajectories represented by a Markovian process throughout a set of protein-like structures artificially constructed to be topological intermediates between the structure of two natural occurring proteins. The number and type of intermediate steps defines how constrained the total evolutionary process is. By using a coarse-grained representation for the protein structures, we derive an analytic formulation of the transition rates between each of the intermediate structures. The results indicate that compact structures with a high number of hydrogen bonds are more probable and have a higher likelihood to arise during evolution. Knowledge of the transition rates allows for the study of complex evolutionary pathways represented by trajectories through a set of intermediate structures.

Protein misinteraction avoidance causes highly expressed proteins to evolve slowly

The tempo and mode of protein evolution have been central questions in biology. Genomic data have shown a strong influence of the expression level of a protein on its rate of sequence evolution (E-R anticorrelation), which is currently explained by the protein misfolding avoidance hypothesis. Here, we show that this hypothesis does not fully explain the E-R anticorrelation, especially for protein surface residues. We propose that natural selection against protein-protein misinteraction, which wastes functional molecules and is potentially toxic, constrains the evolution of surface residues. Because highly expressed proteins are under stronger pressures to avoid misinteraction, surface residues are expected to show an E-R anticorrelation. Our molecular-level evolutionary simulation and yeast genomic analysis confirm multiple predictions of the hypothesis. These findings show a pluralistic origin of the E-R anticorrelation and reveal the role of protein misinteraction, an inherent property of complex cellular systems, in constraining protein evolution.

Protein-binding dynamics imaged in a living cell

Historically, rate constants were determined in vitro and it was unknown whether they were valid for in vivo biological processes. Here, we bridge this gap by measuring binding dynamics between a pair of proteins in living HeLa cells. Binding of a β-lactamase to its protein inhibitor was initiated by microinjection and monitored by Förster resonance energy transfer. Association rate constants for the wild-type and an electrostatically optimized mutant were only 25% and 50% lower than in vitro values, whereas no change in the rate constant was observed for a slower binding mutant. These changes are much smaller than might be anticipated considering the high macromolecular crowding within the cell. Single-cell analyses of association rate constants and fluorescence recovery after photobleaching reveals a naturally occurring variation in cell density, which is translated to an up to a twofold effect on binding rate constants. The data show that for this model protein interaction the intracellular environment had only a small effect on the association kinetics, justifying the extrapolation of in vitro data to processes in the cell.

Topology of protein interaction network shapes protein abundances and strengths of their functional and nonspecific interactions

How do living cells achieve sufficient abundances of functional protein complexes while minimizing promiscuous nonfunctional interactions? Here we study this problem using a first-principle model of the cell whose phenotypic traits are directly determined from its genome through biophysical properties of protein structures and binding interactions in a crowded cellular environment. The model cell includes three independent prototypical pathways, whose topologies of protein-protein interaction (PPI) subnetworks are different, but whose contributions to the cell fitness are equal. Model cells evolve through genotypic mutations and phenotypic protein copy number variations. We found a strong relationship between evolved physical-chemical properties of protein interactions and their abundances due to a "frustration" effect: Strengthening of functional interactions brings about hydrophobic interfaces, which make proteins prone to promiscuous binding. The balancing act is achieved by lowering concentrations of hub proteins while raising solubilities and abundances of functional monomers. On the basis of these principles we generated and analyzed a possible realization of the proteome-wide PPI network in yeast. In this simulation we found that high-throughput affinity capture-mass spectroscopy experiments can detect functional interactions with high fidelity only for high-abundance proteins while missing most interactions for low-abundance proteins.

The evolution of protein interaction networks

The availability of high-throughput methods to detect protein interactions made construction of comprehensive protein interaction networks for several important model organisms possible. Many studies have since focused on uncovering the structural principles of these networks and relating these structures to biological processes. On a global scale, there are striking similarities in the structure of different protein interaction networks, even when distantly related species, such as the yeast Saccharomyces cerevisiae and the fruit fly Drosophila melanogaster, are compared. However, there is also considerable variance in network structures caused by the gain and loss of genes and mutations which alter the interaction behavior of the encoded proteins. Here, we focus on the current state of knowledge on the structure of protein interaction networks and the evolutionary processes that shaped these structures.