Cryo-EM structure of the octameric pore of Clostridium perfringens ß-toxin.

Clostridium perfringens is one of the most widely distributed and successful pathogens producing an impressive arsenal of toxins. One of the most potent toxins produced is the C. perfringens ß-toxin (CPB). This toxin is the main virulence factor of type C strains. We describe the cryo-electron microscopy (EM) structure of CPB oligomer. We show that CPB forms homo-octameric pores like the hetero-oligomeric pores of the bi-component leukocidins, with important differences in the receptor binding region and the N-terminal latch domain. Intriguingly, the octameric CPB pore complex contains a second 16-stranded ß-barrel protrusion atop of the cap domain that is formed by the N-termini of the eight protomers. We propose that CPB, together with the newly identified Epx toxins, is a member a new subclass of the hemolysin-like family. In addition, we show that the ß-barrel protrusion domain can be modified without affecting the pore-forming ability, thus making the pore particularly attractive for macromolecule sensing and nanotechnology. The cryo-EM structure of the octameric pore of CPB will facilitate future developments in both nanotechnology and basic research.

Biobox: a toolbox for biomolecular modelling.

MOTIVATION: The implementation of biomolecular modelling methods and analyses can be cumbersome, often carried out with in-house software reimplementing common tasks, and requiring the integration of diverse software libraries. RESULTS: We present Biobox, a Python-based toolbox facilitating the implementation of biomolecular modelling methods. AVAILABILITY AND IMPLEMENTATION: Biobox is freely available on https://github.com/degiacom/biobox, along with its API and interactive Jupyter notebook tutorials.

The epilepsy-associated protein TBC1D24 is required for normal development, survival and vesicle trafficking in mammalian neurons.

Mutations in the Tre2/Bub2/Cdc16 (TBC)1 domain family member 24 (TBC1D24) gene are associated with a range of inherited neurological disorders, from drug-refractory lethal epileptic encephalopathy and DOORS syndrome (deafness, onychodystrophy, osteodystrophy, mental retardation, seizures) to non-syndromic hearing loss. TBC1D24 has been implicated in neuronal transmission and maturation, although the molecular function of the gene and the cause of the apparently complex disease spectrum remain unclear. Importantly, heterozygous TBC1D24 mutation carriers have also been reported with seizures, suggesting that haploinsufficiency for TBC1D24 is significant clinically. Here we have systematically investigated an allelic series of disease-associated mutations in neurons alongside a new mouse model to investigate the consequences of TBC1D24 haploinsufficiency to mammalian neurodevelopment and synaptic physiology. The cellular studies reveal that disease-causing mutations that disrupt either of the conserved protein domains in TBC1D24 are implicated in neuronal development and survival and are likely acting as loss-of-function alleles. We then further investigated TBC1D24 haploinsufficiency in vivo and demonstrate that TBC1D24 is also crucial for normal presynaptic function: genetic disruption of Tbc1d24 expression in the mouse leads to an impairment of endocytosis and an enlarged endosomal compartment in neurons with a decrease in spontaneous neurotransmission. These data reveal the essential role for TBC1D24 at the mammalian synapse and help to define common synaptic mechanisms that could underlie the varied effects of TBC1D24 mutations in neurological disease.

Transmembrane Protein Docking with JabberDock.

Transmembrane proteins act as an intermediary for a broad range of biological process. Making up 20% to 30% of the proteome, their ubiquitous nature has resulted in them comprising 50% of all targets in drug design. Despite their importance, they make up only 4% of all structures in the PDB database, primarily owing to difficulties associated with isolating and characterizing them. Membrane protein docking algorithms could help to fill this knowledge gap, yet only few exist. Moreover, these existing methods achieve success rates lower than the current best soluble proteins docking software. We present and test a pipeline using our software, JabberDock, to dock membrane proteins. JabberDock docks shapes representative of membrane protein structure and dynamics in their biphasic environment. We verify JabberDock's ability to yield accurate predictions by applying it to a benchmark of 20 transmembrane dimers, returning a success rate of 75.0%. This makes our software very competitive among available membrane protein-protein docking tools.

