Predicting the rate enhancement of protein complex formation from the electrostatic energy of interaction

Residue interaction network and molecular dynamics simulation study on the binding of S239D/I332E Fc variant with enhanced affinity to FcγRIIIa receptor

Engineering of Fc has been adapted as an efficient method for enhanced or reduced affinity towards Fc receptors in the development of therapeutic antibodies. S239D/I332E mutation of Fc induces approximately two logs greater affinity to the FcγRIIIa receptor and has been extensively employed in various Fc engineering studies. It is known that the mutation gives rise to the formation of salt bridges between the mutated residues of Fc and FcγRIIIa, but the overall effect of the mutation in the binding interface of the Fc-FcγRIIIa complex is still unclear. In this study, the molecular interactions in the binding interface of mutant Fc and FcγRIIIa were analyzed and compared with those of wild-type Fc binding through residue interaction network (RIN) analysis and molecular dynamics (MD) simulation. RIN analysis identified specific molecular interactions and Hub residues in the interfaces, and their numbers were increased by introducing the mutation, with maintaining most of the molecular interactions in the wild-type complex. MD simulation study revealed that the numbers of stable electrostatic interactions and stable Hub residues in the mutant complex were higher than those in the wild-type complex. The introduced mutations were shown to form further charge-charge attractive interactions in addition to the identified salt bridges without generating any repulsive interactions. These results are expected to provide further structural insight into Fc variants' design based on the S239D/I332E mutation.

Stress-Responsive Periplasmic Chaperones in Bacteria

Periplasmic proteins are involved in a wide range of bacterial functions, including motility, biofilm formation, sensing environmental cues, and small-molecule transport. In addition, a wide range of outer membrane proteins and proteins that are secreted into the media must travel through the periplasm to reach their final destinations. Since the porous outer membrane allows for the free diffusion of small molecules, periplasmic proteins and those that travel through this compartment are more vulnerable to external environmental changes, including those that result in protein unfolding, than cytoplasmic proteins are. To enable bacterial survival under various stress conditions, a robust protein quality control system is required in the periplasm. In this review, we focus on several periplasmic chaperones that are stress responsive, including Spy, which responds to envelope-stress, DegP, which responds to temperature to modulate chaperone/protease activity, HdeA and HdeB, which respond to acid stress, and UgpB, which functions as a bile-responsive chaperone.

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

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

Co-localization and confinement of ecto-nucleotidases modulate extracellular adenosine nucleotide distributions

Nucleotides comprise small molecules that perform critical signaling roles in biological systems. Adenosine-based nucleotides, including adenosine tri-, di-, and mono-phosphate, are controlled through their rapid degradation by diphosphohydrolases and ecto-nucleotidases (NDAs). The interplay between nucleotide signaling and degradation is especially important in synapses formed between cells, which create signaling 'nanodomains'. Within these 'nanodomains', charged nucleotides interact with densely-packed membranes and biomolecules. While the contributions of electrostatic and steric interactions within such nanodomains are known to shape diffusion-limited reaction rates, less is understood about how these factors control the kinetics of nucleotidase activity. To quantify these factors, we utilized reaction-diffusion numerical simulations of 1) adenosine triphosphate (ATP) hydrolysis into adenosine monophosphate (AMP) and 2) AMP into adenosine (Ado) via two representative nucleotidases, CD39 and CD73. We evaluate these sequentially-coupled reactions in nanodomain geometries representative of extracellular synapses, within which we localize the nucleotidases. With this model, we find that 1) nucleotidase confinement reduces reaction rates relative to an open (bulk) system, 2) the rates of AMP and ADO formation are accelerated by restricting the diffusion of substrates away from the enzymes, and 3) nucleotidase co-localization and the presence of complementary (positive) charges to ATP enhance reaction rates, though the impact of these contributions on nucleotide pools depends on the degree to which the membrane competes for substrates. As a result, these contributions integratively control the relative concentrations and distributions of ATP and its metabolites within the junctional space. Altogether, our studies suggest that CD39 and CD73 nucleotidase activity within junctional spaces can exploit their confinement and favorable electrostatic interactions to finely control nucleotide signaling.

Electrostatic interactions are important for chaperone-client interaction in vivo

It has long been thought that chaperones are primarily attracted to their clients through the hydrophobic effect. However, in in vitro studies on the interaction between the chaperone Spy and its substrate Im7, we recently showed that long-range electrostatic interactions also play a key role. Spy functions in the periplasm of Gram-negative bacteria, which is surrounded by a permeable outer membrane. The ionic conditions in the periplasm therefore closely mimic those in the media, which allowed us to vary the ionic strength of the in vivo folding environment. Using folding biosensors that link protein folding to antibiotic resistance, we were able to monitor Spy chaperone activity in Escherichia coli in vivo as a function of media salt concentration. The chaperone activity of Spy decreased when the ionic strength of the media was increased, strongly suggesting that electrostatic forces play a vital role in the action of Spy in vivo.

BLIP-II Employs Differential Hotspot Residues To Bind Structurally Similar Staphylococcus aureus PBP2a and Class A β-Lactamases

The interaction of β-lactamase inhibitory protein II (BLIP-II) with β-lactamases serves as a model system to investigate the principles underlying protein-protein interactions. Previous studies have focused on identifying the determinants of binding affinity and specificity between BLIP-II and class A β-lactamases. However, interactions between BLIP-II and other bacterial proteins have yet to be explored. Here, we provide evidence that BLIP-II binds penicillin binding protein 2a (PBP2a) from methicillin-resistant Staphylococcus aureus (MRSA) with a K_D in the low micromolar range. In comparison to the binding constants for the potent interaction between BLIP-II and TEM-1 β-lactamase (K_D = 0.5 pM), the on-rate for BLIP-II binding PBP2a is 44 000 times slower and the off-rate is 170 times faster. Therefore, a slow association rate is a limiting factor for the potency of the interaction between BLIP-II and PBP2a. Results from alanine scanning mutagenesis of the predicted interface residues of BLIP-II indicate that charged residues on the periphery of the BLIP-II interface play a critical role for binding PBP2a, in contrast to previous findings that aromatic residues at the center of the BLIP-II interface are critical for the interaction with β-lactamases. Interestingly, many of the alanine mutants at the BLIP-II interface increase k_on for binding PBP2a, consistent with the association rate being a limiting factor for affinity. In summary, the results of the study reveal that BLIP-II binds PBP2a, although weakly compared to binding of β-lactamases, and provides insights into the different binding strategies used for these targets.

VP24-Karyopherin Alpha Binding Affinities Differ between Ebolavirus Species, Influencing Interferon Inhibition and VP24 Stability

