Molecular mechanisms of protein-protein recognition: whether the surface placed charged residues determine the recognition process?

We studied the structure and composition of contact areas in 812 different kind dimeric protein-protein complexes from Brookhaven data base (PDB ) in order to reveal their pecularities with regard to protein-protein recognition. We have found, that the large portion of complexes (approximately 70%) have oppositely charged residues in the contact areas (interfaces) on the subunits surfaces, which form electrostatic contacts - R:E, R:D, K:E, K:D, H:E, H:D. These results are consistent with the current view that high rate complex formation may be driven by the long-range electrostatic interaction between charged AA residues of subunits surfaces. However, there are many complexes among the studied ones (approximately 30%), which have no electrostatic contacts at all in their contact area. Thus a question arises: what forces account for high complex formation rates (i.e. for the distant orienting of subunits before encounter) by forming complexes where the surface contact areas lack electrostatic contacts? We believe that the long-range orienting electrostatic interaction of subunits may account for all cases of efficient complex formation if one drops the traditional view that protein subunits interact mainly through their surfaces. We suggest that the distant orienting being due to the electrostatic interaction between the whole aggregates of partial electric charges of atoms of each complex subunits. Our preliminary model calculations (unpublished) made for ribonuclease dimer (does not have electrostatic contacts) conform this suggestion.     

What forces can determine the formation of highly specific protein-protein complexes?

A software package was designed and used in a detailed study of the contact regions (interfaces) of a large number of protein-protein complexes using the PDB data. It appeared that for about 75% of the complexes the amino acid composition of the subunit surface in the contact region is not essential. Thus one may suggest that, along with the amino acid residues at the interface, the residues in the interior of the globules substantially contribute to protein-protein recognition. Such interactions between quite remote residues are most probably of electrical nature, and are involved in recognition by contributing to the overall electric field created by the protein molecule; the configuration of this field is perhaps the definitive factor of recognition. The overall field of the protein molecule is additively built of the fields created by each constituent residue, and it can be calculated as a sum of the fields created by the protein multipole (aggregate of 'partial' electric charges assigned to every atom of the protein molecule). Preliminary calculations of the remote electrostatic interaction have been performed for ribonuclease subunits in vacuum. The results are indicative of a real possibility that the electric field created by the protein multipole can strongly influence the mutual orientation of molecules before Brownian collisions.     

About factors providing the fast protein-protein recognition in processes of complex formation

A package of programs for the examination of areas of subunit contacts (interface) in protein-protein (PP) complexes has been created and used for a detailed study of amino acid (AA) composition and interface structure in a large number of PP complexes from Brookhaven database (PBD). It appeared that in about 75% of the complexes, the AA composition of the subunit surface is not important. This suggests that, along with the surface AA composition, interactions between AA from the inner parts of protein globules may play a significant role in PP recognition. Such interactions between relatively distant AA residues can only be of electrostatic nature and contribute to the total electric field of the protein molecule. The configuration of the electric field itself appears to determine the PP recognition. The total electric field created by protein molecules can be calculated as a result of superimposition of the fields created by the protein multipole (i.e. by the totality of partial electric charges assigned to each atom of the molecule). We performed preliminary calculations for the distant electrostatic interaction of ribonuclease subunits in a vacuum. The results reveal that the effect of the electric fields of the protein multipole is strong enough to orient protein molecules prior to their Brown collision.     

How optimal are the binding energetics of barnase and barstar?

The extracellular ribonuclease barnase and its intracellular inhibitor barstar bind fast and with high affinity. Although extensive experimental and theoretical studies have been carried out on this system, it is unclear what the relative importance of different contributions to the high affinity is and whether binding can be improved through point mutations. In this work, we first applied Poisson-Boltzmann electrostatic calculations to 65 barnase-barstar complexes with mutations in both barnase and barstar. The continuum electrostatic calculations with a van der Waals surface dielectric boundary definition result in the electrostatic interaction free energy providing the dominant contribution favoring barnase-barstar binding. The results show that the computed electrostatic binding free energy can be improved through mutations at W44/barstar and E73/barnase. Furthermore, the determinants of binding affinity were quantified by applying COMparative BINding Energy (COMBINE) analysis to derive quantitative structure-activity relationships (QSARs) for the 65 complexes. The COMBINE QSAR model highlights approximately 20 interfacial residue pairs as responsible for most of the differences in binding affinity between the mutant complexes, mainly due to electrostatic interactions. Based on the COMBINE model, together with Brownian dynamics simulations to compute diffusional association rate constants, several mutants were designed to have higher binding affinities than the wild-type proteins.     

