Search

Although the αC-β4 loop is a stable feature of all protein kinases, the importance of this motif as a conserved element of secondary structure, as well as its links to the hydrophobic architecture of the kinase core, has been underappreciated. We first review the motif and then describe how it is linked to the hydrophobic spine architecture of the kinase core, which we first discovered using a computational tool, local spatial Pattern (LSP) alignment. Based on NMR predictions that a mutation in this motif abolishes the synergistic high-affinity binding of ATP and a pseudo substrate inhibitor, we used LSP to interrogate the F100A mutant. This comparison highlights the importance of the αC-β4 loop and key residues at the interface between the N- and C-lobes. In addition, we delved more deeply into the structure of the apo C-subunit, which lacks ATP. While apo C-subunit showed no significant changes in backbone dynamics of the αC-β4 loop, we found significant differences in the side chain dynamics of K105. The LSP analysis suggests disruption of communication between the N- and C-lobes in the F100A mutant, which would be consistent with the structural changes predicted by the NMR spectroscopy.

Contactin 2 (CNTN2) is a cell adhesion molecule involved in axon guidance, neuronal migration, and fasciculation. The ectodomains of CNTN1-CNTN6 are composed of six Ig domains (Ig1-Ig6) and four FN domains. Here, we show that CNTN2 forms transient homophilic interactions (KD ∼200 nM). Cryo-EM structures of full-length CNTN2 and CNTN2_Ig1-Ig6 reveal a T-shaped homodimer formed by intertwined, parallel monomers. Unexpectedly, the horseshoe-shaped Ig1-Ig4 headpieces extend their Ig2-Ig3 tips outwards on either side of the homodimer, while Ig4, Ig5, Ig6, and the FN domains form a central stalk. Cross-linking mass spectrometry and cell-based binding assays confirm the 3D assembly of the CNTN2 homodimer. The interface mediating homodimer formation differs between CNTNs, as do the homophilic versus heterophilic interaction mechanisms. The CNTN family thus encodes a versatile molecular platform that supports a very diverse portfolio of protein interactions and that can be leveraged to strategically guide neural circuit development.

Dynamic protein structures are crucial for deciphering their diverse biological functions. Two-dimensional infrared (2DIR) spectroscopy stands as an ideal tool for tracing rapid conformational evolutions in proteins. However, linking spectral characteristics to dynamic structures poses a formidable challenge. Here, we present a pretrained machine learning model based on 2DIR spectra analysis. This model has learned signal features from approximately 204,300 spectra to establish a spectrum-structure correlation, thereby tracing the dynamic conformations of proteins. It excels in accurately predicting the dynamic content changes of various secondary structures and demonstrates universal transferability on real folding trajectories spanning timescales from microseconds to milliseconds. Beyond exceptional predictive performance, the model offers attention-based spectral explanations of dynamic conformational changes. Our 2DIR-based pretrained model is anticipated to provide unique insights into the dynamic structural information of proteins in their native environments.

Cells employ many large macromolecular machines for the execution and regulation of processes that are vital for cell and organismal viability. Interestingly, cells cannot synthesize these machines as functioning units. Instead, cells synthesize the molecular parts that must then assemble into the functional complex. Many important machines, including chaperones such as GroEL and proteases such as the proteasome, comprise protein rings that are stacked on top of one another. While there is some experimental data regarding how stacked-ring complexes such as the proteasome self-assemble, a comprehensive understanding of the dynamics of stacked-ring assembly is currently lacking. Here, we developed a mathematical model of stacked-trimer assembly and performed an analysis of the assembly of the stacked homomeric trimer, which is the simplest stacked-ring architecture. We found that stacked rings are particularly susceptible to a form of kinetic trapping that we term deadlock, in which the system gets stuck in a state where there are many large intermediates that are not the fully assembled structure but that cannot productively react. When interaction affinities are uniformly strong, deadlock severely limits assembly yield. We thus predicted that stacked rings would avoid situations where all interfaces in the structure have high affinity. Analysis of available crystal structures indicated that indeed the majority-if not all-of stacked trimers do not contain uniformly strong interactions. Finally, to better understand the origins of deadlock, we developed a formal pathway analysis and showed that, when all the binding affinities are strong, many of the possible pathways are utilized. In contrast, optimal assembly strategies utilize only a small number of pathways. Our work suggests that deadlock is a critical factor influencing the evolution of macromolecular machines and provides general principles for understanding the self-assembly efficiency of existing machines.

Thermostable clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR-associated (Cas9) enzymes could improve genome-editing efficiency and delivery due to extended protein lifetimes. However, initial experimentation demonstrated Geobacillus stearothermophilus Cas9 (GeoCas9) to be virtually inactive when used in cultured human cells. Laboratory-evolved variants of GeoCas9 overcome this natural limitation by acquiring mutations in the wedge (WED) domain that produce >100-fold-higher genome-editing levels. Cryoelectron microscopy (cryo-EM) structures of the wild-type and improved GeoCas9 (iGeoCas9) enzymes reveal extended contacts between the WED domain of iGeoCas9 and DNA substrates. Biochemical analysis shows that iGeoCas9 accelerates DNA unwinding to capture substrates under the magnesium-restricted conditions typical of mammalian but not bacterial cells. These findings enabled rational engineering of other Cas9 orthologs to enhance genome-editing levels, pointing to a general strategy for editing enzyme improvement. Together, these results uncover a new role for the Cas9 WED domain in DNA unwinding and demonstrate how accelerated target unwinding dramatically improves Cas9-induced genome-editing activity.

The transcription factor and cell cycle regulator p53 is marked for degradation by the ubiquitin ligase MDM2. The interaction between these 2 proteins is mediated by a conserved binding motif in the disordered p53 transactivation domain (p53TAD) and the folded SWIB domain in MDM2. The conserved motif in p53TAD from zebrafish displays a 20-fold weaker interaction with MDM2, compared to the interaction in human and chicken. To investigate this apparent difference, we tracked the molecular evolution of the p53TAD/MDM2 interaction among ray-finned fishes (Actinopterygii), the largest vertebrate clade. Intriguingly, phylogenetic analyses, ancestral sequence reconstructions, and binding experiments showed that different loss-of-affinity changes in the canonical binding motif within p53TAD have occurred repeatedly and convergently in different fish lineages, resulting in relatively low extant affinities (KD = 0.5 to 5 μM). However, for 11 different fish p53TAD/MDM2 interactions, nonconserved regions flanking the canonical motif increased the affinity 4- to 73-fold to be on par with the human interaction. Our findings suggest that compensating changes at conserved and nonconserved positions within the motif, as well as in flanking regions of low conservation, underlie a stabilizing selection of functional affinity in the p53TAD/MDM2 interaction. Such interplay complicates bioinformatic prediction of binding and calls for experimental validation. Motif-mediated protein-protein interactions involving short binding motifs and folded interaction domains are very common across multicellular life. It is likely that the evolution of affinity in motif-mediated interactions often involves an interplay between specific interactions made by conserved motif residues and nonspecific interactions by nonconserved disordered regions.

