A novel strategy for NMR resonance assignment and protein structure determination

Large-scale discovery, analysis, and design of protein energy landscapes

All folded proteins continuously fluctuate between their low-energy native structures and higher energy conformations that can be partially or fully unfolded. These rare states influence protein function, interactions, aggregation, and immunogenicity, yet they remain far less understood than protein native states. Although native protein structures are now often predictable with impressive accuracy, conformational fluctuations and their energies remain largely invisible and unpredictable, and experimental challenges have prevented large-scale measurements that could improve machine learning and physics-based modeling. Here, we introduce a multiplexed experimental approach to analyze the energies of conformational fluctuations for hundreds of protein domains in parallel using intact protein hydrogen-deuterium exchange mass spectrometry. We analyzed 5,778 domains 28-64 amino acids in length, revealing hidden variation in conformational fluctuations even between sequences sharing the same fold and global folding stability. Site-resolved hydrogen exchange NMR analysis of 13 domains showed that these fluctuations often involve entire secondary structural elements with lower stability than the overall fold. Computational modeling of our domains identified structural features that correlated with the experimentally observed fluctuations, enabling us to design mutations that stabilized low-stability structural segments. Our dataset enables new machine learning-based analysis of protein energy landscapes, and our experimental approach promises to reveal these landscapes at unprecedented scale.

Heat-induced structural and chemical changes to a computationally designed miniprotein

The de novo design of miniprotein inhibitors has recently emerged as a new technology to create proteins that bind with high affinity to specific therapeutic targets. Their size, ease of expression, and apparent high stability makes them excellent candidates for a new class of protein drugs. However, beyond circular dichroism melts and hydrogen/deuterium exchange experiments, little is known about their dynamics, especially at the elevated temperatures they seemingly tolerate quite well. To address that and gain insight for future designs, we have focused on identifying unintended and previously overlooked heat-induced structural and chemical changes in a particularly stable model miniprotein, EHEE_rd2_0005. Nuclear magnetic resonance (NMR) studies suggest the presence of dynamics on multiple time and temperature scales. Transiently elevating the temperature results in spontaneous chemical deamidation visible in the NMR spectra, which we validate using both capillary electrophoresis and mass spectrometry (MS) experiments. High temperatures also result in greatly accelerated intrinsic rates of hydrogen exchange and signal loss in NMR heteronuclear single quantum coherence spectra from local unfolding. These losses are in excellent agreement with both room temperature hydrogen exchange experiments and hydrogen bond disruption in replica exchange molecular dynamics simulations. Our analysis reveals important principles for future miniprotein designs and the potential for high stability to result in long-lived alternate conformational states.

The RavA-ViaA chaperone complex modulates bacterial persistence through its association with the fumarate reductase enzyme

Regulatory ATPase variant A (RavA) is a MoxR AAA+ protein that functions together with a partner protein termed von Willebrand factor type A interacting with AAA+ ATPase (ViaA). RavA-ViaA are functionally associated with anaerobic respiration in Escherichia coli through interactions with the fumarate reductase (Frd) electron transport complex. Through this association, RavA and ViaA modulate the activity of the Frd complex and, hence, are proposed to have chaperone-like activity. However, the functional role of RavA-ViaA in the cell is not yet well established. We had demonstrated that RavA-ViaA can sensitize E. coli cells to sublethal concentrations of the aminoglycoside class of antibiotics. Since Frd has been associated with bacterial persistence against antibiotics, the relationship of RavA-ViaA and Frd was explored within this context. Experiments performed here reveal a function of RavA-ViaA in bacterial persistence upon treatment with antibiotics through the association of the chaperone complex with Frd. As part of this work, the NMR structure of the N-terminal domain of ViaA was solved. The structure reveals a novel alpha helical fold, which we name the VAN fold, that has not been observed before. We show that this domain is required for the function of the chaperone complex. We propose that modulating the levels of RavA-ViaA could enhance the susceptibility of Gram-negative bacteria to antibiotics.

De novo design of immunoglobulin-like domains

Antibodies, and antibody derivatives such as nanobodies, contain immunoglobulin-like (Ig) β-sandwich scaffolds which anchor the hypervariable antigen-binding loops and constitute the largest growing class of drugs. Current engineering strategies for this class of compounds rely on naturally existing Ig frameworks, which can be hard to modify and have limitations in manufacturability, designability and range of action. Here, we develop design rules for the central feature of the Ig fold architecture—the non-local cross-β structure connecting the two β-sheets—and use these to design highly stable Ig domains de novo, confirm their structures through X-ray crystallography, and show they can correctly scaffold functional loops. Our approach opens the door to the design of antibody-like scaffolds with tailored structures and superior biophysical properties.

The immunoglobulin domain framework of antibodies has been a long standing design challenge. Here, the authors describe design rules for tailoring these domains and show they can be accurately designed, de novo, with high stability and the ability to scaffold functional loops.

The MYC oncoprotein directly interacts with its chromatin cofactor PNUTS to recruit PP1 phosphatase

