Comprehensive automation for NMR structure determination of proteins

NMR Assignment through Linear Programming

Nuclear Magnetic Resonance (NMR) Spectroscopy is the second most used technique (after X-ray crystallography) for structural determination of proteins. A computational challenge in this technique involves solving a discrete optimization problem that assigns the resonance frequency to each atom in the protein. This paper introduces LIAN (LInear programming Assignment for NMR), a novel linear programming formulation of the problem which yields state-of-the-art results in simulated and experimental datasets.

An engineered construct of cFLIP provides insight into DED1 structure and interactions

Cellular FLICE-like inhibitory protein (cFLIP) is a member of the Death Domain superfamily with pivotal roles in many cellular processes and disease states, including cancer and autoimmune disorders. In the context of the death-inducing signaling complex (DISC), cFLIP isoforms regulate extrinsic apoptosis by controlling procaspase-8 activation. The function of cFLIP is mediated through a series of protein-protein interactions, engaging the two N-terminal death effector domains (DEDs). Here, we solve the structure of an engineered DED1 domain of cFLIP using solution nuclear magnetic resonance (NMR) and we define the interaction with FADD and calmodulin, protein-protein interactions that regulate the function of cFLIP in the DISC. cFLIP DED1 assumes a canonical DED fold characterized by six α helices and is able to bind calmodulin and FADD through two separate interfaces. Our results clearly demonstrate the role of DED1 in the cFLIP/FADD association and contribute to the understanding of the assembly of DISC filaments.

Assessment of prediction methods for protein structures determined by NMR in CASP14: Impact of AlphaFold2

NMR studies can provide unique information about protein conformations in solution. In CASP14, three reference structures provided by solution NMR methods were available (T1027, T1029, and T1055), as well as a fourth data set of NMR‐derived contacts for an integral membrane protein (T1088). For the three targets with NMR‐based structures, the best prediction results ranged from very good (GDT_TS = 0.90, for T1055) to poor (GDT_TS = 0.47, for T1029). We explored the basis of these results by comparing all CASP14 prediction models against experimental NMR data. For T1027, NMR data reveal extensive internal dynamics, presenting a unique challenge for protein structure prediction methods. The analysis of T1029 motivated exploration of a novel method of “inverse structure determination,” in which an AlphaFold2 model was used to guide NMR data analysis. NMR data provided to CASP predictor groups for target T1088, a 238‐residue integral membrane porin, was also used to assess several NMR‐assisted prediction methods. Most groups involved in this exercise generated similar beta‐barrel models, with good agreement with the experimental data. However, as was also observed in CASP13, some pure prediction groups that did not use any NMR data generated models for T1088 that better fit the NMR data than the models generated using these experimental data. These results demonstrate the remarkable power of modern methods to predict structures of proteins with accuracies rivaling solution NMR structures, and that it is now possible to reliably use prediction models to guide and complement experimental NMR data analysis.

Structural Insights into the Mechanism of Human T-cell Leukemia Virus Type 1 Gag Targeting to the Plasma Membrane for Assembly

Retroviral Gag targeting to the plasma membrane (PM) for assembly is mediated by the N-terminal matrix (MA) domain. For many retroviruses, Gag-PM interaction is dependent on phosphatidylinositol 4,5-bisphosphate (PI(4,5)P₂). However, it has been shown that for human T-cell leukemia virus type 1 (HTLV-1), Gag binding to membranes is less dependent on PI(4,5)P₂ than HIV-1, suggesting that other factors may modulate Gag assembly. To elucidate the mechanism by which HTLV-1 Gag binds to the PM, we employed NMR techniques to determine the structure of unmyristoylated MA (myr(-)MA) and to characterize its interactions with lipids and liposomes. The MA structure consists of four α-helices and unstructured N- and C-termini. We show that myr(-)MA binds to PI(4,5)P₂ via the polar head and that binding to inositol phosphates (IPs) is significantly enhanced by increasing the number of phosphate groups on the inositol ring, indicating that the MA-IP binding is governed by charge-charge interactions. The IP binding site was mapped to a well-defined basic patch formed by lysine and arginine residues. Using an NMR-based liposome binding assay, we show that PI(4,5)P₂and phosphatidylserine enhance myr(-)MA binding in a synergistic fashion. Confocal microscopy data revealed formation of puncta on the PM of Gag expressing cells. However, G2A-Gag mutant, lacking myristoylation, is diffuse and cytoplasmic. These results suggest that although myr(-)MA binds to membranes, myristoylation appears to be key for formation of HTLV-1 Gag puncta on the PM. Altogether, these findings advance our understanding of a key mechanism in retroviral assembly.

Structures of Pathological and Functional Amyloids and Prions, a Solid-State NMR Perspective

Infectious proteins or prions are a remarkable class of pathogens, where pathogenicity and infectious state correspond to conformational transition of a protein fold. The conformational change translates into the formation by the protein of insoluble amyloid aggregates, associated in humans with various neurodegenerative disorders and systemic protein-deposition diseases. The prion principle, however, is not limited to pathogenicity. While pathological amyloids (and prions) emerge from protein misfolding, a class of functional amyloids has been defined, consisting of amyloid-forming domains under natural selection and with diverse biological roles. Although of great importance, prion amyloid structures remain challenging for conventional structural biology techniques. Solid-state nuclear magnetic resonance (SSNMR) has been preferentially used to investigate these insoluble, morphologically heterogeneous aggregates with poor crystallinity. SSNMR methods have yielded a wealth of knowledge regarding the fundamentals of prion biology and have helped to solve the structures of several prion and prion-like fibrils. Here, we will review pathological and functional amyloid structures and will discuss some of the obtained structural models. We will finish the review with a perspective on integrative approaches combining solid-state NMR, electron paramagnetic resonance and cryo-electron microscopy, which can complement and extend our toolkit to structurally explore various facets of prion biology.

