Interaction between photosystem I and ferredoxin. Identification by chemical cross-linking of the polypeptide which binds ferredoxin

Ferredoxin has been effectively cross-linked to photosystem I complex by treatment of purified particles or thylakoids with N-ethyl-3-(3-dimethylaminopropyl)carbodiimide, a zero-length cross-linker which stabilizes protein-protein electrostatic interactions. Analysis of photosystem I polypeptide composition after such a treatment showed a specific decrease of the 20-kDa subunit and the appearance of a new component of about 42 kDa which was recognized by the anti-ferredoxin antibody. Cross-linking of ferredoxin to thylakoids allowed the membrane preparation to photoreduce cytochrome c without requiring exogenous ferredoxin, whereas photosystem I particles purified from treated thylakoids were inactivated in the NADP+ photoreduction activity. From these results, it can be inferred that the polypeptide of 20 kDa is the photosystem I subunit which interacts with ferredoxin during the photosynthetic electron transport.     

Voltammetric measurement of the plastoquinone redox state in isolated thylakoids

The oxidation of plastoquinol by the cytochrome bf complex is commonly believed to be the rate limiting step in photosynthetic electron transport. When input of electrons from PS II exceeds electron flow through the cytochrome bf complex the plastoquinone pool becomes reduced. A voltammetric technique previously used to measure the redox state of the ubiquinone pool in plant mitochondria, was modified to measure the redox state of the plastoquinone pool in thylakoids. The presence or absence of a proton gradient strongly influenced the relationship between the redox state of the plastoquinone pool and other photosynthetic parameters. A linear relationship between the rate of electron transport and the reduction of plastoquinone was found. The slope of this relationship was greater in coupled than in uncoupled thylakoids, indicating that under coupled conditions the plastoquinone pool is more reduced at any given rate of electron flow. A complex relationship was found between QA reduction, calculated as 1 – q_p, and the redox state of the plastoquinone pool. The extent of Q_A reduction was similar in coupled and uncoupled thylakoids, but at any given level of Q_A reduction, PQ was always more reduced in coupled thylakoids. These results suggest that the presence of a proton gradient changes the equilibrium constant between Q_A and PQ.     

Molecular architecture of the thylakoid membrane: lipid diffusion space for plastoquinone

We have determined the stoichiometric composition of membrane components (lipids and proteins) in spinach thylakoids and have derived the molecular area occupied by these components. From this analysis, the lipid phase diffusion space, the fraction of lipids located in the first protein solvation shell (boundary lipids), and the plastoquinone (PQ) concentration are derived. On the basis of these stoichiometric data, we have analyzed the motion of PQ between photosystem (PS) II and cytochrome (cyt.) bf complexes in this highly protein obstructed membrane (protein area about 70%) using percolation theory. This analysis reveals an inefficient diffusion process. We propose that distinct structural features of the thylakoid membrane (grana formation, microdomains) could help to minimize these inefficiencies and ensure a non-rate limiting PQ diffusion process. A large amount of published evidence supports the idea that higher protein associations exist, especially in grana thylakoids. From the quantification of the boundary lipid fraction (about 60%), we conclude that protein complexes involved in these associations should be spaced by lipids. Lipid-spaced protein aggregations in thylakoids are qualitatively different to previously characterized associations (multisubunit complexes, supercomplexes). We derive a hierarchy of protein and lipid interactions in the thylakoid membrane.     

Reduced levels of cytochrome bf complex in transgenic tobacco leads to marked photochemical reduction of the plastoquinone pool,without significant change in acclimation to irradiance