Zaire ebolavirus (EBOV), Bundibugyo ebolavirus (BDBV), and Reston ebolavirus (RESTV) belong to the same genus but exhibit different virulence properties. VP24 protein, a structural protein present in all family members, blocks interferon (IFN) signaling and likely contributes to virulence. Inhibition of IFN signaling by EBOV VP24 (eVP24) involves its interaction with the NPI-1 subfamily of karyopherin alpha (KPNA) nuclear transporters. Here, we evaluated eVP24, BDBV VP24 (bVP24), and RESTV VP24 (rVP24) interactions with three NPI-1 subfamily KPNAs (KPNA1, KPNA5, and KPNA6). Using purified proteins, we demonstrated that each VP24 binds to each of the three NPI-1 KPNAs. bVP24, however, exhibited approximately 10-fold-lower KPNA binding affinity than either eVP24 or rVP24. Cell-based assays also indicate that bVP24 exhibits decreased KPNA interaction, decreased suppression of IFN induced gene expression, and a decreased half-life in transfected cells compared to eVP24 or rVP24. Amino acid sequence alignments between bVP24 and eVP24 also identified residues within and surrounding the previously defined eVP24-KPNA5 binding interface that decrease eVP24-KPNA affinity or bVP24-KPNA affinity. VP24 mutations that lead to reduced KPNA binding affinity also decrease IFN inhibition and shorten VP24 half-lives. These data identify novel functional differences in VP24-KPNA interaction and reveal a novel impact of the VP24-KPNA interaction on VP24 stability.

IMPORTANCE

The interaction of Ebola virus (EBOV) VP24 protein with host karyopherin alpha (KPNA) proteins blocks type I interferon (IFN) signaling, which is a central component of the host innate immune response to viral infection. Here, we quantitatively compared the interactions of VP24 proteins from EBOV and two members of the Ebolavirus genus, Bundibugyo virus (BDBV) and Reston virus (RESTV). The data reveal lower binding affinity of the BDBV VP24 (bVP24) for KPNAs and demonstrate that the interaction with KPNA modulates inhibition of IFN signaling and VP24 stability. The effect of KPNA interaction on VP24 stability is a novel functional consequence of this virus-host interaction, and the differences identified between viral species may contribute to differences in pathogenesis.

Forces Driving Chaperone Action

It is still unclear what molecular forces drive chaperone-mediated protein folding. Here, we obtain a detailed mechanistic understanding of the forces that dictate the four key steps of chaperone-client interaction: initial binding, complex stabilization, folding, and release. Contrary to the common belief that chaperones recognize unfolding intermediates by their hydrophobic nature, we discover that the model chaperone Spy uses long-range electrostatic interactions to rapidly bind to its unfolded client protein Im7. Short-range hydrophobic interactions follow, which serve to stabilize the complex. Hydrophobic collapse of the client protein then drives its folding. By burying hydrophobic residues in its core, the client's affinity to Spy decreases, which causes client release. By allowing the client to fold itself, Spy circumvents the need for client-specific folding instructions. This mechanism might help explain how chaperones can facilitate the folding of various unrelated proteins.

Subetta increases phosphorylation of insulin receptor β-subunit alone and in the presence of insulin

It has been previously shown that Subetta (a drug containing released-active forms of antibodies to the insulin receptor β-subunit and antibodies to endothelial nitric oxide synthase) stimulated insulin-induced adiponectin production by mature human adipocytes in the absence of insulin. Therefore, it was assumed that Subetta could activate the insulin receptor. To confirm this hypothesis, the capacity of Subetta to activate the insulin receptor in mature human adipocytes in the absence or presence of the insulin was investigated. Cells were incubated either with Subetta or with vehicle, or with basal medium for 3 days. Then, adipocytes were treated with water or insulin (100 nm) for 15 min. Following treatment, lysates were prepared and phosphorylation of insulin receptor β-subunits was analyzed by western blot analysis. It was shown that Subetta significantly increased (P<0.001) the ‘phosphorylated-insulin receptor β-subunit/total insulin receptor β-subunit' ratios in both the presence and the absence of insulin. These results support previously published data and indicate that Subetta could activate the insulin receptor through the effect on its β-subunits, whose conformational state is essential for insulin receptor activation. This action might serve as one of the primary mechanisms of the drug's antidiabetic effect.

Role of β-lactamase residues in a common interface for binding the structurally unrelated inhibitory proteins BLIP and BLIP-II

The β-lactamase inhibitory proteins (BLIPs) are a model system for examining molecular recognition in protein-protein interactions. BLIP and BLIP-II are structurally unrelated proteins that bind and inhibit TEM-1 β-lactamase. Both BLIPs share a common binding interface on TEM-1 and make contacts with many of the same TEM-1 surface residues. BLIP-II, however, binds TEM-1 over 150-fold tighter than BLIP despite the fact that it has fewer contact residues and a smaller binding interface. The role of eleven TEM-1 amino acid residues that contact both BLIP and BLIP-II was examined by alanine mutagenesis and determination of the association (k on) and dissociation (k off) rate constants for binding each partner. The substitutions had little impact on association rates and resulted in a wide range of dissociation rates as previously observed for substitutions on the BLIP side of the interface. The substitutions also had less effect on binding affinity for BLIP than BLIP-II. This is consistent with the high affinity and small binding interface of the TEM-1-BLIP-II complex, which predicts per residue contributions should be higher for TEM-1 binding to BLIP-II versus BLIP. Two TEM-1 residues (E104 and M129) were found to be hotspots for binding BLIP while five (L102, Y105, P107, K111, and M129) are hotspots for binding BLIP-II with only M129 as a common hotspot for both. Thus, although the same TEM-1 surface binds to both BLIP and BLIP-II, the distribution of binding energy on the surface is different for the two target proteins, that is, different binding strategies are employed.

pH-Dependent aggregation and disaggregation of native β-lactoglobulin in low salt

The aggregation of β-lactoglobulin (BLG) near its isoelectric point was studied as a function of ionic strength and pH. We compared the behavior of native BLG with those of its two isoforms, BLG-A and BLG-B, and with that of a protein with a very similar pI, bovine serum albumin (BSA). Rates of aggregation were obtained through a highly precise and convenient pH/turbidimetric titration that measures transmittance to ±0.05 %T. A comparison of BLG and BSA suggests that the difference between pHmax (the pH of the maximum aggregation rate) and pI is systematically related to the nature of protein charge asymmetry, as further supported by the effect of localized charge density on the dramatically different aggregation rates of the two BLG isoforms. Kinetic measurements including very short time periods show well-differentiated first and second steps. BLG was analyzed by light scattering under conditions corresponding to maxima in the first and second steps. Dynamic light scattering (DLS) was used to monitor the kinetics, and static light scattering (SLS) was used to evaluate the aggregate structure fractal dimensions at different quench points. The rate of the first step is relatively symmetrical around pHmax and is attributed to the local charges within the negative domain of the free protein. In contrast, the remarkably linear pH dependence of the second step is related to the uniform reduction in global protein charge with increasing pH below pI, accompanied by an attractive force due to surface charge fluctuations.

