scholarly journals Mechanistic insights from replica exchange molecular dynamics simulations into mutation induced disordered-to-ordered transition in Hahellin, a βγ-crystallin

2018 ◽  
Author(s):  
Sunita Patel ◽  
Bal Krishnan ◽  
Ramakrishna V. Hosur ◽  
Kandala V. R. Chary

AbstractIntrinsically disordered proteins (IDPs) form a special category because they lack a unique well-folded 3D structure under physiological conditions. They play crucial role in cell signaling, regulatory functions and responsible for several diseases. Although, they are abundant in nature, only a small fraction of it has been characterized till date. Such proteins adopt a range of conformations and can undergo transformation from disordered-to-ordered state or vice-versa upon binding to ligand. Insights of such conformational transition is perplexing in several cases. In the present study, we characterized disordered as well as ordered states and the factors contributing the transitions through a mutational study by employing replica exchange molecular dynamics simulation on a βγ-crystallin. Most of the proteins within this superfamily are inherently ordered. However, Hahellin, although a member of βγ-crystallin, it is intrinsically disordered in its apo-form which takes a well-ordered βγ-crystallin fold upon binding to Ca2+. It is intriguing that the mutation at the 5th position of the canonical motif to Arg increases the domain stability in several ordered microbial βγ-crystallins with concomitant loss in Ca2+ binding affinity. We carried out similar Ser to Arg mutations at 5th position of the canonical motif for the first time in an intrinsically disordered protein to understand the mechanistic insights of conformational transition. Our study revealed that newly formed ionic and hydrogen bonding interactions at the canonical Ca2+ binding sites play crucial role in transforming the disordered conformation into ordered βγ-crystallin.Author summaryIntrinsically disordered proteins lack a unique ordered 3D structure under physiological condition. Although, they are abundant in nature, only a small fraction of these proteins has been characterized till date due to adaptation of multiple conformations and methodological limitation. βγ-crystallins are inherently ordered, however recently a small number of proteins within this superfamily have been identified as intrinsically disordered protein. Hahellin is one such protein which is intrinsically disordered in its apo-form but takes a well-ordered βγ-crystallin fold upon binding to Ca2+. In the present study, we decipher the underlying mechanism of disordered-to-ordered transition in Hahellin by mutations, employing replica exchange molecular dynamics simulations. Earlier experimental studies reported an increase in stabilization of the ordered βγ-crystallion upon mutation to Arg at 5th position of the canonical Ca2+ binding motifs, N/D-N/D-X1-X2-S/T-S. We performed similar Ser to Arg mutation in an intrinsically disordered Hahellin to get the mechanistic insights of the conformational transition in the absence of Ca2+. Our study revealed that several newly formed ionic and hydrogen bonding interactions contributed by the mutant residues are responsible for both intra- and inter-motif rigidification, resulting in overall stability of βγ-crystallin domain.

2019 ◽  
Vol 116 (41) ◽  
pp. 20446-20452 ◽  
Author(s):  
Utsab R. Shrestha ◽  
Puneet Juneja ◽  
Qiu Zhang ◽  
Viswanathan Gurumoorthy ◽  
Jose M. Borreguero ◽  
...  

Intrinsically disordered proteins (IDPs) are abundant in eukaryotic proteomes, play a major role in cell signaling, and are associated with human diseases. To understand IDP function it is critical to determine their configurational ensemble, i.e., the collection of 3-dimensional structures they adopt, and this remains an immense challenge in structural biology. Attempts to determine this ensemble computationally have been hitherto hampered by the necessity of reweighting molecular dynamics (MD) results or biasing simulation in order to match ensemble-averaged experimental observables, operations that reduce the precision of the generated model because different structural ensembles may yield the same experimental observable. Here, by employing enhanced sampling MD we reproduce the experimental small-angle neutron and X-ray scattering profiles and the NMR chemical shifts of the disordered N terminal (SH4UD) of c-Src kinase without reweighting or constraining the simulations. The unbiased simulation results reveal a weakly funneled and rugged free energy landscape of SH4UD, which gives rise to a heterogeneous ensemble of structures that cannot be described by simple polymer theory. SH4UD adopts transient helices, which are found away from known phosphorylation sites and could play a key role in the stabilization of structural regions necessary for phosphorylation. Our findings indicate that adequately sampled molecular simulations can be performed to provide accurate physical models of flexible biosystems, thus rationalizing their biological function.


