Electrostatic Interactions In Proteins Research Articles

Unraveling the Impact of pH on the Crystallization of Pharmaceutical Proteins: A Case Study of Human Insulin.

One of the most crucial parameters in protein crystallization is pH, as it governs the protein’s electrostatic interactions. However, the fundamental role of pH on crystallization still remains unknown. Here, we systematically investigated the crystallization of human insulin (isoelectric point 5.3) at various pHs between 6.0 and 6.7 at different supersaturation ratios, up to 20.9. Our results demonstrate that the pH has an opposing effect on solubility and nucleation rate as a shift in pH toward a more basic milieu increases the solubility by 5-fold while the onset of nucleation was accelerated by a maximum of 8.6-fold. To shed light on this opposing effect, we evaluated the protein–protein interactions as a function of pH by measuring the second virial coefficient and hydrodynamic radius and showed that a change in pH of less than one unit has no significant impact on the protein–protein interactions. As it is widely understood that the increase in protein solubility as a function of pH is due to the increase in the repulsive electrostatic interactions, we have demonstrated that the increase in insulin solubility and decrease in the onset of nucleation are independent of the protein–protein interactions. We hypothesize that it is the electrostatic interactions between both ions and solvent molecules and the protein residues that are governing the crystallization of human insulin. The findings of this study will be of crucial importance for the design of novel crystallization pathways.

Molecular dynamic investigate the affection of EGFR by Tubemoside

Tubemoside as a common traditional Chinese medicine is playing an important role in the field of prevention and treatment of lung cancer without any side effects. However, the reason and its mechanism remain unclear. In our study, the molecular dynamic simulation was used to investigate the mechanism at the molecular level. We found that the hydrogen bond network of proteins (three states of EGFR) was affected by Tubemoside. The movement and opening/closing state of protein was changed when combine with Tubemoside. The results of principal component analysis were used to prove the transform of proteins and the change of its movement. Electrostatic interactions of proteins also were studied. The numbers of active interaction sites will decrease while Tubemoside emerged in the protein, which will cause the activity change of EGFR for forming asymmetric dimers required for activation.

Loss of a conserved salt bridge in bacterial glycosyl hydrolase BgIM-G1 improves substrate binding in temperate environments

Salt bridges are the strongest electrostatic interactions in proteins. They substantially contribute to a protein’s structural stability. Thus, mutations of salt bridges are typically selected against. Here, we report on the evolutionary loss of a highly conserved salt bridge in the GH1 family glycosyl hydrolase BglM-G1. BglM-G1’s gene was found in the bacterial metagenome of a temperate, seasonally cold marine habitat. In BglM-G1, arginine 75 is replaced by a histidine. While fully retaining β-glucosidase activity, BglM-G1 is less heat stable than an H75R variant, in which the salt bridge was artificially re-introduced. However, the Km toward its substrates was lower in wild type, leading to an overall higher catalytic efficiency. Our results indicate that this loss of the salt bridge leads to higher flexibility in BglM-G1’s active site, trading structural stability at high temperatures, a trait not needed in a temperate, seasonally cold habitat, for a more effective catalytic activity.

Electrostatic Interactions in Protein Structure, Folding, Binding, and Condensation.

Charged and polar groups, through forming ion pairs, hydrogen bonds, and other less specific electrostatic interactions, impart important properties to proteins. Modulation of the charges on the amino acids, e.g., by pH and by phosphorylation and dephosphorylation, have significant effects such as protein denaturation and switch-like response of signal transduction networks. This review aims to present a unifying theme among the various effects of protein charges and polar groups. Simple models will be used to illustrate basic ideas about electrostatic interactions in proteins, and these ideas in turn will be used to elucidate the roles of electrostatic interactions in protein structure, folding, binding, condensation, and related biological functions. In particular, we will examine how charged side chains are spatially distributed in various types of proteins and how electrostatic interactions affect thermodynamic and kinetic properties of proteins. Our hope is to capture both important historical developments and recent experimental and theoretical advances in quantifying electrostatic contributions of proteins.

Endocrine Imbalance Associated With Proteome Changes in Diabetes.

