O-fucosylation stabilizes the TSR3 motif in thrombospondin-1 by interacting with nearby amino acids and protecting a disulfide bond.

Thrombospondin type-1 repeats (TSRs) are small protein motifs containing six conserved cysteines forming three disulfide bonds that can be modified with an O-linked fucose. Protein O-fucosyltransferase 2 (POFUT2) catalyzes the addition of O-fucose to TSRs containing the appropriate consensus sequence, and the O-fucose modification can be elongated to a Glucose-Fucose disaccharide with the addition of glucose by ß3-glucosyltransferase (B3GLCT). Elimination of Pofut2 in mice results in embryonic lethality in mice, highlighting the biological significance of O-fucose modification on TSRs. Knockout of POFUT2 in HEK293T cells has been shown to cause complete or partial loss of secretion of many proteins containing O-fucosylated TSRs. In addition, POFUT2 is localized to the endoplasmic reticulum (ER) and only modifies folded TSRs, stabilizing their structures. These observations suggest that POFUT2 is involved in an ER quality control mechanism for TSR folding and that B3GLCT also participates in quality control by providing additional stabilization to TSRs. However, the mechanisms by which addition of these sugars result in stabilization are poorly understood. Here, we conducted molecular dynamics (MD) simulations and provide crystallographic and NMR evidence that the Glucose-Fucose disaccharide interacts with specific amino acids in the TSR3 domain in thrombospondin-1 that are within proximity to the O-fucosylation modification site resulting in protection of a nearby disulfide bond. We also show that mutation of these amino acids reduces the stabilizing effect of the sugars in vitro. These data provide mechanistic details regarding the importance of O-fucosylation and how it participates in quality control mechanisms inside the ER.

Accurate de novo design of hyperstable constrained peptides.

Naturally occurring, pharmacologically active peptides constrained with covalent crosslinks generally have shapes that have evolved to fit precisely into binding pockets on their targets. Such peptides can have excellent pharmaceutical properties, combining the stability and tissue penetration of small-molecule drugs with the specificity of much larger protein therapeutics. The ability to design constrained peptides with precisely specified tertiary structures would enable the design of shape-complementary inhibitors of arbitrary targets. Here we describe the development of computational methods for accurate de novo design of conformationally restricted peptides, and the use of these methods to design 18-47 residue, disulfide-crosslinked peptides, a subset of which are heterochiral and/or N-C backbone-cyclized. Both genetically encodable and non-canonical peptides are exceptionally stable to thermal and chemical denaturation, and 12 experimentally determined X-ray and NMR structures are nearly identical to the computational design models. The computational design methods and stable scaffolds presented here provide the basis for development of a new generation of peptide-based drugs.

Skp1 Dimerization Conceals Its F-Box Protein Binding Site.

Skp1 is an adapter that links F-box proteins to cullin-1 in the Skp1/cullin-1/F-box (SCF) protein family of E3 ubiquitin ligases that targets specific proteins for polyubiquitination and subsequent protein degradation. Skp1 from the amoebozoan Dictyostelium forms a stable homodimer in vitro with a Kd of 2.5 µM as determined by sedimentation velocity studies yet is monomeric in crystal complexes with F-box proteins. To investigate the molecular basis for the difference, we determined the solution NMR structure of a doubly truncated Skp1 homodimer (Skp1ΔΔ). The solution structure of the Skp1ΔΔ dimer reveals a 2-fold symmetry with an interface that buries ∼750 Å2 of predominantly hydrophobic surface. The dimer interface overlaps with subsite 1 of the F-box interaction area, explaining why only the Skp1 monomer binds F-box proteins (FBPs). To confirm the model, Rosetta was used to predict amino acid substitutions that might disrupt the dimer interface, and the F97E substitution was chosen to potentially minimize interference with F-box interactions. A nearly full-length version of Skp1 with this substitution (Skp1ΔF97E) behaved as a stable monomer at concentrations of ≤500 µM and actively bound a model FBP, mammalian Fbs1, which suggests that the dimeric state is not required for Skp1 to carry out a basic biochemical function. Finally, Skp1ΔF97E is expected to serve as a monomer model for high-resolution NMR studies previously hindered by dimerization.

The copBL operon protects Staphylococcus aureus from copper toxicity: CopL is an extracellular membrane-associated copper-binding protein.

