X‐ray scattering reveals disordered linkers and dynamic interfaces in complexes and mechanisms for DNA double‐strand break repair impacting cell and cancer biology

Evolutionary selection ensures specificity and efficiency in dynamic metastable macromolecular machines that repair DNA damage without releasing toxic and mutagenic intermediates. Here we examine non‐homologous end joining (NHEJ) as the primary conserved DNA double‐strand break (DSB) repair process in human cells. NHEJ has exemplary key roles in networks determining the development, outcome of cancer treatments by DSB‐inducing agents, generation of antibody and T‐cell receptor diversity, and innate immune response for RNA viruses. We determine mechanistic insights into NHEJ structural biochemistry focusing upon advanced small angle X‐ray scattering (SAXS) results combined with X‐ray crystallography (MX) and cryo‐electron microscopy (cryo‐EM). SAXS coupled to atomic structures enables integrated structural biology for objective quantitative assessment of conformational ensembles and assemblies in solution, intra‐molecular distances, structural similarity, functional disorder, conformational switching, and flexibility. Importantly, NHEJ complexes in solution undergo larger allosteric transitions than seen in their cryo‐EM or MX structures. In the long‐range synaptic complex, X‐ray repair cross‐complementing 4 (XRCC4) plus XRCC4‐like‐factor (XLF) form a flexible bridge and linchpin for DNA ends bound to KU heterodimer (Ku70/80) and DNA‐PKcs (DNA‐dependent protein kinase catalytic subunit). Upon binding two DNA ends, auto‐phosphorylation opens DNA‐PKcs dimer licensing NHEJ via concerted conformational transformations of XLF‐XRCC4, XLF–Ku80, and LigIV^BRCT–Ku70 interfaces. Integrated structures reveal multifunctional roles for disordered linkers and modular dynamic interfaces promoting DSB end processing and alignment into the short‐range complex for ligation by LigIV. Integrated findings define dynamic assemblies fundamental to designing separation‐of‐function mutants and allosteric inhibitors targeting conformational transitions in multifunctional complexes.     

Cryo‐EM structure of metazoan TRAPPIII,the multi‐subunit complex that activates the GTPase Rab1

The TRAPP complexes are nucleotide exchange factors that play essential roles in membrane traffic and autophagy. TRAPPII activates Rab11, and TRAPPIII activates Rab1, with the two complexes sharing a core of small subunits that affect nucleotide exchange but being distinguished by specific large subunits that are essential for activity in vivo. Crystal structures of core subunits have revealed the mechanism of Rab activation, but how the core and the large subunits assemble to form the complexes is unknown. We report a cryo‐EM structure of the entire Drosophila TRAPPIII complex. The TRAPPIII‐specific subunits TRAPPC8 and TRAPPC11 hold the catalytic core like a pair of tongs, with TRAPPC12 and TRAPPC13 positioned at the joint between them. TRAPPC2 and TRAPPC2L link the core to the two large arms, with the interfaces containing residues affected by disease‐causing mutations. The TRAPPC8 arm is positioned such that it would contact Rab1 that is bound to the core, indicating how the arm could determine the specificity of the complex. A lower resolution structure of TRAPPII shows a similar architecture and suggests that the TRAPP complexes evolved from a single ur‐TRAPP.     

Cryo‐EM reveals the complex architecture of dynactin's shoulder region and pointed end

Dynactin is a 1.1 MDa complex that activates the molecular motor dynein for ultra‐processive transport along microtubules. In order to do this, it forms a tripartite complex with dynein and a coiled‐coil adaptor. Dynactin consists of an actin‐related filament whose length is defined by its flexible shoulder domain. Despite previous cryo‐EM structures, the molecular architecture of the shoulder and pointed end of the filament is still poorly understood due to the lack of high‐resolution information in these regions. Here we combine multiple cryo‐EM datasets and define precise masking strategies for particle signal subtraction and 3D classification. This overcomes domain flexibility and results in high‐resolution maps into which we can build the shoulder and pointed end. The unique architecture of the shoulder securely houses the p150 subunit and positions the four identical p50 subunits in different conformations to bind dynactin’s filament. The pointed end map allows us to build the first structure of p62 and reveals the molecular basis for cargo adaptor binding to different sites at the pointed end.     

