Solution structure and backbone dynamics of the AF-6 PDZ domain/Bcr peptide complex

The human AF-6, a scaffold protein between cell membrane-associated proteins and the actin cytoskeleton, plays an important role in special cell-cell junctions and signal transduction. It can be phosphorylated by the protein kinase Bcr, which allows efficient binding of the C terminus of Bcr to the PDZ domain of AF-6 and consequently enhances the binding affinity of AF-6 to Ras. Formation of the AF-6, Bcr, and Ras ternary complex results in down-regulation of the Ras-mediated signal transduction pathway. To better understand the molecular basis for the recognition of the AF-6 PDZ domain and Bcr, we solve the solution structure of the AF-6 PDZ domain complexed with the C-terminal peptide of Bcr and explore the interactions between them in detail. Compared with previously reported structures, the complex exhibits a noncanonical binding mode of PDZ/peptide. Owing to the distinct residues involved in the AF-6 PDZ domain and Bcr peptide interaction, the interaction mode does not adapt to the existing classification rules that have been put forward, based on the ligand or the PDZ domain specificity. Furthermore, the PDZ domain of AF-6 can bind to the C terminus of Bcr efficiently after phosphorylation of AF-6 by the Bcr kinase. The phosphorylation may induce a conformational change of AF-6, which makes the binding surface on the PDZ domain accessible to Bcr for efficient binding. This study not only characterizes the structural details of the AF-6 PDZ/Bcr peptide complex, but also provides a potential target for future drug design and disease therapy.     

Solution structure and backbone dynamics of the DNA-binding domain of mouse Sox-5

The fold of the murine Sox-5 (mSox-5) HMG box in free solution has been determined by multidimensional NMR using (15)N-labeled protein and has been found to adopt the characteristic twisted L-shape made up of two wings: the major wing comprising helix 1 (F10--F25) and helix 2 (N32--A43), the minor wing comprising helix 3 (P51--Y67) in weak antiparallel association with the N-terminal extended segment. (15)N relaxation measurements show considerable mobility (reduced order parameter, S(2)) in the minor wing that increases toward the amino and carboxy termini of the chain. The mobility of residues C-terminal to Q62 is significantly greater than the equivalent residues of non-sequence-specific boxes, and these residues show a weaker association with the extended N-terminal segment than in non-sequence boxes. Comparison with previously determined structures of HMG boxes both in free solution and complexed with DNA shows close similarity in the packing of the hydrophobic cores and the relative disposition of the three helices. Only in hSRY/DNA does the arrangement of aromatic sidechains differ significantly from that of mSox-5, and only in rHMG1 box 1 bound to cisplatinated DNA does helix 1 have no kink. Helix 3 in mSox-5 is terminated by P68, a conserved residue in DNA sequence-specific HMG boxes, which results in the chain turning through approximately 90 degrees.     

Solution structure and backbone dynamics of the TGFbeta type II receptor extracellular domain

Isoforms of transforming growth factor beta (TGFbeta) are 25 kDa homodimeric polypeptides that signal by binding and bringing together two related, functionally distinct cell surface receptors designated as TbetaR1 and TbetaR2. Here, we report the solution structure of the 13.8 kDa extracellular domain of human TbetaR2 (ecTbetaR2) as calculated from N(N)-H(N), C(alpha)-H(alpha), and C(alpha)-C(O) residual dipolar coupling restraints in conjunction with NOE distance, dihedral angle, and scalar coupling restraints. Comparison of the free ecTbetaR2 solution structure with the TGFbeta3-bound ecTbetaR2 crystal structure reveals backbone conformations that superimpose with RMSDs of 1.0 A over the regions of regular secondary structure and 1.4 A overall. The differences in structure fall mainly in loop regions that are either poorly defined by the available NMR data or are involved in crystal contacts. The noted similarities between the NMR structure of the free form and the crystal structure of the TGFbeta-bound form are also consistent with the close correspondence, 0.16 A RMSD for regions of secondary structure and 0.51 A RMSD overall, for the crystal structure of free ecTbetaR2 as compared to the crystal structure of TGFbeta3-bound ecTbetaR2. Despite the apparent similarities between the free and the bound forms, there appears to be small but significant differences in structure involving the interfacial contact region of the receptor. Measurements of backbone (15)N relaxation times and interpretation of these by the model-free formalism with axial diffusional anisotropy further reveal significant ms to micros time scale motions centered about two of the conserved disulfide bonds and in several residues that comprise the TGFbeta binding surface. Together, these observations indicate that binding likely occurs through a mechanism with a small component of induced fit character, whereby flexibility within the receptor facilitates the transition to the TGFbeta-bound state.     

