Electrostatic interactions in leucine zippers: thermodynamic analysis of the contributions of Glu and His residues and the effect of mutating salt bridges

Biomimetic Charge-Neutral Anion Receptors for Reversible Binding and Release of Highly Hydrated Phosphate in Water

Control of phosphate capture and release is vital in environmental, biological, and pharmaceutical contexts. However, the binding of trivalent phosphate (PO4 3-) in water is exceptionally difficult due to its high hydration energy. Based on the anion coordination chemistry of phosphate, in this study, four charge-neutral tripodal hexaurea receptors (L1-L4), which were equipped with morpholine and polyethylene glycol terminal groups to enhance their solubility in water, were synthesized to enable the pH-triggered phosphate binding and release in aqueous solutions. Encouragingly, the receptors were found to bind PO4 3- anion in a 1 : 1 ratio via hydrogen bonds in 100 % water solutions, with L1 exhibiting the highest binding constant (1.2×103 M-1). These represent the first neutral anion ligands to bind phosphate in 100 % water and demonstrate the potential for phosphate capture and release in water through pH-triggered mechanisms, mimicking native phosphate binding proteins. Furthermore, L1 can also bind multiple bioavailable phosphate species, which may serve as model systems for probing and modulating phosphate homeostasis in biological and biomedical researches.

Plausible blockers of Spike RBD in SARS-CoV2-molecular design and underlying interaction dynamics from high-level structural descriptors

Abstract

COVID-19 is characterized by an unprecedented abrupt increase in the viral transmission rate (SARS-CoV-2) relative to its pandemic evolutionary ancestor, SARS-CoV (2003). The complex molecular cascade of events related to the viral pathogenicity is triggered by the Spike protein upon interacting with the ACE2 receptor on human lung cells through its receptor binding domain (RBD_Spike). One potential therapeutic strategy to combat COVID-19 could thus be limiting the infection by blocking this key interaction. In this current study, we adopt a protein design approach to predict and propose non-virulent structural mimics of the RBD_Spike which can potentially serve as its competitive inhibitors in binding to ACE2. The RBD_Spike is an independently foldable protein domain, resilient to conformational changes upon mutations and therefore an attractive target for strategic re-design. Interestingly, in spite of displaying an optimal shape fit between their interacting surfaces (attributed to a consequently high mutual affinity), the RBD_Spike–ACE2 interaction appears to have a quasi-stable character due to a poor electrostatic match at their interface. Structural analyses of homologous protein complexes reveal that the ACE2 binding site of RBD_Spike has an unusually high degree of solvent-exposed hydrophobic residues, attributed to key evolutionary changes, making it inherently “reaction-prone.” The designed mimics aimed to block the viral entry by occupying the available binding sites on ACE2, are tested to have signatures of stable high-affinity binding with ACE2 (cross-validated by appropriate free energy estimates), overriding the native quasi-stable feature. The results show the apt of directly adapting natural examples in rational protein design, wherein, homology-based threading coupled with strategic “hydrophobic ↔ polar” mutations serve as a potential breakthrough.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1007/s00894-021-04779-0.

A possible founder mutation in FZD6 gene in a Turkish family with autosomal recessive nail dysplasia

BACKGROUND

Autosomal recessive nail dysplasia is characterized by thick and hard nails with a very slow growth on the hands and feet. Mutations in FZD6 gene were found to be associated with autosomal recessive nail dysplasia in 2011. Presently, only seven mutations have been reported in FZD6 gene; five mutations are clustered in the C-terminus, one is at the seventh transmembrane domain, and another is at the very beginning of third extracellular loop.

METHODS

Whole exome sequencing (WES) was applied to the index case, her one affected sister and her healthy consanguineous parents. The mutation was verified via Sanger sequencing. Molecular dynamics simulations of the predicted structures of native and mutant proteins were compared to gain insight into the pathogenicity mechanism of the mutation.

RESULTS

Here, we report a homozygous 8 bp deletion mutation, p.Gly559Aspfs*16; c.1676_1683delGAACCAGC, in FZD6 gene which causes a frameshift and creates a premature stop codon at position 16 of the new reading frame. Our molecular dynamics calculations predict that the pathogenicity of this frameshift mutation may be caused by the change in entropy of the protein with negative manner, disturbing the C-terminal domain structure, and hence interaction partners of FZD6.

CONCLUSION

We identified a homozygous deletion mutation in FZD6 in a consanguineous Turkish family with nail dysplasia. We also provide a molecular mechanism about the effects of the deletion on the protein structure and its possible motions. This study provides a pathogenicity mechanism for this mutation in nail dysplasia for the first time.

Conformational Dynamics of Asparagine at Coiled-Coil Interfaces

Coiled coils (CCs) are among the best-understood protein folds. Nonetheless, there are gaps in our knowledge of CCs. Notably, CCs are likely to be structurally more dynamic than often considered. Here, we explore this in an abundant class of CCs, parallel dimers, focusing on polar asparagine (Asn) residues in the hydrophobic interface. It is well documented that such inclusions discriminate between different CC oligomers, which has been rationalized in terms of whether the Asn can make side-chain hydrogen bonds. Analysis of parallel CC dimers in the Protein Data Bank reveals a variety of Asn side-chain conformations, but not all of these make the expected inter-side-chain hydrogen bond. We probe the structure and dynamics of a de novo-designed coiled-coil homodimer, CC-Di, by multidimensional nuclear magnetic resonance spectroscopy, including model-free dynamical analysis and relaxation-dispersion experiments. We find dynamic exchange on the millisecond time scale between Asn conformers with the side chains pointing into and out of the core. We perform molecular dynamics simulations that are consistent with this, revealing that the side chains are highly dynamic, exchanging between hydrogen-bonded-paired conformations in picoseconds to nanoseconds. Combined, our data present a more dynamic view for Asn at CC interfaces. Although inter-side-chain hydrogen bonding states are the most abundant, Asn is not always buried or engaged in such interactions. Because interfacial Asn residues are key design features for modulating CC stability and recognition, these further insights into how they are accommodated within CC structures will aid their predictive modeling, engineering, and design.

