Smooth statistical torsion angle potential derived from a large conformational database via adaptive kernel density estimation improves the quality of NMR protein structures

Xplor-NIH: Better parameters and protocols for NMR protein structure determination

The present work describes an update to the protein covalent geometry and atomic radii parameters in the Xplor-NIH biomolecular structure determination package. In combination with an improved treatment of selected non-bonded interactions between atoms three bonds apart, such as those involving methyl hydrogens, and a previously developed term that affects the system's gyration volume, the new parameters are tested using structure calculations on 30 proteins with restraints derived from nuclear magnetic resonance data. Using modern structure validation criteria, including several formally adopted by the Protein Data Bank, and a clear measure of structural accuracy, the results show superior performance relative to previous Xplor-NIH implementations. Additionally, the Xplor-NIH structures compare favorably against originally determined NMR models.

Structural basis of protein condensation on microtubules underlying branching microtubule nucleation

Targeting protein for Xklp2 (TPX2) is a key factor that stimulates branching microtubule nucleation during cell division. Upon binding to microtubules (MTs), TPX2 forms condensates via liquid-liquid phase separation, which facilitates recruitment of microtubule nucleation factors and tubulin. We report the structure of the TPX2 C-terminal minimal active domain (TPX2α5-α7) on the microtubule lattice determined by magic-angle-spinning NMR. We demonstrate that TPX2α5-α7 forms a co-condensate with soluble tubulin on microtubules and binds to MTs between two adjacent protofilaments and at the intersection of four tubulin heterodimers. These interactions stabilize the microtubules and promote the recruitment of tubulin. Our results reveal that TPX2α5-α7 is disordered in solution and adopts a folded structure on MTs, indicating that TPX2α5-α7 undergoes structural changes from unfolded to folded states upon binding to microtubules. The aromatic residues form dense interactions in the core, which stabilize folding of TPX2α5-α7 on microtubules. This work informs on how the phase-separated TPX2α5-α7 behaves on microtubules and represents an atomic-level structural characterization of a protein that is involved in a condensate on cytoskeletal filaments.

Structural basis of HIV-1 maturation inhibitor binding and activity

HIV-1 maturation inhibitors (MIs), Bevirimat (BVM) and its analogs interfere with the catalytic cleavage of spacer peptide 1 (SP1) from the capsid protein C-terminal domain (CACTD), by binding to and stabilizing the CACTD-SP1 region. MIs are under development as alternative drugs to augment current antiretroviral therapies. Although promising, their mechanism of action and associated virus resistance pathways remain poorly understood at the molecular, biochemical, and structural levels. We report atomic-resolution magic-angle-spinning NMR structures of microcrystalline assemblies of CACTD-SP1 complexed with BVM and/or the assembly cofactor inositol hexakisphosphate (IP6). Our results reveal a mechanism by which BVM disrupts maturation, tightening the 6-helix bundle pore and quenching the motions of SP1 and the simultaneously bound IP6. In addition, BVM-resistant SP1-A1V and SP1-V7A variants exhibit distinct conformational and binding characteristics. Taken together, our study provides a structural explanation for BVM resistance as well as guidance for the design of new MIs.

Divergent self-association properties of paralogous proteins TRIM2 and TRIM3 regulate their E3 ligase activity

Tripartite motif (TRIM) proteins constitute a large family of RING-type E3 ligases that share a conserved domain architecture. TRIM2 and TRIM3 are paralogous class VII TRIM members that are expressed mainly in the brain and regulate different neuronal functions. Here we present a detailed structure-function analysis of TRIM2 and TRIM3, which despite high sequence identity, exhibit markedly different self-association and activity profiles. We show that the isolated RING domain of human TRIM3 is monomeric and inactive, and that this lack of activity is due to a few placental mammal-specific amino acid changes adjacent to the core RING domain that prevent self-association but not E2 recognition. We demonstrate that the activity of human TRIM3 RING can be restored by substitution with the relevant region of human TRIM2 or by hetero-dimerization with human TRIM2, establishing that subtle amino acid changes can profoundly affect TRIM protein activity. Finally, we show that TRIM2 and TRIM3 interact in a cellular context via their filamin and coiled-coil domains, respectively.

Magic-angle-spinning NMR structure of the kinesin-1 motor domain assembled with microtubules reveals the elusive neck linker orientation

Microtubules (MTs) and their associated proteins play essential roles in maintaining cell structure, organelle transport, cell motility, and cell division. Two motors, kinesin and cytoplasmic dynein link the MT network to transported cargos using ATP for force generation. Here, we report an all-atom NMR structure of nucleotide-free kinesin-1 motor domain (apo-KIF5B) in complex with paclitaxel-stabilized microtubules using magic-angle-spinning (MAS) NMR spectroscopy. The structure reveals the position and orientation of the functionally important neck linker and how ADP induces structural and dynamic changes that ensue in the neck linker. These results demonstrate that the neck linker is in the undocked conformation and oriented in the direction opposite to the KIF5B movement. Chemical shift perturbations and intensity changes indicate that a significant portion of ADP-KIF5B is in the neck linker docked state. This study also highlights the unique capability of MAS NMR to provide atomic-level information on dynamic regions of biological assemblies.

A kink in DWORF helical structure controls the activation of the sarcoplasmic reticulum Ca2+-ATPase

SERCA is a P-type ATPase embedded in the sarcoplasmic reticulum and plays a central role in muscle relaxation. SERCA's function is regulated by single-pass membrane proteins called regulins. Unlike other regulins, dwarf open reading frame (DWORF) expressed in cardiac muscle has a unique activating effect. Here, we determine the structure and topology of DWORF in lipid bilayers using a combination of oriented sample solid-state NMR spectroscopy and replica-averaged orientationally restrained molecular dynamics. We found that DWORF's structural topology consists of a dynamic N-terminal domain, an amphipathic juxtamembrane helix that crosses the lipid groups at an angle of 64°, and a transmembrane C-terminal helix with an angle of 32°. A kink induced by Pro15, unique to DWORF, separates the two helical domains. A single Pro15Ala mutant significantly decreases the kink and eliminates DWORF's activating effect on SERCA. Overall, our findings directly link DWORF's structural topology to its activating effect on SERCA.

The structure of a minimum amyloid fibril core formed by necroptosis-mediating RHIM of human RIPK3

