Predicting the Sequence-Dependent Backbone Dynamics of Intrinsically Disordered Proteins

How the sequences of intrinsically disordered proteins (IDPs) code for functions is still an enigma. Dynamics, in particular residue-specific dynamics, holds crucial clues. Enormous efforts have been spent to characterize residue-specific dynamics of IDPs, mainly through NMR spin relaxation experiments. Here we present a sequence-based method, SeqDYN, for predicting residue-specific backbone dynamics of IDPs. SeqDYN employs a mathematical model with 21 parameters: one is a correlation length and 20 are the contributions of the amino acids to slow dynamics. Training on a set of 45 IDPs reveals aromatic, Arg, and long-branched aliphatic amino acids as the most active in slow dynamics whereas Gly and short polar amino acids as the least active. SeqDYN predictions not only provide an accurate and insightful characterization of sequence-dependent IDP dynamics but may also serve as indicators in a host of biophysical processes, including the propensities of IDP sequences to undergo phase separation.

Charge Relaying within a Phospho-Motif Rescue Binding Competency of a Disordered Transcription Factor

Structural disorder in proteins is central to cellular signaling, where conformational plasticity equips molecules to promiscuously interact with different partners. By engaging with multiple binding partners via the rearrangement of its three helices, the nuclear coactivator binding domain (NCBD) of the CBP/p300 transcription factor is a paradigmatic example of promiscuity. Recently, molecular simulations and experiments revealed that, through the establishment of long-range electrostatic interactions, intended as salt-bridges formed between the post-translationally inserted phosphate and positively charged residues in helix H3 of NCBD, phosphorylation triggers NCBD compaction, lowering its affinity for binding partners. By means of extensive molecular simulations, we here investigated the effect of short-range electrostatics on the conformational ensemble of NCBD, by monitoring the interactions between a phosphorylated serine and conserved positively charged residues within the NCBD phospho-motif. We found that empowering proximal electrostatic interactions, as opposed to long-range electrostatics, can reshape the NCBD ensemble rescuing the binding competency of phosphorylated NCBD. Given the conservation of positive charges in phospho-motifs, proximal electrostatic interactions might dampen the effects of phosphorylation and act as a relay to regulate phosphorylated intrinsically disordered proteins, ultimately tuning the binding affinity for different cellular partners.

Evaluation of implicit solvent models in molecular dynamics simulation of α-Synuclein

We report conventional and accelerated molecular dynamics simulations of α-Synuclein, designed to assess performance of using different starting conformation, solvation environment and force field combination. Backbone and sidechain chemical shifts, radius of gyration, presence of β-hairpin structures in KTK(E/Q)GV repeats and secondary structure percentages were used to evaluate how variations in forcefield, solvation model and simulation protocol provide results that correlate with experimental findings. We show that with suitable choice of forcefield and solvent, ff03ws and OBC implicit model, respectively, acceptable reproduction of experimental data on size and secondary structure is obtained by both conventional and accelerated MD. In contrast to the implicit solvent model, simulations in explicit TIP4P/2005 solvent do not properly represent size or secondary structure of α-Synuclein.Communicated by Ramaswamy H. Sarma.

Replica exchange molecular dynamics simulations reveal self-association sites in M-crystallin caused by mutations provide insights of cataract

Crystallins are ubiquitous, however, prevalence is seen in eye lens. Eye lens crystallins are long-lived and structural intactness is required for maintaining lens transparency and protein solubility. Mutations in crystallins often lead to cataract. In this study, we performed mutations at specific sites of M-crystallin, a close homologue of eye lens crystallin and studied by using replica exchange molecular dynamics simulation with generalized Born implicit solvent model. Mutations were made on the Ca²⁺ binding residues (K34D and S77D) and in the hydrophobic core (W45R) which is known to cause congenital cataract in homologous γD-crystallin. The chosen mutations caused large motion of the N-terminal Greek key, concomitantly broke the interlocking Greek keys interactions and perturbed the compact core resulting in several folded and partially unfolded states. Partially unfolded states exposed large hydrophobic patches that could act as precursors for self-aggregation. Accumulation of such aggregates is the potential cause of cataract in homologous eye lens crystallins.

Exploring the sequence-structure-function relationship for the intrinsically disordered βγ-crystallin Hahellin

βγ-Crystallins are a superfamily of proteins containing crystallin-type Greek key motifs. Some βγ-crystallin domains have been shown to bind Ca²⁺. Hahellin is a newly identified intrinsically disordered βγ-crystallin domain from Hahella chejuensis. It folds into a typical βγ-crystallin structure upon Ca²⁺ binding and acts as a Ca²⁺-regulated conformational switch. Besides Hahellin, another two putative βγ-crystallins from Caulobacter crescentus and Yersinia pestis are shown to be partially disordered in their apo-form and undergo large conformational changes upon Ca²⁺ binding, although whether they acquire a βγ-crystallin fold is not known. The extent of conformational disorder/order of a protein is determined by its amino acid sequence. To date how this sequence-structure relationship is reflected in the βγ-crystallin superfamily has not been investigated. In this work, we comparatively studied the sequence and structure of Hahellin with those of Protein S, an ordered βγ-crystallin, via various computational biophysical techniques. We found that several factors, including presence of a C-terminal disorder prone region, high content of energetic frustrations, and low contact density, may promote the formation of the disordered state of apo-Hahellin. We also analyzed the disorder propensities for other putative disordered βγ-crystallin domains. This study provides new clues for further understanding the sequence-structure-function relationship of βγ-crystallins.