Phosphorylation of a PDZ domain extension modulates binding affinity and interdomain interactions in postsynaptic density-95 (PSD-95) protein, a membrane-associated guanylate kinase (MAGUK)

Postsynaptic density-95 is a multidomain scaffolding protein that recruits glutamate receptors to postsynaptic sites and facilitates signal processing and connection to the cytoskeleton. It is the leading member of the membrane-associated guanylate kinase family of proteins, which are defined by the PSD-95/Discs large/ZO-1 (PDZ)-Src homology 3 (SH3)-guanylate kinase domain sequence. We used NMR to show that phosphorylation of conserved tyrosine 397, which occurs in vivo and is located in an atypical helical extension (α3), initiates a rapid equilibrium of docked and undocked conformations. Undocking reduced ligand binding affinity allosterically and weakened the interaction of PDZ3 with SH3 even though these domains are separated by a ~25-residue linker. Additional phosphorylation at two linker sites further disrupted the interaction, implicating α3 and the linker in tuning interdomain communication. These experiments revealed a novel mode of regulation by a detachable PDZ element and offer a first glimpse at the dynamic interaction of PDZ and SH3-guanylate kinase domains in membrane-associated guanylate kinases.     

Solution structure and backbone dynamics of the defunct domain of calcium vector protein.

CaVP (calcium vector protein) is a Ca(2+) sensor of the EF-hand protein family which is highly abundant in the muscle of Amphioxus. Its three-dimensional structure is not known, but according to the sequence analysis, the protein is composed of two domains, each containing a pair of EF-hand motifs. We determined recently the solution structure of the C-terminal domain (Trp81-Ser161) and characterized the large conformational and dynamic changes induced by Ca(2+) binding. In contrast, the N-terminal domain (Ala1-Asp86) has lost the capacity to bind the metal ion due to critical mutations and insertions in the two calcium loops. In this paper, we report the solution structure of the N-terminal domain and its backbone dynamics based on NMR spectroscopy, nuclear relaxation, and molecular modeling. The well-resolved three-dimensional structure is typical of a pair of EF-hand motifs, joined together by a short antiparallel beta-sheet. The tertiary arrangement of the two EF-hands results in a closed-type conformation, with near-antiparallel alpha-helices, similar to other EF-hand pairs in the absence of calcium ions. To characterize the internal dynamics of the protein, we measured the (15)N nuclear relaxation rates and the heteronuclear NOE effect in (15)N-labeled N-CaVP at a magnetic field of 11.74 T and 298 K. The domain is mainly monomeric in solution and undergoes an isotropic Brownian rotational diffusion with a correlation time of 7.1 ns, in good agreement with the fluorescence anisotropy decay measurements. Data analysis using a model-free procedure showed that the amide backbone groups in the alpha-helices and beta-strands undergo highly restricted movements on a picosecond to nanosecond time scale. The amide groups in Ca(2+) binding loops and in the linker fragment also display rapid fluctuations with slightly increased amplitudes.     

Solution structure and backbone dynamics of the human DNA ligase IIIalpha BRCT domain.