Despite MYC dysregulation in most human cancers, strategies to target this potent oncogenic driver remain an urgent unmet need. Recent evidence shows the PP1 phosphatase and its regulatory subunit PNUTS control MYC phosphorylation, chromatin occupancy, and stability, however the molecular basis remains unclear. Here we demonstrate that MYC interacts directly with PNUTS through the MYC homology Box 0 (MB0), a highly conserved region recently shown to be important for MYC oncogenic activity. By NMR we identified a distinct peptide motif within MB0 that interacts with PNUTS residues 1–148, a functional unit, here termed PNUTS amino-terminal domain (PAD). Using NMR spectroscopy we determined the solution structure of PAD, and characterised its MYC-binding patch. Point mutations of residues at the MYC-PNUTS interface significantly weaken their interaction both in vitro and in vivo, leading to elevated MYC phosphorylation. These data demonstrate that the MB0 region of MYC directly interacts with the PAD of PNUTS, which provides new insight into the control mechanisms of MYC as a regulator of gene transcription and a pervasive cancer driver.

Alternative splicing and allosteric regulation modulate the chromatin binding of UHRF1

UHRF1 is an important epigenetic regulator associated with apoptosis and tumour development. It is a multidomain protein that integrates readout of different histone modification states and DNA methylation with enzymatic histone ubiquitylation activity. Emerging evidence indicates that the chromatin-binding and enzymatic modules of UHRF1 do not act in isolation but interplay in a coordinated and regulated manner. Here, we compared two splicing variants (V1, V2) of murine UHRF1 (mUHRF1) with human UHRF1 (hUHRF1). We show that insertion of nine amino acids in a linker region connecting the different TTD and PHD histone modification-binding domains causes distinct H3K9me3-binding behaviour of mUHRF1 V1. Structural analysis suggests that in mUHRF1 V1, in contrast to V2 and hUHRF1, the linker is anchored in a surface groove of the TTD domain, resulting in creation of a coupled TTD-PHD module. This establishes multivalent, synergistic H3-tail binding causing distinct cellular localization and enhanced H3K9me3-nucleosome ubiquitylation activity. In contrast to hUHRF1, H3K9me3-binding of the murine proteins is not allosterically regulated by phosphatidylinositol 5-phosphate that interacts with a separate less-conserved polybasic linker region of the protein. Our results highlight the importance of flexible linkers in regulating multidomain chromatin binding proteins and point to divergent evolution of their regulation.

Structure and physiological function of the human KCNQ1 channel voltage sensor intermediate state

Voltage-gated ion channels feature voltage sensor domains (VSDs) that exist in three distinct conformations during activation: resting, intermediate, and activated. Experimental determination of the structure of a potassium channel VSD in the intermediate state has previously proven elusive. Here, we report and validate the experimental three-dimensional structure of the human KCNQ1 voltage-gated potassium channel VSD in the intermediate state. We also used mutagenesis and electrophysiology in Xenopus laevisoocytes to functionally map the determinants of S4 helix motion during voltage-dependent transition from the intermediate to the activated state. Finally, the physiological relevance of the intermediate state KCNQ1 conductance is demonstrated using voltage-clamp fluorometry. This work illuminates the structure of the VSD intermediate state and demonstrates that intermediate state conductivity contributes to the unusual versatility of KCNQ1, which can function either as the slow delayed rectifier current (IKs) of the cardiac action potential or as a constitutively active epithelial leak current.

The MLL1 trimeric catalytic complex is a dynamic conformational ensemble stabilized by multiple weak interactions

Histone H3K4 methylation is an epigenetic mark associated with actively transcribed genes. This modification is catalyzed by the mixed lineage leukaemia (MLL) family of histone methyltransferases including MLL1, MLL2, MLL3, MLL4, SET1A and SET1B. The catalytic activity of this family is dependent on interactions with additional conserved proteins, but the structural basis for subunit assembly and the mechanism of regulation is not well understood. We used a hybrid methods approach to study the assembly and biochemical function of the minimally active MLL1 complex (MLL1, WDR5 and RbBP5). A combination of small angle X-ray scattering, cross-linking mass spectrometry, nuclear magnetic resonance spectroscopy and computational modeling were used to generate a dynamic ensemble model in which subunits are assembled via multiple weak interaction sites. We identified a new interaction site between the MLL1 SET domain and the WD40 β-propeller domain of RbBP5, and demonstrate the susceptibility of the catalytic function of the complex to disruption of individual interaction sites.

NMR Structural Analysis of Isolated Shaker Voltage-Sensing Domain in LPPG Micelles