Receptor-interacting protein kinases 3 (RIPK3), a central node in necroptosis, polymerizes in response to the upstream signals and then activates its downstream mediator to induce cell death. The active polymeric form of RIPK3 has been indicated as the form of amyloid fibrils assembled via its RIP homotypic interaction motif (RHIM). In this study, we combine cryogenic electron microscopy and solid-state NMR to determine the amyloid fibril structure of RIPK3 RHIM-containing C-terminal domain (CTD). The structure reveals a single protofilament composed of the RHIM domain. RHIM forms three β-strands (referred to as strands 1 through 3) folding into an S shape, a distinct fold from that in complex with RIPK1. The consensus tetrapeptide VQVG of RHIM forms strand 2, which zips up strands 1 and 3 via heterozipper-like interfaces. Notably, the RIPK3-CTD fibril, as a physiological fibril, exhibits distinctive assembly compared with pathological fibrils. It has an exceptionally small fibril core and twists in both handedness with the smallest pitch known so far. These traits may contribute to a favorable spatial arrangement of RIPK3 kinase domain for efficient phosphorylation.

The amyloid structure of mouse RIPK3 (receptor interacting protein kinase 3) in cell necroptosis

RIPK3 amyloid complex plays crucial roles during TNF-induced necroptosis and in response to immune defense in both human and mouse. Here, we have structurally characterized mouse RIPK3 homogeneous self-assembly using solid-state NMR, revealing a well-ordered N-shaped amyloid core structure featured with 3 parallel in-register β-sheets. This structure differs from previously published human RIPK1/RIPK3 hetero-amyloid complex structure, which adopted a serpentine fold. Functional studies indicate both RIPK1-RIPK3 binding and RIPK3 amyloid formation are essential but not sufficient for TNF-induced necroptosis. The structural integrity of RIPK3 fibril with three β-strands is necessary for signaling. Molecular dynamics simulations with a mouse RIPK1/RIPK3 model indicate that the hetero-amyloid is less stable when adopting the RIPK3 fibril conformation, suggesting a structural transformation of RIPK3 from RIPK1-RIPK3 binding to RIPK3 amyloid formation. This structural transformation would provide the missing link connecting RIPK1-RIPK3 binding to RIPK3 homo-oligomer formation in the signal transduction.

Receptor Interacting Protein Kinase 3 (RIPK3) has a key role in TNF-induced necroptosis. Here, the authors combine solid state NMR measurements, MD simulations and cell based assays to characterize mouse RIPK3 and they present the structure of the RIPK3 amyloid core.

Structural Flexibility of Cyclosporine A Is Mediated by Amide Cis-Trans Isomerization and the Chameleonic Roles of Calcium

Falling outside of Lipinski's rule of five, macrocyclic drugs have accessed unique binding sites of their target receptors unreachable by traditional small molecules. Cyclosporin(e) A (CycA), an extensively studied macrocyclic natural product, is an immunosuppressant with undesirable side effects such as electrolytic imbalances. In this work, a comprehensive view on the conformational landscape of CycA, its interactions with Ca²⁺, and host-guest interactions with cyclophilin A (CypA) is reported through exhaustive analyses that combine ion-mobility spectrometry-mass spectrometry (IMS-MS), nuclear magnetic resonance (NMR) spectroscopy, distance-geometry modeling, and NMR-driven molecular dynamics. Our IMS-MS data show that CycA can adopt extremely compact conformations with significantly smaller collisional cross sections than the closed conformation observed in CDCl₃. To adopt these conformations, the macrocyclic ring has to twist and bend via cis-trans isomerization of backbone amides, and thus, we termed this family of structures the "bent" conformation. Furthermore, NMR measurements indicate that the closed conformation exists at 19% in CD₃OD/H₂O and 55% in CD₃CN. However, upon interacting with Ca²⁺, in addition to the bent and previously reported closed conformations of free CycA, the CycA:Ca²⁺ complex is open and has all-trans peptide bonds. Previous NMR studies using calcium perchlorate reported only the closed conformation of CycA (which contains one cis peptide bond). Here, calcium chloride, a more biologically relevant salt, was used, and interestingly, it helps converting the cis-MeLeu9-MeLeu10 peptide bond into a trans bond. Last, we were able to capture the native binding of CycA and CypA to give forth evidence that IMS-MS is able to probe the solution-phase structures of the complexes and that the Ca²⁺:CycA complex may play an essential role in the binding of CycA to CypA.

Atomic-resolution structure of HIV-1 capsid tubes by magic-angle spinning NMR

HIV-1 capsid plays multiple key roles in viral replication, and inhibition of capsid assembly is an attractive target for therapeutic intervention. Here, we report the atomic-resolution structure of the capsid protein (CA) tubes, determined by magic-angle-spinning NMR and data-guided molecular dynamics simulations. Functionally important regions, including the NTD β-hairpin, the cyclophilin A loop, residues in the hexamer center pore, and the NTD-CTD linker region, are well defined. The structure of individual CA chains, their arrangement in the pseudo-hexameric units of the tube and the inter-hexamer interfaces are consistent with those in intact capsid cores and substantially different from the organization in crystal structures, which featured flat hexamers. The inherent curvature in the CA tubes is controlled by conformational variability of residues in the linker region and of dimer and trimer interfaces. The present structure reveals atomic-level detail into capsid architecture and provides important guidance for the design of novel capsid inhibitors.

Tightening the Crosslinking Distance Restraints for Better Resolution of Protein Structure and Dynamics

Chemical crosslinking coupled with mass spectrometry (CXMS) has been increasingly used in structural biology. CXMS distance restraints are usually applied to Cα or Cβ atoms of the crosslinked residues, with upper bounds typically over 20 Å. The incorporation of loose CXMS restraints only marginally improves the resolution of the calculated structures. Here, we present a revised format of CXMS distance restraints, which works by first modifying the crosslinked residue with a rigid extension derived from the crosslinker. With the flexible side chain explicitly represented, the reformatted restraint can be applied to the modification group instead, with an upper bound of 6 Å or less. The short distance restraint can be represented and back-calculated simply with a straight line. The use of tighter restraints not only afford better-resolved structures but also uncover protein dynamics. Together, our approach enables more information extracted from the CXMS data.

Inter-domain dynamics in the chaperone SurA and multi-site binding to its outer membrane protein clients

