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

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.