scholarly journals Heterogeneous dynamics in partially disordered proteins

2020 ◽  
Vol 22 (37) ◽  
pp. 21185-21196 ◽  
Author(s):  
Salla I. Virtanen ◽  
Anne M. Kiirikki ◽  
Kornelia M. Mikula ◽  
Hideo Iwaï ◽  
O. H. Samuli Ollila

Combination of novel isotopic labeling, NMR experiments and MD simulations reveal heterogeneous dynamics in partially disordered proteins.

2019 ◽  
Vol 9 (1) ◽  
Author(s):  
Tobias Schneider ◽  
Andrej Berg ◽  
Zeynel Ulusoy ◽  
Martin Gamerdinger ◽  
Christine Peter ◽  
...  

AbstractUbiquitylation is an eminent posttranslational modification referring to the covalent attachment of single ubiquitin molecules or polyubiquitin chains to a target protein dictating the fate of such labeled polypeptide chains. Here, we have biochemically produced artificially Lys11-, and Lys27-, and Lys63-linked ubiquitin dimers based on click-chemistry generating milligram quantities in high purity. We show that the artificial linkage used for the conjugation of two ubiquitin moieties represents a fully reliable surrogate of the natural isopeptide bond by acquiring highly resolved nuclear magnetic resonance (NMR) spectroscopic data including ligand binding studies. Extensive coarse grained and atomistic molecular dynamics (MD) simulations allow to extract structures representing the ensemble of domain-domain conformations used to verify the experimental data. Advantageously, this methodology does not require individual isotopic labeling of both ubiquitin moieties as NMR data have been acquired on the isotopically labeled proximal moiety and complementary MD simulations have been used to fully interpret the experimental data in terms of domain-domain conformation. This combined approach intertwining NMR spectroscopy with MD simulations makes it possible to describe the conformational space non-canonically Lys11-, and Lys27-linked ubiquitin dimers occupy in a solution averaged ensemble by taking atomically resolved information representing all residues in ubiquitin dimers into account.


2018 ◽  
Vol 35 (7) ◽  
pp. 1234-1236
Author(s):  
Yi He ◽  
Suhani Nagpal ◽  
Mourad Sadqi ◽  
Eva de Alba ◽  
Victor Muñoz

Abstract Motivation Many proteins are partially disordered in physiological conditions and only fold, fully or partially, upon binding. Their structural analysis is challenging because the accessible information, typically chemical shifts (CS) from nuclear magnetic resonance experiments, are averages over broad ensembles of conformations. We aim to develop a database for the analysis of such data in terms of conformational distributions of the protein backbone rather than of individual high-resolution structures. Results Glutton is the largest available database linking CS and protein 3D structures (5270 entries organized in three levels) and is searchable via a python script. It generates statistical distributions of ϕ−ψ dihedral angles based on CS or vice versa. Such ϕ−ψ distributions are used to calculate structural ensembles of partially disordered proteins from their CS. For folded proteins, such ensembles are excellent starting points for further refinement with additional experimental restraints (structure determination) or computational methods (structure prediction). Availability and implementation Glutton is freely available at https://github.com/YeeHo/Glutton. Supplementary information Supplementary data are available at Bioinformatics online.


2020 ◽  
Author(s):  
Alan Hicks ◽  
Cristian A. Escobar ◽  
Timothy A. Cross ◽  
Huan-Xiang Zhou

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.


2019 ◽  
Author(s):  
Joao Victor de Souza Cunha ◽  
Francesc Sabanes Zariquiey ◽  
Agnieszka K. Bronowska

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.


2020 ◽  
Author(s):  
Suman Samantray ◽  
Feng Yin ◽  
Batuhan Kav ◽  
Birgit Strodel

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.


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>


2020 ◽  
Vol 21 (17) ◽  
pp. 6166
Author(s):  
Joao V. de Souza ◽  
Francesc Sabanés Zariquiey ◽  
Agnieszka K. Bronowska

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.


Biomolecules ◽  
2020 ◽  
Vol 10 (6) ◽  
pp. 946 ◽  
Author(s):  
Alan Hicks ◽  
Cristian Escobar ◽  
Timothy Cross ◽  
Huan-Xiang Zhou

How sequences of intrinsically disordered proteins (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. MD simulations, first validated by SAXS and secondary chemical shift data, found scant α-helices or β-strands but a considerable propensity for polyproline II (PPII) torsion angles. Importantly, several blocks of residues (e.g., 11–29) emerge as “correlated segments”, identified by their 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 the correlated segment as a defining feature for the conformations and dynamics of IDPs.


Sign in / Sign up

Export Citation Format

Share Document