The periplasmic chaperone SurA plays a key role in outer membrane protein (OMP) biogenesis. E. coli SurA comprises a core domain and two peptidylprolyl isomerase domains (P1 and P2), but its mechanisms of client binding and chaperone function have remained unclear. Here, we use chemical cross-linking, hydrogen-deuterium exchange mass spectrometry, single-molecule FRET and molecular dynamics simulations to map the client binding site(s) on SurA and interrogate the role of conformational dynamics in OMP recognition. We demonstrate that SurA samples an array of conformations in solution in which P2 primarily lies closer to the core/P1 domains than suggested in the SurA crystal structure. OMP binding sites are located primarily in the core domain, and OMP binding results in conformational changes between the core/P1 domains. Together, the results suggest that unfolded OMP substrates bind in a cradle formed between the SurA domains, with structural flexibility between domains assisting OMP recognition, binding and release.

Pro-islet amyloid polypeptide in micelles contains a helical prohormone segment

Pro-islet amyloid polypeptide (proIAPP) is the prohormone precursor molecule to IAPP, also known as amylin. IAPP is a calcitonin family peptide hormone that is cosecreted with insulin, and largely responsible for hunger satiation and metabolic homeostasis. Amyloid plaques containing mixtures of mature IAPP and misprocessed proIAPP deposit on, and destroy pancreatic β-cell membranes, and they are recognized as a clinical hallmark of type 2 diabetes mellitus. In order to better understand the interaction with cellular membranes, we solved the solution NMR structure of proIAPP bound to dodecylphosphocholine micelles at pH 4.5. We show that proIAPP is a dynamic molecule with four α-helices. The first two helices are contained within the mature IAPP sequence, while the second two helices are part of the C-terminal prohormone segment (Cpro). We mapped the membrane topology of the amphipathic helices by paramagnetic relaxation enhancement, and we used CD and diffusion-ordered spectroscopy to identify environmental factors that impact proIAPP membrane affinity. We discuss how our structural results relate to prohormone processing based on the varied pH environments and lipid compositions of organelle membranes within the regulated secretory pathway, and the likelihood of Cpro survival for cosecretion with IAPP. DATABASE: The assigned resonances have been deposited in the Biological Magnetic Resonance Bank (BMRB) with accession numbers 50007 and 50019 for proIAPP and Cpro, respectively. The lowest energy structures have been deposited in the Protein Data Bank (PDB) with access codes 6UCJ and 6UCK.

On the necessity of an integrative approach to understand protein structural dynamics

Proteins are dynamic, fluctuating between multiple conformational states. Protein dynamics, spanning orders of magnitude in time and space, allow proteins to perform specific functions. Moreover, under certain conditions, proteins can morph into a different set of conformations. Thus, a complete understanding of protein structural dynamics can provide mechanistic insights into protein function. Here, we review the latest developments in methods used to determine protein ensemble structures and to characterize protein dynamics. Techniques including X-ray crystallography, cryogenic electron microscopy, and small angle scattering can provide structural information on specific conformational states or on the averaged shape of the protein, whereas techniques including nuclear magnetic resonance, fluorescence resonance energy transfer (FRET), and chemical cross-linking coupled with mass spectrometry provide information on the fluctuation of the distances between protein domains, residues, and atoms for the multiple conformational states of the protein. In particular, FRET measurements at the single-molecule level allow rapid resolution of protein conformational states, where information is otherwise obscured in bulk measurements. Taken together, the different techniques complement each other and their integrated use can offer a clear picture of protein structure and dynamics.

A three-dimensional potential of mean force to improve backbone and sidechain hydrogen bond geometry in Xplor-NIH protein structure determination

We introduce a new hydrogen bonding potential of mean force generated from high-quality crystal structures for use in Xplor-NIH structure calculations. This term applies to hydrogen bonds involving both backbone and sidechain atoms. When used in structure refinement calculations of 10 example protein systems with experimental distance, dihedral and residual dipolar coupling restraints, we demonstrate that the new term has superior performance to the previously developed hydrogen bonding potential of mean force used in Xplor-NIH.

Accuracy and precision of protein structures determined by magic angle spinning NMR spectroscopy: for some 'with a little help from a friend'

We present a systematic investigation into the attainable accuracy and precision of protein structures determined by heteronuclear magic angle spinning solid-state NMR for a set of four proteins of varied size and secondary structure content. Structures were calculated using synthetically generated random sets of C-C distances up to 7 Å at different degrees of completeness. For single-domain proteins, 9-15 restraints per residue are sufficient to derive an accurate model structure, while maximum accuracy and precision are reached with over 15 restraints per residue. For multi-domain proteins and protein assemblies, additional information on domain orientations, quaternary structure and/or protein shape is needed. As demonstrated for the HIV-1 capsid protein assembly, this can be accomplished by integrating MAS NMR with cryoEM data. In all cases, inclusion of TALOS-derived backbone torsion angles improves the accuracy for small number of restraints, while no further increases are noted for restraint completeness above 40%. In contrast, inclusion of TALOS-derived torsion angle restraints consistently increases the precision of the structural ensemble at all degrees of distance restraint completeness.

Solution NMR Spectroscopy for the Determination of Structures of Membrane Proteins in a Lipid Environment

NMR spectroscopy has harnessed the recent technical advances to emerge as a competitive, elegant, and eminently viable technique for determining the solution structures of membrane proteins at the level of atomic resolution. Once a good level of cell-based or cell-free expression and purification of a suitably sized membrane protein has been achieved, then NMR offers a combination of several versatile strategies, for example choice of appropriate deuterated or nondeuterated detergents, temperature, and ionic strength; isotope labeling with ²H, ¹³C, ¹⁵N, with or without protonation of Ile (δ1), Leu, and Val methyl protons; combinatorial labeling or unlabeling of specific amino acids; TROSY based-, nonuniform sampling (NUS) based-, and other NMR experiments; measurement of residual dipolar couplings using stretched polyacrylamide gels or DNA nanotubes; spin labeling and paramagnetic relaxation enhancements (PRE). Strategic combinations of these advancements together with availability of highly sensitive cryogenically cooled-probes equipped high-field NMR spectrometers (up to 1 GHz ¹H frequency) have allowed the perseverant investigator to successfully overcome several of the conventional pitfalls associated with the NMR technique and membrane proteins, viz., low sensitivity, poor sample stability, spectral crowding, and a limited number of NOEs and other constraints for structure calculations. This has resulted in an unprecedented growth in the number of successfully determined NMR structures of large and complex membrane proteins over the last two decades, and this technique now holds great promise for the structure determination of an ever larger body of membrane proteins.

