skip to main content

Explore Research Products in the NSF-PAR

Title: How a highly acidic SH3 domain folds in the absence of its charged peptide target
Abstract

Charged residues on the surface of proteins are critical for both protein stability and interactions. However, many proteins contain binding regions with a high net charge that may destabilize the protein but are useful for binding to oppositely charged targets. We hypothesized that these domains would be marginally stable, as electrostatic repulsion would compete with favorable hydrophobic collapse during folding. Furthermore, by increasing the salt concentration, we predict that these protein folds would be stabilized by mimicking some of the favorable electrostatic interactions that take place during target binding. We varied the salt and urea concentrations to probe the contributions of electrostatic and hydrophobic interactions for the folding of the yeast SH3 domain found in Abp1p. The SH3 domain was significantly stabilized with increased salt concentrations due to Debye–Huckel screening and a nonspecific territorial ion‐binding effect. Molecular dynamics and NMR show that sodium ions interact with all 15 acidic residues but do little to change backbone dynamics or overall structure. Folding kinetics experiments show that the addition of urea or salt primarily affects the folding rate, indicating that almost all the hydrophobic collapse and electrostatic repulsion occur in the transition state. After the transition state formation, modest yet favorable short‐range salt bridges are formed along with hydrogen bonds, as the native state fully folds. Thus, hydrophobic collapse offsets electrostatic repulsion to ensure this highly charged binding domain can still fold and be ready to bind to its charged peptide targets, a property that is likely evolutionarily conserved over 1 billion years.

  more » « less
