Thermodynamics and folding kinetics analysis of the SH3 domain form discrete molecular dynamics

We perform a detailed analysis of the thermodynamics and folding kinetics of the SH3 domain fold with discrete molecular dynamic simulations. We propose a protein model that reproduces some of the experimentally observed thermodynamic and folding kinetic properties of proteins. Specifically, we use our model to study the transition state ensemble of the SH3 fold family of proteins, a set of unstable conformations that fold to the protein native state with probability 1/2. We analyze the participation of each secondary structure element formed at the transition state ensemble. We also identify the folding nucleus of the SH3 fold and test extensively its importance for folding kinetics. We predict that a set of amino acid contacts between the RT-loop and the distal hairpin are the critical folding nucleus of the SH3 fold and propose a hypothesis that explains this result.     

The effects of nonnative interactions on protein folding rates: theory and simulation

Proteins are minimally frustrated polymers. However, for realistic protein models, nonnative interactions must be taken into account. In this paper, we analyze the effect of nonnative interactions on the folding rate and on the folding free energy barrier. We present an analytic theory to account for the modification on the free energy landscape upon introduction of nonnative contacts, added as a perturbation to the strong native interactions driving folding. Our theory predicts a rate-enhancement regime at fixed temperature, under the introduction of weak, nonnative interactions. We have thoroughly tested this theoretical prediction with simulations of a coarse-grained protein model, by using an off-lattice C(alpha)model of the src-SH3 domain. The strong agreement between results from simulations and theory confirm the nontrivial result that a relatively small amount of nonnative interaction energy can actually assist the folding to the native structure.     

Conserved patterns and interactions in the unfolding transition state across SH3 domain structural homologues

Proteins with similar structures are generally assumed to arise from similar sequences. However, there are more cases than not where this is not true. The dogma is that sequence determines structure; how, then, can very different sequences fold to the same structure? Here, we employ high temperature unfolding simulations to probe the pathways and specific interactions that direct the folding and unfolding of the SH3 domain. The SH3 metafold in the Dynameomics Database consists of 753 proteins with the same structure, but varied sequences and functions. To investigate the relationship between sequence and structure, we selected 17 targets from the SH3 metafold with high sequence variability. Six unfolding simulations were performed for each target, transition states were identified, revealing two general folding/unfolding pathways at the transition state. Transition states were also expressed as mathematical graphs of connected chemical nodes, and it was found that three positions within the structure, independent of sequence, were consistently more connected within the graph than any other nearby positions in the sequence. These positions represent a hub connecting different portions of the structure. Multiple sequence alignment and covariation analyses also revealed certain positions that were more conserved due to packing constraints and stabilizing long‐range contacts. This study demonstrates that members of the SH3 domain with different sequences can unfold through two main pathways, but certain characteristics are conserved regardless of the sequence or unfolding pathway. While sequence determines structure, we show that disparate sequences can provide similar interactions that influence folding and lead to similar structures.     

The folding pathway of spectrin R17 from experiment and simulation: using experimentally validated MD simulations to characterize States hinted at by experiment

We present an experimental and computational analysis of the folding pathway of the 17th domain of chicken brain alpha-spectrin, R17. Wild-type R17 folds in a two-state manner and the chevron plot (plot of the logarithm of the observed rate constant against concentration of urea) shows essentially linear folding and unfolding arms. A number of mutant proteins, however, show a change in slope of the unfolding arm at high concentration of denaturant, hinting at complexity in the folding landscape. Through a combination of mutational studies and high temperature molecular dynamics simulations we show that the folding of R17 can be described by a model with two sequential transition states separated by an intermediate species. The rate limiting transition state for folding in water has been characterized both through experimental Phi-value analysis and by simulation. In contrast, a detailed analysis of the transition state predicted to dominate under highly denaturing conditions is only possible by simulation.     

Design and structural thermodynamic studies of the chimeric protein derived from spectrin SH3 domain

