Calculation of weak protein-protein interactions: the pH dependence of the second virial coefficient

Atomically detailed simulations of concentrated protein solutions: the effects of salt, pH, point mutations, and protein concentration in simulations of 1000-molecule systems

An ability to accurately simulate the dynamic behavior of concentrated macromolecular solutions would be of considerable utility in studies of a wide range of biological systems. With this goal in mind, a Brownian dynamics (BD) simulation method is reported here that allows systems to be modeled that comprise in excess of 1000 protein molecules, all of which are treated in atomic detail. Intermolecular forces are described in the method using an energy function that incorporates electrostatic and hydrophobic interactions and that is calibrated to reproduce experimental thermodynamic information with a single adjustable parameter. Using the method, BD simulations have been performed over a wide range of pH and ionic strengths for three proteins: hen egg white lysozyme (HEWL), chymotrypsinogen, and T4 lysozyme. The simulations reproduce experimental trends in second virial coefficients (B(22)) and translational diffusion coefficients, correctly capture changes in B(22) values due to single amino acid substitutions, and reveal a new explanation for the difficulties reported previously in the literature in reproducing B(22) values for protein solutions of very low ionic strength. In addition, a strong correlation is found between a residue's probability of being involved in a protein-protein contact in the simulations and its probability of being involved in an experimental crystal contact. Finally, exploratory simulations of HEWL indicate that the simulation model also gives a promising description of behavior at very high protein concentrations (approximately 250 g/L), suggesting that it may provide a suitable computational framework for modeling the complex behavior exhibited by macromolecules in cellular conditions.

Why are proteins charged? Networks of charge-charge interactions in proteins measured by charge ladders and capillary electrophoresis

Almost all proteins contain charged amino acids. While the function in catalysis or binding of individual charges in the active site can often be identified, it is less clear how to assign function to charges beyond this region. Are they necessary for solubility? For reasons other than solubility? Can manipulating these charges change the properties of proteins? A combination of capillary electrophoresis (CE) and protein charge ladders makes it possible to study the roles of charged residues on the surface of proteins outside the active site. This method involves chemical modification of those residues to generate a large number of derivatives of the protein that differ in charge. CE separates those derivatives into groups with the same number of modified charged groups. By studying the influence of charge on the properties of proteins using charge ladders, it is possible to estimate the net charge and hydrodynamic radius and to infer the role of charged residues in ligand binding and protein folding.

Ultrasonic storage modulus as a novel parameter for analyzing protein-protein interactions in high protein concentration solutions: correlation with static and dynamic light scattering measurements