Selective cross‐linking of coinciding protein assemblies by in‐gel cross‐linking mass spectrometry

Cross‐linking mass spectrometry has developed into an important method to study protein structures and interactions. The in‐solution cross‐linking workflows involve time and sample consuming steps and do not provide sensible solutions for differentiating cross‐links obtained from co‐occurring protein oligomers, complexes, or conformers. Here we developed a cross‐linking workflow combining blue native PAGE with in‐gel cross‐linking mass spectrometry (IGX‐MS). This workflow circumvents steps, such as buffer exchange and cross‐linker concentration optimization. Additionally, IGX‐MS enables the parallel analysis of co‐occurring protein complexes using only small amounts of sample. Another benefit of IGX‐MS, demonstrated by experiments on GroEL and purified bovine heart mitochondria, is the substantial reduction of undesired over‐length cross‐links compared to in‐solution cross‐linking. We next used IGX‐MS to investigate the complement components C5, C6, and their hetero‐dimeric C5b6 complex. The obtained cross‐links were used to generate a refined structural model of the complement component C6, resembling C6 in its inactivated state. This finding shows that IGX‐MS can provide new insights into the initial stages of the terminal complement pathway.     

SARS‐CoV‐2 nucleocapsid protein phase‐separates with RNA and with human hnRNPs

Tightly packed complexes of nucleocapsid protein and genomic RNA form the core of viruses and assemble within viral factories, dynamic compartments formed within the host cells associated with human stress granules. Here, we test the possibility that the multivalent RNA‐binding nucleocapsid protein (N) from severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2) condenses with RNA via liquid–liquid phase separation (LLPS) and that N protein can be recruited in phase‐separated forms of human RNA‐binding proteins associated with SG formation. Robust LLPS with RNA requires two intrinsically disordered regions (IDRs), the N‐terminal IDR and central‐linker IDR, as well as the folded C‐terminal oligomerization domain, while the folded N‐terminal domain and the C‐terminal IDR are not required. N protein phase separation is induced by addition of non‐specific RNA. In addition, N partitions in vitro into phase‐separated forms of full‐length human hnRNPs (TDP‐43, FUS, hnRNPA2) and their low‐complexity domains (LCs). These results provide a potential mechanism for the role of N in SARS‐CoV‐2 viral genome packing and in host‐protein co‐opting necessary for viral replication and infectivity.     

Co‐occurrence models fail to infer underlying patterns of avoidance and aggregation when closure is violated

Advances in multi‐species monitoring have prompted an increase in the use of multi‐species occupancy analyses to assess patterns of co‐occurrence among species, even when data were collected at scales likely violating the assumption that sites were closed to changes in the occupancy state for the target species. Violating the closure assumption may lead to erroneous conclusions related to patterns of co‐occurrence among species. Occurrence for two hypothetical species was simulated under patterns of avoidance, aggregation, or independence, when the closure assumption was either met or not. Simulated populations were sampled at two levels (N = 250 or 100 sites) and two scales of temporal resolution for surveys. Sample data were analyzed with conditional two‐species occupancy models, and performance was assessed based on the proportion of simulations recovering the true pattern of co‐occurrence. Estimates of occupancy were unbiased when closure was met, but biased when closure violations occurred; bias increased when sample size was small and encounter histories were collapsed to a large‐scale temporal resolution. When closure was met and patterns of avoidance and aggregation were simulated, conditional two‐species models tended to correctly find support for non‐independence, and estimated species interaction factors (SIF) aligned with predicted values. By contrast, when closure was violated, models tended to incorrectly infer a pattern of independence and power to detect simulated patterns of avoidance or aggregation that decreased with smaller sample size. Results suggest that when the closure assumption is violated, co‐occurrence models often fail to detect underlying patterns of avoidance or aggregation, and incorrectly identify a pattern of independence among species, which could have negative consequences for our understanding of species interactions and conservation efforts. Thus, when closure is violated, inferred patterns of independence from multi‐species occupancy should be interpreted cautiously, and evidence of avoidance or aggregation is likely a conservative estimate of true pattern or interaction.     

The atomic portrait of SARS‐CoV‐2 as captured by cryo‐electron microscopy

