Sequence specificity despite intrinsic disorder: how a disease-associated Val/Met polymorphism rearranges tertiary interactions in a long disordered protein

Tertiary Structure ◽

Protein A ◽

Md Simulations ◽

Disordered Proteins ◽

Sequence Specificity ◽

Specific Effects ◽

Non Local

<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>

Sequence specificity despite intrinsic disorder: how a disease-associated Val/Met polymorphism rearranges tertiary interactions in a long disordered protein

10.26434/chemrxiv.8135777.v3 ◽

2019 ◽

Author(s):

Ruchi Lohia ◽

Reza Salari ◽

Grace Brannigan

Keyword(s):

Tertiary Structure ◽

Protein A ◽

Md Simulations ◽

Disordered Proteins ◽

Sequence Specificity ◽

Specific Effects ◽

Non Local

<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>

Sequence Specificity Despite Intrinsic Disorder: How a Disease-Associated Val/Met Polymorphism Shifts Tertiary Interactions in a Long Disordered Protein

10.26434/chemrxiv.8135777.v1 ◽

2019 ◽

Author(s):

Ruchi Lohia ◽

Reza Salari ◽

Grace Brannigan

Keyword(s):

Secondary Structure ◽

Disordered Proteins ◽

Sequence Specificity ◽

Specific Effects ◽

Heterogeneous Polymer ◽

Non Local ◽

Dynamics Simulations

<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>

Genetically tunable frustration controls allostery in an intrinsically disordered transcription factor

eLife ◽

10.7554/elife.30688 ◽

2017 ◽

Vol 6 ◽

Cited By ~ 37

Author(s):

Jing Li ◽

Jordan T White ◽

Harry Saavedra ◽

James O Wrabl ◽

Hesam N Motlagh ◽

...

Keyword(s):

Transcriptional Activation ◽

Structural Changes ◽

Tertiary Structure ◽

Control Function ◽

Disordered Proteins ◽

Novel Strategy ◽

Structured Proteins

Intrinsically disordered proteins (IDPs) present a functional paradox because they lack stable tertiary structure, but nonetheless play a central role in signaling, utilizing a process known as allostery. Historically, allostery in structured proteins has been interpreted in terms of propagated structural changes that are induced by effector binding. Thus, it is not clear how IDPs, lacking such well-defined structures, can allosterically affect function. Here, we show a mechanism by which an IDP can allosterically control function by simultaneously tuning transcriptional activation and repression, using a novel strategy that relies on the principle of ‘energetic frustration’. We demonstrate that human glucocorticoid receptor tunes this signaling in vivo by producing translational isoforms differing only in the length of the disordered region, which modulates the degree of frustration. We expect this frustration-based model of allostery will prove to be generally important in explaining signaling in other IDPs.

Charge fluctuation effects on the shape of flexible polyampholytes with applications to Intrinsically disordered proteins

10.1101/301911 ◽

2018 ◽

Author(s):

Himadri S. Samanta ◽

Debayan Chakraborty ◽

D. Thirumalai

Keyword(s):

Debye Length ◽

Disordered Proteins ◽

Radius Of Gyration ◽

Charge Fluctuation ◽

Net Charge ◽

X Ray Scattering ◽

Theoretical Predictions

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.

Sequence-dependent correlated segments in the intrinsically disordered region of ChiZ

10.1101/2020.04.22.055590 ◽

2020 ◽

Author(s):

Alan Hicks ◽

Cristian A. Escobar ◽

Timothy A. Cross ◽

Huan-Xiang Zhou

Keyword(s):

Conformational Dynamics ◽

Md Simulations ◽

Proton Exchange ◽

Disordered Proteins ◽

Amide Proton ◽

Relaxation Rates ◽

Amide Proton Exchange ◽

Sequence Dependent