BRCT (BRCA1 carboxyl terminus) domains are found in a number of DNA repair enzymes and cell cycle regulators and are believed to mediate important protein-protein interactions. The DNA ligase IIIalpha BRCT domain partners with the distal BRCT domain of the DNA repair protein XRCC1 (X1BRCTb) in the DNA base excision repair (BER) pathway. To elucidate the mechanisms by which these two domains can interact, we have determined the solution structure of human ligase IIIalpha BRCT (L3[86], residues 837-922). The structure of L3[86] consists of a beta2beta1beta3beta4 parallel sheet with a two-alpha-helix bundle packed against one face of the sheet. This fold is conserved in several proteins having a wide range of activities, including X1BRCTb [Zhang, X. D., et al. (1998) EMBO J. 17, 6404-6411]. L3[86] exists as a dimer in solution, but an insufficient number of NOE restraints precluded the determination of the homodimer structure. However, 13C isotope-filtered and hydrogen-deuterium exchange experiments indicate that the N-terminus, alpha1, the alpha1-beta2 loop, and the three residues following alpha2 are involved in forming the dimer interface, as similarly observed in the structure of X1BRCTb. NOE and dynamic data indicate that several residues (837-844) in the N-terminal region appear to interconvert between helix and random coil conformations. Further studies of other BRCT domains and of their complexes are needed to address how these proteins interact with one another, and to shed light on how mutations can lead to disruption of function and ultimately disease.     

Solution structure and backbone dynamics of the non-receptor protein-tyrosine kinase-6 Src homology 2 domain

Human protein-tyrosine kinase-6 (PTK6, also known as breast tumor kinase (Brk)) is a member of the non-receptor protein-tyrosine kinase family and is expressed in two-thirds of all breast tumors. To understand the structural basis of PTK6 function, we have determined the solution structure and backbone dynamics of the PTK6-Src homology 2 (SH2) domain using multidimensional NMR spectroscopy. The solution structure clearly indicates that the SH2 domain of human PTK6 contains a consensus alpha/beta-fold and a Tyr(P) peptide binding surface, which are common to other SH2 domains. However, two of the alpha-helices (alphaA and alphaB) are located on opposite faces of the central beta-sheet. In addition, the topological arrangement of a central four-stranded antiparallel beta-sheet (strands betaA, betaB, betaC, and betaD) differs from that of other Src family members. Backbone dynamics and Tyr(P) peptide titration experiments revealed that the putative ligand binding sites of the PTK6-SH2 domain undergo distinctive internal motions when compared with other regions of the protein. Surface plasmon resonance analysis showed that the Tyr(P) peptide had a dissociation constant of about 60 microm, which is substantially weaker binding than previously reported for Src family members. The solution structure together with data from the ligand binding mode of PTK6-SH2 provides insight into the molecular basis of the autoinhibitory role of PTK6.     

Structure, dynamics and binding characteristics of the second PDZ domain of PTP-BL

The PDZ domains of the protein tyrosine phosphatase PTP-BL mediate interactions by binding to specific amino acid sequences in target proteins. The solution structure of the second PDZ domain of PTP-BL, PDZ2, displays a compact fold with six β strands and two α-helices. A unique feature of this domain compared to the canonical PDZ fold is an extended flexible loop at the base of the binding pocket, termed L1, that folds back onto the protein backbone, a feature that is shared by both the murine and human orthologues. The structure of PDZ2 differs significantly from the orthologous human structure. A comparison of structural quality indicators clearly demonstrates that the PDZ2 ensemble is statistically more reasonable than that of the human orthologue. The analysis of ¹⁵N relaxation data for PDZ2 shows a normal pattern, with more rigid secondary structures and more flexible loop structures. Close to the binding pocket, Leu85 and Thr88 display greater mobility when compared to surrounding residues. Peptide binding studies demonstrated a lack of interaction between murine PDZ2 and the C terminus of the murine Fas/CD95 receptor, suggesting that the Fas/CD95 receptor is not an in vivo target for PDZ2. In addition, PDZ2 specifically binds the C termini of both human Fas/CD95 receptor and the RIL protein, despite RIL containing a non-canonical PDZ-interacting sequence of E-x-V. A model of PDZ2 with the RIL peptide reveals that the PDZ2 binding pocket is able to accommodate the bulkier side-chain of glutamic acid while maintaining crucial protein to peptide hydrogen bond interactions.     

Solution structure and backbone dynamics of an omega-conotoxin precursor

