Hsc70 contacts helix III of the J domain from polyomavirus T antigens: addressing a dilemma in the chaperone hypothesis of how they release E2F from pRb

Hsc70's expected binding site on helix II of the J domain of T antigens appears to be blocked in its structure bound to tumor suppressor pRb. We used NMR to map where mammalian Hsc70 binds the J domain of murine polyomavirus T antigens (PyJ). The ATPase domain of Hsc70 unexpectedly has its biggest effects on the NMR peak positions of the C-terminal end of helix III of PyJ. The Hsc70 ATPase domain protects the C-terminal end of helix III of PyJ from an uncharged paramagnetic probe of chelated Gd(III), clearly suggesting the interface. Effects on the conserved HPD loop and helix II of PyJ are smaller. The NMR results are supported by a novel assay of Hsc70's ATP hydrolysis showing that mutations of surface residues in PyJ helix III impair PyJ-dependent stimulation of Hsc70 activity. Evolutionary trace analysis of J domains suggests that helix III usually may join helix II in contributing specificities for cognate hsp70s. Our novel evidence implicating helix III differs from evidence that Escherichia coli DnaK primarily affects helix II and the HPD loop of DnaJ. We find the pRb-binding fragment of E2F1 to be intrinsically unfolded and a good substrate for Hsc70 in vitro. This suggests that E2F1 could be a substrate for Hsc70 recruited by T antigen to an Rb family member. Importantly, our results strengthen the chaperone hypothesis for E2F release from an Rb family member by Hsc70 recruited by large T antigen. That is, it now appears that Hsc70 can freely access helix III and the HPD motif of large T antigen bound to an Rb family member.     

A functional role for correlated motion in the N-terminal RNA-binding domain of human U1A protein

The N-terminal RNA-binding domain of the human U1A protein (RBD1) undergoes local conformational changes upon binding to its target RNA. Here, the wild-type RBD1 and two mutants are examined with molecular dynamics simulations that are analyzed using the reorientational eigenmode dynamics (RED) formalism. The results reveal changes in the magnitude and extent of coupled intra-domain motions resulting from single amino acid substitutions. Interpretation of the novel RED results and corresponding NMR relaxation data suggests that the loss of collective motions in the mutants could account for their weak RNA-binding.     

Topology and dynamics of the 10 kDa C-terminal domain of DnaK in solution

Hsp70 molecular chaperones contain three distinct structural domains, a 44 kDa N-terminal ATPase domain, a 17 kDa peptide-binding domain, and a 10 kDa C-terminal domain. The ATPase and peptide binding domains are conserved in sequence and are functionally well characterized. The function of the 10 kDa variable C-terminal domain is less well understood. We have characterized the secondary structure and dynamics of the C-terminal domain from the Escherichia coli Hsp70, DnaK, in solution by high-resolution NMR. The domain was shown to be comprised of a rigid structure consisting of four helices and a flexible C-terminal subdomain of approximately 33 amino acids. The mobility of the flexible region is maintained in the context of the full-length protein and does not appear to be modulated by the nucleotide state. The flexibility of this region appears to be a conserved feature of Hsp70 architecture and may have important functional implications. We also developed a method to analyze 15N nuclear spin relaxation data, which allows us to extract amide bond vector directions relative to a unique diffusion axis. The extracted angles and rotational correlation times indicate that the helices form an elongated, bundle-like structure in solution.     

The response of internal dynamics to hydrophobic core mutations in the SH3 domain from the Fyn tyrosine kinase

We have used (15)N- and (2)H-NMR spin relaxation experiments to study the response of backbone and side-chain dynamics when a leucine or valine is substituted for a completely buried phenylalanine residue in the SH3 domain from the Fyn tyrosine kinase. Several residues show differences in the time scales and temperature dependences of internal motions when data for the three proteins are compared. Changes were also observed in the magnitude of dynamics, with the valine, and to a lesser extent leucine mutant, showing enhanced flexibility compared to the wild-type (WT) protein. The motions of many of the same amide and methyl groups are affected by both mutations, identifying a set of loci where dynamics are sensitive to interactions involving the targeted side chain. These results show that contacts within the hydrophobic core affect many aspects of internal mobility throughout the Fyn SH3 domain.     

Solution structure of GOPC PDZ domain and its interaction with the C-terminal motif of neuroligin

