scholarly journals Salt-bridge dynamics in intrinsically disordered proteins: A trade-off between electrostatic interactions and structural flexibility

2018 ◽  
Vol 1866 (5-6) ◽  
pp. 624-641 ◽  
Author(s):  
Sankar Basu ◽  
Parbati Biswas
2017 ◽  
Author(s):  
Sankar Basu ◽  
Parbati Biswas

AbstractIntrinsically Disordered Proteins (IDPs) are enriched in charged and polar residues; and, therefore, electrostatic interactions play a predominant role in their dynamics. In order to remain multi-functional and exhibit their characteristic binding promiscuity, they need to retain considerable dynamic flexibility. At the same time, they also need to accommodate a large number of oppositely charged residues, which eventually lead to the formation of salt-bridges, imparting local rigidity. The formation of salt-bridges therefore oppose the desired dynamic flexibility. Hence, there appears to be a meticulous trade-off between the two mechanisms which the current study attempts to unravel. With this objective, we identify and analyze salt-bridges, both as isolated as well as composite ionic bond motifs, in the molecular dynamic trajectories of a set of appropriately chosen IDPs. Time evolved structural properties of these salt-bridges like persistence, associated secondary structural ′order-disorder′ transitions, correlated atomic movements, contribution in the overall electrostatic balance of the proteins have been studied in necessary detail. The results suggest that the key to maintain such a trade-off over time is the continuous formation and dissolution of salt-bridges with a wide range of persistence. Also, the continuous dynamic interchange of charged-atom-pairs (coming from a variety of oppositely charged side-chains) in the transient ionic bonds supports a model of dynamic flexibility concomitant with the well characterized stochastic conformational switching in these proteins. The results and conclusions should facilitate the future design of salt-bridges as a mean to further explore the disordered-globular interface in proteins.


2020 ◽  
Author(s):  
Abhirup Bandyopadhyay ◽  
Sankar Basu

AbstractIntrinsically disordered proteins (IDP) serve as one of the key components in the global proteome. In contrast to the dominant class of cytosolic globular proteins, they harbor an enormous amount of physical flexibility and structural plasticity enforcing them to be retained in conformational ensembles rather than well defined stable folds. Previous studies in an aligned direction have revealed the importance of transient dynamical phenomena like that of saltbridge formation in IDPs to support their physical flexibility and have further highlighted their functional relevance. For this characteristic flexibility, IDPs remain amenable and accessible to different ordered binding partners, supporting their potential multi-functionality. The current study further addresses this complex structure-functional interplay in IDPs using phase transition dynamics to conceptualize the underlying (avalanche type) mechanism of their being distributed across and hopping around degenerate structural states (conformational ensembles). For this purpose, extensive molecular dynamics simulations have been done and the data analyzed from a statistical physics perspective. Investigation of the plausible scope ‘selforganized criticality’ (SOC) to fit into the complex dynamics of IDPs was found to be assertive, relating the conformational degeneracy of these proteins to their multi-functionality. In accordance with the transient nature of ‘salt-bridge dynamics’, the study further uses it as a probe to explain the structural basis of the proposed criticality in the conformational phase transition among self-similar groups in IDPs. The analysis reveal scale-invariant self-similar fractal geometries in structural conformations of different IDPs. Also, as discussed in the conclusion, the study has the potential to benefit structural tinkering of bio-medically relevant IDPs in the design of biotherapeutics against them.


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>


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>


2018 ◽  
Author(s):  
Bin Sun ◽  
Eric C. Cook ◽  
Trevor P. Creamer ◽  
Peter M. Kekenes-Huskey

calcineurin (CaN) is a serine/threonine phosphatase that regulates a variety of physiological and pathophysiological processes in mammalian tissue. The CaN regulatory domain (RD) is responsible for regulating the enzyme’s phosphatase activity, and is believed to be highly-disordered when inhibiting CaN, but undergoes a disorderto-order transition upon diffusion-limited binding with the regulatory protein calmodulin (CaM). The prevalence of polar and charged amino acids in the regulatory domain (RD) suggests electrostatic interactions are involved in mediating CaM binding, yet the lack of atomistic-resolution data for the bound complex has stymied efforts to probe how the RD sequence controls its conformational ensemble and long-range attractions contribute to target protein binding. In the present study, we investigated via computational modeling the extent to which electrostatics and structural disorder cofacilitate or hinder CaM/CaN association kinetics. Specifically, we examined several RD constructs that contain the CaM binding region (CAMBR) to characterize the roles of electrostatics versus conformational diversity in controlling diffusion-limited association rates, via microsecond-scale molecular dynamics (MD) and Brownian dynamic (BD) simulations. Our results indicate that the RD amino acid composition and sequence length influence both the dynamic availability of conformations amenable to CaM binding, as well as long-range electrostatic interactions to steer association. These findings provide intriguing insight into the interplay between conformational diversity and electrostatically-driven protein-protein association involving CaN, which are likely to extend to wide-ranging diffusion-limited processes regulated by intrinsically-disordered proteins.


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