A number of the chimeric constructs with spectrin SH3 domain were designed for structural and thermodynamic studies of protein self-assembly and protein-ligand interactions. SH3 domains, components of many regulatory proteins, operate through weak interactions with proline-rich regions of polypeptide chains. The recombinant construct (WT-CIIA) studied in this work was constructed by linking the peptide ligand PPPVPPYSAG to the domain C-terminus via a long 12-residue linker to increase the affinity of this ligand for the spectrin domain, thereby ensuring a stable positioning of the polyproline helix to the conserved ligand-binding site in orientation II, which is regarded as untypical of the interaction between this domain and oligopeptides. A comparison of fluorescence spectra of the initial domain and the recombinant protein WT-CIIA suggests that the ligand sticks to the conservative binding site. However, analysis of the equilibrium urea-induced unfolding has demonstrated that this is an unstable contact, which leads to a two-stage unfolding of the chimeric protein. The protein WT-CIIA was crystallized; a set of X-ray diffraction data with a resolution of 1.75 Å was recorded from individual crystals. A preliminary analysis of these diffraction data has demonstrated that the crystals belong to space group P32 with the following unit cell parameters: a = b = 36.39, c = 112.17 Å, a = β = 90.0, and γ = 120.0.     

The adapter protein Nck: Role of individual SH3 and SH2 binding modules for protein interactions in T lymphocytes

Nck is a ubiquitously expressed, primarily cytosolic adapter protein consisting of one SH2 domain and three SH3 domains. It links receptor and nonreceptor tyrosine kinases to actin cytoskeleton reorganizing proteins. In T lymphocytes, Nck is a crucial component of signaling pathways for T cell activation and effector function. It recruits actin remodeling proteins to T cell receptor (TCR)‐associated activation clusters and thereby initiates changes in cell polarity and morphology. Moreover, Nck is crucial for the TCR‐induced mobilization of secretory vesicles to the cytotoxic immunological synapse. To identify the interactome of Nck in human T cells, we performed a systematic screen for interaction partners in untreated or pervanadate‐treated cells. We used GST fusion proteins containing full length Nck, the combined SH3 domains or the individual SH3 and SH2 domains to precipitate putative Nck interactors from cellular lysates. Protein bands were excised from gels, processed by tryptic in‐gel digestion and analyzed by mass spectrometry. Using this approach, we confirmed previously established interactions (e.g., with Slp76, CD3ε, WASP, and WIPF1) and identified several novel putative Nck‐binding proteins. We subsequently verified the SH2 domain binding to the actin‐binding protein HIP55 and to FYB/ADAP, and the SH3‐mediated binding to the nuclear proteins SFPQ/NONO. Using laser scanning microscopy, we provide new evidence for a nuclear localization of Nck in human T cells. Our data highlight the fundamental role of Nck in the TCR‐to‐cytoskeleton crosstalk and point to yet unknown nuclear functions of Nck also in T lymphocytes.     

Autophagosome biogenesis in plants: Roles of SH3P2

The autophagy core machinery is essentially conserved in eukaryotic cells for autophagy regulation. However, the underlying mechanisms for autophagosome formation in plant cells remain elusive. We have recently demonstrated that SH3 domain-containing protein 2 (SH3P2), a BAR (Bin-Amphiphysin-Rvs) domain protein, functions as a novel regulator for autophagosome biogenesis in Arabidopsis thaliana. Using SH3P2 and its GFP fusion as probes, we have characterized the dynamics and structures of autophagosome formation in plant cells. The phagophore assembly site, marked by SH3P2, is identified as having a close connection with the ER. SH3P2 also binds to phosphatidylinositol 3-phosphate (PtdIns3P) and functions downstream of the phosphatidylinositol 3-kinase (PtdIns3K) complex. Thus, SH3P2 serves as a novel membrane-associated protein in regulating autophagosome formation in Arabidopsis thaliana.     

The L49F mutation in alpha erythroid spectrin induces local disorder in the tetramer association region: Fluorescence and molecular dynamics studies of free and bound alpha spectrin

