Effects of pH on an IDP conformational ensemble explored by molecular dynamics simulation

Conformational Analyses of the AHD1-UBAN Region of TNIP1 Highlight Key Amino Acids for Interaction with Ubiquitin

Tumor necrosis factor ɑ (TNFɑ)-induced protein 3 (TNFAIP3)-interacting protein 1 (TNIP1) is genetically and functionally linked to limiting auto-immune and inflammatory responses. We have shown that TNIP1 (alias A20-binding inhibitor of NF-κB 1, ABIN1), functioning as a hub location to coordinate other proteins in repressing inflammatory signaling, aligns with biophysical traits indicative of its being an intrinsically disordered protein (IDP). IDPs move through a repertoire of three-dimensional structures rather than being in one set conformation. Here we employed bioinformatic analysis and biophysical interventions via amino acid mutations to assess and alter, respectively, conformational flexibility along a crucial region of TNIP1, encompassing the ABIN homology domain 1 and ubiquitin-binding domain in ABIN proteins and NEMO (AHD1-UBAN), by purposeful replacement of key residues. In vitro secondary structure measurements were mostly in line with, but not necessarily to the same degree as, expected results from in silico assessments. Notably, changes in single amino acids outside of the ubiquitin-binding region for gain-of-order effects had consequences along the length of the AHD1-UBAN propagating to that region. This is evidenced by differences in recognition of the partner protein polyubiquitin ≥ 28 residues away, depending on the mutation site, from the previously identified key binding site. These findings serve to demonstrate the role of conformational flexibility in protein partner recognition by TNIP1, thus identifying key amino acids likely to impact the molecular dynamics involved in TNIP1 repression of inflammatory signaling at large.

Electric field-induced conformational dynamics of CA9: a potential biomarker for glioblastoma multiforme

GBM, a malignant brain tumor prevalent in adults, can be treated using Electric field (EF) therapy. However, the underlying mechanism of EF-based GBM therapy is not well understood. In this study, we used bioinformatics and MD analysis to explore CA9 in EF therapy for GBM. CA9 was identified as a differentially expressed gene (DEG) sensitive to EF stimulation in GBM using GEO and TCGA for integrated analysis. Elevated CA9 expression was associated with reduced overall survival in GBM patients, indicating that CA9 was an adverse prognostic factor. Single-cell data demonstrated that CA9 expression was significantly higher in GBM cells than in normal cells, suggesting that CA9 could be an EF-sensitive biomarker for GBM. GSVA analysis suggested that CA9 was related to hypoxia and glucose metabolism in glioblastoma. MD simulations were employed to examine the impact of EF (0 V/nm ≤ E ≤ 0.5 V/nm) on the conformation of the CA9 protein, including RMSF, RMSD, Rg, secondary structure, and dipole moment. The CA9 protein structure was altered with different EF intensities, affecting the motion of protein atoms in an EF intensity-dependent manner. The number of hydrogen bonds was significantly reduced as the EF intensity increased, indicating that EF disrupted the hydrogen bonds. Additionally, the EF intensity affected the dipole moment and characteristic time. Besides, the CA9 gene family analysis suggested that this gene family was highly conserved. Overall, CA9 showed potential as a GBM biomarker sensitive to EF, presenting a prospective target for therapeutic interventions in EF-mediated GBM treatment.

Visualizing gaussian-chain like structural models of human α-synuclein in monomeric pre-fibrillar state: Solution SAXS data and modeling analysis

Here, using small angle X-ray scattering (SAXS) data profile as reference, we attempted to visualize conformational ensemble accessible prefibrillar monomeric state of α-synuclein in solution. In agreement with previous reports, our analysis also confirmed that α-synuclein molecules adopted disordered shape profile under non-associating conditions. Chain-ensemble modeling protocol with dummy residues provided two weighted averaged clusters of semi-extended shapes. Further, Ensemble Optimization Method (EOM) computed mole fractions of semi-extended "twisted" conformations which might co-exist in solution. Since these were only Cα traces of the models, ALPHAFOLD2 server was used to search for all-atom models. Comparison with experimental data showed all predicted models disagreed equally, as individuals. Finally, we employed molecular dynamics simulations and normal mode analysis-based search coupled with SAXS data to seek better agreeing models. Overall, our analysis concludes that a shifting equilibrium of curved models with low α-helical content best-represents non-associating monomeric α-synuclein.