Tobacco transgenics with decreased amounts of the FeS apoprotein were generated using an antisense RNA construct targeted against the nuclear-encoded Rieske FeS protein of cytochrome bf complex [Price et al. (1995) Aust J Plant Physiol 22: 285–297]. FeS phenotypes ranging from intermediate to low were obtained which had 69% and 26% of the Rieske FeS protein of wild type. Similar reductions in the other subunits of cytochrome bf complex, cytochrome f, cytochrome b₅₆₃and the 17 kDa subunit, were demonstrated in the thylakoids of intermediate and low FeS phenotypes. Confirmation that the levels of assembled cytochrome bf in leaves matched the levels of the FeS protein was demonstrated by laser flash-induced redox absorbance changes in leaves, with the extents of cytochrome f oxidation and cytochrome b₅₆₃reduction being equivalent to the decreased amounts of the subunits in isolated thylakoids of the antisense plants. Despite greatly enhanced photochemical reduction of Q_Aand the plastoquinone pool in the antisense plants, light acclimation of the FeS phenotypes to irradiance did not occur. Furthermore, the state 1–state 2 transitions were identical in wild type and antisense plants. Our results suggest that neither Q_Anor the plastoquinone pool acts alone in either the redox control of gene expression or the regulation of light energy distribution between the photosystems. We suggest rather that reduced plastoquinone acting at the inner Q_psite of cytochrome bf complex is involved in molecular redox signalling.     

Effect of training datasets on support vector machine prediction of protein-protein interactions

Knowledge of protein-protein interaction is useful for elucidating protein function via the concept of 'guilt-by-association'. A statistical learning method, Support Vector Machine (SVM), has recently been explored for the prediction of protein-protein interactions using artificial shuffled sequences as hypothetical noninteracting proteins and it has shown promising results (Bock, J. R., Gough, D. A., Bioinformatics 2001, 17, 455-460). It remains unclear however, how the prediction accuracy is affected if real protein sequences are used to represent noninteracting proteins. In this work, this effect is assessed by comparison of the results derived from the use of real protein sequences with that derived from the use of shuffled sequences. The real protein sequences of hypothetical noninteracting proteins are generated from an exclusion analysis in combination with subcellular localization information of interacting proteins found in the Database of Interacting Proteins. Prediction accuracy using real protein sequences is 76.9% compared to 94.1% using artificial shuffled sequences. The discrepancy likely arises from the expected higher level of difficulty for separating two sets of real protein sequences than that for separating a set of real protein sequences from a set of artificial sequences. The use of real protein sequences for training a SVM classification system is expected to give better prediction results in practical cases. This is tested by using both SVM systems for predicting putative protein partners of a set of thioredoxin related proteins. The prediction results are consistent with observations, suggesting that real sequence is more practically useful in development of SVM classification system for facilitating protein-protein interaction prediction.     

Mechanistic Scrutiny Identifies a Kinetic Role for Cytochrome b5 Regulation of Human Cytochrome P450c17 (CYP17A1, P450 17A1)

Cytochrome P450c17 (P450 17A1, CYP17A1) is a critical enzyme in the synthesis of androgens and is now a target enzyme for the treatment of prostate cancer. Cytochrome P450c17 can exhibit either one or two physiological enzymatic activities differentially regulated by cytochrome b5. How this is achieved remains unknown. Here, comprehensive in silico, in vivo and in vitro analyses were undertaken. Fluorescence Resonance Energy Transfer analysis showed close interactions within living cells between cytochrome P450c17 and cytochrome b5. In silico modeling identified the sites of interaction and confirmed that E48 and E49 residues in cytochrome b5 are essential for activity. Quartz crystal microbalance studies identified specific protein-protein interactions in a lipid membrane. Voltammetric analysis revealed that the wild type cytochrome b5, but not a mutated, E48G/E49G cyt b5, altered the kinetics of electron transfer between the electrode and the P450c17. We conclude that cytochrome b5 can influence the electronic conductivity of cytochrome P450c17 via allosteric, protein-protein interactions.     

Diffusion-limited rates for monoclonal antibody binding to cytochrome c.

