Revisiting Schistosoma mansoni Micro-Exon Gene (MEG) Protein Family: A Tour into Conserved Motifs and Annotation

Genome sequencing of the human parasite Schistosoma mansoni revealed an interesting gene superfamily, called micro-exon gene (meg), that encodes secreted MEG proteins. The genes are composed of short exons (3-81 base pairs) regularly interspersed with long introns (up to 5 kbp). This article recollects 35 S. mansoni specific meg genes that are distributed over 7 autosomes and one pair of sex chromosomes and that code for at least 87 verified MEG proteins. We used various bioinformatics tools to produce an optimal alignment and propose a phylogenetic analysis. This work highlighted intriguing conserved patterns/motifs in the sequences of the highly variable MEG proteins. Based on the analyses, we were able to classify the verified MEG proteins into two subfamilies and to hypothesize their duplication and colonization of all the chromosomes. Together with motif identification, we also proposed to revisit MEGs' common names and annotation in order to avoid duplication, to help the reproducibility of research results and to avoid possible misunderstandings.

Divide, conquer and reconstruct: How to solve the 3D structure of recalcitrant Micro-Exon Gene (MEG) protein from Schistosoma mansoni

Micro-Exon Genes are a widespread class of genes known for their high variability, widespread in the genome of parasitic trematodes such as Schistosoma mansoni. In this study, we present a strategy that allowed us to solve the structures of three alternatively spliced isoforms from the Schistoma mansoni MEG 2.1 family for the first time. All isoforms are hydrophobic, intrinsically disordered, and recalcitrant to be expressed in high yield in heterologous hosts. We resorted to the chemical synthesis of shorter pieces, before reconstructing the entire sequence. Here, we show that isoform 1 partially folds in a-helix in the presence of trifluoroethanol while isoform 2 features two rigid elbows, that maintain the peptide as disordered, preventing any structuring. Finally, isoform 3 is dominated by the signal peptide, which folds into a-helix. We demonstrated that combining biophysical techniques, like circular dichroism and nuclear magnetic resonance at natural abundance, with in silico molecular dynamics simulation for isoform 1 only, was the key to solve the structure of MEG 2.1. Our results provide a crucial piece to the puzzle of this elusive and highly variable class of proteins.

The battle for silver binding: How the interplay between the SilE, SilF, and SilB proteins contributes to the silver efflux pump mechanism

The resistance of gram-negative bacteria to silver ions is mediated by a silver efflux pump, which mainly relies on a tripartite efflux complex SilCBA, a metallochaperone SilF and an intrinsically disordered protein SilE. However, the precise mechanism by which silver ions are extruded from the cell and the different roles of SilB, SilF, and SilE remain poorly understood. To address these questions, we employed nuclear magnetic resonance and mass spectrometry to investigate the interplay between these proteins. We first solved the solution structures of SilF in its free and Ag+-bound forms, and we demonstrated that SilB exhibits two silver binding sites in its N and C termini. Conversely to the homologous Cus system, we determined that SilF and SilB interact without the presence of silver ions and that the rate of silver dissociation is eight times faster when SilF is bound to SilB, indicating the formation of a SilF-Ag-SilB intermediate complex. Finally, we have shown that SilE does not bind to either SilF or SilB, regardless of the presence or absence of silver ions, further corroborating that it merely acts as a regulator that prevents the cell from being overloaded with silver. Collectively, we have provided further insights into protein interactions within the sil system that contribute to bacterial resistance to silver ions.

Structural and dynamical insights into SilE silver binding from combined analytical probes

Silver has been used for its antimicrobial properties to fight infection for thousands of years. Unfortunately, some Gram-negative bacteria have developed silver resistance causing the death of patients in a burn unit. The genes responsible for silver resistance have been designated as the sil operon. Among the proteins of the sil operon, SilE has been shown to play a key role in bacterial silver resistance. Based on the limited information available, it has been depicted as an intrinsically disordered protein that folds into helices upon silver ion binding. Herein, this work demonstrates that SilE is composed of 4 clearly identified helical segments in the presence of several silver ions. The combination of analytical and biophysical techniques (NMR spectroscopy, CD, SAXS, HRMS, CE-ICP-MS, and IM-MS) reveals that SilE harbors four strong silver binding sites among the eight sites available. We have also further evidenced that SilE does not adopt a globular structure but rather samples a large conformational space from elongated to more compact structures. This particular structural organization facilitates silver binding through much higher accessibility of the involved His and Met residues. These valuable results will advance our current understanding of the role of SilE in the silver efflux pump complex mechanism and will help in the future rational design of inhibitors to fight bacterial silver resistance.