The accommodation index measures the perturbation associated with insertions and deletions in coiled-coils: Application to understand signaling in histidine kinases

Coiled-coils are essential components of many protein complexes. First discovered in structural proteins such as keratins, they have since been found to figure largely in the assembly and dynamics required for diverse functions, including membrane fusion, signal transduction and motors. Coiled-coils have a characteristic repeating seven-residue geometric and sequence motif, which is sometimes interrupted by the insertion of one or more residues. Such insertions are often highly conserved and critical to interdomain communication in signaling proteins such as bacterial histidine kinases. Here we develop the "accommodation index" as a parameter that allows automatic detection and classification of insertions based on the three dimensional structure of a protein. This method allows precise identification of the type of insertion and the "accommodation length" over which the insertion is structurally accommodated. A simple theory is presented that predicts the structural perturbations of 1, 3, 4 residue insertions as a function of the length over which the insertion is accommodated. Analysis of experimental structures is in good agreement with theory, and shows that short accommodation lengths give rise to greater perturbation of helix packing angles, changes in local helical phase, and increased structural asymmetry relative to long accommodation lengths. Cytoplasmic domains of histidine kinases in different signaling states display large changes in their accommodation lengths, which can now be seen to underlie diverse structural transitions including symmetry/asymmetry and local variations in helical phase that accompany signal transduction.

Effects of localized interactions and surface properties on stability of protein-based therapeutics

OBJECTIVES

Protein-based therapeutics garner significant attention because of exquisite specificity and limited side effects and are now being used to accomplish targeted delivery of small-molecule drugs. This review identifies and highlights individual chemical attributes and categorizes how site-specific changes affect protein stability based on published high-resolution molecular analyses.

KEY FINDINGS

Because it is challenging to determine the mechanisms by which the stability of large, complex molecules is altered and data are sparse, smaller, therapeutic proteins (insulin, erythropoietin, interferons) are examined alongside antibody data. Integrating this large pool of information with the limited available studies on antibodies reveals common mechanisms by which specific alterations affect protein structure and stability.

SUMMARY

Physical and chemical stability of therapeutic proteins and antibody drug conjugates (ADCs) is of critical importance because insufficient stability prevents molecules from making it to market. Individual moieties on/near the surface of proteins have substantial influence on structure and stability. Seemingly small, superficial modification may have far-reaching consequences on structure, conformational dynamics, and solubility of the protein, and hence physical stability of the molecule. Chemical modifications, whether spontaneous (e.g. oxidation, deamidation) or intentional, as with ADCs, may adversely impact stability by disrupting local surface properties or higher order protein structure.

Contributions of pocket depth and electrostatic interactions to affinity and selectivity of receptors for methylated lysine in water

Dynamic combinatorial chemistry was used to generate a set of receptors for peptides containing methylated lysine (KMen, n = 0-3) and study the contribution of electrostatic effects and pocket depth to binding affinity and selectivity. We found that changing the location of a carboxylate resulted in an increase in preference for KMe2, presumably based on ability to form a salt bridge with KMe2. The number of charged groups on either the receptor or peptide guest systematically varied the binding affinities to all guests by approximately 1-1.5 kcal mol(-1), with little influence on selectivity. Lastly, formation of a deeper pocket led to both increased affinity and selectivity for KMe3 over the lower methylation states. From these studies, we identified that the tightest binder was a receptor with greater net charge, with a Kd of 0.2 μM, and the receptor with the highest selectivity was the one with the deepest pocket, providing 14-fold selectivity between KMe3 and KMe2 and a Kd for KMe3 of 0.3 μM. This work provides key insights into approaches to improve binding affinity and selectivity in water, while also demonstrating the versatility of dynamic combinatorial chemistry for rapidly exploring the impact of subtle changes in receptor functionality on molecular recognition in water.

Environmentally responsive histidine-carboxylate zipper formation between proteins and nanoparticles

Interfacing synthetic materials with biomacromolecules provides new systems for biological applications. We report the creation of a reversible multivalent supramolecular "zipper" recognition motif between gold nanoparticles and proteins. In this assembly, carboxylate-functionalized nanoparticles interact strongly with oligohistidine tags. This interaction can be tuned through His-tag length, and offers unique binding profiles based on the pH and electrolyte concentration of the medium.

Salt-bridge energetics in halophilic proteins

Halophilic proteins have greater abundance of acidic over basic and very low bulky hydrophobic residues. Classical electrostatic stabilization was suggested as the key determinant for halophilic adaptation of protein. However, contribution of specific electrostatic interactions (i.e. salt-bridges) to overall stability of halophilic proteins is yet to be understood. To understand this, we use Adaptive-Poison-Boltzmann-Solver Methods along with our home-built automation to workout net as well as associated component energy terms such as desolvation energy, bridge energy and background energy for 275 salt-bridges from 20 extremely halophilic proteins. We then perform extensive statistical analysis on general and energetic attributes on these salt-bridges. On average, 8 salt-bridges per 150 residues protein were observed which is almost twice than earlier report. Overall contributions of salt-bridges are −3.0 kcal mol⁻¹. Majority (78%) of salt-bridges in our dataset are stable and conserved in nature. Although, average contributions of component energy terms are equal, their individual details vary greatly from one another indicating their sensitivity to local micro-environment. Notably, 35% of salt-bridges in our database are buried and stable. Greater desolvation penalty of these buried salt-bridges are counteracted by stable network salt-bridges apart from favorable equal contributions of bridge and background terms. Recruitment of extensive network salt-bridges (46%) with a net contribution of −5.0 kcal mol⁻¹ per salt-bridge, seems to be a halophilic design wherein favorable average contribution of background term (−10 kcal mol⁻¹) exceeds than that of bridge term (−7 kcal mol⁻¹). Interiors of proteins from halophiles are seen to possess relatively higher abundance of charge and polar side chains than that of mesophiles which seems to be satisfied by cooperative network salt-bridges. Overall, our theoretical analyses provide insight into halophilic signature in its specific electrostatic interactions which we hope would help in protein engineering and bioinformatics studies.

