Role of Electrostatic Interactions in Binding of Peptides and Intrinsically Disordered Proteins to Their Folded Targets: 2. The Model of Encounter Complex Involving the Double Mutant of the c-Crk N-SH3 Domain and Peptide Sos

Biochemistry ◽  
2016 ◽  
Vol 55 (12) ◽  
pp. 1784-1800 ◽  
Author(s):  
Tairan Yuwen ◽  
Yi Xue ◽  
Nikolai R. Skrynnikov
2020 ◽  
Author(s):  
Gabriella J. Gerlach ◽  
Rachel Carrock ◽  
Robyn Stix ◽  
Elliott J. Stollar ◽  
K. Aurelia Ball

AbstractProtein-protein interactions are involved in a wide range of cellular processes. These interactions often involve intrinsically disordered proteins (IDPs) and protein binding domains. However, the details of IDP binding pathways are hard to characterize using experimental approaches, which can rarely capture intermediate states present at low populations. SH3 domains are common protein interaction domains that typically bind proline-rich disordered segments and are involved in cell signaling, regulation, and assembly. We hypothesized, given the flexibility of SH3 binding peptides, that their binding pathways include multiple steps important for function. Molecular dynamics simulations were used to characterize the steps of binding between the yeast Abp1p SH3 domain (AbpSH3) and a proline-rich IDP, ArkA. Before binding, the N-terminal segment 1 of ArkA is pre-structured and adopts a polyproline II helix, while segment 2 of ArkA (C-terminal) adopts a 310 helix, but is far less structured than segment 1. As segment 2 interacts with AbpSH3, it becomes more structured, but retains flexibility even in the fully engaged state. Binding simulations reveal that ArkA enters a flexible encounter complex before forming the fully engaged bound complex. In the encounter complex, transient nonspecific hydrophobic and long-range electrostatic contacts form between ArkA and the binding surface of SH3. The encounter complex ensemble includes conformations with segment 1 in both the forward and reverse orientation, suggesting that segment 2 may play a role in stabilizing the correct binding orientation. While the encounter complex forms quickly, the slow step of binding is the transition from the disordered encounter ensemble to the fully engaged state. In this transition, ArkA makes specific contacts with AbpSH3 and buries more hydrophobic surface. Simulating the binding between ApbSH3 and ArkA provides insight into the role of encounter complex intermediates and nonnative hydrophobic interactions for other SH3 domains and IDPs in general.Author SummaryComplex cellular processes are mediated by interactions between proteins, and to determine how these interactions affect cellular function and binding kinetics we often must understand the protein binding pathway. Many protein interaction domains, such as the SH3 domain, bind to intrinsically disordered proteins in a coupled folding and binding process. Using molecular dynamics simulations, we find that the binding of the disordered ArkA peptide to the yeast Abp1p SH3 domain proceeds through a flexible, disordered encounter complex before reaching a stable fully bound state. The encounter complex is stabilized by nonspecific long-range electrostatic interactions and nonspecific hydrophobic interactions between the peptide and domain. Our simulations highlight the important role of hydrophobic interactions in the entire SH3 binding process: both nonspecific hydrophobic contacts in the encounter complex and specific hydrophobic contacts in the fully bound complex. The encounter complex could be key to understanding the functional behavior of SH3 domain interactions because the encounter complex forms very quickly and the transition to the fully bound state is slower. In cells, an SH3 domain may form an encounter complex quickly and nonspecifically with many potential binding partners, allowing it to search for the correct recognition sequence before completing the binding process.


2019 ◽  
Author(s):  
Ruchi Lohia ◽  
Reza Salari ◽  
Grace Brannigan