Solution Structure of the C-terminal Domain of A20, the Missing Brick for the Characterization of the Interface between Vaccinia Virus DNA Polymerase and its Processivity Factor

Poxviruses are enveloped viruses with a linear, double-stranded DNA genome. Viral DNA synthesis is achieved by a functional DNA polymerase holoenzyme composed of three essential proteins. For vaccinia virus (VACV) these are E9, the catalytic subunit, a family B DNA polymerase, and the heterodimeric processivity factor formed by D4 and A20. The A20 protein links D4 to the catalytic subunit. High-resolution structures have been obtained for the VACV D4 protein in complex with an N-terminal fragment of A20 as well as for E9. In addition, biochemical studies provided evidence that a poxvirus-specific insertion (insert 3) in E9 interacts with the C-terminal residues of A20. Here, we provide solution structures of two different VACV A20 C-terminal constructs containing residues 304-426, fused at their C-terminus to either a BAP (Biotin Acceptor Peptide)-tag or a short peptide containing the helix of E9 insert 3. Together with results from titration studies, these structures shed light on the molecular interface between the catalytic subunit and the processivity factor component A20. The interface comprises hydrophobic residues conserved within the Chordopoxvirinae subfamily. Finally, we constructed a HADDOCK model of the VACV A20_304-426-E9 complex, which is in excellent accordance with previous experimental data.

A pleurocidin analogue with greater conformational flexibility, enhanced antimicrobial potency and in vivo therapeutic efficacy

Antimicrobial peptides (AMPs) are a potential alternative to classical antibiotics that are yet to achieve a therapeutic breakthrough for treatment of systemic infections. The antibacterial potency of pleurocidin, an AMP from Winter Flounder, is linked to its ability to cross bacterial plasma membranes and seek intracellular targets while also causing membrane damage. Here we describe modification strategies that generate pleurocidin analogues with substantially improved, broad spectrum, antibacterial properties, which are effective in murine models of bacterial lung infection. Increasing peptide-lipid intermolecular hydrogen bonding capabilities enhances conformational flexibility, associated with membrane translocation, but also membrane damage and potency, most notably against Gram-positive bacteria. This negates their ability to metabolically adapt to the AMP threat. An analogue comprising D-amino acids was well tolerated at an intravenous dose of 15 mg/kg and similarly effective as vancomycin in reducing EMRSA-15 lung CFU. This highlights the therapeutic potential of systemically delivered, bactericidal AMPs.

ARIAweb: a server for automated NMR structure calculation

Nuclear magnetic resonance (NMR) spectroscopy is a method of choice to study the dynamics and determine the atomic structure of macromolecules in solution. The standalone program ARIA (Ambiguous Restraints for Iterative Assignment) for automated assignment of nuclear Overhauser enhancement (NOE) data and structure calculation is well established in the NMR community. To ultimately provide a perfectly transparent and easy to use service, we designed an online user interface to ARIA with additional functionalities. Data conversion, structure calculation setup and execution, followed by interactive visualization of the generated 3D structures are all integrated in ARIAweb and freely accessible at https://ariaweb.pasteur.fr.

Revealing the mechanism of repressor inactivation during switching of a temperate bacteriophage

Temperate bacteriophages can enter one of two life cycles following infection of a sensitive host: the lysogenic or the lytic life cycle. The choice between the two alternative life cycles is dependent upon a tight regulation of promoters and their cognate regulatory proteins within the phage genome. We investigated the genetic switch of TP901-1, a bacteriophage of Lactococcus lactis, controlled by the CI repressor and the modulator of repression (MOR) antirepressor and their interactions with DNA. We determined the solution structure of MOR, and we solved the crystal structure of MOR in complex with the N-terminal domain of CI, revealing the structural basis of MOR inhibition of CI binding to the DNA operator sites. ¹⁵N NMR Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion and rotating frame R _1ρ measurements demonstrate that MOR displays molecular recognition dynamics on two different time scales involving a repacking of aromatic residues at the interface with CI. Mutations in the CI:MOR binding interface impair complex formation in vitro, and when introduced in vivo, the bacteriophage switch is unable to choose the lytic life cycle showing that the CI:MOR complex is essential for proper functioning of the genetic switch. On the basis of sequence alignments, we show that the structural features of the MOR:CI complex are likely conserved among a larger family of bacteriophages from human pathogens implicated in transfer of antibiotic resistance.

A protocol to automatically calculate homo-oligomeric protein structures through the integration of evolutionary constraints and NMR ambiguous contacts

Protein assemblies are involved in many important biological processes. Solid-state NMR (SSNMR) spectroscopy is a technique suitable for the structural characterization of samples with high molecular weight and thus can be applied to such assemblies. A significant bottleneck in terms of both effort and time required is the manual identification of unambiguous intermolecular contacts. This is particularly challenging for homo-oligomeric complexes, where simple uniform labeling may not be effective. We tackled this challenge by exploiting coevolution analysis to extract information on homo-oligomeric interfaces from NMR-derived ambiguous contacts. After removing the evolutionary couplings (ECs) that are already satisfied by the 3D structure of the monomer, the predicted ECs are matched with the automatically generated list of experimental contacts. This approach provides a selection of potential interface residues that is used directly in monomer-monomer docking calculations. We validated the protocol on tetrameric L-asparaginase II and dimeric Sod1.

Intrinsic Instability of BOK Enables Membrane Permeabilization in Apoptosis