AbstractIntrinsically disordered proteins (IDPs) account for a significant fraction of any proteome and are central to numerous cellular functions. Yet how sequences of IDPs code for their conformational dynamics is poorly understood. Here we combined NMR spectroscopy, small-angle X-ray scattering (SAXS), and molecular dynamics (MD) simulations to characterize the conformations and dynamics of ChiZ1-64. This IDP is the N-terminal fragment (residues 1-64) of the transmembrane protein ChiZ, a component of the cell division machinery in Mycobacterium tuberculosis. Its N-half contains most of the prolines and all of the anionic residues while the C-half most of the glycines and cationic residues. MD simulations, first validated by SAXS and secondary chemical shift data, found scant α-helices or β-strands but considerable propensity for polyproline II (PPII) torsion angles. Importantly, several blocks of residues (e.g., 11-29) emerge as “correlated segments”, identified by frequent formation of PPII stretches, salt bridges, cation-π interactions, and sidechain-backbone hydrogen bonds. NMR relaxation experiments showed non-uniform transverse relaxation rates (R2s) and nuclear Overhauser enhancements (NOEs) along the sequence (e.g., high R2s and NOEs for residues 11-14 and 23-28). MD simulations further revealed that the extent of segmental correlation is sequence-dependent: segments where internal interactions are more prevalent manifest elevated “collective” motions on the 5-10 ns timescale and suppressed local motions on the sub-ns timescale. Amide proton exchange rates provides corroboration, with residues in the most correlated segment exhibiting the highest protection factors. We propose correlated segment as a defining feature for the conformation and dynamics of IDPs.

Development of Charge-Augmented Three-Point Water Model (CAIPi3P) for Accurate Simulations of Intrinsically Disordered Proteins

10.26434/chemrxiv.7706867.v1 ◽

2019 ◽

Author(s):

Joao Victor de Souza Cunha ◽

Francesc Sabanes Zariquiey ◽

Agnieszka K. Bronowska

Keyword(s):

Experimental Data ◽

Molecular Dynamics ◽

Md Simulations ◽

Disordered Proteins ◽

Water Model ◽

High Plasticity ◽

Solvation Model ◽

Intrinsically Disordered Regions

Intrinsically disordered proteins (IDPs) are molecules without a fixed tertiary structure, exerting crucial roles in cellular signalling, growth and molecular recognition events. Due to their high plasticity, IDPs are very challenging in experimental and computational structural studies. To provide detailed atomic insight in IDPs dynamics governing its functional mechanisms, all-atom molecular dynamics (MD) simulations are widely employed. However, the current generalist force fields and solvent models are unable to generate satisfactory ensembles for IDPs when compared to existing experimental data. In this work, we present a new solvation model, denoted as Charge-Augmented 3 Point Water model for Intrinsically-disordered Proteins (CAIPi3P). CAIPi3P has been generated by performing a systematic scanning of atomic partial charges assigned to the widely popular molecular scaffold of the three-point TIP3P water model. We found that explicit solvent MD simulations employing CAIPi3P solvation considerably improved the SAXS scattering profiles for three different IDPs. Not surprisingly, this improvement was further enhanced by using CAIPi3P water in combination with the protein force field parametrized for IDPs. We have also demonstrated applicability of CAIPi3P to molecular systems containing structured as well as intrinsically disordered regions/domains. Our results highlight the crucial importance of solvent effects for generating molecular ensembles of IDPs which reproduce the experimental data available. Hence, we conclude that our newly developed CAIPi3P solvation model is a valuable tool assisting molecular simulations of intrinsically disordered proteins and assessing their molecular dynamics.

Aptamers Selected for Recognizing Amyloid β-protein—A Case for Cautious Optimism

10.20944/preprints201802.0047.v1 ◽

2018 ◽

Author(s):

Farid Rahimi

Keyword(s):

Scientific Evidence ◽

Protein A ◽

Intrinsically Disordered Protein ◽

Disordered Proteins ◽

Full Potential ◽

Amyloid Beta Protein ◽

Metastable Structures

Aptamers are versatile oligonucleotide ligands used for molecular recognition of diverse targets. However, application of aptamers to the field of amyloid β-protein (Aβ) has been limited so far. Aβ is an intrinsically disordered protein that exists in a dynamic conformational equilibrium, presenting time-dependent ensembles of short-lived, metastable structures and assemblies that have been generally difficult to isolate and characterize. Moreover, despite understanding of potential physiological roles of Aβ, this peptide has been linked to the pathogenesis of Alzheimer disease, and its pathogenic roles remain controversial. Accumulated scientific evidence thus far highlights undesirable or nonspecific interactions between selected aptamers and different Aβ assemblies likely due to metastable nature of Aβ or inherent affinity of RNA oligonucleotides to β-sheet-rich fibrillar structures of amyloidogenic proteins. Accordingly, lessons drawn from Aβ–aptamer studies emphasize that purity and uniformity of the protein target and rigorous characterization of aptamers’ specificity are important for realizing and garnering the full potential of aptamers selected for recongizing Aβ or other intrinsically disordered proteins. This review summarizes studies of aptamers selected for recognizing different Aβ assemblies and highlights controversies, difficulties, and limitations of such studies.