Molecular crowding and amyloidogenic self-assembly: Emergent perspectives from modern computations

In recent decades, the conventional protein folding paradigm has been challenged by intriguing properties of disordered peptide sequences that do not adopt stably folded conformations. Such intrinsically disordered proteins and protein regions (IDPs and IDRs) are poised uniquely in biology due to their propensity for self-aggregation, amyloidogenesis, and correlations with a cluster of debilitating diseases. Complexities underlying their structural and functional manifestations are enhanced in the presence of molecular crowding via non-specific protein-protein and protein-solvent contacts. Enabled by technological advances, physics-based algorithms, and data science, modern computer simulations provide unprecedented insights into the structure, function, dynamics, and thermodynamics of complex macromolecular systems. These characteristics are frequently correlated and manifest into unique observables. This chapter presents an overview of how such methodologies can lend insights and drive investigations into the molecular trifecta of crowding, protein self-aggregation, and amyloidogenesis. It begins with a general overview of disordered proteins in relation to biological function and of a suite of relevant experimental methods. Specific examples are showcased in the biological context. This is followed by a description of the computational approaches that supplant experimental efforts, with an elaboration on enhanced molecular simulation methods. The chapter concludes by alluding to expanded possibilities in disease amelioration.

Molecular simulations integrated with experiments for probing the interaction dynamics and binding mechanisms of intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) exploit their plasticity to deploy a rich panoply of soft interactions and binding phenomena. Advances in tailoring molecular simulations for IDPs combined with experimental cross-validation offer an atomistic view of the mechanisms that control IDP binding, function, and dysfunction. The emerging theme is that unbound IDPs autonomously form transient local structures and self-interactions that determine their binding behavior. Recent results have shed light on whether and how IDPs fold, stay disordered or drive condensation upon binding; how they achieve binding specificity and select among competing partners. The disorder-binding paradigm is now being proactively used by researchers to target IDPs for rational drug design and engineer molecular responsive elements for biosensing applications.

Novel Covalent Modifier-Induced Local Conformational Changes within the Intrinsically Disordered Region of the Androgen Receptor

Intrinsically disordered regions (IDRs) of transcription factors play an important biological role in liquid condensate formation and gene regulation. It is thus desirable to investigate the druggability of IDRs and how small-molecule binders can alter their conformational stability. For the androgen receptor (AR), certain covalent ligands induce important changes, such as the neutralization of the condensate. To understand the specificity of ligand-IDR interaction and potential implications for the mechanism of neutralizing liquid-liquid phase separation (LLPS), we modeled and performed computer simulations of ligand-bound peptide segments obtained from the human AR. We analyzed how different covalent ligands affect local secondary structure, protein contact map, and protein-ligand contacts for these protein systems. We find that effective neutralizers make specific interactions (such as those between cyanopyrazole and tryptophan) that alter the helical propensity of the peptide segments. These findings on the mechanism of action can be useful for designing molecules that influence IDR structure and condensate of the AR in the future.

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.

Enhancing the Conformational Stability of the cl-Par-4 Tumor Suppressor via Site-Directed Mutagenesis

Intrinsically disordered proteins play important roles in cell signaling, and dysregulation of these proteins is associated with several diseases. Prostate apoptosis response-4 (Par-4), an approximately 40 kilodalton proapoptotic tumor suppressor, is a predominantly intrinsically disordered protein whose downregulation has been observed in various cancers. The caspase-cleaved fragment of Par-4 (cl-Par-4) is active and plays a role in tumor suppression by inhibiting cell survival pathways. Here, we employed site-directed mutagenesis to create a cl-Par-4 point mutant (D313K). The expressed and purified D313K protein was characterized using biophysical techniques, and the results were compared to that of the wild-type (WT). We have previously demonstrated that WT cl-Par-4 attains a stable, compact, and helical conformation in the presence of a high level of salt at physiological pH. Here, we show that the D313K protein attains a similar conformation as the WT in the presence of salt, but at an approximately two times lower salt concentration. This establishes that the substitution of a basic residue for an acidic residue at position 313 alleviates inter-helical charge repulsion between dimer partners and helps to stabilize the structural conformation.