The kinetic and spectroscopic changes accompanying the binding of two monoclonal antibodies to the oxidized form of horse heart cytochrome c have been investigated. The two epitopes recognized by the antibodies are distinct and noninteracting: antibody 2B5 binds to native cytochrome c near a type II turn (residue 44) while antibody 5F8 binds on the opposite face of the protein near the amino terminus of an alpha-helical segment (residue 60). Antibody-cytochrome c binding obeys a simple bimolecular reaction mechanism with second-order rate constants approaching those expected for diffusion-limited protein-protein interactions. The association rate constants have small activation enthalpies and are inversely dependent on solvent viscosity, as expected for diffusion-controlled reactions. There is a moderate ionic strength dependence of the rate of association between the 2B5 antibody and cytochrome c, with the rate constant increasing about 4-fold as the ionic strength is varied between 0.14 and 0 M. Comparison of the rates for antibody-cytochrome c complex formation for binding to the reduced-native, oxidized-native, and alkaline conformations shows that for MAb 2B5 the forward rate constant depends slightly on cytochrome c conformation. Investigation of the pH-induced transition between the native and alkaline conformational states for free cytochrome c and for antibody-cytochrome c complexes shows that antibody binding stabilizes the native form of the protein.(ABSTRACT TRUNCATED AT 250 WORDS)     

Interactions between Rho GTPases and Rho GDP dissociation inhibitor (Rho-GDI)

The Rho-GDP dissociation inhibitor (Rho-GDI) was used as bait in a two-hybrid screen of a human leucocyte cDNA library. Most of the isolated cDNAs encoded GTPases of the Rho subfamily: RhoA, B, C, Rac1, 2, CDC42 and RhoG. The newly discovered RhoH interacted very poorly with Rho-GDI. Another protein partner shared a homology with RhoA that points to Asp67(RhoA)-Arg68(RhoA)-Leu69(RhoA) as critical for interaction with Rho-GDI. A second screen with RhoA as bait led to the isolation of GDI only. In order to investigate the relative role of protein-protein and protein-lipid interactions between Rho GTPases and Rho-GDI, CAAX box mutants of RhoA were produced. They were found to interact with Rho-GDI as efficiently as wild type RhoA, indicating that protein-protein interactions alone lead to strong binding of the two proteins. The C-terminal polybasic region of RhoA was also shown to be a site of protein-protein interaction with Rho-GDI.     

Protein-protein interactions in concentrated electrolyte solutions

Protein-protein interactions were measured for ovalbumin and for lysozyme in aqueous salt solutions. Protein-protein interactions are correlated with a proposed potential of mean force equal to the free energy to desolvate the protein surface that is made inaccessible to the solvent due to the protein-protein interaction. This energy is calculated from the surface free energy of the protein that is determined from protein-salt preferential-interaction parameter measurements. In classical salting-out behavior, the protein-salt preferential interaction is unfavorable. Because addition of salt raises the surface free energy of the protein according to the surface-tension increment of the salt, protein-protein attraction increases, leading to a reduction in solubility. When the surface chemistry of proteins is altered by binding of a specific ion, salting-in is observed when the interactions between (kosmotrope) ion-protein complexes are more repulsive than those between the uncomplexed proteins. However, salting-out is observed when interactions between (chaotrope) ion-protein complexes are more attractive than those of the uncomplexed proteins.     

A de novo protein binding pair by computational design and directed evolution

The de novo design of protein-protein interfaces is a stringent test of our understanding of the principles underlying protein-protein interactions and would enable unique approaches to biological and medical challenges. Here we describe a motif-based method to computationally design protein-protein complexes with native-like interface composition and interaction density. Using this method we designed a pair of proteins, Prb and Pdar, that heterodimerize with a Kd of 130 nM, 1000-fold tighter than any previously designed de novo protein-protein complex. Directed evolution identified two point mutations that improve affinity to 180 pM. Crystal structures of an affinity-matured complex reveal binding is entirely through the designed interface residues. Surprisingly, in the in vitro evolved complex one of the partners is rotated 180° relative to the original design model, yet still maintains the central computationally designed hotspot interaction and preserves the character of many peripheral interactions. This work demonstrates that high-affinity protein interfaces can be created by designing complementary interaction surfaces on two noninteracting partners and underscores remaining challenges.     

