Constructing ensembles for intrinsically disordered proteins

The relatively flat energy landscapes associated with intrinsically disordered proteins makes modeling these systems especially problematic. A comprehensive model for these proteins requires one to build an ensemble consisting of a finite collection of structures, and their corresponding relative stabilities, which adequately capture the range of accessible states of the protein. In this regard, methods that use computational techniques to interpret experimental data in terms of such ensembles are an essential part of the modeling process. In this review, we critically assess the advantages and limitations of current techniques and discuss new methods for the validation of these ensembles

Comparative Studies of Disordered Proteins with Similar Sequences: Application to Aβ40 and Aβ42

Quantitative comparisons of intrinsically disordered proteins (IDPs) with similar sequences, such as mutant forms of the same protein, may provide insights into IDP aggregation—a process that plays a role in several neurodegenerative disorders. Here we describe an approach for modeling IDPs with similar sequences that simplifies the comparison of the ensembles by utilizing a single library of structures. The relative population weights of the structures are estimated using a Bayesian formalism, which provides measures of uncertainty in the resulting ensembles. We applied this approach to the comparison of ensembles for Aβ40 and Aβ42. Bayesian hypothesis testing finds that although both Aβ species sample β-rich conformations in solution that may represent prefibrillar intermediates, the probability that Aβ42 samples these prefibrillar states is roughly an order of magnitude larger than the frequency in which Aβ40 samples such structures. Moreover, the structure of the soluble prefibrillar state in our ensembles is similar to the experimentally determined structure of Aβ that has been implicated as an intermediate in the aggregation pathway. Overall, our approach for comparative studies of IDPs with similar sequences provides a platform for future studies on the effect of mutations on the structure and function of disordered proteins

Intrinsically Disordered Proteins: Where Computation Meets Experiment

