Interactions between proteins and metal cations are central to biochemical processes and shape protein structures. SilE, an intrinsically disordered protein involved in bacterial silver resistance, folds into α-helices upon binding Ag+ ions. Focusing on the B1 peptide fragment from SilE, we investigate the mechanism of Ag+-induced folding with atomistic simulations and experiments. Guided by Mass Spectrometry and NMR, we prepare a structural model of Ag+-bound B1, which we parametrize using DFT. Then, with replica-exchange simulations and deep learning, we map B1's folding landscape and how it is shaped by Ag+. Specifically, Ag+ binding promotes folding by lowering the entropy of the disordered state and stabilizing the folded state. We also describe how Ag+ alters the folding pathways. Overall, we improve our understanding of metal-induced protein folding and lay the groundwork for further computational investigations of the bacterial silver-resistance machinery.

The emergence of very high NMR magnetic fields will certainly encourage the study of larger biological systems with their dynamics and interactions. NMR spin relaxation allows probing the dynamical properties of proteins where the 15N longitudinal (R1) and transverse (R2) relaxation rates in addition to the 1H-15N heteronuclear NOE describe the ps-ns time scale. Their analytical representation involves the chemical shift anisotropy (CSA) effect that represents the major contribution at a very high magnetic field above 18.8 T. An accurate analysis of the latter parameters in terms of model free (MF) requires considering its effect. Until now, a uniform value of -160 ppm for the CSA has been widely used to derive the backbone order parameters (S2), giving rise to a large fluctuation of its value at very high magnetic fields. Conversely, the use of a site-specific CSA improves the accurate analysis of protein dynamics but requires a cost-effective experimental multi-field approach. In the present paper, we show how the CSA mainly contributes to the relaxation parameters at 28.2 T compared to lower magnetic fields and may bias the determination of S2. We propose to replace the time-consuming measurement of spin relaxation at multiple fields by a combination of molecular dynamics (MD) and the measurement of spin relaxation at one very high magnetic field only. We applied this strategy to three well-folded proteins (ubiquitin, GB3 and ribonuclease H) to show that the determined order parameters are in good agreement with the ones obtained by means of experimental data only.

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.

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.

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.

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.

AbstractGenome sequencing of the human parasiteSchistosoma mansonirevealed an interesting gene superfamily calledmicro-exon gene(MEG) that encodes MEG secreted proteins. The genes are composed of short exons (3-81 base pairs) with symmetrically inserted long introns (up to 5 kbp). This article recollects 35S. mansonispecificmeggenes 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, propose a phylogenetic analysis and highlight intriguing conserved patterns/motifs in the sequences of these MEG proteins. Based on the analyses, we were able to classify the 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 reproducibility of research results and to avoid possible misunderstandings.Author AbstractSchistosoma mansoniis a parasitic worm, the etiological agent of schistosomiasis or bilharzia, a chronic tropical disease. It is a vector-borne parasite with a complex life cycle and an equally complex genome, assembled in 7 autosomes and a pair of sexual chromosomes. Within the gene products, one superfamily is particularly interesting, since it is specific toSchistosomatidae, highly variable and redundant: the micro-exon gene (MEG) family. As the name implies, these genes are made by short coding exons (3 to 81 base pairs), symmetrically interspersed by long introns (from 0.2 to 5 kbp). There are 35megsallover the chromosomes, which code for at least 87 MEG proteins. We have aligned all of them, constructed a phylogenetic tree and proposed a theory for their duplication and genome colonization. Based on that, we propose a rational nomenclature to help the community to study MEG’s elusive role. We also propose to help WormBaseParaSite to adopt this new nomenclature to avoid giving the same acronym to different protein sequences.

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.

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.

Among the various biophysical methods available to investigate protein dynamics, NMR present the ability to scrutinize protein motions on a broad range of time scales. 1H-15N NMR spin relaxation experiments can reveal the extent of protein motions across the ps-ns dynamics probed by the fundamental parameters 15N-R1, 15N-R2 and 1H-15N NOE that can be well sampled by molecular dynamics simulations (MD). An accurate prediction of these parameters is subjected to a proper description of the rotational diffusion and anisotropy. Indeed, a strong rotational anisotropy has a profound effect on the various relaxation parameters and could be mistaken for conformational exchange. Although the principle of NMR spin relaxation predictions from MD is now well established, numerous NMR/MD comparisons have hitherto focused on proteins that show low to moderate anisotropy and make use of a scaling factor to remove artifacts arising from water model-dependence of the rotational diffusion. In the present work, we have used NMR to characterize the rotational diffusion of the α-helical STAM2-UIM domain by measuring the 15N-R1, 15N-R2 and 1H-15N NOE relaxation parameters. We therefore highlight the use of the polarizable AMOEBA FF and show that it improves the prediction of the rotational diffusion in the particular case of strong rotational anisotropy, which in turn enhances the prediction of the 15N-R1, 15N-R2 and 1H-15N NOE relaxation parameters without requirement to a scaling factor. Our findings suggest that the use of polarizable FFs could potentially enrich our understanding of protein dynamics in situations where charges distribution or protein shape is remodeled over time like in the case of multidomain proteins or intrinsically disordered proteins.

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.

Multidomain proteins represent a broad spectrum of the protein landscape and are involved in various interactions. They could be considered as modular building blocks assembled in distinct fashion and connected by linkers of varying lengths and sequences. Due to their intrinsic flexibility, these linkers provide proteins a subtle way to modulate interactions and explore a wide range of conformational space. In the present study, we are seeking to understand the effect of the flexibility and dynamics of the linker involved in the STAM2 UIM-SH3 dual domain protein with respect to molecular recognition. We have engineered several constructs of UIM-SH3 with different length linkers or domain deletion. By means of SAXS and NMR experiments, we have shown that the modification of the linker modifies the flexibility and the dynamics of UIM-SH3. Indeed, the global tumbling of both the UIM and SH3 domain is different but not independent from each other while the length of the linker has an impact on the ps-ns time scale dynamics of the respective domains. Finally, the modification of the flexibility and dynamics of the linker has a drastic effect on the interaction of UIM-SH3 with Lys63-linked diubiquitin with a roughly eight-time weaker dissociation constant.

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.

1H-15N 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 R1, R2, 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 S2 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.