A compact RNA tertiary structure contains a buried backbone-K+ complex

The structure of a 58 nucleotide ribosomal RNA fragment buries several phosphate groups of a hairpin loop within a large tertiary core. During refinement of an X-ray crystal structure containing this RNA, a potassium ion was found to be contacted by six oxygen atoms from the buried phosphate groups; the ion is contained completely within the solvent-accessible surface of the RNA. The electrostatic potential at the ion chelation site is unusually large, and more than compensates for the substantial energetic penalties associated with partial dehydration of the ion and displacement of delocalized ions. The very large predicted binding free energy, approximately -30 kcal/mol, implies that the site must be occupied for the RNA to fold. These findings agree with previous studies of the ion-dependent folding of tertiary structure in this RNA, which concluded that a monovalent ion was bound in a partially dehydrated environment where Mg2+ could not easily compete for binding. By compensating the unfavorable free energy of buried phosphate groups with a chelated ion, the RNA is able to create a larger and more complex tertiary fold than would be possible otherwise.     

Electrostatic interactions in the GCN4 leucine zipper: substantial contributions arise from intramolecular interactions enhanced on binding.

The GCN4 leucine zipper is a peptide homodimer that has been the subject of a number of experimental and theoretical investigations into the determinants of affinity and specificity. Here, we utilize this model system to investigate electrostatic effects in protein binding using continuum calculations. A particularly novel feature of the computations made here is that they provide an interaction-by-interaction breakdown of the electrostatic contributions to the free energy of docking that includes changes in the interaction of each functional group with solvent and changes in interactions between all pairs of functional groups on binding. The results show that (1) electrostatic effects disfavor binding by roughly 15 kcal/mol due to desolvation effects that are incompletely compensated in the bound state, (2) while no groups strongly stabilize binding, the groups that are most destabilizing are charged and polar side chains at the interface that have been implicated in determining binding specificity, and (3) attractive intramolecular interactions (e.g., backbone hydrogen bonds) that are enhanced on binding due to reduced solvent screening in the bound state contribute significantly to affinity and are likely to be a general effect in other complexes. A comparison is made between the results obtained in an electrostatic analysis carried out calculationally and simulated results corresponding to idealized data from a scanning mutagenesis experiment. It is shown that scanning experiments provide incomplete information on interactions and, if overinterpreted, tend to overestimate the energetic effect of individual side chains that make attractive interactions. Finally, a comparison is made between the results available from a continuum electrostatic model and from a simpler surface-area dependent solvation model. In this case, although the simpler model neglects certain interactions, on average it performs rather well.     

Energetics of charge-charge interactions in proteins

Electrostatic interactions between pairs of atoms in proteins are calculated with a model based on the linearized Poisson-Boltzmann equation. The equation is solved accurately by a method that takes into account the detailed shape of the protein. This paper presents applications to several systems. Experimental data for the interaction of ionized residues with an active site histidine in subtilisin BPN' allow the model to be tested, using various assumptions for the electrical properties of the protein and solvent. The electrostatic stabilization of the active site thiolate of rhodanese is analyzed, with attention to the influence of alpha-helices. Finally, relationships between electrostatic potential and charge-charge distance are reported for large and small globular proteins. The above results are compared with those of simpler electrostatic models, including Coulomb's law with both a distance-dependent dielectric constant (epsilon = R) and a fixed dielectric constant (epsilon = 2), and Tanford-Kirkwood theory. The primary conclusions are as follows: 1) The Poisson-Boltzmann model agrees with the subtilisin data over a range of ionic strengths; 2) two alpha-helices generate a large potential in the active site of rhodanese; 3) epsilon = R overestimates weak electrostatic interactions but yields relatively good results for strong ones; 4) Tanford-Kirkwood theory is a useful approximation to detailed solutions of the linearized Poisson-Boltzmann equation in globular proteins; and 5) the modified Tanford-Kirkwood theory over-screens the measured electrostatic interactions in subtilisin.     

On the calculation of binding free energies using continuum methods: Application to MHC class I protein-peptide interactions