On the binding affinity of macromolecular interactions: daring to ask why proteins interact

Interactions between proteins are orchestrated in a precise and time-dependent manner, underlying cellular function. The binding affinity, defined as the strength of these interactions, is translated into physico-chemical terms in the dissociation constant (K_d), the latter being an experimental measure that determines whether an interaction will be formed in solution or not. Predicting binding affinity from structural models has been a matter of active research for more than 40 years because of its fundamental role in drug development. However, all available approaches are incapable of predicting the binding affinity of protein–protein complexes from coordinates alone. Here, we examine both theoretical and experimental limitations that complicate the derivation of structure–affinity relationships. Most work so far has concentrated on binary interactions. Systems of increased complexity are far from being understood. The main physico-chemical measure that relates to binding affinity is the buried surface area, but it does not hold for flexible complexes. For the latter, there must be a significant entropic contribution that will have to be approximated in the future. We foresee that any theoretical modelling of these interactions will have to follow an integrative approach considering the biology, chemistry and physics that underlie protein–protein recognition.

Complete thermodynamic and kinetic characterization of the isomer-specific interaction between Pin1-WW domain and the amyloid precursor protein cytoplasmic tail phosphorylated at Thr668

Peptidyl prolyl cis-trans isomerization acts as an effective molecular timer that plays significant roles in biological and pathological processes. Enzymes such as Pin1 catalyze cis-trans isomerization, accelerating the otherwise slow isomerization rate into time scales relevant for cellular signaling. Here we have combined NMR line shape analysis, fluorescence spectroscopy, and isothermal titration calorimetry to determine the kinetic and thermodynamic parameters describing the trans-specific interaction between the binding domain of Pin1 (WW domain) and a key cis-trans molecular switch in the amyloid precursor protein cytoplasmic tail. A three-state model, in which the cis-trans isomerization equilibrium is coupled to the binding equilibrium through the trans isomer, was found to fit the data well. The trans isomer binds the WW domain with ∼22 μM affinity via very fast association (approaching the diffusion limit) and dissociation rates. The common structural and electrostatic characteristics of Pin1 substrates, which contain a phosphorylated serine/threonine-proline motif, suggest that very rapid binding kinetics are a general feature of Pin1 interactions with other substrates. The fast binding kinetics of the WW domain allows rapid response of Pin1 to the dynamic events of phosphorylation and dephosphorylation in the cell that alter the relative populations of diverse Pin1 substrates. Furthermore, our results also highlight the vastly different rates at which slow uncatalyzed cis-trans isomerization and fast isomer-specific binding events occur. These results, along with the experimental methods presented herein, should guide future experiments aimed at the thermodynamic and kinetic characterization of cis-trans molecular switches and isomer-specific interactions involved in various biological processes.

Bovine β-lactoglobulin is dimeric under imitative physiological conditions: dissociation equilibrium and rate constants over the pH range of 2.5-7.5

The oligomerization of β-lactoglobulin (βLg) has been studied extensively, but with somewhat contradictory results. Using analytical ultracentrifugation in both sedimentation equilibrium and sedimentation velocity modes, we studied the oligomerization of βLg variants A and B over a pH range of 2.5-7.5 in 100 mM NaCl at 25°C. For the first time, to our knowledge, we were able to estimate rate constants (k(off)) for βLg dimer dissociation. At pH 2.5 k(off) is low (0.008 and 0.009 s(-1)), but at higher pH (6.5 and 7.5) k(off) is considerably greater (>0.1 s(-1)). We analyzed the sedimentation velocity data using the van Holde-Weischet method, and the results were consistent with a monomer-dimer reversible self-association at pH 2.5, 3.5, 6.5, and 7.5. Dimer dissociation constants K(D)(2-1) fell close to or within the protein concentration range of ∼5 to ∼45 μM, and at ∼45 μM the dimer predominated. No species larger than the dimer could be detected. The K(D)(2-1) increased as |pH-pI| increased, indicating that the hydrophobic effect is the major factor stabilizing the dimer, and suggesting that, especially at low pH, electrostatic repulsion destabilizes the dimer. Therefore, through Poisson-Boltzmann calculations, we determined the electrostatic dimerization energy and the ionic charge distribution as a function of ionic strength at pH above (pH 7.5) and below (pH 2.5) the isoelectric point (pI∼5.3). We propose a mechanism for dimer stabilization whereby the added ionic species screen and neutralize charges in the vicinity of the dimer interface. The electrostatic forces of the ion cloud surrounding βLg play a key role in the thermodynamics and kinetics of dimer association/dissociation.

Temperature-dependent antibody kinetics as a tool in antibody lead selection

Antibody-antigen interactions can principally be classified into three different temperature-dependent kinetic rate profiles. The affinity K (D) can persist, decrease, or increase in the temperature gradient. Today, the impact of temperature-dependent antibody kinetics is recognized, especially as part of the development of best in class monoclonal antibodies. Here, a robust surface plasmon resonance-based protocol is presented, which describes a sensitive temperature-dependent kinetic measurement and evaluation method.

Clustering of HIV-1 Subtypes Based on gp120 V3 Loop electrostatic properties

BACKGROUND

The V3 loop of the glycoprotein gp120 of HIV-1 plays an important role in viral entry into cells by utilizing as coreceptor CCR5 or CXCR4, and is implicated in the phenotypic tropisms of HIV viruses. It has been hypothesized that the interaction between the V3 loop and CCR5 or CXCR4 is mediated by electrostatics. We have performed hierarchical clustering analysis of the spatial distributions of electrostatic potentials and charges of V3 loop structures containing consensus sequences of HIV-1 subtypes.

RESULTS

Although the majority of consensus sequences have a net charge of +3, the spatial distribution of their electrostatic potentials and charges may be a discriminating factor for binding and infectivity. This is demonstrated by the formation of several small subclusters, within major clusters, which indicates common origin but distinct spatial details of electrostatic properties. Some of this information may be present, in a coarse manner, in clustering of sequences, but the spatial details are largely lost. We show the effect of ionic strength on clustering of electrostatic potentials, information that is not present in clustering of charges or sequences. We also make correlations between clustering of electrostatic potentials and net charge, coreceptor selectivity, global prevalence, and geographic distribution. Finally, we interpret coreceptor selectivity based on the N6X7T8|S8X9 sequence glycosylation motif, the specific positive charge location according to the 11/24/25 rule, and the overall charge and electrostatic potential distribution.

CONCLUSIONS

We propose that in addition to the sequence and the net charge of the V3 loop of each subtype, the spatial distributions of electrostatic potentials and charges may also be important factors for receptor recognition and binding and subsequent viral entry into cells. This implies that the overall electrostatic potential is responsible for long-range recognition of the V3 loop with coreceptors CCR5/CXCR4, whereas the charge distribution contributes to the specific short-range interactions responsible for the formation of the bound complex. We also propose a scheme for coreceptor selectivity based on the sequence glycosylation motif, the 11/24/25 rule, and net charge.