The voltage-sensing domain (VSD) is a conserved structural module that regulates the gating of voltage-dependent ion channels in response to a change in membrane potential. Although the structures of many VSD-containing ion channels are now available, our understanding of the structural dynamics associated with gating transitions remains limited. To probe dynamics with site-specific resolution, we utilized NMR spectroscopy to characterize the VSD derived from Shaker potassium channel in 1-palmitoyl-2-hydroxy-sn-glycero-3-phospho-(1'-rac-glycerol) (LPPG) micelles. The backbone dihedral angles predicted based on secondary chemical shifts using torsion angle likeliness obtained from shift (TALOS+) showed that the Shaker-VSD shares many structural features with the homologous Kv1.2/2.1 chimera, including a transition from α-helix to 3₁₀ helix in the C-terminal portion of the fourth transmembrane helix. Nevertheless, there are clear differences between the Shaker-VSD and Kv1.2/2.1 chimera in the S2-S3 linker and S3 transmembrane region, where the organization of secondary structure elements in Shaker-VSD appears to more closely resemble the KvAP-VSD. Comparison of microsecond-long molecular dynamics simulations of Kv 1.2-VSD in LPPG micelles and a 1-palmitoyl-2-oleoyl-glycero-3-phosphocholine (POPC) bilayer showed that LPPG micelles do not induce significant structural distortion in the isolated voltage sensor. To assess the integrity of the tertiary fold, we directly probed the binding of BrMT analog 2-[2-({[3-(2-amino-ethyl)-6-bromo-1H-indol-2-yl]methoxy}k7methyl)-6-bromo-1H-indol-3-yl]ethan-1-amine (BrET), a gating modifier toxin, and identified the location of the putative binding site. Our results suggest that the Shaker-VSD in LPPG micelles is in a native-like fold and is likely to provide valuable insights into the dynamics of voltage-gating and its regulation.

Global analysis of protein folding using massively parallel design, synthesis, and testing

Proteins fold into unique native structures stabilized by thousands of weak interactions that collectively overcome the entropic cost of folding. Although these forces are "encoded" in the thousands of known protein structures, "decoding" them is challenging because of the complexity of natural proteins that have evolved for function, not stability. We combined computational protein design, next-generation gene synthesis, and a high-throughput protease susceptibility assay to measure folding and stability for more than 15,000 de novo designed miniproteins, 1000 natural proteins, 10,000 point mutants, and 30,000 negative control sequences. This analysis identified more than 2500 stable designed proteins in four basic folds-a number sufficient to enable us to systematically examine how sequence determines folding and stability in uncharted protein space. Iteration between design and experiment increased the design success rate from 6% to 47%, produced stable proteins unlike those found in nature for topologies where design was initially unsuccessful, and revealed subtle contributions to stability as designs became increasingly optimized. Our approach achieves the long-standing goal of a tight feedback cycle between computation and experiment and has the potential to transform computational protein design into a data-driven science.

Characterization of DNA Binding by the Isolated N-Terminal Domain of Vaccinia Virus DNA Topoisomerase IB

Vaccinia TopIB (vTopIB), a 314-amino acid eukaryal-type IB topoisomerase, recognizes and transesterifies at the DNA sequence 5'-(T/C)CCTT↓, leading to the formation of a covalent DNA-(3'-phosphotyrosyl²⁷⁴)-enzyme intermediate in the supercoil relaxation reaction. The C-terminal segment of vTopIB (amino acids 81-314), which engages the DNA minor groove at the scissile phosphodiester, comprises an autonomous catalytic domain that retains cleavage specificity, albeit with a cleavage site affinity lower than that of the full-length enzyme. The N-terminal domain (amino acids 1-80) engages the major groove on the DNA face opposite the scissile phosphodiester. Whereas DNA contacts of the N-terminal domain have been implicated in the DNA site affinity of vTopIB, it was not known whether the N-terminal domain per se could bind DNA. Here, using isothermal titration calorimetry, we demonstrate the ability of the isolated N-terminal domain to bind a CCCTT-containing 24-mer duplex with an apparent affinity that is ∼2.2-fold higher than that for an otherwise identical duplex in which the pentapyrimidine sequence is changed to ACGTG. Analyses of the interactions of the isolated N-terminal domain with duplex DNA via solution nuclear magnetic resonance methods are consistent with its DNA contacts observed in DNA-bound crystal structures of full-length vTopIB. The chemical shift perturbations and changes in hydrodynamic properties triggered by CCCTT DNA versus non-CCCTT DNA suggest differences in DNA binding dynamics. The importance of key N-terminal domain contacts in the context of full-length vTopIB is underscored by assessing the effects of double-alanine mutations on DNA transesterification and its sensitivity to ionic strength.

Staphylokinase has distinct modes of interaction with antimicrobial peptides, modulating its plasminogen-activation properties

Staphylokinase (Sak) is a plasminogen activator protein that is secreted by many Staphylococcus aureus strains. Sak also offers protection by binding and inhibiting specific antimicrobial peptides (AMPs). Here, we evaluate Sak as a more general interaction partner for AMPs. Studies with melittin, mCRAMP, tritrpticin and bovine lactoferricin indicate that the truncation of the first ten residues of Sak (SakΔN10), which occurs in vivo and uncovers important residues in a bulge region, improves its affinity for AMPs. Melittin and mCRAMP have a lower affinity for SakΔN10, and in docking studies, they bind to the N-terminal segment and bulge region of SakΔN10. By comparison, lactoferricin and tritrpticin form moderately high affinity 1:1 complexes with SakΔN10 and their cationic residues form several electrostatic interactions with the protein's α-helix. Overall, our work identifies two distinct AMP binding surfaces on SakΔN10 whose occupation would lead to either inhibition or promotion of its plasminogen activating properties.

Structure of the Single-lobe Myosin Light Chain C in Complex with the Light Chain-binding Domains of Myosin-1C Provides Insights into Divergent IQ Motif Recognition