This paper describes a methodology to calculate the binding free energy (ΔG) of a protein-ligand complex using a continuum model of the solvent. A formal thermodynamic cycle is used to decompose the binding free energy into electrostatic and non-electrostatic contributions. In this cycle, the reactants are discharged in water, associated as purely nonpolar entities, and the final complex is then recharged. The total electrostatic free energies of the protein, the ligand, and the complex in water are calculated with the finite difference Poisson-Boltzmann (FDPB) method. The nonpolar (hydrophobic) binding free energy is calculated using a free energy-surface area relationship, with a single alkane/water surface tension coefficient (γ_aw). The loss in backbone and side-chain configurational entropy upon binding is estimated and added to the electrostatic and the nonpolar components of ΔG. The methodology is applied to the binding of the murine MHC class I protein H-2K^b with three distinct peptides, and to the human MHC class I protein HLA-A2 in complex with five different peptides. Despite significant differences in the amino acid sequences of the different peptides, the experimental binding free energy differences (ΔΔG_exp) are quite small (<0.3 and <2.7 kcal/mol for the H-2K^b and HLA-A2 complexes, respectively). For each protein, the calculations are successful in reproducing a fairly small range of values for ΔΔG_calc (<4.4 and <5.2 kcal/mol, respectively) although the relative peptide binding affinities of H-2K^b and HLA-A2 are not reproduced. For all protein-peptide complexes that were treated, it was found that electrostatic interactions oppose binding whereas nonpolar interactions drive complex formation. The two types of interactions appear to be correlated in that larger nonpolar contributions to binding are generally opposed by increased electrostatic contributions favoring dissociation. The factors that drive the binding of peptides to MHC proteins are discussed in light of our results.     

Do salt bridges stabilize proteins? A continuum electrostatic analysis

The electrostatic contribution to the free energy of folding was calculated for 21 salt bridges in 9 protein X-ray crystal structures using a continuum electrostatic approach with the DELPHI computer-program package. The majority (17) were found to be electrostatically destabilizing; the average free energy change, which is analogous to mutation of salt bridging side chains to hydrophobic isosteres, was calculated to be 3.5 kcal/mol. This is fundamentally different from stability measurements using pKa shifts, which effectively measure the strength of a salt bridge relative to 1 or more charged hydrogen bonds. The calculated effect was due to a large, unfavorable desolvation contribution that was not fully compensated by favorable interactions within the salt bridge and between salt-bridge partners and other polar and charged groups in the folded protein. Some of the salt bridges were studied in further detail to determine the effect of the choice of values for atomic radii, internal protein dielectric constant, and ionic strength used in the calculations. Increased ionic strength resulted in little or no change in calculated stability for 3 of 4 salt bridges over a range of 0.1-0.9 M. The results suggest that mutation of salt bridges, particularly those that are buried, to "hydrophobic bridges" (that pack at least as well as wild type) can result in proteins with increased stability. Due to the large penalty for burying uncompensated ionizable groups, salt bridges could help to limit the number of low free energy conformations of a molecule or complex and thus play a role in determining specificity (i.e., the uniqueness of a protein fold or protein-ligand binding geometry).     

Electrostatic contribution to the binding stability of protein-protein complexes

To investigate roles of electrostatic interactions in protein binding stability, electrostatic calculations were carried out on a set of 64 mutations over six protein-protein complexes. These mutations alter polar interactions across the interface and were selected for putative dominance of electrostatic contributions to the binding stability. Three protocols of implementing the Poisson-Boltzmann model were tested. In vdW4 the dielectric boundary between the protein low dielectric and the solvent high dielectric is defined as the protein van der Waals surface and the protein dielectric constant is set to 4. In SE4 and SE20, the dielectric boundary is defined as the surface of the protein interior inaccessible to a 1.4-A solvent probe, and the protein dielectric constant is set to 4 and 20, respectively. In line with earlier studies on the barnase-barstar complex, the vdW4 results on the large set of mutations showed the closest agreement with experimental data. The agreement between vdW4 and experiment supports the contention of dominant electrostatic contributions for the mutations, but their differences also suggest van der Waals and hydrophobic contributions. The results presented here will serve as a guide for future refinement in electrostatic calculation and inclusion of nonelectrostatic effects.     