Different force fields give rise to different amyloid aggregation pathways in molecular dynamics simulations

10.1101/2020.09.09.290320 ◽

2020 ◽

Author(s):

Suman Samantray ◽

Feng Yin ◽

Batuhan Kav ◽

Birgit Strodel

Keyword(s):

Molecular Dynamics ◽

Force Field ◽

Force Fields ◽

Md Simulations ◽

Peptide Sequence ◽

Disordered Proteins ◽

Test Case ◽

Peptide Aggregation ◽

Intrinsically Disordered

AbstractThe progress towards understanding the molecular basis of Alzheimers’s disease is strongly connected to elucidating the early aggregation events of the amyloid-β (Aβ) peptide. Molecular dynamics (MD) simulations provide a viable technique to study the aggregation of Aβ into oligomers with high spatial and temporal resolution. However, the results of an MD simulation can only be as good as the underlying force field. A recent study by our group showed that none of the force fields tested can distinguish between aggregation-prone and non-aggregating peptide sequences, producing the same and in most cases too fast aggregation kinetics for all peptides. Since then, new force fields specially designed for intrinsically disordered proteins such as Aβ were developed. Here, we assess the applicability of these new force fields to studying peptide aggregation using the Aβ16−22 peptide and mutations of it as test case. We investigate their performance in modeling the monomeric state, the aggregation into oligomers, and the stability of the aggregation end product, i.e., the fibrillar state. A main finding is that changing the force field has a stronger effect on the simulated aggregation pathway than changing the peptide sequence. Also the new force fields are not able to reproduce the experimental aggregation propensity order of the peptides. Dissecting the various energy contributions shows that AMBER99SB-disp overestimates the interactions between the peptides and water, thereby inhibiting peptide aggregation. More promising results are obtained with CHARMM36m and especially its version with increased protein–water interactions. It is thus recommended to use this force field for peptide aggregation simulations and base future reparameterizations on it.

Development of Charge-Augmented Three-Point Water Model (CAIPi3P) for Accurate Simulations of Intrinsically Disordered Proteins

International Journal of Molecular Sciences ◽

10.3390/ijms21176166 ◽

2020 ◽

Vol 21 (17) ◽

pp. 6166

Author(s):

Joao V. de Souza ◽

Francesc Sabanés Zariquiey ◽

Agnieszka K. Bronowska

Keyword(s):

Experimental Data ◽

Molecular Dynamics ◽

Md Simulations ◽

Disordered Proteins ◽

Water Model ◽

High Plasticity ◽

Solvation Model ◽

Intrinsically Disordered Regions

Intrinsically disordered proteins (IDPs) are molecules without a fixed tertiary structure, exerting crucial roles in cellular signalling, growth and molecular recognition events. Due to their high plasticity, IDPs are very challenging in experimental and computational structural studies. To provide detailed atomic insight in IDPs’ dynamics governing their functional mechanisms, all-atom molecular dynamics (MD) simulations are widely employed. However, the current generalist force fields and solvent models are unable to generate satisfactory ensembles for IDPs when compared to existing experimental data. In this work, we present a new solvation model, denoted as the Charge-Augmented Three-Point Water Model for Intrinsically Disordered Proteins (CAIPi3P). CAIPi3P has been generated by performing a systematic scan of atomic partial charges assigned to the widely popular molecular scaffold of the three-point TIP3P water model. We found that explicit solvent MD simulations employing CAIPi3P solvation considerably improved the small-angle X-ray scattering (SAXS) scattering profiles for three different IDPs. Not surprisingly, this improvement was further enhanced by using CAIPi3P water in combination with the protein force field parametrized for IDPs. We also demonstrated the applicability of CAIPi3P to molecular systems containing structured as well as intrinsically disordered regions/domains. Our results highlight the crucial importance of solvent effects for generating molecular ensembles of IDPs which reproduce the experimental data available. Hence, we conclude that our newly developed CAIPi3P solvation model is a valuable tool for molecular simulations of intrinsically disordered proteins and assessing their molecular dynamics.