Crystal structure of the Marburg virus GP2 core domain in its postfusion conformation

Marburg virus (MARV) and Ebola virus (EBOV) are members of the family Filoviridae ("filoviruses") and cause severe hemorrhagic fever with human case fatality rates of up to 90%. Filovirus infection requires fusion of the host cell and virus membranes, a process that is mediated by the envelope glycoprotein (GP). GP contains two subunits, the surface subunit (GP1), which is responsible for cell attachment, and the transmembrane subunit (GP2), which catalyzes membrane fusion. The GP2 ectodomain contains two heptad repeat regions, N-terminal and C-terminal (NHR and CHR, respectively), that adopt a six-helix bundle during the fusion process. The refolding of this six-helix bundle provides the thermodynamic driving force to overcome barriers associated with membrane fusion. Here we report the crystal structure of the MARV GP2 core domain in its postfusion (six-helix bundle) conformation at 1.9 Å resolution. The MARV GP2 core domain backbone conformation is virtually identical to that of EBOV GP2 (reported previously), and consists of a central NHR core trimeric coiled coil packed against peripheral CHR α-helices and an intervening loop and helix-turn-helix segments. We previously reported that the stability of the MARV GP2 postfusion structure is highly pH-dependent, with increasing stability at lower pH [Harrison, J. S., Koellhoffer, J. K., Chandran, K., and Lai, J. R. (2012) Biochemistry51, 2515-2525]. We hypothesized that this pH-dependent stability provides a mechanism for conformational control such that the postfusion six-helix bundle is promoted in the environments of appropriately mature endosomes. In this report, a structural rationale for this pH-dependent stability is described and involves a high-density array of core and surface acidic side chains at the midsection of the structure, termed the "anion stripe". In addition, many surface-exposed salt bridges likely contribute to the stabilization of the postfusion structure at low pH. These results provide structural insights into the mechanism of MARV GP2-mediated membrane fusion.

Detergent-induced activation of the hepatitis C virus genotype 1b RNA polymerase

Recently, we found that sphingomyelin bound and activated hepatitis C virus (HCV) 1b RNA polymerase (RdRp), thereby recruiting the HCV replication complex into lipid raft structures. Detergents are commonly used for resolving lipids and purifying proteins, including HCV RdRp. Here, we tested the effect of detergents on HCV RdRp activity in vitro and found that non-ionic (Triton X-100, NP-40, Tween 20, Tween 80, and Brij 35) and twitterionic (CHAPS) detergents activated HCV 1b RdRps by 8-16.6 folds, but did not affect 1a or 2a RdRps. The maximum effect of these detergents was observed at around their critical micelle concentrations. On the other hand, ionic detergents (SDS and DOC) completely inactivated polymerase activity at 0.01%. In the presence of Triton X-100, HCV 1b RdRp did not form oligomers, but recruited more template RNA and increased the speed of polymerization. Comparison of polymerase and RNA-binding activity between JFH1 RdRp and Triton X-100-activated 1b RdRp indicated that monomer RdRp showed high activity because JFH1 RdRp was a monomer in physiological conditions of transcription. Besides, 502H plays a key role on oligomerization of 1b RdRp, while 2a RdRps which have the amino acid S at position 502 are monomers. This oligomer formed by 502H was disrupted both by high salt and Triton X-100. On the contrary, HCV 1b RdRp completely lost fidelity in the presence of 0.02% Triton X-100, which suggests that caution should be exercised while using Triton X-100 in anti-HCV RdRp drug screening tests.

A hamiltonian replica exchange method for building protein-protein interfaces applied to a leucine zipper

Leucine zippers consist of alpha helical monomers dimerized (or oligomerized) into alpha superhelical structures known as coiled coils. Forming the correct interface of a dimer from its monomers requires an exploration of configuration space focused on the side chains of one monomer that must interdigitate with sites on the other monomer. The aim of this work is to generate good interfaces in short simulations starting from separated monomers. Methods are developed to accomplish this goal based on an extension of a previously introduced [Su and Cukier, J. Phys. Chem. B 113, 9595, (2009)] hamiltonian temperature replica exchange method (HTREM), which scales the hamiltonian in both potential and kinetic energies that was used for the simulation of dimer melting curves. The new method, HTREM_MS (MS designates mean square), focused on interface formation, adds restraints to the hamiltonians for all but the physical system, which is characterized by the normal molecular dynamics force field at the desired temperature. The restraints in the nonphysical systems serve to prevent the monomers from separating too far, and have the dual aims of enhancing the sampling of close in configurations and breaking unwanted correlations in the restrained systems. The method is applied to a 31-residue truncation of the 33-residue leucine zipper (GCN4-p1) of the yeast transcriptional activator GCN4. The monomers are initially separated by a distance that is beyond their capture length. HTREM simulations show that the monomers oscillate between dimerlike and monomerlike configurations, but do not form a stable interface. HTREM_MS simulations result in the dimer interface being faithfully reconstructed on a 2 ns time scale. A small number of systems (one physical and two restrained with modified potentials and higher effective temperatures) are sufficient. An in silico mutant that should not dimerize because it lacks charged residues that provide electrostatic stabilization of the dimer does not with HTREM_MS, giving confidence in the method. The interface formation time scale is sufficiently short that using HTREM_MS as a screening tool to validate leucine zipper design methods may be feasible.

The effects of pK(a) tuning on the thermodynamics and kinetics of folding: design of a solvent-shielded carboxylate pair at the a-position of a coiled-coil