Effect of Electric Field on α-Synuclein Fibrils: Revealed by Molecular Dynamics Simulations

The self-association of amylogenic proteins to the fibril form is considered a pivotal factor in the pathogenesis of neurodegenerative diseases, including Parkinson’s disease (PD). PD causes unintended or uncontrollable movements in its common symptoms. α-synuclein is the major cause of PD development and thus has been the main target of numerous studies to suppress and sequester its expression or effectively degrade it. Nonetheless, to date, there are no efficient and proven ways to prevent pathological protein aggregation. Recent investigations proposed applying an external electric field to interrupt the fibrils. This method is a non-invasive approach that has a certain benefit over others. We performed molecular dynamics (MD) simulations by applying an electric field on highly toxic fibrils of α-synuclein to gain a molecular-level insight into fibril disruption mechanisms. The results revealed that the applied external electric field induces substantial changes in the conformation of the α-synuclein fibrils. Furthermore, we show the threshold value for electric field strength required to completely disrupt the α-synuclein fibrils by opening the hydrophobic core of the fibril. Thus, our findings might serve as a valuable foundation to better understand molecular-level mechanisms of the α-synuclein fibrils disaggregation process under an applied external electric field.

pH Induced Switch in the Conformational Ensemble of Intrinsically Disordered Protein Prothymosin-α and Its Implications for Amyloid Fibril Formation

Aggregation of intrinsically disordered proteins (IDPs) can lead to neurodegenerative diseases. Although there is experimental evidence that acidic pH promotes IDP monomer compaction leading to aggregation, the general mechanism is unclear. We studied the pH effect on the conformational ensemble of prothymosin-α (proTα), which is involved in multiple essential functions, and probed its role in aggregation using computer simulations. We show that compaction in the proTα dimension at low pH is due to the protein's collapse in the intermediate region (E41-D80) rich in glutamic acid residues, enhancing its β-sheet content. We observed by performing dimer simulations that the conformations with high β-sheet content could act as aggregation-prone (N*) states and nucleate the aggregation process. The simulations initiated using N* states form dimers within a microsecond time scale, whereas the non-N* states do not form dimers within this time scale. This study contributes to understanding the general principles of pH-induced IDP aggregation.

Conformational Dynamics of Glucagon-like Peptide-2 with Different Electric Field

Molecular dynamics (MD) simulation was used to study the influence of electric field on Glucagon-like Peptide-2 (GLP-2). Different electric field strengths (0 V/nm ≤ E ≤ 1 V/nm) were mainly carried out on GLP-2. The structural changes in GLP-2 were analyzed by the Root Mean Square Deviation (RMSD), Root Mean Square Fluctuation (RMSF), Radius of Gyration (Rg), Solvent Accessible Surface Area (SASA), Secondary Structure and the number of hydrogen bonds. The stable α—helix structure of GLP-2 was unwound and transformed into an unstable Turn and Coil structure since the stability of the GLP-2 protein structure was reduced under the electric field. Our results show that the degree of unwinding of the GLP-2 structure was not linearly related to the electric field intensity. E = 0.5 V/nm was a special point where the degree of unwinding of the GLP-2 structure reached the maximum at this electric field strength. Under a weak electric field, E < 0.5 V/nm, the secondary structure of GLP-2 becomes loose, and the entropy of the chain increases. When E reaches a certain value (E > 0.5 V/nm), the electric force of the charged residues reaches equilibrium, along the z-direction. Considering the confinement of moving along another direction, the residue is less free. Thus, entropy decreases and enthalpy increases, which enhance the interaction of adjacent residues. It is of benefit to recover hydrogen bonds in the middle region of the protein. These investigations, about the effect of an electric field on the structure of GLP-2, can provide some theoretical basis for the biological function of GLP-2 in vivo.

Statistical Analysis of Protein-Ligand Interaction Patterns in Nuclear Receptor RORγ