Control of the photosynthetic electron transport by PQ diffusion microdomains in thylakoids of higher plants

We investigate the role of plastoquinone (PQ) diffusion in the control of the photosynthetic electron transport. A control analysis reveals an unexpected flux control of the whole chain electron transport by photosystem (PS) II. The contribution of PSII to the flux control of whole chain electron transport was high in stacked thylakoids (control coefficient, CJ(PSII) =0.85), but decreased after destacking (CJ(PSII)=0.25). From an 'electron storage' experiment, we conclude that in stacked thylakoids only about 50 to 60% of photoreducable PQ is involved in the light-saturated linear electron transport. No redox equilibration throughout the membrane between fixed redox groups at PSII and cytochrome (cyt) bf complexes, and the diffusable carrier PQ is achieved. The data support the PQ diffusion microdomain concept by Lavergne et al. [J. Lavergne, J.-P. Bouchaud, P. Joliot, Biochim. Biophys. Acta 1101 (1992) 13-22], but we come to different conclusions about size, structure and size distribution of domains. From an analysis of cyt b6 reduction, as a function of PSII inhibition, we conclude that in stacked thylakoids about 70% of PSII is located in small domains, where only 1 to 2 PSII share a local pool of a few PQ molecules. Thirty percent of PSII is located in larger domains. No small domains were found in destacked thylakoids. We present a structural model assuming a hierarchy of specific, strong and weak interactions between PSII core, light harvesting complexes (LHC) II and cyt bf. Peripheral LHCII's may serve to connect PSII-LHCII supercomplexes to a flexible protein network, by which small closed lipid diffusion compartments are formed. Within each domain, PQ moves rapidly and shuttles electrons between PSII and cyt bf complexes in the close vicinity. At the same time, long range diffusion is slow. We conclude, that in high light, cyt bfcomplexes located in distant stromal lamellae (20 to 30%) are not involved in the linear electron transport.     

Prediction-based fingerprints of protein-protein interactions

The recognition of protein interaction sites is an important intermediate step toward identification of functionally relevant residues and understanding protein function, facilitating experimental efforts in that regard. Toward that goal, the authors propose a novel representation for the recognition of protein-protein interaction sites that integrates enhanced relative solvent accessibility (RSA) predictions with high resolution structural data. An observation that RSA predictions are biased toward the level of surface exposure consistent with protein complexes led the authors to investigate the difference between the predicted and actual (i.e., observed in an unbound structure) RSA of an amino acid residue as a fingerprint of interaction sites. The authors demonstrate that RSA prediction-based fingerprints of protein interactions significantly improve the discrimination between interacting and noninteracting sites, compared with evolutionary conservation, physicochemical characteristics, structure-derived and other features considered before. On the basis of these observations, the authors developed a new method for the prediction of protein-protein interaction sites, using machine learning approaches to combine the most informative features into the final predictor. For training and validation, the authors used several large sets of protein complexes and derived from them nonredundant representative chains, with interaction sites mapped from multiple complexes. Alternative machine learning techniques are used, including Support Vector Machines and Neural Networks, so as to evaluate the relative effects of the choice of a representation and a specific learning algorithm. The effects of induced fit and uncertainty of the negative (noninteracting) class assignment are also evaluated. Several representative methods from the literature are reimplemented to enable direct comparison of the results. Using rigorous validation protocols, the authors estimated that the new method yields the overall classification accuracy of about 74% and Matthews correlation coefficients of 0.42, as opposed to up to 70% classification accuracy and up to 0.3 Matthews correlation coefficient for methods that do not utilize RSA prediction-based fingerprints. The new method is available at http://sppider.cchmc.org.     