The tuning of the pK(a) of ionizable residues plays a critical role in various protein functions, such as ligand-binding, catalysis, and allostery. Proteins harness the free energy of folding to position ionizable groups in highly specific environments that strongly affect their pK(a) values. To investigate the interplay among protein folding kinetics, thermodynamics, and pK(a) modulation, we introduced a pair of Asp residues at neighboring interior positions of a coiled-coil. A single Asp residue was replaced for an Asn side chain at the a-position of the coiled-coil from GCN4, which was also crosslinked at the C-terminus via a flexible disulfide bond. The thermodynamic and kinetic stability of the system was measured by circular dichroism and stopped-flow fluorescence as a function of pH and concentration of guanidine HCl. Both sets of data are consistent with a two-state equilibrium between fully folded and unfolded forms. Distinct pK(a) values of 6.3 and 5.35 are assigned to the first and second protonation of the Asp pair; together they represent an energetic difference of 5 kcal/mol relative to the protonation of two Asp residues with unperturbed pK(a) values. Analysis of the rate data as a function of pH and denaturant concentration allowed calculation of the kinetic constants for the conformational transitions of the peptide with the Asp residues in the doubly protonated, singly protonated, and unprotonated forms. The doubly and singly protonated forms fold rapidly, and a ϕ-value analysis shows that their contribution to folding occurs subsequent to the transition state ensemble for folding. By contrast, the doubly charged state shows a reduced rate of folding and a ϕ-value near 0.5 indicative of a repulsive interaction, and possibly also heterogeneity in the transition state ensemble.

Two intracellular helices of G-protein coupling receptors could generally support oligomerization and coupling with transducers

For many G-protein coupling receptors (GPCRs), the upkeep of receptor dimers could depend on association with functional Gi α subunits. This is known for Y1, Y2 and Y4 neuropeptide Y receptors [presented in the companion paper (Estes et al., Amino Acids, doi: 10.1007/s00726-010-0642-z , 2010)]. Interactions with transducers use mainly intracellular domains of the receptors. Intracellular loops 1 and 2 in GPCRs are short and lack extensive helicity that could support transducer anchoring. Interaction with G-proteins is known to use the juxtamembrane Helix 8 in the fourth intracellular domain, for which we document a helix-stabilizing n/(n + 4) pattern of large hydrophobic sidechains. Another intracellular helix located in the C-terminal portion of the third intracellular loop does not display a strong stabilizing pattern, and is found in many studies to serve dynamically in association and activation of transducers and effectors. We show that these tracts share features across metazoan phyla not only in opsins and opsin-like receptors (including the Y receptors), but also in Taste-2 and Frizzled receptors. Similarities of these helices across GPCR groups could have both phylogenetic and functional roots.

Titration_DB: storage and analysis of NMR-monitored protein pH titration curves

NMR-monitored pH titration experiments are routinely used to measure site-specific protein pKa values. Accurate experimental pKa values are essential in dissecting enzyme catalysis, in studying the pH-dependence of protein stability and ligand binding, in benchmarking pKa prediction algorithms, and ultimately in understanding electrostatic effects in proteins. However, due to the complex ways in which pH-dependent electrostatic and structural changes manifest themselves in NMR spectra, reported apparent pKa values are often dependent on the way that NMR pH-titration curves are analyzed. It is therefore important to retain the raw NMR spectroscopic data to allow for documentation and possible re-interpretation. We have constructed a database of primary NMR pH-titration data, which is accessible via a web interface. Here, we report statistics of the database contents and analyze the data with a global perspective to provide guidelines on best practice for fitting NMR titration curves. Titration_DB is available at http://enzyme.ucd.ie/Titration_DB. Proteins 2010. (c) 2009 Wiley-Liss, Inc.

Crystal structure of a trimeric form of the K(V)7.1 (KCNQ1) A-domain tail coiled-coil reveals structural plasticity and context dependent changes in a putative coiled-coil trimerization motif

Coiled-coils are widespread protein-protein interaction motifs typified by the heptad repeat (abcdefg)(n) in which "a" and "d" positions are hydrophobic residues. Although identification of likely coiled-coil sequences is robust, prediction of strand order remains elusive. We present the X-ray crystal structure of a short form (residues 583-611), "Q1-short," of the coiled-coil assembly specificity domain from the voltage-gated potassium channel Kv7.1 (KCNQ1) determined at 1.7 A resolution. Q1-short lacks one and half heptads present in a previously studied tetrameric coiled-coil construct, Kv7.1 585-621, "Q1-long." Surprisingly, Q1-short crystallizes as a trimer. In solution, Q1-short self-assembles more poorly than Q1-long and depends on an R-h-x-x-h-E motif common to trimeric coiled-coils. Addition of native sequences that include "a" and "d" positions C-terminal to Q1-short overrides the R-h-x-x-h-E motif influence and changes assembly state from a weakly associated trimer to a strongly associated tetramer. These data provide a striking example of a naturally occurring amino sequence that exhibits context-dependent folding into different oligomerization states, a three-stranded versus a four-stranded coiled-coil. The results emphasize the degenerate nature of coiled-coil energy landscapes in which small changes can have drastic effects on oligomerization. Discovery of these properties in an ion channel assembly domain and prevalence of the R-h-x-x-h-E motif in coiled-coil assembly domains of a number of different channels that are thought to function as tetrameric assemblies raises the possibility that such sequence features may be important for facilitating the assembly of intermediates en route to the final native state.

Electrostatic contributions to the stabilities of native proteins and amyloid complexes