Many peptide hormones form an α-helix on binding their receptors1-4, and sensitive methods for their detection could contribute to better clinical management of disease5. De novo protein design can now generate binders with high affinity and specificity to structured proteins6,7. However, the design of interactions between proteins and short peptides with helical propensity is an unmet challenge. Here we describe parametric generation and deep learning-based methods for designing proteins to address this challenge. We show that by extending RFdiffusion8 to enable binder design to flexible targets, and to refining input structure models by successive noising and denoising (partial diffusion), picomolar-affinity binders can be generated to helical peptide targets by either refining designs generated with other methods, or completely de novo starting from random noise distributions without any subsequent experimental optimization. The RFdiffusion designs enable the enrichment and subsequent detection of parathyroid hormone and glucagon by mass spectrometry, and the construction of bioluminescence-based protein biosensors. The ability to design binders to conformationally variable targets, and to optimize by partial diffusion both natural and designed proteins, should be broadly useful.

Complex II (CII) activity controls phenomena that require crosstalk between metabolism and signaling, including neurodegeneration, cancer metabolism, immune activation, and ischemia-reperfusion injury. CII activity can be regulated at the level of assembly, a process that leverages metastable assembly intermediates. The nature of these intermediates and how CII subunits transfer between metastable complexes remains unclear. In this work, we identify metastable species containing the SDHA subunit and its assembly factors, and we assign a preferred temporal sequence of appearance of these species during CII assembly. Structures of two species show that the assembly factors undergo disordered-to-ordered transitions without the appearance of significant secondary structure. The findings identify that intrinsically disordered regions are critical in regulating CII assembly, an observation that has implications for the control of assembly in other biomolecular complexes.

ImportanceInactivation of EP300/CREBB paralogous cellular lysine acetyltransferases (KATs) during the early phase of infection is a consistent feature of DNA viruses. The cell responds by stabilizing transcription factor IRF3 which activates transcription of scores of interferon-stimulated genes (ISGs), inhibiting viral replication. Human respiratory adenoviruses counter this by assembling a CUL4-based ubiquitin ligase complex that polyubiquitinylates RUVBL1 and 2 inducing their proteasomal degradation. This inhibits accumulation of active IRF3 and the expression of anti-viral ISGs, allowing replication of the respiratory HAdVs in the face of inhibition of EP300/CBEBBP KAT activity by the N-terminal region of E1A.

Iron delivery to the plasma is closely coupled to erythropoiesis, the production of red blood cells, as this process consumes most of the circulating plasma iron. In response to hemorrhage and other erythropoietic stresses, increased erythropoietin stimulates the production of the hormone erythroferrone (ERFE) by erythrocyte precursors (erythroblasts) developing in erythropoietic tissues. ERFE acts on the liver to inhibit bone morphogenetic protein (BMP) signaling and thereby decrease hepcidin production. Decreased circulating hepcidin concentrations then allow the release of iron from stores and increase iron absorption from the diet. Guided by evolutionary analysis and Alphafold2 protein complex modeling, we used targeted ERFE mutations, deletions, and synthetic ERFE segments together with cell-based bioassays and surface plasmon resonance to probe the structural features required for bioactivity and BMP binding. We define the ERFE active domain and multiple structural features that act together to entrap BMP ligands. In particular, the hydrophobic helical segment 81 to 86 and specifically the highly conserved tryptophan W82 in the N-terminal region are essential for ERFE bioactivity and Alphafold2 modeling places W82 between two tryptophans in its ligands BMP2, BMP6, and the BMP2/6 heterodimer, an interaction similar to those that bind BMPs to their cognate receptors. Finally, we identify the cationic region 96-107 and the globular TNFα-like domain 186-354 as structural determinants of ERFE multimerization that increase the avidity of ERFE for BMP ligands. Collectively, our results provide further insight into the ERFE-mediated inhibition of BMP signaling in response to erythropoietic stress.

CASP assessments primarily rely on comparing predicted coordinates with experimental reference structures. However, experimental structures by their nature are only models themselves-their construction involves a certain degree of subjectivity in interpreting density maps and translating them to atomic coordinates. Here, we directly utilized density maps to evaluate the predictions by employing a method for ranking the quality of protein chain predictions based on their fit into the experimental density. The fit-based ranking was found to correlate well with the CASP assessment scores. Overall, the evaluation against the density map indicated that the models are of high accuracy, and occasionally even better than the reference structure in some regions of the model. Local assessment of predicted side chains in a 1.52 Å resolution map showed that side-chains are sometimes poorly positioned. Additionally, the top 118 predictions associated with 9 protein target reference structures were selected for automated refinement, in addition to the top 40 predictions for 11 RNA targets. For both proteins and RNA, the refinement of CASP15 predictions resulted in structures that are close to the reference target structure. This refinement was successful despite large conformational changes often being required, showing that predictions from CASP-assessed methods could serve as a good starting point for building atomic models in cryo-EM maps for both proteins and RNA. Loop modeling continued to pose a challenge for predictors, and together with the lack of consensus amongst models in these regions suggests that modeling, in combination with model-fit to the density, holds the potential for identifying more flexible regions within the structure.

p97, also known as valosin-containing protein, is an essential cytosolic AAA+ (ATPases associated with diverse cellular activities) hexamer that unfolds substrate polypeptides to support protein homeostasis and macromolecular disassembly. Distinct sets of p97 adaptors guide cellular functions but their roles in direct control of the hexamer are unclear. The UBXD1 adaptor localizes with p97 in critical mitochondria and lysosome clearance pathways and contains multiple p97-interacting domains. Here we identify UBXD1 as a potent p97 ATPase inhibitor and report structures of intact human p97-UBXD1 complexes that reveal extensive UBXD1 contacts across p97 and an asymmetric remodeling of the hexamer. Conserved VIM, UBX and PUB domains tether adjacent protomers while a connecting strand forms an N-terminal domain lariat with a helix wedged at the interprotomer interface. An additional VIM-connecting helix binds along the second (D2) AAA+ domain. Together, these contacts split the hexamer into a ring-open conformation. Structures, mutagenesis and comparisons to other adaptors further reveal how adaptors containing conserved p97-remodeling motifs regulate p97 ATPase activity and structure.

For the first time, the 2022 CASP (Critical Assessment of Structure Prediction) community experiment included a section on computing multiple conformations for protein and RNA structures. There was full or partial success in reproducing the ensembles for four of the nine targets, an encouraging result. For protein structures, enhanced sampling with variations of the AlphaFold2 deep learning method was by far the most effective approach. One substantial conformational change caused by a single mutation across a complex interface was accurately reproduced. In two other assembly modeling cases, methods succeeded in sampling conformations near to the experimental ones even though environmental factors were not included in the calculations. An experimentally derived flexibility ensemble allowed a single accurate RNA structure model to be identified. Difficulties included how to handle sparse or low-resolution experimental data and the current lack of effective methods for modeling RNA/protein complexes. However, these and other obstacles appear addressable.

Prediction categories in the Critical Assessment of Structure Prediction (CASP) experiments change with the need to address specific problems in structure modeling. In CASP15, four new prediction categories were introduced: RNA structure, ligand-protein complexes, accuracy of oligomeric structures and their interfaces, and ensembles of alternative conformations. This paper lists technical specifications for these categories and describes their integration in the CASP data management system.