The receptor RORγ belongs to the nuclear receptor superfamily that senses small signaling molecules and regulates at the gene transcription level. Since RORγ has a high basal activity and plays an important role in immune responses, inhibitors targeting this receptor have been a focus for many studies. The receptor-ligand interaction is complex, and often subtle differences in ligand structure can determine its role as an inverse agonist or an agonist. We examined more than 130 existing RORγ crystal structures that have the same receptor complexed with different ligands. We reported the features of receptor-ligand interaction patterns and the differences between agonist and inverse agonist binding. Specific changes in the contact interaction map are identified to distinguish active and inactive conformations. Further statistical analysis of the contact interaction patterns using principal component analysis reveals a dominant mode which separates allosteric binding vs. canonical binding and a second mode which may indicate active vs. inactive structures. We also studied the nature of constitutive activity by performing a 100-ns computer simulation of apo RORγ. Using constitutively active nuclear receptor CAR as a comparison, we identified a group of conserved contacts that have similar contact strength between the two receptors. These conserved contact interactions, especially a couple key contacts in H11–H12 interaction, can be considered essential to the constitutive activity of RORγ. These protein-ligand and internal protein contact interactions can be useful in the development of new drugs that direct receptor activity.

A modular approach to map out the conformational landscapes of unbound intrinsically disordered proteins

SignificanceIntrinsically disordered proteins have the unique ability to morph in response to multiple partners and thereby process sophisticated inputs and outputs. It is, however, a mystery whether their response is passive, that is, entirely determined by the partner, or controlled via an internal, yet unknown, folding mechanism. Here we introduce an approach to examine this key question and demonstrate its potential by dissecting the conformational properties of the partially disordered protein NCBD and obtaining important clues about how it performs its biological function.

Insights into the Interactions that Trigger the Primary Nucleation of Polymorphic α-Synuclein Dimers

Parkinson's disease is associated with the accumulation of α-synuclein (AS) aggregates that include polymorphic AS oligomers and polymorphic fibrils. There have been advances in solving the polymorphic state of AS fibrils, both by experimental techniques and molecular modeling tools. Yet, the polymorphic AS oligomers are now considered as the neurotoxic species, thus current and future studies making efforts to solve their structures at the molecular level. Importantly, it is crucial to explore the specific interactions between AS monomers within the dimer that stabilize the dimer and yield nucleation. Herein, we present a first work that probes at the molecular level the specific interactions between monomers in polymorphic AS dimers are derived from AS fibrils by applying molecular modeling tools. Our work reveals that both N-terminal and the non-amyloidogenic component domains play a role in the dimerization of all polymorphic AS dimers. In addition, helices along the N-terminal of AS monomers impede the contacts between AS monomers, thus preventing the nucleation or the dimerization of AS. This work provides insights into several mechanisms of the production of polymorphic AS dimers. Thus, the findings obtained in this work may assist in developing new therapeutic strategies for inhibiting the formation of the early-stage neurotoxic AS dimers.

Disordered regions tune order in chromatin organization and function

Intrinsically disordered proteins or hybrid proteins with ordered domains and disordered regions (both collectively designated as IDP(R)s) defy the well-established structure-function paradigm due to their ability to perform multiple biological functions even in the absence of a well-defined 3D structure. IDP(R)s have a unique ability to exist as a functional heterogeneous ensemble, where they adopt multiple thermodynamically stable conformations with low energy barriers between states. The resultant structural plasticity or conformational adaptability provides them with a high functional diversity and ease of regulation. Hence, IDP(R)s are highly efficient biological machinery to mediate intricate cellular functions such as signaling, gene expression, and assembly of complex structures. One such structure is the nucleoprotein complex known as Chromatin. Interestingly, the proteins involved in shaping up the structure and function of chromatin are abundant in disordered regions, which serve more than just as mere flexible linkers. The disordered regions are involved in crucial processes such as gene expression regulation, chromatin architecture maintenance, and liquid-liquid phase separation initiation. This review is an attempt to explore the advantages and the functional and regulatory roles of intrinsic disorder in several Chromatin Associated Proteins from a mechanistic standpoint.