Transmission electron microscopy has historically been indispensable for virology research, as it offers unique insight into virus function. In the past decade, as cryo‐electron microscopy (cryo‐EM) has matured and become more accessible, we have been able to peer into the structure of viruses at the atomic level and understand how they interact with the host cell, with drugs or with antibodies. Perhaps, there was no time in recent history where cryo‐EM was more needed, as SARS‐CoV‐2 has spread around the globe, causing millions of deaths and almost unquantifiable economic devastation. In this concise review, we aim to mark the most important contributions of cryo‐EM to understanding the structure and function of SARS‐CoV‐2 proteins, from surface spikes to the virus core and from virus‐receptor interactions to antibody binding.     

IL‐23 plays a significant role in the augmentation of particulate matter‐mediated allergic airway inflammation

It has been recently that particulate matter (PM) exposure increases the risk and exacerbation of allergic asthma. However, the underlying mechanisms and factors associated with increased allergic responses remain elusive. We evaluated IL‐23 and IL‐23R (receptor) expression, as well as changes in the asthmatic phenotype in mice administered PM and a low dose of house dust mite (HDM). Next, changes in the phenotype and immune responses were evaluated after intranasal administration of anti‐IL‐23 antibody during co‐exposure to PM and low‐dose HDM. We also performed in vitro experiments to investigate the effect of IL‐23. IL‐23 expression was significantly increased in Epcam+CD45− and CD11c+ cells, while that of IL‐23R was increased in Epcam+CD45− cells only in mice administered PM and low‐dose HDM. Administration of anti‐IL‐23 antibody led to decreased airway hyperresponsiveness, eosinophils, and activation of dendritic cells, reduced populations of Th2 Th17, ILC2, the level of IL‐33 and granulocyte‐macrophage colony‐stimulating factor (GM‐CSF). Inhibition of IL‐23 in PM and low‐dose HDM stimulated airway epithelial cell line resulted in decreased IL‐33, GM‐CSF and affected ILC2 and the activation of BMDCs. PM augmented the phenotypes and immunologic responses of asthma even at low doses of HDM. Interestingly, IL‐23 affected immunological changes in airway epithelial cells.     

Membranome 3.0: Database of single‐pass membrane proteins with AlphaFold models

The Membranome database provides comprehensive structural information on single‐pass (i.e., bitopic) membrane proteins from six evolutionarily distant organisms, including protein–protein interactions, complexes, mutations, experimental structures, and models of transmembrane α‐helical dimers. We present a new version of this database, Membranome 3.0, which was significantly updated by revising the set of 5,758 bitopic proteins and incorporating models generated by AlphaFold 2 in the database. The AlphaFold models were parsed into structural domains located at the different membrane sides, modified to exclude low‐confidence unstructured terminal regions and signal sequences, validated through comparison with available experimental structures, and positioned with respect to membrane boundaries. Membranome 3.0 was re‐developed to facilitate visualization and comparative analysis of multiple 3D structures of proteins that belong to a specified family, complex, biological pathway, or membrane type. New tools for advanced search and analysis of proteins, their interactions, complexes, and mutations were included. The database is freely accessible at https://membranome.org.     

NMR hawk‐eyed view of AlphaFold2 structures