Characterisation of class VI TRIM RING domains: linking RING activity to C-terminal domain identity

TRIM E3 ubiquitin ligases regulate multiple cellular processes, and their dysfunction is linked to disease. They are characterised by a conserved N-terminal tripartite motif comprising a RING, B-box domains, and a coiled-coil region, with C-terminal domains often mediating substrate recruitment. TRIM proteins are grouped into 11 classes based on C-terminal domain identity. Class VI TRIMs, TRIM24, TRIM33, and TRIM28, have been described as transcriptional regulators, a function linked to their C-terminal plant homeodomain and bromodomain, and independent of their ubiquitination activity. It is unclear whether E3 ligase activity is regulated in family members where the C-terminal domains function independently. Here, we provide a detailed biochemical characterisation of the RING domains of class VI TRIMs and describe the solution structure of the TRIM28 RING. Our study reveals a lack of activity of the isolated RING domains, which may be linked to the absence of self-association. We propose that class VI TRIMs exist in an inactive state and require additional regulatory events to stimulate E3 ligase activity, ensuring that associated chromatin-remodelling factors are not injudiciously degraded.

Structural basis for the recognition of K48-linked Ub chain by proteasomal receptor Rpn13

The interaction between K48-linked ubiquitin (Ub) chain and Rpn13 is important for proteasomal degradation of ubiquitinated substrate proteins. Only the complex structure between the N-terminal domain of Rpn13 (Rpn13^NTD) and Ub monomer has been characterized, while it remains unclear how Rpn13 specifically recognizes K48-linked Ub chain. Using single-molecule FRET, here we show that K48-linked diubiquitin (K48-diUb) fluctuates among distinct conformational states, and a preexisting compact state is selectively enriched by Rpn13^NTD. The same binding mode is observed for full-length Rpn13 and longer K48-linked Ub chain. Using solution NMR spectroscopy, we have determined the complex structure between Rpn13^NTD and K48-diUb. In this structure, Rpn13^NTD simultaneously interacts with proximal and distal Ub subunits of K48-diUb that remain associated in the complex, thus corroborating smFRET findings. The proximal Ub interacts with Rpn13^NTD similarly as the Ub monomer in the known Rpn13^NTD:Ub structure, while the distal Ub binds to a largely electrostatic surface of Rpn13^NTD. Thus, a charge-reversal mutation in Rpn13^NTD weakens the interaction between Rpn13 and K48-linked Ub chain, causing accumulation of ubiquitinated proteins. Moreover, physical blockage of the access of the distal Ub to Rpn13^NTD with a proximity-attached Ub monomer can disrupt the interaction between Rpn13 and K48-diUb. Taken together, the bivalent interaction of K48-linked Ub chain with Rpn13 provides the structural basis for Rpn13 linkage selectivity, which opens a new window for modulating proteasomal function.

The C-terminal domain of HpDprA is a DNA-binding winged helix domain that does not bind double-stranded DNA

DNA-processing protein A, a ubiquitous multidomain DNA-binding protein, plays a crucial role during natural transformation in bacteria. Here, we carried out the structural analysis of DprA from the human pathogen Helicobacter pylori by combining data issued from the 1.8-Å resolution X-ray structure of the Pfam02481 domain dimer (RF), the NMR structure of the carboxy terminal domain (CTD), and the low-resolution structure of the full-length DprA dimer obtained in solution by SAXS. In particular, we sought a molecular function for the CTD, a domain that we show here is essential for transformation in H. pylori. Albeit its structural homology to winged helix DNA-binding motifs, we confirmed that the isolated CTD does not interact with ssDNA nor with dsDNA. The key R52 and K137 residues of RF are crucial for these two interactions. Search for sequences harboring homology to either HpDprA or Rhodopseudomonas palustris DprA CTDs led to the identification of conserved patches in the two CTD. Our structural study revealed the similarity of the structures adopted by these residues in RpDprA CTD and HpDprA CTD. This argues for a conserved, but yet to be defined, CTD function, distinct from DNA binding.

Theory and practice of using solvent paramagnetic relaxation enhancement to characterize protein conformational dynamics

Paramagnetic relaxation enhancement (PRE) has been established as a powerful tool in NMR for investigating protein structure and dynamics. The PRE is usually measured with a paramagnetic probe covalently attached at a specific site of an otherwise diamagnetic protein. The present work provides the numerical formulation for probing protein structure and conformational dynamics based on the solvent PRE (sPRE) measurement, using two alternative approaches. An inert paramagnetic cosolute randomly collides with the protein, and the resulting sPRE manifests the relative solvent exposure of protein nuclei. To make the back-calculated sPRE values most consistent with the observed values, the protein structure is either refined against the sPRE, or an ensemble of conformers is selected from a pre-generated library using a Monte Carlo algorithm. The ensemble structure comprises either N conformers of equal occupancy, or two conformers with different relative populations. We demonstrate the sPRE method using GB1, a structurally rigid protein, and calmodulin, a protein comprising two domains and existing in open and closed states. The sPRE can be computed with a stand-alone program for rapid evaluation, or with the invocation of a module in the latest release of the structure calculation software Xplor-NIH. As a label-free method, the sPRE measurement can be readily integrated with other biophysical techniques. The current limitations of the sPRE method are also discussed, regarding accurate measurement and theoretical calculation, model selection and suitable timescale.

Protein Structure Elucidation from NMR Data with the Program Xplor-NIH

Xplor-NIH is a popular software package for biomolecular structure determination from NMR (and other) experimental data. This chapter illustrates its use with the de novo structure determination of the B1 domain of streptococcal protein G (GB1), based on distances from nuclear Overhauser effects, torsion angles from scalar couplings, and bond-vector orientations from residual dipolar couplings. Including Xplor-NIH's latest developments, a complete structure calculation script is discussed in detail, and is intended to serve as a basis for other applications.

Xplor-NIH for molecular structure determination from NMR and other data sources