The ability to predict electrostatic contributions to protein stability from structure has been a long-standing goal of experimentalists and theorists. With recent advances in NMR spectroscopy, it is possible to determine pK(a) values of all ionizable residues for at least small proteins, and to use the pK(a) shift between the folded and unfolded states to calculate the thermodynamic contribution from a change in charge to the change in free energy of unfolding. Results for globular proteins and for α-helical coiled coils show that electrostatic contributions to stability are typically small on an individual basis, particularly for surface-exposed residues. We discuss why NMR often suggests smaller electrostatic contributions to stability than X-ray crystallography or site-directed mutagenesis, and discuss the type of information needed to improve structure-based modeling of electrostatic forces. Large pK(a) shifts from random coil values are observed for proteins bound to negatively charged sodium dodecyl sulfate micelles. The results suggest that electrostatic interactions between proteins and charges on the surfaces of membrane lipid bilayers could be a major driving force in stabilizing the structures of peripheral membrane proteins. Finally, we discuss how changes in ionization states affect amyloid-β fibril formation and suggest that electrostatic repulsion may be a common destabilizing force in amyloid fibrils.

Role of conserved salt bridges in homeodomain stability and DNA binding

The sequence information available for homeodomains reveals that salt bridges connecting pairs 19/30, 31/42, and 17/52 are frequent, whereas aliphatic residues at these sites are rare and mainly restricted to proteins from homeotherms. We have analyzed the influence of salt and hydrophobic bridges at these sites on the stability and DNA binding properties of human Hesx-1 homeodomain. Regarding the protein stability, our analysis shows that hydrophobic side chains are clearly preferred at positions 19/30 and 31/42. This stabilizing influence results from the more favorable packing of the aliphatic side chains with the protein core, as illustrated by the three-dimensional solution structure of a thermostable variant, herein reported. In contrast only polar side chains seem to be tolerated at positions 17/52. Interestingly, despite the significant influence of pairs 19/30 and 31/42 on the stability of the homeodomain, their effect on DNA binding ranges from modest to negligible. The observed lack of correlation between binding strength and conformational stability in the analyzed variants suggests that salt/hydrophobic bridges at these specific positions might have been employed by evolution to independently modulate both properties.

Salt bridges in the hyperthermophilic protein Ssh10b are resilient to temperature increases

A double mutant cycle (DMC) approach was employed to estimate the effect of temperature on the contribution of two highly conserved salt bridges to protein stability in the hyperthermophilic protein Ssh10b. The coupling free energy were 2.4 +/- 0.4 kJ/mol at 298 K and 2.2 +/- 0.4 kJ/mol at 353 K for Glu-54/Arg-57, and 6.0 +/- 0.2 kJ/mol at 298 K and 5.9 +/- 0.6 kJ/mol at 353 K for Glu-36/Lys-68. The stability free energy of Ssh10b decrease greatly with increasing temperature, while the direct contribution of these two salt bridges to protein stability remain almost constant, providing evidence supporting the theoretical prediction that salt bridges are extremely resilient to temperature increases and thus are specially suited to improving protein stability at high temperatures. The reason for the difference in coupling free energy between salt bridges Glu-54/Arg-57 and Glu-36/Lys-68 is discussed. Comparing our results with published DMC data for the contribution of salt bridges to stability in other proteins, we found that the energy contribution of a salt bridge formed by two charged residues far apart in the primary sequence is higher than that of those formed between two very close ones. Implications of this finding are useful for engineering proteins with enhanced thermostability.

The effects of NaCl concentration and pH on the stability of hyperthermophilic protein Ssh10b

Background

Hyperthermophiles constitute a group of microorganisms with an optimum growth temperature of between 80°C and 100°C. Although the molecular underpinnings of protein thermostabilization have been the focus of many theoretical and experimental efforts, the properties leading to the higher denaturation temperature of hyperthermophilic proteins are still controversial. Among the large number of factors identified as responsible for the thermostability of hyperthermophilic proteins, the electrostatic interactions are thought to be a universally important factor.

Results

In this study, we report the effects of pH and salt concentration on the urea-induced denaturation of the protein Ssh10b from a hyperthermophile in low ionic strength buffer. In the absence of NaCl, the unfolding ΔG of the protein increased from about 33 kJ/mol at pH 3 to about 78 kJ/mol at pH 10. At all values of pH, the ΔG increased with increasing NaCl concentration, indicating that salt stabilizes the protein significantly.

Conclusion

These findings suggests that the increased number of charged residues and ion pairs in the protein Ssh10b from hyperthermophiles does not contribute to the stabilization of the folded protein, but may play a role in determining the denatured state ensemble and also in increasing the denaturation temperature.

Electrostatic contributions to the stability of the GCN4 leucine zipper structure

Ion pairs are ubiquitous in X-ray structures of coiled coils, and mutagenesis of charged residues can result in large stability losses. By contrast, pK(a) values determined by NMR in solution often predict only small contributions to stability from charge interactions. To help reconcile these results we used triple-resonance NMR to determine pK(a) values for all groups that ionize between pH 1 and 13 in the 33 residue leucine zipper fragment, GCN4p. In addition to the native state we also determined comprehensive pK(a) values for two models of the GCN4p denatured state: the protein in 6 M urea, and unfolded peptide fragments of the protein in water. Only residues that form ion pairs in multiple X-ray structures of GCN4p gave large pK(a) differences between the native and denatured states. Moreover, electrostatic contributions to stability were not equivalent for oppositely charged partners in ion pairs, suggesting that the interactions between a charge and its environment are as important as those within the ion pair. The pH dependence of protein stability calculated from NMR-derived pK(a) values agreed with the stability profile measured from equilibrium urea-unfolding experiments as a function of pH. The stability profile was also reproduced with structure-based continuum electrostatic calculations, although contributions to stability were overestimated at the extremes of pH. We consider potential sources of errors in the calculations, and how pK(a) predictions could be improved. Our results show that although hydrophobic packing and hydrogen bonding have dominant roles, electrostatic interactions also make significant contributions to the stability of the coiled coil.

Energetic determinants of oligomeric state specificity in coiled coils

