Using NMR Chemical Shifts to Determine Residue-Specific Secondary Structure Populations for Intrinsically Disordered Proteins

CoVAMPnet: Comparative Markov State Analysis for Studying Effects of Drug Candidates on Disordered Biomolecules

Computational study of the effect of drug candidates on intrinsically disordered biomolecules is challenging due to their vast and complex conformational space. Here, we developed a comparative Markov state analysis (CoVAMPnet) framework to quantify changes in the conformational distribution and dynamics of a disordered biomolecule in the presence and absence of small organic drug candidate molecules. First, molecular dynamics trajectories are generated using enhanced sampling, in the presence and absence of small molecule drug candidates, and ensembles of soft Markov state models (MSMs) are learned for each system using unsupervised machine learning. Second, these ensembles of learned MSMs are aligned across different systems based on a solution to an optimal transport problem. Third, the directional importance of inter-residue distances for the assignment to different conformational states is assessed by a discriminative analysis of aggregated neural network gradients. This final step provides interpretability and biophysical context to the learned MSMs. We applied this novel computational framework to assess the effects of ongoing phase 3 therapeutics tramiprosate (TMP) and its metabolite 3-sulfopropanoic acid (SPA) on the disordered Aβ42 peptide involved in Alzheimer's disease. Based on adaptive sampling molecular dynamics and CoVAMPnet analysis, we observed that both TMP and SPA preserved more structured conformations of Aβ42 by interacting nonspecifically with charged residues. SPA impacted Aβ42 more than TMP, protecting α-helices and suppressing the formation of aggregation-prone β-strands. Experimental biophysical analyses showed only mild effects of TMP/SPA on Aβ42 and activity enhancement by the endogenous metabolization of TMP into SPA. Our data suggest that TMP/SPA may also target biomolecules other than Aβ peptides. The CoVAMPnet method is broadly applicable to study the effects of drug candidates on the conformational behavior of intrinsically disordered biomolecules.

Disordered-to-ordered transitions in assembly factors allow the complex II catalytic subunit to switch binding partners

Complex II (CII) activity controls phenomena that require crosstalk between metabolism and signaling, including neurodegeneration, cancer metabolism, immune activation, and ischemia-reperfusion injury. CII activity can be regulated at the level of assembly, a process that leverages metastable assembly intermediates. The nature of these intermediates and how CII subunits transfer between metastable complexes remains unclear. In this work, we identify metastable species containing the SDHA subunit and its assembly factors, and we assign a preferred temporal sequence of appearance of these species during CII assembly. Structures of two species show that the assembly factors undergo disordered-to-ordered transitions without the appearance of significant secondary structure. The findings identify that intrinsically disordered regions are critical in regulating CII assembly, an observation that has implications for the control of assembly in other biomolecular complexes.

Protein disorder and autoinhibition: The role of multivalency and effective concentration

Regulation of protein binding through autoinhibition commonly occurs via interactions involving intrinsically disordered regions (IDRs). These intramolecular interactions can directly or allosterically inhibit intermolecular protein or DNA binding, regulate enzymatic activity, and control the assembly of large macromolecular complexes. Autoinhibitory interactions mediated by protein disorder are inherently transient, making their identification and characterization challenging. In this review, we explore the structural and functional diversity of disorder-mediated autoinhibition for a variety of biological mechanisms, with a focus on the role of multivalency and effective concentration. We also discuss the evolution of disordered motifs that participate in autoinhibition using examples where sequence conservation varies from high to low. In some cases, identifiable motifs that are essential for autoinhibition remain intact within a rapidly evolving sequence, over long evolutionary distances. Finally, we examine the potential of AlphaFold2 to predict autoinhibitory intramolecular interactions involving IDRs.

Complex Conformational Space of the RNA Polymerase II C-Terminal Domain upon Phosphorylation