Proteins are heteropolymers that play important roles in virtually every biological reaction. While many proteins have well-defined three-dimensional structures that are inextricably coupled to their function, intrinsically disordered proteins (IDPs) do not have a well-defined structure, and it is this lack of structure that facilitates their function. As many IDPs are involved in essential cellular processes, various diseases have been linked to their malfunction, thereby making them important drug targets. In this review we discuss methods for studying IDPs and provide examples of how computational methods can improve our understanding of IDPs. We focus on two intensely studied IDPs that have been implicated in very different pathologic pathways. The first, p53, has been linked to over 50% of human cancers, and the second, Amyloid-β (Aβ), forms neurotoxic aggregates in the brains of patients with Alzheimer’s disease. We use these representative proteins to illustrate some of the challenges associated with studying IDPs and demonstrate how computational tools can be fruitfully applied to arrive at a more comprehensive understanding of these fascinating heteropolymers.National Science Foundation (U.S.). Directorate for Biological Sciences. Postdoctoral Research Fellowship (Grant 1309247

ECG Morphological Variability in Beat Space for Risk Stratification After Acute Coronary Syndrome

Background: Identification of patients who are at high risk of adverse cardiovascular events after an acute coronary syndrome (ACS) remains a major challenge in clinical cardiology. We hypothesized that quantifying variability in electrocardiogram (ECG) morphology may improve risk stratification post‐ACS. Methods and Results: We developed a new metric to quantify beat‐to‐beat morphologic changes in the ECG: morphologic variability in beat space (MVB), and compared our metric to published ECG metrics (heart rate variability [HRV], deceleration capacity [DC], T‐wave alternans, heart rate turbulence, and severe autonomic failure). We tested the ability of these metrics to identify patients at high risk of cardiovascular death (CVD) using 1082 patients (1‐year CVD rate, 4.5%) from the MERLIN‐TIMI 36 (Metabolic Efficiency with Ranolazine for Less Ischemia in Non‐ST‐Elevation Acute Coronary Syndrome—Thrombolysis in Myocardial Infarction 36) clinical trial. DC, HRV/low frequency–high frequency, and MVB were all associated with CVD (hazard ratios [HRs] from 2.1 to 2.3 [P<0.05 for all] after adjusting for the TIMI risk score [TRS], left ventricular ejection fraction [LVEF], and B‐type natriuretic peptide [BNP]). In a cohort with low‐to‐moderate TRS (N=864; 1‐year CVD rate, 2.7%), only MVB was significantly associated with CVD (HR, 3.0; P=0.01, after adjusting for LVEF and BNP). Conclusions: ECG morphological variability in beat space contains prognostic information complementary to the clinical variables, LVEF and BNP, in patients with low‐to‐moderate TRS. ECG metrics could help to risk stratify patients who might not otherwise be considered at high risk of CVD post‐ACS

Hidden States within Disordered Regions of the CcdA Antitoxin Protein

The bacterial toxin–antitoxin system CcdB–CcdA provides a mechanism for the control of cell death and quiescence. The antitoxin protein CcdA is a homodimer composed of two monomers that each contain a folded N-terminal region and an intrinsically disordered C-terminal arm. Binding of the intrinsically disordered C-terminal arm of CcdA to the toxin CcdB prevents CcdB from inhibiting DNA gyrase and thereby averts cell death. Accurate models of the unfolded state of the partially disordered CcdA antitoxin can therefore provide insight into general mechanisms whereby protein disorder regulates events that are crucial to cell survival. Previous structural studies were able to model only two of three distinct structural states, a closed state and an open state, that are adopted by the C-terminal arm of CcdA. Using a combination of free energy simulations, single-pair Förster resonance energy transfer experiments, and existing NMR data, we developed structural models for all three states of the protein. Contrary to prior studies, we find that CcdA samples a previously unknown state where only one of the disordered C-terminal arms makes extensive contacts with the folded N-terminal domain. Moreover, our data suggest that previously unobserved conformational states play a role in regulating antitoxin concentrations and the activity of CcdA’s cognate toxin. These data demonstrate that intrinsic disorder in CcdA provides a mechanism for regulating cell fate

Crystallographic and Molecular Dynamics Analysis of Loop Motions Unmasking the Peptidoglycan-Binding Site in Stator Protein MotB of Flagellar Motor

Background: The C-terminal domain of MotB (MotB-C) shows high sequence similarity to outer membrane protein A and related peptidoglycan (PG)-binding proteins. It is believed to anchor the power-generating MotA/MotB stator unit of the bacterial flagellar motor to the peptidoglycan layer of the cell wall. We previously reported the first crystal structure of this domain and made a puzzling observation that all conserved residues that are thought to be essential for PG recognition are buried and inaccessible in the crystal structure. In this study, we tested a hypothesis that peptidoglycan binding is preceded by, or accompanied by, some structural reorganization that exposes the key conserved residues. Methodology/Principal Findings: We determined the structure of a new crystalline form (Form B) of Helicobacter pylori MotB-C. Comparisons with the existing Form A revealed conformational variations in the petal-like loops around the carbohydrate binding site near one end of the b-sheet. These variations are thought to reflect natural flexibility at this site required for insertion into the peptidoglycan mesh. In order to understand the nature of this flexibility we have performed molecular dynamics simulations of the MotB-C dimer. The results are consistent with the crystallographic data and provide evidence that the three loops move in a concerted fashion, exposing conserved MotB residues that have previously been implicated in binding of the peptide moiety of peptidoglycan. Conclusion/Significance: Our structural analysis provides a new insight into the mechanism by which MotB inserts into th

The Effect of a ΔK280 Mutation on the Unfolded State of a Microtubule-Binding Repeat in Tau

Tau is a natively unfolded protein that forms intracellular aggregates in the brains of patients with Alzheimer's disease. To decipher the mechanism underlying the formation of tau aggregates, we developed a novel approach for constructing models of natively unfolded proteins. The method, energy-minima mapping and weighting (EMW), samples local energy minima of subsequences within a natively unfolded protein and then constructs ensembles from these energetically favorable conformations that are consistent with a given set of experimental data. A unique feature of the method is that it does not strive to generate a single ensemble that represents the unfolded state. Instead we construct a number of candidate ensembles, each of which agrees with a given set of experimental constraints, and focus our analysis on local structural features that are present in all of the independently generated ensembles. Using EMW we generated ensembles that are consistent with chemical shift measurements obtained on tau constructs. Thirty models were constructed for the second microtubule binding repeat (MTBR2) in wild-type (WT) tau and a ΔK280 mutant, which is found in some forms of frontotemporal dementia. By focusing on structural features that are preserved across all ensembles, we find that the aggregation-initiating sequence, PHF6*, prefers an extended conformation in both the WT and ΔK280 sequences. In addition, we find that residue K280 can adopt a loop/turn conformation in WT MTBR2 and that deletion of this residue, which can adopt nonextended states, leads to an increase in locally extended conformations near the C-terminus of PHF6*. As an increased preference for extended states near the C-terminus of PHF6* may facilitate the propagation of β-structure downstream from PHF6*, these results explain how a deletion at position 280 can promote the formation of tau aggregates

The Assembly of Individual Chaplin Peptides from Streptomyces coelicolor into Functional Amyloid Fibrils

The self-association of proteins into amyloid fibrils offers an alternative to the natively folded state of many polypeptides. Although commonly associated with disease, amyloid fibrils represent the natural functional state of some proteins, such as the chaplins from the soil-dwelling bacterium Streptomyces coelicolor, which coat the aerial mycelium and spores rendering them hydrophobic. We have undertaken a biophysical characterisation of the five short chaplin peptides ChpD-H to probe the mechanism by which these peptides self-assemble in solution to form fibrils. Each of the five chaplin peptides produced synthetically or isolated from the cell wall is individually surface-active and capable of forming fibrils under a range of solution conditions in vitro. These fibrils contain a highly similar cross-β core structure and a secondary structure that resembles fibrils formed in vivo on the spore and mycelium surface. They can also restore the growth of aerial hyphae to a chaplin mutant strain. We show that cysteine residues are not required for fibril formation in vitro and propose a role for the cysteine residues conserved in four of the five short chaplin peptides

Elimination of Endogenous Toxin, Creatinine from Blood Plasma Depends on Albumin Conformation: Site Specific Uremic Toxicity & Impaired Drug Binding

Uremic syndrome results from malfunctioning of various organ systems due to the retention of uremic toxins which, under normal conditions, would be excreted into the urine and/or metabolized by the kidneys. The aim of this study was to elucidate the mechanisms underlying the renal elimination of uremic toxin creatinine that accumulate in chronic renal failure. Quantitative investigation of the plausible correlations was performed by spectroscopy, calorimetry, molecular docking and accessibility of surface area. Alkalinization of normal plasma from pH 7.0 to 9.0 modifies the distribution of toxin in the body and therefore may affect both the accumulation and the rate of toxin elimination. The ligand loading of HSA with uremic toxin predicts several key side chain interactions of site I that presumably have the potential to impact the specificity and impaired drug binding. These findings provide useful information for elucidating the complicated mechanism of toxin disposition in renal disease state