Ligand binding remodels protein side chain conformational heterogeneity

  1. Stephanie A Wankowicz
  2. Saulo H de Oliveira
  3. Daniel W Hogan
  4. Henry van den Bedem
  5. James S Fraser  Is a corresponding author
  1. University of California, San Francisco, United States
  2. Atomwise, Inc, United States

Abstract

While protein conformational heterogeneity plays an important role in many aspects of biological function, including ligand binding, its impact has been difficult to quantify. Macromolecular X-ray diffraction is commonly interpreted with a static structure, but it can provide information on both the anharmonic and harmonic contributions to conformational heterogeneity. Here, through multiconformer modeling of time- and space-averaged electron density, we measure conformational heterogeneity of 743 stringently matched pairs of crystallographic datasets that reflect unbound/apo and ligand-bound/holo states. When comparing the conformational heterogeneity of side chains, we observe that when binding site residues become more rigid upon ligand binding, distant residues tend to become more flexible, especially in non-solvent exposed regions. Among ligand properties, we observe increased protein flexibility as the number of hydrogen bonds decrease and relative hydrophobicity increases. Across a series of 13 inhibitor bound structures of CDK2, we find that conformational heterogeneity is correlated with inhibitor features and identify how conformational changes propagate differences in conformational heterogeneity away from the binding site. Collectively, our findings agree with models emerging from NMR studies suggesting that residual side chain entropy can modulate affinity and point to the need to integrate both static conformational changes and conformational heterogeneity in models of ligand binding.

Data availability

Refined models are available here: https://zenodo.org/record/5533006#.YVJr2Z5KgUsCode can be found in the following repositories:-Dataset selection: https://github.com/stephaniewankowicz/PDB_selection_pipeline-Refinement/qFit pipeline: https://github.com/stephaniewankowicz/refinement_qFit-Analysis/Figures: https://github.com/fraser-lab/Apo_Holo_Analysis-qFit: https://github.com/ExcitedStates/qfit-3.0.

The following data sets were generated

Article and author information

Author details

  1. Stephanie A Wankowicz

    Department of Bioengineering and Therapeutic Sciences, University of California, San Francisco, San Francisco, United States
    Competing interests
    No competing interests declared.
  2. Saulo H de Oliveira

    Atomwise, Inc, San Francisco, United States
    Competing interests
    Saulo H de Oliveira, is an employee of Atomwise Inc..
  3. Daniel W Hogan

    Department of Bioengineering and Therapeutic Sciences, University of California, San Francisco, San Francisco, United States
    Competing interests
    No competing interests declared.
  4. Henry van den Bedem

    Atomwise, Inc, San Francisco, United States
    Competing interests
    Henry van den Bedem, is an employee of Atomwise Inc..
  5. James S Fraser

    Department of Bioengineering and Therapeutic Sciences, University of California, San Francisco, San Francisco, United States
    For correspondence
    jfraser@fraserlab.com
    Competing interests
    James S Fraser, has equity, has received consulting fees, and has sponsored research agreements with Relay Therapeutics..
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-5080-2859

Funding

National Science Foundation (GRFP 2034836)

  • Stephanie A Wankowicz

National Institutes of Health (GM123159)

  • James S Fraser

National Institutes of Health (GM124149)

  • James S Fraser

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Copyright

© 2022, Wankowicz et al.

This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 4,060
    views
  • 582
    downloads
  • 49
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Stephanie A Wankowicz
  2. Saulo H de Oliveira
  3. Daniel W Hogan
  4. Henry van den Bedem
  5. James S Fraser
(2022)
Ligand binding remodels protein side chain conformational heterogeneity
eLife 11:e74114.
https://doi.org/10.7554/eLife.74114

Share this article

https://doi.org/10.7554/eLife.74114

Further reading

    1. Neuroscience
    2. Structural Biology and Molecular Biophysics
    Amy N Shore, Keyong Li ... Matthew C Weston
    Research Article

    More than 20 recurrent missense gain-of-function (GOF) mutations have been identified in the sodium-activated potassium (KNa) channel gene KCNT1 in patients with severe developmental and epileptic encephalopathies (DEEs), most of which are resistant to current therapies. Defining the neuron types most vulnerable to KCNT1 GOF will advance our understanding of disease mechanisms and provide refined targets for precision therapy efforts. Here, we assessed the effects of heterozygous expression of a Kcnt1 GOF variant (Kcnt1Y777H) on KNa currents and neuronal physiology among cortical glutamatergic and GABAergic neurons in mice, including those expressing vasoactive intestinal polypeptide (VIP), somatostatin (SST), and parvalbumin (PV), to identify and model the pathogenic mechanisms of autosomal dominant KCNT1 GOF variants in DEEs. Although the Kcnt1Y777H variant had no effects on glutamatergic or VIP neuron function, it increased subthreshold KNa currents in both SST and PV neurons but with opposite effects on neuronal output; SST neurons became hypoexcitable with a higher rheobase current and lower action potential (AP) firing frequency, whereas PV neurons became hyperexcitable with a lower rheobase current and higher AP firing frequency. Further neurophysiological and computational modeling experiments showed that the differential effects of the Kcnt1Y777H variant on SST and PV neurons are not likely due to inherent differences in these neuron types, but to an increased persistent sodium current in PV, but not SST, neurons. The Kcnt1Y777H variant also increased excitatory input onto, and chemical and electrical synaptic connectivity between, SST neurons. Together, these data suggest differential pathogenic mechanisms, both direct and compensatory, contribute to disease phenotypes, and provide a salient example of how a pathogenic ion channel variant can cause opposite functional effects in closely related neuron subtypes due to interactions with other ionic conductances.

    1. Biochemistry and Chemical Biology
    2. Structural Biology and Molecular Biophysics
    Marina Dajka, Tobias Rath ... Benesh Joseph
    Research Article

    Lipopolysaccharides (LPS) confer resistance against harsh conditions, including antibiotics, in Gram-negative bacteria. The lipopolysaccharide transport (Lpt) complex, consisting of seven proteins (A-G), exports LPS across the cellular envelope. LptB2FG forms an ATP-binding cassette transporter that transfers LPS to LptC. How LptB2FG couples ATP binding and hydrolysis with LPS transport to LptC remains unclear. We observed the conformational heterogeneity of LptB2FG and LptB2FGC in micelles and/or proteoliposomes using pulsed dipolar electron spin resonance spectroscopy. Additionally, we monitored LPS binding and release using laser-induced liquid bead ion desorption mass spectrometry. The β-jellyroll domain of LptF stably interacts with the LptG and LptC β-jellyrolls in both the apo and vanadate-trapped states. ATP binding at the cytoplasmic side is allosterically coupled to the selective opening of the periplasmic LptF β-jellyroll domain. In LptB2FG, ATP binding closes the nucleotide binding domains, causing a collapse of the first lateral gate as observed in structures. However, the second lateral gate, which forms the putative entry site for LPS, exhibits a heterogeneous conformation. LptC binding limits the flexibility of this gate to two conformations, likely representing the helix of LptC as either released from or inserted into the transmembrane domains. Our results reveal the regulation of the LPS entry gate through the dynamic behavior of the LptC transmembrane helix, while its β-jellyroll domain is anchored in the periplasm. This, combined with long-range ATP-dependent allosteric gating of the LptF β-jellyroll domain, may ensure efficient and unidirectional transport of LPS across the periplasm.