NLRP3 phosphorylation in its LRR domain critically regulates inflammasome assembly

NLRP3 controls the secretion of inflammatory cytokines IL-1β/18 and pyroptosis by assembling the inflammasome. Upon coordinated priming and activation stimuli, NLRP3 recruits NEK7 within hetero-oligomers that nucleate ASC and caspase-1 filaments, but the apical molecular mechanisms underlying inflammasome assembly remain elusive. Here we show that NEK7 recruitment to NLRP3 is controlled by the phosphorylation status of NLRP3 S803 located within the interaction surface, in which NLRP3 S803 is phosphorylated upon priming and later dephosphorylated upon activation. Phosphomimetic substitutions of S803 abolish NEK7 recruitment and inflammasome activity in macrophages in vitro and in vivo. In addition, NLRP3-NEK7 binding is also essential for NLRP3 deubiquitination by BRCC3 and subsequently inflammasome assembly, with NLRP3 phosphomimetic mutants showing enhanced ubiquitination and degradation than wildtype NLRP3. Finally, we identify CSNK1A1 as the kinase targeting NLRP3 S803. Our findings thus reveal NLRP3 S803 phosphorylation status as a druggable apical molecular mechanism controlling inflammasome assembly.

Nlrp3 inflammasome activation requires Nek7 recruitment to drive ASC speck formation. Here the authors show how Nlrp3 phosphorylation events control this Nek7 recruitment.

NMR reveals the interplay between SilE and SilB model peptides in the context of silver resistance

SilE and SilB are both proteins involved in the silver efflux pump found in Gram-negative bacteria such as S. typhimurium. Using model peptides along with NMR and CD experiments, we show how SilE may store silver ions prior to delivery and we hypothesize for the first time the interplay between SilB and SilE.

Illustrative sketch representative of the interplay between the SilB and SilE proteins in the context of silver bacteria resistance. In case of an excess of silver, SilE plays the role of a regulator and avoids a collapse of the silver efflux pump.

Regulation of measles virus gene expression by P protein coiled-coil properties

The polymerase of negative-stranded RNA viruses consists of the large protein (L) and the phosphoprotein (P), the latter serving both as a chaperon and a cofactor for L. We mapped within measles virus (MeV) P the regions responsible for binding and stabilizing L and showed that the coiled-coil multimerization domain (MD) of P is required for gene expression. MeV MD is kinked as a result of the presence of a stammer. Both restoration of the heptad regularity and displacement of the stammer strongly decrease or abrogate activity in a minigenome assay. By contrast, P activity is rather tolerant of substitutions within the stammer. Single substitutions at the "a" or "d" hydrophobic anchor positions with residues of variable hydrophobicity revealed that P functionality requires a narrow range of cohesiveness of its MD. Results collectively indicate that, beyond merely ensuring P oligomerization, the MD finely tunes viral gene expression through its cohesiveness.

Alpha-helical folding of SilE models upon Ag(His)(Met) motif formation

The SilE protein is suspected to have a prominent role in Ag+ detoxification of silver resistant bacteria. Using model peptides, we elucidated both qualitative and quantitative aspects of the Ag+-induced α-helical structuring role of His- and Met-rich sequences of SilE, improving our understanding of its function within the Sil system.

Ab Initio Prediction of NMR Spin Relaxation Parameters from Molecular Dynamics Simulations

¹H-¹⁵N NMR spin relaxation and relaxation dispersion experiments can reveal the time scale and extent of protein motions across the ps-ms range, where the ps-ns dynamics revealed by fundamental quantities R₁, R₂, and heteronuclear NOE can be well-sampled by molecular dynamics simulations (MD). Although the principles of relaxation prediction from simulations are well-established, numerous NMR-MD comparisons have hitherto focused upon the aspect of order parameters S² due to common artifacts in the prediction of transient dynamics. We therefore summarize here all necessary components and highlight existing and proposed solutions, such as the inclusion of quantum mechanical zero-point vibrational corrections and separate MD convergence of global and local motions in coarse-grained and all-atom force fields, respectively. For the accuracy of the MD prediction to be tested, two model proteins GB3 and Ubiquitin are used to validate five atomistic force fields against published NMR data supplemented by the coarse-grained force field MARTINI+EN. In Amber and CHARMM-type force fields, quantitative agreement was achieved for structured elements with minimum adjustment of global parameters. Deviations from experiment occur in flexible loops and termini, indicating differences in both the extent and time scale of backbone motions. The lack of systematic patterns and water model dependence suggests that modeling of the local environment limits prediction accuracy. Nevertheless, qualitative accuracy in a 2 μs CHARMM36m Stam2 VHS domain simulation demonstrates the potential of MD-based interpretation in combination with NMR-measured dynamics, increasing the utility of spin relaxation in integrative structural biology.