The effector B cell lymphoma-2 (BCL-2) protein BCL-2 ovarian killer (BOK) induces mitochondrial outer membrane permeabilization (MOMP) to initiate apoptosis upon inhibition of the proteasome. How BOK mediates MOMP is mechanistically unknown. The NMR structure of the BCL-2 core of human BOK reveals a conserved architecture with an atypical hydrophobic groove that undergoes conformational exchange. Remarkably, the BCL-2 core of BOK spontaneously associates with purified mitochondria to release cytochrome c in MOMP assays. Alanine substitution of a unique glycine in helix α1 stabilizes BOK, as shown by thermal shift and urea denaturation analyses, and significantly inhibits MOMP, liposome permeabilization, and cell death. Activated BID does not activate WT BOK or the stabilized alanine mutant to promote cell death. We propose that BOK-mediated membrane permeabilization is governed in part by its unique metastability of the hydrophobic groove and helix α1 and not through activation by BH3 ligands.

Systematic Exploration of Protein Conformational Space Using a Distance Geometry Approach

The optimization approaches classically used during the determination of protein structure encounter various difficulties, especially when the size of the conformational space is large. Indeed, in such a case, algorithmic convergence criteria are more difficult to set up. Moreover, the size of the search space makes it difficult to achieve a complete exploration. The interval branch-and-prune (iBP) approach, based on the reformulation of the distance geometry problem (DGP) provides a theoretical frame for the generation of protein conformations, by systematically sampling the conformational space. When an appropriate subset of interatomic distances is known exactly, this worst-case exponential-time algorithm is provably complete and fixed-parameter tractable. These guarantees, however, immediately disappear as distance measurement errors are introduced. Here we propose an improvement of this approach: threading-augmented interval branch-and-prune (TAiBP), where the combinatorial explosion of the original iBP approach arising from its exponential complexity is alleviated by partitioning the input instances into consecutive peptide fragments and by using self-organizing maps (SOMs) to obtain clusters of similar solutions. A validation of the TAiBP approach is presented here on a set of proteins of various sizes and structures. The calculation inputs are a uniform covalent geometry extracted from force field covalent terms, the backbone dihedral angles with error intervals, and a few long-range distances. For most of the proteins smaller than 50 residues and interval widths of 20°, the TAiBP approach yielded solutions with RMSD values smaller than 3 Å with respect to the initial protein conformation. The efficiency of the TAiBP approach for proteins larger than 50 residues will require the use of nonuniform covalent geometry and may have benefits from the recent development of residue-specific force-fields.

Structural characterization of bacteriophage viruses by NMR

Magic-angle spinning (MAS) solid-state NMR has provided structural insights into various bacteriophage systems including filamentous, spherical, and tailed bacteriophage viruses. A variety of methodologies have been utilized including elementary two and three-dimensional assignment experiments, proton-detection techniques at fast spinning speeds, non-uniform sampling, structure determination protocols, conformational dynamics revealed by recoupling of anisotropic interactions, and enhancement by dynamic nuclear polarization. This review summarizes most of the studies performed during the last decade by MAS techniques and makes comparisons with prior knowledge obtained from static and solution NMR techniques. Chemical shifts for the capsids of the various systems are reported and analyzed, and DNA shifts are reported and discussed in the context of general high molecular-weight DNA molecules. Chemical shift and torsion angle prediction techniques are compared and applied to the various phage systems. The structures of the intact M13 filamentous bacteriophage and that of the Acinetobacter phage AP205 capsid, determined using MAS-based experimental data, are presented. Finally, filamentous phages, which are highly rigid systems, show interesting dynamics at the interface of the capsid and DNA, and their mutual electrostatic interactions are shown to be mediated by highly mobile positively charged residues. Novel results obtained from recoupling the chemical shift anisotropy of a single arginine in IKe phage, which is in contact with its DNA, further demonstrate this point. MAS NMR thus provides many new insights into phage structure, and on the other hand the richness, complexity and variety of bacteriophage systems provide opportunities for new NMR methodologies and technique developments.

Temporin L and aurein 2.5 have identical conformations but subtly distinct membrane and antibacterial activities

Frogs such as Rana temporaria and Litoria aurea secrete numerous closely related antimicrobial peptides (AMPs) as an effective chemical dermal defence. Damage or penetration of the bacterial plasma membrane is considered essential for AMP activity and such properties are commonly ascribed to their ability to form secondary amphipathic, α-helix conformations in membrane mimicking milieu. Nevertheless, despite the high similarity in physical properties and preference for adopting such conformations, the spectrum of activity and potency of AMPs often varies considerably. Hence distinguishing apparently similar AMPs according to their behaviour in, and effects on, model membranes will inform understanding of primary-sequence-specific antimicrobial mechanisms. Here we use a combination of molecular dynamics simulations, circular dichroism and patch-clamp to investigate the basis for differing anti-bacterial activities in representative AMPs from each species; temporin L and aurein 2.5. Despite adopting near identical, α-helix conformations in the steady-state in a variety of membrane models, these two AMPs can be distinguished both in vitro and in silico based on their dynamic interactions with model membranes, notably their differing conformational flexibility at the N-terminus, ability to form higher order aggregates and the characteristics of induced ion conductance. Taken together, these differences provide an explanation of the greater potency and broader antibacterial spectrum of activity of temporin L over aurein 2.5. Consequently, while the secondary amphipathic, α-helix conformation is a key determinant of the ability of a cationic AMP to penetrate and disrupt the bacterial plasma membrane, the exact mechanism, potency and spectrum of activity is determined by precise structural and dynamic contributions from specific residues in each AMP sequence.

Minor sequence modifications in temporin B cause drastic changes in antibacterial potency and selectivity by fundamentally altering membrane activity