The bundling of the N‐terminal, partial domain helix (Helix C′) of human erythroid α‐spectrin (αI) with the C‐terminal, partial domain helices (Helices A′ and B′) of erythroid β‐spectrin (βI) to give a spectrin pseudo structural domain (triple helical bundle A′B′C′) has long been recognized as a crucial step in forming functional spectrin tetramers in erythrocytes. We have used apparent polarity and Stern–Volmer quenching constants of Helix C′ of αI bound to Helices A′ and B′ of βI, along with previous NMR and EPR results, to propose a model for the triple helical bundle. This model was used as the input structure for molecular dynamics simulations for both wild type (WT) and αI mutant L49F. The simulation output structures show a stable helical bundle for WT, but not for L49F. In WT, four critical interactions were identified: two hydrophobic clusters and two salt bridges. However, in L49F, the region downstream of Helix C′ was unable to assume a helical conformation and one critical hydrophobic cluster was disrupted. Other molecular interactions critical to the WT helical bundle were also weakened in L49F, possibly leading to the lower tetramer levels observed in patients with this mutation‐induced blood disorder.     

Hierarchy of simulation models in predicting molecular recognition mechanisms from the binding energy landscapes: structural analysis of the peptide complexes with SH2 domains.

Computer simulations using the simplified energy function and simulated tempering dynamics have accurately determined the native structure of the pYVPML, SVLpYTAVQPNE, and SPGEpYVNIEF peptides in the complexes with SH2 domains. Structural and equilibrium aspects of the peptide binding with SH2 domains have been studied by generating temperature-dependent binding free energy landscapes. Once some native peptide-SH2 domain contacts are constrained, the underlying binding free energy profile has the funnel-like shape that leads to a rapid and consistent acquisition of the native structure. The dominant native topology of the peptide-SH2 domain complexes represents an extended peptide conformation with strong specific interactions in the phosphotyrosine pocket and hydrophobic interactions of the peptide residues C-terminal to the pTyr group. The topological features of the peptide-protein interface are primarily determined by the thermodynamically stable phosphotyrosyl group. A diversity of structurally different binding orientations has been observed for the amino-terminal residues to the phosphotyrosine. The dominant native topology for the peptide residues carboxy-terminal to the phosphotyrosine is tolerant to flexibility in this region of the peptide-SH2 domain interface observed in equilibrium simulations. The energy landscape analysis has revealed a broad, entropically favorable topology of the native binding mode for the bound peptides, which is robust to structural perturbations. This could provide an additional positive mechanism underlying tolerance of the SH2 domains to hydrophobic conservative substitutions in the peptide specificity region.     

The origins of asymmetry in the folding transition states of protein L and protein G

Topology has been shown to be an important determinant of many features of protein folding; however, the delineation of sequence effects on folding remains obscure. Furthermore, differentiation between the two influences proves difficult due to their intimate relationship. To investigate the effect of sequence in the absence of significant topological differences, we examined the folding mechanisms of segment B1 peptostreptococcal protein L and segment B1 of streptococcal protein G. These proteins share the same highly symmetrical topology. Despite this symmetry, neither protein folds through a symmetrical transition state. We analyzed the origins of this difference using theoretical models. We found that the strength of the interactions present in the N-terminal hairpin of protein L causes this hairpin to form ahead of the C-terminal hairpin. The difference in chain entropy associated with the formation of the hairpins of protein G proves sufficient to beget initiation of folding at the shorter C-terminal hairpin. Our findings suggest that the mechanism of folding may be understood by examination of the free energy associated with the formation of partially folded microstates.     

Multicolor BiFC analysis of competition among G protein β and γ subunit interactions