Quantitative, directional measurement of electric field heterogeneity in the active site of ketosteroid isomerase

Understanding the electrostatic forces and features within highly heterogeneous, anisotropic, and chemically complex enzyme active sites and their connection to biological catalysis remains a longstanding challenge, in part due to the paucity of incisive experimental probes of electrostatic properties within proteins. To quantitatively assess the landscape of electrostatic fields at discrete locations and orientations within an enzyme active site, we have incorporated site-specific thiocyanate vibrational probes into multiple positions within bacterial ketosteroid isomerase. A battery of X-ray crystallographic, vibrational Stark spectroscopy, and NMR studies revealed electrostatic field heterogeneity of 8 MV/cm between active site probe locations and widely differing sensitivities of discrete probes to common electrostatic perturbations from mutation, ligand binding, and pH changes. Electrostatic calculations based on active site ionization states assigned by literature precedent and computational pK(a) prediction were unable to quantitatively account for the observed vibrational band shifts. However, electrostatic models of the D40N mutant gave qualitative agreement with the observed vibrational effects when an unusual ionization of an active site tyrosine with a pK(a) near 7 was included. UV-absorbance and (13)C NMR experiments confirmed the presence of a tyrosinate in the active site, in agreement with electrostatic models. This work provides the most direct measure of the heterogeneous and anisotropic nature of the electrostatic environment within an enzyme active site, and these measurements provide incisive benchmarks for further developing accurate computational models and a foundation for future tests of electrostatics in enzymatic catalysis.

On the role of electrostatics in protein-protein interactions

The role of electrostatics in protein-protein interactions and binding is reviewed in this paper. A brief outline of the computational modeling, in the framework of continuum electrostatics, is presented and the basic electrostatic effects occurring upon the formation of the complex are discussed. The effect of the salt concentration and pH of the water phase on protein-protein binding free energy is demonstrated which indicates that the increase of the salt concentration tends to weaken the binding, an observation that is attributed to the optimization of the charge-charge interactions across the interface. It is pointed out that the pH-optimum (pH of optimal binding affinity) varies among the protein-protein complexes, and perhaps is a result of their adaptation to particular subcellular compartments. The similarities and differences between hetero- and homo-complexes are outlined and discussed with respect to the binding mode and charge complementarity.

Triathlon for energy functions: who is the winner for design of protein-protein interactions?

Computational prediction of stabilizing mutations into monomeric proteins has become an almost ordinary task. Yet, computational stabilization of protein–protein complexes remains a challenge. Design of protein–protein interactions (PPIs) is impeded by the absence of an energy function that could reliably reproduce all favorable interactions between the binding partners. In this work, we present three energy functions: one function that was trained on monomeric proteins, while the other two were optimized by different techniques to predict side-chain conformations in a dataset of PPIs. The performances of these energy functions are evaluated in three different tasks related to design of PPIs: predicting side-chain conformations in PPIs, recovering native binding-interface sequences, and predicting changes in free energy of binding due to mutations. Our findings show that both functions optimized on side-chain repacking in PPIs are more suitable for PPI design compared to the function trained on monomeric proteins. Yet, no function performs best at all three tasks. Comparison of the three energy functions and their performances revealed that (1) burial of polar atoms should not be penalized significantly in PPI design as in single-protein design and (2) contribution of electrostatic interactions should be increased several-fold when switching from single-protein to PPI design. In addition, the use of a softer van der Waals potential is beneficial in cases when backbone flexibility is important. All things considered, we define an energy function that captures most of the nuances of the binding energetics and hence, should be used in future for design of PPIs.

Automated computational framework for the analysis of electrostatic similarities of proteins

Charge plays an important role in protein-protein interactions. In the case of excessively charged proteins, their electrostatic potentials contribute to the processes of recognition and binding with other proteins or ligands. We present an automated computational framework for determining the contribution of each charged amino acid to the electrostatic properties of proteins, at atomic resolution level. This framework involves computational alanine scans, calculation of Poisson-Boltzmann electrostatic potentials, calculation of electrostatic similarity distances (ESDs), hierarchical clustering analysis of ESDs, calculation of solvation free energies of association, and visualization of the spatial distributions of electrostatic potentials. The framework is useful to classify families of mutants with similar electrostatic properties and to compare them with the parent proteins in the complex. The alanine scan mutants introduce perturbations in the local electrostatic properties of the proteins and aim in delineating the contribution of each mutated amino acid in the spatial distribution of electrostatic potential, and in biological function when electrostatics is a dominant contributing factor in protein-protein interactions. The framework can be used to design new proteins with tailored electrostatic properties, such as immune system regulators, inhibitors, and vaccines, and in guiding experimental studies. We present an example for the interaction of the immune system protein C3d (the d-fragment of complement protein C3) with its receptor CR2, and we discuss our data in view of a binding site controversy.

Identification of a β-lactamase inhibitory protein variant that is a potent inhibitor of Staphylococcus PC1 β-lactamase

β-Lactamase inhibitory protein (BLIP) binds and inhibits a diverse collection of class A β-lactamases. Widespread resistance to β-lactam antibiotics currently limits the treatment strategies for Staphylococcus infections. The goals of this study were to determine the binding affinity of BLIP for Staphylococcus aureus PC1 β-lactamase and to identify mutants that alter binding affinity. The BLIP inhibition constant (K(i)) for PC1 β-lactamase was measured at 350 nM, and isothermal titration calorimetry experiments indicated a binding constant (K(d)) of 380 nM. Twenty-three residue positions in BLIP that contact β-lactamase were randomized, and phage display was used to sort the libraries for tight binders to immobilized PC1 β-lactamase. The BLIP(K74G) mutant was the dominant clone selected, and it was found to inhibit the PC1 β-lactamase with a K(i) of 42 nM, while calorimetry indicated a K(d) of 26 nM. Molecular modeling studies suggested that BLIP binds weakly to the PC1 β-lactamase due to the presence of alanine at position 104 of PC1. This position is occupied by glutamate in the TEM-1 enzyme, where it forms a salt bridge with the BLIP residue Lys74 that is important for the stability of the complex. This hypothesis was confirmed by showing that the PC1(A104E) enzyme binds BLIP with 15-fold greater affinity than wild-type PC1 β-lactamase. Kinetic measurements indicated similar association rates for all complexes with variation in affinity due to altered dissociation rate constants, suggesting that changes in short-range interactions are responsible for the altered binding properties of the mutants.

A biosensor study indicating that entropy, electrostatics, and receptor glycosylation drive the binding interaction between interleukin-7 and its receptor