GOPC (Golgi-associated PDZ and coiled-coil motif-containing protein) represents a PDZ domain-containing protein associated with the Golgi apparatus, which plays important roles in vesicular trafficking in secretory and endocytic pathways. GOPC interacts with many other proteins, such as the Wnt receptors Frizzled 8 and neuroligin via its PDZ domain. Neuroligin is a neural cell-adhesion molecule of the post-synapse, which binds to the presynapse molecule neurexin to form a heterotypic intercellular junction. Here we report the solution structure of the GOPC PDZ domain by NMR. Our results show that it is a canonical class I PDZ domain, which contains two alpha-helices and six beta-strands. Using chemical shift perturbation experiments, we further studied the binding properties of the GOPC PDZ domain with the C-terminal motif of neuroligin. The observations showed that the ensemble of the interaction belongs to fast exchange with low affinity. The 3D model of the GOPC PDZ domain/neuroligin C-terminal peptide complex was constructed with the aid of the molecular dynamics simulation method. Our discoveries provide insight into the specific interaction of the GOPC PDZ domain with the C-terminal peptide of Nlg and also provide a general insight about the possible binding mode of the interaction of Nlg with other PDZ domain-containing proteins.     

The role of backbone conformational heat capacity in protein stability: temperature dependent dynamics of the B1 domain of Streptococcal protein G

The contributions of backbone NH group dynamics to the conformational heat capacity of the B1 domain of Streptococcal protein G have been estimated from the temperature dependence of 15N NMR-derived order parameters. Longitudinal (R1) and transverse (R2) relaxation rates, transverse cross-relaxation rates (eta(xy)), and steady state [1H]-15N nuclear Overhauser effects were measured at temperatures of 0, 10, 20, 30, 40, and 50 degrees C for 89-100% of the backbone secondary amide nitrogen nuclei in the B1 domain. The ratio R2/eta(xy) was used to identify nuclei for which conformational exchange makes a significant contribution to R2. Relaxation data were fit to the extended model-free dynamics formalism, incorporating an axially symmetric molecular rotational diffusion tensor. The temperature dependence of the order parameter (S2) was used to calculate the contribution of each NH group to conformational heat capacity (Cp) and a characteristic temperature (T*), representing the density of conformational energy states accessible to each NH group. The heat capacities of the secondary structure regions of the B1 domain are significantly higher than those of comparable regions of other proteins, whereas the heat capacities of less structured regions are similar to those in other proteins. The higher local heat capacities are estimated to contribute up to approximately 0.8 kJ/mol K to the total heat capacity of the B1 domain, without which the denaturation temperature would be approximately 9 degrees C lower (78 degrees C rather than 87 degrees C). Thus, variation of backbone conformational heat capacity of native proteins may be a novel mechanism that contributes to high temperature stabilization of proteins.     

Ligand-induced changes in dynamics in the RT loop of the C-terminal SH3 domain of Sem-5 indicate cooperative conformational coupling

We report the effects of peptide binding on the (15)N relaxation rates and chemical shifts of the C-SH3 of Sem-5. (15)N spin-lattice relaxation time (T(1)), spin-spin relaxation time (T(2)), and ((1)H)-(15)N NOE were obtained from heteronuclear 2D NMR experiments. These parameters were then analyzed using the Lipari-Szabo model free formalism to obtain parameters that describe the internal motions of the protein. High-order parameters (S(2) > 0.8) are found in elements of regular secondary structure, whereas some residues in the loop regions show relatively low-order parameters, notably the RT loop. Peptide binding is characterized by a significant decrease in the (15)N relaxation in the RT loop. Concomitant with the change in dynamics is a cooperative change in chemical shifts. The agreement between the binding constants calculated from chemical shift differences and that obtained from ITC indicates that the binding of Sem-5 C-SH3 to its putative peptide ligand is coupled to a cooperative conformational change in which a portion of the binding site undergoes a significant reduction in conformational heterogeneity.     

The N-terminal domain (IF2N) of bacterial translation initiation factor IF2 is connected to the conserved C-terminal domains by a flexible linker