Electrostatic interactions play an essential role in DNA repair and cold-adaptation of Uracil DNA glycosylase

Life has adapted to most environments on earth, including low and high temperature niches. The increased catalytic efficiency and thermoliability observed for enzymes from organisms living in constantly cold regions when compared to their mesophilic and thermophilic cousins are poorly understood at the molecular level. Uracil DNA glycosylase (UNG) from cod (cUNG) catalyzes removal of uracil from DNA with an increased k_cat and reduced K_m relative to its warm-active human (hUNG) counterpart. Specific issues related to DNA repair and substrate binding/recognition (K_m) are here investigated by continuum electrostatics calculations, MD simulations and free energy calculations. Continuum electrostatic calculations reveal that cUNG has surface potentials that are more complementary to the DNA potential at and around the catalytic site when compared to hUNG, indicating improved substrate binding. Comparative MD simulations combined with free energy calculations using the molecular mechanics-Poisson Boltzmann surface area (MM-PBSA) method show that large opposing energies are involved when forming the enzyme-substrate complexes. Furthermore, the binding free energies obtained reveal that the Michaelis-Menten complex is more stable for cUNG, primarily due to enhanced electrostatic properties, suggesting that energetic fine-tuning of electrostatics can be utilized for enzymatic temperature adaptation. Energy decomposition pinpoints the residual determinants responsible for this adaptation. Figure Electrostatic isosurfaces of cod uracil DNA glycosylase in complex with double stranded DNA     

Distance dependence and salt sensitivity of pairwise, coulombic interactions in a protein

Histidine pK(a) values were measured in charge-reversal (K78E, K97E, K127E, and K97E/K127E) and charge-neutralization (E10A, E101A, and R35A) mutants of staphylococcal nuclease (SNase) by (1)H-NMR spectroscopy. Energies of interaction between pairs of charges (DeltaG(ij)) were obtained from the shifts in pK(a) values relative to wild-type values. The data describe the distance dependence and salt sensitivity of pairwise coulombic interactions. Calculations with a continuum electrostatics method captured the experimental DeltaG(ij) when static structures were used and when the protein interior was treated empirically with a dielectric constant of 20. The DeltaG(ij) when r(ij) < or = 10 A were exaggerated slightly in the calculations. Coulomb's law with a dielectric constant near 80 and a Debye-Hückel term to account for screening by the ionic strength reproduced the salt sensitivity and distance dependence of DeltaG(ij) as well as the structure-based method. In their interactions with each other, surface charges behave as if immersed in water; the Debye length describes realistically the distance where interactions become negligible at a given ionic strength. On average, charges separated by distances (r(ij)) approximately 5 A interacted with DeltaG(ij) approximately 0.6 kcal/mole in 0.01 M KCl, but DeltaG(ij) decayed to < or =0.10 kcal/mole when r(ij) = 20 A. In 0.10 M KCl, DeltaG(ij) approximately 0.10 kcal/mole when r(ij) = 10 A. In 1.5 M KCl, only short-range interactions with r(ij) < or = 5 A persisted. Although at physiological ionic strengths the interactions between charges separated by more than 10 A are extremely weak, in situations where charge imbalance exists many weak interactions can cumulatively produce substantial effects.     

Hydrophobic vs coulombic interactions in the binding of steroidal polyamines to DNA

Seven new steroidal polyamines derived from bile acids, either lithocholic or deoxycholic acid, have been studied as DNA-binding agents using four complimentary methods: an ethidium displacement assay, observed changes in the thermal denaturation of poly[d(AT)], effects on hyperchromicity of DNA, and circular dichroism. In addition, modelling studies were conducted to examine the electrostatic surface potential of the polycations. The results point to a key role for a large hydrophobic surface area on the steroid in addition to the Coulombic attraction by ammonium and guanidinium groups on the steroid interacting with the polyphosphate backbone.     

Three-dimensional structure of a fluorescein-Fab complex crystallized in 2-methyl-2,4-pentanediol