The purpose of this work was to establish ultrasonic storage modulus (G') as a novel parameter for characterizing protein-protein interactions (PPI) in high concentration protein solutions. Using an indigenously developed ultrasonic shear rheometer, G' for 20-120 mg/ml solutions of a monoclonal antibody (IgG(2)), between pH 3.0 and 9.0 at 4 mM ionic strength, was measured at frequency of 10 MHz. Our understanding of ultrasonic rheology indicated decrease in repulsive and increase in attractive PPI with increasing solution pH. To confirm this behavior, dynamic (DLS) and static (SLS) light scattering measurements were conducted in dilute solutions. Due to technical limitations, light scattering measurements could not be conducted in concentrated solutions. Mutual-diffusion coefficient, measured by DLS, increased with IgG(2) concentration at pH 4.0 and this trend reversed as pH was increased to 9.0. Second virial coefficient, measured by SLS, decreased with increasing pH. These observations were consistent with the nature of PPI understood from G' measurements. Ultrasonic rheology, DLS, and SLS measurements were also conducted under conditions of increased ionic strength. The consistency between rheology and light scattering analysis under various solution conditions established the utility of ultrasonic G' measurements as a novel tool for analyzing PPI in high protein concentration systems.

Application of high‐frequency rheology measurements for analyzing protein–protein interactions in high protein concentration solutions using a model monoclonal antibody (IgG2)

The purpose of this work was to explore the utilization of high-frequency rheology analysis for assessing protein-protein interactions in high protein concentration solutions. Rheology analysis of a model monoclonal immunoglobulin G2 solutions was conducted on indigenously developed ultrasonic shear rheometer at frequency of 10 MHz. Solutions at pH 9.0 behaved as most viscous and viscoelastic whereas those at pH 4.0 and 5.4 exhibited lower viscosity and viscoelasticity, respectively. Intrinsic viscosity, hydrophobicity, and conformational analysis could not account for the rheological behavior of IgG2 solutions. Zeta potential and light scattering measurements showed the significance of electroviscous and specific protein-protein interactions in governing rheology of IgG2 solutions. Specific protein-protein interactions resulted in formation of reversible higher order species of monomer. Solution storage modulus (G'), and not loss modulus or complex viscosity, was the more reliable parameter for predicting protein-protein interactions. Predictions about the nature of protein-protein interactions made on the basis of solution G' were found to be consistent with observed effect of pH and ionic strength on zeta potential and scattered intensity of IgG2 solutions. Results demonstrated the potential of high-frequency storage modulus measurements for understanding behavior of proteins in solutions and predicting the nature of protein-protein interactions.

Interaction energy decomposition in protein-protein association: a quantum mechanical study of barnase-barstar complex

Protein-protein interactions are very important in the function of a cell. Computational studies of these interactions have been of interest, but often they have utilized classical modelling techniques. In recent years, quantum mechanical (QM) treatment of entire proteins has emerged as a powerful approach to study biomolecular systems. Herein, we apply a semi-empirical divide and conquer (DC) methodology coupled with a dielectric continuum model for the solvent, to explore the contribution of electrostatics, polarization and charge transfer to the interaction energy between barnase and barstar in their complex form. Molecular dynamic (MD) simulation was performed to account for the dynamic behavior of the complex. The results show that electrostatics, charge transfer and polarization favor the formation of the complex. Our study shows that electrostatics dominates the interaction between barnase and barstar ( approximately 73%), while charge transfer and polarization are approximately 21% and approximately 6%, respectively. Close inspection of the polarization and charge-transfer effects on the charge distribution of the complex reveals the existence of two, well localized, regions in barstar. The first region includes the residues between P27 and Y47 and the second region is between N65 and D83. Since no such regions could be detected in barnase clearly suggests that barstar is well optimized for efficiently binding barnase. Furthermore, using our interaction energy decomposition scheme, we were able to identify all residues that have been experimentally determined to be important for the complex formation and to suggest other residues never have been investigated. This suggests that our approach will be useful as an aid in further understanding protein-protein contacts for the ultimate goal to produce successful inhibitors for protein complexes.

Investigation of interaction between enolase and phosphoglycerate mutase using molecular dynamics simulation

Two glycolytic enzymes, phosphoglycerate mutase (PGM) and enolase from Saccharomyces cerevisiae have been chosen to detect complex formation between active centers (a/c), using molecular dynamics simulation. Enzymes have been separated by 10 A distance and placed in a water box of size 173 x 173 x 173 A. Three different orientations where a/c of PGM and enolase were positioned toward each other have been used for investigation. The two initial 3-phosphoglycerate substrates at near active centers of initial structure of PGM have been replaced with final 2-phosphoglycerate products. 150mM of NaCl have been added to the system to observe binding activity in the near physiological conditions. Analysis of interaction energies and conformation changes for 3ns simulation indicates that PGM and enolase do show binding affinity between their near active regions. Moreover the similarity between final conformations of the first two orientations with the initial conformation of the third orientation suggests that complex formation between a/c of enzymes is not confined only by discussed orientations. Clear interaction of enolase with C-terminal tail of PGM has been recorded. These results suggest that substrate direct transfer mechanism may exist between enzymes.

Oligomerization of nitrophorins

Rhodnius prolixus is a blood-sucking insect that uses a mixture of nitrophorin (NP) proteins to deliver nitric oxide (NO) from the insect saliva to the hosts via a ferric heme coordinated to the protein, causing vasodilatation and anticoagulation to support their feeding. R. prolixus NPs 1-4 are very similar proteins ( approximately 20 kDa) with different NO affinities for stepwise NO release triggered by pH increase and histamine binding in hosts. Ultra-high-resolution X-ray structures of native and mutant NPs and their kinetic analysis already have revealed the fundamental steps of NO binding and release. In this study, we found that NPs can exist in multiple oligomerization states at higher concentrations. The oligomers are characterized by a combination of multiple biophysical methods. The intrinsic features of the oligomerization revealed here led us to propose that this intensive, moderately pH- and ligand-dependent oligomerization of NPs has physiological implications in the facilitation of the efficient storage and release of the highly reactive NO in the insect saliva and the victim, respectively.

The Combined Simulation Approach of Atomistic and Continuum Models for the Thermodynamics of Lysozyme Crystals

We have studied the thermodynamic properties of hen egg white lysozyme crystals using a novel simulation method combining atomistic Monte Carlo simulation to calculate van der Waals interactions and the boundary element method to solve the Poisson-Boltzmann equation for the electrostatic interactions. For computational simplicity, we treat the protein as a rigid body, using the crystallographic coordinates of all non-hydrogen atoms of the protein to describe the detailed shape. NVT Monte Carlo simulations are carried out for tetragonal and orthorhombic crystals to obtain the van der Waals energy, incorporating an implicit solvation effect. For crystal phases, an optimally linearized Poisson-Boltzmann equation is used to include the effect of the Donnan equilibrium of the salt ions. The Helmholtz energy is obtained from expanded ensemble Monte Carlo simulations. By using the force field parameters that had previously been tuned for the solution properties, reasonable agreement with experiment is found for the crystallization energy of the tetragonal form. The prediction of the entropy is also reasonable with a slight underestimation suggesting the release of a few water molecules per protein during the crystallization. However, the predictions of the properties of the orthorhombic crystal are poor, probably due to differences in the solvation structure as indicated by experiments, and also as a result of the approximate force field used.

Lysozyme dimerization: Brownian dynamics simulation

The lysozyme dimerization reaction has been studied within the framework of encounter-complex (EC) formation theory using the MacroDox software package. Two types of energetically favorite ECs were determined. In the first of them, active-center amino acids of lysozyme take part in the complex formation or the second molecule blocks accessibility to active center sterically. Epitope amino-acid residues are involved in the complex of type II. The existence of both types of complexes does not contradict experimental data. Dimer-formation rate constants for different kinds of EC were calculated. Increasing the pH from 2.0 to 10.0 decreases the total positive lysozyme charge and eliminates the unfavorable repulsive electrostatic interaction. The rate constant of EC formation is inversely proportional to the protein total charge. The association rate constant was also enhanced by an increase of ionic strength that screened repulsive electrostatic interaction between positively charged proteins.

Potentials of mean force for the interaction of blocked alanine dipeptide molecules in water and gas phase from MD simulations

We calculate potentials of mean force (PMFs) for the intermolecular interaction of two blocked alanine dipeptide (AcAlaNHMe) molecules in water and gas phase at two temperatures, 278 and 300 K, from all-atom molecular dynamics simulations. Simple models based on buried solvent accessible surface and one-dimensional potentials derived from distance-based radial distribution functions are not capable of expressing the short- and long-range complexity of the solute-solute interactions in water. Instead, radial and angular variations in the PMFs are observed with the two-dimensional potentials. The strength of the interactions for specific relative orientations of the molecules in the two-dimensional PMFs is more than double that observed in the one-dimensional PMFs. The populations of specific blocked alanine dipeptide conformations in water, such as alpha(R) and PPII, vary with temperature, and most significantly, with the distance between the centers of mass. A preference for helical conformations is observed at close encounter between molecules.

Reversal of some viral IL-6 electrostatic properties compared to IL-6 contributes to a loss of alpha receptor component recruitment

Human interleukin-6 (hIL-6) is a pleiotropic mediator of activation and proliferation across a large number of different cell types. Human herpesvirus-8 (HHV-8) has been associated with classical and AIDS-related Kaposi's sarcoma (KS). HHV-8 encodes viral IL-6 (vIL-6), a functional homolog of human interleukin-6, that promotes the growth of KS and of some lymphoma cells. Signaling induced by human IL-6 requires recruitment of the glycoprotein gp130, which acts as the signal transducing chain, and of IL-6Ralpha, which is necessary for cognate recognition and high affinity receptor complex formation. In contrast, the formation of a functional complex between vIL-6 and gp130 does not require the presence of IL-6Ralpha. The physico-chemical properties of vIL-6 have been analyzed and compared to those of hIL-6 and of the receptor chains, gp130 and IL-6Ralpha. Interaction sites on vIL-6 involve more hydrophobic residues than those of hIL-6. The electrostatic fields induced by vIL-6 and IL-6Ralpha are repulsive and prevent interaction between vIL-6 and IL-6Ralpha, whereas the electrostatic field induced by hIL-6 steers the complex formation with IL-6Ralpha. Subsequently, electrostatic binding free energy in the vIL-6/IL-6Ralpha complex is destabilizing, whereas it is stabilizing in the complex comprising hIL-6. These properties result from charge reversals between viral and human IL-6, an unusual phenomenon of amino acid substitutions within a homologous protein family. This suggests a selection pressure for vIL-6 to by-pass the IL-6Ralpha control of host defense against virus infection. This selection pressure has yielded the reversal of electrostatic properties of vIL-6 when compared to hIL-6.

Light-scattering studies of protein solutions: role of hydration in weak protein-protein interactions

We model the hydration contribution to short-range electrostatic/dispersion protein interactions embodied in the osmotic second virial coefficient, B(2), by adopting a quasi-chemical description in which water molecules associated with the protein are identified through explicit molecular dynamics simulations. These water molecules reduce the surface complementarity of highly favorable short-range interactions, and therefore can play an important role in mediating protein-protein interactions. Here we examine this quasi-chemical view of hydration by predicting the interaction part of B(2) and comparing our results with those derived from light-scattering measurements of B(2) for staphylococcal nuclease, lysozyme, and chymotrypsinogen at 25 degrees C as a function of solution pH and ionic strength. We find that short-range protein interactions are influenced by water molecules strongly associated with a relatively small fraction of the protein surface. However, the effect of these strongly associated water molecules on the surface complementarity of short-range protein interactions is significant, and must be taken into account for an accurate description of B(2). We also observe remarkably similar hydration behavior for these proteins despite substantial differences in their three-dimensional structures and spatial charge distributions, suggesting a general characterization of protein hydration.

Electrostatic interactions between double layers: influence of surface roughness, regulation, and chemical heterogeneities

Electrostatic interactions between two surfaces as measured by atomic force microscopy (AFM) are usually analyzed in terms of DLVO theory. The discrepancies often observed between the experimental and theoretical behavior are usually ascribed to the occurrence of chemical regulation processes and/or to the presence of surface chemical or morphological heterogeneities (roughness). In this paper, a two-gradient mean-field lattice analysis is elaborated to quantifying double layer interactions between nonplanar surfaces. It allows for the implementation of the aforementioned sources of deviation from DLVO predictions. Two types of ion-surface interaction ensure the adjustment of charges and potentials upon double layer overlap, i.e., specific ionic adsorption at the surfaces and/or the presence of charge-determining ions for the surfaces considered. Upon double layer overlap, charges and potentials are adjusted via reequilibrium of the different ion adsorption processes. Roughness is modeled by grafting asperities on supporting planar surfaces, with their respective positions, shapes, and chemical properties being assigned at will. Local potential and charge distributions are derived by numerically solving the nonlinear Poisson-Boltzmann equation under the boundary conditions imposed by the surface profiles and regulation mechanism chosen. Finite size of the ions is taken into account. A number of characteristic situations are briefly discussed. It is shown how the surface irregularities are reflected in the Gibbs energy of interaction.

Prediction of RNA-binding proteins from primary sequence by a support vector machine approach

Elucidation of the interaction of proteins with different molecules is of significance in the understanding of cellular processes. Computational methods have been developed for the prediction of protein-protein interactions. But insufficient attention has been paid to the prediction of protein-RNA interactions, which play central roles in regulating gene expression and certain RNA-mediated enzymatic processes. This work explored the use of a machine learning method, support vector machines (SVM), for the prediction of RNA-binding proteins directly from their primary sequence. Based on the knowledge of known RNA-binding and non-RNA-binding proteins, an SVM system was trained to recognize RNA-binding proteins. A total of 4011 RNA-binding and 9781 non-RNA-binding proteins was used to train and test the SVM classification system, and an independent set of 447 RNA-binding and 4881 non-RNA-binding proteins was used to evaluate the classification accuracy. Testing results using this independent evaluation set show a prediction accuracy of 94.1%, 79.3%, and 94.1% for rRNA-, mRNA-, and tRNA-binding proteins, and 98.7%, 96.5%, and 99.9% for non-rRNA-, non-mRNA-, and non-tRNA-binding proteins, respectively. The SVM classification system was further tested on a small class of snRNA-binding proteins with only 60 available sequences. The prediction accuracy is 40.0% and 99.9% for snRNA-binding and non-snRNA-binding proteins, indicating a need for a sufficient number of proteins to train SVM. The SVM classification systems trained in this work were added to our Web-based protein functional classification software SVMProt, at http://jing.cz3.nus.edu.sg/cgi-bin/svmprot.cgi. Our study suggests the potential of SVM as a useful tool for facilitating the prediction of protein-RNA interactions.

Mitochondrial Ca2+ uptake is important over low [Ca2+]i range in arterial smooth muscle

Mitochondrial Ca(2+) uptake is usually thought to occur only when intracellular Ca(2+) concentration ([Ca(2+)](i)) is high. We investigated whether mitochondrial Ca(2+) removal participates in shaping [Ca(2+)](i) signals in arterial smooth muscle over a low [Ca(2+)](i) range. [Ca(2+)](i) was measured using fura 2-loaded, voltage-clamped cells from rat femoral arteries. Both diazoxide and carbonyl cyanide m-chlorophenylhydrazone (CCCP) depolarized the mitochondria. Diazoxide application increased resting [Ca(2+)](i), suggesting that Ca(2+) is sequestered in mitochondria. Over a low [Ca(2+)](i) range, diazoxide and CCCP slowed Ca(2+) removal rate, determined after a brief depolarization. When [Ca(2+)](i) was measured during sustained depolarization to -30 mV, CCCP application increased [Ca(2+)](i). When Ca(2+) transients were repeatedly evoked by caffeine applications, CCCP application elevated resting [Ca(2+)](i). Caffeine-induced Ca(2+) transients were compared before and after CCCP application using the half decay time, or time required to reduce increase in [Ca(2+)](i) by 50% (t((1/2))). CCCP treatment significantly increased t((1/2)). These results suggest that Ca(2+) removal to mitochondria in arterial smooth muscle cells may be important at a low [Ca(2+)](i).

Protein-protein interactions in concentrated electrolyte solutions

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

Hydrophobic forces between protein molecules in aqueous solutions of concentrated electrolyte

Protein-protein interactions have been measured for a mutant (D101F) lysozyme and for native lysozyme in concentrated solutions of ammonium sulfate at pH 7 and sodium chloride at pH 4.5. In the mutant lysozyme, a surface aspartate residue has been replaced with a hydrophobic phenylalanine residue. The protein-protein interactions of D101F lysozyme are more attractive than those of native lysozyme for all conditions studied. The salt-induced attraction is correlated with a solvation potential of mean force given by the work required to desolvate the part of the protein surfaces that is buried by the protein-protein interaction. This work is proportional to the aqueous surface-tension increment of the salt and the fractional non-polar surface coverage of the protein. Experimental measurements of osmotic second virial coefficients validate a proposed potential of mean force that ascribes the salt-induced attraction between protein molecules to an enhancement of the hydrophobic attraction. This model provides a first approximation for predicting the protein-protein potential of mean force in concentrated aqueous electrolyte solutions; this potential is useful for determining solution conditions favorable for protein crystallization.

Prediction of protein-protein interactions by docking methods

Recently, developments have been made in predicting the structure of docked complexes when the coordinates of the components are known. The process generally consists of a stage during which the components are combined rigidly and then a refinement stage. Several rapid new algorithms have been introduced in the rigid docking problem and promising refinement techniques have been developed, based on modified molecular mechanics force fields and empirical measures of desolvation, combined with minimisations that switch on the short-range interactions gradually. There has also been progress in developing a benchmark set of targets for docking and a blind trial, similar to the trials of protein structure prediction, has taken place.