As complications associated with antibiotic resistance have intensified, copper (Cu) is attracting attention as an antimicrobial agent. Recent studies have shown that copper surfaces decrease microbial burden, and host macrophages use Cu to increase bacterial killing. Not surprisingly, microbes have evolved mechanisms to tightly control intracellular Cu pools and protect against Cu toxicity. Here, we identified two genes (copB and copL) encoded within the Staphylococcus aureus arginine-catabolic mobile element (ACME) that we hypothesized function in Cu homeostasis. Supporting this hypothesis, mutational inactivation of copB or copL increased copper sensitivity. We found that copBL are co-transcribed and that their transcription is increased during copper stress and in a strain in which csoR, encoding a Cu-responsive transcriptional repressor, was mutated. Moreover, copB displayed genetic synergy with copA, suggesting that CopB functions in Cu export. We further observed that CopL functions independently of CopB or CopA in Cu toxicity protection and that CopL from the S. aureus clone USA300 is a membrane-bound and surface-exposed lipoprotein that binds up to four Cu+ ions. Solution NMR structures of the homologous Bacillus subtilis CopL, together with phylogenetic analysis and chemical-shift perturbation experiments, identified conserved residues potentially involved in Cu+ coordination. The solution NMR structure also revealed a novel Cu-binding architecture. Of note, a CopL variant with defective Cu+ binding did not protect against Cu toxicity in vivo Taken together, these findings indicate that the ACME-encoded CopB and CopL proteins are additional factors utilized by the highly successful S. aureus USA300 clone to suppress copper toxicity.

Glycosylation Promotes the Random Coil to Helix Transition in a Region of a Protist Skp1 Associated with F-Box Binding.

Cullin-ring-ligases mediate protein polyubiquitination, a signal for degradation in the 26S proteasome. The CRL1 class consists of Skp1/cullin-1/F-box protein/Rbx1 (SCF) complexes that cyclically associate with ubiquitin-E2 to build the polyubiquitin chain. Within the SCF complex, the 162-amino acid DdSkp1 from Dictyostelium bridges cullin-1 with an F-box protein (FBP), the specificity factor for substrate selection. The hydroxylation-dependent glycosylation of Pro143 of DdSkp1 by a pentasaccharide forms the basis of a novel O2-sensing mechanism in the social amoeba Dictyostelium and other protists. Previous evidence indicated that glycosylation promotes increased α-helical content correlating with enhanced interaction with three F-box proteins. To localize these differences, we used nuclear magnetic resonance (NMR) methods to compare nonglycosylated DdSkp1 and a glycoform with a single GlcNAc sugar (Gn-DdSkp1). We report NMR assignments of backbone 1HN, 15N, 13Cα, and 13CO nuclei as well as side-chain 13Cß and methyl 13C/1H nuclei of Ile(δ1), Leu, and Val in both unmodified DdSkp1 and Gn-DdSkp1. The random coil index and 15N{1H} HNOE indicate that the C-terminal region, which forms a helix-loop-helix motif centered on Pro143 at the crystallographically defined binding interface with F-box domains, remains dynamic in both DdSkp1 and Gn-DdSkp1. Chemical shifts indicate that the variation of conformation in Gn-DdSkp1, relative to DdSkp1, is limited to this region and characterized by increased helical fold. Extension of the glycan chain results in further changes, also limited to this region. Thus, glycosylation may control F-box protein interactions via a local effect on DdSkp1 conformation, by a mechanism that may be general to many unicellular eukaryotes.

Solution NMR Experiment for Measurement of (15)N-(1)H Residual Dipolar Couplings in Large Proteins and Supramolecular Complexes.

NMR residual dipolar couplings (RDCs) are exquisite probes of protein structure and dynamics. A new solution NMR experiment named 2D SE2 J-TROSY is presented to measure N-H RDCs for proteins and supramolecular complexes in excess of 200 kDa. This enables validation and refinement of their X-ray crystal and solution NMR structures and the characterization of structural and dynamic changes occurring upon complex formation. Accurate N-H RDCs were measured at 750 MHz (1)H resonance frequency for 11-mer 93 kDa (2)H,(15)N-labeled Trp RNA-binding attenuator protein tumbling with a correlation time τc of 120 ns. This is about twice as long as that for the most slowly tumbling system, for which N-H RDCs could be measured, so far, and corresponds to molecular weights of ∼200 kDa at 25 °C. Furthermore, due to the robustness of SE2 J-TROSY with respect to residual (1)H density from exchangeable protons, increased sensitivity at (1)H resonance frequencies around 1 GHz promises to enable N-H RDC measurement for even larger systems.

Solution NMR structures of homeodomains from human proteins ALX4, ZHX1, and CASP8AP2 contribute to the structural coverage of the Human Cancer Protein Interaction Network.