The interaction between interleukin-7 (IL-7) and its α-receptor, IL-7Rα, plays fundamental roles in the development, survival, and homeostasis of B- and T-cells. N-Linked glycosylation of human IL-7Rα enhances its binding affinity for human IL-7 300-fold versus that of the nonglycosylated receptor through an allosteric mechanism. The N-glycans of IL-7Rα do not participate directly in the binding interface with IL-7. This biophysical study involves dissection of the properties of binding of IL-7 to both nonglycosylated and glycosylated forms of the IL-7Rα extracellular domain (ECD) as functions of salt, pH, and temperature using surface plasmon resonance (SPR) spectroscopy. Interactions of IL-7 with both IL-7Rα variants display weaker binding affinities with increasing salt concentrations primarily reflected by changes in the first on rates of a two-step reaction pathway. The electrostatic parameter of the IL-7-IL-7Rα interaction is not driven by complementary charge interactions through residues at the binding interface or N-glycan composition of IL-7Rα, but presumably by favorable global charges of the two proteins. van't Hoff analysis indicates both IL-7-IL-7Rα interactions are driven by large favorable entropy changes and smaller unfavorable (nonglycosylated complex) and favorable (glycosylated complex) enthalpy changes. Eyring analysis of the IL-7-IL-7Rα interactions reveals different reaction pathways and barriers for the transition-state thermodynamics with the enthalpy and entropy changes of IL-7 binding to nonglycosylated and glycosylated IL-7Rα. There were no discernible heat capacity changes for the equilibrium or transition-state binding thermodynamics of the IL-7-IL-7Rα interactions. The results suggest that the unbound nonglycosylated IL-7Rα samples an extensive conformational landscape relative to the unbound glycosylated IL-7Rα, potentially explaining the switch from a "conformationally controlled" reaction (k(1) ∼ 10(2) M(-1) s(-1)) for the nonglycosylated interaction to a "diffusion-controlled" reaction (k(1) ∼ 10(6) M(-1) s(-1)) for the glycosylated interaction. Thus, a large favorable entropy change, a global favorable electrostatic component, and glycosylation of the receptor, albeit not at the interface, contribute significantly to the interaction between IL-7 and the IL-7Rα ECD.

Glycosaminoglycans as polyelectrolytes

One of the barriers to understanding structure-property relations for glycosaminoglycans has been the lack of constructive interplay between the principles and methodologies of the life sciences (molecular biology, biochemistry and cell biology) and the physical sciences, particularly in the field of polyelectrolytes. To address this, we first review the similarities and differences between the physicochemical properties of GAGs and other statistical chain polyelectrolytes of both natural and abioitic origin. Since the biofunctionality and regulation of the structures of GAGs is intimately connected with interactions with their cognate proteins, we particularly compare and contrast aspects of protein binding, i.e. effects of both GAGs and other polyelectrolytes on protein stability, protein aggregation and phase behavior. The protein binding affinities and their dependences on pH and ionic strength for the two groups are discussed not only in terms of observable differences, but also with regard to contrasting descriptions of the bound state and the role of electrostatics. We conclude that early studies of the heparin-Antithromin system, proceeding to a large extent through the methods and models of protein chemistry and drug discovery, established not only many enabling precedents but also constraining paradigms. Current studies on heparan sulfate and chondroitin sulfate seem to reflect a more ecumenical view likely to be more compatible with concepts from physical and polymer chemistry.

Assembly dynamics and stability of the pneumococcal epsilon zeta antitoxin toxin (PezAT) system from Streptococcus pneumoniae

The pneumococcal epsilon zeta antitoxin toxin (PezAT) system is a chromosomally encoded, class II toxin antitoxin system from the human pathogen Streptococcus pneumnoniae. Neutralization of the bacteriotoxic protein PezT is carried out by complex formation with its cognate antitoxin PezA. Here we study the stability of the inhibitory complex in vivo and in vitro. We found that toxin release is impeded in Escherichia coli and Bacillus subtilis due to the proteolytic resistance of PezA once bound to PezT. These findings are supported by in vitro experiments demonstrating a strong thermodynamic stabilization of both proteins upon binding. A detailed kinetic analysis of PezAT assembly revealed that these particular features of PezAT are based on a strong, electrostatically guided binding mechanism leading to a stable toxin antitoxin complex with femtomolar affinity. Our data show that PezAT complex formation is distinct to all other conventional toxin antitoxin modules and a controlled mode of toxin release is required for activation.

pH-dependent association of proteins. The test case of monoclonal antibody HyHEL-5 and its antigen hen egg white lysozyme

We describe a method for determining diffusion-controlled rate constants for protein-protein association that explicitly includes the solution pH. The method combines the transient-complex theory for computing electrostatically enhanced association rates with an approach based on a rigorous thermodynamic cycle and partition functions for energy levels characterizing protonation states of associating proteins and their complexes. To test our method, we determine the pH-dependent kinetics of association of the HyHEL-5 antibody with its antigen hen egg white lysozyme. It was shown experimentally that their association rate constant depends on pH, increasing linearly in the pH range 6-8 and saturating or even exhibiting a flat maximum in the pH range 8-10. The presented methodology leads to a qualitative agreement with the experimental data. Our approach allows one to study diffusion-controlled protein-protein association under different pH conditions by taking into account the ensembles of protonation states rather than just the most probable protonation state of each protein.

The effect of different force applications on the protein-protein complex Barnase-Barstar

Steered molecular dynamics simulations are a tool to examine the energy landscape of protein-protein complexes by applying external forces. Here, we analyze the influence of the velocity and geometry of the probing forces on a protein complex using this tool. With steered molecular dynamics, we probe the stability of the protein-protein complex Barnase-Barstar. The individual proteins are mechanically labile. The Barnase-Barstar binding site is more stable than the folds of the individual proteins. By using different force protocols, we observe a variety of responses of the system to the applied tension.

Structure of Escherichia coli Hfq bound to polyriboadenylate RNA

Hfq is a small, highly abundant hexameric protein that is found in many bacteria and plays a critical role in mRNA expression and RNA stability. As an "RNA chaperone," Hfq binds AU-rich sequences and facilitates the trans annealing of small RNAs (sRNAs) to their target mRNAs, typically resulting in the down-regulation of gene expression. Hfq also plays a key role in bacterial RNA decay by binding tightly to polyadenylate [poly(A)] tracts. The structural mechanism by which Hfq recognizes and binds poly(A) is unknown. Here, we report the crystal structure of Escherichia coli Hfq bound to the poly(A) RNA, A(15). The structure reveals a unique RNA binding mechanism. Unlike uridine-containing sequences, which bind to the "proximal" face, the poly(A) tract binds to the "distal" face of Hfq using 6 tripartite binding motifs. Each motif consists of an adenosine specificity site (A site), which is effected by peptide backbone hydrogen bonds, a purine nucleotide selectivity site (R site), and a sequence-nondiscriminating RNA entrance/exit site (E site). The resulting implication that Hfq can bind poly(A-R-N) triplets, where R is a purine nucleotide and N is any nucleotide, was confirmed by binding studies. Indeed, Hfq bound to the oligoribonucleotides (AGG)(8), (AGC)(8), and the shorter (A-R-N)(4) sequence, AACAACAAGAAG, with nanomolar affinities. The abundance of (A-R-N)(4) and (A-R-N)(5) triplet repeats in the E. coli genome suggests additional RNA targets for Hfq. Further, the structure provides insight into Hfq-mediated sRNA-mRNA annealing and the role of Hfq in RNA decay.

