skip to main content


This content will become publicly available on August 13, 2024

Title: Engineering gain‐of‐function mutants of a WW domain by dynamics and structural analysis
Abstract

Proteins gain optimal fitness such as foldability and function through evolutionary selection. However, classical studies have found that evolutionarily designed protein sequences alone cannot guarantee foldability, or at least not without considering local contacts associated with the initial folding steps. We previously showed that foldability and function can be restored by removing frustration in the folding energy landscape of a model WW domain protein, CC16, which was designed based on Statistical Coupling Analysis (SCA). Substitutions ensuring the formation of five local contacts identified as “on‐path” were selected using the closest homolog native folded sequence, N21. Surprisingly, the resulting sequence, CC16‐N21, bound to Group I peptides, while N21 did not. Here, we identified single‐point mutations that enable N21 to bind a Group I peptide ligand through structure and dynamic‐based computational design. Comparison of the docked position of the CC16‐N21/ligand complex with the N21 structure showed that residues at positions 9 and 19 are important for peptide binding, whereas the dynamic profiles identified position 10 as allosterically coupled to the binding site and exhibiting different dynamics between N21 and CC16‐N21. We found that swapping these positions in N21 with matched residues from CC16‐N21 recovers nature‐like binding affinity to N21. This study validates the use of dynamic profiles as guiding principles for affecting the binding affinity of small proteins.

 
more » « less
Award ID(s):
1901709
NSF-PAR ID:
10451573
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
Protein Science
Volume:
32
Issue:
9
ISSN:
0961-8368
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    AlphaFold2 has revolutionized protein structure prediction from amino‐acid sequence. In addition to protein structures, high‐resolution dynamics information about various protein regions is important for understanding protein function. Although AlphaFold2 has neither been designed nor trained to predict protein dynamics, it is shown here how the information returned by AlphaFold2 can be used to predict dynamic protein regions at the individual residue level. The approach, which is termed cdsAF2, uses the 3D protein structure returned by AlphaFold2 to predict backbone NMR NHS2order parameters using a local contact model that takes into account the contacts made by each peptide plane along the backbone with its environment. By combining for each residue AlphaFold2's pLDDT confidence score for the structure prediction accuracy with the predictedS2value using the local contact model, an estimator is obtained that semi‐quantitatively captures many of the dynamics features observed in experimental backbone NMR NHS2order parameter profiles. The method is demonstrated for a set nine proteins of different sizes and variable amounts of dynamics and disorder.

     
    more » « less
  2. Abstract

    Human acyl protein thioesterases (APTs) catalyze the depalmitoylation ofS‐acylated proteins attached to the plasma membrane, facilitating reversible cycles of membrane anchoring and detachment. We previously showed that a bacterial APT homologue, FTT258 from the gram‐negative pathogenFrancisella tularensis, exists in equilibrium between a closed and open state based on the structural dynamics of a flexible loop overlapping its active site. Although the structural dynamics of this loop are not conserved in human APTs, the amino acid sequence of this loop is highly conserved, indicating essential but divergent functions for this loop in human APTs. Herein, we investigated the role of this loop in regulating the catalytic activity, ligand binding, and protein folding of human APT1, a depalmitoylase connected with cancer, immune, and neurological signaling. Using a combination of substitutional analysis with kinetic, structural, and biophysical characterization, we show that even in its divergent structural location in human APT1 that this loop still regulates the catalytic activity of APT1 through contributions to ligand binding and substrate positioning. We confirmed previously known roles for multiple residues (Phe72 and Ile74) in substrate binding and catalysis while adding new roles in substrate selectivity (Pro69), in catalytic stabilization (Asp73 and Ile75), and in transitioning between the membrane binding β‐tongue and substrate‐binding loops (Trp71). Even conservative substitution of this tryptophan (Trp71) fulcrum led to complete loss of catalytic activity, a 13°C decrease in total protein stability, and drastic drops in ligand affinity, indicating that the combination of the size, shape, and aromaticity of Trp71 are essential to the proper structure of APT1. Mixing buried hydrophobic surface area with contributions to an exposed secondary surface pocket, Trp71 represents a previously unidentified class of essential tryptophans within α/β hydrolase structure and a potential allosteric binding site within human APTs.

     
    more » « less
  3. Abstract

    Structure and functions of S100 proteins are regulated by two distinct calcium binding EF hand motifs. In this work, we used solution‐state NMR spectroscopy to investigate the cooperativity between the two calcium binding sites and map the allosteric changes at the target binding site. To parse the contribution of the individual calcium binding events, variants of S100A12 were designed to selectively bind calcium to either the EF‐I (N63A) or EF‐II (E31A) loop, respectively. Detailed analysis of the backbone chemical shifts for wildtype protein and its mutants indicates that calcium binding to the canonical EF‐II loop is the principal trigger for the conformational switch between ‘closed’ apo to the ‘open’ Ca2+‐bound conformation of the protein. Elimination of binding in S100‐specific EF‐I loop has limited impact on the calcium binding affinity of the EF‐II loop and the concomitant structural rearrangement. In contrast, deletion of binding in the EF‐II loop significantly attenuates calcium affinity in the EF‐I loop and the structure adopts a ‘closed’ apo‐like conformation. Analysis of experimental amide nitrogen (15N) relaxation rates (R1,R2, and15N–{1H} NOE) and molecular dynamics (MD) simulations demonstrate that the calcium bound state is relatively floppy with pico–nanosecond motions induced in functionally relevant domains responsible for target recognition such as the hinge domain and the C‐terminal residues. Experimental relaxation studies combined with MD simulations show that while calcium binding in the EF‐I loop alone does not induce significant motions in the polypeptide chain, EF‐I regulates fluctuations in the polypeptide in the presence of bound calcium in the EF‐II loop. These results offer novel insights into the dynamic regulation of target recognition by calcium binding and unravels the role of cooperativity between the two calcium binding events in S100A12.

     
    more » « less
  4. Abstract

    Charged residues on the surface of proteins are critical for both protein stability and interactions. However, many proteins contain binding regions with a high net charge that may destabilize the protein but are useful for binding to oppositely charged targets. We hypothesized that these domains would be marginally stable, as electrostatic repulsion would compete with favorable hydrophobic collapse during folding. Furthermore, by increasing the salt concentration, we predict that these protein folds would be stabilized by mimicking some of the favorable electrostatic interactions that take place during target binding. We varied the salt and urea concentrations to probe the contributions of electrostatic and hydrophobic interactions for the folding of the yeast SH3 domain found in Abp1p. The SH3 domain was significantly stabilized with increased salt concentrations due to Debye–Huckel screening and a nonspecific territorial ion‐binding effect. Molecular dynamics and NMR show that sodium ions interact with all 15 acidic residues but do little to change backbone dynamics or overall structure. Folding kinetics experiments show that the addition of urea or salt primarily affects the folding rate, indicating that almost all the hydrophobic collapse and electrostatic repulsion occur in the transition state. After the transition state formation, modest yet favorable short‐range salt bridges are formed along with hydrogen bonds, as the native state fully folds. Thus, hydrophobic collapse offsets electrostatic repulsion to ensure this highly charged binding domain can still fold and be ready to bind to its charged peptide targets, a property that is likely evolutionarily conserved over 1 billion years.

     
    more » « less
  5. Middle East Respiratory Syndrome Coronavirus (MERS-CoV) causes severe pneumonia-like symptoms and is still pose a significant threat to global public health. A key component in the virulence of MERS-CoV is the Spike (S) protein, which binds with the host membrane receptor dipeptidyl peptidase 4 (DPP4). The goal of the present investigation is to examine the effects of missense mutations in the MERS-CoV S protein on protein stability and binding affinity with DPP4 to provide insight that is useful in developing vaccines to prevent coronavirus infection. We utilized a saturation mutagenesis approach to simulate all possible mutations in the MERS-CoV full-length S, S Receptor Binding Domain (RBD) and DPP4. We found the mutations in MERS-CoV S protein residues, G552, C503, C526, N468, G570, S532, S451, S419, S465, and S435, affect protein stability. We identified key residues, G538, E513, V555, S557, L506, L507, R511, M452, D537, and S454 in the S protein RBD region are important in the binding of MERS-CoV S protein to the DPP4 receptor. We investigated the effects of MERS-CoV S protein viral mutations on protein stability and binding affinity. In addition, we studied all DPP4 mutations and found the functional substitution R336T weakens both DPP4 protein stability and S-DPP4 binding affinity. We compared the S protein structures of MERS-CoV, SARS-CoV, and SARS-CoV-2 viruses and identified the residues like C526, C383, and N468 located in equivalent positions of these viruses have effects on S protein structure. These findings provide further information on how mutations in coronavirus S proteins effect protein function. 
    more » « less