High-quality solution NMR structures of three homeodomains from human proteins ALX4, ZHX1 and CASP8AP2 were solved. These domains were chosen as targets of a biomedical theme project pursued by the Northeast Structural Genomics Consortium. This project focuses on increasing the structural coverage of human proteins associated with cancer.

Solution NMR structures provide first structural coverage of the large protein domain family PF08369 and complementary structural coverage of dark operative protochlorophyllide oxidoreductase complexes.

High-quality NMR structures of the C-terminal domain comprising residues 484-537 of the 537-residue protein Bacterial chlorophyll subunit B (BchB) from Chlorobium tepidum and residues 9-61 of 61-residue Asr4154 from Nostoc sp. (strain PCC 7120) exhibit a mixed α/ß fold comprised of three α-helices and a small ß-sheet packed against second α-helix. These two proteins share 29% sequence similarity and their structures are globally quite similar. The structures of BchB(484-537) and Asr4154(9-61) are the first representative structures for the large protein family (Pfam) PF08369, a family of unknown function currently containing 610 members in bacteria and eukaryotes. Furthermore, BchB(484-537) complements the structural coverage of the dark-operating protochlorophyllide oxidoreductase.

Solution NMR structure of CD1104B from pathogenic Clostridium difficile reveals a distinct α-helical architecture and provides first structural representative of protein domain family PF14203.

A high-quality structure of the 68-residue protein CD1104B from Clostridium difficile strain 630 exhibits a distinct all α-helical fold. The structure presented here is the first representative of bacterial protein domain family PF14203 (currently 180 members) of unknown function (DUF4319) and reveals that the side-chains of the only two strictly conserved residues (Glu 8 and Lys 48) form a salt bridge. Moreover, these two residues are located in the vicinity of the largest surface cleft which is predicted to contribute to a surface area involved in protein-protein interactions. This, along with its coding in transposon CTn4, suggests that CD1104B (and very likely all members of Pfam 14203) functions by interacting with other proteins required for the transfer of transposons between different bacterial species.

Solution NMR structures reveal a distinct architecture and provide first structures for protein domain family PF04536.

The protein family (Pfam) PF04536 is a broadly conserved domain family of unknown function (DUF477), with more than 1,350 members in prokaryotic and eukaryotic proteins. High-quality NMR structures of the N-terminal domain comprising residues 41-180 of the 684-residue protein CG2496 from Corynebacterium glutamicum and the N-terminal domain comprising residues 35-182 of the 435-residue protein PG0361 from Porphyromonas gingivalis both exhibit an α/ß fold comprised of a four-stranded ß-sheet, three α-helices packed against one side of the sheet, and a fourth α-helix attached to the other side. In spite of low sequence similarity (18%) assessed by structure-based sequence alignment, the two structures are globally quite similar. However, moderate structural differences are observed for the relative orientation of two of the four helices. Comparison with known protein structures reveals that the α/ß architecture of CG2496(41-180) and PG0361(35-182) has previously not been characterized. Moreover, calculation of surface charge potential and identification of surface clefts indicate that the two domains very likely have different functions.

Solution NMR structures reveal unique homodimer formation by a winged helix-turn-helix motif and provide first structures for protein domain family PF10771.

High-quality NMR structures of the homo-dimeric proteins Bvu3908 (69-residues in monomeric unit) from Bacteroides vulgatus and Bt2368 (74-residues) from Bacteroides thetaiotaomicron reveal the presence of winged helix-turn-helix (wHTH) motifs mediating tight complex formation. Such homo-dimer formation by winged HTH motifs is otherwise found only in two DNA-binding proteins with known structure: the C-terminal wHTH domain of transcriptional activator FadR from E. coli and protein TubR from B. thurigensis, which is involved in plasmid DNA segregation. However, the relative orientation of the wHTH motifs is different and residues involved in DNA-binding are not conserved in Bvu3908 and Bt2368. Hence, the proteins of the present study are not very likely to bind DNA, but are likely to exhibit a function that has thus far not been ascribed to homo-dimers formed by winged HTH motifs. The structures of Bvu3908 and Bt2368 are the first atomic resolution structures for PFAM family PF10771, a family of unknown function (DUF2582) currently containing 128 members.

Solution NMR structure of yeast succinate dehydrogenase flavinylation factor Sdh5 reveals a putative Sdh1 binding site.

The yeast mitochondrial protein Sdh5 is required for the covalent attachment of flavin adenine dinucleotide (FAD) to protein Sdh1, a subunit of the heterotetrameric enzyme succinate dehydrogenase. The NMR structure of Sdh5 represents the first eukaryotic structure of Pfam family PF03937 and reveals a conserved surface region, which likely represents a putative Sdh1-Sdh5 interaction interface. Point mutations in this region result in the loss of covalent flavinylation of Sdh1. Moreover, chemical shift perturbation measurements showed that Sdh5 does not bind FAD in vitro, indicating that it is not a simple cofactor transporter in vivo.