Improved binding of raf to Ras.GDP is correlated with biological activity

The GTP-binding protein Ras plays a central role in the regulation of various cellular processes, acting as a molecular switch that triggers signaling cascades. Only Ras bound to GTP is able to interact strongly with effector proteins like Raf kinase, phosphatidylinositol 3-kinase, and RalGDS, whereas in the GDP-bound state, the stability of the complex is strongly decreased, and signaling is interrupted. To determine whether this process is only controlled by the stability of the complex, we used computer-aided protein design to improve the interaction between Ras and effector. We challenged the Ras.Raf complex in this study because Raf among all effectors shows the highest Ras affinity and the fastest association kinetics. The proposed mutations were characterized as to their changes in dynamics and binding strength. We demonstrate that Ras-Raf interaction can only be improved at the cost of a loss in specificity of Ras.GTP versus Ras.GDP. As shown by NMR spectroscopy, the Raf mutation A85K leads to a shift of Ras switch I in the GTP-bound as well as in the GDP-bound state, thereby increasing the complex stability. In a luciferase-based reporter gene assay, Raf A85K is associated with higher signaling activity, which appears to be a mere matter of Ras-Raf affinity.

Cross-species analysis of the glycolytic pathway by comparison of molecular interaction fields

The electrostatic potential of an enzyme is a key determinant of its substrate interactions and catalytic turnover. Here we invoke comparative analysis of protein electrostatic potentials, along with sequence and structural analysis, to classify and characterize all the enzymes in an entire pathway across a set of different organisms. The electrostatic potentials of the enzymes from the glycolytic pathway of 11 eukaryotes were analyzed by qPIPSA (quantitative protein interaction property similarity analysis). The comparison allows the functional assignment of neuron-specific isoforms of triosephosphate isomerase from zebrafish, the identification of unusual protein surface interaction properties of the mosquito glucose-6-phosphate isomerase and the functional annotation of ATP-dependent phosphofructokinases and cofactor-dependent phosphoglycerate mutases from plants. We here show that plants possess two parallel pathways to convert glucose. One is similar to glycolysis in humans, the other is specialized to let plants adapt to their environmental conditions. We use differences in electrostatic potentials to estimate kinetic parameters for the triosephosphate isomerases from nine species for which published parameters are not available. Along the core glycolytic pathway, phosphoglycerate mutase displays the most conserved electrostatic potential. The largest cross-species variations are found for glucose-6-phosphate isomerase, enolase and fructose-1,6-bisphosphate aldolase. The extent of conservation of electrostatic potentials along the pathway is consistent with the absence of a single rate-limiting step in glycolysis.

Fruitful and futile encounters along the association reaction between proteins

The association reaction between pairs of proteins proceeds through an encounter complex that develops into the final complex. Here, we combined Brownian dynamics simulations with experimental studies to analyze the structures of the encounter complexes along the association reaction between TEM1-beta-lactamase and its inhibitor, beta-lactamase-inhibitor protein. The encounter complex can be considered as an ensemble of short-lived low free-energy states that are stabilized primarily by electrostatic forces and desolvation. For the wild-type, the simulation showed two main encounter regions located outside the physical binding site. One of these regions was located near the experimentally determined transition state. To validate whether these encounters are fruitful or futile, we examined three groups of mutations that altered the encounter. The first group consisted of mutations that increased the experimental rate of association through electrostatic optimization. This resulted in an increase in the size of the encounter region located near the experimentally determined transition state, as well as a decrease in the energy of this region and an increase in the number of successful trajectories (i.e., encounters that develop into complex). A second group of mutations was specifically designed to either increase or decrease the size and energy of the second encounter complex, but either way it did not affect k(on). A third group of mutations consisted of residues that increased k(on) without significantly affecting the encounter complexes. These results indicate that the size and energy of the encounter regions are only two of several parameters that lead to fruitful association, and that electrostatic optimization is a major driving force in fast association.

Binding mechanism of an SH3 domain studied by NMR and ITC

Complexes between Src-homology 3 domains and proline-rich target peptides can have lifetimes on the order of milliseconds, making them too short-lived for kinetic characterization by conventional methods. Nuclear magnetic resonance (NMR) dynamics experiments are ideally suited to study such rapid binding equilibria, and additionally provide information on partly bound intermediate states. We used NMR together with isothermal titration calorimetry (ITC) to characterize the interaction of the SH3 domain from the Fyn tyrosine kinase with a 12-residue peptide at temperatures between 10 and 50 degrees C. NMR data at all temperatures are consistent with an effectively two-state binding reaction, such that any intermediates are either very weakly populated or exchange extremely rapidly with the free or bound forms. Dissociation rate constants, determined by CPMG and ZZ-exchange NMR experiments, range from k(off)(10 degrees C) = 4.5 s(-1) to k(off)(50 degrees C) = 331 s(-1). ITC data at all temperatures follow a simple two-state interaction model. Binding is favored enthalpically, with a dissociation enthalpy, DeltaH(D)(30 degrees C) = 15.4 kcal mol(-1), and disfavored entropically, with a dissociation entropy, DeltaS(D)(30 degrees C) = 20.0 cal mol(-1) K(-1). The free protein and peptide have significantly higher heat capacity than the bound complex, DeltaC(p) = 352 cal mol(-1) K(-1), which is consistent with the largely hydrophobic character of the binding interface. An Eyring plot of k(off) values gives an activation enthalpy of dissociation, DeltaH(D)(double dagger)(30 degrees C) = 19.3 kcal mol(-1) and exhibits slight curvature consistent with the ITC-derived value of DeltaC(p). The curvature suggests that nonpolar residues of the hydrophobic interface are solvated in the transition state for dissociation. Association rate constants were calculated as k(on) = k(off)/K(D), and range from k(on)(10 degrees C) = 1.03 x 10(8) M(-1) s(-1) to k(on)(50 degrees C) = 2.0 x 10(8) M(-1) s(-1), with an apparent activation enthalpy, DeltaH(A)(double dagger) = 3.4 kcal mol(-1). Both the magnitudes and temperature dependence of k(on) values are consistent with a diffusion-limited association mechanism. The combination of NMR and ITC data sheds light on how the Fyn tyrosine kinase is activated by binding to proline-rich targets, and represents a powerful approach for characterizing transient protein/ligand interactions.