Processing of CASP15 targets into evaluation units (EUs) and assigning them to evolutionary-based prediction classes is presented in this study. The targets were first split into structural domains based on compactness and similarity to other proteins. Models were then evaluated against these domains and their combinations. The domains were joined into larger EUs if predictors' performance on the combined units was similar to that on individual domains. Alternatively, if most predictors performed better on the individual domains, then they were retained as EUs. As a result, 112 evaluation units were created from 77 tertiary structure prediction targets. The EUs were assigned to four prediction classes roughly corresponding to target difficulty categories in previous CASPs: TBM (template-based modeling, easy or hard), FM (free modeling), and the TBM/FM overlap category. More than a third of CASP15 EUs were attributed to the historically most challenging FM class, where homology or structural analogy to proteins of known fold cannot be detected.

Targeted intracellular delivery via receptor-mediated endocytosis requires the delivered cargo to escape the endosome to prevent lysosomal degradation. This can in principle be achieved by membrane lysis tightly restricted to endosomal membranes upon internalization to avoid general membrane insertion and lysis. Here, we describe the design of small monomeric proteins with buried histidine containing pH-responsive hydrogen bond networks and membrane permeating amphipathic helices. Of the 30 designs that were experimentally tested, all expressed in Escherichia coli, 13 were monomeric with the expected secondary structure, and 4 designs disrupted artificial liposomes in a pH-dependent manner. Mutational analysis showed that the buried histidine hydrogen bond networks mediate pH-responsiveness and control lysis of model membranes within a very narrow range of pH (6.0-5.5) with almost no lysis occurring at neutral pH. These tightly controlled lytic monomers could help mediate endosomal escape in designed targeted delivery platforms.

β-arrestins play a key role in G protein-coupled receptor (GPCR) internalization, trafficking, and signaling. Whether β-arrestins act independently of G protein-mediated signaling has not been fully elucidated. Studies using genome-editing approaches revealed that whereas G proteins are essential for mitogen-activated protein kinase activation by GPCRs., β-arrestins play a more prominent role in signal compartmentalization. However, in the absence of G proteins, GPCRs may not activate β-arrestins, thereby limiting the ability to distinguish G protein from β-arrestin-mediated signaling events. We used β2-adrenergic receptor (β2AR) and its β2AR-C tail mutant expressed in human embryonic kidney 293 cells wildtype or CRISPR-Cas9 gene edited for Gαs, β-arrestin1/2, or GPCR kinases 2/3/5/6 in combination with arrestin conformational sensors to elucidate the interplay between Gαs and β-arrestins in controlling gene expression. We found that Gαs is not required for β2AR and β-arrestin conformational changes, β-arrestin recruitment, and receptor internalization, but that Gαs dictates the GPCR kinase isoforms involved in β-arrestin recruitment. By RNA-Seq analysis, we found that protein kinase A and mitogen-activated protein kinase gene signatures were activated by stimulation of β2AR in wildtype and β-arrestin1/2-KO cells but absent in Gαs-KO cells. These results were validated by re-expressing Gαs in the corresponding KO cells and silencing β-arrestins in wildtype cells. These findings were extended to cellular systems expressing endogenous levels of β2AR. Overall, our results support that Gs is essential for β2AR-promoted protein kinase A and mitogen-activated protein kinase gene expression signatures, whereas β-arrestins initiate signaling events modulating Gαs-driven nuclear transcriptional activity.

p97/valosin-containing protein is an essential eukaryotic AAA+ ATPase with diverse functions including protein homeostasis, membrane remodeling, and chromatin regulation. Dysregulation of p97 function causes severe neurodegenerative disease and is associated with cancer, making this protein a significant therapeutic target. p97 extracts polypeptide substrates from macromolecular assemblies by hydrolysis-driven translocation through its central pore. Growing evidence indicates that this activity is highly coordinated by adapter partner proteins, of which more than 30 have been identified and are commonly described to facilitate translocation through substrate recruitment or modification. In so doing, these adapters enable critical p97-dependent functions such as extraction of misfolded proteins from the endoplasmic reticulum or mitochondria, and are likely the reason for the extreme functional diversity of p97 relative to other AAA+ translocases. Here, we review the known functions of adapter proteins and highlight recent structural and biochemical advances that have begun to reveal the diverse molecular bases for adapter-mediated regulation of p97 function. These studies suggest that the range of mechanisms by which p97 activity is controlled is vastly underexplored with significant advances possible for understanding p97 regulation by the most known adapters.

FAM86A is a class I lysine methyltransferase (KMT) that generates trimethylation on the eukaryotic translation elongation factor 2 (EEF2) at Lys525. Publicly available data from The Cancer Dependency Map project indicate high dependence of hundreds of human cancer cell lines on FAM86A expression. This classifies FAM86A among numerous other KMTs as potential targets for future anticancer therapies. However, selective inhibition of KMTs by small molecules can be challenging due to high conservation within the S-adenosyl methionine (SAM) cofactor binding domain among KMT subfamilies. Therefore, understanding the unique interactions within each KMT-substrate pair can facilitate developing highly specific inhibitors. The FAM86A gene encodes an N-terminal FAM86 domain of unknown function in addition to its C-terminal methyltransferase domain. Here, we used a combination of X-ray crystallography, the AlphaFold algorithms, and experimental biochemistry to identify an essential role of the FAM86 domain in mediating EEF2 methylation by FAM86A. To facilitate our studies, we also generated a selective EEF2K525 methyl antibody. Overall, this is the first report of a biological function for the FAM86 structural domain in any species and an example of a noncatalytic domain participating in protein lysine methylation. The interaction between the FAM86 domain and EEF2 provides a new strategy for developing a specific FAM86A small molecule inhibitor, and our results provide an example in which modeling a protein-protein interaction with AlphaFold expedites experimental biology.

We have investigated the equilibrium properties and rapid-reaction kinetics of the isolated butyryl-CoA dehydrogenase (bcd) component of the electron-bifurcating crotonyl-CoA-dependent NADH:ferredoxin oxidoreductase (EtfAB-bcd) from Megasphaera elsdenii. We find that a neutral FADH• semiquinone accumulates transiently during both reduction with sodium dithionite and with NADH in the presence of catalytic concentrations of EtfAB. In both cases full reduction of bcd to the hydroquinone is eventually observed, but the accumulation of FADH• indicates that a substantial portion of reduction occurs in sequential one-electron processes rather than a single two-electron event. In rapid-reaction experiments following the reaction of reduced bcd with crotonyl-CoA and oxidized bcd with butyryl-CoA, long-wavelength-absorbing intermediates are observed that are assigned to bcdred:crotonyl-CoA and bcdox:butyryl-CoA charge-transfer complexes, demonstrating their kinetic competence in the course of the reaction. In the presence of crotonyl-CoA there is an accumulation of semiquinone that is unequivocally the anionic FAD•- rather than the neutral FADH• seen in the absence of substrate, indicating that binding of substrate/product results in ionization of the bcd semiquinone. In addition to fully characterizing the rapid-reaction kinetics of both the oxidative and reductive half-reactions, our results demonstrate that one-electron processes play an important role in the reduction of bcd in EtfAB-bcd.

