Functional profiling of asymmetrically-organized human CCT/TRiC chaperonin.

Molecular organization of the eukaryote chaperonin known as CCT/TRiC complex was recently clarified. Eight distinct subunits are uniquely organized, providing a favorable folding cavity for specific client proteins such as tubulin and actin. Because of its heterogeneous subunit composition, CCT complex has polarized inner faces, which may underlie an essential part of its chaperonin function. In this study, we structurally characterized the closed and open states of CCT complex, using molecular dynamics analyses. Our results showed that the inter-subunit interaction energies were asymmetrically distributed and were remodeled during conformational changes of CCT complex. In addition, exploration of redox related characteristics indicated changes in inner surface properties, including electrostatic potential, pKa and exposure of inner cysteine thiol groups, between the closed and open states. Cysteine activation events were experimentally verified by interaction analyses, using tubulin as a model substrate. Our data highlighted the importance of dynamics-based structural profiling of asymmetrically oriented chaperonin function.

All-atom molecular dynamics analysis of multi-peptide systems reproduces peptide solubility in line with experimental observations.

In order to investigate the contribution of individual amino acids to protein and peptide solubility, we carried out 100 ns molecular dynamics (MD) simulations of 10(6) Å(3) cubic boxes containing ~3 × 10(4) water molecules and 27 tetra-peptides regularly positioned at 23 Å from each other and composed of a single amino acid type for all natural amino acids but cysteine and glycine. The calculations were performed using Amber with a standard force field on a special purpose MDGRAPE-3 computer, without introducing any "artificial" hydrophobic interactions. Tetra-peptides composed of I, V, L, M, N, Q, F, W, Y, and H formed large amorphous clusters, and those containing A, P, S, and T formed smaller ones. Tetra-peptides made of D, E, K, and R did not cluster at all. These observations correlated well with experimental solubility tendencies as well as hydrophobicity scales with correlation coefficients of 0.5 to > 0.9. Repulsive Coulomb interactions were dominant in ensuring high solubility, whereas both Coulomb and van der Waals (vdW) energies contributed to the aggregations of low solubility amino acids. Overall, this very first all-atom molecular dynamics simulation of a multi-peptide system appears to reproduce the basic properties of peptide solubility, essentially in line with experimental observations.

Calculating the Na⁺ translocating V-ATPase catalytic site affinity for substrate binding by homology modeled NtpA monomer using molecular dynamics/free energy calculation.

Vacuolar ATPase (V-ATPase) of Enterococcus hirae is composed of a soluble catalytic domain (V1; NtpA3-B3-D-G) and an integral membrane domain (V0; NtpI-K10) connected by a central and two peripheral stalks (NtpC, NtpD-G and NtpE-F). Recently nucleotide binding of catalytic NtpA monomer has been reported (Arai et al.). In the present study, we calculated the nucleotide binding affinity of NtpA by molecular dynamics (MD) simulation/free energy calculation using MM-GBSA approach based on homology modeled structure of NtpA monomer docked with ATP analogue, adenosine 5'-[ß, γ-imido] triphosphate (AMP-PNP). The calculated binding free energies showed qualitatively good agreement with experimental data. The calculation was cross-validated further by the rigorous method, thermodynamic integration (TI) simulation. Finally, the interaction between NtpA and nucleotides at the atomic level was investigated by the analyses of components of free energy and the optimized model structures obtained from MD simulations, suggesting that electrostatic contribution is responsible for the difference in nucleotide binding to NtpA monomer. This is the first observation and suggestion to explain the difference of nucleotide binding properties in V-ATPase NtpA subunit, and our method can be a valuable primary step to predict nucleotide binding affinity to other subunits (NtpAB, NtpA3B3) and to explore subunit interactions and eventually may help to understand energy transduction mechanism of E. hirae V-ATPase.

Blockade of inflammatory responses by a small-molecule inhibitor of the Rac activator DOCK2.