Electrostatic interactions of Hsp-organizing protein tetratricopeptide domains with Hsp70 and Hsp90: computational analysis and protein engineering

The Hsp-organizing protein (HOP) binds to the C termini of the chaperones Hsp70 and Hsp90, thus bringing them together so that substrate proteins can be passed from Hsp70 to Hsp90. Because Hsp90 is essential for the correct folding and maturation of many oncogenic proteins, it has become a significant target for anti-cancer drug design. HOP binds to Hsp70 and Hsp90 via two independent tetratricopeptide (TPR) domains, TPR1 and TPR2A, respectively. We have analyzed ligand binding using Poisson-Boltzmann continuum electrostatic calculations, free energy perturbation, molecular dynamics simulations, and site-directed mutagenesis to delineate the contribution of different interactions to the affinity and specificity of the TPR-peptide interactions. We found that continuum electrostatic calculations could be used to guide protein design by removing unfavorable interactions to increase binding affinity, with an 80-fold increase in affinity for TPR2A. Contributions at buried charged residues, however, were better predicted by free energy perturbation calculations. We suggest using a combination of the two approaches for increasing the accuracy of results, with free energy perturbation calculations used only at selected buried residues of the ligand binding pocket. Finally we present the crystal structure of TPR2A in complex with its non-cognate Hsp70 ligand, which provides insight on the origins of specificity in TPR domain-peptide recognition.

Docking of antizyme to ornithine decarboxylase and antizyme inhibitor using experimental mutant and double-mutant cycle data

Antizyme (Az) is a highly conserved key regulatory protein bearing a major role in regulating polyamine levels in the cell. It has the ability to bind and inhibit ornithine decarboxylase (ODC), targeting it for degradation. Az inhibitor (AzI) impairs the activity of Az. In this study, we mapped the binding sites of ODC and AzI on Az using Ala scan mutagenesis and generated models of the two complexes by constrained computational docking. In order to scan a large number of mutants in a short time, we developed a workflow combining high-throughput mutagenesis, small-scale parallel partial purification of His-tagged proteins and their immobilization on a tris-nitrilotriacetic-acid-coated surface plasmon resonance chip. This combination of techniques resulted in a significant reduction in time for production and measurement of large numbers of mutant proteins. The data-driven docking results suggest that both proteins occupy the same binding site on Az, with Az binding within a large groove in AzI and ODC. However, single-mutant data provide information concerning the location of the binding sites only, not on their relative orientations. Therefore, we generated a large number of double-mutant cycles between residues on Az and ODC and used the resulting interaction energies to restrict docking. The model of the complex is well defined and accounts for the mutant data generated here, and previously determined biochemical data for this system. Insights on the structure and function of the complexes, as well as general aspects of the method, are discussed.

Co-chaperonin GroES as a modulator of proteasomal activity

The proteasome has a crucial part in the degradation of normal, damaged, mutant or misfolded proteins within both the ubiquitin ATP-dependent and the ubiquitin ATP-independent pathways. Proteasome-mediated proteolysis is modulated by diverse factors, and in this regard, chaperonins have been attracting great interest. The investigation on the role of a co-chaperonin, namely GroES, in the modulation of proteasomal activity was the focus of this work. Our study reports on an analytical approach based on combined fluorimetric, chromatographic (applied to the enzymatic activity evaluation), surface plasmon resonance techniques and molecular modelling, addressed to the assessment and characterization of the interaction. Globally, we described a high affinity interaction between GroES and two different 20 S (immuno- and constitutive) proteasomes, uncovering new scenarios on their possible physio-pathological role, specifically on the ability of proteasomes to interact both with unfolding and folding- assisting macromolecules.

Fast predictions of thermodynamics and kinetics of protein-protein recognition from structures: from molecular design to systems biology

The increasing call for an overall picture of the interactions between the components of a biological system that give rise to the observed function is often summarized by the expression systems biology. Both the interpretative and predictive capabilities of holistic models of biochemical systems, however, depend to a large extent on the level of physico-chemical knowledge of the individual molecular interactions making up the network. This review is focused on the structure-based quantitative characterization of protein-protein interactions, ubiquitous in any biochemical pathway. Recently developed, fast and effective computational methods are reviewed, which allow the assessment of kinetic and thermodynamic features of the association-dissociation processes of protein complexes, both in water soluble and membrane environments. The performance and the accuracy of fast and semi-empirical structure-based methods have reached comparable levels with respect to the classical and more elegant molecular simulations. Nevertheless, the broad accessibility and lower computational cost provide the former methods with the advantageous possibility to perform systems-level analyses including extensive in silico mutagenesis screenings and large-scale structural predictions of multiprotein complexes.

Interaction of onconase with the human ribonuclease inhibitor protein

One of the tightest known protein-protein interactions in biology is that between members of the ribonuclease A superfamily and the ribonuclease inhibitor protein (RI). Some members of this superfamily are able to kill cancer cells, and the ability to evade RI is a major determinant of whether a ribonuclease will be cytotoxic. The archetypal cytotoxic ribonuclease, onconase (ONC), is in late-stage clinical trials for the treatment of malignant mesothelioma. We present here the first measurement of the inhibition of the ribonucleolytic activity of ONC by RI. This inhibition occurs with K(i)=0.15muM in a solution of low salt concentration.

Conformational stability and folding mechanisms of dimeric proteins

The folding of multisubunit proteins is of tremendous biological significance since the large majority of proteins exist as protein-protein complexes. Extensive experimental and computational studies have provided fundamental insights into the principles of folding of small monomeric proteins. Recently, important advances have been made in extending folding studies to multisubunit proteins, in particular homodimeric proteins. This review summarizes the equilibrium and kinetic theory and models underlying the quantitative analysis of dimeric protein folding using chemical denaturation, as well as the experimental results that have been obtained. Although various principles identified for monomer folding also apply to the folding of dimeric proteins, the effects of subunit association can manifest in complex ways, and are frequently overlooked. Changes in molecularity typically give rise to very different overall folding behaviour than is observed for monomeric proteins. The results obtained for dimers have provided key insights pertinent to understanding biological assembly and regulation of multisubunit proteins. These advances have set the stage for future advances in folding involving protein-protein interactions for natural multisubunit proteins and unnatural assemblies involved in disease.

Dual GPCR and GAG mimicry by the M3 chemokine decoy receptor

Viruses have evolved a myriad of evasion strategies focused on undermining chemokine-mediated immune surveillance, exemplified by the mouse γ-herpesvirus 68 M3 decoy receptor. Crystal structures of M3 in complex with C chemokine ligand 1/lymphotactin and CC chemokine ligand 2/monocyte chemoattractant protein 1 reveal that invariant chemokine features associated with G protein–coupled receptor binding are primarily recognized by the decoy C-terminal domain, whereas the N-terminal domain (NTD) reconfigures to engage divergent basic residue clusters on the surface of chemokines. Favorable electrostatic forces dramatically enhance the association kinetics of chemokine binding by M3, with a primary role ascribed to acidic NTD regions that effectively mimic glycosaminoglycan interactions. Thus, M3 employs two distinct mechanisms of chemical imitation to potently sequester chemokines, thereby inhibiting chemokine receptor binding events as well as the formation of chemotactic gradients necessary for directed leukocyte trafficking.