Nuclear magnetic resonance spectroscopy was used to characterize the solution structure and backbone dynamics of a putative precursor form of omega-conotoxin MVIIA, a 25-amino-acid residue peptide antagonist of voltage-gated Ca(2+) channels. The mature peptide is found in the venom of a fish-hunting marine snail Conus magus and contains an amidated carboxyl terminus that is generated by oxidative cleavage of a Gly residue. The form examined in this study is identical to the mature peptide except for the presence of the unmodified carboxy-terminal Gly. This form, referred to as omega-MVIIA-Gly, has previously been shown to refold and form its disulfides more efficiently than the mature form, suggesting that the presence of the terminal Gly may favor folding in vivo. The nuclear magnetic resonance (NMR) structure determination indicated that the fold of omega-MVIIA-Gly is very similar to that previously determined for the mature form, but revealed that the terminal Gly residue participates in a network of hydrogen bonds involving both backbone and side chain atoms, very likely accounting for the enhanced stability and folding efficiency. (15)N relaxation experiments indicated that the backbone is well ordered on the nanosecond time scale but that residues 9-15 undergo a conformational exchange processes with a time constant of approximately 35 microseconds. Other studies have implicated this segment in the binding of the peptide to its physiological target, and the observed motions may play a role in allowing the peptide to enter the binding site     

Formation of nNOS/PSD-95 PDZ dimer requires a preformed beta-finger structure from the nNOS PDZ domain

PDZ domains are modular protein units that play important roles in organizing signal transduction complexes. PDZ domains mediate interactions with both C-terminal peptide ligands and other PDZ domains. Here, we used PDZ domains from neuronal nitric oxide synthase (nNOS) and postsynaptic density protein-95 (PSD-95) to explore the mechanism for PDZ-dimer formation. The nNOS PDZ domain terminates with a approximately 30 residue amino acid beta-finger peptide that is shown to be required for nNOS/PSD-95 PDZ dimer formation. In addition, formation of the PDZ dimer requires this beta-finger peptide to be physically anchored to the main body of the canonical nNOS PDZ domain. A buried salt bridge between the beta-finger and the PDZ domain induces and stabilizes the beta-hairpin structure of the nNOS PDZ domain. In apo-nNOS, the beta-finger peptide is partially flexible and adopts a transient beta-strand like structure that is stabilized upon PDZ dimer formation. The flexibility of the NOS PDZ beta-finger is likely to play a critical role in supporting the formation of nNOS/PSD-95 complex. The experimental data also suggest that nNOS PDZ and the second PDZ domain of PSD-95 form a "head-to-tail" dimer similar to the nNOS/syntrophin complex characterized by X-ray crystallography.     

Solution structure and backbone dynamics of the XPC-binding domain of the human DNA repair protein hHR23B

Human cells contain two homologs of the yeast RAD23 protein, hHR23A and hHR23B, which participate in the DNA repair process. hHR23B houses a domain (residues 277-332, called XPCB) that binds specifically and directly to the xeroderma pigmentosum group C protein (XPC) to initiate nucleotide excision repair (NER). This domain shares sequence homology with a heat shock chaperonin-binding motif that is also found in the stress-inducible yeast phosphoprotein STI1. We determined the solution structure of a protein fragment containing amino acids 275-342 of hHR23B (termed XPCB-hHR23B) and compared it with the previously reported solution structures of the corresponding domain of hHR23A. The periodic positioning of proline residues in XPCB-hHR23B produced kinked alpha helices and assisted in the formation of a compact domain. Although the overall structure of the XPCB domain was similar in both XPCB-hHR23B and XPCB-hHR23A, the N-terminal part (residues 275-283) of XPCB-hHR23B was more flexible than the corresponding part of hHR23A. We tried to infer the characteristics of this flexibility through (15)N-relaxation studies. The hydrophobic surface of XPCB-hHR23B, which results from the diverse distribution of N-terminal region, might give rise to the functional pleiotropy observed in vivo for hHR23B, but not for hHR23A.     

Solution structure and backbone dynamics of the Cu(I) and apo forms of the second metal-binding domain of the Menkes protein ATP7A

The second domain of the human Menkes protein (MNK2), formed by 72 residues, has been expressed in Escherichia coli, and its structure has been determined by NMR in both the apo and copper-loaded forms. The structures, obtained with (13)C- and (15)N-labeled samples, are of high quality with backbone rmsd values of 0.51 and 0.41 A and CYANA target functions of 0.39 and 0.38 A(2), respectively. The loop involved in copper binding is part of a hydrophobic patch, which is maintained in both forms. Conformational mobility is observed in the apo form in the same loop. A comparison with metallochaperones and soluble domains of P-type ATPases allows us to relate the primary structure to the occurrence of structural rearrangements upon copper binding.     