We have applied multicolor BiFC to study the association preferences of G protein β and γ subunits in living cells. Cells co-express multiple isoforms of β and γ subunits, most of which can form complexes. Although many βγ complexes exhibit similar properties when assayed in reconstituted systems, knockout experiments in vivo suggest that individual isoforms have unique functions. BiFC makes it possible to correlate βγ complex formation with functionality in intact cells by comparing the amounts of fluorescent βγ complexes with their abilities to modulate effector proteins. The relative predominance of specific βγ complexes in vivo is not known. To address this issue, multicolor BiFC can determine the association preferences of β and γ subunits by simultaneously visualizing the two fluorescent complexes formed when β or γ subunits fused to amino terminal fragments of yellow fluorescent protein (YFP-N) and cyan fluorescent protein (CFP-N) compete to interact with limiting amounts of a common γ or β subunit, respectively, fused to a carboxyl terminal fragment of CFP (CFP-C). Multicolor BiFC also makes it possible to determine the roles of interacting proteins in the subcellular targeting of complexes, study the formation of protein complexes that are unstable under isolation conditions, determine the roles of co-expressed proteins in regulating the association preferences of interacting proteins, and visualize dynamic events affecting multiple protein complexes. These approaches can be applied to studying the assembly and functions of a wide variety of protein complexes in the context of a living cell.     

Artificial proteins as allosteric modulators of PDZ3 and SH3 in two‐domain constructs: A computational characterization of novel chimeric proteins

Artificial multidomain proteins with enhanced structural and functional properties can be utilized in a broad spectrum of applications. The design of chimeric fusion proteins utilizing protein domains or one‐domain miniproteins as building blocks is an important advancement for the creation of new biomolecules for biotechnology and medical applications. However, computational studies to describe in detail the dynamics and geometry properties of two‐domain constructs made from structurally and functionally different proteins are lacking. Here, we tested an in silico design strategy using all‐atom explicit solvent molecular dynamics simulations. The well‐characterized PDZ3 and SH3 domains of human zonula occludens (ZO‐1) (3TSZ), along with 5 artificial domains and 2 types of molecular linkers, were selected to construct chimeric two‐domain molecules. The influence of the artificial domains on the structure and dynamics of the PDZ3 and SH3 domains was determined using a range of analyses. We conclude that the artificial domains can function as allosteric modulators of the PDZ3 and SH3 domains. Proteins 2016; 84:1358–1374. © 2016 Wiley Periodicals, Inc.     

NMR structure and regulated expression in APL cell of human SH3BGRL3

SH3 domain binding glutamic acid-rich protein like 3 (SH3BGRL3) is the new member of thioredoxin (TRX) super family, whose posttranslational modified form was identified as tumor necrosis factor alpha (TNF-alpha) inhibitory protein, TIP-B1. In this paper, we determined its solution structure by multi-dimensional nuclear magnetic resonance spectroscopy. The overall structure of human SH3BGRL3 conformed to a TRX-like fold. To understand its function in vivo, the upregulated expression in acute promyelocytic leukemia cell line NB4 at both mRNA and protein level was elucidated. Immunofluorescence and immunohistochemistry staining with monoclonal antibody against SH3BGRL3 demonstrated that it was a cytoplasmic protein in both NB4 cell and human tissues. These results, as a whole, indicate that SH3BGRL3 may function as a regulator in all-trans retinoic acid-induced pathway.     

Structural investigation of CDCA3‐Cdh1 protein–protein interactions using in vitro studies and molecular dynamics simulation

The anaphase‐promoting complex/cyclosome (APC/C) ubiquitin ligase and its cofactor, Cdh1, regulate the expression of several cell‐cycle proteins and their functions during mitosis. Levels of the protein cell division cycle‐associated protein 3 (CDCA3), which is functionally required for mitotic entry, are regulated by APC/C^Cdh1. CDCA3 is an intrinsically disordered protein and contains both C‐terminal KEN box and D‐box recognition motifs, enabling binding to Cdh1. Our previous findings demonstrate that CDCA3 has a phosphorylation‐dependent non‐canonical ABBA‐like motif within the linker region bridging these two recognition motifs and is required for efficient binding to Cdh1. Here, we sought to identify and further characterize additional residues that participate within this ABBA‐like motif using detailed in vitro experiments and in silico modeling studies. We identified the role of H‐bonds, hydrophobic and ionic interactions across the CDCA3 ABBA‐like motif in the linker region between KEN and D‐box motifs. This linker region adopts a well‐defined structure when bound to Cdh1 in the presence of phosphorylation. Upon alanine mutation, the structure of this region is lost, leading to higher flexibility, and alteration in affinities due to binding to alternate sites on Cdh1. Our findings identify roles for the anchoring residues in the non‐canonical ABBA‐like motif to promote binding to the APC/C^Cdh1 and regulation of CDCA3 protein levels.     