Xplor-NIH is a popular software package for biomolecular structure determination from nuclear magnetic resonance (NMR) and other data sources. Here, some of Xplor-NIH's most useful data-associated energy terms are reviewed, including newer alternative options for using residual dipolar coupling data in structure calculations. Further, we discuss new developments in the implementation of strict symmetry for the calculation of symmetric homo-oligomers, and in the representation of the system as an ensemble of structures to account for motional effects. Finally, the different available force fields are presented, among other Xplor-NIH capabilities.

Structural Insights into the Yersinia pestis Outer Membrane Protein Ail in Lipid Bilayers

Yersinia pestis the causative agent of plague, is highly pathogenic and poses very high risk to public health. The outer membrane protein Ail (Adhesion invasion locus) is one of the most highly expressed proteins on the cell surface of Y. pestis, and a major target for the development of medical countermeasures. Ail is essential for microbial virulence and is critical for promoting the survival of Y. pestis in serum. Structures of Ail have been determined by X-ray diffraction and solution NMR spectroscopy, but the protein's activity is influenced by the detergents in these samples, underscoring the importance of the surrounding environment for structure-activity studies. Here we describe the backbone structure of Ail, determined in lipid bilayer nanodiscs, using solution NMR spectroscopy. We also present solid-state NMR data obtained for Ail in membranes containing lipopolysaccharide (LPS), a major component of the bacterial outer membranes. The protein in lipid bilayers, adopts the same eight-stranded β-barrel fold observed in the crystalline and micellar states. The membrane composition, however, appears to have a marked effect on protein dynamics, with LPS enhancing conformational order and slowing down the ¹⁵N transverse relaxation rate. The results provide information about the way in which an outer membrane protein inserts and functions in the bacterial membrane.

Ubiquitin S65 phosphorylation engenders a pH-sensitive conformational switch

Ubiquitin (Ub) is an important signaling protein. Recent studies have shown that Ub can be enzymatically phosphorylated at S65, and that the resulting pUb exhibits two conformational states-a relaxed state and a retracted state. However, crystallization efforts have yielded only the structure for the relaxed state, which was found similar to that of unmodified Ub. Here we present the solution structures of pUb in both states obtained through refinement against state-specific NMR restraints. We show that the retracted state differs from the relaxed state by the retraction of the last β-strand and by the extension of the second α-helix. Further, we show that at 7.2, the pK_a value for the phosphoryl group in the relaxed state is higher by 1.4 units than that in the retracted state. Consequently, pUb exists in equilibrium between protonated and deprotonated forms and between retracted and relaxed states, with protonated/relaxed species enriched at slightly acidic pH and deprotonated/retracted species enriched at slightly basic pH. The heterogeneity of pUb explains the inability of phosphomimetic mutants to fully mimic pUb. The pH-sensitive conformational switch is likely preserved for polyubiquitin, as single-molecule FRET data indicate that pH change leads to quaternary rearrangement of a phosphorylated K63-linked diubiquitin. Because cellular pH varies among compartments and changes upon pathophysiological insults, our finding suggests that pH and Ub phosphorylation confer additional target specificities and enable an additional layer of modulation for Ub signals.

Lanthanoid tagging via an unnatural amino acid for protein structure characterization

Lanthanoid pseudo-contact shift (PCS) provides long-range structural information between a paramagnetic tag and protein nuclei. However, for proteins with native cysteines, site-specific attachment may only utilize functional groups orthogonal to sulfhydryl chemistry. Here we report two lanthanoid probes, DTTA-C3-yne and DTTA-C4-yne, which can be conjugated to an unnatural amino acid pAzF in the target protein via azide-alkyne cycloaddition. Demonstrated with ubiquitin and cysteine-containing enzyme EIIB, we show that large PCSs of distinct profiles can be generated for each tag/lanthanoid combination. The DTTA-based lanthanoid tags are associated with large magnetic susceptibility tensors owing to the rigidity of the tags. In particular, introduction of the DTTA-C3 tag affords intermolecular PCSs and enables structural characterization of a transient protein complex between ubiquitin and a UBA domain. Together, we have expanded the repertoire of paramagnetic tags and the applicability of paramagnetic NMR.

Modeling Protein Excited-state Structures from "Over-length" Chemical Cross-links

Chemical cross-linking coupled with mass spectroscopy (CXMS) provides proximity information for the cross-linked residues and is used increasingly for modeling protein structures. However, experimentally identified cross-links are sometimes incompatible with the known structure of a protein, as the distance calculated between the cross-linked residues far exceeds the maximum length of the cross-linker. The discrepancies may persist even after eliminating potentially false cross-links and excluding intermolecular ones. Thus the "over-length" cross-links may arise from alternative excited-state conformation of the protein. Here we present a method and associated software DynaXL for visualizing the ensemble structures of multidomain proteins based on intramolecular cross-links identified by mass spectrometry with high confidence. Representing the cross-linkers and cross-linking reactions explicitly, we show that the protein excited-state structure can be modeled with as few as two over-length cross-links. We demonstrate the generality of our method with three systems: calmodulin, enzyme I, and glutamine-binding protein, and we show that these proteins alternate between different conformations for interacting with other proteins and ligands. Taken together, the over-length chemical cross-links contain valuable information about protein dynamics, and our findings here illustrate the relationship between dynamic domain movement and protein function.

High quality NMR structures: a new force field with implicit water and membrane solvation for Xplor-NIH

Structure determination of proteins by NMR is unique in its ability to measure restraints, very accurately, in environments and under conditions that closely mimic those encountered in vivo. For example, advances in solid-state NMR methods enable structure determination of membrane proteins in detergent-free lipid bilayers, and of large soluble proteins prepared by sedimentation, while parallel advances in solution NMR methods and optimization of detergent-free lipid nanodiscs are rapidly pushing the envelope of the size limit for both soluble and membrane proteins. These experimental advantages, however, are partially squandered during structure calculation, because the commonly used force fields are purely repulsive and neglect solvation, Van der Waals forces and electrostatic energy. Here we describe a new force field, and updated energy functions, for protein structure calculations with EEFx implicit solvation, electrostatics, and Van der Waals Lennard-Jones forces, in the widely used program Xplor-NIH. The new force field is based primarily on CHARMM22, facilitating calculations with a wider range of biomolecules. The new EEFx energy function has been rewritten to enable OpenMP parallelism, and optimized to enhance computation efficiency. It implements solvation, electrostatics, and Van der Waals energy terms together, thus ensuring more consistent and efficient computation of the complete nonbonded energy lists. Updates in the related python module allow detailed analysis of the interaction energies and associated parameters. The new force field and energy function work with both soluble proteins and membrane proteins, including those with cofactors or engineered tags, and are very effective in situations where there are sparse experimental restraints. Results obtained for NMR-restrained calculations with a set of five soluble proteins and five membrane proteins show that structures calculated with EEFx have significant improvements in accuracy, precision, and conformation, and that structure refinement can be obtained by short relaxation with EEFx to obtain improvements in these key metrics. These developments broaden the range of biomolecular structures that can be calculated with high fidelity from NMR restraints.