The dynamics of cellular metabolism involves rapid interactions between proteins and nucleic acids, proteins and proteins, and signaling. These involve the interactions with respect to the sulfur bond, noncovalent electrostatic interactions, protein structure stabilization and protein-ligand binding, weak electrostatic interactions in proteins, oxygen radicals that initiate a change in conformation and a chain of events. We review a development in molecular medicine that is a very promising work in progress. We also review the current and future research methods involving mitochondria. Long-term effects of diabetes include glycation of proteins, for example, glycohemoglobin (HbA1c), increased risk of cardiovascular diseases, atherosclerosis, retinopathy, nephropathy, and neurological dysfunctions. Tissues are exposed to significant quantities of highly reactive chemical species including nitric oxide • NO and reactive oxygen species ROS over months to years, to an extent generated by mitochondrial activities. The reactions of • NO can be broadly discussed with reference to three main processes which control their fate in biological systems: (1) diffusion and intra-cellular consumption; (2) autooxidation to form nitrous anhydride N2 O3 ; and (3) reaction with superoxide O2• - to form peroxynitrite ONOO-. Reactive nitrogen species produced by macrophages and neutrophils in the interstitial space, with emphasis on • NO, N2 O3 , ONOO-, and nitrogen dioxide radicals • NO2 generate protein and DNA damage. Serum thiol (-SH) groups act as an important extracellular scavenger of peroxides and are therefore helpful in protecting the surrounding tissues. The events described here are a homeostatic endocrine imbalance that is associated with proteostasis. The advances we have seen in untangling this web of interactions are sure to continue at a breathtaking pace. J. Cell. Biochem. 118: 3569-3576, 2017. © 2017 Wiley Periodicals, Inc.

Pick and Choose the Spectroscopic Method to Calibrate the Local Electric Field inside Proteins.

Electrostatic interactions in proteins play a crucial role in determining the structure-function relation in biomolecules. In recent years, fluorescent probes have been extensively employed to interrogate the polarity in biological cavities through dielectric constants or semiempirical polarity scales. A choice of multiple spectroscopic methods, not limited by fluorophores, along with a molecular level description of electrostatics involving solute-solvent interactions, would allow more flexibility to pick and choose the experimental technique to determine the local electrostatics within protein interiors. In this work we report that ultraviolet/visible-absorption, infrared-absorption, or (13)C NMR can be used to calibrate the local electric field in both hydrogen bonded and non-hydrogen bonded protein environments. The local electric field at the binding site of a serum protein has been determined using the absorption wavelength as well as the carbonyl stretching frequency of its natural steroid substrate, testosterone. Excellent agreement is observed in the results obtained from two independent spectroscopic techniques.

Statistical analysis of protein structures suggests that buried ionizable residues in proteins are hydrogen bonded or form salt bridges

It is well known that nonpolar residues are largely buried in the interior of proteins, whereas polar and ionizable residues tend to be more localized on the protein surface where they are solvent exposed. Such a distribution of residues between surface and interior is well understood from a thermodynamic point: nonpolar side chains are excluded from the contact with the solvent water, whereas polar and ionizable groups have favorable interactions with the water and thus are preferred at the protein surface. However, there is an increasing amount of information suggesting that polar and ionizable residues do occur in the protein core, including at positions that have no known functional importance. This is inconsistent with the observations that dehydration of polar and in particular ionizable groups is very energetically unfavorable. To resolve this, we performed a detailed analysis of the distribution of fractional burial of polar and ionizable residues using a large set of ˜2600 nonhomologous protein structures. We show that when ionizable residues are fully buried, the vast majority of them form hydrogen bonds and/or salt bridges with other polar/ionizable groups. This observation resolves an apparent contradiction: the energetic penalty of dehydration of polar/ionizable groups is paid off by favorable energy of hydrogen bonding and/or salt bridge formation in the protein interior. Our conclusion agrees well with the previous findings based on the continuum models for electrostatic interactions in proteins.

Effects of surface curvature and surface chemistry on the structure and activity of proteins adsorbed in nanopores