Helix-helix interactions in lipid bilayers.

Using a continuum model, we calculated the electrostatic interaction free energy between two alpha-helices in three environments: the aqueous phase, a low dielectric alkane phase, and a simple representation of a lipid bilayer. As was found in previous work, helix-helix interactions in the aqueous phase are quite weak, because of solvent screening, and slightly repulsive, because of desolvation effects that accompany helix assembly. In contrast, the interactions can be quite strong in a hypothetical alkane phase because desolvation effects are essentially nonexistent and because helix-helix interactions are not well screened. In this type of environment, the antiparallel helix orientation is strongly favored over the parallel orientation. In previous work we found that the free energy penalty associated with burying helix termini in a bilayer is quite high, which is why the termini tend to protrude into the solvent. Under these conditions the electrostatic interaction is strongly screened by solvent; indeed, it is sufficient for the termini to protrude a few angstroms from the two surfaces of the bilayer for their interaction to diminish almost completely. The effect is consistent with the classical model of the helix dipole in which the dipole moment is represented by point charges located at either terminus. Our results suggest, in agreement with previous models, that there is no significant nonspecific driving force for helix aggregation and, hence, that membrane protein folding must be driven by specific interactions such as close packing and salt-bridge and hydrogen bond formation.     

Helix–coil transitions in a simple polypeptide model

A simplified model of a polypeptide chain is described. Each residue is represented by a single interaction center. The energy of the chain and the force acting on each residue are given as a function of the residue coordinates. Terms to approximate the effect of solvent and the stabilization energy of helix formation are included. The model is used to study equilibrium and dynamical aspects of the helix–coil transition. The equilibrium properties examined include helix–coil equilibrium constants and their dependence on chain position. Dynamical properties are examined by a stochastic simulation of the Brownian motion of the chain in its solvent surroundings. Correlations in the motions of the residues are found to have an important influence on the helix–coil transition rates.     

Significant role of electrostatic interactions for stabilization of protein assemblies

Contribution of electrostatic interactions to stability of BPTI orthorhombic, pig-insulin cubic crystals, and horse L ferritin crystals was evaluated with numerical calculation of Poisson-Boltzmann equation based on a dielectric model. The stability of a ferritin molecule (24-mer) composed of 24 subunits was also evaluated. It was found that the surface charge-charge interactions at separation distances (< 5 Å) were insensitive to variations in the ionic strength, and thus stabilized assembled states of the proteins (i.e., crystalline state and oligomeric state). It was also revealed that the charge density and the packing of the protein crystals were largely responsible for the ionic strength dependence of the crystal stability. The stability of the 5PTI crystalline state with a high charge density drastically increased as the concentration of the solvent ions increased. In contrast, that of the insulin crystal with a low charge density and large solvent region was insensitive to changes in the ionic concentration. The electrostatic interaction between ferritin 24-mers was attributed to two salt bridges mediated by Cd ion. For the stability of the ferritin 24-mer, which is evolutionally designed, the electrostatic stabilization between the subunits was attributed to polar bonds such as buried salt bridges or hydrogen bonds, which occasionally yielded more than 5 kcal/mol and were numerous and very strong compared with the bonds between molecules in the 5PTI and 9INS crystals.By analyzing the atomic charge-charge interactions in detail, it was found that charge pairs separated by less than 3 Å, such as hydrogen bonds, dominantly stabilize the assembled states, and that pairs 3 to 5 Å apart were also important. The stability of the assembled states evaluated by the total EET was determined by the fine balance between the two competing contributions arising from the stabilizing atoms and the destabilizing atoms.Changes of the ASA and hydration free energy were also evaluated in accordance with the process of the subunit assembly. The change of hydration free energy, which was very large (i.e., ~+ 100 kcal/mol/subunit) and unfavorable for the assembly, was proportional to the electrostatic hydration energy (i.e., Born energy change in the hydration process). Hydrophobic groups were likely to appear more frequently than hydrophilic groups at the interfaces.This study offers a method which can improve the stability of protein crystals by introducing polar or charged residues that are properly designed to form specific hydrogen bonds or salt bridges between neighboring protein molecules. This method is also applicable to crystallography, because it improves refinement of protein structures in crystals by taking the inter-protein interactions into account.     

Release factors eRF1 and RF2: a universal mechanism controls the large conformational changes