Myosin light chains are key regulators of class 1 myosins and typically comprise two domains, with calmodulin being the archetypal example. They bind IQ motifs within the myosin neck region and amplify conformational changes in the motor domain. A single lobe light chain, myosin light chain C (MlcC), was recently identified and shown to specifically bind to two sequentially divergent IQ motifs of the Dictyostelium myosin-1C. To provide a molecular basis of this interaction, the structures of apo-MlcC and a 2:1 MlcC·myosin-1C neck complex were determined. The two non-functional EF-hand motifs of MlcC pack together to form a globular four-helix bundle that opens up to expose a central hydrophobic groove, which interacts with the N-terminal portion of the divergent IQ1 and IQ2 motifs. The N- and C-terminal regions of MlcC make critical contacts that contribute to its specific interactions with the myosin-1C divergent IQ motifs, which are contacts that deviate from the traditional mode of calmodulin-IQ recognition.

The molecular mechanism of nuclear transport revealed by atomic-scale measurements

Nuclear pore complexes (NPCs) form a selective filter that allows the rapid passage of transport factors (TFs) and their cargoes across the nuclear envelope, while blocking the passage of other macromolecules. Intrinsically disordered proteins (IDPs) containing phenylalanyl-glycyl (FG)-rich repeats line the pore and interact with TFs. However, the reason that transport can be both fast and specific remains undetermined, through lack of atomic-scale information on the behavior of FGs and their interaction with TFs. We used nuclear magnetic resonance spectroscopy to address these issues. We show that FG repeats are highly dynamic IDPs, stabilized by the cellular environment. Fast transport of TFs is supported because the rapid motion of FG motifs allows them to exchange on and off TFs extremely quickly through transient interactions. Because TFs uniquely carry multiple pockets for FG repeats, only they can form the many frequent interactions needed for specific passage between FG repeats to cross the NPC.

DOI:http://dx.doi.org/10.7554/eLife.10027.001

Eukaryotic cells have a nucleus that contains most of the organism's genetic material. Two layers of membrane form an envelope around the nucleus and protect its contents from the rest of the cell's interior. However, this protective barrier must also allow certain proteins and nucleic acids(collectively called ‘cargo’) to move in and out of the nucleus.

Cargo molecules can pass through channel-like structures called nuclear pore complexes, which are embedded in the nuclear envelope. However, transport across this barrier is highly selective. While small molecules can pass freely through nuclear pore complexes, larger cargo can only be transported when they are bound to so-called transport factors. The nuclear pore complex is a large structure made up of more than 30 different proteins called nucleoporins. Like all proteins, nucleoporins are built from amino acids. Many nucleoporins contain repeating units of two amino acids, namely phenylalanine (which is often referred to as ‘F’) and glycine (or ‘G’). These ‘FG nucleoporins’ are found on the inside of the nuclear pore complex and interact with transport factors to allow them to transit across the nuclear envelope.

Several models have been put forward to explain how FG nucleoporins block the passage of most molecules. But it was unclear from these models how these nucleoporins could do this while simultaneously allowing the selective and fast transport of nuclear transport receptors. There was also a major lack of experimental data that probed the behavior of FG nucleoporins in detail.

Hough, Dutta et al. have now used a technique called nuclear magnetic resonance spectroscopy (or NMR for short) to address this issue. NMR can be used to analyze the structure of proteins and how they interact with other molecules. This analysis revealed that FG nucleoporins never adopt an ordered three-dimensional shape, even briefly; instead they remain unfolded or disordered, moving constantly. Nevertheless, and unlike many other unfolded proteins, FG nucleoporins do not aggregate into clumps. This is because they are constantly changing and continuously interacting with other molecules present inside the cell, which prevents them from aggregating.

Hough, Dutta et al. also observed that the repeating units in the FG nucleoporins engaged briefly with a large number of sites or pockets present on the transport factors. These FG repeats can bind and then release the transport factors at unusually high speeds, which enables the transport factors to move quickly through the nuclear pore complex. This transit is specific because only transport factors have a high capacity for interacting with the FG repeats. These findings provide an explanation for how the nuclear pore complex achieves fast and selective transport. Further work is needed to see whether certain FG nucleoporins specifically interact with a particular type of transport factor, to provide preferred transport routes through the nuclear pore complex.

DOI:http://dx.doi.org/10.7554/eLife.10027.002

The second round of Critical Assessment of Automated Structure Determination of Proteins by NMR: CASD-NMR-2013

The second round of the community-wide initiative Critical Assessment of automated Structure Determination of Proteins by NMR (CASD-NMR-2013) comprised ten blind target datasets, consisting of unprocessed spectral data, assigned chemical shift lists and unassigned NOESY peak and RDC lists, that were made available in both curated (i.e. manually refined) or un-curated (i.e. automatically generated) form. Ten structure calculation programs, using fully automated protocols only, generated a total of 164 three-dimensional structures (entries) for the ten targets, sometimes using both curated and un-curated lists to generate multiple entries for a single target. The accuracy of the entries could be established by comparing them to the corresponding manually solved structure of each target, which was not available at the time the data were provided. Across the entire data set, 71 % of all entries submitted achieved an accuracy relative to the reference NMR structure better than 1.5 Å. Methods based on NOESY peak lists achieved even better results with up to 100% of the entries within the 1.5 Å threshold for some programs. However, some methods did not converge for some targets using un-curated NOESY peak lists. Over 90% of the entries achieved an accuracy better than the more relaxed threshold of 2.5 Å that was used in the previous CASD-NMR-2010 round. Comparisons between entries generated with un-curated versus curated peaks show only marginal improvements for the latter in those cases where both calculations converged.