Structure of the tandem MA-3 region of Pdcd4 protein and characterization of its interactions with eIF4A and eIF4G: molecular mechanisms of a tumor suppressor

One of the key regulatory points of translation initiation is recruitment of the 43S preinitation complex to the 5' mRNA cap by the eIF4F complex (eIF4A, eIF4E, and eIF4G). The tumor suppressor protein Pdcd4 has been shown to inhibit cap-dependent translation by interacting tightly with the RNA helicase eIF4A via its tandem MA-3 domains. The NMR studies reported here reveal a fairly extensive and well defined interface between the two MA-3 domains in solution, which appears to be stabilized by a network of interdomain salt bridges and hydrogen bonds, and reveals a unique orientation of the two domains. Characterization of the stoichiometry of the Pdcd4-eIF4A complex suggests that under physiological conditions Pdcd4 binds to a single molecule of eIF4A, which involves contacts with both Pdcd4 MA-3 domains. We also show that contacts mediated by a conserved acidic patch on the middle MA-3 domain of Pdcd4 are essential for forming a tight complex with eIF4A in vivo, whereas the equivalent region of the C-terminal MA-3 domain appears to have no role in complex formation in vivo. The formation of a 1:1 eIF4A-Pdcd4 complex in solution is consistent with the reported presence in vivo of only one molecule of eIF4A in the eIF4F complex. Pdcd4 has also been reported to interact directly with the middle region of eIF4G, however, we were unable to obtain any evidence for even a weak, transient direct interaction.     

Roles of G beta gamma in membrane recruitment and activation of p110 gamma/p101 phosphoinositide 3-kinase gamma

Receptor-regulated class I phosphoinositide 3-kinases (PI3K) phosphorylate the membrane lipid phosphatidylinositol (PtdIns)-4,5-P2 to PtdIns-3,4,5-P3. This, in turn, recruits and activates cytosolic effectors with PtdIns-3,4,5-P3-binding pleckstrin homology (PH) domains, thereby controlling important cellular functions such as proliferation, survival, or chemotaxis. The class IB p110 gamma/p101 PI3K gamma is activated by G beta gamma on stimulation of G protein-coupled receptors. It is currently unknown whether in living cells G beta gamma acts as a membrane anchor or an allosteric activator of PI3K gamma, and which role its noncatalytic p101 subunit plays in its activation by G beta gamma. Using GFP-tagged PI3K gamma subunits expressed in HEK cells, we show that G beta gamma recruits the enzyme from the cytosol to the membrane by interaction with its p101 subunit. Accordingly, p101 was found to be required for G protein-mediated activation of PI3K gamma in living cells, as assessed by use of GFP-tagged PtdIns-3,4,5-P3-binding PH domains. Furthermore, membrane-targeted p110 gamma displayed basal enzymatic activity, but was further stimulated by G beta gamma, even in the absence of p101. Therefore, we conclude that in vivo, G beta gamma activates PI3K gamma by a mechanism assigning specific roles for both PI3K gamma subunits, i.e., membrane recruitment is mediated via the noncatalytic p101 subunit, and direct stimulation of G beta gamma with p110 gamma contributes to activation of PI3K gamma.     

Structural insights into the interactions of phorbol ester and bryostatin complexed with protein kinase C: a comparative molecular dynamics simulation study