Class I release factors 1 and 2 (RF1 and RF2) terminate protein synthesis by recognizing stop codons on the mRNA via their conserved amino acid motifs (NIKS in eRF1 and SPF in RF2) and by the conserved tripeptide (GGQ) interactions with the ribosomal peptidyltransferase center. Crystal structures of eRF1 and RF2 do not fit their ribosomal binding pocket (approximately 73 angstroms). Cryoelectron microscopy indicates large conformational changes in the ribosome-bound RF2. Here, we investigate the conformational dynamics of the eRF1 and RF2 using molecular dynamics simulation, structural alignment, and electrostatic analysis of domain interactions. We show that relaxed eRF1 has a shape remarkably similar to the ribosome-bound RF2 observed by cryoelectron microscopy. The similarity between the two release factors is as good as between elongation factor G and elongation factor Tu-guanosine-5'(beta,gamma-imido)triphosphate-tRNA. Further, the conformational transitions and dynamics of eRF1 and RF2 between the free and ribosome-bound states are most likely controlled by protonation of conserved histidines. For eRF1, the distance between the NIKS and GGQ motifs shrinks from 97.5 angstroms in the crystal to 70-80 angstroms. For RF2, the separation between SPF and GGQ elongates from 32 angstroms in the crystal to 50 angstroms. Coulombic interaction strongly favors the open conformation of eRF1; however, solvation and histidine protonation modulate the domain interactions, making the closed conformation of eRF1 more accessible. Thus, RF1 and RF2 function like molecular machines, most likely fueled by histidine protonation. The unified conformational control and the shapes of eRF1 and RF2 support the proposition that the termination of protein synthesis involves similar mechanisms across species.     

Charge-charge interactions influence the denatured state ensemble and contribute to protein stability

Several recent studies have shown that it is possible to increase protein stability by improving electrostatic interactions among charged groups on the surface of the folded protein. However, the stability increases are considerably smaller than predicted by a simple Coulomb's law calculation, and in some cases, a charge reversal on the surface leads to a decrease in stability when an increase was predicted. These results suggest that favorable charge-charge interactions are important in determining the denatured state ensemble, and that the free energy of the denatured state may be decreased more than that of the native state by reversing the charge of a side chain. We suggest that when the hydrophobic and hydrogen bonding interactions that stabilize the folded state are disrupted, the unfolded polypeptide chain rearranges to compact conformations with favorable long-range electrostatic interactions. These charge-charge interactions in the denatured state will reduce the net contribution of electrostatic interactions to protein stability and will help determine the denatured state ensemble. To support this idea, we show that the denatured state ensemble of ribonuclease Sa is considerably more compact at pH 7 where favorable charge-charge interactions are possible than at pH 3, where unfavorable electrostatic repulsion among the positive charges causes an expansion of the denatured state ensemble. Further support is provided by studies of the ionic strength dependence of the stability of charge-reversal mutants of ribonuclease Sa. These results may have important implications for the mechanism of protein folding.     

Estimating the contribution of engineered surface electrostatic interactions to protein stability by using double-mutant cycles

Coulombic interactions between charges on the surface of proteins contribute to stability. It is difficult, however, to estimate their importance by protein engineering methods because mutation of one residue in an ion pair alters the energetics of many interactions in addition to the coulombic energy between the two components. We have estimated the interaction energy between two charged residues, Asp-12 and Arg-16, in an alpha-helix on the surface of a barnase mutant by invoking a double-mutant cycle involving wild-type enzyme (Asp-12, Thr-16), the single mutants Thr----Arg-16 and Asp----Ala-12, and the double mutant Asp----Ala-12, Thr----Arg-16. The changes in free energy of unfolding of the single mutants are not additive because of the coulombic interaction energy. Additivity is restored at high concentrations of salt that shield electrostatic interactions. The geometry of the ion pair in the mutant was assumed to be the same as that in the highly homologous ribonuclease from Bacillus intermedius, binase, which has Asp-12 and Arg-16 in the native enzyme. The ion pair does not form a hydrogen-bonded salt bridge, but the charges are separated by 5-6 A. The mutant barnase containing the ion pair Asp-12/Arg-16 is more stable than wild type by 0.5 kcal/mol, but only a part of the increased stability is attributable to the electrostatic interaction. We present a formal analysis of how double-mutant cycles can be used to measure the energetics of pairwise interactions.     

What Forces Can Determine the Formation of Highly Specific Protein–Protein Complexes?