KCMF1 (potassium channel modulatory factor 1) Links RAD6 to UBR4 (ubiquitin N-recognin domain-containing E3 ligase 4) and lysosome-mediated degradation

RAD6 is a ubiquitin E2 protein with roles in a number of different biological processes. Here, using affinity purification coupled with mass spectrometry, we identify a number of new RAD6 binding partners, including the poorly characterized ubiquitin E3 ligases KCMF1 (potassium channel modulatory factor 1) and UBR4 (ubiquitin N-recognin domain-containing E3 ligase 4), a protein that can bind N-end rule substrates, and which was recently linked to lysosome-mediated degradation and autophagy. NMR, combined with in vivo and in vitro interaction mapping, demonstrate that the KCMF1 C terminus binds directly to RAD6, whereas N-terminal domains interact with UBR4 and other intracellular vesicle- and mitochondria-associated proteins. KCMF1 and RAD6 colocalize at late endosomes and lysosomes, and cells disrupted for KCMF1 or RAD6 function display defects in late endosome vesicle dynamics. Notably, we also find that two different RAD6A point mutants (R7W and R11Q) found in X-linked intellectual disability (XLID) patients specifically lose the interaction with KCMF1 and UBR4, but not with other previously identified RAD6 interactors. We propose that RAD6-KCMF1-UBR4 represents a unique new E2-E3 complex that targets unknown N-end rule substrates for lysosome-mediated degradation, and that disruption of this complex via RAD6A mutations could negatively affect neuronal function in XLID patients.

Mechanisms mediating enhanced neutralization efficacy of staphylococcal enterotoxin B by combinations of monoclonal antibodies

Staphylococcal enterotoxin B (SEB) is a superantigen that cross-links the major histocompatibility complex class II and specific V-β chains of the T-cell receptor, thus forming a ternary complex. Developing neutralizing mAb to disrupt the ternary complex and abrogate the resulting toxicity is a major therapeutic challenge because SEB is effective at very low concentrations. We show that combining two SEB-specific mAbs enhances their efficacy, even though one of the two mAbs by itself has no effect on neutralization. Crystallography was employed for fine-mapping conformational epitopes in binary and ternary complexes between SEB and Fab fragments. NMR spectroscopy was used to validate and identify subtle allosteric changes induced by mAbs binding to SEB. The mapping of epitopes established that a combination of different mAbs can enhance efficacy of mAb-mediated protection from SEB induced lethal shock by two different mechanisms: one mAb mixture promoted clearance of the toxin both in vitro and in vivo by FcR-mediated cross-linking and clearance, whereas the other mAb mixture induced subtle allosteric conformational changes in SEB that perturbed formation of the SEB·T-cell receptor·major histocompatibility complex class II trimer. Finally structural information accurately predicted mAb binding to other superantigens that share conformational epitopes with SEB. Fine mapping of conformational epitopes is a powerful tool to establish the mechanism and optimize the action of synergistic mAb combinations.

Structural characterization of a flexible two-domain protein in solution using small angle X-ray scattering and NMR data

Multidomain proteins in which individual domains are connected by linkers often possess inherent interdomain flexibility that significantly complicates their structural characterization in solution using either nuclear magnetic resonance (NMR) spectroscopy or small-angle X-ray scattering (SAXS) alone. Here, we report a protocol for joint refinement of flexible multidomain protein structures against NMR distance and angular restraints, residual dipolar couplings, and SAXS data. The protocol is based on the ensemble optimization method principle (Bernadó et al., 2007) and is compared with different refinement strategies for the structural characterization of the flexible two-domain protein sf3636 from Shigella flexneri 2a. The results of our refinement suggest the existence of a dominant population of configurational states in solution possessing an overall elongated shape and restricted relative twisting of the two domains.

Structural and functional characterization of DUF1471 domains of Salmonella proteins SrfN, YdgH/SssB, and YahO

Bacterial species in the Enterobacteriaceae typically contain multiple paralogues of a small domain of unknown function (DUF1471) from a family of conserved proteins also known as YhcN or BhsA/McbA. Proteins containing DUF1471 may have a single or three copies of this domain. Representatives of this family have been demonstrated to play roles in several cellular processes including stress response, biofilm formation, and pathogenesis. We have conducted NMR and X-ray crystallographic studies of four DUF1471 domains from Salmonella representing three different paralogous DUF1471 subfamilies: SrfN, YahO, and SssB/YdgH (two of its three DUF1471 domains: the N-terminal domain I (residues 21–91), and the C-terminal domain III (residues 244–314)). Notably, SrfN has been shown to have a role in intracellular infection by Salmonella Typhimurium. These domains share less than 35% pairwise sequence identity. Structures of all four domains show a mixed α+β fold that is most similar to that of bacterial lipoprotein RcsF. However, all four DUF1471 sequences lack the redox sensitive cysteine residues essential for RcsF activity in a phospho-relay pathway, suggesting that DUF1471 domains perform a different function(s). SrfN forms a dimer in contrast to YahO and SssB domains I and III, which are monomers in solution. A putative binding site for oxyanions such as phosphate and sulfate was identified in SrfN, and an interaction between the SrfN dimer and sulfated polysaccharides was demonstrated, suggesting a direct role for this DUF1471 domain at the host-pathogen interface.

An overview of tools for the validation of protein NMR structures

