1
|
LeMaster DM, Bashir Q, Hernández G. Propagation of conformational instability in FK506-binding protein FKBP12. BIOCHIMICA ET BIOPHYSICA ACTA. PROTEINS AND PROTEOMICS 2024; 1872:140990. [PMID: 38142946 PMCID: PMC10939819 DOI: 10.1016/j.bbapap.2023.140990] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/27/2023] [Revised: 12/12/2023] [Accepted: 12/18/2023] [Indexed: 12/26/2023]
Abstract
FKBP12 is the archetype of the FK506 binding domains that define the family of FKBP proteins which participate in the regulation of various distinct physiological signaling processes. As the drugs FK506 and rapamycin inhibit many of these FKBP proteins, there is need to develop therapeutics which exhibit selectivity within this family. The long β4-β5 loop of the FKBP domain is known to regulate transcriptional activity for the steroid hormone receptors and appears to participate in regulating calcium channel activity for the cardiac and skeletal muscle ryanodine receptors. The β4-β5 loop of FKBP12 has been shown to undergo extensive conformational dynamics, and here we report hydrogen exchange measurements for a series of mutational variants in that loop which indicate deviations from a two-state kinetics for those dynamics. In addition to a previously characterized local transition near the tip of this loop, evidence is presented for a second site of conformational dynamics in the stem of this loop. These mutation-dependent hydrogen exchange effects extend beyond the β4-β5 loop, primarily by disrupting the hydrogen bond between the Gly 58 amide and the Tyr 80 carbonyl oxygen which links the two halves of the structural rim that surrounds the active site cleft. Mutationally-induced opening of the cleft between Gly 58 and Tyr 80 not only modulates the global stability of the protein, it promotes a conformational transition in the distant β2-β3a hairpin that modulates the binding affinity for a FKBP51-selective inhibitor previously designed to exploit a localized conformational transition at the homologous site.
Collapse
Affiliation(s)
- David M LeMaster
- Biggs Laboratory Wadsworth Center, NYS Department of Health, Empire State Plaza, Albany, NY 12237, United States of America
| | - Qamar Bashir
- Biggs Laboratory Wadsworth Center, NYS Department of Health, Empire State Plaza, Albany, NY 12237, United States of America
| | - Griselda Hernández
- Biggs Laboratory Wadsworth Center, NYS Department of Health, Empire State Plaza, Albany, NY 12237, United States of America.
| |
Collapse
|
2
|
Anderson JS, LeMaster DM, Hernández G. Transient conformations in the unliganded FK506 binding domain of FKBP51 correspond to two distinct inhibitor-bound states. J Biol Chem 2023; 299:105159. [PMID: 37579948 PMCID: PMC10514456 DOI: 10.1016/j.jbc.2023.105159] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2023] [Revised: 08/08/2023] [Accepted: 08/10/2023] [Indexed: 08/16/2023] Open
Abstract
Members of the FK506-binding protein (FKBP) family regulate a range of important physiological processes. Unfortunately, current therapeutics such as FK506 and rapamycin exhibit only modest selectivity among these functionally distinct proteins. Recent progress in developing selective inhibitors has been reported for FKBP51 and FKBP52, which act as mutual antagonists in the regulation of steroid hormone signaling. Two structurally similar inhibitors yield distinct protein conformations at the binding site. Localized conformational transition in the binding site of the unliganded FK1 domain of FKBP51 is suppressed by a K58T mutation that also suppresses the binding of these inhibitors. Here, it is shown that the changes in amide hydrogen exchange kinetics arising from this K58T substitution are largely localized to this structural region. Accurate determination of the hydroxide-catalyzed exchange rate constants in both the wildtype and K58T variant proteins impose strong constraints upon the pattern of amide exchange reactivities within either a single or a pair of transient conformations that could give rise to the differences between these two sets of measured rate constants. Poisson-Boltzmann continuum dielectric calculations provide moderately accurate predictions of the structure-dependent hydrogen exchange reactivity for solvent-exposed protein backbone amides. Applying such calculations to the local protein conformations observed in the two inhibitor-bound FKBP51 domains demonstrated that the experimentally determined exchange rate constants for the wildtype domain are robustly predicted by a population-weighted sum of the experimental hydrogen exchange reactivity of the K58T variant and the predicted exchange reactivities in model conformations derived from the two inhibitor-bound protein structures.
Collapse
Affiliation(s)
- Janet S Anderson
- Department of Chemistry, Union College, Schenectady, New York, USA
| | - David M LeMaster
- New York State Department of Health, Wadsworth Center, Albany, New York, USA
| | - Griselda Hernández
- New York State Department of Health, Wadsworth Center, Albany, New York, USA.
| |
Collapse
|
3
|
Devaurs D, Antunes DA, Borysik AJ. Computational Modeling of Molecular Structures Guided by Hydrogen-Exchange Data. JOURNAL OF THE AMERICAN SOCIETY FOR MASS SPECTROMETRY 2022; 33:215-237. [PMID: 35077179 DOI: 10.1021/jasms.1c00328] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
Data produced by hydrogen-exchange monitoring experiments have been used in structural studies of molecules for several decades. Despite uncertainties about the structural determinants of hydrogen exchange itself, such data have successfully helped guide the structural modeling of challenging molecular systems, such as membrane proteins or large macromolecular complexes. As hydrogen-exchange monitoring provides information on the dynamics of molecules in solution, it can complement other experimental techniques in so-called integrative modeling approaches. However, hydrogen-exchange data have often only been used to qualitatively assess molecular structures produced by computational modeling tools. In this paper, we look beyond qualitative approaches and survey the various paradigms under which hydrogen-exchange data have been used to quantitatively guide the computational modeling of molecular structures. Although numerous prediction models have been proposed to link molecular structure and hydrogen exchange, none of them has been widely accepted by the structural biology community. Here, we present as many hydrogen-exchange prediction models as we could find in the literature, with the aim of providing the first exhaustive list of its kind. From purely structure-based models to so-called fractional-population models or knowledge-based models, the field is quite vast. We aspire for this paper to become a resource for practitioners to gain a broader perspective on the field and guide research toward the definition of better prediction models. This will eventually improve synergies between hydrogen-exchange monitoring and molecular modeling.
Collapse
Affiliation(s)
- Didier Devaurs
- MRC Institute of Genetics and Cancer, University of Edinburgh, Edinburgh EH4 2XU, U.K
| | - Dinler A Antunes
- Department of Biology and Biochemistry, University of Houston, Houston, Texas 77005, United States
| | - Antoni J Borysik
- Department of Chemistry, King's College London, London SE1 1DB, U.K
| |
Collapse
|
4
|
Claesen J, Politis A. POPPeT: a New Method to Predict the Protection Factor of Backbone Amide Hydrogens. JOURNAL OF THE AMERICAN SOCIETY FOR MASS SPECTROMETRY 2019; 30:67-76. [PMID: 30338451 PMCID: PMC6318252 DOI: 10.1007/s13361-018-2068-x] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/13/2018] [Revised: 08/13/2018] [Accepted: 08/28/2018] [Indexed: 05/29/2023]
Abstract
Hydrogen exchange (HX) has become an important tool to monitor protein structure and dynamics. The interpretation of HX data with respect to protein structure requires understanding of the factors that influence exchange. Simulated protein structures can be validated by comparing experimental deuteration profiles with the profiles derived from the modeled protein structure. To do this, we propose here a new method, POPPeT, for protection factor prediction based on protein motions that enable HX. By comparing POPPeT with two existing methods, the phenomenological approximation and COREX, we show enhanced predictability measured at both protection factor and deuteration level. This method can be subsequently used by modeling strategies for protein structure prediction. Graphical Abstract ᅟ.
Collapse
Affiliation(s)
| | - Argyris Politis
- Department of Chemistry, King's College London, 7 Trinity Street, London, SE1 1DB, UK.
| |
Collapse
|
5
|
Devaurs D, Antunes DA, Papanastasiou M, Moll M, Ricklin D, Lambris JD, Kavraki LE. Coarse-Grained Conformational Sampling of Protein Structure Improves the Fit to Experimental Hydrogen-Exchange Data. Front Mol Biosci 2017; 4:13. [PMID: 28344973 PMCID: PMC5344923 DOI: 10.3389/fmolb.2017.00013] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/11/2016] [Accepted: 02/24/2017] [Indexed: 11/13/2022] Open
Abstract
Monitoring hydrogen/deuterium exchange (HDX) undergone by a protein in solution produces experimental data that translates into valuable information about the protein's structure. Data produced by HDX experiments is often interpreted using a crystal structure of the protein, when available. However, it has been shown that the correspondence between experimental HDX data and crystal structures is often not satisfactory. This creates difficulties when trying to perform a structural analysis of the HDX data. In this paper, we evaluate several strategies to obtain a conformation providing a good fit to the experimental HDX data, which is a premise of an accurate structural analysis. We show that performing molecular dynamics simulations can be inadequate to obtain such conformations, and we propose a novel methodology involving a coarse-grained conformational sampling approach instead. By extensively exploring the intrinsic flexibility of a protein with this approach, we produce a conformational ensemble from which we extract a single conformation providing a good fit to the experimental HDX data. We successfully demonstrate the applicability of our method to four small and medium-sized proteins.
Collapse
Affiliation(s)
- Didier Devaurs
- Department of Computer Science, Rice UniversityHouston, TX, USA
| | | | - Malvina Papanastasiou
- Department of Pathology and Laboratory Medicine, University of PennsylvaniaPhiladelphia, PA, USA
- Broad Institute of MIT & HarvardCambridge, MA, USA
| | - Mark Moll
- Department of Computer Science, Rice UniversityHouston, TX, USA
| | - Daniel Ricklin
- Department of Pathology and Laboratory Medicine, University of PennsylvaniaPhiladelphia, PA, USA
- Department of Pharmaceutical Sciences, University of BaselBasel, Switzerland
| | - John D. Lambris
- Department of Pathology and Laboratory Medicine, University of PennsylvaniaPhiladelphia, PA, USA
| | | |
Collapse
|
6
|
LeMaster DM, Hernandez G. Conformational Dynamics in FKBP Domains: Relevance to Molecular Signaling and Drug Design. Curr Mol Pharmacol 2016; 9:5-26. [PMID: 25986571 DOI: 10.2174/1874467208666150519113146] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2015] [Revised: 02/25/2015] [Accepted: 05/17/2015] [Indexed: 01/05/2023]
Abstract
Among the 22 FKBP domains in the human genome, FKBP12.6 and the first FKBP domains (FK1) of FKBP51 and FKBP52 are evolutionarily and structurally most similar to the archetypical FKBP12. As such, the development of inhibitors with selectivity among these four FKBP domains poses a significant challenge for structure-based design. The pleiotropic effects of these FKBP domains in a range of signaling processes such as the regulation of ryanodine receptor calcium channels by FKBP12 and FKBP12.6 and steroid receptor regulation by the FK1 domains of FKBP51 and FKBP52 amply justify the efforts to develop selective therapies. In contrast to their close structural similarities, these four FKBP domains exhibit a substantial diversity in their conformational flexibility. A number of distinct conformational transitions have been characterized for FKBP12 spanning timeframes from 20 s to 10 ns and in each case these dynamics have been shown to markedly differ from the conformational behavior for one or more of the other three FKBP domains. Protein flexibilitybased inhibitor design could draw upon the transitions that are significantly populated in only one of the targeted proteins. Both the similarities and differences among these four proteins valuably inform the understanding of how dynamical effects propagate across the FKBP domains as well as potentially how such intramolecular transitions might couple to the larger scale transitions that are central to the signaling complexes in which these FKBP domains function.
Collapse
Affiliation(s)
| | - Griselda Hernandez
- Wadsworth Center, New York State Department of Health, Empire State Plaza, Albany, New York, 12201, USA; Department of Biomedical Sciences, School of Public Health, University at Albany - SUNY, Empire State Plaza, Albany, New York, 12201, USA.
| |
Collapse
|
7
|
McAllister RG, Konermann L. Challenges in the Interpretation of Protein H/D Exchange Data: A Molecular Dynamics Simulation Perspective. Biochemistry 2015; 54:2683-92. [DOI: 10.1021/acs.biochem.5b00215] [Citation(s) in RCA: 52] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Robert G. McAllister
- Department
of Chemistry, The University of Western Ontario, London, Ontario N6A 5B7, Canada
| | - Lars Konermann
- Department
of Chemistry, The University of Western Ontario, London, Ontario N6A 5B7, Canada
| |
Collapse
|
8
|
Smith AE, Zhou LZ, Pielak GJ. Hydrogen exchange of disordered proteins in Escherichia coli. Protein Sci 2015; 24:706-13. [PMID: 25611326 DOI: 10.1002/pro.2643] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2015] [Accepted: 01/16/2015] [Indexed: 01/16/2023]
Abstract
A truly disordered protein lacks a stable fold and its backbone amide protons exchange with solvent at rates predicted from studies of unstructured peptides. We have measured the exchange rates of two model disordered proteins, FlgM and α-synuclein, in buffer and in Escherichia coli using the NMR experiment, SOLEXSY. The rates are similar in buffer and cells and are close to the rates predicted from data on small, unstructured peptides. This result indicates that true disorder can persist inside the crowded cellular interior and that weak interactions between proteins and macromolecules in cells do not necessarily affect intrinsic rates of exchange.
Collapse
Affiliation(s)
- Austin E Smith
- Department of Chemistry, University of North Carolina, Chapel Hill, North Carolina, 27599
| | | | | |
Collapse
|
9
|
Differential conformational dynamics in the closely homologous FK506-binding domains of FKBP51 and FKBP52. Biochem J 2014; 461:115-23. [PMID: 24749623 PMCID: PMC4060953 DOI: 10.1042/bj20140232] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Abstract
As co-chaperones of Hsp90 (heat-shock protein 90), FKBP51 (FK506-binding protein of 51 kDa) and FKBP52 (FK506-binding protein of 52 kDa) act as antagonists in regulating the hormone affinity and nuclear transport of steroid receptor complexes. Exchange of Leu119 in FKBP51 for Pro119 in FKBP52 has been shown to largely reverse the steroid receptor activities of FKBP51 and FKBP52. To examine whether differences in conformational dynamics/plasticity might correlate with changes in the reported receptor activities, 15N-NMR relaxation measurements were carried out on the N-terminal FKBP domains of FKBP51 and FKBP52 as well as their residue-swapped variants. Both proteins exhibit a similar pattern of motion in the picosecond–nanosecond timeframe as well as a small degree of 15N line-broadening, indicative of motion in the microsecond–millisecond timeframe, in the β3a strand of the central sheet. Only the FKBP51 domain exhibits much larger line-broadening in the adjacent β3 bulge (40′s loop of FKBP12) and throughout the long β4–β5 loop (80′s loop of FKBP12). The L119P mutation at the tip of the β4–β5 loop completely suppressed the line-broadening in this loop while partially suppressing the line-broadening in the neighbouring β2 and β3a strands. The complementary P119L and P119L/P124S variants of FKBP52 yielded similar patterns of line-broadening for the β4–β5 loop as that for FKBP51, although only 20% and 60% as intense respectively. However, despite the close structural similarity in the packing interactions between the β4–β5 loop and the β3a strand for FKBP51 and FKBP52, the line-broadening in the β3a strand is unaffected by the P119L or P119L/P124S mutations in FKBP52. Unlike FKBP52, the FK1 domain of FKBP51 exhibits microsecond–millisecond conformational dynamics in the β3 bulge and the β4–β5 loop, known sites of protein signalling interactions. Swapping residue 119 yields altered conformational dynamics in a pattern reminiscent of reported modulations in steroid receptor activity.
Collapse
|
10
|
Anderson JS, Mustafi SM, Hernández G, LeMaster DM. Statistical allosteric coupling to the active site indole ring flip equilibria in the FK506-binding domain. Biophys Chem 2014; 192:41-8. [PMID: 25016286 DOI: 10.1016/j.bpc.2014.06.004] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2014] [Accepted: 06/16/2014] [Indexed: 10/25/2022]
Abstract
In solution, the Trp 59 indole ring at the base of the active site cleft in the FKBP domain protein FKBP12 is rotated by ~90° at a population level of 20%, relative to its canonical crystallographic orientation. NMR measurements on the homologous FK1 domains of human FKBP51 and FKBP52 indicate no observable indole ring flip conformation, while the V101I variant of FKBP12 decreases the population having a perpendicular indole orientation by 10-fold. A set of three parallel 400 ns CHARMM27 molecular simulations for both wild type FKBP12 and the V101I variant examined how this ring flip might be energetically coupled to a transition of the Glu 60 sidechain which interacts with the backbone of the 50's loop located ~12 Å from the indole nitrogen. Analysis of the transition matrix for the local dynamics of the Glu 60 sidechain, the Trp 59 sidechain, and of the structurally interposed α-helix hydrogen bonding pattern yielded a statistical allosteric coupling of 10 kJ/mol with negligible concerted dynamical coupling for the transitions of the two sidechains.
Collapse
Affiliation(s)
- Janet S Anderson
- Department of Chemistry, Union College, Schenectady, NY 12308, United States
| | - Sourajit M Mustafi
- Wadsworth Center, New York State Department of Health, Empire State Plaza, Albany, NY 12201, United States
| | - Griselda Hernández
- Wadsworth Center, New York State Department of Health, Empire State Plaza, Albany, NY 12201, United States; Department of Biomedical Sciences, School of Public Health, University at Albany - SUNY, Empire State Plaza, Albany, NY 12201, United States
| | - David M LeMaster
- Wadsworth Center, New York State Department of Health, Empire State Plaza, Albany, NY 12201, United States; Department of Biomedical Sciences, School of Public Health, University at Albany - SUNY, Empire State Plaza, Albany, NY 12201, United States.
| |
Collapse
|
11
|
Structural basis of conformational transitions in the active site and 80's loop in the FK506-binding protein FKBP12. Biochem J 2014; 458:525-36. [PMID: 24405377 PMCID: PMC3940039 DOI: 10.1042/bj20131429] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
The extensive set of NMR doublings exhibited by the immunophilin FKBP12 (FK506-binding protein 12) arose from a slow transition to the cis-peptide configuration at Gly89 near the tip of the 80′s loop, the site for numerous protein-recognition interactions for both FKBP12 and other FKBP domain proteins. The 80′s loop also exhibited linebroadening, indicative of microsecond to millisecond conformational dynamics, but only in the trans-peptide state. The G89A variant shifted the trans–cis peptide equilibrium from 88:12 to 33:67, whereas a proline residue substitution induced fully the cis-peptide configuration. The 80′s loop conformation in the G89P crystal structure at 1.50 Å resolution differed from wild-type FKBP12 primarily at residues 88, 89 and 90, and it closely resembled that reported for FKBP52. Structure-based chemical-shift predictions indicated that the microsecond to millisecond dynamics in the 80′s loop probably arose from a concerted main chain (ψ88 and ϕ89) torsion angle transition. The indole side chain of Trp59 at the base of the active-site cleft was reoriented ~90o and the adjacent backbone was shifted in the G89P crystal structure. NOE analysis of wild-type FKBP12 demonstrated that this indole populates the perpendicular orientation at 20%. The 15N relaxation analysis was consistent with the indole reorientation occurring in the nanosecond timeframe. Recollection of the G89P crystal data at 1.20 Å resolution revealed a weaker wild-type-like orientation for the indole ring. Differences in the residues that underlie the Trp59 indole ring and altered interactions linking the 50′s loop to the active site suggested that reorientation of this ring may be disfavoured in the other six members of the FKBP domain family that bear this active-site tryptophan residue. Extensive resonance doubling arises from a cis–trans peptide transition at Gly89, whereas linebroadening appears due to a concerted shift in the neighbouring torsion angles. The active site Trp59 ring adopts a perpendicular orientation at a population of 20%.
Collapse
|
12
|
Chen H, Mustafi SM, LeMaster DM, Li Z, Héroux A, Li H, Hernández G. Crystal structure and conformational flexibility of the unligated FK506-binding protein FKBP12.6. ACTA CRYSTALLOGRAPHICA. SECTION D, BIOLOGICAL CRYSTALLOGRAPHY 2014; 70:636-46. [PMID: 24598733 PMCID: PMC3949516 DOI: 10.1107/s1399004713032112] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/09/2013] [Accepted: 11/25/2013] [Indexed: 12/15/2022]
Abstract
The primary known physiological function of FKBP12.6 involves its role in regulating the RyR2 isoform of ryanodine receptor Ca(2+) channels in cardiac muscle, pancreatic β islets and the central nervous system. With only a single previously reported X-ray structure of FKBP12.6, bound to the immunosuppressant rapamycin, structural inferences for this protein have been drawn from the more extensive studies of the homologous FKBP12. X-ray structures at 1.70 and 1.90 Å resolution from P2₁ and P3₁21 crystal forms are reported for an unligated cysteine-free variant of FKBP12.6 which exhibit a notable diversity of conformations. In one monomer from the P3₁21 crystal form, the aromatic ring of Phe59 at the base of the active site is rotated perpendicular to its typical orientation, generating a steric conflict for the immunosuppressant-binding mode. The peptide unit linking Gly89 and Val90 at the tip of the protein-recognition `80s loop' is flipped in the P2₁ crystal form. Unlike the >30 reported FKBP12 structures, the backbone conformation of this loop closely follows that of the first FKBP domain of FKBP51. The NMR resonances for 21 backbone amides of FKBP12.6 are doubled, corresponding to a slow conformational transition centered near the tip of the 80s loop, as recently reported for 31 amides of FKBP12. The comparative absence of doubling for residues along the opposite face of the active-site pocket in FKBP12.6 may in part reflect attenuated structural coupling owing to increased conformational plasticity around the Phe59 ring.
Collapse
Affiliation(s)
- Hui Chen
- Wadsworth Center, New York State Department of Health, Empire State Plaza, Albany, NY 12201, USA
| | - Sourajit M. Mustafi
- Wadsworth Center, New York State Department of Health, Empire State Plaza, Albany, NY 12201, USA
| | - David M. LeMaster
- Wadsworth Center, New York State Department of Health, Empire State Plaza, Albany, NY 12201, USA
- Department of Biomedical Sciences, School of Public Health, University at Albany – SUNY, Empire State Plaza, Albany, NY 12201, USA
| | - Zhong Li
- Wadsworth Center, New York State Department of Health, Empire State Plaza, Albany, NY 12201, USA
| | - Annie Héroux
- Department of Biology, Brookhaven National Laboratory, Upton, NY 11973, USA
| | - Hongmin Li
- Wadsworth Center, New York State Department of Health, Empire State Plaza, Albany, NY 12201, USA
- Department of Biomedical Sciences, School of Public Health, University at Albany – SUNY, Empire State Plaza, Albany, NY 12201, USA
| | - Griselda Hernández
- Wadsworth Center, New York State Department of Health, Empire State Plaza, Albany, NY 12201, USA
- Department of Biomedical Sciences, School of Public Health, University at Albany – SUNY, Empire State Plaza, Albany, NY 12201, USA
| |
Collapse
|
13
|
Kukic P, Farrell D, McIntosh LP, García-Moreno E B, Jensen KS, Toleikis Z, Teilum K, Nielsen JE. Protein dielectric constants determined from NMR chemical shift perturbations. J Am Chem Soc 2013; 135:16968-76. [PMID: 24124752 DOI: 10.1021/ja406995j] [Citation(s) in RCA: 68] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Understanding the connection between protein structure and function requires a quantitative understanding of electrostatic effects. Structure-based electrostatic calculations are essential for this purpose, but their use has been limited by a long-standing discussion on which value to use for the dielectric constants (ε(eff) and ε(p)) required in Coulombic and Poisson-Boltzmann models. The currently used values for ε(eff) and ε(p) are essentially empirical parameters calibrated against thermodynamic properties that are indirect measurements of protein electric fields. We determine optimal values for ε(eff) and ε(p) by measuring protein electric fields in solution using direct detection of NMR chemical shift perturbations (CSPs). We measured CSPs in 14 proteins to get a broad and general characterization of electric fields. Coulomb's law reproduces the measured CSPs optimally with a protein dielectric constant (ε(eff)) from 3 to 13, with an optimal value across all proteins of 6.5. However, when the water-protein interface is treated with finite difference Poisson-Boltzmann calculations, the optimal protein dielectric constant (ε(p)) ranged from 2 to 5 with an optimum of 3. It is striking how similar this value is to the dielectric constant of 2-4 measured for protein powders and how different it is from the ε(p) of 6-20 used in models based on the Poisson-Boltzmann equation when calculating thermodynamic parameters. Because the value of ε(p) = 3 is obtained by analysis of NMR chemical shift perturbations instead of thermodynamic parameters such as pK(a) values, it is likely to describe only the electric field and thus represent a more general, intrinsic, and transferable ε(p) common to most folded proteins.
Collapse
Affiliation(s)
- Predrag Kukic
- School of Biomolecular and Biomedical Science, Centre for Synthesis and Chemical Biology, UCD Conway Institute, University College Dublin , Belfield, Dublin 4, Ireland
| | | | | | | | | | | | | | | |
Collapse
|
14
|
Abstract
The 1H-15N 2D NMR correlation spectrum of the widely studied FK506-binding protein FKBP12 (FK506-binding protein of 12 kDa) contains previously unreported peak doublings for at least 31 residues that arise from a minor conformational state (12% of total) which exchanges with the major conformation with a time constant of 3.0 s at 43°C. The largest differences in chemical shift occur for the 80′s loop that forms critical recognition interactions with many of the protein partners for the FKBP family. The residues exhibiting doubling extend into the adjacent strands of the β-sheet, across the active site to the α-helix and into the 50′s loop. Each of the seven proline residues adopts a trans-peptide linkage in both the major and minor conformations, indicating that this slow transition is not the result of prolyl isomerization. Many of the residues exhibiting resonance doubling also participate in conformational line-broadening transition(s) that occur ~105-fold more rapidly, proposed previously to arise from a single global process. The 1.70 Å (1 Å=0.1 nm) resolution X-ray structure of the H87V variant is strikingly similar to that of FKBP12, yet this substitution quenches the slow conformational transition throughout the protein while quenching the line-broadening transition for residues near the 80′s loop. Line-broadening was also decreased for the residues in the α-helix and 50′s loop, whereas line-broadening in the 40′s loop was unaffected. The K44V mutation selectively reduces the line-broadening in the 40′s loop, verifying that at least three distinct conformational transitions underlie the line-broadening processes of FKBP12.
Collapse
|
15
|
Smith AE, Sarkar M, Young GB, Pielak GJ. Amide proton exchange of a dynamic loop in cell extracts. Protein Sci 2013; 22:1313-9. [PMID: 23904228 DOI: 10.1002/pro.2318] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2013] [Revised: 07/15/2013] [Accepted: 07/16/2013] [Indexed: 11/07/2022]
Abstract
Intrinsic rates of exchange are essential parameters for obtaining protein stabilities from amide (1) H exchange data. To understand the influence of the intracellular environment on stability, one must know the effect of the cytoplasm on these rates. We probed exchange rates in buffer and in Escherichia coli lysates for the dynamic loop in the small globular protein chymotrypsin inhibitor 2 using a modified form of the nuclear magnetic resonance experiment, SOLEXSY. No significant changes were observed, even in 100 g dry weight L(-1) lysate. Our results suggest that intrinsic rates from studies conducted in buffers are applicable to studies conducted under cellular conditions.
Collapse
Affiliation(s)
- Austin E Smith
- Department of Chemistry, University of North Carolina, Chapel Hill, North Carolina, 27599
| | | | | | | |
Collapse
|
16
|
Heitzer HM, Marks TJ, Ratner MA. First-Principles Calculation of Dielectric Response in Molecule-Based Materials. J Am Chem Soc 2013; 135:9753-9. [PMID: 23734640 DOI: 10.1021/ja401904d] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Affiliation(s)
- Henry M. Heitzer
- Department
of Chemistry and the Materials Research
Center, Northwestern University, 2145 Sheridan
Road, Evanston Illinois 60208, United States
| | - Tobin J. Marks
- Department
of Chemistry and the Materials Research
Center, Northwestern University, 2145 Sheridan
Road, Evanston Illinois 60208, United States
| | - Mark A. Ratner
- Department
of Chemistry and the Materials Research
Center, Northwestern University, 2145 Sheridan
Road, Evanston Illinois 60208, United States
| |
Collapse
|
17
|
Anderson JS, Hernández G, LeMaster DM. Assessing the chemical accuracy of protein structures via peptide acidity. Biophys Chem 2012. [PMID: 23182463 DOI: 10.1016/j.bpc.2012.10.005] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Although the protein native state is a Boltzmann conformational ensemble, practical applications often require a representative model from the most populated region of that distribution. The acidity of the backbone amides, as reflected in hydrogen exchange rates, is exquisitely sensitive to the surrounding charge and dielectric volume distribution. For each of four proteins, three independently determined X-ray structures of differing crystallographic resolution were used to predict exchange for the static solvent-exposed amide hydrogens. The average correlation coefficients range from 0.74 for ubiquitin to 0.93 for Pyrococcus furiosus rubredoxin, reflecting the larger range of experimental exchange rates exhibited by the latter protein. The exchange prediction errors modestly correlate with the crystallographic resolution. MODELLER 9v6-derived homology models at ~60% sequence identity (36% identity for chymotrypsin inhibitor CI2) yielded correlation coefficients that are ~0.1 smaller than for the cognate X-ray structures. The most recently deposited NOE-based ubiquitin structure and the original NMR structure of CI2 fail to provide statistically significant predictions of hydrogen exchange. However, the more recent RECOORD refinement study of CI2 yielded predictions comparable to the X-ray and homology model-based analyses.
Collapse
Affiliation(s)
- Janet S Anderson
- Department of Chemistry, Union College, Schenectady, New York 12308, USA
| | | | | |
Collapse
|
18
|
Skinner JJ, Lim WK, Bédard S, Black BE, Englander SW. Protein hydrogen exchange: testing current models. Protein Sci 2012; 21:987-95. [PMID: 22544567 DOI: 10.1002/pro.2082] [Citation(s) in RCA: 110] [Impact Index Per Article: 9.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2012] [Accepted: 04/09/2012] [Indexed: 11/06/2022]
Abstract
To investigate the determinants of protein hydrogen exchange (HX), HX rates of most of the backbone amide hydrogens of Staphylococcal nuclease were measured by NMR methods. A modified analysis was used to improve accuracy for the faster hydrogens. HX rates of both near surface and well buried hydrogens are spread over more than 7 orders of magnitude. These results were compared with previous hypotheses for HX rate determination. Contrary to a common assumption, proximity to the surface of the native protein does not usually produce fast exchange. The slow HX rates for unprotected surface hydrogens are not well explained by local electrostatic field. The ability of buried hydrogens to exchange is not explained by a solvent penetration mechanism. The exchange rates of structurally protected hydrogens are not well predicted by algorithms that depend only on local interactions or only on transient unfolding reactions. These observations identify some of the present difficulties of HX rate prediction and suggest the need for returning to a detailed hydrogen by hydrogen analysis to examine the bases of structure-rate relationships, as described in the companion paper (Skinner et al., Protein Sci 2012;21:996-1005).
Collapse
Affiliation(s)
- John J Skinner
- Johnson Research Foundation, Department of Biochemistry and Biophysics, Perelman School of Medicine, University of Pennsylvania, Philadelphia, Pennsylvania 19104-6059, USA.
| | | | | | | | | |
Collapse
|
19
|
Hernández G, Anderson JS, LeMaster DM. Experimentally assessing molecular dynamics sampling of the protein native state conformational distribution. Biophys Chem 2012; 163-164:21-34. [PMID: 22425325 DOI: 10.1016/j.bpc.2012.02.002] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2012] [Revised: 02/05/2012] [Accepted: 02/06/2012] [Indexed: 11/19/2022]
Abstract
The acute sensitivity to conformation exhibited by amide hydrogen exchange reactivity provides a valuable test for the physical accuracy of model ensembles developed to represent the Boltzmann distribution of the protein native state. A number of molecular dynamics studies of ubiquitin have predicted a well-populated transition in the tight turn immediately preceding the primary site of proteasome-directed polyubiquitylation Lys 48. Amide exchange reactivity analysis demonstrates that this transition is 10(3)-fold rarer than these predictions. More strikingly, for the most populated novel conformational basin predicted from a recent 1 ms MD simulation of bovine pancreatic trypsin inhibitor (at 13% of total), experimental hydrogen exchange data indicates a population below 10(-6). The most sophisticated efforts to directly incorporate experimental constraints into the derivation of model protein ensembles have been applied to ubiquitin, as illustrated by three recently deposited studies (PDB codes 2NR2, 2K39 and 2KOX2K392KOX). Utilizing the extensive set of experimental NOE constraints, each of these three ensembles yields a modestly more accurate prediction of the exchange rates for the highly exposed amides than does a standard unconstrained molecular simulation. However, for the less frequently exposed amide hydrogens, the 2NR2 ensemble offers no improvement in rate predictions as compared to the unconstrained MD ensemble. The other two NMR-constrained ensembles performed markedly worse, either underestimating (2KOX) or overestimating (2K39) the extent of conformational diversity.
Collapse
Affiliation(s)
- Griselda Hernández
- Wadsworth Center, New York State Department of Health, School of Public Health, University at Albany-SUNY, Empire State Plaza, Albany, NY 12201, USA
| | | | | |
Collapse
|
20
|
Hernández G, Anderson JS, Lemaster DM. Electrostatics of hydrogen exchange for analyzing protein flexibility. Methods Mol Biol 2012; 831:369-405. [PMID: 22167684 DOI: 10.1007/978-1-61779-480-3_20] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/31/2023]
Abstract
Electrostatic interactions at the protein-aqueous interface modulate the reactivity of solvent-exposed backbone amides by a factor of at least a billion fold. The brief (∼10 ps) lifetime of the peptide anion formed during the hydroxide-catalyzed exchange reaction helps enable the experimental rates to be robustly predictable by continuum dielectric methods. Since this ability to predict the structural dependence of exchange reactivity also applies to the protein amide hydrogens that are only rarely exposed to the bulk solvent phase, electrostatic analysis of the experimental exchange rates provides an effective assessment of whether a given model ensemble is consistent with the properly weighted Boltzmann conformational distribution of the protein native state.
Collapse
Affiliation(s)
- Griselda Hernández
- Department of Health and Department of Biomedical Sciences, Wadsworth Center, School of Public Health, University at Albany - SUNY, Albany, NY, USA
| | | | | |
Collapse
|
21
|
Zhu X, Lopes PE, MacKerell AD. Recent Developments and Applications of the CHARMM force fields. WILEY INTERDISCIPLINARY REVIEWS. COMPUTATIONAL MOLECULAR SCIENCE 2012; 2:167-185. [PMID: 23066428 PMCID: PMC3468154 DOI: 10.1002/wcms.74] [Citation(s) in RCA: 127] [Impact Index Per Article: 10.6] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Empirical force fields commonly used to describe the condensed phase properties of complex systems such as biological macromolecules are continuously being updated. Improvements in quantum mechanical (QM) methods used to generate target data, availability of new experimental target data, incorporation of new classes of compounds and new theoretical developments (eg. polarizable methods) make force-field development a dynamic domain of research. Accordingly, a number of improvements and extensions of the CHARMM force fields have occurred over the years. The objective of the present review is to provide an up-to-date overview of the CHARMM force fields. A limited presentation on the historical aspects of force fields will be given, including underlying methodologies and principles, along with a brief description of the strategies used for parameter development. This is followed by information on the CHARMM additive and polarizable force fields, including examples of recent applications of those force fields.
Collapse
Affiliation(s)
- Xiao Zhu
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, MD 21201
| | - Pedro E.M. Lopes
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, MD 21201
| | - Alexander D. MacKerell
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, MD 21201
| |
Collapse
|
22
|
Craig PO, Lätzer J, Weinkam P, Hoffman RMB, Ferreiro DU, Komives EA, Wolynes PG. Prediction of native-state hydrogen exchange from perfectly funneled energy landscapes. J Am Chem Soc 2011; 133:17463-72. [PMID: 21913704 DOI: 10.1021/ja207506z] [Citation(s) in RCA: 63] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Simulations based on perfectly funneled energy landscapes often capture many of the kinetic features of protein folding. We examined whether simulations based on funneled energy functions can also describe fluctuations in native-state protein ensembles. We quantitatively compared the site-specific local stability determined from structure-based folding simulations, with hydrogen exchange protection factors measured experimentally for ubiquitin, chymotrypsin inhibitor 2, and staphylococcal nuclease. Different structural definitions for the open and closed states based on the number of native contacts for each residue, as well as the hydrogen-bonding state, or a combination of both criteria were evaluated. The predicted exchange patterns agree with the experiments under native conditions, indicating that protein topology indeed has a dominant effect on the exchange kinetics. Insights into the simplest mechanistic interpretation of the amide exchange process were thus obtained.
Collapse
Affiliation(s)
- Patricio O Craig
- Department of Chemistry and Biochemistry, University of California at San Diego, La Jolla, California 92093-0374, USA
| | | | | | | | | | | | | |
Collapse
|
23
|
Ponomarev SY, Kaminski GA. Polarizable Simulations with Second order Interaction Model (POSSIM) force field: Developing parameters for alanine peptides and protein backbone. J Chem Theory Comput 2011; 7:1415-1427. [PMID: 21743799 DOI: 10.1021/ct1007197] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Abstract
A previously introduced POSSIM (POlarizable Simulations with Second order Interaction Model) force field has been extended to include parameters for alanine peptides and protein backbones. New features were introduced into the fitting protocol, as compared to the previous generation of the polarizable force field for proteins. A reduced amount of quantum mechanical data was employed in fitting the electrostatic parameters. Transferability of the electrostatics between our recently developed NMA model and the protein backbone was confirmed. Binding energy and geometry for complexes of alanine dipeptide with a water molecule were estimated and found in a good agreement with high-level quantum mechanical results (for example, the intermolecular distances agreeing within ca. 0.06Å). Following the previously devised procedure, we calculated average errors in alanine di- and tetra-peptide conformational energies and backbone angles and found the agreement to be adequate (for example, the alanine tetrapeptide extended-globular conformational energy gap was calculated to be 3.09 kcal/mol quantim mechanically and 3.14 kcal/mol with the POSSIM force field). However, we have now also included simulation of a simple alpha-helix in both gas-phase and water as the ultimate test of the backbone conformational behavior. The resulting alanine and protein backbone force field is currently being employed in further development of the POSSIM fast polarizable force field for proteins.
Collapse
Affiliation(s)
- Sergei Y Ponomarev
- Department of Chemistry and Biochemistry, Worcester Polytechnic Institute, Worcester, MA 01609
| | | |
Collapse
|
24
|
Shaw BF, Arthanari H, Narovlyansky M, Durazo A, Frueh DP, Pollastri MP, Lee A, Bilgicer B, Gygi SP, Wagner G, Whitesides GM. Neutralizing positive charges at the surface of a protein lowers its rate of amide hydrogen exchange without altering its structure or increasing its thermostability. J Am Chem Soc 2010; 132:17411-25. [PMID: 21090618 DOI: 10.1021/ja9067035] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
This paper combines two techniques--mass spectrometry and protein charge ladders--to examine the relationship between the surface charge and hydrophobicity of a representative globular protein (bovine carbonic anhydrase II; BCA II) and its rate of amide hydrogen-deuterium (H/D) exchange. Mass spectrometric analysis indicated that the sequential acetylation of surface lysine-ε-NH3(+) groups--a type of modification that increases the net negative charge and hydrophobicity of the surface of BCA II without affecting its secondary or tertiary structure--resulted in a linear decrease in the aggregate rate of amide H/D exchange at pD 7.4, 15 °C. According to analysis with MS, the acetylation of each additional lysine generated between 1.4 and 0.9 additional hydrogens that are protected from H/D exchange during the 2 h exchange experiment at 15 °C, pD 7.4. NMR spectroscopy demonstrated that none of the hydrogen atoms which became protected upon acetylation were located on the side chain of the acetylated lysine residues (i.e., lys-ε-NHCOCH3) but were instead located on amide NHCO moieties in the backbone. The decrease in rate of exchange associated with acetylation paralleled a decrease in thermostability: the most slowly exchanging rungs of the charge ladder were the least thermostable (as measured by differential scanning calorimetry). This observation--that faster rates of exchange are associated with slower rates of denaturation--is contrary to the usual assumptions in protein chemistry. The fact that the rates of H/D exchange were similar for perbutyrated BCA II (e.g., [lys-ε-NHCO(CH2)2CH3]18) and peracetylated BCA II (e.g., [lys-ε-NHCOCH3]18) suggests that the electrostatic charge is more important than the hydrophobicity of surface groups in determining the rate of H/D exchange. These electrostatic effects on the kinetics of H/D exchange could complicate (or aid) the interpretation of experiments in which H/D exchange methods are used to probe the structural effects of non-isoelectric perturbations to proteins (i.e., phosphorylation, acetylation, or the binding of the protein to an oligonucleotide or to another charged ligand or protein).
Collapse
Affiliation(s)
- Bryan F Shaw
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, USA.
| | | | | | | | | | | | | | | | | | | | | |
Collapse
|
25
|
Hernández G, Anderson JS, LeMaster DM. Assessing the native state conformational distribution of ubiquitin by peptide acidity. Biophys Chem 2010; 153:70-82. [PMID: 21055867 DOI: 10.1016/j.bpc.2010.10.006] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2010] [Revised: 10/08/2010] [Accepted: 10/10/2010] [Indexed: 11/18/2022]
Abstract
At equilibrium, every energetically feasible conformation of a protein occurs with a non-zero probability. Quantitative analysis of protein flexibility is thus synonymous with determining the proper Boltzmann-weighting of this conformational distribution. The exchange reactivity of solvent-exposed amide hydrogens greatly varies with conformation, while the short-lived peptide anion intermediate implies an insensitivity to the dynamics of conformational motion. Amides that are well-exposed in model conformational ensembles of ubiquitin vary a million-fold in exchange rates which continuum dielectric methods can predict with an rmsd of 3. However, the exchange rates for many of the more rarely exposed amides are markedly overestimated in the PDB-deposited 2K39 and 2KN5 ubiquitin ensembles, while the 2NR2 ensemble predictions are largely consistent with those of the Boltzmann-weighted conformational distribution sampled at the level of 1%. The correlation between the fraction of solvent-accessible conformations for a given amide hydrogen and the exchange rate constant for that residue provides a useful monitor of the degree of completeness with which a given ensemble has sampled the energetically accessible conformational space. These exchange predictions correlate with the degree to which each ensemble deviates from a set of 46 ubiquitin X-ray structures. Kolmogorov-Smirnov analysis for the distribution of intra- and inter-ensemble pairwise structural rmsd values assisted the identification of a subensemble of 2K39 that eliminates the overestimations of hydrogen exchange rates observed for the full ensemble. The relative merits of incorporating experimental restraints into the conformational sampling process are compared to using these restraints as filters to select subpopulations consistent with the experimental data.
Collapse
Affiliation(s)
- Griselda Hernández
- Wadsworth Center, New York State Department of Health, University at Albany-SUNY, 12201, USA
| | | | | |
Collapse
|
26
|
Aleksandrov A, Polydorides S, Archontis G, Simonson T. Predicting the Acid/Base Behavior of Proteins: A Constant-pH Monte Carlo Approach with Generalized Born Solvent. J Phys Chem B 2010; 114:10634-48. [DOI: 10.1021/jp104406x] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Alexey Aleksandrov
- Laboratoire de Biochimie (CNRS UMR7654), Department of Biology, Ecole Polytechnique, 91128 Palaiseau, France, and Department of Physics, University of Cyprus, PO20537, CY1678, Nicosia, Cyprus
| | - Savvas Polydorides
- Laboratoire de Biochimie (CNRS UMR7654), Department of Biology, Ecole Polytechnique, 91128 Palaiseau, France, and Department of Physics, University of Cyprus, PO20537, CY1678, Nicosia, Cyprus
| | - Georgios Archontis
- Laboratoire de Biochimie (CNRS UMR7654), Department of Biology, Ecole Polytechnique, 91128 Palaiseau, France, and Department of Physics, University of Cyprus, PO20537, CY1678, Nicosia, Cyprus
| | - Thomas Simonson
- Laboratoire de Biochimie (CNRS UMR7654), Department of Biology, Ecole Polytechnique, 91128 Palaiseau, France, and Department of Physics, University of Cyprus, PO20537, CY1678, Nicosia, Cyprus
| |
Collapse
|
27
|
Anderson JS, Hernández G, LeMaster DM. Sidechain conformational dependence of hydrogen exchange in model peptides. Biophys Chem 2010; 151:61-70. [PMID: 20627534 DOI: 10.1016/j.bpc.2010.05.006] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2010] [Revised: 05/10/2010] [Accepted: 05/12/2010] [Indexed: 10/19/2022]
Abstract
Peptide hydrogens that are exposed to solvent in protein X-ray structures exhibit a billion-fold range in hydroxide-catalyzed exchange rates, and these rates have previously been shown to be predictable by continuum dielectric methods to within a factor of 7, based on single protein conformations. When using a protein coil library to model the Boltzmann-weighted conformational distribution for the various N-acetyl-[X-Ala]-N-methylamides and N-acetyl-[Ala-Y]-N-methylamides, the acidity of the central amide in the individual conformers of each peptide spans nearly a million-fold range. Nevertheless, population averaging of these conformer acidities predicts the standard sidechain-dependent hydrogen exchange correction factors for nonpolar model peptides to within a factor of 30% (10(0.11)) with a correlation coefficient r=0.91. Comparison with the analogous continuum dielectric calculations for the other N-acetyl-[X-Y]-N-methylamides indicates that deviations from the isolated residue hypothesis of classical polymer theory predict appreciable errors in the exchange rates for conformationally disordered peptides when the standard sidechain-dependent hydrogen exchange rate correction factors are assumed to be independently additive. Although electronic polarizability generally dominates the dielectric shielding for the approximately 10ps lifetime of peptide ionization, evidence is presented for modest contributions from rapid intrarotamer conformational reorganization of Asn and Gln sidechains.
Collapse
Affiliation(s)
- Janet S Anderson
- Department of Chemistry, Union College, Schenectady, NY 12308, USA
| | | | | |
Collapse
|
28
|
LeMaster DM, Anderson JS, Hernández G. Peptide conformer acidity analysis of protein flexibility monitored by hydrogen exchange. Biochemistry 2009; 48:9256-65. [PMID: 19722680 PMCID: PMC2754664 DOI: 10.1021/bi901219x] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
The amide hydrogens that are exposed to solvent in the high-resolution X-ray structures of ubiquitin, FK506-binding protein, chymotrypsin inhibitor 2, and rubredoxin span a billion-fold range in hydroxide-catalyzed exchange rates which are predictable by continuum dielectric methods. To facilitate analysis of transiently accessible amides, the hydroxide-catalyzed rate constants for every backbone amide of ubiquitin were determined under near physiological conditions. With the previously reported NMR-restrained molecular dynamics ensembles of ubiquitin (PDB codes 2NR2 and 2K39) used as representations of the Boltzmann-weighted conformational distribution, nearly all of the exchange rates for the highly exposed amides were more accurately predicted than by use of the high-resolution X-ray structure. More strikingly, predictions for the amide hydrogens of the NMR relaxation-restrained ensemble that become exposed to solvent in more than one but less than half of the 144 protein conformations in this ensemble were almost as accurate. In marked contrast, the exchange rates for many of the analogous amides in the residual dipolar coupling-restrained ubiquitin ensemble are substantially overestimated, as was particularly evident for the Ile 44 to Lys 48 segment which constitutes the primary interaction site for the proteasome targeting enzymes involved in polyubiquitylation. For both ensembles, “excited state” conformers in this active site region having markedly elevated peptide acidities are represented at a population level that is 102 to 103 above what can exist in the Boltzmann distribution of protein conformations. These results indicate how a chemically consistent interpretation of amide hydrogen exchange can provide insight into both the population and the detailed structure of transient protein conformations.
Collapse
Affiliation(s)
- David M LeMaster
- Wadsworth Center, New York State Department of Health, School of Public Health, University at Albany-SUNY, Empire State Plaza, Albany, New York 12201, USA
| | | | | |
Collapse
|