Linkage between substrate recognition and catalysis during cleavage of sarcin/ricin loop RNA by restrictocin

Restrictocin is a site-specific endoribonuclease that inactivates ribosomes by cleaving the sarcin/ricin loop (SRL) of 23S-28S rRNA. Here we present a kinetic and thermodynamic analysis of the SRL cleavage reaction based on monitoring the cleavage of RNA oligonucleotides (2-27-mers). Restrictocin binds to a 27-mer SRL model substrate (designated wild-type SRL) via electrostatic interactions to form a nonspecific ground state complex E:S. At pH 6.7, physical steps govern the reaction rate: the wild-type substrate reacts at a partially diffusion-limited rate, and a faster-reacting SRL, containing a 3'-sulfur atom at the scissile phosphate, reacts at a fully diffusion-limited rate (k2/K1/2 = 1.1 x 10(9) M-1 s-1). At pH 7.4, the chemical step apparently limits the SRL cleavage rate. After the nonspecific binding step, restrictocin recognizes the SRL structure, which imparts 4.3 kcal/mol transition state stabilization relative to a single-stranded RNA. The two conserved SRL modules, bulged-G motif and GAGA tetraloop, contribute at least 2.4 and 1.9 kcal/mol, respectively, to the recognition. These findings suggest a model of SRL recognition in which restrictocin contacts the GAGA tetraloop and the bulged guanosine of the bulged-G motif to progress from the nonspecific ground state complex (E:S) to the higher-energy-specific complex (E.S) en route to the chemical transition state. Comparison of restrictocin with other ribonucleases revealed that restrictocin exhibits a 10(3)-10(6)-fold smaller ribonuclease activity against single-stranded RNA than do the restrictocin homologues, non-structure-specific ribonucleases T1 and U2. Together, these findings show how structural features of the SRL substrate facilitate catalysis and provide a mechanism for distinguishing between cognate and noncognate RNA.

qPIPSA: relating enzymatic kinetic parameters and interaction fields

Background

The simulation of metabolic networks in quantitative systems biology requires the assignment of enzymatic kinetic parameters. Experimentally determined values are often not available and therefore computational methods to estimate these parameters are needed. It is possible to use the three-dimensional structure of an enzyme to perform simulations of a reaction and derive kinetic parameters. However, this is computationally demanding and requires detailed knowledge of the enzyme mechanism. We have therefore sought to develop a general, simple and computationally efficient procedure to relate protein structural information to enzymatic kinetic parameters that allows consistency between the kinetic and structural information to be checked and estimation of kinetic constants for structurally and mechanistically similar enzymes.

Results

We describe qPIPSA: quantitative Protein Interaction Property Similarity Analysis. In this analysis, molecular interaction fields, for example, electrostatic potentials, are computed from the enzyme structures. Differences in molecular interaction fields between enzymes are then related to the ratios of their kinetic parameters. This procedure can be used to estimate unknown kinetic parameters when enzyme structural information is available and kinetic parameters have been measured for related enzymes or were obtained under different conditions. The detailed interaction of the enzyme with substrate or cofactors is not modeled and is assumed to be similar for all the proteins compared. The protein structure modeling protocol employed ensures that differences between models reflect genuine differences between the protein sequences, rather than random fluctuations in protein structure.

Conclusion

Provided that the experimental conditions and the protein structural models refer to the same protein state or conformation, correlations between interaction fields and kinetic parameters can be established for sets of related enzymes. Outliers may arise due to variation in the importance of different contributions to the kinetic parameters, such as protein stability and conformational changes. The qPIPSA approach can assist in the validation as well as estimation of kinetic parameters, and provide insights into enzyme mechanism.

Evidence for the importance of electrostatics in the function of two distinct families of ribosome inactivating toxins

Alpha-sarcin and ricin represent two structurally and mechanistically distinct families of site-specific enzymes that block translation by irreversibly modifying the sarcin/ricin loop (SRL) of 23S-28S rRNA. alpha-Sarcin family enzymes are designated as ribotoxins and act as endonucleases. Ricin family enzymes are designated as ribosome inactivating proteins (RIP) and act as N-glycosidases. Recently, we demonstrated that basic surface residues of the ribotoxin restrictocin promote rapid and specific ribosome targeting by this endonuclease. Here, we report that three RIP: ricin A, saporin, and gypsophilin depurinate the ribosome with strong salt sensitivity and achieve unusually fast kcat/Km approximately 10(9)-10(10) M(-1) s(-1), implying that RIP share with ribotoxins a common mechanism of electrostatically facilitated ribosome targeting. Bioinformatics analysis of RIP revealed that surface charge properties correlate with the presence of the transport chain in the RIP molecule, suggesting a second role for the surface charge in RIP transport. These findings put forward surface electrostatics as an important determinant of RIP activity.

On the dynamic nature of the transition state for protein-protein association as determined by double-mutant cycle analysis and simulation

The process of protein-protein association starts with their random collision, which may develop into an encounter complex followed by a transition state and final complex formation. Here we aim to experimentally characterize the nature of the transition state of protein-protein association for three different protein-protein interactions; Barnase-Barstar, TEM1-BLIP and IFNalpha2-IFNAR2, and use the data to model the transition state structures. To model the transition state, we determined inter-protein distance-constraints of the activated complex by using double mutant cycles (DMC) assuming that interacting residues are spatially close. Significant DeltaDeltaG(double dagger)(int) values were obtained only between residues on Barnase and Barstar. However, introducing specific mutations that optimize the charge complementarity between BLIP and TEM1 resulted in the introduction of significant DeltaDeltaG(double dagger)(int) values also between residues of these two proteins. While electrostatic interactions make major contributions towards stabilizing the transition state, we show two examples where steric hindrance exerts an effect on the transition state as well. To model the transition-state structures from the experimental DeltaDeltaG(double dagger)(int) values, we introduced a method for structure perturbation, searching for those inter-protein orientations that best support the experimental DeltaDeltaG(double dagger)(int) values. Two types of transition states were found, specific as observed for Barnase-Barstar and the electrostatically optimized TEM1-BLIP mutants, and diffusive as shown for wild-type TEM1-BLIP and IFNalpha2-IFNAR2. The specific transition states are characterized by defined inter-protein orientations, which cannot be modeled for the diffusive transition states. Mutations introduced through rational design can change the transition state from diffusive to specific. Together, these data provide a structural view of the mechanism allowing rates of association to differ by five orders of magnitude between different protein complexes.