Biomolecular structures at atomic resolution present a valuable resource for the understanding of biology. NMR spectroscopy accounts for 11% of all structures in the PDB repository. In response to serious problems with the accuracy of some of the NMR-derived structures and in order to facilitate proper analysis of the experimental models, a number of program suites are available. We discuss nine of these tools in this review: PROCHECK-NMR, PSVS, GLM-RMSD, CING, Molprobity, Vivaldi, ResProx, NMR constraints analyzer and QMEAN. We evaluate these programs for their ability to assess the structural quality, restraints and their violations, chemical shifts, peaks and the handling of multi-model NMR ensembles. We document both the input required by the programs and output they generate. To discuss their relative merits we have applied the tools to two representative examples from the PDB: a small, globular monomeric protein (Staphylococcal nuclease from S. aureus, PDB entry 2kq3) and a small, symmetric homodimeric protein (a region of human myosin-X, PDB entry 2lw9).

Basic Tilted Helix Bundle - a new protein fold in human FKBP25/FKBP3 and HectD1

In this paper, we describe the structure of a N-terminal domain motif in nuclear-localized FKBP251-73, a member of the FKBP family, together with the structure of a sequence-related subdomain of the E3 ubiquitin ligase HectD1 that we show belongs to the same fold. This motif adopts a compact 5-helix bundle which we name the Basic Tilted Helix Bundle (BTHB) domain. A positively charged surface patch, structurally centered around the tilted helix H4, is present in both FKBP25 and HectD1 and is conserved in both proteins, suggesting a conserved functional role. We provide detailed comparative analysis of the structures of the two proteins and their sequence similarities, and analysis of the interaction of the proposed FKBP25 binding protein YY1. We suggest that the basic motif in BTHB is involved in the observed DNA binding of FKBP25, and that the function of this domain can be affected by regulatory YY1 binding and/or interactions with adjacent domains.

Time-resolved multidimensional NMR with non-uniform sampling

Time-resolved experiments demand high resolution both in spectral dimensions and in time of the studied kinetic process. The latter requirement traditionally prohibits applications of the multidimensional experiments, which, although capable of providing invaluable information about structure and dynamics and almost unlimited spectral resolution, require too lengthy data collection. Our work shows that the problem has a solution in using modern methods of NMR data collection and signal processing. A continuous fast pulsing three-dimensional experiment is acquired using non-uniform sampling during full time of the studied reaction. High sensitivity and time-resolution of a few minutes is achieved by simultaneous processing of the full data set with the multi-dimensional decomposition. The method is verified and illustrated in realistic simulations and by measuring deuterium exchange rates of amide protons in ubiquitin. We applied the method for characterizing kinetics of in vitro phosphorylation of two tyrosine residues in an intrinsically disordered cytosolic domain of the B cell receptor protein CD79b. Signals of many residues including tyrosines in both phosphorylated and unmodified forms of CD79b are found in a heavily crowded region of 2D ¹H-¹⁵N correlation spectrum and the significantly enhanced spectral resolution provided by the 3D time-resolved approach was essential for the quantitative site-specific analysis.

Screening proteins for NMR suitability

NMR spectroscopy is a valuable tool in structural genomics. Identification of protein samples that are amenable to structure determination by NMR spectroscopy requires efficient screening. The preparation of multiple samples in parallel and screening by NMR is described. The method described is applicable to large structural genomics projects but can easily be scaled down for application to small structural biology projects. All the equipment used is commonly found in any NMR structural biology laboratory.

Solution NMR structure and histone binding of the PHD domain of human MLL5

Mixed Lineage Leukemia 5 (MLL5) is a histone methyltransferase that plays a key role in hematopoiesis, spermatogenesis and cell cycle progression. In addition to its catalytic domain, MLL5 contains a PHD finger domain, a protein module that is often involved in binding to the N-terminus of histone H3. Here we report the NMR solution structure of the MLL5 PHD domain showing a variant of the canonical PHD fold that combines conserved H3 binding features from several classes of other PHD domains (including an aromatic cage) along with a novel C-terminal α-helix, not previously seen. We further demonstrate that the PHD domain binds with similar affinity to histone H3 tail peptides di- and tri-methylated at lysine 4 (H3K4me2 and H3K4me3), the former being the putative product of the MLL5 catalytic reaction. This work establishes the PHD domain of MLL5 as a bone fide ‘reader’ domain of H3K4 methyl marks suggesting that it may guide the spreading or further methylation of this site on chromatin.

Effects of NMR spectral resolution on protein structure calculation

Adequate digital resolution and signal sensitivity are two critical factors for protein structure determinations by solution NMR spectroscopy. The prime objective for obtaining high digital resolution is to resolve peak overlap, especially in NOESY spectra with thousands of signals where the signal analysis needs to be performed on a large scale. Achieving maximum digital resolution is usually limited by the practically available measurement time. We developed a method utilizing non-uniform sampling for balancing digital resolution and signal sensitivity, and performed a large-scale analysis of the effect of the digital resolution on the accuracy of the resulting protein structures. Structure calculations were performed as a function of digital resolution for about 400 proteins with molecular sizes ranging between 5 and 33 kDa. The structural accuracy was assessed by atomic coordinate RMSD values from the reference structures of the proteins. In addition, we monitored also the number of assigned NOESY cross peaks, the average signal sensitivity, and the chemical shift spectral overlap. We show that high resolution is equally important for proteins of every molecular size. The chemical shift spectral overlap depends strongly on the corresponding spectral digital resolution. Thus, knowing the extent of overlap can be a predictor of the resulting structural accuracy. Our results show that for every molecular size a minimal digital resolution, corresponding to the natural linewidth, needs to be achieved for obtaining the highest accuracy possible for the given protein size using state-of-the-art automated NOESY assignment and structure calculation methods.