A software package was designed and used in a detailed study of the contact regions (interfaces) of a large number of protein–protein complexes using the PDB data. It appeared that for about 75% of the complexes the amino acid composition of the subunit surface in the contact region is not essential. Thus one may suggest that, along with the amino acid residues at the interface, the residues in the interior of the globules substantially contribute to protein–protein recognition. Such interactions between quite remote residues are most probably of electrical nature, and are involved in recognition by contributing to the overall electric field created by the protein molecule; the configuration of this field is perhaps the definitive factor of recognition. The overall field of the protein molecule is additively built of the fields created by each constituent residue, and it can be calculated as a sum of the fields created by the protein multipole (aggregate of partial electric charges assigned to every atom of the protein molecule). Preliminary assessment of the remote electrostatic interaction has been performed for ribonuclease subunits in vacuum. The results are indicative of a real possibility that the electric field created by the protein multipole can strongly influence the mutual orientation of molecules before Brownian collision.     

Calculation of the energy difference between the quaternary structures of deoxy- and oxyhemoglobin.

The energy difference between the quaternary structures of deoxy- and oxyhemoglobin is evaluated on the basis of the atomic coordinates determined by X-ray diffraction analysis. Calculation of the van der Waals interaction between subunits shows that in a hemoglobin molecule as a whole, the interaction is more attractive in the oxy form than in the deoxy form by about 8 kcal/mol, and that in each pair of two subunits except the pair alpha1alpha2, the interaction energy varies by about 15 kcal/mol. The electrostatic interactions originating in the partial charges on all constituent atoms of hemoglobin and in the polar residues on the surface of hemoglobin make only a small contribution to the energy difference between the quaternary structures of deoxy- and oxyhemoglobin. Thus, the contribution of the clusters of the polar residues in the internal cavity between like subunits and also of the freedom of rotation of the C-terminal of each subunit in oxyhemoglobin may be important energetically in the transition from deoxy to oxy quaternary structure. In this point, the present calculation supports Perutz' model, but suggests necessity of further investigations on the transitional characteristics of the quaternary structure in the intermediate steps of oxygenation. The discussion on the transitional characteristics is given in the last section.     

Structural basis of heteromeric smad protein assembly in TGF-beta signaling

The formation of protein complexes between phosphorylated R-Smads and Smad4 is a central event in the TGF-beta signaling pathway. We have determined the crystal structure of two R-Smad/Smad4 complexes, Smad3/Smad4 to 2.5 angstroms, and Smad2/Smad4 to 2.7 angstroms. Both complexes are heterotrimers, comprising two phosphorylated R-Smad subunits and one Smad4 subunit, a finding that was corroborated by isothermal titration calorimetry and mutational studies. Preferential formation of the R-Smad/Smad4 heterotrimer over the R-Smad homotrimer is largely enthalpy driven, contributed by the unique presence of strong electrostatic interactions within the heterotrimeric interfaces. The study supports a common mechanism of Smad protein assembly in TGF-beta superfamily signaling.     

Pseudosymmetry near a magnetic surface in a plasma confinement system

The possibility is demonstrated of finding vacuum equilibrium magnetic configurations with an exactly pseudosymmetric nonparaxial boundary magnetic surface in the vicinity of which the pseudosymmetry condition is satisfied approximately. Equations are derived for calculating the boundary surface from a prescribed particular dependence of the magnetic field strength in special magnetic flux coordinates. In calculations, magnetic coordinates serve as ordinary angular coordinates, while their “magnetic” character is specified by additional integral conditions. As an example, a “tubular” orthogonal magnetic surface is calculated analytically.     

On the energetics of conformational changes and pH dependent redox behaviour of electron transfer proteins

Calculations of the electrostatic interaction energies for four metalloproteins that carry out electron transfer are reported. Each protein has a pH dependent redox potential from which the measured electrostatic interaction energy is obtained. The calculations were made using the X-ray structure coordinates and a semimacroscopic model of the interactions. For cytochrome c-551 and HIPIP the calculated and observed interaction energies were found to be approximately the same, in agreement with the fact that significant conformational changes do not accompany the ionisations. For cytochrome c2 and azurin, however, major differences were found between the calculated and observed values. These are accounted for primarily by the occurrence of significant conformational changes accompanying the ionisations. The reorganisation energies for these conformational changes are approximately 7.0 and approximately 11.1 kJ.mol-1, respectively.     

Implementation of an experimental and computational tool set to study protein-mAb interactions