Flavonoids from Radix Scutellariae as potential stroke therapeutic agents by targeting the second postsynaptic density 95 (PSD-95)/disc large/zonula occludens-1 (PDZ) domain of PSD-95.

Excessive activation of N-methyl-D-aspartate receptors (NMDARs) and subsequent production of nitric oxide by neuronal nitric oxide synthase (nNOS) contribute to neuronal damage resulting from hypoxic and ischemic insults. NMDARs and nNOS are coupled together at the postsynaptic membrane through their interaction with postsynaptic density protein (PSD) 95 via PSD-95/disc large/zonula occludens-1 (PDZ) domains. We used NMR (nuclear magnetic resonance) spectroscopy to screen medicinal herbs used in traditional Chinese medicine (TCM) stroke therapy for compounds binding to the second PDZ domain (PDZ2) of PSD-95, the domain linking nNOS and PSD-95. Aqueous extract of Huangqin, the root of Scutellaria baicalensis Georgi (Labiatae), showed significant binding to PDZ2 of PSD-95. The binding site of the active components in the extract overlapped with the nNOS/NR2B-binding pocket of PDZ2 of PSD-95. Four flavones, baicalin, norwogonoside, oroxylin A-glucuronide (oroxyloside), and wogonoside were isolated and found to account for the PDZ-binding activity of the extract. NMR titration experiments showed that baicalin and norwogonoside displayed the highest PDZ2 binding affinity, while oroxylin A-glucuronide and wogonoside showed 4-5 fold less potency in binding to the PDZ domain. Identification of the PDZ binding activity of these compounds will allow investigating whether or not it contributes to the observed clinical effects of Radix Scutellariae. Furthermore, these molecules might provide leads for the development of drugs targeting the signaling pathways mediated by PDZ domains.     

Solution structure and backbone dynamics of an N-terminal ubiquitin-like domain in the GLUT4-regulating protein, TUG

The GLUT4-regulating protein, TUG, functions to retain GLUT4-containing membrane vesicles intracellularly and, in response to insulin stimulation, releases these vesicles to the cellular exocytic machinery for translocation to the plasma membrane. As part of our on going effort to describe the molecular basis for TUG function, we have determined the tertiary structure and characterized the backbone dynamics for an N-terminal ubiquitin-like domain (TUG-UBL1) using NMR spectroscopy. A well-ordered conformation is observed for residues 10-83 of full-length TUG, and confirms a beta-grasp or ubiquitin-like topology. Although not required for in vitro association with GLUT4, the functional role of the TUG-UBL1 domain has not yet been described. We undertook a limited literature review of similar N-terminal UBL domains and noted that a majority participate in protein-protein interactions, generally functioning as adaptor modules to physically associate the over all activity of the protein with a specific cellular process, such as the ubiquitin-proteasome pathway. In consistent fashion, TUG-UBL1 is not expected to participate in a covalent protein modification reaction as it lacks the characteristic C-terminal diglycine ("GG") motif required for conjugation to an acceptor lysine, and also lacks the three most common acceptor lysine residues involved in polyubiquitination. Additionally, analysis of the TUG-UBL1 molecular surface reveals a lack of conservation of the "Ile-44 hydrophobic face" typically involved in ubiquitin recognition. Instead, we speculate on the possible significance of a concentrated area of negative electrostatic potential with increased backbone mobility, both of which are features suggestive of a potential protein-protein interaction site.     

Solution structure and backbone dynamics of an antigen-free heavy chain variable domain (VHH) from Llama