Regulatory interactions between a bacterial tyrosine kinase and its cognate phosphatase

The cyclic process of autophosphorylation of the C-terminal tyrosine cluster (YC) of a bacterial tyrosine kinase and its subsequent dephosphorylation following interactions with a counteracting tyrosine phosphatase regulates diverse physiological processes, including the biosynthesis and export of polysaccharides responsible for the formation of biofilms or virulence-determining capsules. We provide here the first detailed insight into this hitherto uncharacterized regulatory interaction at residue-specific resolution using Escherichia coli Wzc, a canonical bacterial tyrosine kinase, and its opposing tyrosine phosphatase, Wzb. The phosphatase Wzb utilizes a surface distal to the catalytic elements of the kinase, Wzc, to dock onto its catalytic domain (WzcCD). WzcCD binds in a largely YC-independent fashion near the Wzb catalytic site, inducing allosteric changes therein. YC dephosphorylation is proximity-mediated and reliant on the elevated concentration of phosphorylated YC near the Wzb active site resulting from WzcCD docking. Wzb principally recognizes the phosphate of its phosphotyrosine substrate and further stabilizes the tyrosine moiety through ring stacking interactions with a conserved active site tyrosine.

Vivaldi: visualization and validation of biomacromolecular NMR structures from the PDB