The interactions of proteins with the surface of cylindrical nanopores are systematically investigated to elucidate how surface curvature and surface chemistry affect the conformation and activity of confined proteins in an aqueous, buffered environment. Two globular proteins, lysozyme and myoglobin, with different catalytic functions, were used as model proteins to analyze structural changes in proteins after adsorption on ordered mesoporous silica SBA-15 and propyl-functionalized SBA-15 (C(3)SBA-15) with carefully controlled pore size. Liquid phase ATR-FTIR spectroscopy was used to study the amide I and II bands of the adsorbed proteins. The amide I bands showed that the secondary structures of free and adsorbed protein molecules differ, and that the secondary structure of the adsorbed protein is influenced by the local geometry as well as by the surface chemistry of the nanopores. The conformation of the adsorbed proteins inside the nanopores of SBA-15 and C(3)SBA-15 is strongly correlated with the local geometry and the surface properties of the nanoporous materials, which results in different catalytic activities. Adsorption by electrostatic interaction of proteins in nanopores of an optimal size provides a favorably confining and protecting environment, which may lead to considerably enhanced structural stability and catalytic activity.

A demonstration of the inhomogeneity of the local dielectric response of proteins by molecular dynamics simulations

The dielectric properties of proteins are important in determining responses to changes in pH, responses to external electric fields, and the energy associated with charge injection or charge transfer. Although the importance of electrostatic interactions in proteins is well established and they are routinely included in simulations, the dielectric properties of proteins are less frequently discussed. We use a combination of atomistic molecular dynamics simulation and continuum electrostatics models to calculate the local dielectric constant within regions of the model protein hen egg white lysozyme. The local dielectric response of the protein was found to be highly inhomogeneous and strongly dependent on the local chemical environment. In light of these findings, we discuss the use of models using a single effective dielectric constant for the whole protein and the relationship between variations in local protein dielectric response and the activity of protein sites.

Residue-Specific pKa Determination of Lysine and Arginine Side Chains by Indirect 15N and 13C NMR Spectroscopy: Application to apo Calmodulin

Electrostatic interactions in proteins can be probed experimentally through determination of residue-specific acidity constants. We describe here triple-resonance NMR techniques for direct determination of lysine and arginine side-chain protonation states in proteins. The experiments are based on detection of nonexchangeable protons over the full range of pH and temperature and therefore are well suited for pKa determination of individual amino acid side chains. The experiments follow the side-chain 15Nzeta (lysine) and 15Nepsilon or 13Czeta (arginine) chemical shift, which changes due to sizable changes in the heteronuclear electron distribution upon (de)protonation. Since heteronuclear chemical shifts are overwhelmed by the charge state of the amino acid side chain itself, these methods supersede 1H-based NMR in terms of accuracy, sensitivity, and selectivity. Moreover, the 15Nzeta and 15Nepsilon nuclei may be used to probe changes in the local electrostatic environment. Applications to three proteins are described: apo calmodulin, calbindin D9k, and FKBP12. For apo calmodulin, residue-specific pKa values of lysine side chains were determined to fall between 10.7 and 11.2 as a result of the high net negative charge on the protein surface. Ideal two-state titration behavior observed for all lysines indicates the absence of significant direct charge interactions between the basic residues. These results are compared with earlier studies based on chemical modification.

Peptide 17O Chemical Shielding and Electric Field Gradient Tensors

Complete (17)O chemical shielding (CS) and quadrupole coupling (QC) tensors and their molecular orientations were determined for the central residues in two tripeptides Gly-Gly-Val (GGV) and Ala-Gly-Gly (AGG) by single-crystal NMR methods. Tensor orientations in the two peptides are very similar, however, principal components are different. The most shielded CS and smallest magnitude QC components are normal to the peptide plane, while the most deshielded CS and largest QC components are in the peptide plane either at an angle of 17 degrees (CS) or perpendicular (QC) to the C=O bond. Comparisons of principal components from experiment and DFT calculations indicate that the smaller shielding tensor span in GGV (549 ppm) compared to AGG (606 ppm) is likely due to two factors: a shorter "direct" H-bond distance to the peptide carbonyl oxygen and an "indirect" H bond of the peptide NH to a carboxylate rather than a carbonyl. We anticipate that (17)O NMR should be generally useful for probing H-bonding and local electrostatic interactions in proteins and polypeptides. Using the single-crystal data as an accurate reference, we show that a useful subset of the NMR parameters, QC and CS principal components and their relative orientation, can be obtained with reasonable accuracy from a very high-field (21.2 T), stationary sample powder spectrum.