High-Resolution NMR Determination of the Dynamic Structure of Membrane Proteins

¹⁵N spin‐relaxation rates are demonstrated to provide critical information about the long‐range structure and internal motions of membrane proteins. Combined with an improved calculation method, the relaxation‐rate‐derived structure of the 283‐residue human voltage‐dependent anion channel revealed an anisotropically shaped barrel with a rigidly attached N‐terminal helix. Our study thus establishes an NMR spectroscopic approach to determine the structure and dynamics of mammalian membrane proteins at high accuracy and resolution.

Solid-state NMR structure of a pathogenic fibril of full-length human α-synuclein

Misfolded α-synuclein amyloid fibrils are the principal components of Lewy bodies and neurites, hallmarks of Parkinson's disease (PD). We present a high-resolution structure of an α-synuclein fibril, in a form that induces robust pathology in primary neuronal culture, determined by solid-state NMR spectroscopy and validated by EM and X-ray fiber diffraction. Over 200 unique long-range distance restraints define a consensus structure with common amyloid features including parallel, in-register β-sheets and hydrophobic-core residues, and with substantial complexity arising from diverse structural features including an intermolecular salt bridge, a glutamine ladder, close backbone interactions involving small residues, and several steric zippers stabilizing a new orthogonal Greek-key topology. These characteristics contribute to the robust propagation of this fibril form, as supported by the structural similarity of early-onset-PD mutants. The structure provides a framework for understanding the interactions of α-synuclein with other proteins and small molecules, to aid in PD diagnosis and treatment.

Improving NMR Structures of RNA

Here, we show that modern solution nuclear magnetic resonance (NMR) structures of RNA exhibit more steric clashes and conformational ambiguities than their crystallographic X-ray counterparts. To tackle these issues, we developed RNA-ff1, a new force field for structure calculation with Xplor-NIH. Using seven published NMR datasets, RNA-ff1 improves covalent geometry and MolProbity validation criteria for clashes and backbone conformation in most cases, relative to both the previous Xplor-NIH force field and the original structures associated with the experimental data. In addition, with smaller base-pair step rises in helical stems, RNA-ff1 structures enjoy more favorable base stacking. Finally, structural accuracy improves in the majority of cases, as supported by complete residual dipolar coupling cross-validation. Thus, the reported advances show great promise in bridging the quality gap that separates NMR and X-ray structures of RNA.

Structural insights into the mechanism of activation of the TRPV1 channel by a membrane-bound tarantula toxin

Venom toxins are invaluable tools for exploring the structure and mechanisms of ion channels. Here, we solve the structure of double-knot toxin (DkTx), a tarantula toxin that activates the heat-activated TRPV1 channel. We also provide improved structures of TRPV1 with and without the toxin bound, and investigate the interactions of DkTx with the channel and membranes. We find that DkTx binds to the outer edge of the external pore of TRPV1 in a counterclockwise configuration, using a limited protein-protein interface and inserting hydrophobic residues into the bilayer. We also show that DkTx partitions naturally into membranes, with the two lobes exhibiting opposing energetics for membrane partitioning and channel activation. Finally, we find that the toxin disrupts a cluster of hydrophobic residues behind the selectivity filter that are critical for channel activation. Collectively, our findings reveal a novel mode of toxin-channel recognition that has important implications for the mechanism of thermosensation.

Dynamic equilibrium between closed and partially closed states of the bacterial Enzyme I unveiled by solution NMR and X-ray scattering

Enzyme I (EI) is the first component in the bacterial phosphotransferase system, a signal transduction pathway in which phosphoryl transfer through a series of bimolecular protein-protein interactions is coupled to sugar transport across the membrane. EI is a multidomain, 128-kDa homodimer that has been shown to exist in two conformational states related to one another by two large (50-90°) rigid body domain reorientations. The open conformation of apo EI allows phosphoryl transfer from His189 located in the N-terminal domain α/β (EIN(α/β)) subdomain to the downstream protein partner bound to the EIN(α) subdomain. The closed conformation, observed in a trapped phosphoryl transfer intermediate, brings the EIN(α/β) subdomain into close proximity to the C-terminal dimerization domain (EIC), thereby permitting in-line phosphoryl transfer from phosphoenolpyruvate (PEP) bound to EIC to His189. Here, we investigate the solution conformation of a complex of an active site mutant of EI (H189A) with PEP. Simulated annealing refinement driven simultaneously by solution small angle X-ray scattering and NMR residual dipolar coupling data demonstrates unambiguously that the EI(H189A)-PEP complex exists in a dynamic equilibrium between two approximately equally populated conformational states, one corresponding to the closed structure and the other to a partially closed species. The latter likely represents an intermediate in the open-to-closed transition.

Backbone structure of Yersinia pestis Ail determined in micelles by NMR-restrained simulated annealing with implicit membrane solvation

The outer membrane protein Ail (attachment invasion locus) is a virulence factor of Yersinia pestis that mediates cell invasion, cell attachment and complement resistance. Here we describe its three-dimensional backbone structure determined in decyl-phosphocholine (DePC) micelles by NMR spectroscopy. The NMR structure was calculated using the membrane function of the implicit solvation potential, eefxPot, which we have developed to facilitate NMR structure calculations in a physically realistic environment. We show that the eefxPot force field guides the protein towards its native fold. The resulting structures provide information about the membrane-embedded global position of Ail, and have higher accuracy, higher precision and improved conformational properties, compared to the structures calculated with the standard repulsive potential.

A Practical Implicit Membrane Potential for NMR Structure Calculations of Membrane Proteins