Antimicrobial peptides (AMPs) are a potential source of new molecules to counter the increase in antimicrobial resistant infections but a better understanding of their properties is required to understand their native function and for effective translation as therapeutics. Details of the mechanism of their interaction with the bacterial plasma membrane are desired since damage or penetration of this structure is considered essential for AMPs activity. Relatively modest modifications to AMPs primary sequence can induce substantial changes in potency and/or spectrum of activity but, hitherto, have not been predicted to substantially alter the mechanism of interaction with the bacterial plasma membrane. Here we use a combination of molecular dynamics simulations, circular dichroism, solid-state NMR and patch clamp to investigate the extent to which temporin B and its analogues can be distinguished both in vitro and in silico on the basis of their interactions with model membranes. Enhancing the hydrophobicity of the N-terminus and cationicity of the C-terminus in temporin B improves its membrane activity and potency against both Gram-negative and Gram-positive bacteria. In contrast, enhancing the cationicity of the N-terminus abrogates its ability to trigger channel conductance and renders it ineffective against Gram-positive bacteria while nevertheless enhancing its potency against Escherichia coli. Our findings suggest even closely related AMPs may target the same bacterium with fundamentally differing mechanisms of action.

Induced conformational changes activate the peptidoglycan synthase PBP1B

Bacteria surround their cytoplasmic membrane with an essential, stress‐bearing peptidoglycan (PG) layer consisting of glycan chains linked by short peptides into a mesh‐like structure. Growing and dividing cells expand their PG layer using inner‐membrane anchored PG synthases, including Penicillin‐binding proteins (PBPs), which participate in dynamic protein complexes to facilitate cell wall growth. In Escherichia coli, and presumably other Gram‐negative bacteria, growth of the mainly single layered PG is regulated by outer membrane‐anchored lipoproteins. The lipoprotein LpoB is required to activate PBP1B, which is a major, bi‐functional PG synthase with glycan chain polymerising (glycosyltransferase) and peptide cross‐linking (transpeptidase) activities. In this work we show how the binding of LpoB to the regulatory UB2H domain of PBP1B activates both activities. Binding induces structural changes in the UB2H domain, which transduce to the two catalytic domains by distinct allosteric pathways. We also show how an additional regulator protein, CpoB, is able to selectively modulate the TPase activation by LpoB without interfering with GTase activation.

Structural basis for targeting avian sarcoma virus Gag polyprotein to the plasma membrane for virus assembly

For most retroviruses, including HIV-1, binding of the Gag polyprotein to the plasma membrane (PM) is mediated by interactions between Gag's N-terminal myristoylated matrix (MA) domain and phosphatidylinositol 4,5-bisphosphate (PI(4,5)P₂) in the PM. The Gag protein of avian sarcoma virus (ASV) lacks the N-myristoylation signal but contains structural domains having functions similar to those of HIV-1 Gag. The molecular mechanism by which ASV Gag binds to the PM is incompletely understood. Here, we employed NMR techniques to elucidate the molecular determinants of the membrane-binding domain of ASV MA (MA₈₇) to lipids and liposomes. We report that MA₈₇ binds to the polar head of phosphoinositides such as PI(4,5)P₂ We found that MA₈₇ binding to inositol phosphates (IPs) is significantly enhanced by increasing the number of phosphate groups, indicating that the MA₈₇-IP binding is governed by charge-charge interactions. Using a sensitive NMR-based liposome-binding assay, we show that binding of MA₈₇ to liposomes is enhanced by incorporation of PI(4,5)P₂ and phosphatidylserine. We also show that membrane binding is mediated by a basic surface formed by Lys-6, Lys-13, Lys-23, and Lys-24. Substitution of these residues to glutamate abolished binding of MA₈₇ to both IPs and liposomes. In an accompanying paper, we further report that mutation of these lysine residues diminishes Gag assembly on the PM and inhibits ASV particle release. These findings provide a molecular basis for ASV Gag binding to the inner leaflet of the PM and advance our understanding of the basic mechanisms of retroviral assembly.

A histidine-rich Pseudomonas metallothionein with a disordered tail displays higher binding capacity for cadmium than zinc

The NMR solution structure of a Pseudomonas metallothionein reveals a different binding capacity for Zn^II and Cd^II ions that results in two novel metal-cluster topologies. Replacement of a non-coordinating residue by histidine decreases the kinetic lability of the cluster. All three structures reported show an identical protein fold.

Atomic Scale Structural Studies of Macromolecular Assemblies by Solid-state Nuclear Magnetic Resonance Spectroscopy

Supramolecular protein assemblies play fundamental roles in biological processes ranging from host-pathogen interaction, viral infection to the propagation of neurodegenerative disorders. Such assemblies consist in multiple protein subunits organized in a non-covalent way to form large macromolecular objects that can execute a variety of cellular functions or cause detrimental consequences. Atomic insights into the assembly mechanisms and the functioning of those macromolecular assemblies remain often scarce since their inherent insolubility and non-crystallinity often drastically reduces the quality of the data obtained from most techniques used in structural biology, such as X-ray crystallography and solution Nuclear Magnetic Resonance (NMR). We here present magic-angle spinning solid-state NMR spectroscopy (SSNMR) as a powerful method to investigate structures of macromolecular assemblies at atomic resolution. SSNMR can reveal atomic details on the assembled complex without size and solubility limitations. The protocol presented here describes the essential steps from the production of 13C/15N isotope-labeled macromolecular protein assemblies to the acquisition of standard SSNMR spectra and their analysis and interpretation. As an example, we show the pipeline of a SSNMR structural analysis of a filamentous protein assembly.

A Residue Specific Insight into the Arkadia E3 Ubiquitin Ligase Activity and Conformational Plasticity