A beta-galactosidase-based bacterial two-hybrid system to assess protein-protein interactions in the correct cellular environment

The vast majority of proteins functions in complex with one or more of the same or other proteins, indicating that protein-protein interactions play crucial roles in biology. Here, we present a beta-galactosidase reconstitution-based bacterial two-hybrid system in which two proteins of interest are fused to two non-functional but complementing beta-galactosidase truncations (Delta alpha and Delta omega). The level of complemented beta-galactosidase activity, driven by the protein-protein recognition between both non-beta-galactosidase parts of the chimeras, reflects whether or not the proteins of interest interact. Our approach was validated by reconfirming some well-established Escherichia coli cytoplasmic and membranous interactions, including well-chosen mutants, and providing the first in vivo evidence for the transient periplasmic interaction between Rhodobacter capsulatus cytochrome c2 and cytochrome c peroxidase. We demonstrated the major advantages of this in vivo two-hybrid technique: i) analyses of interactions are not limited to particular cellular compartments, ii) the potential of using the system in mutation-driven structure-function studies, and iii) the possibility of its application to transiently interacting proteins. These benefits demonstrate the relevance of the method as a powerful new tool in the broad spectrum of interaction assessment methods.     

Aib-based peptide backbone as scaffolds for helical peptide mimics.

Helical peptides that can intervene and disrupt therapeutically important protein-protein interactions are attractive drug targets. In order to develop a general strategy for developing such helical peptide mimics, we have studied the effect of incorporating alpha-amino isobutyric acid (Aib), an amino acid with strong preference for helical backbone, as the sole helix promoter in designed peptides. Specifically, we focus on the hdm2-p53 interaction, which is central to development of many types of cancer. The peptide corresponding to the hdm2 interacting part of p53, helical in bound state but devoid of structure in solution, served as the starting point for peptide design that involved replacement of noninteracting residues by Aib. Incorporation of Aib, while preserving the interacting residues, led to significant increase in helical structure, particularly at the C-terminal region as judged by nuclear magnetic resonance and circular dichroism. The interaction with hdm2 was also found to be enhanced. Most interestingly, trypsin cleavage was found to be retarded by several orders of magnitude. We conclude that incorporation of Aib is a feasible strategy to create peptide helical mimics with enhanced receptor binding and lower protease cleavage rate.     

Isothermal titration calorimetry of protein-protein interactions.

The interaction of biologicalmacromolecules, whether protein-DNA, antibody-antigen, hormone-receptor, etc., illustrates the complexity and diversity of molecular recognition. The importance of such interactions in the immune response, signal transduction cascades, and gene expression cannot be overstated. It is of great interest to determine the nature of the forces that stabilize the interaction. The thermodynamics of association are characterized by the stoichiometry of the interaction (n), the association constant (K(a)), the free energy (DeltaG(b)), enthalpy (DeltaH(b)), entropy (DeltaS(b)), and heat capacity of binding (DeltaC(p)). In combination with structural information, the energetics of binding can provide a complete dissection of the interaction and aid in identifying the most important regions of the interface and the energetic contributions. Various indirect methods (ELISA, RIA, surface plasmon resonance, etc.) are routinely used to characterize biologically important interactions. Here we describe the use of isothermal titration calorimetry (ITC) in the study of protein-protein interactions. ITC is the most quantitative means available for measuring the thermodynamic properties of a protein-protein interaction. ITC measures the binding equilibrium directly by determining the heat evolved on association of a ligand with its binding partner. In a single experiment, the values of the binding constant (K(a)), the stoichiometry (n), and the enthalpy of binding (DeltaH(b)) are determined. The free energy and entropy of binding are determined from the association constant. The temperature dependence of the DeltaH(b) parameter, measured by performing the titration at varying temperatures, describes the DeltaC(p) term. As a practical application of the method, we describe the use of ITC to study the interaction between cytochrome c and two monoclonal antibodies.     