Sprouty-related EVH-1 domain-containing (SPRED) proteins are a family of proteins that negatively regulate the RAS-Mitogen-Activated Protein Kinase (MAPK) pathway, which is involved in the regulation of the mitogenic response and cell proliferation. However, the mechanism by which these proteins affect RAS-MAPK signaling has not been elucidated. Patients with mutations in SPRED give rise to unique disease phenotypes; thus, we hypothesized that distinct interactions across SPRED proteins may account for alternative nodes of regulation. To characterize the SPRED interactome and evaluate how members of the SPRED family function through unique binding partners, we performed affinity purification mass spectrometry. We identified 90-kDa ribosomal S6 kinase 2 (RSK2) as a specific interactor of SPRED2 but not SPRED1 or SPRED3. We identified that the N-terminal kinase domain of RSK2 mediates the interaction between amino acids 123 to 201 of SPRED2. Using X-ray crystallography, we determined the structure of the SPRED2-RSK2 complex and identified the SPRED2 motif, F145A, as critical for interaction. We found that the formation of this interaction is regulated by MAPK signaling events. We also find that this interaction between SPRED2 and RSK2 has functional consequences, whereby the knockdown of SPRED2 resulted in increased phosphorylation of RSK substrates, YB1 and CREB. Furthermore, SPRED2 knockdown hindered phospho-RSK membrane and nuclear subcellular localization. We report that disruption of the SPRED2-RSK complex has effects on RAS-MAPK signaling dynamics. Our analysis reveals that members of the SPRED family have unique protein binding partners and describes the molecular and functional determinants of SPRED2-RSK2 complex dynamics.

MDGAs (MAM domain-containing glycosylphosphatidylinositol anchors) are synaptic cell surface molecules that regulate the formation of trans-synaptic bridges between neurexins (NRXNs) and neuroligins (NLGNs), which promote synaptic development. Mutations in MDGAs are implicated in various neuropsychiatric diseases. MDGAs bind NLGNs in cis on the postsynaptic membrane and physically block NLGNs from binding to NRXNs. In crystal structures, the six immunoglobulin (Ig) and single fibronectin III domains of MDGA1 reveal a striking compact, triangular shape, both alone and in complex with NLGNs. Whether this unusual domain arrangement is required for biological function or other arrangements occur with different functional outcomes is unknown. Here, we show that WT MDGA1 can adopt both compact and extended 3D conformations that bind NLGN2. Designer mutants targeting strategic molecular elbows in MDGA1 alter the distribution of 3D conformations while leaving the binding affinity between soluble ectodomains of MDGA1 and NLGN2 intact. In contrast, in a cellular context, these mutants result in unique combinations of functional consequences, including altered binding to NLGN2, decreased capacity to conceal NLGN2 from NRXN1β, and/or suppressed NLGN2-mediated inhibitory presynaptic differentiation, despite the mutations being located far from the MDGA1-NLGN2 interaction site. Thus, the 3D conformation of the entire MDGA1 ectodomain appears critical for its function, and its NLGN-binding site on Ig1-Ig2 is not independent of the rest of the molecule. As a result, global 3D conformational changes to the MDGA1 ectodomain via strategic elbows may form a molecular mechanism to regulate MDGA1 action within the synaptic cleft.

We present a deep-learning-based platform, MIND-S, for protein post-translational modification (PTM) predictions. MIND-S employs a multi-head attention and graph neural network and assembles a 15-fold ensemble model in a multi-label strategy to enable simultaneous prediction of multiple PTMs with high performance and computation efficiency. MIND-S also features an interpretation module, which provides the relevance of each amino acid for making the predictions and is validated with known motifs. The interpretation module also captures PTM patterns without any supervision. Furthermore, MIND-S enables examination of mutation effects on PTMs. We document a workflow, its applications to 26 types of PTMs of two datasets consisting of ∼50,000 proteins, and an example of MIND-S identifying a PTM-interrupting SNP with validation from biological data. We also include use case analyses of targeted proteins. Taken together, we have demonstrated that MIND-S is accurate, interpretable, and efficient to elucidate PTM-relevant biological processes in health and diseases.

Amyloid protein aggregation is commonly associated with progressive neurodegenerative diseases, however not all amyloid fibrils are pathogenic. The neuronal cytoplasmic polyadenylation element binding protein is a regulator of synaptic mRNA translation and has been shown to form functional amyloid aggregates that stabilize long-term memory. In adult Drosophila neurons, the cytoplasmic polyadenylation element binding homolog Orb2 is expressed as 2 isoforms, of which the Orb2B isoform is far more abundant, but the rarer Orb2A isoform is required to initiate Orb2 aggregation. The N terminus is a distinctive feature of the Orb2A isoform and is critical for its aggregation. Intriguingly, replacement of phenylalanine in the fifth position of Orb2A with tyrosine (F5Y) in Drosophila impairs stabilization of long-term memory. The structure of endogenous Orb2B fibers was recently determined by cryo-EM, but the structure adopted by fibrillar Orb2A is less certain. Here we use micro-electron diffraction to determine the structure of the first 9 N-terminal residues of Orb2A, at a resolution of 1.05 Å. We find that this segment (which we term M9I) forms an amyloid-like array of parallel in-register β-sheets, which interact through side chain interdigitation of aromatic and hydrophobic residues. Our structure provides an explanation for the decreased aggregation observed for the F5Y mutant and offers a hypothesis for how the addition of a single atom (the tyrosyl oxygen) affects long-term memory. We also propose a structural model of Orb2A that integrates our structure of the M9I segment with the published Orb2B cryo-EM structure.

The power of structural information for informing biological mechanisms is clear for stable folded macromolecules, but similar structure-function insight is more difficult to obtain for highly dynamic systems such as intrinsically disordered proteins (IDPs) which must be described as structural ensembles. Here, we present IDPConformerGenerator, a flexible, modular open-source software platform for generating large and diverse ensembles of disordered protein states that builds conformers that obey geometric, steric, and other physical restraints on the input sequence. IDPConformerGenerator samples backbone phi (φ), psi (ψ), and omega (ω) torsion angles of relevant sequence fragments from loops and secondary structure elements extracted from folded protein structures in the RCSB Protein Data Bank and builds side chains from robust Monte Carlo algorithms using expanded rotamer libraries. IDPConformerGenerator has many user-defined options enabling variable fractional sampling of secondary structures, supports Bayesian models for assessing the agreement of IDP ensembles for consistency with experimental data, and introduces a machine learning approach to transform between internal and Cartesian coordinates with reduced error. IDPConformerGenerator will facilitate the characterization of disordered proteins to ultimately provide structural insights into these states that have key biological functions.

Antibiotics that use novel mechanisms are needed to combat antimicrobial resistance1-3. Teixobactin4 represents a new class of antibiotics with a unique chemical scaffold and lack of detectable resistance. Teixobactin targets lipid II, a precursor of peptidoglycan5. Here we unravel the mechanism of teixobactin at the atomic level using a combination of solid-state NMR, microscopy, in vivo assays and molecular dynamics simulations. The unique enduracididine C-terminal headgroup of teixobactin specifically binds to the pyrophosphate-sugar moiety of lipid II, whereas the N terminus coordinates the pyrophosphate of another lipid II molecule. This configuration favours the formation of a β-sheet of teixobactins bound to the target, creating a supramolecular fibrillar structure. Specific binding to the conserved pyrophosphate-sugar moiety accounts for the lack of resistance to teixobactin4. The supramolecular structure compromises membrane integrity. Atomic force microscopy and molecular dynamics simulations show that the supramolecular structure displaces phospholipids, thinning the membrane. The long hydrophobic tails of lipid II concentrated within the supramolecular structure apparently contribute to membrane disruption. Teixobactin hijacks lipid II to help destroy the membrane. Known membrane-acting antibiotics also damage human cells, producing undesirable side effects. Teixobactin damages only membranes that contain lipid II, which is absent in eukaryotes, elegantly resolving the toxicity problem. The two-pronged action against cell wall synthesis and cytoplasmic membrane produces a highly effective compound targeting the bacterial cell envelope. Structural knowledge of the mechanism of teixobactin will enable the rational design of improved drug candidates.