New model of metalloantibiotic: synthesis, structure and biological activity of a zinc(ii) mononuclear complex carrying two enrofloxacin and sulfadiazine antibiotics

Synthesis of a new model of the Zn-based complex with two complementary antibiotics (sulfonamide and quinolone).

Model peptide studies of Ag+ binding sites from the silver resistance protein SilE

A model peptide study characterizes several Ag⁺-binding sites of the bacterial silver resistant protein SilE, providing new insights into its physiological role.

NMR Reveals the Interplay among the AMSH SH3 Binding Motif, STAM2, and Lys63-Linked Diubiquitin

AMSH [associated molecule with a Src homology 3 domain of signal transducing adaptor molecule (STAM)] is one of the deubiquitinating enzymes associated in the regulation of endocytic cargo trafficking. It shows an exquisite selectivity for Lys63-linked polyubiquitin chains that are the main chains involved in cargo sorting. The first step requires the ESCRT-0 complex that comprises the STAM and hepatocyte growth factor-regulated substrate (Hrs) proteins. Previous studies have shown that the presence of the STAM protein increases the efficiency of Lys63-linked polyubiquitin chain cleavage by AMSH, one of the deubiquitinating enzyme involved in lysosomal degradation. In the present study, we are seeking to understand if a particular structural organization among these three key players is responsible for the stimulation of the catalytic activity of AMSH. To address this question, we first monitored the interaction between the ubiquitin interacting motif (UIM)-SH3 construct of STAM2 and the Lys63-linked diubiquitin (Lys63-Ub₂) chains by means of NMR. We show that Lys63-Ub₂ is able to bind either the UIM or the SH3 domain without any selectivity. We further demonstrate that the SH3 binding motif (SBM) of AMSH (AMSH-SBM) outcompetes Lys63-Ub₂ for binding SH3. Additionally, we show how different AMSH-SBM variants, modified by their sequence and length, exhibit similar equilibrium dissociation constants when binding SH3 but significantly differ in their dissociation rate constants. Finally, we report the solution NMR structure of the AMSH-SBM/SH3 complex and propose a structural organization where the AMSH-SBM interacts with the STAM2-SH3 domain and contributes to the correct positioning of AMSH prior to polyubiquitin chains' cleavage.

Francisella tularensis IglG Belongs to a Novel Family of PAAR-Like T6SS Proteins and Harbors a Unique N-terminal Extension Required for Virulence

The virulence of Francisella tularensis, the etiological agent of tularemia, relies on an atypical type VI secretion system (T6SS) encoded by a genomic island termed the Francisella Pathogenicity Island (FPI). While the importance of the FPI in F. tularensis virulence is clearly established, the precise role of most of the FPI-encoded proteins remains to be deciphered. In this study, using highly virulent F. tularensis strains and the closely related species F. novicida, IglG was characterized as a protein featuring a unique α-helical N-terminal extension and a domain of unknown function (DUF4280), present in more than 250 bacterial species. Three dimensional modeling of IglG and of the DUF4280 consensus protein sequence indicates that these proteins adopt a PAAR-like fold, suggesting they could cap the T6SS in a similar way as the recently described PAAR proteins. The newly identified PAAR-like motif is characterized by four conserved cysteine residues, also present in IglG, which may bind a metal atom. We demonstrate that IglG binds metal ions and that each individual cysteine is required for T6SS-dependent secretion of IglG and of the Hcp homologue, IglC and for the F. novicida intracellular life cycle. In contrast, the Francisella-specific N-terminal α-helical extension is not required for IglG secretion, but is critical for F. novicida virulence and for the interaction of IglG with another FPI-encoded protein, IglF. Altogether, our data suggest that IglG is a PAAR-like protein acting as a bi-modal protein that may connect the tip of the Francisella T6SS with a putative T6SS effector, IglF.