The interaction of cytochrome c with monolayers of phosphatidylethanolamine

1. The interaction between [(14)C]carboxymethylated cytochrome c and monolayers of egg phosphatidylethanolamine at the air/water interface has been investigated by measurements of surface radioactivity, pressure and potential. 2. On adding (14)C-labelled cytochrome c to the subphase under monolayers with a surface pressure below 24dynes/cm. there was an initial surface pressure increment as the protein penetrated, followed by an adsorption that could be detected only by a continued increase in the surface radioactivity. 3. Above film pressures of 24dynes/cm. only adsorption was observed, i.e. an increment in surface radioactivity with none in surface pressure. 4. The changes in surface parameters with penetration of cytochrome c added to the subphase were indirectly proportional to the initial pressure of the monolayer. With hydrogenated phosphatidylethanolamine the constant of proportionality was increased but penetration again ceased at 24dynes/cm. 5. On compressing a phosphatidylethanolamine film containing penetrated cytochrome c to 40dynes/cm. only a proportion of the protein was ejected on a subphase of 10mm-sodium chloride, whereas on a subphase of m-sodium chloride nearly all the protein was lost. 6. With both penetration and adsorption only a small proportion of the added cytochrome c interacted with the phospholipid films, and initially the amount bound was proportional to the added protein concentration. There was no evidence of a stoicheiometric relationship between the protein and phospholipid or the build-up of multilayers. The bonded protein was not released by removing cytochrome c from the subphase. 7. The addition of m-sodium chloride to the subphase delays the rate of protein penetration into low-pressure films, but the final surface-pressure increment is not appreciably decreased. In contrast, m-sodium chloride almost completely stops adsorption on to films at all pressures. 8. When sodium chloride is added to the subphase below cytochrome c adsorbed to monolayers at high pressures, so that the final concentration is 1m, only a proportion of the protein is desorbed and this decreases as the time of the interaction increases. This indicates that adsorption is initially electrostatic, followed by the formation of non-ionic bonds. 9. Alteration of the subphase pH under a high-pressure film leads to a steady increase in adsorption from pH3 to 8.5 followed by a rapid fall to zero adsorption at pH11. 10. The penetration into phospholipid monolayers at 10dynes/cm. shows a rate that is consistent with the relative electrostatic status of the two components of the interaction as the subphase pH is varied between 3 and 10.5. The final equilibrium penetration shows a pronounced peak in the increments of surface pressure at pH9.0 although a similar peak is not observed in the surface radioactivity. This indicates that more residues of the protein are penetrating into the film at about this pH. 11. Determinations were made of the electrophoretic mobilities of phosphatidylethanolamine particles both alone and after interaction with cytochrome c. 12. The electrophoretic mobilities of cytochrome c adsorbed on lipid particles showed an isoelectric point below that of cytochrome c. This and the observations on the monolayers suggest that, with cytochrome c, protein-protein interactions are weak compared with other proteins.     

Protein microarrays, biosensors, and cell-based methods for secretome-wide extracellular protein-protein interaction mapping

Approximately one quarter of all human genes encode proteins that function in the extracellular space or serve to bridge the extracellular and intracellular environments. Physical associations between these secretome proteins serve to regulate a wide range of biological activities and consequently represent important therapeutic targets. Moreover, some extracellular proteins are targeted by pathogens to allow host access or immune evasion. Despite the importance of extracellular protein-protein interactions, our knowledge in this area has remained sparse. Weak affinities and low abundance have often hindered efforts to identify these interactions using traditional methods such as biochemical purification and cDNA library expression cloning. Moreover, current large-scale protein-protein interaction mapping techniques largely under represent extracellular protein-protein interactions. This review highlights emerging biosensor and protein microarray technology, along with more traditional cell-based techniques, that are compatible with secretome-wide screens for extracellular protein-protein interaction discovery. A combination of these approaches will serve to rapidly expand our knowledge of the extracellular protein-protein interactome.     