Lipid binding attenuates channel closure of the outer membrane protein OmpF.

Strong interactions between lipids and proteins occur primarily through association of charged headgroups and amino acid side chains, rendering the protonation status of both partners important. Here we use native mass spectrometry to explore lipid binding as a function of charge of the outer membrane porin F (OmpF). We find that binding of anionic phosphatidylglycerol (POPG) or zwitterionic phosphatidylcholine (POPC) to OmpF is sensitive to electrospray polarity while the effects of charge are less pronounced for other proteins in outer or mitochondrial membranes: the ferripyoverdine receptor (FpvA) or the voltage-dependent anion channel (VDAC). Only marginal charge-induced differences were observed for inner membrane proteins: the ammonia channel (AmtB) or the mechanosensitive channel. To understand these different sensitivities, we performed an extensive bioinformatics analysis of membrane protein structures and found that OmpF, and to a lesser extent FpvA and VDAC, have atypically high local densities of basic and acidic residues in their lipid headgroup-binding regions. Coarse-grained molecular dynamics simulations, in mixed lipid bilayers, further implicate changes in charge by demonstrating preferential binding of anionic POPG over zwitterionic POPC to protonated OmpF, an effect not observed to the same extent for AmtB. Moreover, electrophysiology and mass-spectrometry-based ligand-binding experiments, at low pH, show that POPG can maintain OmpF channels in open conformations for extended time periods. Since the outer membrane is composed almost entirely of anionic lipopolysaccharide, with similar headgroup properties to POPG, such anionic lipid binding could prevent closure of OmpF channels, thereby increasing access of antibiotics that use porin-mediated pathways.

Structural and functional consequences of age-related isomerization in α-crystallins.

Long-lived proteins are subject to spontaneous degradation and may accumulate a range of modifications over time, including subtle alterations such as side-chain isomerization. Recently, tandem MS has enabled identification and characterization of such peptide isomers, including those differing only in chirality. However, the structural and functional consequences of these perturbations remain largely unexplored. Here, we examined the impact of isomerization of aspartic acid or epimerization of serine at four sites mapping to crucial oligomeric interfaces in human αA- and αB-crystallin, the most abundant chaperone proteins in the eye lens. To characterize the effect of isomerization on quaternary assembly, we utilized synthetic peptide mimics, enzyme assays, molecular dynamics calculations, and native MS experiments. The oligomerization of recombinant forms of αA- and αB-crystallin that mimic isomerized residues deviated from native behavior in all cases. Isomerization also perturbs recognition of peptide substrates, either enhancing or inhibiting kinase activity. Specifically, epimerization of serine (αASer-162) dramatically weakened inter-subunit binding. Furthermore, phosphorylation of αBSer-59, known to play an important regulatory role in oligomerization, was severely inhibited by serine epimerization and altered by isomerization of nearby αBAsp-62. Similarly, isomerization of αBAsp-109 disrupted a vital salt bridge with αBArg-120, a contact that when broken has previously been shown to yield aberrant oligomerization and aggregation in several disease-associated variants. Our results illustrate how isomerization of amino acid residues, which may seem to be only a minor structural perturbation, can disrupt native structural interactions with profound consequences for protein assembly and activity.

Predicting the Shapes of Protein Complexes through Collision Cross Section Measurements and Database Searches.

In structural biology, collision cross sections (CCSs) from ion mobility mass spectrometry (IM-MS) measurements are routinely compared to computationally or experimentally derived protein structures. Here, we investigate whether CCS data can inform about the shape of a protein in the absence of specific reference structures. Analysis of the proteins in the CCS database shows that protein complexes with low apparent densities are structurally more diverse than those with a high apparent density. Although assigning protein shapes purely on CCS data is not possible, we find that we can distinguish oblate- and prolate-shaped protein complexes by using the CCS, molecular weight, and oligomeric states to mine the Protein Data Bank (PDB) for potentially similar protein structures. Furthermore, comparing the CCS of a ferritin cage to the solution structures in the PDB reveals significant deviations caused by structural collapse in the gas phase. We then apply the strategy to an integral membrane protein by comparing the shapes of a prokaryotic and a eukaryotic sodium/proton antiporter homologue. We conclude that mining the PDB with IM-MS data is a time-effective way to derive low-resolution structural models.