Francisella tularensis is a highly pathogenic bacterium causing tularemia. Its ability to cause disease is linked to its ability to replicate in the macrophage cytosol. The intracellular life cycle of Francisella is controlled by a type VI secretion system (T6SS), which is thought to inject effectors into the host cell to allow bacterial escape into the host cytosol. The molecular mechanisms behind this process are still largely unclear. In this work, we identify IglG as a protein with two important domains, one conserved in proteins from more than 250 bacterial species (DUF4280, renamed here as PAAR-like domain) and one specific for the Francisella genus. Using protein sequence analysis and three-dimensional structure predictions, comparative modeling and biochemistry approaches, our data demonstrate that IglG is a metal-binding protein that based on its PAAR-like domain might cap the VgrG spike of the T6SS and act as a membrane-puncturing protein. Furthermore, we identified that the Francisella-specific domain is directly involved in forming a protein complex with another virulence protein, IglF. This work, in addition to enhancing the molecular understanding of the Francisella T6SS, defines the features of the conserved DUF4280, a novel PAAR-like domain involved in type VI secretion (T6S) of many bacterial species.

Competitive binding of UBPY and ubiquitin to the STAM2 SH3 domain revealed by NMR

To date, the signal transducing adaptor molecule 2 (STAM2) was shown to harbour two ubiquitin binding domains (UBDs) known as the VHS and UIM domains, while the SH3 domain of STAM2 was reported to interact with deubiquitinating enzymes (DUBs) like UBPY and AMSH. In the present study, NMR evidences the interaction of the STAM2 SH3 domain with ubiquitin, demonstrating that SH3 constitutes the third UBD of STAM2. Furthermore, we show that a UBPY-derived peptide can outcompete ubiquitin for SH3 binding and vice versa. These results suggest that the SH3 domain of STAM2 plays versatile roles in the context of ubiquitin mediated receptor sorting.

Dynamical properties of the loop 320s of substrate-free and substrate-bound muscle creatine kinase by NMR: evidence for independent subunits

Muscle creatine kinase (MCK; EC2.7.3.2) is a 86 kDa homodimer that belongs to the family of guanidino kinases. MCK has been intensively studied for several decades, but it is still not known why it is a dimer because this quaternary structure does not translate into obvious structural or functional advantages over the homologous monomeric arginine kinase. In particular, it remains to be demonstrated whether MCK subunits are independent. Here, we describe NMR chemical-shift perturbation and relaxation experiments designed to study the active site 320s flexible loop of this enzyme. The analysis was performed with the enzyme in its ligand-free and MgADP-complexed forms, as well as with the transition-state analogue abortive complex (MCK-Mg-ADP-creatine-nitrate ion). Our data indicate that each subunit can bind substrates independently.

Interaction domain on thioredoxin for Pseudomonas aeruginosa 5'-adenylylsulfate reductase

NMR spectroscopy has been used to map the interaction domain on Escherichia coli thioredoxin for the thioredoxin- dependent 5'-adenylylsulfate reductase from Pseudomonas aeruginosa (PaAPR). Seventeen thioredoxin amino acids, all clustered around Cys-32 (the more surface-exposed of the two active-site cysteines), have been located at the PaAPR binding site. The center of the binding domain is dominated by nonpolar amino acids, with a smaller number of charged and polar amino acids located on the periphery of the site. Twelve of the amino acids detected by NMR have non-polar, hydrophobic side chains, including one aromatic amino acid (Trp-31). Four of the thioredoxin amino acids at the PaAPR binding site have polar side chains (Lys-36, Asp-61, Gln-62 and Arg-73), with three of the four having charged side chains. Site-directed mutagenesis experiments have shown that replacement of Lys-36, Asp-61, and Arg-73 and of the absolutely conserved Trp-31 significantly decreases the V(max) for the PaAPR-catalyzed reduction of 5'-adenylylsulfate, with E. coli thioredoxin serving as the electron donor. The most dramatic effect was observed with the W31A variant, which showed no activity as a donor to PaAPR. Although the thiol of the active-site Cys-256 of PaAPR is the point of the initial nucleophilic attack by reduced thioredoxin, mutagenic replacement of Cys-256 by serine has no effect on thioredoxin binding to PaAPR.

Protein Side-Chain Dynamics Observed by Solution- and Solid-State NMR: Comparative Analysis of Methyl 2H Relaxation Data