This work focused on the development of a combined experimental and computational tool set to study protein-mAb interactions. A model protein library was first screened using cross interaction chromatography to identify proteins with the strongest retention. Fluorescence polarization was then employed to study the interactions and thermodynamics of the selected proteins—lactoferrin, pyruvate kinase, and ribonuclease B with the mAb. Binding affinities of lactoferrin and pyruvate kinase to the mAb were seen to be relatively salt insensitive in the range examined. Further, a strong entropic contribution was observed, suggesting the importance of hydrophobic interactions. On the other hand, ribonuclease B-mAb binding was seen to be enthalpically driven and salt sensitive, indicating the importance of electrostatic interactions. Protein–protein docking was then carried out and the results identified the CDR region on the mAb as an important binding site for all three proteins. The binding interfaces identified for the mAb-lactoferrin and mAb-pyruvate kinase systems were found to contain complementary hydrophobic and oppositely charged clusters on the interacting regions which were indicative of both hydrophobic and electrostatic interactions. On the other hand, the binding site on ribonuclease B was predominantly positively charged with minimal hydrophobicity. This resulted in an alignment with negatively charged clusters on the mAb, supporting the contention that these interactions were primarily electrostatic in nature. Importantly, these computational results were found to be consistent with the fluorescence polarization data and this combined approach may have utility in examining mAb-HCP interactions which can often complicate the downstream processing of biologics. © 2019 American Institute of Chemical Engineers     

Brownian dynamics simulation of the lateral distribution of charged membrane components

Brownian dynamics simulations were performed to study the contribution of electric interactions between charged membrane components to their lateral distribution in a two-dimensional viscous liquid (bilayer lipid membrane). The electrostatic interaction potential was derived from an analytical solution of the linearized Poisson-Boltzmann equation for point charges in an electrolyte solution — membrane — electrolyte solution system. Equilibrium as well as dynamic quantities were investigated. The lateral organization of membrane particles, modelled by mobile cylinders in a homogeneous membrane separating two electrolyte solutions was described by spatial distribution functions, diffusion coefficients and cluster statistics. Disorder, local order and crystal-like arrangements were observed as a function of the particle charge, the closest possible distances between the charges and the particle density. The simulations revealed that the system is very sensitive to the position of the charges with respect to the electrolyte solution — membrane interface. Electrostatic interactions of charges placed directly on the membrane surface were almost negligible, whereas deeper charges demonstrated pronounced interaction. Biologically relevant parameters corresponded at most to local and transient ordering. It was found that lateral electric forces can give rise to a preferred formation of clusters with an even number of constituents provided that the closest possible charge-charge distances are small. It is concluded that lateral electrostatic interactions can account for local particle aggregations, but their impact on the global arrangement and movement of membrane components is limited. Correspondence to: D. Walther     

Long range RNA-RNA interactions in the 30 S ribosomal subunit of E. coli.

We have attempted to identify long-range interactions in the tertiary structure of RNA in the E. coli 30 S ribosome. Native subunits were cleaved with ribonuclease and separated into nucleoprotein fragments which were deproteinized and fractionated into multi-oligonucleotide complexes under conditions intended to preserve RNA-RNA interactions. The final products were denatured by urea and heat and their constituent oligonucleotides resolved and sequenced. Many complexes contained complementary sequences known to be bound together in the RNA secondary structure, attesting to the validity of the technique. Other co-migrating oligonucleotides, not joined in the secondary structure, contained mutually complementary sequences in locations that allow base-pairing interaction without disrupting pre-existing secondary structure. In seven instances the complementary relationship was found to have been preserved during phylogenetic diversification.     

Incorporation of lone pairs into the electrostatic term in the empirical intra- and intermolecular energy calculations

The method hitherto used for estimating the electrostatic term in empirical intramolecular calculations of stable conformations of biologically important molecules and macromolecules and intermolecular calculations of molecular associations or packing energy in molecular crystals had been analyzed. It has been shown that the contribution of atomic hybridization moments is omitted in the calculation of electrostatic interactions from net atomic charges localized on nuclei which have been determined by standard quantum-chemical methods. This contribution plays an important part in determining electrostatic interactions, mainly in molecules containing atoms with lone pairs. Simultaneously, a modified method for calculating the electrostatic term comprising the interaction of the lone pairs, which are represented by atomic hybridization moments, has been proposed. The relationship between the atomic hybridization moment and the bond angle has been expressed for some typical configurations occurring in biologically important molecules. Finally, this new approach is illustrated by results of the conformational analysis of some model compounds for biomolecules and compared with the approach used so far for the estimation of the electrostatic interaction in empirical methods of calculation of the intra- and intermolecular energy.