Protein Interactions and Fluctuations in a Proteomic Network using an Elastic Network Model

Abstract

A set of protein conformations are analyzed by normal mode analysis. An elastic network model is used to obtain fluctuation and cooperativity of residues with low amplitude fluctuations across different species. Slow modes that are associated with the function of proteins have common features among different protein structures. We show that the degree of flexibility of the protein is important for proteins to interact with other proteins and as the species gets more complex its proteins become more flexible. In the complex organism, higher cooperativity arises due to protein structure and connectivity.     

Probing possible downhill folding: native contact topology likely places a significant constraint on the folding cooperativity of proteins with approximately 40 residues

Experiments point to appreciable variations in folding cooperativity among natural proteins with approximately 40 residues, indicating that the behaviors of these proteins are valuable for delineating the contributing factors to cooperative folding. To explore the role of native topology in a protein's propensity to fold cooperatively and how native topology might constrain the degree of cooperativity achievable by a given set of physical interactions, we compared folding/unfolding kinetics simulated using three classes of native-centric C^α chain models with different interaction schemes. The approach was applied to two homologous 45-residue fragments from the peripheral subunit-binding domain family and a 39-residue fragment of the N-terminal domain of ribosomal protein L9. Free-energy profiles as functions of native contact number were computed to assess the heights of thermodynamic barriers to folding. In addition, chevron plots of folding/unfolding rates were constructed as functions of native stability to facilitate comparison with available experimental data. Although common Gō-like models with pairwise Lennard-Jones-type interactions generally fold less cooperatively than real proteins, the rank ordering of cooperativity predicted by these models is consistent with experiment for the proteins investigated, showing increasing folding cooperativity with increasing nonlocality of a protein's native contacts. Models that account for water-expulsion (desolvation) barriers and models with many-body (nonadditive) interactions generally entail higher degrees of folding cooperativity indicated by more linear model chevron plots, but the rank ordering of cooperativity remains unchanged. A robust, experimentally valid rank ordering of model folding cooperativity independent of the multiple native-centric interaction schemes tested here argues that native topology places significant constraints on how cooperatively a protein can fold.     

Sequential hydration of dry proteins: a direct difference IR investigation of sequence homologs lysozyme and alpha- lactalbumin

Direct difference ir spectra are presented as a function of hydration for lysozyme and α-lactalbumin, and detailed sequential hydration molecular events identified. Despite the strong sequence homology between the two proteins, and their expected conformational similarity, the hydration behaviour of the polar groups is different for the two proteins. Using a Hill-type analysis, we conclude that the acid groups ionize and hydrate rapidly and noncooperatively in both proteins, consistent with the known (lysozyme) and postulated (α-lactalbumin) surface chemistry. The polar group hydration shows a clear cooperativity, which is quantitatively different in the two proteins. Complementary work suggests this cooperativity relates to a hydration-induced “loosening up” of the lysozyme conformation at about 55 mol water/mol protein. α-Lactalbumin appears to “open up” more easily for hydration than does lysozyme, consistent with its lower stability against thermal and acid denaturation.     

Model for cooperativity of biological membranes

We present a mathematical model for the complex cooperativity observed in biological membranes. In our model, it is assumed that the proteins bound on the membrane are noncooperative and possess a Bohr proton. It is further assumed that the net charge of the unliganded state of the protein is different from that of the liganded state owing to the structural change upon binding the ligand. With this model, we show how an all-or-none response, a graded response, and a noncooperative response arise in the binding curve of such biological membranes. In addition, we show how an effector, which can alter the pK_a involved in the binding site, induces a complex cooperativity.     

Highly sequential binding of protein kinase C and related proteins to membranes

Protein kinase C belongs to a class of proteins that displays simultaneous interaction with calcium and phospholipids. Other members of this class include two proteins (Mr 64K and 32K) isolated from bovine brain. The association of these proteins with membranes exhibited highly unusual properties that were not consistent with a simple equilibrium. Titration of protein-phospholipid binding as a function of calcium showed an apparently normal curve with a low degree of cooperativity. The binding was rapid and quickly adjusted to changes in the calcium concentration. Calcium was readily exchanged from the protein-phospholipid complex. However, at each calcium concentration, membrane-bound protein was not in rapid equilibrium with free protein in solution; the half-time for dissociation exceeded 24 h. Titration of phospholipid vesicles with proteins showed different saturation levels of bound protein at different calcium concentrations. The amount of protein bound was almost entirely determined by the concentration of calcium and was virtually unaffected by the free protein concentration. These properties suggested that protein-phospholipid binding involved a sequence of steps that were each irreversible upon completion. These binding properties were consistent with high-affinity interaction between protein and phospholipid, high cooperativity with respect to calcium (N greater than or equal to 10), clustering of acidic phospholipids, and negative cooperativity with respect to protein density on the membrane. A major apparent problem with the complete titration of PKC-membrane interaction was a requirement for calcium in excess of intracellular levels. However, a highly sequential binding process showed that a number of protein-binding sites on the membrane would be saturated with calcium at physiological levels.(ABSTRACT TRUNCATED AT 250 WORDS)     