Arkadia (Rnf111) is an E3 ubiquitin ligase that plays a central role in the amplification of transforming growth factor beta (TGF-β) signaling responses by targeting for degradation the negative regulators of the pathway, Smad6 and Smad7, and the nuclear co-repressors Ski and Skil (SnoN). Arkadia's function in vivo depends on the really interesting new gene (RING)-H2 interaction with the E2 enzyme UbcH5b in order to ligate ubiquitin chains on its substrates. A conserved tryptophan (W972) in the C-terminal α-helix is widely accepted as essential for E2 recruitment and interaction and thus also for E3 enzymatic activity. The present NMR-driven study provides an atomic-level investigation of the structural and dynamical properties of two W972 Arkadia RING mutants, attempting to illuminate for the first time the differences between a functional and a nonfunctional mutant W972A and W972R, respectively. A TGF-β-responsive promoter driving luciferase was used to assay for Arkadia function in vivo. These experiments showed that the Arkadia W972A mutant has the same activity as wild-type (WT) Arkadia in enhancing TGF-β signaling responses, while W972R does not. Only minor structural differences exist between the W972A RING domain and WT-RING. In contrast, the W972R mutant hardly interacts with E2. The loss of function correlates with structural changes in the C-terminal α-helix and an increase in the distance between the Zn(II) ions. Our data show that the position occupied by W972 within WT Arkadia is critical for the function of RING and that it depends on the nature of the residue at this position.

Bacterial Filamentous Appendages Investigated by Solid-State NMR Spectroscopy

The assembly of filamentous appendages at the surface of bacteria is essential in many infection mechanisms. The extent of mechanical, dynamical, and functional properties of such appendages is very diverse, ranging from a structural scaffold of the pathogen-host cell interaction to cell motility, surface adhesion, or the export of virulence effectors. In particular, the architectures of several bacterial secretion systems have revealed the presence of filamentous architectures, known as pili, fimbriae, andneedles. At the macroscopic level, filamentous bacterial appendages appear as thin extracellular filaments of several nanometers in diameter and up to several microns in length. The structural characterization of these appendages at atomic-scale resolution represents an extremely challenging task because of their inherent noncrystallinity and very poor solubility. Here, we describe protocols based on recent advances in solid-state NMR spectroscopy to investigate the secondary structure, subunit-subunit protein interactions, symmetry parameters, and atomic architecture of bacterial filaments.

Structure of fully protonated proteins by proton-detected magic-angle spinning NMR

Protein structure determination by proton-detected magic-angle spinning (MAS) NMR has focused on highly deuterated samples, in which only a small number of protons are introduced and observation of signals from side chains is extremely limited. Here, we show in two fully protonated proteins that, at 100-kHz MAS and above, spectral resolution is high enough to detect resolved correlations from amide and side-chain protons of all residue types, and to reliably measure a dense network of (1)H-(1)H proximities that define a protein structure. The high data quality allowed the correct identification of internuclear distance restraints encoded in 3D spectra with automated data analysis, resulting in accurate, unbiased, and fast structure determination. Additionally, we find that narrower proton resonance lines, longer coherence lifetimes, and improved magnetization transfer offset the reduced sample size at 100-kHz spinning and above. Less than 2 weeks of experiment time and a single 0.5-mg sample was sufficient for the acquisition of all data necessary for backbone and side-chain resonance assignment and unsupervised structure determination. We expect the technique to pave the way for atomic-resolution structure analysis applicable to a wide range of proteins.

Structure-function analysis of the extracellular domain of the pneumococcal cell division site positioning protein MapZ

Accurate placement of the bacterial division site is a prerequisite for the generation of two viable and identical daughter cells. In Streptococcus pneumoniae, the positive regulatory mechanism involving the membrane protein MapZ positions precisely the conserved cell division protein FtsZ at the cell centre. Here we characterize the structure of the extracellular domain of MapZ and show that it displays a bi-modular structure composed of two subdomains separated by a flexible serine-rich linker. We further demonstrate in vivo that the N-terminal subdomain serves as a pedestal for the C-terminal subdomain, which determines the ability of MapZ to mark the division site. The C-terminal subdomain displays a patch of conserved amino acids and we show that this patch defines a structural motif crucial for MapZ function. Altogether, this structure–function analysis of MapZ provides the first molecular characterization of a positive regulatory process of bacterial cell division.

Placement of the bacterial division site is crucial for the creation of identical daughter cells. Here, the authors solve the structure of the MapZ protein, which helps to position the cell division protein FtsZ at the cell centre, and further analyse the function of the protein in vivo.

Structural biology computing: Lessons for the biomedical research sciences

The field of structural biology, whose aim is to elucidate the molecular and atomic structures of biological macromolecules, has long been at the forefront of biomedical sciences in adopting and developing computational research methods. Operating at the intersection between biophysics, biochemistry, and molecular biology, structural biology's growth into a foundational framework on which many concepts and findings of molecular biology are interpreted1 has depended largely on parallel advancements in computational tools and techniques. Without these computing advances, modern structural biology would likely have remained an exclusive pursuit practiced by few, and not become the widely practiced, foundational field it is today. As other areas of biomedical research increasingly embrace research computing techniques, the successes, failures and lessons of structural biology computing can serve as a useful guide to progress in other biomedically related research fields.

Structure determination of helical filaments by solid-state NMR spectroscopy