The highly anisotropic environment of the lipid bilayer membrane imposes significant constraints on the structures and functions of membrane proteins. However, NMR structure calculations typically use a simple repulsive potential that neglects the effects of solvation and electrostatics, because explicit atomic representation of the solvent and lipid molecules is computationally expensive and impractical for routine NMR-restrained calculations that start from completely extended polypeptide templates. Here, we describe the extension of a previously described implicit solvation potential, eefxPot, to include a membrane model for NMR-restrained calculations of membrane protein structures in XPLOR-NIH. The key components of eefxPot are an energy term for solvation free energy that works together with other nonbonded energy functions, a dedicated force field for conformational and nonbonded protein interaction parameters, and a membrane function that modulates the solvation free energy and dielectric screening as a function of the atomic distance from the membrane center, relative to the membrane thickness. Initial results obtained for membrane proteins with structures determined experimentally in lipid bilayer membranes show that eefxPot affords significant improvements in structural quality, accuracy, and precision. Calculations with eefxPot are straightforward to implement and can be used to both fold and refine structures, as well as to run unrestrained molecular-dynamics simulations. The potential is entirely compatible with the full range of experimental restraints measured by various techniques. Overall, it provides a useful and practical way to calculate membrane protein structures in a physically realistic environment.

Application of information theory to a three-body coarse-grained representation of proteins in the PDB: insights into the structural and evolutionary roles of residues in protein structure

Knowledge-based methods for analyzing protein structures, such as statistical potentials, primarily consider the distances between pairs of bodies (atoms or groups of atoms). Considerations of several bodies simultaneously are generally used to characterize bonded structural elements or those in close contact with each other, but historically do not consider atoms that are not in direct contact with each other. In this report, we introduce an information-theoretic method for detecting and quantifying distance-dependent through-space multibody relationships between the sidechains of three residues. The technique introduced is capable of producing convergent and consistent results when applied to a sufficiently large database of randomly chosen, experimentally solved protein structures. The results of our study can be shown to reproduce established physico-chemical properties of residues as well as more recently discovered properties and interactions. These results offer insight into the numerous roles that residues play in protein structure, as well as relationships between residue function, protein structure, and evolution. The techniques and insights presented in this work should be useful in the future development of novel knowledge-based tools for the evaluation of protein structure.

PONDEROSA-C/S: client-server based software package for automated protein 3D structure determination

Peak-picking Of Noe Data Enabled by Restriction Of Shift Assignments-Client Server (PONDEROSA-C/S) builds on the original PONDEROSA software (Lee et al. in Bioinformatics 27:1727-1728. doi: 10.1093/bioinformatics/btr200, 2011) and includes improved features for structure calculation and refinement. PONDEROSA-C/S consists of three programs: Ponderosa Server, Ponderosa Client, and Ponderosa Analyzer. PONDEROSA-C/S takes as input the protein sequence, a list of assigned chemical shifts, and nuclear Overhauser data sets ((13)C- and/or (15)N-NOESY). The output is a set of assigned NOEs and 3D structural models for the protein. Ponderosa Analyzer supports the visualization, validation, and refinement of the results from Ponderosa Server. These tools enable semi-automated NMR-based structure determination of proteins in a rapid and robust fashion. We present examples showing the use of PONDEROSA-C/S in solving structures of four proteins: two that enable comparison with the original PONDEROSA package, and two from the Critical Assessment of automated Structure Determination by NMR (Rosato et al. in Nat Methods 6:625-626. doi: 10.1038/nmeth0909-625 , 2009) competition. The software package can be downloaded freely in binary format from http://pine.nmrfam.wisc.edu/download_packages.html. Registered users of the National Magnetic Resonance Facility at Madison can submit jobs to the PONDEROSA-C/S server at http://ponderosa.nmrfam.wisc.edu, where instructions, tutorials, and instructions can be found. Structures are normally returned within 1-2 days.

Structure of the membrane protein MerF, a bacterial mercury transporter, improved by the inclusion of chemical shift anisotropy constraints

MerF is a mercury transport membrane protein from the bacterial mercury detoxification system. By performing a solid-state INEPT experiment and measuring chemical shift anisotropy frequencies in aligned samples, we are able to improve on the accuracy and precision of the initial structure that we presented. MerF has four N-terminal and eleven C-terminal residues that are mobile and unstructured in phospholipid bilayers. The structure presented here has average pairwise RMSDs of 1.78 Å for heavy atoms and 0.92 Å for backbone atoms.

Using small angle solution scattering data in Xplor-NIH structure calculations

This contribution describes the use of small and wide angle X-ray and small angle neutron scattering for biomolecular structure calculation using the program Xplor-NIH, both with and without NMR data. The current algorithms used for calculating scattering curves are described, and the use of scattering data as a structural restraint is given concrete form as a fragment of an Xplor-NIH structure calculation script. We review five examples of the use of scattering data in structure calculation, including the treatment of single domain proteins, nucleic acids, structure determination of large proteins, and the use of ensemble representations to characterize small and large amplitude motions.

A practical implicit solvent potential for NMR structure calculation

The benefits of protein structure refinement in water are well documented. However, performing structure refinement with explicit atomic representation of the solvent molecules is computationally expensive and impractical for NMR-restrained structure calculations that start from completely extended polypeptide templates. Here we describe a new implicit solvation potential, EEFx (Effective Energy Function for XPLOR-NIH), for NMR-restrained structure calculations of proteins in XPLOR-NIH. The key components of EEFx are an energy term for solvation energy that works together with other nonbonded energy functions, and a dedicated force field for conformational and nonbonded protein interaction parameters. The initial results obtained with EEFx show that significant improvements in structural quality can be obtained. EEFx is computationally efficient and can be used both to fold and refine structures. Overall, EEFx improves the quality of protein conformation and nonbonded atomic interactions. Moreover, such benefits are accompanied by enhanced structural precision and enhanced structural accuracy, reflected in improved agreement with the cross-validated dipolar coupling data. Finally, implementation of EEFx calculations is straightforward and computationally efficient. Overall, EEFx provides a useful method for the practical calculation of experimental protein structures in a physically realistic environment.

Structure of the Na,K-ATPase regulatory protein FXYD2b in micelles: implications for membrane-water interfacial arginines