Protein complex evolution does not involve extensive network rewiring

The formation of proteins into stable protein complexes plays a fundamental role in the operation of the cell. The study of the degree of evolutionary conservation of protein complexes between species and the evolution of protein-protein interactions has been hampered by lack of comprehensive coverage of the high-throughput (HTP) technologies that measure the interactome. We show that new high-throughput datasets on protein co-purification in yeast have a substantially lower false negative rate than previous datasets when compared to known complexes. These datasets are therefore more suitable to estimate the conservation of protein complex membership than hitherto possible. We perform comparative genomics between curated protein complexes from human and the HTP data in Saccharomyces cerevisiae to study the evolution of co-complex memberships. This analysis revealed that out of the 5,960 protein pairs that are part of the same complex in human, 2,216 are absent because both proteins lack an ortholog in S. cerevisiae, while for 1,828 the co-complex membership is disrupted because one of the two proteins lacks an ortholog. For the remaining 1,916 protein pairs, only 10% were never co-purified in the large-scale experiments. This implies a conservation level of co-complex membership of 90% when the genes coding for the protein pairs that participate in the same protein complex are also conserved. We conclude that the evolutionary dynamics of protein complexes are, by and large, not the result of network rewiring (i.e. acquisition or loss of co-complex memberships), but mainly due to genomic acquisition or loss of genes coding for subunits. We thus reveal evidence for the tight interrelation of genomic and network evolution.     

Determination of binding constants for cooperative site-specific protein-DNA interactions using the gel mobility-shift assay

We have investigated the question of whether the gel mobility-shift assay can provide data that are useful to the demonstration of cooperativity in the site-specific binding of proteins to DNA. Three common patterns of protein-DNA interaction were considered: (i) the cooperative binding of a protein to two sites (illustrated by the Escherichia coli Gal repressor); (ii) the cooperative binding of a bidentate protein to two sites (illustrated by the E. coli Lac repressor); and (iii) the cooperative binding of a protein to three sites (illustrated by the lambda cI repressor). A simple, rigorous, and easily extendable statistical mechanical approach to the derivation of the binding equations for the different patterns is presented. Both simulated and experimental data for each case are analyzed. The mobility-shift assay provides estimates of the macroscopic binding constants for each step of ligation based on its separation of liganded species by the number of ligands bound. Resolution of the binding constants depends on the precision with which the equilibrium distribution of liganded species is determined over the entire range of titration of each of the sites. However, the evaluation of cooperativity from the macroscopic binding constants is meaningful only for data that are also accurate. Some criteria that are useful in evaluating accuracy are introduced and illustrated. Resolution of cooperative effects is robust only for the simplest case, in which there are two identical protein binding sites. In this case, cooperative effects of up to 1,000-fold are precisely determined. For heterogeneous sites, cooperative effects of greater than 1,000-fold are resolvable, but weak cooperativity is masked by the heterogeneity. For three-site systems, only averaged pair-wise cooperative effects are resolvable.     

Hierarchical regulation of WASP/WAVE proteins

Members of the Wiskott-Aldrich syndrome protein (WASP) family control actin dynamics in eukaryotic cells by stimulating the actin nucleating activity of the Arp2/3 complex. The prevailing paradigm for WASP regulation invokes allosteric relief of autoinhibition by diverse upstream activators. Here we demonstrate an additional level of regulation that is superimposed upon allostery: dimerization increases the affinity of active WASP species for Arp2/3 complex by up to 180-fold, greatly enhancing actin assembly by this system. This finding explains a large and apparently disparate set of observations under a common mechanistic framework. These include WASP activation by the bacterial effector EspFu and a large number of SH3 domain proteins, the effects on WASP of membrane localization/clustering and assembly into large complexes, and cooperativity between different family members. Allostery and dimerization act in hierarchical fashion, enabling WASP/WAVE proteins to integrate different classes of inputs to produce a wide range of cellular actin responses.     