Life ◽  
2020 ◽  
Vol 10 (7) ◽  
pp. 109 ◽  
Author(s):  
Sandra S. Sullivan ◽  
Robert O.J. Weinzierl

Many of the proteins involved in key cellular regulatory events contain extensive intrinsically disordered regions that are not readily amenable to conventional structure/function dissection. The oncoprotein c-MYC plays a key role in controlling cell proliferation and apoptosis and more than 70% of the primary sequence is disordered. Computational approaches that shed light on the range of secondary and tertiary structural conformations therefore provide the only realistic chance to study such proteins. Here, we describe the results of several tests of force fields and water models employed in molecular dynamics simulations for the N-terminal 88 amino acids of c-MYC. Comparisons of the simulation data with experimental secondary structure assignments obtained by NMR establish a particular implicit solvation approach as highly congruent. The results provide insights into the structural dynamics of c-MYC1-88, which will be useful for guiding future experimental approaches. The protocols for trajectory analysis described here will be applicable for the analysis of a variety of computational simulations of intrinsically disordered proteins.


2021 ◽  
Author(s):  
Antonio B. Oliveira ◽  
Xingcheng Lin ◽  
Prakash Kulkarni ◽  
José N. Onuchic ◽  
Susmita Roy ◽  
...  

AbstractIntrinsically disordered proteins (IDPs) lack a rigid 3D structure and populate a polymorphic ensemble of conformations. Because of the lack of a reference conformation, their energy landscape representation in terms of reaction coordinates presents a daunting challenge. Here, our newly developed Energy Landscape Visualization Method (ELViM), a reaction coordinate-free approach, shows its prime application to explore frustrated energy landscapes of an intrinsically disordered protein, Prostate-Associated Gene 4 (PAGE4). PAGE4 is a transcriptional coactivator that potentiates the oncogene c-Jun. Two kinases, namely HIPK1 and CLK2, phosphorylate PAGE4 generating variants phosphorylated at different serine/threonine residues (HIPK1-PAGE4 and CLK2-PAGE4, respectively) with opposing functions. While HIPK1-PAGE4 predominantly phosphorylates Thr51 and potentiates c-Jun, CLK2-PAGE4 hyper-phosphorylates PAGE4 and attenuates transactivation. To understand the underlying mechanisms of conformational diversity among different phosphoforms, we have analyzed their atomistic trajectories simulated using AWSEM forcefield and the energy landscapes were elucidated using ELViM. This method allows us to identify and compare the population distributions of different conformational ensembles of PAGE4 phosphoforms using the same effective phase space. The results reveal a predominant conformational ensemble with an extended C-terminal segment of WT PAGE4, which exposes a functional residue Thr51, implying its potential of undertaking a fly-casting mechanism while binding to its cognate partner. In contrast, for HIPK1-PAGE4, a compact conformational ensemble enhances its population sequestering phosphorylated-Thr51. This clearly explains the experimentally observed weaker affinity of HIPK1-PAGE4 for c-Jun. ELViM appears as a powerful tool especially to analyze the highly-frustrated energy landscape representation of IDPs where appropriate reaction coordinates are hard to apprehend.


2020 ◽  
Vol 10 (1) ◽  
Author(s):  
Pablo Herrera-Nieto ◽  
Adrià Pérez ◽  
Gianni De Fabritiis

Abstract The exploration of intrinsically disordered proteins in isolation is a crucial step to understand their complex dynamical behavior. In particular, the emergence of partially ordered states has not been explored in depth. The experimental characterization of such partially ordered states remains elusive due to their transient nature. Molecular dynamics mitigates this limitation thanks to its capability to explore biologically relevant timescales while retaining atomistic resolution. Here, millisecond unbiased molecular dynamics simulations were performed in the exemplar N-terminal region of p53. In combination with state-of-the-art Markov state models, simulations revealed the existence of several partially ordered states accounting for $$\sim $$ ∼ 40% of the equilibrium population. Some of the most relevant states feature helical conformations similar to the bound structure of p53 to Mdm2, as well as novel $$\beta $$ β -sheet elements. This highlights the potential complexity underlying the energy surface of intrinsically disordered proteins.


Sign in / Sign up

Export Citation Format

Share Document