Highly precise measurement of kinetic isotope effects using 1H-detected 2D [13C,1H]-HSQC NMR spectroscopy.

A new method is presented for measuring kinetic isotope effects (KIEs) by (1)H-detected 2D [(13)C,(1)H]-heteronuclear single quantum coherence (HSQC) NMR spectroscopy. The high accuracy of this approach was exemplified for the reaction catalyzed by glucose-6-phosphate dehydrogenase by comparing the 1-(13)C KIE with the published value obtained using isotope ratio mass spectrometry. High precision was demonstrated for the reaction catalyzed by 1-deoxy-D-xylulose-5-phosphate reductoisomerase from Mycobacterium tuberculosis. 2-, 3-, and 4-(13)C KIEs were found to be 1.0031(4), 1.0303(12), and 1.0148(2), respectively. These KIEs provide evidence for a cleanly rate-limiting retroaldol step during isomerization. The high intrinsic sensitivity and signal dispersion of 2D [(13)C,(1)H]-HSQC offer new avenues to study challenging systems where low substrate concentration and/or signal overlap impedes 1D (13)C NMR data acquisition. Moreover, this approach can take advantage of highest-field spectrometers, which are commonly equipped for (1)H detection with cryogenic probes.

Glycan Conformation in the Heavily Glycosylated Protein, CEACAM1.

Glycans attached to glycoproteins can contribute to stability, mediate interactions with other proteins, and initiate signal transduction. Glycan conformation, which is critical to these processes, is highly variable and often depicted as sampling a multitude of conformers. These conformers can be generated by molecular dynamics simulations, and more inclusively by accelerated molecular dynamics, as well as other extended sampling methods. However, experimental assessments of the contribution that various conformers make to a native ensemble are rare. Here, we use long-range pseudo-contact shifts (PCSs) of NMR resonances from an isotopically labeled glycoprotein to identify preferred conformations of its glycans. The N-terminal domain from human Carcinoembryonic Antigen Cell Adhesion Molecule 1, hCEACAM1-Ig1, was used as the model glycoprotein in this study. It has been engineered to include a lanthanide-ion-binding loop that generates PCSs, as well as a homogeneous set of three 13C-labeled N-glycans. Analysis of the PCSs indicates that preferred glycan conformers have extensive contacts with the protein surface. Factors leading to this preference appear to include interactions between N-acetyl methyls of GlcNAc residues and hydrophobic surface pockets on the protein surface.

AssignSLP_GUI, a software tool exploiting AI for NMR resonance assignment of sparsely labeled proteins.

Not all proteins are amenable to uniform isotopic labeling with 13C and 15N, something needed for the widely used, and largely deductive, triple resonance assignment process. Among them are proteins expressed in mammalian cell culture where native glycosylation can be maintained, and proper formation of disulfide bonds facilitated. Uniform labeling in mammalian cells is prohibitively expensive, but sparse labeling with one or a few isotopically enriched amino acid types is an option for these proteins. However, assignment then relies on accessing the best match between a variety of measured NMR parameters and predictions based on 3D structure, often from X-ray crystallography. Finding this match is a challenging process that has benefitted from many computational tools, including trained neural nets for chemical shift prediction, genetic algorithms for searches through a myriad of assignment possibilities, and now AI-based prediction of high-quality structures for protein targets. AssignSLP_GUI, a new version of a software package for assignment of resonances from sparsely-labeled proteins, uses many of these tools. These tools and new additions to the package are highlighted in an application to a sparsely-labeled domain from a glycoprotein, CEACAM1.

Validated determination of NRG1 Ig-like domain structure by mass spectrometry coupled with computational modeling.

High resolution hydroxyl radical protein footprinting (HR-HRPF) is a mass spectrometry-based method that measures the solvent exposure of multiple amino acids in a single experiment, offering constraints for experimentally informed computational modeling. HR-HRPF-based modeling has previously been used to accurately model the structure of proteins of known structure, but the technique has never been used to determine the structure of a protein of unknown structure. Here, we present the use of HR-HRPF-based modeling to determine the structure of the Ig-like domain of NRG1, a protein with no close homolog of known structure. Independent determination of the protein structure by both HR-HRPF-based modeling and heteronuclear NMR was carried out, with results compared only after both processes were complete. The HR-HRPF-based model was highly similar to the lowest energy NMR model, with a backbone RMSD of 1.6 Å. To our knowledge, this is the first use of HR-HRPF-based modeling to determine a previously uncharacterized protein structure.