Adsorption of globular proteins on locally planar surfaces: models for the effect of excluded surface area and aggregation of adsorbed protein on adsorption equilibria.

Equilibrium statistical-thermodynamic models are presented for the surface adsorption of proteins modeled as regular convex hard particles. The adsorbed phase is treated as a two-dimensional fluid, and the chemical potential of adsorbed protein is obtained from scaled particle theory. Adsorption isotherms are calculated for nonassociating and self-associating adsorbing proteins. Area exclusion broadens adsorption isotherms relative to the Langmuir isotherm (negative cooperativity), whereas self-association steepens them (positive cooperativity). The calculated isotherm for adsorption of hard spheres using scaled particle theory for hard discs agrees well with that calculated from the hard disc virial expansion. As the cross section of the adsorbing protein in the plane of the surface becomes less discoidal, the apparent negative cooperativity manifested in the isotherm becomes more pronounced. The model is extended to the case of simultaneous adsorption of a tracer protein at low saturation and a competitor protein with a different size and/or shape at arbitrary fractional saturation. Area exclusion by competitor for tracer (and vice versa) is shown to substantially enhance the displacement of tracer by competitor and to qualitatively invalidate the standard interpretation of ligand competition experiments, according to which the fractional displacement of tracer by competitor is equal to the fractional saturation by competitor.     

NMR-based modeling and binding studies of a ternary complex between chicken liver bile acid binding protein and bile acids

Chicken liver bile acid binding protein (cL-BABP) is involved in bile acid transport in the liver cytosol. A detailed study of the mechanism of binding and selectivity of bile acids binding proteins towards the physiological pool of bile salts is a key issue for the complete understanding of the role of these proteins and their involvement in cholesterol homeostasis. In the present study, we modeled the ternary complex of cL-BABP with two molecules of bile salts using the data driven docking program HADDOCK on the basis of NMR and mass spectrometry data. Docking resulted in good 3D models, satisfying the majority of experimental restraints. The docking procedure represents a necessary step to help in the structure determination and in functional analysis of such systems, in view of the high complexity of the 3D structure determination of a ternary complex with two identical ligands. HADDOCK models show that residues involved in binding are mainly located in the C-terminal end of the protein, with two loops, CD and EF, playing a major role in ligand binding. A spine, comprising polarresidues pointing toward the protein interior and involved in motion communication, has a prominent role in ligand interaction. The modeling approach has been complemented with NMR interaction and competition studies of cL-BABP with chenodeoxycholic and cholic acids. A higher affinity for chenodeoxycholic acid was observed and a Kd upper limit estimate was obtained. The binding is highly cooperative and no site selectivity was detected for the different bile salts, thus indicating that site selectivity and cooperativity are not correlated. Differences in physiological pathways and bile salt pools in different species is discussed in light of the binding results thus enlarging the body of knowledge of BABPs biological functions.     

Kinetic mechanism of ligand binding in human ileal bile acid binding protein as determined by stopped-flow fluorescence analysis

Cooperative ligand binding to human ileal bile acid binding protein (I-BABP) was studied using the stopped-flow fluorescence technique. The kinetic data obtained for wild-type protein are in agreement with a four-step mechanism where after a fast conformational change on the millisecond time scale, the ligands bind in a sequential manner, followed by another, slow conformational change on the time scale of seconds. This last step is more pronounced in the case of glycocholate (GCA), the bile salt that binds with high positive cooperativity and is absent in mutant I-BABP proteins that lack positive cooperativity in their bile salt binding. These results suggest that positive cooperativity in human I-BABP is related to a slow conformational change of the protein, which occurs after the second binding step. Analogous to that in the intestinal fatty acid binding protein (I-FABP), we hypothesize that ligand binding in I-BABP is linked to a disorder-order transition between an open and a closed form of the protein.     

Determinants of cooperativity and site selectivity in human ileal bile acid binding protein

Human ileal bile acid binding protein (I-BABP) is a member of the family of intracellular lipid-binding proteins and is thought to play a role in the enterohepatic circulation of bile salts. Our group has previously shown that human I-BABP binds two molecules of glycocholate (GCA) with low intrinsic affinity but an extraordinary high degree of positive cooperativity. Besides the strong positive cooperativity, human I-BABP exhibits a high degree of site selectivity in its interactions with GCA and glycochenodeoxycholate (GCDA), the two major bile salts in humans. In this study, on the basis of our first generation nuclear magnetic resonance (NMR) structure of the ternary complex of human I-BABP with GCA and GCDA, we introduced single-residue mutations at certain key positions in the binding pocket that might disrupt a hydrogen-bonding network, a likely way of energetic communication between the two sites. Macroscopic binding parameters were determined using isothermal titration calorimetry, and site selectivity was monitored by NMR spectroscopy of isotopically enriched bile salts. According to our results, cooperativity and site selectivity are not linked in human I-BABP. While cooperativity is governed by a subtle interplay of entropic and enthalpic contributions, site selectivity appears to be determined by more localized enthalpic effects. Possible communication pathways between the two binding sites are discussed.     