<div>The role of electrostatic interactions and mutations that change charge states in intrinsically disordered proteins (IDPs) is well-established, but many disease-associated mutations in IDPs are charge-neutral. The Val66Met single nucleotide polymorphism (SNP) encodes a hydrophobic-to-hydrophobic mutation at the midpoint of the prodomain of precursor brain-derived neurotrophic factor (BDNF), one of the earliest SNPs to be associated with neuropsychiatric disorders, for which the underlying molecular mechanism is unknown. Here we report on over 250 μs of fully-atomistic, explicit solvent, temperature replica exchange molecular dynamics simulations of the 91 residue BDNF prodomain, for both the V66 and M66 sequence.</div><div>The simulations were able to correctly reproduce the location of both local and non-local secondary changes due to the Val66Met mutation when compared with NMR spectroscopy. We find that the local structure change is mediated via entropic and sequence specific effects. We show that the highly disordered prodomain can be meaningfully divided into domains based on sequence alone. Monte Carlo simulations of a self-excluding heterogeneous polymer, with monomers representing each domain, suggest the sequence would be effectively segmented by the long, highly disordered polyampholyte near the sequence midpoint. This is qualitatively consistent with observed interdomain contacts within the BDNF prodomain, although contacts between the two segments are enriched relative to the self-excluding polymer. The Val66Met mutation increases interactions across the boundary between the two segments, due in part to a specific Met-Met interaction with a Methionine in the other segment. This effect propagates to cause the non-local change in secondary structure around the second methionine, previously observed in NMR. The effect is not mediated simply via changes in inter-domain contacts but is also dependent on secondary structure formation around residue 66, indicating a mechanism for secondary structure coupling in disordered proteins. </div>


2021 ◽  
Vol 19 (1) ◽  
Author(s):  
Lasse Staby ◽  
Katrine Bugge ◽  
Rasmus Greve Falbe-Hansen ◽  
Edoardo Salladini ◽  
Karen Skriver ◽  
...  

Abstract Background Signal fidelity depends on protein–protein interaction–‘hubs’ integrating cues from large interactomes. Recently, and based on a common secondary structure motif, the αα-hubs were defined, which are small α-helical domains of large, modular proteins binding intrinsically disordered transcriptional regulators. Methods Comparative structural biology. Results We assign the harmonin-homology-domain (HHD, also named the harmonin N-terminal domain, NTD) present in large proteins such as harmonin, whirlin, cerebral cavernous malformation 2, and regulator of telomere elongation 1 to the αα-hubs. The new member of the αα-hubs expands functionality to include scaffolding of supra-modular complexes mediating sensory perception, neurovascular integrity and telomere regulation, and reveal novel features of the αα-hubs. As a common trait, the αα-hubs bind intrinsically disordered ligands of similar properties integrating similar cellular cues, but without cross-talk. Conclusion The inclusion of the HHD in the αα-hubs has uncovered new features, exemplifying the utility of identifying groups of hub domains, whereby discoveries in one member may cross-fertilize discoveries in others. These features make the αα-hubs unique models for decomposing signal specificity and fidelity. Using these as models, together with other suitable hub domain, we may advance the functional understanding of hub proteins and their role in cellular communication and signaling, as well as the role of intrinsically disordered proteins in signaling networks.


2021 ◽  
Vol 8 ◽  
Author(s):  
George V. Papamokos ◽  
George Tziatzos ◽  
Dimitrios G. Papageorgiou ◽  
Spyros Georgatos ◽  
Efthimios Kaxiras ◽  
...  

Protein phosphorylation is a key regulatory mechanism in eukaryotic cells. In the intrinsically disordered histone tails, phosphorylation is often a part of combinatorial post-translational modifications and an integral part of the “histone code” that regulates gene expression. Here, we study the association between two histone H3 tail peptides modified to different degrees, using fully atomistic molecular dynamics simulations. Assuming that the initial conformations are either α-helical or fully extended, we compare the propensity of the two peptides to associate with one another when both are unmodified, one modified and the other unmodified, or both modified. The simulations lead to the identification of distinct inter- and intramolecular interactions in the peptide dimer, highlighting a prominent role of a fine-tuned phosphorylation rheostat in peptide association. Progressive phosphorylation appears to modulate peptide charge, inducing strong and specific intermolecular interactions between the monomers, which do not result in the formation of amorphous or ordered aggregates, as documented by experimental evidence derived from Circular Dichroism and NMR spectroscopy. However, upon complete saturation of positive charges by phosphate groups, this effect is reversed: intramolecular interactions prevail and dimerization of zero-charge peptides is markedly reduced. These findings underscore the role of phosphorylation thresholds in the dynamics of intrinsically disordered proteins. Phosphorylation rheostats might account for the divergent effects of histone modifications on the modulation of chromatin structure.