Computational Strategies and Challenges for Using Native Ion Mobility Mass Spectrometry in Biophysics and Structural Biology.

Native mass spectrometry (MS) allows the interrogation of structural aspects of macromolecules in the gas phase, under the premise of having initially maintained their solution-phase noncovalent interactions intact. In the more than 25 years since the first reports, the utility of native MS has become well established in the structural biology community. The experimental and technological advances during this time have been rapid, resulting in dramatic increases in sensitivity, mass range, resolution, and complexity of possible experiments. As experimental methods have improved, there have been accompanying developments in computational approaches for analyzing and exploiting the profusion of MS data in a structural and biophysical context. In this perspective, we consider the computational strategies currently being employed by the community, aspects of best practice, and the challenges that remain to be addressed. Our perspective is based on discussions within the European Cooperation in Science and Technology Action on Native Mass Spectrometry and Related Methods for Structural Biology (EU COST Action BM1403), which involved participants from across Europe and North America. It is intended not as an in-depth review but instead to provide an accessible introduction to and overview of the topic-to inform newcomers to the field and stimulate discussions in the community about addressing existing challenges. Our complementary perspective (http://dx.doi.org/10.1021/acs.analchem.9b05792) focuses on software tools available to help researchers tackle some of the challenges enumerated here.

Software Requirements for the Analysis and Interpretation of Native Ion Mobility Mass Spectrometry Data.

The past few years have seen a dramatic increase in applications of native mass and ion mobility spectrometry, especially for the study of proteins and protein complexes. This increase has been catalyzed by the availability of commercial instrumentation capable of carrying out such analyses. As in most fields, however, the software to process the data generated from new instrumentation lags behind. Recently, a number of research groups have started addressing this by developing software, but further improvements are still required in order to realize the full potential of the data sets generated. In this perspective, we describe practical aspects as well as challenges in processing native mass spectrometry (MS) and ion mobility-MS data sets and provide a brief overview of currently available tools. We then set out our vision of future developments that would bring the community together and lead to the development of a common platform to expedite future computational developments, provide standardized processing approaches, and serve as a location for the deposition of data for this emerging field. This perspective has been written by members of the European Cooperation in Science and Technology Action on Native MS and Related Methods for Structural Biology (EU COST Action BM1403) as an introduction to the software tools available in this area. It is intended to serve as an overview for newcomers and to stimulate discussions in the community on further developments in this field, rather than being an in-depth review. Our complementary perspective (http://dx.doi.org/10.1021/acs.analchem.9b05791) focuses on computational approaches used in this field.

Membrane proteins bind lipids selectively to modulate their structure and function.

Previous studies have established that the folding, structure and function of membrane proteins are influenced by their lipid environments and that lipids can bind to specific sites, for example, in potassium channels. Fundamental questions remain however regarding the extent of membrane protein selectivity towards lipids. Here we report a mass spectrometry approach designed to determine the selectivity of lipid binding to membrane protein complexes. We investigate the mechanosensitive channel of large conductance (MscL) from Mycobacterium tuberculosis and aquaporin Z (AqpZ) and the ammonia channel (AmtB) from Escherichia coli, using ion mobility mass spectrometry (IM-MS), which reports gas-phase collision cross-sections. We demonstrate that folded conformations of membrane protein complexes can exist in the gas phase. By resolving lipid-bound states, we then rank bound lipids on the basis of their ability to resist gas phase unfolding and thereby stabilize membrane protein structure. Lipids bind non-selectively and with high avidity to MscL, all imparting comparable stability; however, the highest-ranking lipid is phosphatidylinositol phosphate, in line with its proposed functional role in mechanosensation. AqpZ is also stabilized by many lipids, with cardiolipin imparting the most significant resistance to unfolding. Subsequently, through functional assays we show that cardiolipin modulates AqpZ function. Similar experiments identify AmtB as being highly selective for phosphatidylglycerol, prompting us to obtain an X-ray structure in this lipid membrane-like environment. The 2.3 Å resolution structure, when compared with others obtained without lipid bound, reveals distinct conformational changes that re-position AmtB residues to interact with the lipid bilayer. Our results demonstrate that resistance to unfolding correlates with specific lipid-binding events, enabling a distinction to be made between lipids that merely bind from those that modulate membrane protein structure and/or function. We anticipate that these findings will be important not only for defining the selectivity of membrane proteins towards lipids, but also for understanding the role of lipids in modulating protein function or drug binding.

It takes a dimer to tango: Oligomeric small heat shock proteins dissociate to capture substrate.

Small heat-shock proteins (sHsps) are ubiquitous molecular chaperones, and sHsp mutations or altered expression are linked to multiple human disease states. sHsp monomers assemble into large oligomers with dimeric substructure, and the dynamics of sHsp oligomers has led to major questions about the form that captures substrate, a critical aspect of their mechanism of action. We show here that substructural dimers of two plant dodecameric sHsps, Ta16.9 and homologous Ps18.1, are functional units in the initial encounter with unfolding substrate. We introduced inter-polypeptide disulfide bonds at the two dodecameric interfaces, dimeric and nondimeric, to restrict how their assemblies can dissociate. When disulfide-bonded at the nondimeric interface, mutants of Ta16.9 and Ps18.1 (TaCT-ACD and PsCT-ACD) were inactive but, when reduced, had WT-like chaperone activity, demonstrating that dissociation at nondimeric interfaces is essential for sHsp activity. Moreover, the size of the TaCT-ACD and PsCT-ACD covalent unit defined a new tetrahedral geometry for these sHsps, different from that observed in the Ta16.9 X-ray structure. Importantly, oxidized Tadimer (disulfide bonded at the dimeric interface) exhibited greatly enhanced ability to protect substrate, indicating that strengthening the dimeric interface increases chaperone efficiency. Temperature-induced size and secondary structure changes revealed that folded sHsp dimers interact with substrate and that dimer stability affects chaperone efficiency. These results yield a model in which sHsp dimers capture substrate before assembly into larger, heterogeneous sHsp-substrate complexes for substrate refolding or degradation, and suggest that tuning the strength of the dimer interface can be used to engineer sHsp chaperone efficiency.

Dimer interface of bovine cytochrome c oxidase is influenced by local posttranslational modifications and lipid binding.

Bovine cytochrome c oxidase is an integral membrane protein complex comprising 13 protein subunits and associated lipids. Dimerization of the complex has been proposed; however, definitive evidence for the dimer is lacking. We used advanced mass spectrometry methods to investigate the oligomeric state of cytochrome c oxidase and the potential role of lipids and posttranslational modifications in its subunit interfaces. Mass spectrometry of the intact protein complex revealed that both the monomer and the dimer are stabilized by large lipid entities. We identified these lipid species from the purified protein complex, thus implying that they interact specifically with the enzyme. We further identified phosphorylation and acetylation sites of cytochrome c oxidase, located in the peripheral subunits and in the dimer interface, respectively. Comparing our phosphorylation and acetylation sites with those found in previous studies of bovine, mouse, rat, and human cytochrome c oxidase, we found that whereas some acetylation sites within the dimer interface are conserved, suggesting a role for regulation and stabilization of the dimer, phosphorylation sites were less conserved and more transient. Our results therefore provide insights into the locations and interactions of lipids with acetylated residues within the dimer interface of this enzyme, and thereby contribute to a better understanding of its structure in the natural membrane. Moreover dimeric cytochrome c oxidase, comprising 20 transmembrane, six extramembrane subunits, and associated lipids, represents the largest integral membrane protein complex that has been transferred via electrospray intact into the gas phase of a mass spectrometer, representing a significant technological advance.

Automatic Assignment of Methyl-NMR Spectra of Supramolecular Machines Using Graph Theory.

Methyl groups are powerful probes for the analysis of structure, dynamics and function of supramolecular assemblies, using both solution- and solid-state NMR. Widespread application of the methodology has been limited due to the challenges associated with assigning spectral resonances to specific locations within a biomolecule. Here, we present Methyl Assignment by Graph Matching (MAGMA), for the automatic assignment of methyl resonances. A graph matching protocol examines all possibilities for each resonance in order to determine an exact assignment that includes a complete description of any ambiguity. MAGMA gives 100% accuracy in confident assignments when tested against both synthetic data, and 9 cross-validated examples using both solution- and solid-state NMR data. We show that this remarkable accuracy enables a user to distinguish between alternative protein structures. In a drug discovery application on HSP90, we show the method can rapidly and efficiently distinguish between possible ligand binding modes. By providing an exact and robust solution to methyl resonance assignment, MAGMA can facilitate significantly accelerated studies of supramolecular machines using methyl-based NMR spectroscopy.

Characterisation of Shigella Spa33 and Thermotoga FliM/N reveals a new model for C-ring assembly in T3SS.

Flagellar type III secretion systems (T3SS) contain an essential cytoplasmic-ring (C-ring) largely composed of two proteins FliM and FliN, whereas an analogous substructure for the closely related non-flagellar (NF) T3SS has not been observed in situ. We show that the spa33 gene encoding the putative NF-T3SS C-ring component in Shigella flexneri is alternatively translated to produce both full-length (Spa33-FL) and a short variant (Spa33-C), with both required for secretion. They associate in a 1:2 complex (Spa33-FL/C2) that further oligomerises into elongated arrays in vitro. The structure of Spa33-C2 and identification of an unexpected intramolecular pseudodimer in Spa33-FL reveal a molecular model for their higher order assembly within NF-T3SS. Spa33-FL and Spa33-C are identified as functional counterparts of a FliM-FliN fusion and free FliN respectively. Furthermore, we show that Thermotoga maritima FliM and FliN form a 1:3 complex structurally equivalent to Spa33-FL/C2 , allowing us to propose a unified model for C-ring assembly by NF-T3SS and flagellar-T3SS.

EM∩IM: software for relating ion mobility mass spectrometry and electron microscopy data.

We present EM∩IM, software that allows the calculation of collision cross-sections from electron density maps obtained for example by means of transmission electron microscopy. This allows the assessment of structures other than those described by atomic coordinates with ion mobility mass spectrometry data, and provides a new means for contouring and validating electron density maps. EM∩IM thereby facilitates the use of data obtained in the gas phase within structural biology studies employing diverse experimental methodologies.

Molecular assembly of the aerolysin pore reveals a swirling membrane-insertion mechanism.

Aerolysin is the founding member of a superfamily of ß-pore-forming toxins whose pore structure is unknown. We have combined X-ray crystallography, cryo-EM, molecular dynamics and computational modeling to determine the structures of aerolysin mutants in their monomeric and heptameric forms, trapped at various stages of the pore formation process. A dynamic modeling approach based on swarm intelligence was applied, whereby the intrinsic flexibility of aerolysin extracted from new X-ray structures was used to fully exploit the cryo-EM spatial restraints. Using this integrated strategy, we obtained a radically new arrangement of the prepore conformation and a near-atomistic structure of the aerolysin pore, which is fully consistent with all of the biochemical data available so far. Upon transition from the prepore to pore, the aerolysin heptamer shows a unique concerted swirling movement, accompanied by a vertical collapse of the complex, ultimately leading to the insertion of a transmembrane ß-barrel.

Dual chaperone role of the C-terminal propeptide in folding and oligomerization of the pore-forming toxin aerolysin.

Throughout evolution, one of the most ancient forms of aggression between cells or organisms has been the production of proteins or peptides affecting the permeability of the target cell membrane. This class of virulence factors includes the largest family of bacterial toxins, the pore-forming toxins (PFTs). PFTs are bistable structures that can exist in a soluble and a transmembrane state. It is unclear what drives biosynthetic folding towards the soluble state, a requirement that is essential to protect the PFT-producing cell. Here we have investigated the folding of aerolysin, produced by the human pathogen Aeromonas hydrophila, and more specifically the role of the C-terminal propeptide (CTP). By combining the predictive power of computational techniques with experimental validation using both structural and functional approaches, we show that the CTP prevents aggregation during biosynthetic folding. We identified specific residues that mediate binding of the CTP to the toxin. We show that the CTP is crucial for the control of the aerolysin activity, since it protects individual subunits from aggregation within the bacterium and later controls assembly of the quaternary pore-forming complex at the surface of the target host cell. The CTP is the first example of a C-terminal chain-linked chaperone with dual function.

Identification of Graphene Dispersion Agents through Molecular Fingerprints.

The scalable production and dispersion of 2D materials, like graphene, is critical to enable their use in commercial applications. While liquid exfoliation is commonly used, solvents such as N-methyl-pyrrolidone (NMP) are toxic and difficult to scale up. However, the search for alternative solvents is hindered by the intimidating size of the chemical space. Here, we present a computational pipeline informing the identification of effective exfoliation agents. Classical molecular dynamics simulations provide statistical sampling of interactions, enabling the identification of key molecular descriptors for a successful solvent. The statistically representative configurations from these simulations, studied with quantum mechanical calculations, allow us to gain insights onto the chemophysical interactions at the surface-solvent interface. As an exemplar, through this pipeline we identify a potential graphene exfoliation agent 2-pyrrolidone and experimentally demonstrate it to be as effective as NMP. Our workflow can be generalized to any 2D material and solvent system, enabling the screening of a wide range of compounds and solvents to identify safer and cheaper means of producing dispersions.

Complementing machine learning-based structure predictions with native mass spectrometry.

The advent of machine learning-based structure prediction algorithms such as AlphaFold2 (AF2) and RoseTTa Fold have moved the generation of accurate structural models for the entire cellular protein machinery into the reach of the scientific community. However, structure predictions of protein complexes are based on user-provided input and may require experimental validation. Mass spectrometry (MS) is a versatile, time-effective tool that provides information on post-translational modifications, ligand interactions, conformational changes, and higher-order oligomerization. Using three protein systems, we show that native MS experiments can uncover structural features of ligand interactions, homology models, and point mutations that are undetectable by AF2 alone. We conclude that machine learning can be complemented with MS to yield more accurate structural models on a small and large scale.

Electrospray ionization of native membrane proteins proceeds via a charge equilibration step.

Electrospray ionization mass spectrometry is increasingly applied to study the structures and interactions of membrane protein complexes. However, the charging mechanism is complicated by the presence of detergent micelles during ionization. Here, we show that the final charge of membrane proteins can be predicted by their molecular weight when released from the non-charge reducing saccharide detergents. Our data indicate that PEG detergents lower the charge depending on the number of detergent molecules in the surrounding micelle, whereas fos-choline detergents may additionally participate in ion-ion reactions after desolvation. The supercharging reagent sulfolane, on the other hand, has no discernible effect on the charge of detergent-free membrane proteins. Taking our observations into the context of protein-detergent interactions in the gas phase, we propose a charge equilibration model for the generation of native-like membrane protein ions. During ionization of the protein-detergent complex, the ESI charges are distributed between detergent and protein according to proton affinity of the detergent, number of detergent molecules, and surface area of the protein. Charge equilibration influenced by detergents determines the final charge state of membrane proteins. This process likely contributes to maintaining a native-like fold after detergent release and can be harnessed to stabilize particularly labile membrane protein complexes in the gas phase.