The controlled formation of filamentous protein complexes plays a crucial role in many biological systems and represents an emerging paradigm in signal transduction. The mitochondrial antiviral signaling protein (MAVS) is a central signal transduction hub in innate immunity that is activated by a receptor-induced conversion into helical superstructures (filaments) assembled from its globular caspase activation and recruitment domain. Solid-state NMR (ssNMR) spectroscopy has become one of the most powerful techniques for atomic resolution structures of protein fibrils. However, for helical filaments, the determination of the correct symmetry parameters has remained a significant hurdle for any structural technique and could thus far not be precisely derived from ssNMR data. Here, we solved the atomic resolution structure of helical MAVS(CARD) filaments exclusively from ssNMR data. We present a generally applicable approach that systematically explores the helical symmetry space by efficient modeling of the helical structure restrained by interprotomer ssNMR distance restraints. Together with classical automated NMR structure calculation, this allowed us to faithfully determine the symmetry that defines the entire assembly. To validate our structure, we probed the protomer arrangement by solvent paramagnetic resonance enhancement, analysis of chemical shift differences relative to the solution NMR structure of the monomer, and mutagenesis. We provide detailed information on the atomic contacts that determine filament stability and describe mechanistic details on the formation of signaling-competent MAVS filaments from inactive monomers.

The Solution Structure and Dynamics of Full-length Human Cerebral Dopamine Neurotrophic Factor and Its Neuroprotective Role against α-Synuclein Oligomers

Cerebral dopamine neurotrophic factor (CDNF) is a promising therapeutic agent for Parkinson disease. As such, there has been great interest in studying its mode of action, which remains unknown. The three-dimensional crystal structure of the N terminus (residues 9-107) of CDNF has been determined, but there have been no published structural studies on the full-length protein due to proteolysis of its C-terminal domain, which is considered intrinsically disordered. An improved purification protocol enabled us to obtain active full-length CDNF and to determine its three-dimensional structure in solution. CDNF contains two well folded domains (residues 10-100 and 111-157) that are linked by a loop of intermediate flexibility. We identified two surface patches on the N-terminal domain that were characterized by increased conformational dynamics that should allow them to embrace active sites. One of these patches is formed by residues Ser-33, Leu-34, Ala-66, Lys-68, Ile-69, Leu-70, Ser-71, and Glu-72. The other includes a flexibly disordered N-terminal tail (residues 1-9), followed by the N-terminal portion of α-helix 1 (residues Cys-11, Glu-12, Val-13, Lys-15, and Glu-16) and residue Glu-88. The surface of the C-terminal domain contains two conserved active sites, which have previously been identified in mesencephalic astrocyte-derived neurotrophic factor, a CDNF paralog, which corresponds to its intracellular mode of action. We also showed that CDNF was able to protect dopaminergic neurons against injury caused by α-synuclein oligomers. This advises its use against physiological damages caused by α-synuclein oligomers, as observed in Parkinson disease and several other neurodegenerative diseases.

CASD-NMR 2: robust and accurate unsupervised analysis of raw NOESY spectra and protein structure determination with UNIO

UNIO is a comprehensive software suite for protein NMR structure determination that enables full automation of all NMR data analysis steps involved--including signal identification in NMR spectra, sequence-specific backbone and side-chain resonance assignment, NOE assignment and structure calculation. Within the framework of the second round of the community-wide stringent blind NMR structure determination challenge (CASD-NMR 2), we participated in two categories of CASD-NMR 2, namely using either raw NMR spectra or unrefined NOE peak lists as input. A total of 15 resulting NMR structure bundles were submitted for 9 out of 10 blind protein targets. All submitted UNIO structures accurately coincided with the corresponding blind targets as documented by an average backbone root mean-square deviation to the reference proteins of only 1.2 Å. Also, the precision of the UNIO structure bundles was virtually identical to the ensemble of reference structures. By assessing the quality of all UNIO structures submitted to the two categories, we find throughout that only the UNIO-ATNOS/CANDID approach using raw NMR spectra consistently yielded structure bundles of high quality for direct deposition in the Protein Data Bank. In conclusion, the results obtained in CASD-NMR 2 are another vital proof for robust, accurate and unsupervised NMR data analysis by UNIO for real-world applications.

An algorithm to enumerate all possible protein conformations verifying a set of distance constraints

BACKGROUND

The determination of protein structures satisfying distance constraints is an important problem in structural biology. Whereas the most common method currently employed is simulated annealing, there have been other methods previously proposed in the literature. Most of them, however, are designed to find one solution only.

RESULTS

In order to explore exhaustively the feasible conformational space, we propose here an interval Branch-and-Prune algorithm (iBP) to solve the Distance Geometry Problem (DGP) associated to protein structure determination. This algorithm is based on a discretization of the problem obtained by recursively constructing a search space having the structure of a tree, and by verifying whether the generated atomic positions are feasible or not by making use of pruning devices. The pruning devices used here are directly related to features of protein conformations.

CONCLUSIONS

We described the new algorithm iBP to generate protein conformations satisfying distance constraints, that would potentially allows a systematic exploration of the conformational space. The algorithm iBP has been applied on three α-helical peptides.

NMR-based modeling and refinement of protein 3D structures

NMR is a well-established method to characterize the structure and dynamics of biomolecules in solution. High-quality structures can now be produced thanks to both experimental advances and computational developments that incorporate new NMR parameters and improved protocols and force fields in the structure calculation and refinement process. In this chapter, we give a short overview of the various types of NMR data that can provide structural information, and then focus on the structure calculation methodology itself. We discuss and illustrate with tutorial examples "classical" structure calculation, refinement, and structure validation approaches.

NMR structure validation in relation to dynamics and structure determination

NMR spectroscopy is a key technique for understanding the behaviour of proteins, especially highly dynamic proteins that adopt multiple conformations in solution. Overall, protein structures determined from NMR spectroscopy data constitute just over 10% of the Protein Data Bank archive. This review covers the validation of these NMR protein structures, but rather than describing currently available methodology, it focuses on concepts that are important for understanding where and how validation is most relevant. First, the inherent characteristics of the protein under study have an influence on quality and quantity of the distinct types of data that can be acquired from NMR experiments. Second, these NMR data are necessarily transformed into a model for use in a structure calculation protocol, and the protein structures that result from this reflect the types of NMR data used as well as the protein characteristics. The validation of NMR protein structures should therefore take account, wherever possible, of the inherent behavioural characteristics of the protein, the types of available NMR data, and the calculation protocol. These concepts are discussed in the context of 'knowledge based' and 'model versus data' validation, with suggestions for questions to ask and different validation categories to consider. The principal aim of this review is to stimulate discussion and to help the reader understand the relationships between the above elements in order to make informed decisions on which validation approaches are the most relevant in particular cases.