2018 ◽  
Author(s):  
Himadri S. Samanta ◽  
Debayan Chakraborty ◽  
D. Thirumalai

Random polyampholytes (PAs) contain positively and negatively charged monomers that are distributed randomly along the polymer chain. The interaction between charges is assumed to be given by the Debye-Huckel potential. We show that the size of the PA is determined by an interplay between electrostatic interactions, giving rise to the polyelectrolyte (PE) effect due to net charge per monomer (σ), and an effective attractive PA interaction due to charge fluctuations, δσ. The interplay between these terms gives rise to non-monotonic dependence of the radius of gyration, Rg on the inverse Debye length, κ when PA effects are important . In the opposite limit, Rg decreases monotonically with increasing κ. Simulations of PA chains, using a charged bead-spring model, further corroborates our theoretical predictions. The simulations unambiguously show that conformational heterogeneity manifests itself among sequences that have identical PA parameters. A clear implication is that the phases of PA sequences, and by inference IDPs, cannot be determined using only the bare PA parameters (σ and δσ).The theory is used to calculate the changes in Rg on N, the number of residues for a set of Intrinsically Disordered Proteins (IDPs). For a certain class of IDPs, with N between 24 to 441, the size grows as Rg ~ N0.6, which agrees with data from Small Angle X-ray Scattering (SAXS) experiments.


2019 ◽  
Author(s):  
Ruchi Lohia ◽  
Reza Salari ◽  
Grace Brannigan

<p>The role of electrostatic interactions and mutations that change charge states in intrinsically disordered proteins (IDPs) is well-established, but many disease-associated mutations in IDPs are charge-neutral. The Val66Met single nucleotide polymorphism (SNP) in precursor brain-derived neurotrophic factor (BDNF) is one of the earliest SNPs to be associated with neuropsychiatric disorders, and the underlying molecular mechanism is unknown. Here we report on over 250 μs of fully-atomistic, explicit solvent, temperature replica exchange molecular dynamics (MD) simulations of the 91 residue BDNF prodomain, for both the V66 and M66 sequence. The simulations were able to correctly reproduce the location of both local and non-local secondary changes due to the Val66Met mutation when compared with NMR spectroscopy. We find that the change in local structure is mediated via entropic and sequence specific effects. We developed a hierarchical sequence-based framework for analysis and conceptualization, which first identifies “blobs” of 5-15 residues representing local globular regions or linkers. We use this framework within a novel test for enrichment of higher-order (tertiary) structure in disordered proteins; the size and shape of each blob is extracted from MD simulation of the real protein (RP), and used to parameterize a self-avoiding heterogenous polymer (SAHP). The SAHP version of the BDNF prodomain suggested a protein segmented into three regions, with a central long, highly disordered polyampholyte linker separating two globular regions. This effective segmentation was also observed in full simulations of the RP, but the Val66Met substitution significantly increased interactions across the linker, as well as the number of participating residues. The Val66Met substitution replaces β-bridging between Val66 and Val94 (on either side of the linker) with specific side-chain interactions between Met66 and Met95.The protein backbone in the vicinity of Met95 is then free to form β-bridges with residues 31-41 near the N-terminus, which condenses the protein. A significant role for Met/Met interactions is consistent with previously-observed non-local effects of the Val66Met SNP, as well as established interactions between the Met66 sequence and a Met-rich receptor that initiates neuronal growth cone retraction.</p>


Sign in / Sign up

Export Citation Format

Share Document