We describe Vivaldi (VIsualization and VALidation DIsplay; http://pdbe.org/vivaldi), a web-based service for the analysis, visualization, and validation of NMR structures in the Protein Data Bank (PDB). Vivaldi provides access to model coordinates and several types of experimental NMR data using interactive visualization tools, augmented with structural annotations and model-validation information. The service presents information about the modeled NMR ensemble, validation of experimental chemical shifts, residual dipolar couplings, distance and dihedral angle constraints, as well as validation scores based on empirical knowledge and databases. Vivaldi was designed for both expert NMR spectroscopists and casual non-expert users who wish to obtain a better grasp of the information content and quality of NMR structures in the public archive. © Proteins 2013. © 2012 Wiley Periodicals, Inc.

Structure of a specialized acyl carrier protein essential for lipid A biosynthesis with very long-chain fatty acids in open and closed conformations

The solution nuclear magnetic resonance (NMR) structures and backbone (15)N dynamics of the specialized acyl carrier protein (ACP), RpAcpXL, from Rhodopseudomonas palustris, in both the apo form and holo form modified by covalent attachment of 4'-phosphopantetheine at S37, are virtually identical, monomeric, and correspond to the closed conformation. The structures have an extra α-helix compared to the archetypical ACP from Escherichia coli, which has four helices, resulting in a larger opening to the hydrophobic cavity. Chemical shift differences between apo- and holo-RpAcpXL indicated some differences in the hinge region between α2 and α3 and in the hydrophobic cavity environment, but corresponding changes in nuclear Overhauser effect cross-peak patterns were not detected. In contrast to the NMR structures, apo-RpAcpXL was observed in an open conformation in crystals that diffracted to 2.0 Å resolution, which resulted from movement of α3. On the basis of the crystal structure, the predicted biological assembly is a homodimer. Although the possible biological significance of dimerization is unknown, there is potential that the resulting large shared hydrophobic cavity could accommodate the very long-chain fatty acid (28-30 carbons) that this specialized ACP is known to synthesize and transfer to lipid A. These structures are the first representatives of the AcpXL family and the first to indicate that dimerization may be important for the function of these specialized ACPs.

Transient structure and dynamics in the disordered c-Myc transactivation domain affect Bin1 binding

The crucial role of Myc as an oncoprotein and as a key regulator of cell growth makes it essential to understand the molecular basis of Myc function. The N-terminal region of c-Myc coordinates a wealth of protein interactions involved in transformation, differentiation and apoptosis. We have characterized in detail the intrinsically disordered properties of Myc-1–88, where hierarchical phosphorylation of S62 and T58 regulates activation and destruction of the Myc protein. By nuclear magnetic resonance (NMR) chemical shift analysis, relaxation measurements and NOE analysis, we show that although Myc occupies a very heterogeneous conformational space, we find transiently structured regions in residues 22–33 and in the Myc homology box I (MBI; residues 45–65); both these regions are conserved in other members of the Myc family. Binding of Bin1 to Myc-1–88 as assayed by NMR and surface plasmon resonance (SPR) revealed primary binding to the S62 region in a dynamically disordered and multivalent complex, accompanied by population shifts leading to altered intramolecular conformational dynamics. These findings expand the increasingly recognized concept of intrinsically disordered regions mediating transient interactions to Myc, a key transcriptional regulator of major medical importance, and have important implications for further understanding its multifaceted role in gene regulation.

A software framework for analysing solid-state MAS NMR data

Solid-state magic-angle-spinning (MAS) NMR of proteins has undergone many rapid methodological developments in recent years, enabling detailed studies of protein structure, function and dynamics. Software development, however, has not kept pace with these advances and data analysis is mostly performed using tools developed for solution NMR which do not directly address solid-state specific issues. Here we present additions to the CcpNmr Analysis software package which enable easier identification of spinning side bands, straightforward analysis of double quantum spectra, automatic consideration of non-uniform labelling schemes, as well as extension of other existing features to the needs of solid-state MAS data. To underpin this, we have updated and extended the CCPN data model and experiment descriptions to include transfer types and nomenclature appropriate for solid-state NMR experiments, as well as a set of experiment prototypes covering the experiments commonly employed by solid-sate MAS protein NMR spectroscopists. This work not only improves solid-state MAS NMR data analysis but provides a platform for anyone who uses the CCPN data model for programming, data transfer, or data archival involving solid-state MAS NMR data.

PDBe: Protein Data Bank in Europe

The Protein Data Bank in Europe (PDBe; pdbe.org) is a partner in the Worldwide PDB organization (wwPDB; wwpdb.org) and as such actively involved in managing the single global archive of biomacromolecular structure data, the PDB. In addition, PDBe develops tools, services and resources to make structure-related data more accessible to the biomedical community. Here we describe recently developed, extended or improved services, including an animated structure-presentation widget (PDBportfolio), a widget to graphically display the coverage of any UniProt sequence in the PDB (UniPDB), chemistry- and taxonomy-based PDB-archive browsers (PDBeXplore), and a tool for interactive visualization of NMR structures, corresponding experimental data as well as validation and analysis results (Vivaldi).

Zn-binding AZUL domain of human ubiquitin protein ligase Ube3A

Ube3A (also referred to as E6AP for E6 Associated Protein) is a E3 ubiquitin-protein ligase implicated in the development of Angelman syndrome by controlling degradation of synaptic protein Arc and oncogenic papilloma virus infection by controlling degradation of p53. This article describe the solution NMR structure of the conserved N-terminal domain of human Ube3A (residues 24-87) that contains two residues (Cys44 and Arg62) found to be mutated in patients with Angelman syndrome. The structure of this domain adopts a novel Zn-binding fold we called AZUL (Amino-terminal Zn-finger of Ube3a Ligase). The AZUL domain has a helix-loop-helix architecture with a Zn ion coordinated by four Cys residues arranged in Cys-X(4)-Cys-X(4)-Cys-X(28)-Cys motif. Three of the Zn-bound residues are located in a 23-residue long and well structured loop that connects two α-helicies.

Structure of the BamC two-domain protein obtained by Rosetta with a limited NMR data set

The CS-RDC-NOE Rosetta program was used to generate the solution structure of a 27-kDa fragment of the Escherichia coli BamC protein from a limited set of NMR data. The BamC protein is a component of the essential five-protein β-barrel assembly machine in E. coli. The first 100 residues in BamC were disordered in solution. The Rosetta calculations showed that BamC₁₀₁₋₃₄₄ forms two well-defined domains connected by an ~18-residue linker, where the relative orientation of the domains was not defined. Both domains adopt a helix-grip fold previously observed in the Bet v 1 superfamily. ¹⁵N relaxation data indicated a high degree of conformational flexibility for the linker connecting the N-terminal domain and the C-terminal domain in BamC. The results here show that CS-RDC-NOE Rosetta is robust and has a high tolerance for misassigned nuclear Overhauser effect restraints, greatly simplifying NMR structure determinations.

Protein side-chain resonance assignment and NOE assignment using RDC-defined backbones without TOCSY data

One bottleneck in NMR structure determination lies in the laborious and time-consuming process of side-chain resonance and NOE assignments. Compared to the well-studied backbone resonance assignment problem, automated side-chain resonance and NOE assignments are relatively less explored. Most NOE assignment algorithms require nearly complete side-chain resonance assignments from a series of through-bond experiments such as HCCH-TOCSY or HCCCONH. Unfortunately, these TOCSY experiments perform poorly on large proteins. To overcome this deficiency, we present a novel algorithm, called NASCA: (NOE Assignment and Side-Chain Assignment), to automate both side-chain resonance and NOE assignments and to perform high-resolution protein structure determination in the absence of any explicit through-bond experiment to facilitate side-chain resonance assignment, such as HCCH-TOCSY. After casting the assignment problem into a Markov Random Field (MRF), NASCA: extends and applies combinatorial protein design algorithms to compute optimal assignments that best interpret the NMR data. The MRF captures the contact map information of the protein derived from NOESY spectra, exploits the backbone structural information determined by RDCs, and considers all possible side-chain rotamers. The complexity of the combinatorial search is reduced by using a dead-end elimination (DEE) algorithm, which prunes side-chain resonance assignments that are provably not part of the optimal solution. Then an A* search algorithm is employed to find a set of optimal side-chain resonance assignments that best fit the NMR data. These side-chain resonance assignments are then used to resolve the NOE assignment ambiguity and compute high-resolution protein structures. Tests on five proteins show that NASCA: assigns resonances for more than 90% of side-chain protons, and achieves about 80% correct assignments. The final structures computed using the NOE distance restraints assigned by NASCA: have backbone RMSD 0.8-1.5 Å from the reference structures determined by traditional NMR approaches.