Camelids, (dromedaries, camels, and llamas) produce heavy-chains antibodies, with their antigen recognition sites composed of a single VH-like domain, referred to as VHH. The solution structure of one of these VHHs domains (VHH-H14), raised against the alpha subunit of the human chorionic gonadotropin hormone (hCG), has been determined by (15)N heteronuclear three-dimensional NMR spectroscopy. The framework is well resolved within the set of 20 best-calculated NMR structures and is close to that of classical VH domains from vertebrate antibodies, consisting of two antiparallel beta-sheets organized in a beta-barrel. Loops display a lower precision, especially the Complementarity Determining Regions (CDRs), involved in antigen recognition. Comparison of the three-dimensional VHH-H14 solution structure with its previously solved crystal structure (Spinelli et al., Nature Struct. Biol. 1996;3:752-757) reveals a high similarity to the framework, whereas significant conformational differences occur on CDRs, leading to the assumption that the antigen recognition site is a more mobile part. In order to deepen our insights into the dynamics of VHH-H14 in solution, (15)N relaxation was measured with longitudinal R1 and transverse R2 self-relaxation rates, and (15)N steady-state heteronuclear nuclear Overhauser enhancements (NOE), making it possible to probe picosecond-to-millisecond internal motions. Determination of dynamic parameters (S(2), tau(e), and Rex) through the Lipari-Szabo Model-free approach enables the identification of several regions with enhanced dynamics. Especially, the mobility measurements from NMR confirm that the antigen recognition site is the most mobile part of the VHH-H14 domain on picosecond-to-nanosecond fast time scales. Several residues belonging to the three CDRs are submitted to chemical exchange processes occurring on slow microsecond-to-millisecond time scales, suggesting that the formation of the VHH/antigen complex should be accompanied by structural changes.     

NMR structure and backbone dynamics of the extended second transmembrane domain of the human neuronal glycine receptor alpha1 subunit

The structure and backbone dynamics of an extended second transmembrane segment (TM2e) of the human neuronal glycine receptor alpha(1) subunit in sodium dodecyl sulfate micelles were studied by (1)H and (15)N solution-state NMR. The 28-amino acid segment contained the consensus TM2 domain plus part of the linker between the second and third transmembrane domains. The presence of a well-structured helical region of at least 13 amino acids long and an unstructured region near the linker was evident from the proton chemical shifts and the pattern of midrange nuclear Overhauser effects (NOE). (15)N relaxation rate constants, R(1) and R(2), and (15)N-[(1)H] NOE indicated restricted internal motions in the helical region with NOE values between 0.6 and 0.8. The squared order parameter (S(2)), the effective correlation time for fast internal motions (tau(e)), and the global rotational correlation time (tau(m)) were calculated for all TM2e backbone N-H bonds using the model-free approach. The S(2) values ranged about 0.75-0.86, and the tau(e) values were below 100 ps for most of the residues in the helical region. The tau(m) value, calculated from the dynamics of the helical region, was 5.1 ns. The S(2) values decreased to 0.1, and the tau(e) values sharply increased up to 1.2 ns at the linker near the C-terminus, indicating that the motion of this region is unrestricted. The results suggest a relatively high degree of motional freedom of TM2e in micelles and different propensities of the N- and C-terminal moieties of the transmembrane domain to assume stable helical structures.     

Solution structure and backbone dynamics of the catalytic domain of matrix metalloproteinase-2 complexed with a hydroxamic acid inhibitor

MMP-2 is a member of the matrix metalloproteinase family that has been implicated in tumor cell metastasis and angiogenesis. Here, we describe the solution structure of a catalytic domain of MMP-2 complexed with a hydroxamic acid inhibitor (SC-74020), determined by three-dimensional heteronuclear NMR spectroscopy. The catalytic domain, designated MMP-2C, has a short peptide linker replacing the internal fibronectin-domain insertion and is enzymatically active. Distance geometry-simulated annealing calculations yielded 14 converged structures with atomic root-mean-square deviations (r.m.s.d.) of 1.02 and 1.62 A from the mean coordinate positions for the backbone and for all heavy atoms, respectively, when 11 residues at the N-terminus are excluded. The structure has the same global fold as observed for other MMP catalytic domains and is similar to previously solved crystal structures of MMP-2. Differences observed between the solution and the crystal structures, near the bottom of the S1' specificity loop, appear to be induced by the large inhibitor present in the solution structure. The MMP-2C solution structure is compared with MMP-8 crystal structure bound to the same inhibitor to highlight the differences especially in the S1' specificity loop. The finding provides a structural explanation for the selectivity between MMP-2 and MMP-8 that is achieved by large inhibitors.