Stabilization of Model Beverage Cloud Emulsions Using Protein−Polysaccharide Electrostatic Complexes Formed at the Oil−Water Interface

The potential of utilizing interfacial complexes, formed through the electrostatic interactions of proteins and polysaccharides at oil-water interfaces, to stabilize model beverage cloud emulsions has been examined. These interfacial complexes were formed by mixing charged polysaccharides with oil-in-water emulsions containing oppositely charged protein-coated oil droplets. Model beverage emulsions were prepared that consisted of 0.1 wt % corn oil droplets coated by beta-lactoglobulin (beta-Lg), beta-Lg/alginate, beta-Lg/iota-carrageenan, or beta-Lg/gum arabic interfacial layers (pH 3 or 4). Stable emulsions were formed when the polysaccharide concentration was sufficient to saturate the protein-coated droplets. The emulsions were subjected to variations in pH (from 3 to 7), ionic strength (from 0 to 250 mM NaCl), and thermal processing (from 30 or 90 degrees C), and the influence on their stability was determined. The emulsions containing alginate and carrageenan had the best stability to ionic strength and thermal processing. This study shows that the controlled formation of protein-polysaccharide complexes at droplet surfaces may be used to produce stable beverage emulsions, which may have important implications for industrial applications.

On the factors affecting the enhanced resonance light scattering signals of the interactions between proteins and multiply negatively charged chromophores using water blue as an example

Electrostatic interactions of proteins, including bovine serum albumin (BSA), human serum albumin (HSA), γ-globulin (γ-IgG), α-chymotrypsin (Chy), lysozyme (Lys) and cellulase (Cel), with multiply negatively charged chromophores were investigated based on the measurements of the enhanced resonance light scattering (RLS) signals. Using triply negatively charged water blue (WB) as an example, the factors were discussed that affect the enhanced resonance light scattering signals of the interactions between proteins and the negatively charged chromophores. It was found that the enhanced RLS signals with the maximum light scattering peak at 346.0 nm in these interacting systems are strongly dependent on the isoelectric points of proteins and show adverse linear relationships with increasing ionic strength depending on the positive charges of the inorganic metal ions used to control the ionic strength of the medium, sufficiently disclosing that the electrostatic attraction performs an important role in the combination of proteins with WB. Linear responses were discovered between the enhanced RLS signals and the protein molecular weights ( M w), displaying the dimensions of scattered particles formed by proteins and WB make a key contribution to the RLS enhancements. An empirical equation is proposed which possibly displays the factors affecting the enhanced RLS signals of the interactions between proteins with negatively charged chromophores.

Influence of the Solvent Structure on the Electrostatic Interactions in Proteins

The proper estimation of the influence of the many-body dynamic solvent microstructure on a pairwise electrostatic interaction (PEI) at the protein-solvent interface is very important for solving many biophysical problems. In this work, the PEI energy was calculated for a system that models the interface between a protein and an aqueous solvent. The concept of nonlocal electrostatics for interfacial electrochemical systems was used to evaluate the contribution of a solvent orientational polarization, correlated by the network of hydrogen bonds, into the PEI energy in proteins. The analytical expression for this energy was obtained in the form of Coulomb's law with an effective distance-dependent dielectric function. The asymptotic and numerical analysis carried out for this function revealed several features of dielectric heterogeneity at the protein-solvent interface. For charges located in close proximity to this interface, the values of the dielectric function for the short-distance electrostatic interactions were found to be remarkably smaller than those determined by the classical model, in which the solvent was considered as the uniform dielectric medium of high dielectric constant. Our results have shown that taking into consideration the dynamic solvent microstructure remarkably increases the value of the PEI energy at the protein-solvent interface.

Using charge ladders and capillary electrophoresis to measure the charge, size, and electrostatic interactions of proteins.

This chapter provides an overview of protein charge ladders--collections of protein derivatives that differ in charge--and capillary electrophoresis (CE). The combination of charge ladders and CE is a useful biophysical tool for measuring the net charge of proteins and the role of electrostatics in biochemical processes involving proteins. Methods to synthesize and analyze charge ladders by CE are described. Applications of charge ladders and CE to the simultaneous measurement of net charge and hydrodynamic radius of proteins are presented. Techniques for using charge ladders and CE to measure the role of interactions between charged groups on protein stability and ligand binding are also given. The power of this approach lies in the ability to isolate protein charge as an independent and measurable variable in the study of protein stability and function.

