skip to main content

Title: A threonine zipper that mediates protein–protein interactions: Structure and prediction

We present the structure of an engineered protein–protein interface between two beta barrel proteins, which is mediated by interactions between threonine (Thr) residues. This Thr zipper structure suggests that the protein interface is stabilized by close‐packing of the Thr residues, with only one intermonomer hydrogen bond (H‐bond) between two of the Thr residues. This Thr‐rich interface provides a unique opportunity to study the behavior of Thr in the context of many other Thr residues. In previous work, we have shown that the side chain (χ1) dihedral angles of interface and core Thr residues can be predicted with high accuracy using a hard sphere plus stereochemical constraint (HS) model. Here, we demonstrate that in the Thr‐rich local environment of the Thr zipper structure, we are able to predict theχ1dihedral angles of most of the Thr residues. Some, however, are not well predicted by the HS model. We therefore employed explicitly solvated molecular dynamics (MD) simulations to further investigate the side chain conformations of these residues. The MD simulations illustrate the role that transient H‐bonding to water, in combination with steric constraints, plays in determining the behavior of these Thr side chains.

Broader Audience Statement: Protein–protein interactions are critical to life and the search for ways to disrupt adverse protein–protein interactions involved in disease is an ongoing area of drug discovery. We must better understand protein–protein interfaces, both to be able to disrupt existing ones and to engineer new ones for a variety of biotechnological applications. We have discovered and characterized an artificial Thr‐rich protein–protein interface. This novel interface demonstrates a heretofore unknown property of Thr‐rich surfaces: mediating protein–protein interactions.

more » « less
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
Protein Science
Page Range / eLocation ID:
p. 1969-1977
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Site-directed spin labeling (SDSL) ESR is a valuable tool to probe protein systems that are not amenable to characterization by x-ray crystallography, NMR or EM. While general principles that govern the shape of SDSL ESR spectra are known, its precise relationship with protein structure and dynamics is still not fully understood. To address this problem, we designed seven variants of GB1 domain bearing R1 spin label and recorded the corresponding MD trajectories (combined length 180 μs). The MD data were subsequently used to calculate time evolution of the relevant spin density matrix and thus predict the ESR spectra. The simulated spectra proved to be in good agreement with the experiment. Further analysis confirmed that the spectral shape primarily reflects the degree of steric confinement of the R1 tag and, for the well-folded protein such as GB1, offers little information on local backbone dynamics. The rotameric preferences of R1 side chain are determined by the type of the secondary structure at the attachment site. The rotameric jumps involving dihedral angles χ1and χ2are sufficiently fast to directly influence the ESR lineshapes. However, the jumps involving multiple dihedral angles tend to occur in (anti)correlated manner, causing smaller-than-expected movements of the R1 proxyl ring. Of interest, ESR spectra of GB1 domain with solvent-exposed spin label can be accurately reproduced by means of Redfield theory. In particular, the asymmetric character of the spectra is attributable to Redfield-type cross-correlations. We envisage that the current MD-based, experimentally validated approach should lead to a more definitive, accurate picture of SDSL ESR experiments.

    more » « less
  2. Abstract

    There have been several studies suggesting that protein structures solved by NMR spectroscopy and X‐ray crystallography show significant differences. To understand the origin of these differences, we assembled a database of high‐quality protein structures solved by both methods. We also find significant differences between NMR and crystal structures—in the root‐mean‐square deviations of the Cαatomic positions, identities of core amino acids, backbone, and side‐chain dihedral angles, and packing fraction of core residues. In contrast to prior studies, we identify the physical basis for these differences by modeling protein cores as jammed packings of amino acid‐shaped particles. We find that we can tune the jammed packing fraction by varying the degree of thermalization used to generate the packings. For an athermal protocol, we find that the average jammed packing fraction is identical to that observed in the cores of protein structures solved by X‐ray crystallography. In contrast, highly thermalized packing‐generation protocols yield jammed packing fractions that are even higher than those observed in NMR structures. These results indicate that thermalized systems can pack more densely than athermal systems, which suggests a physical basis for the structural differences between protein structures solved by NMR and X‐ray crystallography.

    more » « less
  3. Knowles, David A ; Mostafavi, Sara (Ed.)
    Accurately modeling protein 3D structure is essential for the design of functional proteins. An important sub-task of structure modeling is protein side-chain packing: predicting the conformation of side-chains (rotamers) given the protein’s backbone structure and amino-acid sequence. Conventional approaches for this task rely on expensive sampling procedures over hand-crafted energy functions and rotamer libraries. Recently, several deep learning methods have been developed to tackle the problem in a data-driven way, albeit with vastly different formulations (from image-to-image translation to directly predicting atomic coordinates). Here, we frame the problem as a joint regression over the side-chains’ true degrees of freedom: the dihedral $\chi$ angles. We carefully study possible objective functions for this task, while accounting for the underlying symmetries of the task. We propose Holographic Packer (H-Packer), a novel two-stage algorithm for side-chain packing built on top of two light-weight rotationally equivariant neural networks. We evaluate our method on CASP13 and CASP14 targets. H-Packer is computationally efficient and shows favorable performance against conventional physics-based algorithms and is competitive against alternative deep learning solutions. 
    more » « less
  4. Abstract

    Hydroxyl radical protein footprinting (HRPF) in combination with mass spectrometry reveals the relative solvent exposure of labeled residues within a protein, thereby providing insight into protein tertiary structure. HRPF labels nineteen residues with varying degrees of reliability and reactivity. Here, we are presenting a dynamics-driven HRPF-guided algorithm for protein structure prediction. In a benchmark test of our algorithm, usage of the dynamics data in a score term resulted in notable improvement of the root-mean-square deviations of the lowest-scoring ab initio models and improved the funnel-like metric Pnearfor all benchmark proteins. We identified models with accurate atomic detail for three of the four benchmark proteins. This work suggests that HRPF data along with side chain dynamics sampled by a Rosetta mover ensemble can be used to accurately predict protein structure.

    more » « less
  5. Abstract

    Regulator of G protein signaling (RGS) proteins play a pivotal role in regulation of G protein‐coupled receptor (GPCR) signaling and are therefore becoming an increasingly important therapeutic target. Recently discovered thiadiazolidinone (TDZD) compounds that target cysteine residues have shown different levels of specificities and potencies for the RGS4 protein, thereby suggesting intrinsic differences in dynamics of this protein upon binding of these compounds. In this work, we investigated using atomistic molecular dynamics (MD) simulations the effect of binding of several small‐molecule inhibitors on perturbations and dynamical motions in RGS4. Specifically, we studied two conformational models of RGS4 in which a buried cysteine residue is solvent‐exposed due to side‐chain motions or due to flexibility in neighboring helices. We found that TDZD compounds with aromatic functional groups perturb the RGS4 structure more than compounds with aliphatic functional groups. Moreover, small‐molecules with aromatic functional groups but lacking sulfur atoms only transiently reside within the protein and spontaneously dissociate to the solvent. We further measured inhibitory effects of TDZD compounds using a protein–protein interaction assay on a single‐cysteine RGS4 protein showing trends in potencies of compounds consistent with our simulation studies. Thermodynamic analyses of RGS4 conformations in the apo‐state and on binding to TDZD compounds revealed links between both conformational models of RGS4. The exposure of cysteine side‐chains appears to facilitate initial binding of TDZD compounds followed by migration of the compound into a bundle of four helices, thereby causing allosteric perturbations in the RGS/Gα protein–protein interface.

    more » « less