The EtfAB components of two bifurcating flavoprotein systems, the crotonyl-CoA-dependent NADH:ferredoxin oxidoreductase from the bacterium Megasphaera elsdenii and the menaquinone-dependent NADH:ferredoxin oxidoreductase from the archaeon Pyrobaculum aerophilum, have been investigated. With both proteins, we find that removal of the electron-transferring flavin adenine dinucleotide (FAD) moiety from both proteins results in an uncrossing of the reduction potentials of the remaining bifurcating FAD; this significantly stabilizes the otherwise very unstable semiquinone state, which accumulates over the course of reductive titrations with sodium dithionite. Furthermore, reduction of both EtfABs depleted of their electron-transferring FAD by NADH was monophasic with a hyperbolic dependence of reaction rate on the concentration of NADH. On the other hand, NADH reduction of the replete proteins containing the electron-transferring FAD was multiphasic, consisting of a fast phase comparable to that seen with the depleted proteins followed by an intermediate phase that involves significant accumulation of FAD⋅-, again reflecting uncrossing of the half-potentials of the bifurcating FAD. This is then followed by a slow phase that represents the slow reduction of the electron-transferring FAD to FADH-, with reduction of the now fully reoxidized bifurcating FAD by a second equivalent of NADH. We suggest that the crossing and uncrossing of the reduction half-potentials of the bifurcating FAD is due to specific conformational changes that have been structurally characterized.

Low-complexity domains (LCDs) of proteins have been shown to self-associate, and pathogenic mutations within these domains often drive the proteins into amyloid aggregation associated with disease. These domains may be especially susceptible to amyloidogenic mutations because they are commonly intrinsically disordered and function in self-association. The question therefore arises whether a search for pathogenic mutations in LCDs of the human proteome can lead to identification of other proteins associated with amyloid disease. Here, we take a computational approach to identify documented pathogenic mutations within LCDs that may favor amyloid formation. Using this approach, we identify numerous known amyloidogenic mutations, including several such mutations within proteins previously unidentified as amyloidogenic. Among the latter group, we focus on two mutations within the TRK-fused gene protein (TFG), known to play roles in protein secretion and innate immunity, which are associated with two different peripheral neuropathies. We show that both mutations increase the propensity of TFG to form amyloid fibrils. We therefore conclude that TFG is a novel amyloid protein and propose that the diseases associated with its mutant forms may be amyloidoses.

Fitness functions map biological sequences to a scalar property of interest. Accurate estimation of these functions yields biological insight and sets the foundation for model-based sequence design. However, the fitness datasets available to learn these functions are typically small relative to the large combinatorial space of sequences; characterizing how much data are needed for accurate estimation remains an open problem. There is a growing body of evidence demonstrating that empirical fitness functions display substantial sparsity when represented in terms of epistatic interactions. Moreover, the theory of Compressed Sensing provides scaling laws for the number of samples required to exactly recover a sparse function. Motivated by these results, we develop a framework to study the sparsity of fitness functions sampled from a generalization of the NK model, a widely used random field model of fitness functions. In particular, we present results that allow us to test the effect of the Generalized NK (GNK) model's interpretable parameters-sequence length, alphabet size, and assumed interactions between sequence positions-on the sparsity of fitness functions sampled from the model and, consequently, the number of measurements required to exactly recover these functions. We validate our framework by demonstrating that GNK models with parameters set according to structural considerations can be used to accurately approximate the number of samples required to recover two empirical protein fitness functions and an RNA fitness function. In addition, we show that these GNK models identify important higher-order epistatic interactions in the empirical fitness functions using only structural information.

The three-dimensional structure of the synthetic lung Surfactant Protein B Peptide Super Mini-B was determined using an integrative experimental approach, including mass spectrometry and isotope enhanced Fourier-transform infrared (FTIR) spectroscopy. Mass spectral analysis of the peptide, oxidized by solvent assisted region-specific disulfide formation, confirmed that the correct folding and disulfide pairing could be facilitated using two different oxidative structure-promoting solvent systems. Residue specific analysis by isotope enhanced FTIR indicated that the N-terminal and C-terminal domains have well defined α-helical amino acid sequences. Using these experimentally derived measures of distance constraints and disulfide connectivity, the ensemble was further refined with molecular dynamics to provide a medium resolution, residue-specific structure for the peptide construct in a simulated synthetic lung surfactant lipid multilayer environment. The disulfide connectivity combined with the α-helical elements stabilize the peptide conformationally to form a helical hairpin structure that resembles critical elements of the Saposin protein fold of the predicted full-length Surfactant Protein B structure.

To understand the dynamic structure-function relationship of soft- and biomolecules, the determination of the three-dimensional (3D) structure of each individual molecule (nonaveraged structure) in its native state is sought-after. Cryo-electron tomography (cryo-ET) is a unique tool for imaging an individual object from a series of tilted views. However, due to radiation damage from the incident electron beam, the tolerable electron dose limits image contrast and the signal-to-noise ratio (SNR) of the data, preventing the 3D structure determination of individual molecules, especially at high-resolution. Although recently developed technologies and techniques, such as the direct electron detector, phase plate, and computational algorithms, can partially improve image contrast/SNR at the same electron dose, the high-resolution structure, such as tertiary structure of individual molecules, has not yet been resolved. Here, we review the cryo-electron microscopy (cryo-EM) and cryo-ET experimental parameters to discuss how these parameters affect the extent of radiation damage. This discussion can guide us in optimizing the experimental strategy to increase the imaging dose or improve image SNR without increasing the radiation damage. With a higher dose, a higher image contrast/SNR can be achieved, which is crucial for individual-molecule 3D structure. With 3D structures determined from an ensemble of individual molecules in different conformations, the molecular mechanism through their biochemical reactions, such as self-folding or synthesis, can be elucidated in a straightforward manner.

Single particle cryogenic-electron microscopy (cryo-EM) is used extensively to determine structures of activated G protein-coupled receptors (GPCRs) in complex with G proteins or arrestins. However, applying it to GPCRs without signaling proteins remains challenging because most receptors lack structural features in their soluble domains to facilitate image alignment. In GPCR crystallography, inserting a fusion protein between transmembrane helices 5 and 6 is a highly successful strategy for crystallization. Although a similar strategy has the potential to broadly facilitate cryo-EM structure determination of GPCRs alone without signaling protein, the critical determinants that make this approach successful are not yet clear. Here, we address this shortcoming by exploring different fusion protein designs, which lead to structures of antagonist bound A2A adenosine receptor at 3.4 Å resolution and unliganded Smoothened at 3.7 Å resolution. The fusion strategies explored here are likely applicable to cryo-EM interrogation of other GPCRs and small integral membrane proteins.