Bacterial translation initiation factor IF2 is a multidomain protein that is an essential component of a system for ensuring that protein synthesis begins at the correct codon within a messenger RNA. Full-length IF2 from Escherichia coli and seven fragments of the protein were expressed, purified, and characterized using nuclear magnetic resonance (NMR) and circular dichroism (CD) methods. Interestingly, resonances of the 6 kD IF2N domain located at the extreme N terminus of IF2 can be clearly identified within the NMR spectra of the full-length 97-kD protein. (15)N NMR relaxation rate data indicate that (1) the IF2N domain is internally well ordered and tumbles in solution in a manner that is independent of the other domains of the IF2 protein, and (2) the IF2N domain is connected to the C-terminal regions of IF2 by a flexible linker. Chemical shifts of resonances within the isolated IF2N domain do not significantly differ from those of the corresponding residues within the context of the full-length 97-kD protein, indicating that IF2N is a structurally independent unit that does not strongly interact with other regions of IF2. CD and NMR data together provide evidence that Domains I-III of IF2 have unstructured and flexible regions as well as substantial helical content; CD data indicate that the helical content of these regions decreases significantly at temperatures above 35 degrees C. The features of structurally well-ordered N- and C-terminal domains connected by a flexible linker with significant helical content are reminiscent of another translation initiation factor, IF3.     

Regulation of endosomal clathrin and retromer‐mediated endosome to Golgi retrograde transport by the J‐domain protein RME‐8

After endocytosis, most cargo enters the pleiomorphic early endosomes in which sorting occurs. As endosomes mature, transmembrane cargo can be sequestered into inwardly budding vesicles for degradation, or can exit the endosome in membrane tubules for recycling to the plasma membrane, the recycling endosome, or the Golgi apparatus. Endosome to Golgi transport requires the retromer complex. Without retromer, recycling cargo such as the MIG‐14/Wntless protein aberrantly enters the degradative pathway and is depleted from the Golgi. Endosome‐associated clathrin also affects the recycling of retrograde cargo and has been shown to function in the formation of endosomal subdomains. Here, we find that the Caemorhabditis elegans endosomal J‐domain protein RME‐8 associates with the retromer component SNX‐1. Loss of SNX‐1, RME‐8, or the clathrin chaperone Hsc70/HSP‐1 leads to over‐accumulation of endosomal clathrin, reduced clathrin dynamics, and missorting of MIG‐14 to the lysosome. Our results indicate a mechanism, whereby retromer can regulate endosomal clathrin dynamics through RME‐8 and Hsc70, promoting the sorting of recycling cargo into the retrograde pathway.     

Contrasting roles of dynamics in protein allostery: NMR and structural studies of CheY and the third PDZ domain from PSD-95

Allosteric regulation is a ubiquitous phenomenon exploited in biological processes to control cells in a myriad of ways. It is also of emerging interest in the design of functional proteins and therapeutics. Even though allostery was proposed over 50 years ago and has been studied intensively from a structural perspective, many key details of allosteric mechanisms remain mysterious. Over the last decade significant attention has been paid to the “dynamic component” of allostery, as opposed to the analysis of rigid structures. Nuclear magnetic resonance spectroscopy and its ability to detect conformationally dynamic processes at atomic resolution have played an important role in expanding our understanding of allosteric mechanisms and opening up new questions. This article focuses on work that highlights how protein dynamics can factor into allosteric processes in distinct ways. Two cases are contrasted. The first considers the “traditionally allosteric” protein CheY, which undergoes a conformational change as a key element of its allostery. The second considers the more rarely observed “dynamic allostery” in a PDZ domain, in which allosteric behavior arises from changes in internal structural dynamics. Interestingly, the dynamic processes in these two contrasting examples occur on different timescales. In the case of the PDZ domain, subsequent experimental and computational work is reviewed to reveal a more complete picture of this interesting case of allostery.     

NMR solution structure and backbone dynamics of domain III of the E protein of tick-borne Langat flavivirus suggests a potential site for molecular recognition

Flaviviruses cause many human diseases, including dengue fever, yellow fever, West Nile viral encephalitis, and hemorrhagic fevers, and are transmitted to their vertebrate hosts by infected mosquitoes and ticks. Domain III of the envelope protein (E-D3) is considered to be the primary viral determinant involved in the virus-host-cell receptor interaction, and thus represents an excellent target for antiviral drug development. Langat (LGT) virus is a naturally attenuated BSL-2 TBE virus and is a model for the pathogenic BSL-3 and BSL-4 viruses in the serogroup. We have determined the solution structure of LGT-E-D3 using heteronuclear NMR spectroscopy. The backbone dynamics of LGT-E-D3 have been investigated using 15N relaxation measurements. A detailed analysis of the solution structure and dynamics of LGT-E-D3 suggests potential residues that could form a surface for molecular recognition, and thereby represent a target site for antiviral therapeutics design.     

NMR structure of the p63 SAM domain and dynamical properties of G534V and T537P pathological mutants, identified in the AEC syndrome