FXYD2 is a membrane protein responsible for regulating the function of the Na,K-ATPase in mammalian kidney epithelial cells. Here we report the structure of FXYD2b, one of two splice variants of the protein, determined by NMR spectroscopy in detergent micelles. Solid-state NMR characterization of the protein embedded in phospholipid bilayers indicates that several arginine side chains may be involved in hydrogen bond interactions with the phospholipid polar head groups. The structure and the NMR data suggest that FXYD2b could regulate the Na,K-ATPase by modulating the effective membrane surface electrostatics near the ion binding sites of the pump.

Characterization of the Raptor/4E-BP1 interaction by chemical cross-linking coupled with mass spectrometry analysis

mTORC1 plays critical roles in the regulation of protein synthesis, growth, and proliferation in response to nutrients, growth factors, and energy conditions. One of the substrates of mTORC1 is 4E-BP1, whose phosphorylation by mTORC1 reverses its inhibitory action on eIF4E, resulting in the promotion of protein synthesis. Raptor in mTOR complex 1 is believed to recruit 4E-BP1, facilitating phosphorylation of 4E-BP1 by the kinase mTOR. We applied chemical cross-linking coupled with mass spectrometry analysis to gain insight into interactions between mTORC1 and 4E-BP1. Using the cross-linking reagent bis[sulfosuccinimidyl] suberate, we showed that Raptor can be cross-linked with 4E-BP1. Mass spectrometric analysis of cross-linked Raptor-4E-BP1 led to the identification of several cross-linked peptide pairs. Compilation of these peptides revealed that the most N-terminal Raptor N-terminal conserved domain (in particular residues from 89 to 180) of Raptor is the major site of interaction with 4E-BP1. On 4E-BP1, we found that cross-links with Raptor were clustered in the central region (amino acid residues 56-72) we call RCR (Raptor cross-linking region). Intramolecular cross-links of Raptor suggest the presence of two structured regions of Raptor: one in the N-terminal region and the other in the C-terminal region. In support of the idea that the Raptor N-terminal conserved domain and the 4E-BP1 central region are closely located, we found that peptides that encompass the RCR of 4E-BP1 inhibit cross-linking and interaction of 4E-BP1 with Raptor. Furthermore, mutations of residues in the RCR decrease the ability of 4E-BP1 to serve as a substrate for mTORC1 in vitro and in vivo.

Structure and dynamics of full-length HIV-1 capsid protein in solution

The HIV-1 capsid protein plays a crucial role in viral infectivity, assembling into a cone that encloses the viral RNA. In the mature virion, the N-terminal domain of the capsid protein forms hexameric and pentameric rings, while C-terminal domain homodimers connect adjacent N-terminal domain rings to one another. Structures of disulfide-linked hexamer and pentamer assemblies, as well as structures of the isolated domains, have been solved previously. The dimer configuration in C-terminal domain constructs differs in solution (residues 144-231) and crystal (residues 146-231) structures by ∼30°, and it has been postulated that the former connects the hexamers while the latter links pentamers to hexamers. Here we study the structure and dynamics of full-length capsid protein in solution, comprising a mixture of monomeric and dimeric forms in dynamic equilibrium, using ensemble simulated annealing driven by experimental NMR residual dipolar couplings and X-ray scattering data. The complexity of the system necessitated the development of a novel computational framework that should be generally applicable to many other challenging systems that currently escape structural characterization by standard application of mainstream techniques of structural biology. We show that the orientation of the C-terminal domains in dimeric full-length capsid and isolated C-terminal domain constructs is the same in solution, and we obtain a quantitative description of the conformational space sampled by the N-terminal domain relative to the C-terminal domain on the nano- to millisecond time scale. The positional distribution of the N-terminal domain relative to the C-terminal domain is large and modulated by the oligomerization state of the C-terminal domain. We also show that a model of the hexamer/pentamer assembly can be readily generated with a single configuration of the C-terminal domain dimer, and that capsid assembly likely proceeds via conformational selection of sparsely populated configurations of the N-terminal domain within the capsid protein dimer.

NMR structure of calmodulin complexed to an N-terminally acetylated α-synuclein peptide

Calmodulin (CaM) is a calcium binding protein that plays numerous roles in Ca-dependent cellular processes, including uptake and release of neurotransmitters in neurons. α-Synuclein (α-syn), one of the most abundant proteins in central nervous system neurons, helps maintain presynaptic vesicles containing neurotransmitters and moderates their Ca-dependent release into the synapse. Ca-Bound CaM interacts with α-syn most strongly at its N-terminus. The N-terminal region of α-syn is important for membrane binding; thus, CaM could modulate membrane association of α-syn in a Ca-dependent manner. In contrast, Ca-free CaM has negligible interaction. The interaction with CaM leads to significant signal broadening in both CaM and α-syn NMR spectra, most likely due to conformational exchange. The broadening is much reduced when binding a peptide consisting of the first 19 residues of α-syn. In neurons, most α-syn is acetylated at the N-terminus, and acetylation leads to a 10-fold increase in binding strength for the α-syn peptide (KD = 35 ± 10 μM). The N-terminally acetylated peptide adopts a helical structure at the N-terminus with the acetyl group contacting the N-terminal domain of CaM and with less ordered helical structure toward the C-terminus of the peptide contacting the CaM C-terminal domain. Comparison with known structures shows that the CaM/α-syn complex most closely resembles Ca-bound CaM in a complex with an IQ motif peptide. However, a search comparing the α-syn peptide sequence with known CaM targets, including IQ motifs, found no homologies; thus, the N-terminal α-syn CaM binding site appears to be a novel CaM target sequence.

Accurate high-throughput structure mapping and prediction with transition metal ion FRET

Mapping the landscape of a protein's conformational space is essential to understanding its functions and regulation. The limitations of many structural methods have made this process challenging for most proteins. Here, we report that transition metal ion FRET (tmFRET) can be used in a rapid, highly parallel screen, to determine distances from multiple locations within a protein at extremely low concentrations. The distances generated through this screen for the protein maltose binding protein (MBP) match distances from the crystal structure to within a few angstroms. Furthermore, energy transfer accurately detects structural changes during ligand binding. Finally, fluorescence-derived distances can be used to guide molecular simulations to find low energy states. Our results open the door to rapid, accurate mapping and prediction of protein structures at low concentrations, in large complex systems, and in living cells.