Elongated structure of the outer-membrane activator of peptidoglycan synthesis LpoA: implications for PBP1A stimulation

The bacterial cell envelope contains the stress-bearing peptidoglycan layer, which is enlarged during cell growth and division by membrane-anchored synthases guided by cytoskeletal elements. In Escherichia coli, the major peptidoglycan synthase PBP1A requires stimulation by the outer-membrane-anchored lipoprotein LpoA. Whereas the C-terminal domain of LpoA interacts with PBP1A to stimulate its peptide crosslinking activity, little is known about the role of the N-terminal domain. Herein we report its NMR structure, which adopts an all-α-helical fold comprising a series of helix-turn-helix tetratricopeptide-repeat (TPR)-like motifs. NMR spectroscopy of full-length LpoA revealed two extended flexible regions in the C-terminal domain and limited, if any, flexibility between the N- and C-terminal domains. Analytical ultracentrifugation and small-angle X-ray scattering results are consistent with LpoA adopting an elongated shape, with dimensions sufficient to span from the outer membrane through the periplasm to interact with the peptidoglycan synthase PBP1A.

•

LpoA’s N-terminal domain features an all-α-helical fold similar to TPR domains

•

The C-terminal domain of E. coli LpoA contains two extensive flexible regions

•

Full-length LpoA adopts an elongated structure with low interdomain flexibility

•

LpoA can span the periplasm to stimulate the peptidoglycan synthase PBP1A

New insights into histidine triad proteins: solution structure of a Streptococcus pneumoniae PhtD domain and zinc transfer to AdcAII

Zinc (Zn²⁺) homeostasis is critical for pathogen host colonization and invasion. Polyhistidine triad (Pht) proteins, located at the surface of various streptococci, have been proposed to be involved in Zn²⁺ homeostasis. The phtD gene, coding for a Zn²⁺-binding protein, is organized in an operon with adcAII coding for the extracellular part of a Zn²⁺ transporter. In the present work, we investigate the relationship between PhtD and AdcAII using biochemical and structural biology approaches. Immuno-precipitation experiments on purified membranes of Streptococcus pneumoniae (S. pneumoniae) demonstrate that native PhtD and AdcAII interact in vivo confirming our previous in vitro observations. NMR was used to demonstrate Zn²⁺ transfer from the Zn²⁺-bound form of a 137 amino acid N-terminal domain of PhtD (t-PhtD) to AdcAII. The high resolution NMR structure of t-PhtD shows that Zn²⁺ is bound in a tetrahedral site by histidines 83, 86, and 88 as well as by glutamate 63. Comparison of the NMR parameters measured for apo- and Zn²⁺-t-PhtD shows that the loss of Zn²⁺ leads to a diminished helical propensity at the C-terminus and increases the local dynamics and overall molecular volume. Structural comparison with the crystal structure of a 55-long fragment of PhtA suggests that Pht proteins are built from short repetitive units formed by three β-strands containing the conserved HxxHxH motif. Taken together, these results support a role for S. pneumoniae PhtD as a Zn²⁺ scavenger for later release to the surface transporter AdcAII, leading to Zn²⁺ uptake.

Iron-sulfur cluster binding by mitochondrial monothiol glutaredoxin-1 of Trypanosoma brucei: molecular basis of iron-sulfur cluster coordination and relevance for parasite infectivity

AIMS

Monothiol glutaredoxins (1-C-Grxs) are small proteins linked to the cellular iron and redox metabolism. Trypanosoma brucei brucei, model organism for human African trypanosomiasis, expresses three 1-C-Grxs. 1-C-Grx1 is a highly abundant mitochondrial protein capable to bind an iron-sulfur cluster (ISC) in vitro using glutathione (GSH) as cofactor. We here report on the functional and structural analysis of 1-C-Grx1 in relation to its ISC-binding properties.

RESULTS

An N-terminal extension unique to 1-C-Grx1 from trypanosomatids affects the oligomeric structure and the ISC-binding capacity of the protein. The active-site Cys104 is essential for ISC binding, and the parasite-specific glutathionylspermidine and trypanothione can replace GSH as the ligands of the ISC. Interestingly, trypanothione forms stable protein-free ISC species that in vitro are incorporated into the dithiol T. brucei 2-C-Grx1, but not 1-C-Grx1. Overexpression of the C104S mutant of 1-C-Grx1 impairs disease progression in a mouse model. The structure of the Grx-domain of 1-C-Grx1 was solved by nuclear magnetic resonance spectroscopy. Despite the fact that several residues--which in other 1-C-Grxs are involved in the noncovalent binding of GSH--are conserved, different physicochemical approaches did not reveal any specific interaction between 1-C-Grx1 and free thiol ligands.

INNOVATION

Parasite Grxs are able to coordinate an ISC formed with trypanothione, suggesting a new mechanism of ISC binding and a novel function for the parasite-specific dithiol. The first 3D structure and in vivo relevance of a 1-C-Grx from a pathogenic protozoan are reported.

CONCLUSION

T. brucei 1-C-Grx1 is indispensable for mammalian parasitism and utilizes a new mechanism for ISC binding.

Mono- and dithiol glutaredoxins in the trypanothione-based redox metabolism of pathogenic trypanosomes

SIGNIFICANCE