The p63 protein is crucial for epidermal development, and its mutations cause the extrodactyly ectodermal dysplasia and cleft lip/palate syndrome. The three-dimensional solution structure of the p63 sterile α-motif (SAM) domain (residues 505–579), a region crucial to explaining the human genetic disease ankyloblepharon-ectodermal dysplasia-clefting syndrome (AEC), has been determined by nuclear magnetic resonance spectroscopy. The structure indicates that the domain is a monomer with the characteristic five-helix bundle topology observed in other SAM domains. It includes five tightly packed helices with an extended hydrophobic core to form a globular and compact structure. The dynamics of the backbone and the global correlation time of the molecule have also been investigated and compared with the dynamical properties obtained through molecular dynamics simulation. Attempts to purify the pathological G534V and T537P mutants, originally identified in AEC, were not successful because of the occurrence of unspecific proteolytic degradation of the mutated SAM domains. Analysis of the structural dynamic properties of the G534V and T537P mutants through molecular dynamics simulation and comparison with the wild type permits detection of differences in the degree of free-dom of individual residues and discussion of the possible causes for the pathology.     

Structure and dynamics in solution of the stop codon decoding N-terminal domain of the human polypeptide chain release factor eRF1

The high-resolution NMR structure of the N-domain of human eRF1, responsible for stop codon recognition, has been determined in solution. The overall fold of the protein is the same as that found in the crystal structure. However, the structures of several loops, including those participating in stop codon decoding, are different. Analysis of the NMR relaxation data reveals that most of the regions with the highest structural discrepancy between the solution and solid states undergo internal motions on the ps-ns and ms time scales. The NMR data show that the N-domain of human eRF1 exists in two conformational states. The distribution of the residues having the largest chemical shift differences between the two forms indicates that helices α2 and α3, with the NIKS loop between them, can switch their orientation relative to the β-core of the protein. Such structural plasticity may be essential for stop codon recognition by human eRF1.     

Characterization of the hinges of the effector loop in the reaction pathway of the activation of ras-proteins. Kinetics of binding of beryllium trifluoride to V29G and I36G mutants of Ha-ras-p21.

This work experimentally confirms the pathway of activation of Ha-ras-p21, which was calculated by the method of Targeted Molecular Dynamics (TMD) (Díaz JF, Wroblowski B, Schlitter J, Engelborghs Y, 1997a, Proteins Struct Funct Genet 28:434-451). The process can be studied experimentally by analyzing the binding of BeF3- to the GDP complex of the active fluorescent mutant Y32W (Díaz JF, Sillen A, Engelborghs Y, 1997b, J Biol Chem 227:23138-23143). Two mutants, V29G and 136G, have been constructed at both sides of the effector loop of the active fluorescent mutant. This was done to check the proposed reaction pathway and to provide further insight into the mechanism of the activation of ras proteins. Both mutations accelerate the conformational isomerization with two orders of magnitude, demonstrating convincingly the role of these residues as hinges of the effector loop in one or more of the transitions of the conformational change. These results provide experimental support to the pathway calculated by TMD analysis.     

The conformation of the extracellular binding domain of Death Receptor 5 in the presence and absence of the activating ligand TRAIL: a molecular dynamics study

The Death Receptor 5 (DR5), a member of tumor necrosis factor receptor (TNFR) superfamily of receptors, triggers apoptosis (programmed cell death) when stimulated by its tridentate ligand TRAIL. Until recently it was generally assumed that the activation of DR5 resulted from the recruitment of three independent receptor units, leading to the trimerization of intracellular domains. However, there is mounting evidence to suggest that, in the absence of ligand, such cytokine receptors primarily reside as preformed complexes. In this work, molecular dynamics simulations of the TRAIL-DR5 complex, the unbound receptor trimer and individual receptor monomers are compared to gain insight in the mechanism of activation. The results suggest that, in the absence of TRAIL, DR5 has a strong propensity to self-associate and that this is primarily mediated through interactions of the membrane proximal domains. The association of the free receptors leads to a loss of the threefold symmetry found within the receptor-ligand complex. The simulations suggest that the primary role of TRAIL is to induce threefold-symmetry within the DR5 complex and to constrain the receptor to a specific conformation. The implications of this in terms of the mechanism by which the receptor switches from an inactive to an active state are discussed.     

Analysis of the binding loops configuration and surface adaptation of different crystallized single‐domain antibodies in response to various antigens