An empirical model for electrostatic interactions in proteins incorporating multiple geometry-dependent dielectric constants.

Here we introduce an electrostatic model that treats the complexity of electrostatic interactions in a heterogeneous protein environment by using multiple parameters that take into account variations in protein geometry, local structure, and the type of interacting residues. The optimal values for these parameters were obtained by fitting the model to a large dataset of 260 experimentally determined pK(a) values distributed over 41 proteins. We obtain fits between the calculated and observed values that are significantly better than the null model. The model performs well on the groups that exhibit large pK(a) shifts from solution values in response to the protein environment and compares favorably with other, successful continuum models. The empirically determined values of the parameters correlate well with experimentally observed contributions of hydrogen bonds and ion pairs as well as theoretically predicted magnitudes of charge-charge and charge-polar interactions. The magnitudes of the dielectric constants assigned to different regions of the protein rank according to the strength of the relaxation effects expected for the core, boundary, and surface. The electrostatic interactions in this model are pairwise decomposable and can be calculated rapidly. This model is therefore well suited for the large computations required for simulating protein properties and especially for prediction of mutations for protein design.

Evaluation of Models of Electrostatic Interactions in Proteins

The conformations of proteins and protein−protein complexes observed in nature must be low in free energy relative to alternative (not observed) conformations, and it is plausible (but not absolutely necessary) that the electrostatic free energies of experimentally observed conformations are also low relative to other conformations. Starting from this assumption, we evaluate alternative models of electrostatic interactions in proteins by comparing the electrostatic free energies of native, nativelike, and non-native structures. We observe that the total electrostatic free energy computed using the Poisson−Boltzmann (PB) equation or the generalized Born (GB) model exhibits free energy gaps that are comparable to, or smaller than, the free energy gaps resulting from Coulomb interactions alone. Detailed characterization of the contributions of different atom types to the total electrostatic free energy showed that, although for most atoms unfavorable solvation energies associated with atom burial are more than compensated by attractive Coulomb interactions, Coulomb interactions do not become more favorable with burial for certain backbone atom types, suggesting inaccuracies in the treatment of backbone electrostatics. Sizable free energy gaps are obtained using simple distance-dependent dielectric models, suggesting their usefulness in approximating the attenuation of long range Coulomb interactions by induced polarization effects. Hydrogen bonding interactions appear to be better modeled with an explicitly orientation-dependent hydrogen bonding potential than with any of the purely electrostatic models of hydrogen bonds, as there are larger free energy gaps with the former. Finally, a combined electrostatics−hydrogen bonding potential is developed that appears to better capture the free energy differences between native, nativelike, and non-native proteins and protein−protein complexes than electrostatic or hydrogen bonding models alone.

Relationship between Ion Pair Geometries and Electrostatic Strengths in Proteins

The electrostatic free energy contribution of an ion pair in a protein depends on two factors, geometrical orientation of the side-chain charged groups with respect to each other and the structural context of the ion pair in the protein. Conformers in NMR ensembles enable studies of the relationship between geometry and electrostatic strengths of ion pairs, because the protein structural contexts are highly similar across different conformers. We have studied this relationship using a dataset of 22 unique ion pairs in 14 NMR conformer ensembles for 11 nonhomologous proteins. In different NMR conformers, the ion pairs are classified as salt bridges, nitrogen–oxygen (N O) bridges and longer-range ion pairs on the basis of geometrical criteria. In salt bridges, centroids of the side-chain charged groups and at least a pair of side-chain nitrogen and oxygen atoms of the ion-pairing residues are within a 4 Å distance. In N O bridges, at least a pair of the side-chain nitrogen and oxygen atoms of the ion-pairing residues are within 4 Å distance, but the distance between the side-chain charged group centroids is greater than 4 Å. In the longer-range ion pairs, the side-chain charged group centroids as well as the side-chain nitrogen and oxygen atoms are more than 4 Å apart. Continuum electrostatic calculations indicate that most of the ion pairs have stabilizing electrostatic contributions when their side-chain charged group centroids are within 5 Å distance. Hence, most (∼92%) of the salt bridges and a majority (68%) of the N O bridges are stabilizing. Most (∼89%) of the destabilizing ion pairs are the longer-range ion pairs. In the NMR conformer ensembles, the electrostatic interaction between side-chain charged groups of the ion-pairing residues is the strongest for salt bridges, considerably weaker for N O bridges, and the weakest for longer-range ion pairs. These results suggest empirical rules for stabilizing electrostatic interactions in proteins.