The prediction of the three‐dimensional (3D) structure of proteins from the amino acid sequence made a stunning breakthrough reaching atomic accuracy. Using the neural network‐based method AlphaFold2, 3D structures of almost the entire human proteome have been predicted and made available (https://www.alphafold.ebi.ac.uk). To gain insight into how well AlphaFold2 structures represent the conformation of proteins in solution, I here compare the AlphaFold2 structures of selected small proteins with their 3D structures that were determined by nuclear magnetic resonance (NMR) spectroscopy. Proteins were selected for which the 3D solution structures were determined on the basis of a very large number of distance restraints and residual dipolar couplings and are thus some of the best‐resolved solution structures of proteins to date. The quality of the backbone conformation of the AlphaFold2 structures is assessed by fitting a large set of experimental residual dipolar couplings (RDCs). The analysis shows that experimental RDCs fit extremely well to the AlphaFold2 structures predicted for GB3, DinI, and ubiquitin. In the case of GB3, the accuracy of the AlphaFold2 structure even surpasses that of a 1.1 Å crystal structure. Fitting of experimental RDCs furthermore allows identification of AlphaFold2 structures that are best representative of the protein''s conformation in solution as seen for the EF hands of the N‐terminal domain of Ca²⁺‐ligated calmodulin. Taken together, the analysis shows that structures predicted by AlphaFold2 can be highly representative of the solution conformation of proteins. The combination of AlphaFold2 structures with RDCs promises to be a powerful approach to study structural changes in proteins.     

Post‐association barrier to host switching maintained despite strong selection in a novel mutualism

Following a host shift, repeated co‐passaging of a mutualistic pair is expected to increase fitness over time in one or both species. Without adaptation, a novel association may be evolutionarily short‐lived as it is likely to be outcompeted by native pairings. Here, we test whether experimental evolution can rescue a low‐fitness novel pairing between two sympatric species of Steinernema nematodes and their symbiotic Xenorhabdus bacteria. Despite low mean fitness in the novel association, considerable variation in nematode reproduction was observed across replicate populations. We selected the most productive infections, co‐passaging this novel mutualism nine times to determine whether selection could improve the fitness of either or both partners. We found that neither partner showed increased fitness over time. Our results suggest that the variation in association success was not heritable and that mutational input was insufficient to allow evolution to facilitate this host shift. Thus, post‐association costs of host switching may represent a formidable barrier to novel partnerships among sympatric mutualists.     

Neutralization of SARS‐CoV‐2 by highly potent,hyperthermostable, and mutation‐tolerant nanobodies

Monoclonal anti‐SARS‐CoV‐2 immunoglobulins represent a treatment option for COVID‐19. However, their production in mammalian cells is not scalable to meet the global demand. Single‐domain (VHH) antibodies (also called nanobodies) provide an alternative suitable for microbial production. Using alpaca immune libraries against the receptor‐binding domain (RBD) of the SARS‐CoV‐2 Spike protein, we isolated 45 infection‐blocking VHH antibodies. These include nanobodies that can withstand 95°C. The most effective VHH antibody neutralizes SARS‐CoV‐2 at 17–50 pM concentration (0.2–0.7 µg per liter), binds the open and closed states of the Spike, and shows a tight RBD interaction in the X‐ray and cryo‐EM structures. The best VHH trimers neutralize even at 40 ng per liter. We constructed nanobody tandems and identified nanobody monomers that tolerate the K417N/T, E484K, N501Y, and L452R immune‐escape mutations found in the Alpha, Beta, Gamma, Epsilon, Iota, and Delta/Kappa lineages. We also demonstrate neutralization of the Beta strain at low‐picomolar VHH concentrations. We further discovered VHH antibodies that enforce native folding of the RBD in the E. coli cytosol, where its folding normally fails. Such “fold‐promoting” nanobodies may allow for simplified production of vaccines and their adaptation to viral escape‐mutations.     

Multiple subunit fitting into a low-resolution density map of a macromolecular complex using a gaussian mixture model

Recently, electron microscopy measurement of single particles has enabled us to reconstruct a low-resolution 3D density map of large biomolecular complexes. If structures of the complex subunits can be solved by x-ray crystallography at atomic resolution, fitting these models into the 3D density map can generate an atomic resolution model of the entire large complex. The fitting of multiple subunits, however, generally requires large computational costs; therefore, development of an efficient algorithm is required. We developed a fast fitting program, “gmfit”, which employs a Gaussian mixture model (GMM) to represent approximated shapes of the 3D density map and the atomic models. A GMM is a distribution function composed by adding together several 3D Gaussian density functions. Because our model analytically provides an integral of a product of two distribution functions, it enables us to quickly calculate the fitness of the density map and the atomic models. Using the integral, two types of potential energy function are introduced: the attraction potential energy between a 3D density map and each subunit, and the repulsion potential energy between subunits. The restraint energy for symmetry is also employed to build symmetrical origomeric complexes. To find the optimal configuration of subunits, we randomly generated initial configurations of subunit models, and performed a steepest-descent method using forces and torques of the three potential energies. Comparison between an original density map and its GMM showed that the required number of Gaussian distribution functions for a given accuracy depended on both resolution and molecular size. We then performed test fitting calculations for simulated low-resolution density maps of atomic models of homodimer, trimer, and hexamer, using different search parameters. The results indicated that our method was able to rebuild atomic models of a complex even for maps of 30 Å resolution if sufficient numbers (eight or more) of Gaussian distribution functions were employed for each subunit, and the symmetric restraints were assigned for complexes with more than three subunits. As a more realistic test, we tried to build an atomic model of the GroEL/ES complex by fitting 21-subunit atomic models into the 3D density map obtained by cryoelectron microscopy using the C7 symmetric restraints. A model with low root mean-square deviations (14.7 Å) was obtained as the lowest-energy model, showing that our fitting method was reasonably accurate. Inclusion of other restraints from biological and biochemical experiments could further enhance the accuracy.     

Quantitative contributions of TNF receptor superfamily members to CD8+ T‐cell responses

T‐cell responses to infections and cancers are regulated by co‐signalling receptors grouped into the binary categories of co‐stimulation or co‐inhibition. The co‐stimulation TNF receptor superfamily (TNFRSF) members 4‐1BB, CD27, GITR and OX40 have similar signalling mechanisms raising the question of whether they have similar impacts on T‐cell responses. Here, we screened for the quantitative impact of these TNFRSFs on primary human CD8⁺ T‐cell cytokine production. Although both 4‐1BB and CD27 increased production, only 4‐1BB was able to prolong the duration over which cytokine was produced, and both had only modest effects on antigen sensitivity. An operational model explained these different phenotypes using shared signalling based on the surface expression of 4‐1BB being regulated through signalling feedback. The model predicted and experiments confirmed that CD27 co‐stimulation increases 4‐1BB expression and subsequent 4‐1BB co‐stimulation. GITR and OX40 displayed only minor effects on their own but, like 4‐1BB, CD27 could enhance GITR expression and subsequent GITR co‐stimulation. Thus, different co‐stimulation receptors can have different quantitative effects allowing for synergy and fine‐tuning of T‐cell responses.     

Gel‐like inclusions of C‐terminal fragments of TDP‐43 sequester stalled proteasomes in neurons

Aggregation of the multifunctional RNA‐binding protein TDP‐43 defines large subgroups of amyotrophic lateral sclerosis and frontotemporal dementia and correlates with neurodegeneration in both diseases. In disease, characteristic C‐terminal fragments of ~25 kDa ("TDP‐25") accumulate in cytoplasmic inclusions. Here, we analyze gain‐of‐function mechanisms of TDP‐25 combining cryo‐electron tomography, proteomics, and functional assays. In neurons, cytoplasmic TDP‐25 inclusions are amorphous, and photobleaching experiments reveal gel‐like biophysical properties that are less dynamic than nuclear TDP‐43. Compared with full‐length TDP‐43, the TDP‐25 interactome is depleted of low‐complexity domain proteins. TDP‐25 inclusions are enriched in 26S proteasomes adopting exclusively substrate‐processing conformations, suggesting that inclusions sequester proteasomes, which are largely stalled and no longer undergo the cyclic conformational changes required for proteolytic activity. Reporter assays confirm that TDP‐25 impairs proteostasis, and this inhibitory function is enhanced by ALS‐causing TDP‐43 mutations. These findings support a patho‐physiological relevance of proteasome dysfunction in ALS/FTD.     

To minimize foraging time,use high‐efficiency,energy‐expensive search and capture methods when food is abundant but low‐efficiency,low‐cost methods during food shortages

Based on a mathematical model, I show that the amount of food in the habitat determines which among alternative methods for search of prey, respectively, for pursuit‐and‐capture give the shortest daily foraging time. The higher the locomotor activity, the higher the rate of energy expenditure and the larger the habitat space a predator can search for prey per time unit. Therefore, I assume that the more efficient a foraging method is, the higher its rate of energy expenditure. Survival selection favors individuals that use foraging methods that cover their energy needs in the shortest possible time. Therefore, I take the optimization criterion to be minimization of the daily foraging time or, equivalently, maximization of the rate of net energy gain. When time is limiting and food is in short supply, as during food bottleneck periods, low‐efficiency, low‐cost foraging methods give shorter daily foraging times than high‐efficiency, energy‐expensive foraging methods. When time is limiting, food is abundant and energy needs are large, as during reproduction, high‐efficiency high‐cost foraging methods give shorter daily foraging times than low‐efficiency low‐cost foraging methods. When time is not limiting, food is abundant, and energy needs are small, the choice of foraging method is not critical. Small animals have lower rates of energy expenditure for locomotion than large animals. At a given food density and with similar diet, small animals are therefore more likely than large ones to minimize foraging time by using high‐efficiency energy‐expansive foraging methods and to exploit patches and sites that require energy‐demanding locomotion modes. Survival selection takes place at food shortages, while low‐efficiency low‐cost foraging methods are used, whereas reproduction selection occurs when food is abundant and high‐efficiency energy‐expensive foraging methods do better. In seasonal environments, selection therefore acts on different foraging methods at different times. Morphological adaptation to one method may oppose adaptation to another. Such conflicts select against foraging and morphological specialization and tend to give species‐poor communities of year‐round resident generalists. But a stable year‐round food supply favors specialization, niche narrowing, and dense species packing.     

Context‐dependency in carnivore co‐occurrence across a multi‐use conservation landscape

Carnivore intraguild dynamics depend on a complex interplay of environmental affinities and interspecific interactions. Context‐dependency is commonly expected with varying suites of interacting species and environmental conditions but seldom empirically described. In South Africa, decentralized approaches to conservation and the resulting multi‐tenure conservation landscapes have markedly altered the environmental stage that shapes the structure of local carnivore assemblages. We explored assemblage‐wide patterns of carnivore spatial (residual occupancy probability) and temporal (diel activity overlap) co‐occurrence across three adjacent wildlife‐oriented management contexts—a provincial protected area, a private ecotourism reserve, and commercial game ranches. We found that carnivores were generally distributed independently across space, but existing spatial dependencies were context‐specific. Spatial overlap was most common in the protected area, where species occur at higher relative abundances, and in game ranches, where predator persecution presumably narrows the scope for spatial asymmetries. In the private reserve, spatial co‐occurrence patterns were more heterogeneous but did not follow a dominance hierarchy associated with higher apex predator densities. Pair‐specific variability suggests that subordinate carnivores may alternate between pre‐emptive behavioral strategies and fine‐scale co‐occurrence with dominant competitors. Consistency in species‐pairs diel activity asynchrony suggested that temporal overlap patterns in our study areas mostly depend on species'' endogenous clock rather than the local context. Collectively, our research highlights the complexity and context‐dependency of guild‐level implications of current management and conservation paradigms; specifically, the unheeded potential for interventions to influence the local network of carnivore interactions with unknown population‐level and cascading effects.     

The STI1‐domain is a flexible alpha‐helical fold with a hydrophobic groove

STI1‐domains are present in a variety of co‐chaperone proteins and are required for the transfer of hydrophobic clients in various cellular processes. The domains were first identified in the yeast Sti1 protein where they were referred to as DP1 and DP2. Based on hidden Markov model searches, this domain had previously been found in other proteins including the mammalian co‐chaperone SGTA, the DNA damage response protein Rad23, and the chloroplast import protein Tic40. Here, we refine the domain definition and carry out structure‐based sequence alignment of STI1‐domains showing conservation of five amphipathic helices. Upon examinations of these identified domains, we identify a preceding helix 0 and unifying sequence properties, determine new molecular models, and recognize that STI1‐domains nearly always occur in pairs. The similarity at the sequence, structure, and molecular levels likely supports a unified functional role.     

Multi-resolution contour-based fitting of macromolecular structures

A novel contour-based matching criterion is presented for the quantitative docking of high-resolution structures of components into low-resolution maps of macromolecular complexes. The proposed Laplacian filter is combined with a six-dimensional search using fast Fourier transforms to rapidly scan the rigid-body degrees of freedom of a probe molecule relative to a fixed target density map. A comparison of the docking performance with the standard cross-correlation criterion demonstrates that contour matching with the Laplacian filter significantly extends the viable resolution range of correlation-based fitting to resolutions as low as 30 A. The gain in docking precision at medium to low resolution (15-30 A) is critical for image reconstructions from electron microscopy (EM). The new algorithm enables for the first time the reliable docking of smaller molecular components into EM densities of large biomolecular assemblies at such low resolutions. As an example of the practical effectiveness of contour-based fitting, a new pseudo-atomic model of a microtubule was constructed from a 20 A resolution EM map and from atomic structures of alpha and beta tubulin subunits.