Modeling and Structure Determination of Homo-Oligomeric Proteins: An Overview of Challenges and Current Approaches

The oligomeric states of dye-decolorizing peroxidases from Streptomyces lividans and their implications for mechanism of substrate oxidation

A common evolutionary mechanism in biology to drive function is protein oligomerization. In prokaryotes, the symmetrical assembly of repeating protein units to form homomers is widespread, yet consideration in vitro of whether such assemblies have functional or mechanistic consequences is often overlooked. Dye-decolorizing peroxidases (DyPs) are one such example, where their dimeric α + β barrel units can form various oligomeric states, but the oligomer influence, if any, on mechanism and function has received little attention. In this work, we have explored the oligomeric state of three DyPs found in Streptomyces lividans, each with very different mechanistic behaviors in their reactions with hydrogen peroxide and organic substrates. Using analytical ultracentrifugation, we reveal that except for one of the A-type DyPs where only a single sedimenting species is detected, oligomer states ranging from homodimers to dodecamers are prevalent in solution. Using cryo-EM on preparations of the B-type DyP, we determined a 3.02 Å resolution structure of a hexamer assembly that corresponds to the dominant oligomeric state in solution as determined by analytical ultracentrifugation. Furthermore, cryo-EM data detected sub-populations of higher-order oligomers, with one of these formed by an arrangement of two B-type DyP hexamers to give a dodecamer assembly. Our solution and structural insights of these oligomer states provide a new framework to consider previous mechanistic studies of these DyP members and are discussed in terms of long-range electron transfer for substrate oxidation and in the "storage" of oxidizable equivalents on the heme until a two-electron donor is available.

Democratizing protein language models with parameter-efficient fine-tuning

Proteomics has been revolutionized by large protein language models (PLMs), which learn unsupervised representations from large corpora of sequences. These models are typically fine-tuned in a supervised setting to adapt the model to specific downstream tasks. However, the computational and memory footprint of fine-tuning (FT) large PLMs presents a barrier for many research groups with limited computational resources. Natural language processing has seen a similar explosion in the size of models, where these challenges have been addressed by methods for parameter-efficient fine-tuning (PEFT). In this work, we introduce this paradigm to proteomics through leveraging the parameter-efficient method LoRA and training new models for two important tasks: predicting protein-protein interactions (PPIs) and predicting the symmetry of homooligomer quaternary structures. We show that these approaches are competitive with traditional FT while requiring reduced memory and substantially fewer parameters. We additionally show that for the PPI prediction task, training only the classification head also remains competitive with full FT, using five orders of magnitude fewer parameters, and that each of these methods outperform state-of-the-art PPI prediction methods with substantially reduced compute. We further perform a comprehensive evaluation of the hyperparameter space, demonstrate that PEFT of PLMs is robust to variations in these hyperparameters, and elucidate where best practices for PEFT in proteomics differ from those in natural language processing. All our model adaptation and evaluation code is available open-source at https://github.com/microsoft/peft_proteomics. Thus, we provide a blueprint to democratize the power of PLM adaptation to groups with limited computational resources.

Rationally Engineering pH Adaptation of Acid-Induced Arginine Decarboxylase from Escherichia coli to Alkaline Environments to Efficiently Biosynthesize Putrescine

Acid-induced arginine decarboxylase AdiA is a typical homo-oligomeric protein biosynthesizing alkaline nylon monomer putrescine. However, upon loss of the AdiA decamer oligomeric state at neutral and alkaline conditions the activity also diminishes, obstructing the whole-cell biosynthesis of alkaline putrescine. Here, a structure cohesion strategy is proposed to change the pH adaptation of AdiA to alkaline environments based on the rational engineering of meridional and latitudinal oligomerization interfaces. After integrating substitutions of E467K at the latitudinal interface and H736E at the meridional channel interface, the structural stability of AdiA decamer and its substrate transport efficiency at neutral and alkaline conditions are improved. Finally, E467K_H736E is well adapted to neutral and alkaline environments (pH 7.0-9.0), and its enzymatic activity is 35-fold higher than that of wild AdiA at pH 8.0. Using E467K_H736E in the putrescine synthesis pathway, the titer of putrescine is up to 128.9 g·L-1 with a conversion of 0.94 mol·mol-1 in whole-cell catalysis. Additionally, the neutral pH adaptation of lysine decarboxylase, with a decamer structure similar to AdiA, is also improved using this cohesion strategy, providing an option for pH-adaptation engineering of other oligomeric decarboxylases.

Structural and functional analysis of two SHMT8 variants associated with soybean cyst nematode resistance