Solution NMR structure of MED25(391-543) comprising the activator-interacting domain (ACID) of human mediator subunit 25.

The solution NMR structure of protein MED25(391-543), comprising the activator interacting domain (ACID) of subunit 25 of the human mediator, is presented along with the measurement of polypeptide backbone heteronuclear 15N-{1H} NOEs to identify fast internal motional modes. This domain interacts with the acidic transactivation domains of Herpes simplex type 1 (HSV-1) protein VP16 and the Varicella-zoster virus (VZV) major transactivator protein IE62, which initiate transcription of viral genes. The structure is similar to the ß-barrel domains of the human protein Ku and the SPOC domain of human protein SHARP, and provides a starting point to understand the structural biology of initiation of HSV-1 and VZV gene activation. Homology models built for the two ACID domains of the prostate tumor overexpressed (PTOV1) protein using the structure of MED25(391-543) as a template suggest that differential biological activities of the ACID domains in MED25 and PTOV1 arise from modulation of quite similar protein-protein interactions by variable residues grouped around highly conserved charged surface areas.

Consistent blind protein structure generation from NMR chemical shift data.

Protein NMR chemical shifts are highly sensitive to local structure. A robust protocol is described that exploits this relation for de novo protein structure generation, using as input experimental parameters the (13)C(alpha), (13)C(beta), (13)C', (15)N, (1)H(alpha) and (1)H(N) NMR chemical shifts. These shifts are generally available at the early stage of the traditional NMR structure determination process, before the collection and analysis of structural restraints. The chemical shift based structure determination protocol uses an empirically optimized procedure to select protein fragments from the Protein Data Bank, in conjunction with the standard ROSETTA Monte Carlo assembly and relaxation methods. Evaluation of 16 proteins, varying in size from 56 to 129 residues, yielded full-atom models that have 0.7-1.8 A root mean square deviations for the backbone atoms relative to the experimentally determined x-ray or NMR structures. The strategy also has been successfully applied in a blind manner to nine protein targets with molecular masses up to 15.4 kDa, whose conventional NMR structure determination was conducted in parallel by the Northeast Structural Genomics Consortium. This protocol potentially provides a new direction for high-throughput NMR structure determination.

Analysis of the varicella-zoster virus IE62 N-terminal acidic transactivating domain and its interaction with the human mediator complex.

The varicella-zoster virus major transactivator, IE62, contains a potent N-terminal acidic transcriptional activation domain (TAD). Our experiments revealed that the minimal IE62 TAD encompasses amino acids (aa) 19 to 67. We showed that the minimal TAD interacts with the human Mediator complex. Site-specific mutations revealed residues throughout the minimal TAD that are important for both activation and Mediator interaction. The TAD interacts directly with aa 402 to 590 of the MED25 subunit, and site-specific TAD mutations abolished this interaction. Two-dimensional nuclear magnetic resonance spectroscopy revealed that the TAD is intrinsically unstructured. Our studies suggest that transactivation may involve the TAD adopting a defined structure upon binding MED25.

NMR and X-RAY structures of human E2-like ubiquitin-fold modifier conjugating enzyme 1 (UFC1) reveal structural and functional conservation in the metazoan UFM1-UBA5-UFC1 ubiquination pathway.

For cell regulation, E2-like ubiquitin-fold modifier conjugating enzyme 1 (Ufc1) is involved in the transfer of ubiquitin-fold modifier 1 (Ufm1), a ubiquitin like protein which is activated by E1-like enzyme Uba5, to various target proteins. Thereby, Ufc1 participates in the very recently discovered Ufm1-Uba5-Ufc1 ubiquination pathway which is found in metazoan organisms. The structure of human Ufc1 was solved by using both NMR spectroscopy and X-ray crystallography. The complementary insights obtained with the two techniques provided a unique basis for understanding the function of Ufc1 at atomic resolution. The Ufc1 structure consists of the catalytic core domain conserved in all E2-like enzymes and an additional N-terminal helix. The active site Cys(116), which forms a thio-ester bond with Ufm1, is located in a flexible loop that is highly solvent accessible. Based on the Ufc1 and Ufm1 NMR structures, a model could be derived for the Ufc1-Ufm1 complex in which the C-terminal Gly(83) of Ufm1 may well form the expected thio-ester with Cys(116), suggesting that Ufm1-Ufc1 functions as described for other E1-E2-E3 machineries. alpha-helix 1 of Ufc1 adopts different conformations in the crystal and in solution, suggesting that this helix plays a key role to mediate specificity.