A complex-centric view of protein network evolution

The recent availability of protein–protein interaction networks for several species makes it possible to study protein complexes in an evolutionary context. In this article, we present a novel network-based framework for reconstructing the evolutionary history of protein complexes. Our analysis is based on generalizing evolutionary measures for single proteins to the level of whole subnetworks, comprehensively considering a broad set of computationally derived complexes and accounting for both sequence and interaction changes. Specifically, we compute sets of orthologous complexes across species, and use these to derive evolutionary rate and age measures for protein complexes. We observe significant correlations between the evolutionary properties of a complex and those of its member proteins, suggesting that protein complexes form early in evolution and evolve as coherent units. Additionally, our approach enables us to directly quantify the extent to which gene duplication has played a role in the evolution of complexes. We find that about one quarter of the sets of orthologous complexes have originated from evolutionary cores of homodimers that underwent duplication and divergence, testifying to the important role of gene duplication in protein complex evolution.     

Cooperativity in the motor activities of the ATP-fueled molecular motors

Kinesin, myosin and F1-ATPase are multi-domain molecular motors with multiple catalytic subunits. The motor mechanochemics are achieved via the conversion of ATP hydrolysis energy into forces and motions. We find that the catalysis of these molecular motors do not follow the simple Michaelis-Menten mechanism. The motor activities, such as the hydrolysis or processive rates, of kinesin, myosin and F1-ATPase have a complex ATP-dependent cooperativity. To understand this complexity in kinetics and mechanochemics, we develop a conformation correlation theory of cooperativity for the ATP-fueled motor proteins. The quantitative analysis and simulations indicate that cooperativity is induced by the conformational coupling of binding states of different subunits and prevails in the motor activities.     

Cooperative folding in a multi-domain protein

Most protein domains are found in multi-domain proteins, yet most studies of protein folding have concentrated on small, single-domain proteins or on isolated domains from larger proteins. Spectrin domains are small (106 amino acid residues), independently folding domains consisting of three long alpha-helices. They are found in multi-domain proteins with a number of spectrin domains in tandem array. Structural studies have shown that in these arrays the last helix of one domain forms a continuous helix with the first helix of the following domain. It has been demonstrated that a number of spectrin domains are stabilised by their neighbours. Here we investigate the molecular basis for cooperativity between adjacent spectrin domains 16 and 17 from chicken brain alpha-spectrin (R16 and R17). We show that whereas the proteins unfold as a single cooperative unit at 25 degrees C, cooperativity is lost at higher temperatures and in the presence of stabilising salts. Mutations in the linker region also cause the cooperativity to be lost. However, the cooperativity does not rely on specific interactions in the linker region alone. Most mutations in the R17 domain cause a decrease in cooperativity, whereas proteins with mutations in the R16 domain still fold cooperatively. We propose a mechanism for this behaviour.     

Unraveling Subunit Cooperativity in Homotetrameric HCN2 Channels

In a multimeric receptor protein, the binding of a ligand can modulate the binding of a succeeding ligand. This phenomenon, called cooperativity, is caused by the interaction of the receptor subunits. By using a complex Markovian model and a set of parameters determined previously, we analyzed how the successive binding of four ligands leads to a complex cooperative interaction of the subunits in homotetrameric HCN2 pacemaker channels. The individual steps in the model were characterized by Gibbs free energies for the equilibria and activation energies, specifying the affinity of the binding sites and the transition rates, respectively. Moreover, cooperative free energies were calculated for each binding step in both the closed and the open channel. We show that the cooperativity sequence positive-negative-positive determined for the binding affinity is generated by the combined effect of very different cooperativity sequences determined for the binding and unbinding rates, which are negative-negative-positive and no-negative-no, respectively. It is concluded that in the ligand-induced activation of HCN2 channels, the sequence of cooperativity based on the binding affinity is caused by two even qualitatively different sequences of cooperativity that are based on the rates of ligand binding and unbinding.