Intrinsically disordered proteins (IDPs) have been closely studied during the past decade due to their importance in many biological processes. The disordered nature of this group of proteins makes it difficult to observe its full span of the conformational space using either experimental or computational studies. In this article, we explored the conformational space of the C-terminal domain (CTD) of RNA polymerase II (Pol II), which is also an intrinsically disordered low complexity domain, using enhanced sampling methods. We provided a detailed conformational analysis of model systems of CTD with different lengths; first with the last 44 residues of the human CTD sequence and finally the CTD model with 2-heptapeptide repeating units. We then investigated the effects of phosphorylation on CTD conformations by performing simulations at different phosphorylated states. We obtained broad conformational spaces in nonphosphorylated CTD models, and phosphorylation has complex effects on the conformations of the CTD. These complex effects depend on the length of the CTD, spacing between the multiple phosphorylation sites, ion coordination, and interactions with the nearby residues.

The influence of random-coil chemical shifts on the assessment of structural propensities in folded proteins and IDPs

In studying secondary structural propensities of proteins by nuclear magnetic resonance (NMR) spectroscopy, secondary chemical shifts (SCSs) serve as the primary atomic scale observables. For SCS calculation, the selection of an appropriate random coil chemical shift (RCCS) dataset is a crucial step, especially when investigating intrinsically disordered proteins (IDPs). The scientific literature is abundant in such datasets, however, the effect of choosing one over all the others in a concrete application has not yet been studied thoroughly and systematically. Hereby, we review the available RCCS prediction methods and to compare them, we conduct statistical inference by means of the nonparametric sum of ranking differences and comparison of ranks to random numbers (SRD-CRRN) method. We try to find the RCCS predictors best representing the general consensus regarding secondary structural propensities. The existence and the magnitude of resulting differences on secondary structure determination under varying sample conditions (temperature, pH) are demonstrated and discussed for globular proteins and especially IDPs.

Semaphorin 3A-Glycosaminoglycans Interaction as Therapeutic Target for Axonal Regeneration

Semaphorin 3A (Sema3A) is a cell-secreted protein that participates in the axonal guidance pathways. Sema3A acts as a canonical repulsive axon guidance molecule, inhibiting CNS regenerative axonal growth and propagation. Therefore, interfering with Sema3A signaling is proposed as a therapeutic target for achieving functional recovery after CNS injuries. It has been shown that Sema3A adheres to the proteoglycan component of the extracellular matrix (ECM) and selectively binds to heparin and chondroitin sulfate-E (CS-E) glycosaminoglycans (GAGs). We hypothesize that the biologically relevant interaction between Sema3A and GAGs takes place at Sema3A C-terminal polybasic region (SCT). The aims of this study were to characterize the interaction of the whole Sema3A C-terminal polybasic region (Sema3A 725–771) with GAGs and to investigate the disruption of this interaction by small molecules. Recombinant Sema3A basic domain was produced and we used a combination of biophysical techniques (NMR, SPR, and heparin affinity chromatography) to gain insight into the interaction of the Sema3A C-terminal domain with GAGs. The results demonstrate that SCT is an intrinsically disordered region, which confirms that SCT binds to GAGs and helps to identify the specific residues involved in the interaction. NMR studies, supported by molecular dynamics simulations, show that a new peptoid molecule (CSIC02) may disrupt the interaction between SCT and heparin. Our structural study paves the way toward the design of new molecules targeting these protein–GAG interactions with potential therapeutic applications.

HACANCOi: a new Hα-detected experiment for backbone resonance assignment of intrinsically disordered proteins

Unidirectional coherence transfer is highly efficient in intrinsically disordered proteins (IDPs). Their elevated ps-ns timescale dynamics ensures long transverse (T₂) relaxation times allowing sophisticated coherence transfer pathway selection in comparison to folded proteins. ¹H^α-detection ensures non-susceptibility to chemical exchange with the solvent and enables chemical shift assignment of consecutive proline residues, typically abundant in IDPs. However, many IDPs undergo a disorder-to-order transition upon interaction with their target protein, which leads to the loss of the favorable relaxation properties. Long coherence transfer routes now result in prohibitively large decrease in sensitivity. We introduce a novel 4D ¹H^α-detected experiment HACANCOi, together with its 3D implementation, which warrant high sensitivity for the assignment of proline-rich regions in IDPs in complex with a globular protein. The experiment correlates ¹H^α_i, ¹³C^α_i, ¹⁵N_i and [Formula: see text] spins by transferring the magnetization concomitantly from ¹³C^α_i to ¹⁵N_i and [Formula: see text]. The B1 domain of protein G (GB1), and the enteropathogenic E. coli EspF in complex with human SNX9 SH3, serve as model systems to demonstrate the attainable sensitivity and successful sequential assignment.