Two amino acid variants in soybean serine hydroxymethyltransferase 8 (SHMT8) are associated with resistance to the soybean cyst nematode (SCN), a devastating agricultural pathogen with worldwide economic impacts on soybean production. SHMT8 is a cytoplasmic enzyme that catalyzes the pyridoxal 5-phosphate-dependent conversion of serine and tetrahydrofolate (THF) to glycine and 5,10-methylenetetrahydrofolate. A previous study of the P130R/N358Y double variant of SHMT8, identified in the SCN-resistant soybean cultivar (cv.) Forrest, showed profound impairment of folate binding affinity and reduced THF-dependent enzyme activity, relative to the highly active SHMT8 in cv. Essex, which is susceptible to SCN. Given the importance of SCN-resistance in soybean agriculture, we report here the biochemical and structural characterization of the P130R and N358Y single variants to elucidate their individual effects on soybean SHMT8. We find that both single variants have reduced THF-dependent catalytic activity relative to Essex SHMT8 (10- to 50-fold decrease in kcat /Km ) but are significantly more active than the P130R/N368Y double variant. The kinetic data also show that the single variants lack THF-substrate inhibition as found in Essex SHMT8, an observation with implications for regulation of the folate cycle. Five crystal structures of the P130R and N358Y variants in complex with various ligands (resolutions from 1.49 to 2.30 Å) reveal distinct structural impacts of the mutations and provide new insights into allosterism. Our results support the notion that the P130R/N358Y double variant in Forrest SHMT8 produces unique and unexpected effects on the enzyme, which cannot be easily predicted from the behavior of the individual variants.

Co-evolution at protein-protein interfaces guides inference of stoichiometry of oligomeric protein complexes by de novo structure prediction

The quaternary structure with specific stoichiometry is pivotal to the specific function of protein complexes. However, determining the structure of many protein complexes experimentally remains a major bottleneck. Structural bioinformatics approaches, such as the deep learning algorithm Alphafold2-multimer (AF2-multimer), leverage the co-evolution of amino acids and sequence-structure relationships for accurate de novo structure and contact prediction. Pseudo-likelihood maximization direct coupling analysis (plmDCA) has been used to detect co-evolving residue pairs by statistical modeling. Here, we provide evidence that combining both methods can be used for de novo prediction of the quaternary structure and stoichiometry of a protein complex. We achieve this by augmenting the existing AF2-multimer confidence metrics with an interpretable score to identify the complex with an optimal fraction of native contacts of co-evolving residue pairs at intermolecular interfaces. We use this strategy to predict the quaternary structure and non-trivial stoichiometries of Bacillus subtilis spore germination protein complexes with unknown structures. Co-evolution at intermolecular interfaces may therefore synergize with AI-based de novo quaternary structure prediction of structurally uncharacterized bacterial protein complexes.

A morpheein equilibrium regulates catalysis in phosphoserine phosphatase SerB2 from Mycobacterium tuberculosis

Mycobacterium tuberculosis phosphoserine phosphatase MtSerB2 is of interest as a new antituberculosis target due to its essential metabolic role in L-serine biosynthesis and effector functions in infected cells. Previous works indicated that MtSerB2 is regulated through an oligomeric transition induced by L-Ser that could serve as a basis for the design of selective allosteric inhibitors. However, the mechanism underlying this transition remains highly elusive due to the lack of experimental structural data. Here we describe a structural, biophysical, and enzymological characterisation of MtSerB2 oligomerisation in the presence and absence of L-Ser. We show that MtSerB2 coexists in dimeric, trimeric, and tetrameric forms of different activity levels interconverting through a conformationally flexible monomeric state, which is not observed in two near-identical mycobacterial orthologs. This morpheein behaviour exhibited by MtSerB2 lays the foundation for future allosteric drug discovery and provides a starting point to the understanding of its peculiar multifunctional moonlighting properties.

Protein Oligomerization

Protein self-association is a biologically remarkable event that involves and affects the structural and functional properties of proteins [...].

A proteome scale study reveals how plastic surfaces and agitation promote protein aggregation