Protein kinase C (PKC) isozymes are important regulatory enzymes that have been implicated in many diseases, including cancer, Alzheimer’s disease, and in the eradication of HIV/AIDS. Given their potential clinical ramifications, PKC modulators, e.g. phorbol esters and bryostatin, are also of great interest in the drug development. However, structural details on the binding between PKC and its modulators, especially bryostatin – the highly potent and non-tumor promoting activator for PKCs, are still lacking. Here, we report the first comparative molecular dynamics study aimed at gaining structural insight into the mechanisms by which the PKC delta cys2 activator domain is used in its binding to phorbol ester and bryostatin-1. As anticipated in the phorbol ester binding, hydrogen bonds are formed through the backbone atoms of Thr242, Leu251, and Gly253 of PKC. However, the opposition of H-bond formation between Thr242 and Gly253 may cause the phorbol ester complex to become less stable when compared with the bryostatin binding. For the PKC delta-bryostatin complex, hydrogen bonds are formed between the Gly253 backbone carbonyl and the C30 carbomethoxy substituent of the ligand. Additionally, the indole Nε1 of the highly homologous Trp252 also forms an H-bond to the C20 ester group on bryostatin. Backbone fluctuations also suggest that this latter H-bond formation may abrogate the transient interaction between Trp252 and His269, thus dampening the fluctuations observed on the nearby Zn²⁺-coordinating residues. This new dynamic fluctuation dampening model can potentially benefit future design of new PKC modulators.     

Characterization of the overall and local dynamics of a protein with intermediate rotational anisotropy: Differentiating between conformational exchange and anisotropic diffusion in the B3 domain of protein G

Because the overall tumbling provides a major contribution to protein spectral densities measured in solution, the choice of a proper model for this motion is critical for accurate analysis of protein dynamics. Here we study the overall and backbone dynamics of the B3 domain of protein G using ¹⁵N relaxation measurements and show that the picture of local motions is markedly dependent on the model of overall tumbling. The main difference is in the interpretation of the elevated R ₂ values in the -helix: the isotropic model results in conformational exchange throughout the entire helix, whereas no exchange is predicted by anisotropic models that place the longitudinal axis of diffusion tensor almost parallel to the helix axis. Due to small size (fast tumbling) of the protein, the T ₁ values have low sensitivity to NH bond orientation. The diffusion tensor derived from orientation dependence of R ₂/R ₁ is anisotropic (D _par/D _perp=1.4), with a small rhombic component. In order to distinguish the correct picture of motion, we apply model-independent methods that are sensitive to conformational exchange and do not require knowledge of protein structure or assumptions about its dynamics. A comparison of the CSA/dipolar cross-correlation rate constants with ¹⁵N relaxation rates and the estimation of R _ex terms from relaxation data at 9.4 and 14.1 T indicate no conformational exchange in the helix, in support of the anisotropic models. The experimentally derived diffusion tensor is in excellent agreement with theoretical predictions from hydrodynamic calculations; a detailed comparison with various hydrodynamic models revealed optimal parameters for hydrodynamic calculations.     

Computational analysis of molecular basis of 1:1 interactions of NRG-1beta wild-type and variants with ErbB3 and ErbB4

The neuregulin/ErbB system is a growth factor/receptor cascade that has been proven to be essential in the development of the heart and the sympathetic nervous system. However, the basis of the specificity of ligand-receptor recognition remains to be elucidated. In this study, the structures of NRG-1beta/ErbB3 and NRG-1beta/ErbB4 complexes were modeled based on the available structures of the homologous proteins. The binding free energies of NRG-1beta to ErbB3 and ErbB4 were calculated using the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) computational method. In addition, computational alanine-scanning mutagenesis was performed in the binding site of NRG-1beta and the difference in the binding free energies between NRG-1beta mutants and the receptors was calculated. The results specify the contribution of each residue at the interaction interfaces to the binding affinity of NRG-1beta with ErbB3 and ErbB4, identifying several important interaction residue pairs that are in agreement with previously acquired experimental data. This indicates that the presented structural models of NRG-1beta/ErbB3 and NRG-1beta/ErbB4 complexes are reliable and could be used to guide future studies, such as performing desirable mutations on NRG-1beta to increase the binding affinity and selectivity to the receptor and discovering new therapeutic agents for the treatment of heart failure.