Tissue infiltration of activated lymphocytes is a hallmark of transplant rejection and organ-specific autoimmune diseases. Migration and activation of lymphocytes depend on DOCK2, an atypical Rac activator predominantly expressed in hematopoietic cells. Although DOCK2 does not contain Dbl homology domain typically found in guanine nucleotide exchange factors, DOCK2 mediates the GTP-GDP exchange reaction for Rac through its DHR-2 domain. Here, we have identified 4-[3'-(2″-chlorophenyl)-2'-propen-1'-ylidene]-1-phenyl-3,5-pyrazolidinedione (CPYPP) as a small-molecule inhibitor of DOCK2. CPYPP bound to DOCK2 DHR-2 domain in a reversible manner and inhibited its catalytic activity in vitro. When lymphocytes were treated with CPYPP, both chemokine receptor- and antigen receptor-mediated Rac activation were blocked, resulting in marked reduction of chemotactic response and T cell activation. These results provide a rational of and a chemical scaffold for development of the DOCK2-targeting immunosuppressant.

Mutation of epidermal growth factor receptor is associated with MIG6 expression.

Controlled activation of epidermal growth factor receptor (EGFR) is systematically guaranteed at the molecular level; however, aberrant activation of EGFR is frequently found in cancer. Transcription induced by EGFR activation often involves the coordinated expression of genes that positively and negatively regulate the original signaling pathway; therefore, alterations in EGFR kinase activity may reflect changes in gene expression associated with the pathway. In the present study, we investigated transcriptional changes after EGF stimulation with or without the EGFR kinase inhibitor Iressa in H1299 human non-small-cell lung cancer cells [parental H1299, H1299 cells that overexpress wild-type EGFR (EGFR-WT) and mutant H1299 cells that overexpress EGFR where Leu858 is substituted with Arg (L858R)]. The results obtained clearly demonstrate differences in transcriptional activity in the absence or presence of EGFR kinase activity, with genes sharing the same molecular functions showing distinct expression dynamics. The results show the particular enrichment of EGFR/ErbB signaling-related genes in a differentially expressed gene set, and significant protein expression of MIG6/RALT(ERRFI1), an EGFR negative regulator, was confirmed in L858R. High MIG6 protein expression was correlated with basal EGFR phosphorylation and inversely correlated with EGF-induced extracellular signal-regulated protein kinase phosphorylation levels. Investigation of the NCI-60 cell lines showed that ERRFI1 expression was correlated with EGFR expression, regardless of tissue type. These results suggest that cells accumulate MIG6 as an inherent negative regulator to suppress excess EGFR activity when basal EGFR kinase activity is considerably high. Taking all the above together, an EGFR mutation can cause transcriptional changes to accommodate the activation potency of the original signaling pathway at the cellular level.

Nanoscale elongating control of the self-assembled protein filament with the cysteine-introduced building blocks.

Self-assembly of artificially designed proteins is extremely desirable for nanomaterials. Here we show a novel strategy for the creation of self-assembling proteins, named "Nanolego." Nanolego consists of "structural elements" of a structurally stable symmetrical homo-oligomeric protein and "binding elements," which are multiple heterointeraction proteins with relatively weak affinity. We have established two key technologies for Nanolego, a stabilization method and a method for terminating the self-assembly process. The stabilization method is mediated by disulfide bonds between Cysteine-residues incorporated into the binding elements, and the termination method uses "capping Nanolegos," in which some of the binding elements in the Nanolego are absent for the self-assembled ends. With these technologies, we successfully constructed timing-controlled and size-regulated filament-shape complexes via Nanolego self-assembly. The Nanolego concept and these technologies should pave the way for regulated nanoarchitecture using designed proteins.

Molecular dynamics simulations reveal that Tyr-317 phosphorylation reduces Shc binding affinity for phosphotyrosyl residues of epidermal growth factor receptor.

The Src homology 2 (SH2) and collagen domain protein Shc plays a pivotal role in signaling via tyrosine kinase receptors, including epidermal growth factor receptor (EGFR). Shc binding to phospho-tyrosine residues on activated receptors is mediated by the SH2 and phospho-tyrosine binding (PTB) domains. Subsequent phosphorylation on Tyr-317 within the Shc linker region induces Shc interactions with Grb2-Son of Sevenless that initiate Ras-mitogen-activated protein kinase signaling. We use molecular dynamics simulations of full-length Shc to examine how Tyr-317 phosphorylation controls Shc conformation and interactions with EGFR. Our simulations reveal that Shc tyrosine phosphorylation results in a significant rearrangement of the relative position of its domains, suggesting a key conformational change. Importantly, computational estimations of binding affinities show that EGFR-derived phosphotyrosyl peptides bind with significantly more strength to unphosphorylated than to phosphorylated Shc. Our results unveil what we believe is a novel structural phenomenon, i.e., tyrosine phosphorylation of Shc within its linker region regulates the binding affinity of SH2 and PTB domains for phosphorylated Shc partners, with important implications for signaling dynamics.

Structural and functional basis of a role for CRKL in a fibroblast growth factor 8-induced feed-forward loop.

The adapter protein CRKL is required for the normal development of multiple tissues that rely on fibroblast growth factor 8 (FGF8). The precise role of CRKL in receptor signaling has been unclear, however. To address this issue, we first modeled the three-dimensional structure of CRKL by molecular dynamics. By taking advantage of structural simulations, we performed in silico analysis of the interactions of the autophosphorylation sites of FGR receptor 1 (FGFR1) with the SH2 domain of CRKL or a highly related protein, CRK. As predicted by simulations, we confirm the specific physical interaction of phosphorylated Y463 (pY463) in FGFR1 with the CRKL SH2 domain at an affinity approximately 30-fold stronger than that of CRK. We also provide evidence that interactions outside of the core YXXP motif have a significant impact on physical association, which is consistent with predictions from molecular-dynamics simulations. Furthermore, we identify CRKL as an essential component of an FGF8-induced feed-forward loop permissive for efficient activation of the mitogen-activated protein kinase Erk1/2, as well as FGF8-induced anchorage-independent cell growth, using Crkl-deficient cells or a pY463 synthetic peptide. Although many cells generally require cell-matrix adhesion, our results demonstrate that CRKL permits cells to bypass the strict need for adhesion in response to FGF8 through direct interaction with receptor.

Crystal structure of the rac activator, Asef, reveals its autoinhibitory mechanism.

The Rac-specific guanine nucleotide exchange factor (GEF) Asef is activated by binding to the tumor suppressor adenomatous polyposis coli mutant, which is found in sporadic and familial colorectal tumors. This activated Asef is involved in the migration of colorectal tumor cells. The GEFs for Rho family GTPases contain the Dbl homology (DH) domain and the pleckstrin homology (PH) domain. When Asef is in the resting state, the GEF activity of the DH-PH module is intramolecularly inhibited by an unidentified mechanism. Asef has a Src homology 3 (SH3) domain in addition to the DH-PH module. In the present study, the three-dimensional structure of Asef was solved in its autoinhibited state. The crystal structure revealed that the SH3 domain binds intramolecularly to the DH domain, thus blocking the Rac-binding site. Furthermore, the RT-loop and the C-terminal region of the SH3 domain interact with the DH domain in a manner completely different from those for the canonical binding to a polyproline-peptide motif. These results demonstrate that the blocking of the Rac-binding site by the SH3 domain is essential for Asef autoinhibition. This may be a common mechanism in other proteins that possess an SH3 domain adjacent to a DH-PH module.

Novel mechanism of interaction of p85 subunit of phosphatidylinositol 3-kinase and ErbB3 receptor-derived phosphotyrosyl peptides.

Ligand-activated and tyrosine-phosphorylated ErbB3 receptor binds to the SH2 domain of the p85 subunit of phosphatidylinositol 3-kinase and initiates intracellular signaling. Here, we studied the interactions between the N- (N-SH2) and C- (C-SH2) terminal SH2 domains of the p85 subunit of the phosphatidylinositol 3-kinase and eight ErbB3 receptor-derived phosphotyrosyl peptides (P-peptides) by using molecular dynamics, free energy, and surface plasmon resonance (SPR) analyses. In SPR analysis, these P-peptides showed no binding to the C-SH2 domain, but P-peptides containing a phospho-YXXM or a non-phospho-YXXM motif did bind to the N-SH2 domain. The N-SH2 domain has two phosphotyrosine binding sites in its N- (N1) and C- (N2) terminal regions. Interestingly, we found that P-peptides of pY1180 and pY1241 favored to bind to the N2 site, although all other P-peptides showed favorable binding to the N1 site. Remarkably, two phosphotyrosines, pY1178 and pY1243, which are just 63 amino acids apart from the pY1241 and pY1180, respectively, showed favorable binding to the N1 site. These findings indicate a possibility that the pair of phosphotyrosines, pY1178-pY1241 or pY1243-pY1180, will fold into an appropriate configuration for binding to the N1 and N2 sites simultaneously. Our model structures of the cytoplasmic C-terminal domain of ErbB3 receptor also strongly supported the speculation. The calculated binding free energies between the N-SH2 domain and P-peptides showed excellent qualitative agreement with SPR data with a correlation coefficient of 0.91. The total electrostatic solvation energy between the N-SH2 domain and P-peptide was the dominant factor for its binding affinity.

Increased rigidity of domain structures enhances the stability of a mutant enzyme created by directed evolution.

A mutant of kanamycin nucleotidyltransferase (KNT) was previously created by directed evolution. This mutant, HTK, has 19 amino acid substitutions, which increase the thermostability by 20 degrees C. In this study, we have examined to what extent each mutation contributes to the increased stability and analyzed how the mutations affect the structure of KNT at 72 degrees C using molecular dynamics simulations. The effects of some mutations on the stability are simply additive, but those of others are cooperative. Mutations with large effects on the stability are introduced into the N-terminal domain, which appears to be less stable than the C-terminal domain. Results of the molecular dynamics simulations have indicated that the rigidity of the domain structures is increased by the mutations: at 72 degrees C, the intradomain fluctuations of HTK are decreased, and in turn, its interdomain motions are pronounced, whereas the structure of the preevolved KNT fluctuates randomly. Chemical modification experiments of cysteine residues have shown that the cysteine residues of HTK are less accessible to an SH reagent than those of the preevolved KNT. The present results suggest that the 19 mutations of HTK stabilize KNT by affecting the dynamic behavior of the structure of this enzyme without significantly changing its static overall structure.

Molecular dynamics, free energy, and SPR analyses of the interactions between the SH2 domain of Grb2 and ErbB phosphotyrosyl peptides.

We studied the interactions between the SH2 domain of growth factor receptor binding protein 2 (Grb2) and ErbB receptor-derived phosphotyrosyl peptides using molecular dynamics, free energy calculations, and surface plasmon resonance (SPR) analysis. Binding free energies for nine phosphotyrosyl peptides were calculated using the MM-PBSA continuum solvent method, and excellent qualitative agreement with the SPR experimental data, with a correlation coefficient of 0.92, was obtained. Consistent with previous experimental findings, phosphotyrosyl peptides with the consensus sequence pYXNX showed favorable binding affinity for the Grb2. Unexpectedly, phosphotyrosyl peptides with the consensus sequence pYQQD, which had not shown any specific binding affinity for the Grb2 in earlier studies, also showed favorable binding affinity for the Grb2 in our experimental and computational analyses. Component analysis of the calculated binding free energies revealed that van der Waals interaction between the Grb2 and the phosphotyrosyl peptide was the dominant factor for specificity and binding affinity. These results indicate that current methods of estimating binding free energies are efficient for obtaining important information about protein-protein interactions, which are essential for the transmission of signals in cellular signaling pathways.