Structures provide a critical breakthrough step for biological analyses, and small angle X-ray scattering (SAXS) is a powerful structural technique to study dynamic DNA repair proteins. As toxic and mutagenic repair intermediates need to be prevented from inadvertently harming the cell, DNA repair proteins often chaperone these intermediates through dynamic conformations, coordinated assemblies, and allosteric regulation. By measuring structural conformations in solution for both proteins, DNA, RNA, and their complexes, SAXS provides insight into initial DNA damage recognition, mechanisms for validation of their substrate, and pathway regulation. Here, we describe exemplary SAXS analyses of a DNA damage response protein spanning from what can be derived directly from the data to obtaining super resolution through the use of SAXS selection of atomic models. We outline strategies and tactics for practical SAXS data collection and analysis. Making these structural experiments in reach of any basic and clinical researchers who have protein, SAXS data can readily be collected at government-funded synchrotrons, typically at no cost for academic researchers. In addition to discussing how SAXS complements and enhances cryo-electron microscopy, X-ray crystallography, NMR, and computational modeling, we furthermore discuss taking advantage of recent advances in protein structure prediction in combination with SAXS analysis.

PurposeTo report novel genotypes and expand the phenotype spectrum of SSBP1-disease and explore potential disease mechanism.MethodsFive families with previously unsolved optic atrophy and retinal dystrophy underwent whole genome sequencing as part of the National Institute for Health Research BioResource Rare-Diseases and the UK's 100,000 Genomes Project. In silico analysis and protein modelling was performed on the identified variants. Deep phenotyping including retinal imaging and International Society for Clinical Electrophysiology of Vision standard visual electrophysiology was performed.ResultsSeven individuals from five unrelated families with bilateral optic atrophy and/or retinal dystrophy with extraocular signs and symptoms in some are described. In total, 6 SSBP1 variants were identified including the previously unreported variants: c.151A>G, p.(Lys51Glu), c.335G>A p.(Gly112Glu), and c.380G>A, p.(Arg127Gln). One individual was found to carry biallelic variants (c.380G>A p.(Arg127Gln); c.394A>G p.(Ile132Val)) associated with likely autosomal recessive SSBP1-disease. In silico analysis predicted all variants to be pathogenic and Three-dimensional protein modelling suggested possible disease mechanisms via decreased single-stranded DNA binding affinity or impaired higher structure formation.ConclusionsSSBP1 is essential for mitochondrial DNA replication and maintenance, with defects leading to a spectrum of disease that includes optic atrophy and/or retinal dystrophy, occurring with or without extraocular features. This study provides evidence of intrafamilial variability and confirms the existence of an autosomal recessive inheritance in SSBP1-disease consequent upon a previously unreported genotype.

In CASP14, 39 research groups submitted more than 2500 3D models on 22 protein complexes. In general, the community performed well in predicting the fold of the assemblies (for 80% of the targets), although it faced significant challenges in reproducing the native contacts. This is especially the case for the complexes without whole-assembly templates. The leading predictor, BAKER-experimental, used a methodology combining classical techniques (template-based modeling, protein docking) with deep learning-based contact predictions and a fold-and-dock approach. The Venclovas team achieved the runner-up position with template-based modeling and docking. By analyzing the target interfaces, we showed that the complexes with depleted charged contacts or dominating hydrophobic interactions were the most challenging ones to predict. We also demonstrated that if AlphaFold2 predictions were at hand, the interface prediction challenge could be alleviated for most of the targets. All in all, it is evident that new approaches are needed for the accurate prediction of assemblies, which undoubtedly will expand on the significant improvements in the tertiary structure prediction field.

An evolutionary-based definition and classification of target evaluation units (EUs) is presented for the 14th round of the critical assessment of structure prediction (CASP14). CASP14 targets included 84 experimental models submitted by various structural groups (designated T1024-T1101). Targets were split into EUs based on the domain organization of available templates and performance of server groups. Several targets required splitting (19 out of 25 multidomain targets) due in part to observed conformation changes. All in all, 96 CASP14 EUs were defined and assigned to tertiary structure assessment categories (Topology-based FM or High Accuracy-based TBM-easy and TBM-hard) considering their evolutionary relationship to existing ECOD fold space: 24 family level, 50 distant homologs (H-group), 12 analogs (X-group), and 10 new folds. Principal component analysis and heatmap visualization of sequence and structure similarity to known templates as well as performance of servers highlighted trends in CASP14 target difficulty. The assigned evolutionary levels (i.e., H-groups) and assessment classes (i.e., FM) displayed overlapping clusters of EUs. Many viral targets diverged considerably from their template homologs and thus were more difficult for prediction than other homology-related targets. On the other hand, some targets did not have sequence-identifiable templates, but were predicted better than expected due to relatively simple arrangements of secondary structural elements. An apparent improvement in overall server performance in CASP14 further complicated traditional classification, which ultimately assigned EUs into high-accuracy modeling (27 TBM-easy and 31 TBM-hard), topology (23 FM), or both (15 FM/TBM).

The high accuracy of some CASP14 models at the domain level prompted a more detailed evaluation of structure predictions on whole targets. For the first time in critical assessment of structure prediction (CASP), we evaluated accuracy of difficult domain assembly in models submitted for multidomain targets where the community predicted individual evaluation units (EUs) with greater accuracy than full-length targets. Ten proteins with domain interactions that did not show evidence of conformational change and were not involved in significant oligomeric contacts were chosen as targets for the domain interaction assessment. Groups were ranked using complementary interaction scores (F1, QS score, and Jaccard coefficient), and their predictions were evaluated for their ability to correctly model inter-domain interfaces and overall protein folds. Target performance was broadly grouped into two clusters. The first consisted primarily of targets containing two EUs wherein predictors more broadly predicted domain positioning and interfacial contacts correctly. The other consisted of complex two-EU and three-EU targets where few predictors performed well. The highest ranked predictor, AlphaFold2, produced high-accuracy models on eight out of 10 targets. Their interdomain scores on three of these targets were significantly higher than all other groups and were responsible for their overall outperformance in the category. We further highlight the performance of AlphaFold2 and the next best group, BAKER-experimental on several interesting targets.

Retinal guanylate cyclases (RetGCs) are regulated by a family of guanylate cyclase-activating proteins (called GCAP1-7). GCAPs form dimers that bind to Ca2+ and confer Ca2+ sensitive activation of RetGC during visual phototransduction. The GCAP5 homologue from zebrafish contains two nonconserved cysteine residues (Cys15 and Cys17) that bind to ferrous ion, which stabilizes GCAP5 dimerization and diminishes its ability to activate RetGC. Here, we present NMR and EPR-DEER structural analysis of a GCAP5 dimer in the Mg2+-bound, Ca2+-free, Fe2+-free activator state. The NMR-derived structure of GCAP5 is similar to the crystal structure of Ca2+-bound GCAP1 (root-mean-square deviation of 2.4 Å), except that the N-terminal helix of GCAP5 is extended by two residues, which allows the sulfhydryl groups of Cys15 and Cys17 to become more solvent exposed in GCAP5 to facilitate Fe2+ binding. Nitroxide spin-label probes were covalently attached to particular cysteine residues engineered in GCAP5: C15, C17, T26C, C28, N56C, C69, C105, N139C, E152C, and S159C. The intermolecular distance of each spin-label probe in dimeric GCAP5 (measured by EPR-DEER) defined restraints for calculating the dimer structure by molecular docking. The GCAP5 dimer possesses intermolecular hydrophobic contacts involving the side chain atoms of H18, Y21, M25, F72, V76, and W93, as well as an intermolecular salt bridge between R22 and D71. The structural model of the GCAP5 dimer was validated by mutations (H18E/Y21E, H18A/Y21A, R22D, R22A, M25E, D71R, F72E, and V76E) at the dimer interface that disrupt dimerization of GCAP5 and affect the activation of RetGC. We propose that GCAP5 dimerization may play a role in the Fe2+-dependent regulation of cyclase activity in zebrafish photoreceptors.