Electrostatic interactions during electron transfer reactions between c-type cytochromes and flavodoxin

The interaction of three different c-type cytochromes with flavodoxin has been studied by computer graphics modelling and computational methods. Flavodoxin and each cytochrome can make similar hypothetical electron transfer complexes that are characterized by nearly coplanar arrangement of the prosthetic groups, close intermolecular contacts at the protein-protein interface, and complementary intermolecular salt linkages. Computation of the electrostatic free energy of each complex showed that all were electrostatically stable. However, both the magnitude and behavior of the electrostatic stabilization as a function of solution ionic strength differed for the three cytochrome c-flavodoxin complexes. Variation in the computed electrostatic stabilization appears to reflect differences in the surface distribution of all charged groups in the complex, rather than differences localized at the site of intermolecular contact. The computed electrostatic association constants for the complexes and the measured kinetic rates of electron transfer in solution show a remarkable similarity in their ionic strength dependence. This correlation suggests electrostatic interactions influence electron transfer rates between protein molecules at the intermolecular association step. Comparative calculations for the three cytochrome c-flavodoxin complexes show that these ionic strength effects also involve all charged groups in both redox partners.     

Competing species in a changing climate: effects of recruitment disturbances on two interacting barnacle species

1. The climate is changing and data-based simulation models can be a valuable tool for predicting population response to such changes and investigate the mechanisms of population change. In this study, a data-based two-species matrix model was constructed to explore the possible effects of elevated sea surface temperature (i.e. climate change) on the interaction between open populations of the south Atlantic barnacle species Chthamalus montagui and the boreal species Semibalanus balanoides in the north-east Atlantic. 2. First, the model was used to perform an elasticity analysis to determine the relative importance of recruitment and survival in the interaction. Further, three scenarios of changes in recruitment, related to climate change, were investigated with model simulations: (i) increased frequencies of low recruitment for S. balanoides; (ii) increased frequencies of high recruitment for C. montagui; (iii) a combination of (i) and (ii). 3. Model simulations showed that in present environmental conditions, S. balanoides occupied most of the space and dominated the interaction through high recruitment and survival. These results matched independent field observations, which validated the model for further analyses. 4. The elasticity analyses showed that although free space was available there was competition for space during recruitment intervals. It was also shown that both populations were sensitive to changes in recruitment. 5. Introducing the three scenarios of recruitment disturbances led to large changes in species abundance and free space. The most significant changes were found when scenario (i) and (ii) were combined, producing a shift in species dynamics towards C. montagui dominance. This demonstrates that recruitment can be an important mechanism in the interaction between populations and that the population response to changes in recruitment depends on the added response of interacting species. 6. In a more general context, this model shows that increased sea surface temperature could rapidly lead to increased competition from southern species at higher latitudes. This might accelerate the effects of climate change on the species distribution at these latitudes and eventually lead to changes in community dynamics on temperate and subarctic shores.     

Potential of mean force for protein-protein interaction studies.

Calculating protein-protein interaction energies is crucial for understanding protein-protein associations. On the basis of the methodology of mean-field potential, we have developed an empirical approach to estimate binding free energy for protein-protein interactions. This knowledge-based approach has been used to derive distance-dependent free energies of protein complexes from a nonredundant training set in the Protein Data Bank (PDB), with a careful treatment of homology. We calculate atom pair potentials for 16 pair interactions, which can reflect the importance of hydrophobic interactions and specific hydrogen-bonding interactions. The derived potentials for hydrogen-bonding interactions show a valley of favorable interactions at a distance of approximately 3 A, corresponding to that of an established hydrogen bond. For the test set of 28 protein complexes, the calculated energies have a correlation coefficient of 0.75 compared with experimental binding free energies. The performance of the method in ranking the binding energies of different protein-protein complexes shows that the energy estimation can be applied to value binding free energies for protein-protein associations.