Award ID(s):
1852677 2018427
NSF-PAR ID:
10409900
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
Protein Science
Volume:
32
Issue:
5
ISSN:
0961-8368
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; ( , Physical Chemistry Chemical Physics)
    null (Ed.)
    Electrostatic protein/DNA interactions arise from the neutralization of the DNA phosphodiester backbone as well as coupled exchanges by charged protein residues as salt bridges or with mobile ions. Much focus has been and continues to be paid to interfacial ion pairs with DNA. The role of extra-interfacial ionic interactions, particularly as dynamic drivers of DNA sequence selectivity, remain poorly known. The ETS family of transcription factors represents an attractive model for addressing this knowledge gap given their diverse ionic composition in primary structures that fold to a tightly conserved DNA-binding motif. To probe the importance of extra-interfacial salt bridges in DNA recognition, we compared the salt-dependent binding by Elk1 with ETV6, two ETS homologs differing markedly in ionic composition. While both proteins exhibit salt-dependent binding with cognate DNA that corresponds to interfacial phosphate contacts, their nonspecific binding diverges from cognate binding as well as each other. Molecular dynamics simulations in explicit solvent, which generated ionic interactions in agreement with the experimental binding data, revealed distinct salt-bridge dynamics in the nonspecific complexes formed by the two proteins. Impaired DNA contact by ETV6 resulted in fewer backbone contacts in the nonspecific complex, while Elk1 exhibited a redistribution of extra-interfacial salt bridges via residues that are non-conserved between the two ETS relatives. Thus, primary structure variation in ionic residues can encode highly differentiated specificity mechanisms in a highly conserved DNA-binding motif. 
    more » « less
  2. ; ; ; ; ; ; ; ( , Frontiers in Molecular Biosciences)
    Disordered proline-rich motifs are common across the proteomes of many species and are often involved in protein-protein interactions. Proline is a unique amino acid due to the covalent bond between the backbone nitrogen and the proline side chain. The resulting five-membered ring allows proline to sample the cis state about its peptide bond, which other residues cannot do as readily. Because proline-rich disordered sequences exist as ensembles that likely include structures with the proline peptide bond in cis , a robust methodology to accurately account for these conformations in the overall ensemble is crucial. Observing the cis conformations of proline in a disordered sequence is challenging both experimentally and computationally. Nitrogen-hydrogen NMR spectroscopy cannot directly observe proline residues, which lack an amide bond, and computational methods struggle to overcome the large kinetic barrier between the cis and trans states, since isomerization usually occurs on the order of seconds. In the current work, Gaussian accelerated molecular dynamics was used to overcome this free energy barrier and simulate proline isomerization in a tetrapeptide (KPTP) and in the 12-residue proline-rich SH3 binding peptide, ArkA. We found that Gaussian accelerated molecular dynamics, when combined with a lowered peptide bond dihedral angle potential energy barrier (15 kcal/mol), allowed sufficient sampling of the proline cis and trans states on a microsecond timescale. All ArkA prolines spend a significant fraction of time in cis , leading to a more compact ensemble with less polyproline II helix structure than an ArkA ensemble with all peptide bonds in trans . The ensemble containing cis prolines also matches more closely to in vitro circular dichroism data than the all- trans ensemble. The ability of the ArkA prolines to isomerize likely affects the peptide’s ability to bind its partner SH3 domain, and should be studied further. This is the first molecular dynamics simulation study of proline isomerization in a biologically relevant proline-rich sequence that we know of, and a similar protocol could be applied to study multi-proline isomerization in other proline-containing proteins to improve conformational diversity and agreement with in vitro data. 
    more » « less
  3. ; ; ; ; ( , Protein Science)
    Abstract

    SurA is thought to be the most important periplasmic chaperone for outer membrane protein (OMP) biogenesis. Its structure is composed of a core region and two peptidylprolyl isomerase domains, termed P1 and P2, connected by flexible linkers. As such these three independent folding units are able to adopt a number of distinct spatial positions with respect to each other. The conformational dynamics of these domains are thought to be functionally important yet are largely unresolved. Here we address this question of the conformational ensemble using sedimentation equilibrium, small‐angle neutron scattering, and folding titrations. This combination of orthogonal methods converges on a SurA population that is monomeric at physiological concentrations. The conformation that dominates this population has the P1 and core domains docked to one another, for example, “P1‐closed” and the P2 domain extended in solution. We discovered that the distribution of domain orientations is defined by modest and favorable interactions between the core domain and either the P1 or the P2 domains. These two peptidylprolyl domains compete with each other for core‐binding but are thermodynamically uncoupled. This arrangement implies two novel insights. Firstly, an open conformation must exist to facilitate P1 and P2 exchange on the core, indicating that the open client‐binding conformation is populated at low levels even in the absence of client unfolded OMPs. Secondly, competition between P1 and P2 binding paradoxically occludes the client binding site on the core, which may serve to preserve the reservoir of binding‐competent apo‐SurA in the periplasm.

     
    more » « less
  4. ; ; ; ; ; ; ( , Proceedings of the National Academy of Sciences)

    Much attention is being paid to conformational biases in the ensembles of intrinsically disordered proteins. However, it is currently unknown whether or how conformational biases within the disordered ensembles of foldable proteins affect function in vivo. Recently, we demonstrated that water can be a good solvent for unfolded polypeptide chains, even those with a hydrophobic and charged sequence composition typical of folded proteins. These results run counter to the generally accepted model that protein folding begins with hydrophobicity-driven chain collapse. Here we investigate what other features, beyond amino acid composition, govern chain collapse. We found that local clustering of hydrophobic and/or charged residues leads to significant collapse of the unfolded ensemble of pertactin, a secreted autotransporter virulence protein fromBordetella pertussis, as measured by small angle X-ray scattering (SAXS). Sequence patterns that lead to collapse also correlate with increased intermolecular polypeptide chain association and aggregation. Crucially, sequence patterns that support an expanded conformational ensemble enhance pertactin secretion to the bacterial cell surface. Similar sequence pattern features are enriched across the large and diverse family of autotransporter virulence proteins, suggesting sequence patterns that favor an expanded conformational ensemble are under selection for efficient autotransporter protein secretion, a necessary prerequisite for virulence. More broadly, we found that sequence patterns that lead to more expanded conformational ensembles are enriched across water-soluble proteins in general, suggesting protein sequences are under selection to regulate collapse and minimize protein aggregation, in addition to their roles in stabilizing folded protein structures.

     
    more » « less
  5. ; ; ; ( , Protein Science)
    Abstract

    Amide−π interactions, in which an amide interacts with an aromatic group, are ubiquitous in biology, yet remain understudied relative to other noncovalent interactions. Recently, we demonstrated that an electrostatically tunable amide−π interaction is key to recognition of histone acyllysine by the AF9 YEATS domain, a reader protein which has emerged as a therapeutic target due to its dysregulation in cancer. Amide isosteres are commonly employed in drug discovery, often to prevent degradation by proteases, and have proven valuable in achieving selectivity when targeting epigenetic proteins. However, like amide−π interactions, interactions of amide isosteres with aromatic rings have not been thoroughly studied despite widespread use. Herein, we evaluate the recognition of a series of amide isosteres by the AF9 YEATS domain using genetic code expansion to evaluate the amide isostere−π interaction. We show that compared to the amide−π interaction with the native ligand, each isostere exhibits similar electrostatic tunability with an aromatic residue in the binding pocket, demonstrating that the isosteres maintain similar interactions with the aromatic residue. We identify a urea‐containing ligand that binds with enhanced affinity for the AF9 YEATS domain, offering a promising starting point for inhibitor development. Furthermore, we demonstrate that carbamate and urea isosteres of crotonyllysine are resistant to enzymatic removal by SIRT1, a protein that cleaves acyl post‐translational modifications, further indicating the potential of amide isosteres in YEATS domain inhibitor development. These results also provide experimental precedent for interactions of these common drug discovery moieties with aromatic rings that can inform computational methods.

     
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email