The coiled coil is one of the simplest and best-studied protein structural motifs, consisting of two to five helices wound around each other. Empirical rules have been established on the tendency of different core sequences to form a certain oligomeric state but the physical forces behind this specificity are unclear. In this work, we model four sequences onto the structures of dimeric, trimeric, tetrameric, and pentameric coiled coils. We first examine the ability of an effective energy function (EEF1.1) to discriminate the correct oligomeric state for a given sequence. We find that inclusion of the translational, rotational, and side-chain conformational entropy is necessary for discriminating the native structures from their misassembled counterparts. The decomposition of the effective energy into residue contributions yields theoretical values for the oligomeric propensity of different residue types at different heptad positions. We find that certain calculated residue propensities are general and consistent with existing rules. For example, leucine at d favors dimers, leucine at a favors tetramers or pentamers, and isoleucine at a favors trimers. Other residue propensities are sequence context dependent. For example, glutamine at d favors trimers in one context and pentamers in another. Charged residues at e and g positions usually destabilize higher oligomers due to higher desolvation. Nonpolar residues at these positions confer pentamer specificity when combined with certain residues at positions a and d. Specifically, the pair Leua-Alag' or the inverse was found to stabilize the pentamer. The small energy gap between the native and misfolded counterparts explains why a few mutations at the core sites are sufficient to induce a change in the oligomeric state of these peptides. A large number of possible experiments are suggested by these results.

Crystal structures of tropomyosin: flexible coiled-coil

Tropomyosin (Tm) is a 400 angstroms long coiled coil protein, and with troponin it regulates contraction in skeletal and cardiac muscles in a [Ca2+]-dependent manner. Tm consists of multiple domains with diverse stabilities in the coiled coil form, thus providing Tm with dynamic flexibility. This flexibility must play important roles in the actin binding and the cooperative transition between the calcium regulated states of the entire muscle thin filament. In order to understand the flexibility of Tm in its entirety, the atomic coordinates of Tm are needed. Here we report the two crystal structures of Tm segments. One is rabbit skeletal muscle alpha-Tm encompassing residues 176-284 with an N-terminal extension of 25 residues from the leucine zipper sequence of GCN4, which includes the region that interacts with the troponin core domain. The other is alpha-Tm encompassing residues 176-273 with N- and C-terminal extensions of the leucine zipper sequences. These two crystal structures imply that this molecule is a flexible coiled coil. First, Tm's are not homogeneous and smooth coiled coils, but instead they undulate, with highly fluctuating local parameters specifying the coiled coil. Independent fluctuating showed by two crystal structures is important. Second, in the first crystal, the coiled coil is bent by 9 degrees in the region centered about Y214-E218-Y221, where the inter-helical distance has its maximum. On the other hand, no bend is observed at the same region in the second crystal even if its inter-helical distance has also its maximum. E218, an unusual negatively charged residue at the a position in the heptad repeat, seems to play the key role in destabilizing the coiled coil with alanine destabilizing clusters.

Antiparallel four-stranded coiled coil specified by a 3-3-1 hydrophobic heptad repeat

Coiled-coil sequences in proteins commonly share a seven-amino acid repeat with nonpolar side chains at the first (a) and fourth (d) positions. We investigate here the role of a 3-3-1 hydrophobic repeat containing nonpolar amino acids at the a, d, and g positions in determining the structures of coiled coils using mutants of the GCN4 leucine zipper dimerization domain. When three charged residues at the g positions in the parental sequence are replaced by nonpolar alanine or valine side chains, stable four-helix structures result. The X-ray crystal structures of the tetramers reveal antiparallel, four-stranded coiled coils in which the a, d, and g side chains interlock in a combination of knobs-into-knobs and knobs-into-holes packing. Interfacial interactions in a coiled coil can therefore be prescribed by hydrophobic-polar patterns beyond the canonical 3-4 heptad repeat. The results suggest that the conserved, charged residues at the g positions in the GCN4 leucine zipper can impart a negative design element to disfavor thermodynamically more stable, antiparallel tetramers.

Mutating His29, His125, His133 or His158 abolishes glycosylphosphatidylinositol-specific phospholipase D catalytic activity

Glycosylphosphatidylinositol (GPI)-specific phospholipase D (GPI-PLD) specifically cleaves GPIs. This phospholipase D is a secreted protein consisting of two domains: an N-terminal catalytic domain and a predicted C-terminal b-propeller. Although the biochemical properties of GPI-PLD have been extensively studied, its catalytic site has not been identified. We hypothesized that a histidine residue(s) may play a critical role in the catalytic activity of GPI-PLD, based on the observations that (i) Zn2+, which utilizes histidine residues for binding, is required for GPI-PLD catalytic activity, (ii) a phosphohistidine intermediate is involved in phospholipase D hydrolysis of phosphatidylcholine, (iii) computer modelling suggests a catalytic site containing histidine residues, and (iv) our observation that diethyl pyrocarbonate, which modifies histidine residues, inhibits GPI-PLD catalytic activity. Individual mutation of the ten histidine residues to asparagine in the catalytic domain of murine GPI-PLD resulted in three general phenotypes: not secreted or retained (His56 or His88), secreted with catalytic activity (His34, His81, His98 or His219) and secreted without catalytic activity (His29, His125, His133 or His158). Changing His133 but not His29, His125 or His158 to Cys resulted in a mutant that retained catalytic activity, suggesting that at least His133 is involved in Zn2+ binding. His133 and His158 also retained the biochemical properties of wild-type GPI-PLD including trypsin cleavage pattern and phosphorylation by protein kinase A. Hence, His29, His125, His133 and His158 are required for GPI-PLD catalytic activity.

Mechanism and Thermodynamics of Guanidinium Chloride-induced Denaturation of ALS-associated Mutant Cu,Zn Superoxide Dismutases