Disorder for Dummies: Functional Mutagenesis of Transient Helical Segments in Disordered Proteins

Most cytosolic eukaryotic proteins contain a mixture of ordered and disordered regions. Disordered regions facilitate cell signaling by concentrating sites for posttranslational modifications and protein-protein interactions into arrays of short linear motifs that can be reorganized by RNA splicing. The evolution of disordered regions looks different from their ordered counterparts. In some cases, selection is focused on maintaining protein binding interfaces and PTM sites, but sequence heterogeneity is common. In other cases, simple properties like charge, length, or end-to-end distance are maintained. Many disordered protein binding sites contain some transient secondary structure that may resemble the structure of the bound state. α-Helical secondary structure is common and a wide range of fractional helicity is observed in different disordered regions. Here we provide a simple protocol to identify transient helical segments and design mutants that can change their structure and function.

Dispersion from Cα or NH: 4D experiments for backbone resonance assignment of intrinsically disordered proteins

Resonance assignment of intrinsically disordered proteins is remarkably challenging due to scant chemical shift dispersion arising from conformational heterogeneity. The challenge is even greater if repeating segments are present in the amino acid sequence. To forward unambiguous resonance assignment of intrinsically disordered proteins, we present iHACANCO, HACACON and (HACA)CONCAHA, three H^α-detected 4D experiments with C^α as an additional dimension. In addition, we present (HACA)CON(CA)NH and (HACA)N(CA)CONH, new 4D H^α-start, H^N-detect experiments which have two N^H dimensions to enhance peak dispersion in a sequential walk through C', N^H and H^N, and provide more accurate N^H/H^N chemical shifts than those that can be obtained from a crowded ¹H, ¹⁵N-HSQC spectrum. Application of these 4D experiments is demonstrated using BilRI (165 aa), an outer-membrane intrinsically disordered protein from the opportunistic oral pathogen Aggregatibacter actinomycetemcomitans. BilRI amino acid sequence encompasses three very similar repeats with a 13-residue identical stretch in two of them.

p53 Phosphomimetics Preserve Transient Secondary Structure but Reduce Binding to Mdm2 and MdmX

The disordered p53 transactivation domain (p53TAD) contains specific levels of transient helical secondary structure that are necessary for its binding to the negative regulators, mouse double minute 2 (Mdm2) and MdmX. The interactions of p53 with Mdm2 and MdmX are also modulated by posttranslational modifications (PTMs) of p53TAD including phosphorylation at S15, T18 and S20 that inhibits p53-Mdm2 binding. It is unclear whether the levels of transient secondary structure in p53TAD are changed by phosphorylation or other PTMs. We used phosphomimetic mutants to determine if adding a negative charge at positions 15 and 18 has any effect on the transient secondary structure of p53TAD and protein-protein binding. Using a combination of biophysical and structural methods, we investigated the effects of single and multisite phosphomimetics on the transient secondary structure of p53TAD and its interaction with Mdm2, MdmX, and the KIX domain. The phosphomimetics reduced Mdm2 and MdmX binding affinity by 3–5-fold, but resulted in minimal changes in transient secondary structure, suggesting that the destabilizing effect of phosphorylation on the p53TAD-Mdm2 interaction is primarily electrostatic. Phosphomimetics had no effect on the p53-KIX interaction, suggesting that increased binding of phosphorylated p53 to KIX may be influenced by decreased competition with its negative regulators.