The crystal structure of a fluorescein-Fab (4-4-20) complex was determined at 2.7 A resolution by molecular replacement methods. The starting model was the refined 2.7 A structure of unliganded Fab from an autoantibody (BV04-01) with specificity for single-stranded DNA. In the 4-4-20 complex fluorescein fits tightly into a relatively deep slot formed by a network of tryptophan and tyrosine side chains. The planar xanthonyl ring of the hapten is accommodated at the bottom of the slot while the phenylcarboxyl group interfaces with solvent. Tyrosine 37 (light chain) and tryptophan 33 (heavy chain) flank the xanthonyl group and tryptophan 101 (light chain) provides the floor of the combining site. Tyrosine 103 (heavy chain) is situated near the phenyl ring of the hapten and tyrosine 102 (heavy chain) forms part of the boundary of the slot. Histidine 31 and arginine 39 of the light chain are located in positions adjacent to the two enolic groups at opposite ends of the xanthonyl ring, and thus account for neutralization of one of two negative charges in the haptenic dianion. Formation of an enol-arginine ion pair in a region of low dielectric constant may account for an incremental increase in affinity of 2-3 orders of magnitude in the 4-4-20 molecule relative to other members of an idiotypic family of monoclonal antifluorescyl antibodies. The phenyl carboxyl group of fluorescein appears to be hydrogen bonded to the phenolic hydroxyl group of tyrosine 37 of the light chain. A molecule of 2-methyl-2,4-pentanediol (MPD), trapped in the interface of the variable domains just below the fluorescein binding site, may be partly responsible for the decrease in affinity for the hapten in MPD.     

Electrostatics of the FeS clusters in respiratory complex I

Respiratory complex I couples the transfer of electrons from NADH to ubiquinone and the translocation of protons across the mitochondrial membrane. A detailed understanding of the midpoint reduction potentials (E_m) of each redox center and the factors which influence those potentials are critical in the elucidation of the mechanism of electron transfer in this enzyme. We present accurate electrostatic interaction energies for the iron-sulfur (FeS) clusters of complex I to facilitate the development of models and the interpretation of experiments in connection to electron transfer (ET) in this enzyme. To calculate redox titration curves for the FeS clusters it is necessary to include interactions between clusters, which in turn can be used to refine E_m values and validate spectroscopic assignments of each cluster. Calculated titration curves for clusters N4, N5, and N6a are discussed. Furthermore, we present some initial findings on the electrostatics of the redox centers of complex I under the influence of externally applied membrane potentials. A means of determining the location of the FeS cofactors within the holo-complex based on electrostatic arguments is proposed. A simple electrostatic model of the protein/membrane system is examined to illustrate the viability of our hypothesis.     

基于分子模拟的离子交换层析中静电相互作用研究

将分子模拟方法引入到蛋白质离子交换层析中的静电相互作用研究。选用蛋清溶菌酶和牛胰凝乳蛋白酶为模型蛋白质,阳离子交换吸附剂SP Sepharose FF等为模型层析介质。从蛋白质数据库(PDB)中获得蛋白质三维结构数据,分析了介质孔径和配基分布,以点电荷模拟离子交换层析介质的功能配基,构筑了蛋白质-介质配基模拟表面体系。采用MCCE、Delphi和GRASP等程序包进行了分子模拟计算,考察了作用方向、作用距离、盐浓度、pH等对蛋白质和模拟配基平面间静电相互作用的影响。结果表明,宏观的层析平衡常数与微观分子模拟计算的相互作用能量参数间存在良好的线性关系。     

Molecular dynamics-solvated interaction energy studies of protein-protein interactions: the MP1-p14 scaffolding complex