Mutations in human copper zinc superoxide dismutase (hSOD) that are associated with amyotrophic lateral sclerosis (ALS) have been proposed to destabilize the protein and thereby enhance toxic protein aggregation. In previous studies, denaturation of metallated (holo) hSODs was found to be irreversible, and complicated by the formation of intermolecular disulfide bonds. Here, ALS-associated mutations (E100G, G93A, G85R and A4V) are introduced into a pseudo wild-type background containing no free cysteine residues. The guanidinium chloride-induced denaturation of the holo proteins is generally found to be highly reversible (except for A4V, which tended to aggregate), enabling quantitative analysis of the effects of the mutations on protein stability. Denaturation and renaturation curves were monitored by tryptophan fluorescence, circular dichroism, enzyme activity, chemical cross-linking and analytical sedimentation, as a function of equilibration time and protein concentration. There is strong kinetic hysteresis, with curves requiring exceptionally long times (many days for pseudo wild-type) to reach equilibrium, and evidence for the formation of kinetic and equilibrium intermediate(s), which are more highly populated at lower protein concentrations. The effects of metal dissociation were included in the data fitting. The full protein concentration dependence is best described using a three-state model involving metallated native dimer, metallated monomeric intermediate and unfolded monomers with no bound metals; however, at high protein concentrations the unfolding approaches a two-state transition with metal binding to both the native dimers and unfolded monomers. We show that the E100G, G93A and G85R mutations decrease overall protein stability, largely by decreasing monomer stability with little effect on dimer dissociation. Comparison of the chemical denaturation data with ALS disease characteristics suggests that aggregation of some mutant hSOD may occur through increased population of partially folded states that are less stable than the monomeric intermediate and accessed from the destabilized holo protein.

Anticooperativity in a Glu-Lys-Glu salt bridge triplet in an isolated alpha-helical peptide

Salt bridges between oppositely charged side chains are well-known to stabilize protein structure, though their contributions vary considerably. Here we study Glu-Lys and Lys-Glu salt bridges, formed when the residues are spaced i, i + 4 surface of an isolated alpha-helix in aqueous solution. Both are stabilizing by -0.60 and -1.02 kcal/mol, respectively, when the interacting residues are fully charged. When the side chains are spaced i, i + 4, i + 8, forming a Glu-Lys-Glu triplet, the second salt bridge provides no additional stabilization to the helix. We attribute this to the inability of the central Lys to form two salt bridges simultaneously. Analysis of these salt bridges in protein structures shows that the Lys-Glu interaction is dominant, with the side chains of the Glu-Lys pair far apart.

Statistical characterization of salt bridges in proteins

The structure and folding mechanism of a given protein are determined by many factors, including the electrostatic interactions between charged residues of protein molecules known in general as salt bridges. In this study, analyses were conducted on 10,370 salt bridges in 2017 proteins and the results compared to previous statistical surveys of 36 protein structures. Although many of the general trends remained consistent with other studies, more detailed information was illuminated by the larger dataset. In particular, it was shown that there is a strong correlation between secondary structure and salt bridge formation, and that salt bridges display preferential formation in an environment of about 30% solvent accessible surface area.

Electrostatic interactions involving lysine make major contributions to collagen triple-helix stability

Important stabilizing features for the collagen triple helix include the presence of Gly as every third residue, a high content of imino acids, and interchain hydrogen bonds. Host-guest peptides have been used previously to characterize triple-helix propensities of individual residues and Gly-X-Y triplets. Here, comparison of the thermal stabilities of host-guest peptides of the form (Gly-Pro-Hyp)3-Gly-X-Y-Gly-X'-Y'-(Gly-Pro-Hyp)3 extends the study to adjacent tripeptide sequences, to encompass the major classes of potential direct intramolecular interactions. Favorable hydrophobic interactions were observed, as well as stabilizing intrachain interactions between residues of opposite charge in the i and i + 3 positions. However, the greatest gain in triple-helix stability was achieved in the presence of Gly-Pro-Lys-Gly-Asp/Glu-Hyp sequences, leading to a T(m) value equal to that seen for a Gly-Pro-Hyp-Gly-Pro-Hyp sequence. This stabilization is seen for Lys but not for Arg and can be assigned to interchain ion pairs, as shown by molecular modeling. Computational analysis shows that Lys-Gly-Asp/Glu sequences are present at a frequency much greater than expected in collagen, suggesting this interaction is biologically important. These results add significantly to the understanding of which surface ion pairs can contribute to protein stability.

A Phosphoserine−Lysine Salt Bridge within an α-Helical Peptide, the Strongest α-Helix Side-Chain Interaction Measured to Date†

Phosphorylation is ubiquitous in control of protein activity, yet its effects on protein structure are poorly understood. Here we investigate the effect of serine phosphorylation in the interior of an alpha-helix when a salt bridge is present between the phosphate group and a positively charged side chain (in this case lysine) at i,i + 4 spacing. The stabilization of the helix is considerable and can overcome the intrinsically low preference of phosphoserine for the interior of the helix. The effect is pH dependent, as both the lysine and phosphate groups are titratable, and so calculations are given for several charge combinations. These results, with our previous work, highlight the different, context-dependent effects of phosphorylation in the alpha-helix. The interaction between the phosphate(2)(-) group and the lysine side chain is the strongest yet recorded in helix-coil studies. The results are of interest both in de novo design of peptides and in understanding the structural modes of control by phosphorylation.

Electrostatic contributions to the kinetics and thermodynamics of protein assembly

The role of electrostatic interactions in the assembly of a native protein structure was studied using fragment complementation. Contributions of salt, pH, or surface charges to the kinetics and equilibrium of calbindin D(9k) reconstitution was measured in the presence of Ca(2+) using surface plasmon resonance and isothermal titration calorimetry. Whereas surface charge substitutions primarily affect the dissociation rate constant, the association rates are correlated with subdomain net charge in a way expected for Coulomb interactions. The affinity is reduced in all mutants, with the largest effect (260-fold) observed for the double mutant K25E+K29E. At low net charge, detailed charge distribution is important, and charges remote from the partner EF-hand have less influence than close ones. The effects of salt and pH on the reconstitution are smaller than mutational effects. The interaction between the wild-type EF-hands occurs with high affinity (K(A) = 1.3 x 10(10) M(-1); K(D) = 80 pM). The enthalpy of association is overall favorable and there appears to be a very large favorable entropic contribution from the desolvation of hydrophobic surfaces that become buried in the complex. Electrostatic interactions contribute significantly to the affinity between the subdomains, but other factors, such as hydrophobic interactions, dominate.