Rapid advances in solid-state MAS NMR made it possible to probe protein dynamics on a per-residue basis, similar to solution experiments. In this work we compare methyl 2H relaxation rates measured in the solid and liquid samples of alpha-spectrin SH3 domain. The solution data are treated using a model-free approach to separate the contributions from the overall molecular tumbling and fast internal motion. The latter part forms the basis for comparison with the solid-state data. Although the accuracy of solid-state measurements is limited by deuterium spin diffusion, the results suggest a significant similarity between methyl dynamics in the two samples. This is a potentially important observation, preparing the ground for combined analysis of the dynamics data by solid- and solution-state NMR.

Characterization of dynamic processes using deuterium in uniformly 2H,13C,15N enriched peptides by MAS solid-state NMR

We present in this paper 2H,13C MAS correlation experiments that are performed on a uniformly 2H,13C,15N labeled sample of Nac-Val, and on the uniformly 2H,15N labeled dipeptide Nac-Val-Leu-OH. The experiments involve the measurement of 2H T1 relaxation times at two different magnetic fields, as well as the measurement of the 2H tensor parameters by evolution of the 2H chemical shift. The data are interpreted quantitatively to differentiate between different side chain motional models.

13C/15N distance determination by CPMAS NMR in uniformly 13C labeled molecules

The REDOR and CPMAS techniques are applied for measuring 13C-15N dipolar coupling constants in glycine. It is shown that the selective CP or SPECIFIC CP technique removes the coherent evolution of the spin system under homonuclear 13C-13C J couplings. While the large coupling constant (approximately 900 Hz) is readily determined because of the presence of large oscillations in the CPMAS dynamics, their absence precludes the measurement of the small coupling constant (approximately 200 Hz). The experimental results and numerical simulations demonstrate that the determination of 13C-15N coupling constants of medium size (<1 kHz) by the CPMAS technique is mainly limited by the strength of the 1H decoupling field and the size of the 13C and 15N chemical shift anisotropies.

Characterization of Dynamics of Perdeuterated Proteins by MAS Solid-State NMR

We show in this communication that dynamic information for uniformly 2H,13C,15N isotopically enriched, crystalline proteins can be obtained by MAS solid-state NMR spectroscopy. The experiments make use of the deuterium quadrupolar tensor, which is the dominant interaction mechanism. Dynamic properties are accessed by measurement of the size of the quadrupolar coupling constant, Cq, and the value of the asymmetry parameter, eta, via evolution of the deuterium chemical shift, as well as by measurement of deuterium T1 relaxation times. Three-dimensional experiments are performed in order to obtain site-specific resolution. Due to proton dilution, no proton decoupling is required in the carbon evolution periods at MAS rotation frequencies of 10 kHz.

1H/31P distance determination by solid state NMR in multiple-spin systems

The results of two techniques of dipolar recoupling, REDOR and CPMAS, are compared in the case of a coupled multiple-spin system. A fundamentally different behavior is observed for these two techniques. In REDOR, the terms associated with each interaction S-I(k) commute with each other and no truncation takes place so that each addition of spin I(k) causes a splitting with its dipolar frequency. In CPMAS, the flip-flop terms of the dipolar Hamiltonian do not commute with the dominant term from the strongly coupled spin pair so that the weak coupling terms from the neighboring spin I(k) are effectively truncated by the dominant pair interaction. Spin dynamics calculations are in agreement with the experimental data in a cubane shaped cluster.

Molecular dynamics as studied by static-powder and magic-angle spinning 2H NMR

The 2H NMR magic-angle spinning (MAS) technique is compared to the static-powder quadrupole echo (QE) and Jeener-Brockaert (JB) pulse sequences for a quantitative investigation of molecular dynamics in solids. The linewidth of individual spinning sidebands of the one-dimensional MAS spectra are observed to be characteristic of the correlation time from approximately 10(-2) to approximately 10(-8) s so that the dynamic range is increased by approximately three orders of magnitude when compared to the QE experiment. As a consequence, MAS 2H NMR is found to be more sensitive to the presence of an inhomogeneous distribution of correlation times than the QE and JB experiments which rely upon lineshape distortions due to anisotropic T2 and T1Q relaxation, respectively. All these results are demonstrated experimentally and numerically using the two-site flip motion of dimethyl sulfone and of the nitrobenzene guest in the alpha-p-tert-butylcalix[4]arene-nitrobenzene inclusion compound.