Using the MP1-p14 scaffolding complex from the mitogen-activated protein kinase signaling pathway as model system, we explored a structure-based computational protocol to probe and characterize binding affinity hot spots at protein-protein interfaces. Hot spots are located by virtual alanine-scanning consensus predictions over three different energy functions and two different single-structure representations of the complex. Refined binding affinity predictions for select hot-spot mutations are carried out by applying first-principle methods such as the molecular mechanics generalized Born surface area (MM-GBSA) and solvated interaction energy (SIE) to the molecular dynamics (MD) trajectories for mutated and wild-type complexes. Here, predicted hot-spot residues were actually mutated to alanine, and crystal structures of the mutated complexes were determined. Two mutated MP1-p14 complexes were investigated, the p14(Y56A)-mutated complex and the MP1(L63A,L65A)-mutated complex. Alternative ways to generate MD ensembles for mutant complexes, not relying on crystal structures for mutated complexes, were also investigated. The SIE function, fitted on protein-ligand binding affinities, gave absolute binding affinity predictions in excellent agreement with experiment and outperformed standard MM-GBSA predictions when tested on the MD ensembles of Ras-Raf and Ras-RalGDS protein-protein complexes. For wild-type and mutant MP1-p14 complexes, SIE predictions of relative binding affinities were supported by a yeast two-hybrid assay that provided semiquantitative relative interaction strengths. Results on the MP1-mutated complex suggested that SIE predictions deteriorate if mutant MD ensembles are approximated by just mutating the wild-type MD trajectory. The SIE data on the p14-mutated complex indicated feasibility for generating mutant MD ensembles from mutated wild-type crystal structure, despite local structural differences observed upon mutation. For energetic considerations, this would circumvent costly needs to produce and crystallize mutated complexes. The sensitized protein-protein interface afforded by the p14(Y56A) mutation identified here has practical applications in screening-based discovery of first-generation small-molecule hits for further development into specific modulators of the mitogen-activated protein kinase signaling pathway.     

Electrostatic interactions contribute to reduced heat capacity change of unfolding in a thermophilic ribosomal protein l30e

The origin of reduced heat capacity change of unfolding (DeltaC(p)) commonly observed in thermophilic proteins is controversial. The established theory that DeltaC(p) is correlated with change of solvent-accessible surface area cannot account for the large differences in DeltaC(p) observed for thermophilic and mesophilic homologous proteins, which are very similar in structures. We have determined the protein stability curves, which describe the temperature dependency of the free energy change of unfolding, for a thermophilic ribosomal protein L30e from Thermococcus celer, and its mesophilic homologue from yeast. Values of DeltaC(p), obtained by fitting the free energy change of unfolding to the Gibbs-Helmholtz equation, were 5.3 kJ mol(-1) K(-1) and 10.5 kJ mol(-1) K(-1) for T.celer and yeast L30e, respectively. We have created six charge-to-neutral mutants of T.celer L30e. Removal of charges at Glu6, Lys9, and Arg92 decreased the melting temperatures of T.celer L30e by approximately 3-9 degrees C, and the differences in melting temperatures were smaller with increasing concentration of salt. These results suggest that these mutations destabilize T.celer L30e by disrupting favorable electrostatic interactions. To determine whether electrostatic interactions contribute to the reduced DeltaC(p) of the thermophilic protein, we have determined DeltaC(p) for wild-type and mutant T.celer L30e by Gibbs-Helmholtz and by van't Hoff analyses. A concomitant increase in DeltaC(p) was observed for those charge-to-neutral mutants that destabilize T.celer L30e by removing favorable electrostatic interactions. The crystal structures of K9A, E90A, and R92A, were determined, and no structural change was observed. Taken together, our results support the conclusion that electrostatic interactions contribute to the reduced DeltaC(p) of T.celer L30e.     

Energetic evaluation of binding modes in the C3d and Factor H (CCP 19-20) complex