Protein structure, stability and solubility in water and other solvents

Proteins carry out the most difficult tasks in living cells. They do so by interacting specifically with other molecules. This requires that they fold to a unique, globular conformation that is only marginally more stable than the large ensemble of unfolded states. The folded state is stabilized mainly by the burial and tight packing of over 80% of the peptide groups and non-polar side chains. If life as we know it is to exist in a solvent other than water, the folded state must be stable and soluble in the new solvent. Our analysis suggests that proteins will be unstable in most polar solvents such as ethanol, extremely stable in non-polar solvents such as cyclohexane, and even more stable in a vacuum. Our solubility studies suggest that protein solubility will be markedly lower in polar solvents such as ethanol and that proteins will be essentially insoluble in non-polar solvents such as cyclohexane. For these and other reasons it seems unlikely that the life we know could exist in any solvent system other than water.

The hyperthermophile protein Sso10a is a dimer of winged helix DNA-binding domains linked by an antiparallel coiled coil rod

Sso10a is a member of a group of DNA-binding proteins thought to be important in chromatin structure and regulation in the hyperthermophilic archaeon Sulfolobus solfataricus. We have determined the structure of Sso10a to 1.47A resolution directly with unlabelled native crystals by a novel approach using sulfur single-wavelength anomalous scattering (SAS) from a chromium X-ray source. The 95 amino acid residue protein contains a winged helix DNA-binding domain with an extended C-terminal alpha-helix that leads to dimerization by forming a two-stranded, antiparallel coiled-coil rod. The winged helix domains are at opposite ends of the extended coiled coil with two putative DNA-recognition helices separated by 55A and rotated by 83 degrees. Formation of stable dimers in solution is demonstrated by both analytical ultracentrifugation and differential scanning calorimetry. With a T0 of 109 degrees C, Sso10a is one of the most stable two-stranded coiled coils known. The coiled coil contains a rare aspartate residue (D69) in the normally hydrophobic d position of the heptad repeat, with two aspartate-lysine (d-g') interhelical ion pairs in the symmetrical dimer. Mutation of D69 to alanine resulted in an increase in thermal stability, indicating that destabilization resulting from the partially buried aspartate residue cannot be offset by ion pair formation. Possible DNA-binding interactions are discussed on the basis of comparisons to other winged helix proteins. The structure of Sso10a provides insight into the structures of the conserved domain represented by COG3432, a group of more than 20 hypothetical transcriptional regulators coded in the genomic sequences of both crenarchaeota and euryarchaeota.

Inverse Electrostatic Effect: Electrostatic Repulsion in the Unfolded State Stabilizes a Leucine Zipper,

The pH-dependent stability of a protein is strongly affected by electrostatic interactions between ionizable residues in the folded as well as unfolded state. Here we characterize the individual contributions of charged Glu and His residues to stability and determine the NMR structure of the designed, heterodimeric leucine zipper AB consisting of an acidic A chain and a basic B chain. Thermodynamic parameters are compared with those of the homologous leucine zipper AB(SS) in which the A and B chains are disulfide-linked. NMR structures of AB based on (1)H NMR data collected at 600 MHz converge, and formation of the same six interchain salt bridges found previously in disulfide-linked AB(SS) [Marti, D. N., and Bosshard, H. R. (2003) J. Mol. Biol. 330, 621-637] is indicated. While the structures of AB and AB(SS) are very similar, their pH-dependent relative stabilities are strikingly different. The stability of AB peaks at pH approximately 4.5 and is higher at pH 8 than at pH 2. In contrast, AB(SS) is most stable at acidic pH where no interhelical salt bridges are formed. The different energetic contributions of charged Glu and His residues to stability of the two coiled coil structures were evaluated from pK(a) shifts induced by folding. The six charged Glu residues involved in salt bridges stabilize leucine zipper AB by 4.5 kJ/mol yet destabilize disulfide-linked AB(SS) by -1.1 kJ/mol. Two non-ion-paired Glu charges destabilize AB by only -1.8 kJ/mol but AB(SS) by -5.6 kJ/mol. The higher relative stability of AB at neutral pH is not caused by more favorable electrostatic interactions in the folded leucine zipper. It is due mainly to unfavorable electrostatic interactions in the unfolded A and B chains and may therefore be called an inverse electrostatic effect. This study illustrates the importance of residual interactions in the unfolded state and how the energetics of the unfolded state affect the stability of the folded protein.

Salt-bridges can stabilize but do not accelerate the folding of the homodimeric coiled-coil peptide GCN4-p1

Double mutant cycle analysis was employed to ascertain the role of intra- and interchain salt-bridges in the folding and stability of the dimeric coiled-coil peptide, GCN4-p1, the 33-residue leucine zipper domain of the transcriptional activator GCN4. Equilibrium circular dichroism studies of the urea-induced unfolding reaction at neutral pH revealed that both types of ionic interactions, localized primarily in the N-terminal portion of the molecule, enhance the stability of the native coiled-coil. By contrast, comparable stopped-flow circular dichroism studies indicate that the salt-bridge interactions, with one possible exception, are not well formed in the transition state for folding. Although the E22Q/R25A double mutant failed to fold, fragmentation studies suggest that the E22/R25 intramolecular salt-bridge may play a critical role in stabilizing C-terminal nascent helices that drive the association reaction. The remaining salt-bridges appear to stabilize the parallel-stranded coiled-coil architecture of GCN4-p1 only after the peptide traverses the rate-limiting, dimeric transition state.