Although β-hairpins are widespread in proteins, there is no tool to coax any small peptide to adopt a β-hairpin conformation, regardless of sequence. Here, we report that δ-linked γ(R)-methyl-ornithine (δ MeOrn) provides an improved β-turn template for inducing a β-hairpin conformation in peptides. We developed a synthesis of protected δ MeOrn as a building block suitable for use in Fmoc-based solid-phase peptide synthesis. The synthesis begins with l-leucine and affords gram quantities of the Nα -Boc-Nδ -Fmoc-γ(R)-methyl-ornithine building block. X-ray crystallography confirms that the δ MeOrn turn unit adopts a folded structure in a macrocyclic β-hairpin peptide. CD and NMR spectroscopy allow comparison of the δ MeOrn turn template to the δ-linked ornithine (δ Orn) turn template that we previously introduced and to the popular d-Pro-Gly turn template. These studies show that the folding of the δ MeOrn turn template is substantially better than that of δ Orn and is comparable to d-Pro-Gly.

We describe PROPER-seq (protein-protein interaction sequencing) to map protein-protein interactions (PPIs) en masse. PROPER-seq first converts transcriptomes of input cells into RNA-barcoded protein libraries, in which all interacting protein pairs are captured through nucleotide barcode ligation, recorded as chimeric DNA sequences, and decoded at once by sequencing and mapping. We applied PROPER-seq to human embryonic kidney cells, T lymphocytes, and endothelial cells and identified 210,518 human PPIs (collected in the PROPER v.1.0 database). Among these, 1,365 and 2,480 PPIs are supported by published co-immunoprecipitation (coIP) and affinity purification-mass spectrometry (AP-MS) data, 17,638 PPIs are predicted by the prePPI algorithm without previous experimental validation, and 100 PPIs overlap human synthetic lethal gene pairs. In addition, four previously uncharacterized interaction partners with poly(ADP-ribose) polymerase 1 (PARP1) (a critical protein in DNA repair) known as XPO1, MATR3, IPO5, and LEO1 are validated in vivo. PROPER-seq presents a time-effective technology to map PPIs at the transcriptome scale, and PROPER v.1.0 provides a rich resource for studying PPIs.

Extant fold-switching proteins remodel their secondary structures and change their functions in response to cellular stimuli, regulating biological processes and affecting human health. Despite their biological importance, these proteins remain understudied. Predictive methods are needed to expedite the process of discovering and characterizing more of these shapeshifting proteins. Most previous approaches require a solved structure or all-atom simulations, greatly constraining their use. Here, we propose a high-throughput sequence-based method for predicting extant fold switchers that transition from α-helix in one conformation to β-strand in the other. This method leverages two previous observations: (a) α-helix ↔ β-strand prediction discrepancies from JPred4 are a robust predictor of fold switching, and (b) the fold-switching regions (FSRs) of some extant fold switchers have different secondary structure propensities when expressed by themselves (isolated FSRs) than when expressed within the context of their parent protein (contextualized FSRs). Combining these two observations, we ran JPred4 on 99-fold-switching proteins and found strong correspondence between predicted and experimentally observed α-helix ↔ β-strand discrepancies. To test the overall robustness of this finding, we randomly selected regions of proteins not expected to switch folds (single-fold proteins) and found significantly fewer predicted α-helix ↔ β-strand discrepancies. Combining these discrepancies with the overall percentage of predicted secondary structure, we developed a classifier to identify extant fold switchers (Matthews correlation coefficient of .71). Although this classifier had a high false-negative rate (7/17), its false-positive rate was very low (2/136), suggesting that it can be used to predict a subset of extant fold switchers from a multitude of available genomic sequences.