Glutaredoxins are ubiquitous small thiol proteins of the thioredoxin-fold superfamily. Two major groups are distinguished based on their active sites: the dithiol (2-C-Grxs) and the monothiol (1-C-Grxs) glutaredoxins with a CXXC and a CXXS active site motif, respectively. Glutaredoxins are involved in cellular redox and/or iron sulfur metabolism. Usually their functions are closely linked to the glutathione system. Trypanosomatids, the causative agents of several tropical diseases, rely on trypanothione as principal low molecular mass thiol, and their glutaredoxins readily react with the unique bis(glutathionyl) spermidine conjugate.

RECENT ADVANCES

Two 2-C-Grxs and three 1-C-Grxs have been identified in pathogenic trypanosomatids. The 2-C-Grxs catalyze the reduction of glutathione disulfide by trypanothione and display reductase activity towards protein disulfides, as well as protein-glutathione mixed disulfides. In vitro, all three 1-C-Grxs as well as the cytosolic 2-C-Grx of Trypanosoma brucei can complex an iron-sulfur cluster. Recently the structure of the 1-C-Grx1 has been solved by NMR spectroscopy. The structure is very similar to those of other 1-C-Grxs, with some differences in the loop containing the conserved cis-Pro and the surface charge distribution.

CRITICAL ISSUES

Although four of the five trypanosomal glutaredoxins proved to coordinate an iron-sulfur cluster in vitro, the physiological role of the mitochondrial and cytosolic proteins, respectively, has only started to be unraveled.

FUTURE DIRECTIONS

The use of trypanothione by the glutaredoxins has established a novel role for this parasite-specific dithiol. Future work should reveal if these differences can be exploited for the development of novel antiparasitic drugs.

MOTOR: model assisted software for NMR structure determination

Eukaryotic proteins with important biological function can be partially unstructured, conformational flexible, or heterogenic. Crystallization trials often fail for such proteins. In NMR spectroscopy, parts of the polypeptide chain undergoing dynamics in unfavorable time regimes cannot be observed. De novo NMR structure determination is seriously hampered when missing signals lead to an incomplete chemical shift assignment resulting in an information content of the NOE data insufficient to determine the structure ab initio. We developed a new protein structure determination strategy for such cases based on a novel NOE assignment strategy utilizing a number of model structures but no explicit reference structure as it is used for bootstrapping like algorithms. The software distinguishes in detail between consistent and mutually exclusive pairs of possible NOE assignments on the basis of different precision levels of measured chemical shifts searching for a set of maximum number of consistent NOE assignments in agreement with 3D space. Validation of the method using the structure of the low molecular-weight-protein tyrosine phosphatase A (MptpA) showed robust results utilizing protein structures with 30-45% sequence identity and 70% of the chemical shift assignments. About 60% of the resonance assignments are sufficient to identify those structural models with highest conformational similarity to the real structure. The software was benchmarked by de novo solution structures of fibroblast growth factor 21 (FGF21) and the extracellular fibroblast growth factor receptor domain FGFR4 D2, which both failed in crystallization trials and in classical NMR structure determination.

BID-induced structural changes in BAK promote apoptosis

The BCL-2 family protein BAK is responsible for mitochondrial outer membrane permeabilization (MOMP), which leads to apoptosis. The BCL-2 homology (BH) 3-only protein BID activates BAK to perform this function. We report the NMR solution structure of the human BID BH3–BAK complex, which identified the activation site at the canonical BH3-binding groove of BAK. Mutating the BAK BH1 in the groove prevented activation and MOMP but not the binding of BID. BAK BH3 mutations allowed BID binding and activation but blunted function by blocking BAK oligomerization. BAK activation follows a “hit-and-run” mechanism whereby BID dissociates from the trigger site allowing BAK oligomerization at an overlapping interface. In contrast, the BH3-only proteins NOXA and BAD are predicted to clash with the trigger site, and are not activators of BAK. These findings provide insights into the early stages of BAK activation.

Structure of Enterococcus faeciuml,d-transpeptidase acylated by ertapenem provides insight into the inactivation mechanism

The maintenance of bacterial cell shape and integrity is largely attributed to peptidoglycan, a biopolymer highly cross-linked through d,d-transpeptidation. Peptidoglycan cross-linking is catalyzed by penicillin-binding proteins (PBPs) that are the essential target of β-lactam antibiotics. PBPs are functionally replaced by l,d-transpeptidases (Ldts) in ampicillin-resistant mutants of Enterococcus faecium and in wild-type Mycobacterium tuberculosis. Ldts are inhibited in vivo by a single class of β-lactams, the carbapenems, which act as a suicide substrate. We present here the first structure of a carbapenem-acylated l,d-transpeptidase, E. faecium Ldtfm acylated by ertapenem, which revealed key contacts between the carbapenem core and residues of the catalytic cavity of the enzyme. Significant reorganization of the antibiotic conformation occurs upon enzyme acylation. These results, together with the analysis of protein-to-carbapenem proton transfers, provide new insights into the mechanism of Ldt acylation by carbapenems.

CING: an integrated residue-based structure validation program suite

We present a suite of programs, named CING for Common Interface for NMR Structure Generation that provides for a residue-based, integrated validation of the structural NMR ensemble in conjunction with the experimental restraints and other input data. External validation programs and new internal validation routines compare the NMR-derived models with empirical data, measured chemical shifts, distance- and dihedral restraints and the results are visualized in a dynamic Web 2.0 report. A red-orange-green score is used for residues and restraints to direct the user to those critiques that warrant further investigation. Overall green scores below ~20 % accompanied by red scores over ~50 % are strongly indicative of poorly modelled structures. The publically accessible, secure iCing webserver ( https://nmr.le.ac.uk ) allows individual users to upload the NMR data and run a CING validation analysis.