Comparison of all atom, continuum, and linear fitting empirical models for charge screening effect of aqueous medium surrounding a protein molecule

To investigate the role hydration plays in the electrostatic interactions of proteins, the time-averaged electrostatic potential of the B1 domain of protein G in an aqueous solution was calculated with full atomic molecular dynamics simulations that explicitly considers every atom (i.e., an all atom model). This all atom calculated potential was compared with the potential obtained from an electrostatic continuum model calculation. In both cases, the charge-screening effect was fairly well formulated with an effective relative dielectric constant which increased linearly with increasing charge–charge distance. This simulated linear dependence agrees with the experimentally determined linear relation proposed by Pickersgill. Cut-off approximations for Coulomb interactions failed to reproduce this linear relation. Correlation between the all atom model and the continuum models was found to be better than the respective correlation calculated for linear fitting to the two models. This confirms that the continuum model is better at treating the complicated shapes of protein conformations than the simple linear fitting empirical model. We have tried a sigmoid fitting empirical model in addition to the linear one. When weights of all data were treated equally, the sigmoid model, which requires two fitting parameters, fits results of both the all atom and the continuum models less accurately than the linear model which requires only one fitting parameter. When potential values are chosen as weighting factors, the fitting error of the sigmoid model became smaller, and the slope of both linear fitting curves became smaller. This suggests the screening effect of an aqueous medium within a short range, where potential values are relatively large, is smaller than that expected from the linear fitting curve whose slope is almost 4. To investigate the linear increase of the effective relative dielectric constant, the Poisson equation of a low-dielectric sphere in a high-dielectric medium was solved and charges distributed near the molecular surface were indicated as leading to the apparent linearity.

Electrostatic effects in highly charged proteins: salt sensitivity of pKa values of histidines in staphylococcal nuclease.

The pK(a) values of most histidines in small peptides and in myoglobin increase on average by 0.30 unit between 0.02 and 1.5 M NaCl [Kao et al. (2000) Biophys. J. 79, 1637]. The DeltapK(a) values reflect primarily the ionic strength dependence of the solvation energy; screening of Coulombic interactions contributes only in a minor way. This implies that Coulombic interactions are weak, or that attractive and repulsive contributions to the pK(a) values are balanced. To distinguish experimentally between these two possibilities, and to further characterize the magnitude and salt sensitivity of surface electrostatic interactions in proteins, the salt dependence of pK(a) values of histidines in staphylococcal nuclease was measured by (1)H NMR spectroscopy. Three of the four histidines titrated with significantly depressed pK(a) values, and the salt sensitivity of all histidine pK(a) values was substantial. In three cases, the pK(a) values increased by a full unit between 0.01 and 1.5 M KCl. Anion-specific effects were found; the pK(a) values measured under equivalent ionic strengths in SCN(-) and SO(4)(2-) were higher than in Cl(-); the order of the sensitivity of pK(a) values to anions was SCN(-) > Cl(-) > SO(4)(2-). Structure-based pK(a) calculations with continuum methods were performed to interpret the measured effects structurally and to test their ability to capture the experimental behavior. Calculations in which the protein interior was treated empirically with a dielectric constant of 20 reproduced the pK(a) values and their dependence on the concentration of Cl(-). According to the calculations, the pK(a) values are depressed because of unfavorable self-energies and repulsive Coulombic interactions. Their striking salt sensitivity reflects screening of weak, repulsive, Coulombic interactions among charges separated by more than 10 A. Long-range Coulombic interactions on the surfaces of proteins are weak, but they can add up to produce substantial electrostatic effects when positive and negative charges are not balanced.