As a part of innate immunity, the complement system relies on activation of the alternative pathway (AP). While feed-forward amplification generates an immune response towards foreign surfaces, the process requires regulation to prevent an immune response on the surface of host cells. Factor H (FH) is a complement protein secreted by native cells to negatively regulate the AP. In terms of structure, FH is composed of 20 complement-control protein (CCP) modules that are structurally homologous but vary in composition and function. Mutations in these CCPs have been linked to states of autoimmunity. In particular, several mutations in CCP 19-20 are correlated to atypical hemolytic uremic syndrome (aHUS). From crystallographic structures there are three putative binding sites of CCP 19-20 on C3d. Since there has been some controversy over the primary mode of binding from experimental studies, we approach characterization of binding using computational methods. Specifically, we compare each binding mode in terms of electrostatic character, structural stability, dissociative and associative properties, and predicted free energy of binding. After a detailed investigation, we found two of the three binding sites to be similarly stable while varying in the number of contacts to C3d and in the energetic barrier to complex dissociation. These sites are likely physiologically relevant and may facilitate multivalent binding of FH CCP 19-20 to C3b and either C3d or host glycosaminoglycans. We propose thermodynamically stable binding with modules 19 and 20, the latter driven by electrostatics, acting synergistically to increase the apparent affinity of FH for host surfaces.     

Electrostatic interactions play a significant role in the affinity of Saccharomyces cerevisiae phosphoenolpyruvate carboxykinase for Mn

Phosphoenolpyruvate (PEP) carboxykinases catalyse the reversible formation of oxaloacetate (OAA) and ATP (or GTP) from PEP, ADP (or GDP) and CO₂. They are activated by Mn²⁺, a metal ion that coordinates to the protein through the ?-amino group of a lysine residue, the N_?-2-imidazole of a histidine residue, and the carboxylate from an aspartic acid residue. Neutrality in the ?-amino group of Lys213 of Saccharomyces cerevisiae PEP carboxykinase is expected to be favoured by the vicinity of ionised Lys212. Glu272 and Glu284, located close to Lys212, should, in turn, electrostatically stabilise its positive charge and hence assist in keeping the ?-amino group of Lys213 in a neutral state. The mutations Glu272Gln, Glu284Gln, and Lys212Met increased the activation constant for Mn²⁺ in the main reaction of the enzyme up to seven-fold. The control mutation Lys213Gln increased this constant by ten-fold, as opposed to control mutation Lys212Arg, which did not affect the Mn²⁺ affinity of the enzyme. These observations indicate a role for Glu272, Glu284, and Lys212 in assisting Lys213 to properly bind Mn²⁺. In an unexpected result, the mutations Glu284Gln, Lys212Met and Lys213Gln changed the nucleotide-independent OAA decarboxylase activity of S. cerevisiae PEP carboxykinase into an ADP-requiring activity, implying an effect on the OAA binding characteristics of PEP carboxykinase.     

13C-n.m.r. studies of the conformational changes in proline oligomers brought about by lithium and calcium salts

A detailed ¹³C-n.m.r. investigation has been carried out on the conformational changes in proline oligomers brought about by interaction with lithium and calcium perchlorates. Interaction of lithium and calcium salts with Piv-(Pro)_n-OMe, n = 2, 4 and 5 results in trans-cis isomerization. In the case of pentaproline, metal salts also give rise to other trans-isomers caused by the rotation about the C^αC(O) bond (Ψ, cis). Calcium salts seem to stabilize cis'-isomers and produce effects somewhat different from those of lithium salts.     

Electrostatics of ion stabilization in a ClC chloride channel homologue from Escherichia coli

The structural determinants of electrostatics of ion stabilization within EcClC, a ClC-type chloride channel homologue from Escherichia coli, are studied using a continuum dielectric approximation. Specifically, the ion occupancy is investigated in the wild-type protein and a mutant thereof, and the contribution to the electrostatic binding free energy of local and non-local interactions is characterized at the single-residue level. This analysis shows that, in spite of the desolvation cost and the strong ion-ion repulsion, all previously reported binding sites can be occupied simultaneously. The stabilizing effect of the protein arises from hydrogen bonding as well as from longer-range favorable interactions, such as with the strictly conserved Lys131 side-chain. The latter is involved in the stabilization of the conserved GSGIP motif that delimits two of the binding sites. Interestingly, an additional low-affinity binding site, mediated by a structurally analogous motif including the side-chain of Arg340, can be identified on the extracellular side of the permeation pathway. Finally, it is shown that, in contrast to K-channels, and in analogy to the SBP/PBP sulfate/phosphate-binding proteins, the contribution of helix macrodipoles to chloride binding in EcClC is only marginal.