Protein aggregation in biotherapeutics can reduce their activity and effectiveness. It may also promote immune reactions responsible for severe adverse effects. The impact of plastic materials on protein destabilization is not totally understood. Here, we propose to deconvolve the effects of material surface, air/liquid interface, and agitation to decipher their respective role in protein destabilization and aggregation. We analyzed the effect of polypropylene, TEFLON, glass and LOBIND surfaces on the stability of purified proteins (bovine serum albumin, hemoglobin and α-synuclein) and on a cell extract composed of 6000 soluble proteins during agitation (P = 0.1-1.2 W/kg). Proteomic analysis revealed that chaperonins, intrinsically disordered proteins and ribosomes were more sensitive to the combined effects of material surfaces and agitation while small metabolic oligomers could be protected in the same conditions. Protein loss observations coupled to Raman microscopy, dynamic light scattering and proteomic allowed us to propose a mechanistic model of protein destabilization by plastics. Our results suggest that protein loss is not primarily due to the nucleation of small aggregates in solution, but to the destabilization of proteins exposed to material surfaces and their subsequent aggregation at the sheared air/liquid interface, an effect that cannot be prevented by using LOBIND tubes. A guidance can be established on how to minimize these adverse effects. Remove one of the components of this combined stress - material, air (even partially), or agitation - and proteins will be preserved.

Predicting the Assembly of the Transmembrane Domains of Viral Channel Forming Proteins and Peptide Drug Screening Using a Docking Approach

A de novo assembly algorithm is provided to propose the assembly of bitopic transmembrane domains (TMDs) of membrane proteins. The algorithm is probed using, in particular, viral channel forming proteins (VCPs) such as M2 of influenza A virus, E protein of severe acute respiratory syndrome corona virus (SARS-CoV), 6K of Chikungunya virus (CHIKV), SH of human respiratory syncytial virus (hRSV), and Vpu of human immunodeficiency virus type 2 (HIV-2). The generation of the structures is based on screening a 7-dimensional space. Assembly of the TMDs can be achieved either by simultaneously docking the individual TMDs or via a sequential docking. Scoring based on estimated binding energies (EBEs) of the oligomeric structures is obtained by the tilt to decipher the handedness of the bundles. The bundles match especially well for all-atom models of M2 referring to an experimentally reported tetrameric bundle. Docking of helical poly-peptides to experimental structures of M2 and E protein identifies improving EBEs for positively charged (K,R,H) and aromatic amino acids (F,Y,W). Data are improved when using polypeptides for which the coordinates of the amino acids are adapted to the Cα coordinates of the respective experimentally derived structures of the TMDs of the target proteins.

Staphylococcus aureus Exfoliative Toxin E, Oligomeric State and Flip of P186: Implications for Its Action Mechanism

Staphylococcal exfoliative toxins (ETs) are glutamyl endopeptidases that specifically cleave the Glu381-Gly382 bond in the ectodomains of desmoglein 1 (Dsg1) via complex action mechanisms. To date, four ETs have been identified in different Staphylococcus aureus strains and ETE is the most recently characterized. The unusual properties of ETs have been attributed to a unique structural feature, i.e., the 180° flip of the carbonyl oxygen (O) of the nonconserved residue 192/186 (ETA/ETE numbering), not conducive to the oxyanion hole formation. We report the crystal structure of ETE determined at 1.61 Å resolution, in which P186(O) adopts two conformations displaying a 180° rotation. This finding, together with free energy calculations, supports the existence of a dynamic transition between the conformations under the tested conditions. Moreover, enzymatic assays showed no significant differences in the esterolytic efficiency of ETE and ETE/P186G, a mutant predicted to possess a functional oxyanion hole, thus downplaying the influence of the flip on the activity. Finally, we observed the formation of ETE homodimers in solution and the predicted homodimeric structure revealed the participation of a characteristic nonconserved loop in the interface and the partial occlusion of the protein active site, suggesting that monomerization is required for enzymatic activity.

A Lumenal Loop Associated with Catalytic Asymmetry in Plant Vacuolar H+-Translocating Pyrophosphatase

Membrane-integral inorganic pyrophosphatases (mPPases) couple pyrophosphate hydrolysis with H⁺ and Na⁺ pumping in plants and microbes. mPPases are homodimeric transporters with two catalytic sites facing the cytoplasm and demonstrating highly different substrate-binding affinities and activities. The structural aspects of the functional asymmetry are still poorly understood because the structure of the physiologically relevant dimer form with only one active site occupied by the substrate is unknown. We addressed this issue by molecular dynamics (MD) simulations of the H⁺-transporting mPPase of Vigna radiata, starting from its crystal structure containing a close substrate analog (imidodiphosphate, IDP) in both active sites. The MD simulations revealed pre-existing subunit asymmetry, which increased upon IDP binding to one subunit and persisted in the fully occupied dimer. The most significant asymmetrical change caused by IDP binding is a ‘rigid body’-like displacement of the lumenal loop connecting α-helices 2 and 3 in the partner subunit and opening its exit channel for water. This highly conserved 14–19-residue loop is found only in plant vacuolar mPPases and may have a regulatory function, such as pH sensing in the vacuole. Our data define the structural link between the loop and active sites and are consistent with the published structural and functional data.