Although most experimentally characterized proteins with similar sequences assume the same folds and perform similar functions, an increasing number of exceptions is emerging. One class of exceptions comprises sequence-similar fold switchers, whose secondary structures shift from α-helix β-sheet through a small number of mutations, a sequence insertion, or a deletion. Predictive methods for identifying sequence-similar fold switchers are desirable because some are associated with disease and/or can perform different functions in cells. Here, we use homology-based secondary structure predictions to identify sequence-similar fold switchers from their amino acid sequences alone. To do this, we predicted the secondary structures of sequence-similar fold switchers using three different homology-based secondary structure predictors: PSIPRED, JPred4, and SPIDER3. We found that α-helix β-strand prediction discrepancies from JPred4 discriminated between the different conformations of sequence-similar fold switchers with high statistical significance (P

TRAAK is a mechanosensitive two-pore domain K+ (K2P) channel localized to nodes of Ranvier in myelinated neurons. TRAAK deletion in mice results in mechanical and thermal allodynia, and gain-of-function mutations cause the human neurodevelopmental disorder FHEIG. TRAAK displays basal and stimulus-gated activities typical of K2Ps, but the mechanistic and structural differences between these modes are unknown. Here, we demonstrate that basal and mechanically gated openings are distinguished by their conductance, kinetics, and structure. Basal openings are low conductance, short duration, and due to a conductive channel conformation with the interior cavity exposed to the surrounding membrane. Mechanically gated openings are high conductance, long duration, and due to a channel conformation in which the interior cavity is sealed to the surrounding membrane. Our results explain how dual modes of activity are produced by a single ion channel and provide a basis for the development of state-selective pharmacology with the potential to treat disease.

Time-resolved crystallography of biomolecules in action has advanced rapidly as methods for serial crystallography have improved, but the large number of crystals and the complex experimental infrastructure that are required remain serious obstacles to its widespread application. Here, millisecond mix-and-quench crystallography (MMQX) has been developed, which yields millisecond time-resolved data using far fewer crystals and routine remote synchrotron data collection. To demonstrate the capabilities of MMQX, the conversion of oxaloacetic acid to phosphoenolpyruvate by phosphoenolpyruvate carboxy-kinase (PEPCK) is observed with a time resolution of 40 ms. By lowering the entry barrier to time-resolved crystallography, MMQX should enable a broad expansion in structural studies of protein dynamics.

Tyrosyl DNA phosphodiesterase 1 (TDP1) and DNA Ligase IIIα (LigIIIα) are key enzymes in single-strand break (SSB) repair. TDP1 removes 3'-tyrosine residues remaining after degradation of DNA topoisomerase (TOP) 1 cleavage complexes trapped by either DNA lesions or TOP1 inhibitors. It is not known how TDP1 is linked to subsequent processing and LigIIIα-catalyzed joining of the SSB. Here we define a direct interaction between the TDP1 catalytic domain and the LigIII DNA-binding domain (DBD) regulated by conformational changes in the unstructured TDP1 N-terminal region induced by phosphorylation and/or alterations in amino acid sequence. Full-length and N-terminally truncated TDP1 are more effective at correcting SSB repair defects in TDP1 null cells compared with full-length TDP1 with amino acid substitutions of an N-terminal serine residue phosphorylated in response to DNA damage. TDP1 forms a stable complex with LigIII170-755, as well as full-length LigIIIα alone or in complex with the DNA repair scaffold protein XRCC1. Small-angle X-ray scattering and negative stain electron microscopy combined with mapping of the interacting regions identified a TDP1/LigIIIα compact dimer of heterodimers in which the two LigIII catalytic cores are positioned in the center, whereas the two TDP1 molecules are located at the edges of the core complex flanked by highly flexible regions that can interact with other repair proteins and SSBs. As TDP1and LigIIIα together repair adducts caused by TOP1 cancer chemotherapy inhibitors, the defined interaction architecture and regulation of this enzyme complex provide insights into a key repair pathway in nonmalignant and cancer cells.

OBJECTIVE: Genetic variants in the GRIN genes that encode N-methyl-D-aspartate receptor (NMDAR) subunits have been identified in various neurodevelopmental disorders, including epilepsy. We identified a GRIN1 variant from an individual with early-onset epileptic encephalopathy, evaluated functional changes to NMDAR properties caused by the variant, and screened FDA-approved therapeutic compounds as potential treatments for the patient. METHODS: Whole exome sequencing identified a missense variant in GRIN1. Electrophysiological recordings were made from Xenopus oocytes and transfected HEK cells to determine the NMDAR biophysical properties as well as the sensitivity to agonists and FDA-approved drugs that inhibit NMDARs. A beta-lactamase reporter assay in transfected HEK cells evaluated the effects of the variant on the NMDAR surface expression. RESULTS: A recurrent de novo missense variant in GRIN1 (c.1923G>A, p.Met641Ile), which encodes the GluN1 subunit, was identified in a pediatric patient with drug-resistant seizures and early-onset epileptic encephalopathy. In vitro analysis indicates that GluN1-M641I containing NMDARs showed enhanced agonist potency and reduced Mg2+ block, which may be associated with the patients phenotype. Results from screening FDA-approved drugs suggested that GluN1-M641I containing NMDARs are more sensitive to the NMDAR channel blockers memantine, ketamine, and dextromethorphan compared to the wild-type receptors. The addition of memantine to the seizure treatment regimen significantly reduced the patients seizure burden. INTERPRETATION: Our finding contributes to the understanding of the phenotype-genotype correlations of patients with GRIN1 gene variants, provides a molecular mechanism underlying the actions of this variant, and explores therapeutic strategies for treating GRIN1-related neurological conditions.

Mutation of C9ORF72 is the most common genetic cause of amyotrophic lateral sclerosis (ALS) and frontal temporal degeneration (FTD), which is attributed to both a gain and loss of function. C9orf72 forms a complex with SMCR8 and WDR41, which was reported to have GTPase activating protein activity toward ARF proteins, RAB8A, and RAB11A. We determined the cryo-EM structure of ARF1-GDP-BeF3- bound to C9orf72:SMCR8:WDR41. The SMCR8longin and C9orf72longin domains form the binding pocket for ARF1. One face of the C9orf72longin domain holds ARF1 in place, while the SMCR8longin positions the catalytic finger Arg147 in the ARF1 active site. Mutations in interfacial residues of ARF1 and C9orf72 reduced or eliminated GAP activity. RAB8A GAP required ~10-fold higher concentrations of the C9orf72 complex than for ARF1. These data support a specific function for the C9orf72 complex as an ARF GAP. The structure also provides a model for the active forms of the longin domain GAPs of FLCN and NPRL2 that regulate the Rag GTPases of the mTORC1 pathway.

The structures of Aspergillus oryzae α-amylase were determined in a tetragonal crystal, having one molecule as asymmetric unit, and a monoclinic crystal with two molecules as asymmetric unit. Both crystal forms were obtained from trace contaminants of an old commercial lipase preparation. Structures were determined and refined to 1.65 Å and 1.43 Å resolution respectively. The latter crystal has a non-crystallographic (NCS) twofold axis within the asymmetric unit. Glycosylation at Asn197 is evident, and in the tetragonal crystal can be seen to include three, partially disordered sugar residues following the initial N-acetyl glucosamine (NAG). Superposition of the tetragonal crystal model on the α-amylases from Bacillus subtilis (PDB:1BAG), pig pancreas (PDB:3L2L), and barley (PDB:1AMY), show a high degree of coincidence, particularly for the (β/α)8-barrel domains, and especially within the active site. Using this structural agreement between amylases, we extrapolated the binding model of a six residue, limit dextrin found in pig pancreas α-amylase to the A. oryzae enzyme model, which predicts substrate interacting amino acid residues.

Infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is initiated by binding of the viral Spike protein to host receptor angiotensin-converting enzyme 2 (ACE2), followed by fusion of viral and host membranes. Although antibodies that block this interaction are in emergency use as early coronavirus disease 2019 (COVID-19) therapies, the precise determinants of neutralization potency remain unknown. We discovered a series of antibodies that potently block ACE2 binding but exhibit divergent neutralization efficacy against the live virus. Strikingly, these neutralizing antibodies can inhibit or enhance Spike-mediated membrane fusion and formation of syncytia, which are associated with chronic tissue damage in individuals with COVID-19. As revealed by cryoelectron microscopy, multiple structures of Spike-antibody complexes have distinct binding modes that not only block ACE2 binding but also alter the Spike protein conformational cycle triggered by ACE2 binding. We show that stabilization of different Spike conformations leads to modulation of Spike-mediated membrane fusion with profound implications for COVID-19 pathology and immunity.

As organelles of the innate immune system, inflammasomes activate caspase-1 and other inflammatory caspases that cleave gasdermin D (GSDMD). Caspase-1 also cleaves inactive precursors of the interleukin (IL)-1 family to generate mature cytokines such as IL-1β and IL-18. Cleaved GSDMD forms transmembrane pores to enable the release of IL-1 and to drive cell lysis through pyroptosis1-9. Here we report cryo-electron microscopy structures of the pore and the prepore of GSDMD. These structures reveal the different conformations of the two states, as well as extensive membrane-binding elements including a hydrophobic anchor and three positively charged patches. The GSDMD pore conduit is predominantly negatively charged. By contrast, IL-1 precursors have an acidic domain that is proteolytically removed by caspase-110. When permeabilized by GSDMD pores, unlysed liposomes release positively charged and neutral cargoes faster than negatively charged cargoes of similar sizes, and the pores favour the passage of IL-1β and IL-18 over that of their precursors. Consistent with these findings, living-but not pyroptotic-macrophages preferentially release mature IL-1β upon perforation by GSDMD. Mutation of the acidic residues of GSDMD compromises this preference, hindering intracellular retention of the precursor and secretion of the mature cytokine. The GSDMD pore therefore mediates IL-1 release by electrostatic filtering, which suggests the importance of charge in addition to size in the transport of cargoes across this large channel.