Monoclonal antibodies have revolutionized the biomedical field through their ubiquitous utilization in different diagnostics and therapeutic applications. Despite this widespread use, their large size and structural complexity have limited their versatility in specific applications. The antibody variable region that is responsible for binding antigen is embodied within domains that can be rescued individually as single‐domain antibody (sdAb) fragments. Because of the unique characteristics of sdAbs, such as low molecular weight, high physicochemical stability, and the ability to bind antigens inaccessible to conventional antibodies, they represent a viable alternative to full‐length antibodies. Consequently, 149 crystal structures of sdAbs, originating from human (VH), camelids (VHH), or sharks (VNAR), were retrieved from the Protein Data Bank, and their structures were compared. The 3 types of sdAbs displayed complementarity determining regions (CDRs) with different lengths and configurations. CDR3 of the VHH and VNAR domains were dominated by pleated and extended orientations, respectively. Although VNAR showed the smallest average molecular weight and molecular surface area compared with VHH and VH antibodies. However, the solvent accessible surface area measurements of the 3 tested sdAbs types were very similar. All the antihapten VHH antibodies showed pleated CDR3, which were sufficient to create a binding pocket to accommodate haptens (methotrexate and azo dyes) in terms of shape and electrostatic potential. The sdAbs that recognized lysozyme showed more diversity in their CDR3 orientation to enable them to recognize various topographies of lysozyme. Subsequently, the three sdAb classes were different in size and surface area and have shown distinguishable ability to optimize their CDR length and orientation to recognize different antigen classes.     

Molecular dynamics study of segment peptides of Bax,Bim, and Mcl-1 BH3 domain of the apoptosis-regulating proteins bound to the anti-apoptotic Mcl-1 protein

Mcl-1 has emerged as a potential therapeutic target in the treatment of several malignancies. Peptides representing BH3 region of pro-apoptotic proteins have been shown to bind the hydrophobic cleft of anti-apoptotic Mcl-1 and this segment is responsible for modulating the apoptotic pathways in living cells. Understanding the molecular basis of protein–peptide interaction is required to develop potent inhibitors specific for Mcl-1. Molecular dynamics simulations were performed for Mcl-1 in complex with three different BH3 peptides derived from Mcl-1, Bax, and Bim. Accordingly, the calculated binding free energies using MM-PBSA method are obtained and comparison with the experimentally determined binding free energies is made. The interactions involving two conserved charged residues (Aspi, and Arg/Lysi-4) and three upstream conserved hydrophobic residues (Leui-5, Ile/Vali-2, and Glyi-1, respectively) of BH3 peptides play the important roles in the structural stability of the complexes. The calculated results exhibit that the interactions of Bim BH3 peptides to Mcl-1 is stronger than the complex with Bax 19BH3 peptides. The hydrophobic residues (position i???9, i???8 and i?+?2) of BH3 peptides can be involved in their inhibitory specificity. The calculated results can be used for designing more effective MCL-1 inhibitors.     

Crystal structure of the streptococcal superantigen SpeI and functional role of a novel loop domain in T cell activation by group V superantigens

Superantigens (SAgs) are potent microbial toxins that bind simultaneously to T cell receptors (TCRs) and class II major histocompatibility complex molecules, resulting in the activation and expansion of large T cell subsets and the onset of numerous human diseases. Within the bacterial SAg family, streptococcal pyrogenic exotoxin I (SpeI) has been classified as belonging to the group V SAg subclass, which are characterized by a unique, relatively conserved ∼15 amino acid extension (amino acid residues 154 to 170 in SpeI; herein referred to as the α3-β8 loop), absent in SAg groups I through IV. Here, we report the crystal structure of SpeI at 1.56 Å resolution. Although the α3-β8 loop in SpeI is several residues shorter than that of another group V SAg, staphylococcal enterotoxin serotype I, the C-terminal portions of these loops, which are located adjacent to the putative TCR binding site, are structurally similar. Mutagenesis and subsequent functional analysis of SpeI indicates that TCR β-chains are likely engaged in a similar general orientation as other characterized SAgs. We show, however, that the α3-β8 loop length, and the presence of key glycine residues, are necessary for optimal activation of T cells. Based on Vβ-skewing analysis of human T cells activated with SpeI and structural models, we propose that the α3-β8 loop is positioned to form productive intermolecular contacts with the TCR β-chain, likely in framework region 3, and that these contacts are required for optimal TCR recognition by SpeI, and likely all other group V SAgs.