ESCPE-1 mediates retrograde endosomal sorting of the SARS-CoV-2 host factor Neuropilin-1.

Endosomal sorting maintains cellular homeostasis by recycling transmembrane proteins and associated proteins and lipids (termed "cargoes") from the endosomal network to multiple subcellular destinations, including retrograde traffic to the trans-Golgi network (TGN). Viral and bacterial pathogens subvert retrograde trafficking machinery to facilitate infectivity. Here, we develop a proteomic screen to identify retrograde cargo proteins of the endosomal SNX-BAR sorting complex promoting exit 1 (ESCPE-1). Using this methodology, we identify Neuropilin-1 (NRP1), a recently characterized host factor for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, as a cargo directly bound and trafficked by ESCPE-1. ESCPE-1 mediates retrograde trafficking of engineered nanoparticles functionalized with the NRP1-interacting peptide of the SARS-CoV-2 spike (S) protein. CRISPR-Cas9 deletion of ESCPE-1 subunits reduces SARS-CoV-2 infection levels in cell culture. ESCPE-1 sorting of NRP1 may therefore play a role in the intracellular membrane trafficking of NRP1-interacting viruses such as SARS-CoV-2.

A potential interaction between the SARS-CoV-2 spike protein and nicotinic acetylcholine receptors.

Changeux et al. (Changeux et al. C. R. Biol. 343:33-39.) recently suggested that the SARS-CoV-2 spike protein may interact with nicotinic acetylcholine receptors (nAChRs) and that such interactions may be involved in pathology and infectivity. This hypothesis is based on the fact that the SARS-CoV-2 spike protein contains a sequence motif similar to known nAChR antagonists. Here, we use molecular simulations of validated atomically detailed structures of nAChRs and of the spike to investigate the possible binding of the Y674-R685 region of the spike to nAChRs. We examine the binding of the Y674-R685 loop to three nAChRs, namely the human α4ß2 and α7 subtypes and the muscle-like αßγδ receptor from Tetronarce californica. Our results predict that Y674-R685 has affinity for nAChRs. The region of the spike responsible for binding contains a PRRA motif, a four-residue insertion not found in other SARS-like coronaviruses. The conformational behavior of the bound Y674-R685 is highly dependent on the receptor subtype; it adopts extended conformations in the α4ß2 and α7 complexes but is more compact when bound to the muscle-like receptor. In the α4ß2 and αßγδ complexes, the interaction of Y674-R685 with the receptors forces the loop C region to adopt an open conformation, similar to other known nAChR antagonists. In contrast, in the α7 complex, Y674-R685 penetrates deeply into the binding pocket in which it forms interactions with the residues lining the aromatic box, namely with TrpB, TyrC1, and TyrC2. Estimates of binding energy suggest that Y674-R685 forms stable complexes with all three nAChR subtypes. Analyses of simulations of the glycosylated spike show that the Y674-R685 region is accessible for binding. We suggest a potential binding orientation of the spike protein with nAChRs, in which they are in a nonparallel arrangement to one another.

De Novo Designed Peptide and Protein Hairpins Self-Assemble into Sheets and Nanoparticles.

The design and assembly of peptide-based materials has advanced considerably, leading to a variety of fibrous, sheet, and nanoparticle structures. A remaining challenge is to account for and control different possible supramolecular outcomes accessible to the same or similar peptide building blocks. Here a de novo peptide system is presented that forms nanoparticles or sheets depending on the strategic placement of a "disulfide pin" between two elements of secondary structure that drive self-assembly. Specifically, homodimerizing and homotrimerizing de novo coiled-coil α-helices are joined with a flexible linker to generate a series of linear peptides. The helices are pinned back-to-back, constraining them as hairpins by a disulfide bond placed either proximal or distal to the linker. Computational modeling indicates, and advanced microscopy shows, that the proximally pinned hairpins self-assemble into nanoparticles, whereas the distally pinned constructs form sheets. These peptides can be made synthetically or recombinantly to allow both chemical modifications and the introduction of whole protein cargoes as required.

Molecular Simulations suggest Vitamins, Retinoids and Steroids as Ligands of the Free Fatty Acid Pocket of the SARS-CoV-2 Spike Protein*.

We investigate binding of linoleate and other potential ligands to the recently discovered fatty acid binding site in the SARS-CoV-2 spike protein, using docking and molecular dynamics simulations. Simulations suggest that linoleate and dexamethasone stabilize the locked spike conformation, thus reducing the opportunity for ACE2 interaction. In contrast, cholesterol may expose the receptor-binding domain by destabilizing the closed structure, preferentially binding to a different site in the hinge region of the open structure. We docked a library of FDA-approved drugs to the fatty acid site using an approach that reproduces the structure of the linoleate complex. Docking identifies steroids (including dexamethasone and vitamin D); retinoids (some known to be active in vitro, and vitamin A); and vitamin K as potential ligands that may stabilize the closed conformation. The SARS-CoV-2 spike fatty acid site may bind a diverse array of ligands, including dietary components, and therefore provides a promising target for therapeutics or prophylaxis.

Emergence of a Thrombospondin Superfamily at the Origin of Metazoans.

Extracellular matrix (ECM) is considered central to the evolution of metazoan multicellularity; however, the repertoire of ECM proteins in nonbilaterians remains unclear. Thrombospondins (TSPs) are known to be well conserved from cnidarians to vertebrates, yet to date have been considered a unique family, principally studied for matricellular functions in vertebrates. Through searches utilizing the highly conserved C-terminal region of TSPs, we identify undisclosed new families of TSP-related proteins in metazoans, designated mega-TSP, sushi-TSP, and poriferan-TSP, each with a distinctive phylogenetic distribution. These proteins share the TSP C-terminal region domain architecture, as determined by domain composition and analysis of molecular models against known structures. Mega-TSPs, the only form identified in ctenophores, are typically >2,700 aa and are also characterized by N-terminal leucine-rich repeats and central cadherin/immunoglobulin domains. In cnidarians, which have a well-defined ECM, Mega-TSP was expressed throughout embryogenesis in Nematostella vectensis, with dynamic endodermal expression in larvae and primary polyps and widespread ectodermal expression in adult Nematostella vectensis and Hydra magnipapillata polyps. Hydra Mega-TSP was also expressed during regeneration and siRNA-silencing of Mega-TSP in Hydra caused specific blockade of head regeneration. Molecular phylogenetic analyses based on the conserved TSP C-terminal region identified each of the TSP-related groups to form clades distinct from the canonical TSPs. We discuss models for the evolution of the newly defined TSP superfamily by gene duplications, radiation, and gene losses from a debut in the last metazoan common ancestor. Together, the data provide new insight into the evolution of ECM and tissue organization in metazoans.

Beyond icosahedral symmetry in packings of proteins in spherical shells.

The formation of quasi-spherical cages from protein building blocks is a remarkable self-assembly process in many natural systems, where a small number of elementary building blocks are assembled to build a highly symmetric icosahedral cage. In turn, this has inspired synthetic biologists to design de novo protein cages. We use simple models, on multiple scales, to investigate the self-assembly of a spherical cage, focusing on the regularity of the packing of protein-like objects on the surface. Using building blocks, which are able to pack with icosahedral symmetry, we examine how stable these highly symmetric structures are to perturbations that may arise from the interplay between flexibility of the interacting blocks and entropic effects. We find that, in the presence of those perturbations, icosahedral packing is not the most stable arrangement for a wide range of parameters; rather disordered structures are found to be the most stable. Our results suggest that (i) many designed, or even natural, protein cages may not be regular in the presence of those perturbations and (ii) optimizing those flexibilities can be a possible design strategy to obtain regular synthetic cages with full control over their surface properties.

Intraring allostery controls the function and assembly of a hetero-oligomeric class II chaperonin.

Class II chaperonins are essential multisubunit complexes that aid the folding of nonnative proteins in the cytosol of archaea and eukarya. They use energy derived from ATP to drive a series of structural rearrangements that enable polypeptides to fold within their central cavity. These events are regulated by an elaborate allosteric mechanism in need of elucidation. We employed mutagenesis and experimental analysis in concert with in silico molecular dynamics simulations and interface-binding energy calculations to investigate the class II chaperonin from Thermoplasma acidophilum. Here we describe the effects on the asymmetric allosteric mechanism and on hetero-oligomeric complex formation in a panel of mutants in the ATP-binding pocket of the α and ß subunits. Our observations reveal a potential model for a nonconcerted folding mechanism optimized for protecting and refolding a range of nonnative substrates under different environmental conditions, starting to unravel the role of subunit heterogeneity in this folding machine and establishing important links with the behavior of the most complex eukaryotic chaperonins.-Shoemark, D. K., Sessions, R. B., Brancaccio, A., Bigotti, M. G. Intraring allostery controls the function and assembly of a hetero-oligomeric class II chaperonin.

The dynamical interplay between a megadalton peptide nanocage and solutes probed by microsecond atomistic MD; implications for design.

Understanding the assembly and dynamics of protein-based supramolecular capsids and cages is of fundamental importance and could lead to applications in synthetic biology and biotechnology. Here we present long and large atomistic molecular dynamics simulations of de novo designed self-assembling protein nanocages (SAGEs) in aqueous media. Microsecond simulations, comprised of ≈42 million atoms for three pre-formed SAGEs of different charges, in the presence of solutes and solvent have been completed. Here, the dynamics, stability and porosity of the peptide networks are explored along with their interactions with ions, small molecules and macromolecular solutes. All assemblies are stable over the µs timescale, and the solutes show a mixture of transport behaviour across or adherence to the fabric of the SAGE particles. Solute proteins largely retained native-like conformation on contact with SAGE. Certain residues of the SAGE peptides are identified as "repeat offenders" for contacting many different solutes, which suggest modifications to reduce non-specific binding. These studies highlight how molecular dynamics can aid the design process of SAGE and similar assemblies for potential applications as diverse as platforms for drug and vaccine delivery and nanoreactors to encapsulate enzyme pathways.

Improved Specificity of Glutamate Decarboxylase 65 Autoantibody Measurement Using Luciferase-Based Immunoprecipitation System Assays.

Autoantibodies to glutamate decarboxylase (GADA) are widely used in the prediction and classification of type 1 diabetes. GADA radiobinding assays (RBAs) using N-terminally truncated antigens offer improved specificity, but radioisotopes limit the high-throughput potential for population screening. Luciferase-based immunoprecipitation system (LIPS) assays are sensitive and specific alternatives to RBAs with the potential to improve risk stratification. The performance of assays using the Nanoluc luciferase (Nluc)-conjugated GAD65 constructs, Nluc-GAD65(96-585) and full length Nluc-GAD65(1-585), were evaluated in 434 well-characterized serum samples from patients with recent-onset type 1 diabetes and first-degree relatives. Nonradioactive, high-throughput LIPS assays are quicker and require less serum than RBAs. Of 171 relatives previously tested single autoantibody positive for autoantibodies to full-length GAD65 by RBA but had not progressed to diabetes, fewer retested positive by LIPS using either truncated (n = 72) or full-length (n = 111) antigen. The Nluc-GAD65(96-585) truncation demonstrated the highest specificity in LIPS assays overall, but in contrast to RBA, N-terminus truncations did not result in a significant increase in disease-specificity compared with the full-length antigen. This suggests that binding of nonspecific antibodies is affected by the conformational changes resulting from addition of the Nluc antigen. Nluc-GAD65(96-585) LIPS assays offer low-blood-volume, high-specificity GADA tests for screening and diagnostics.

Identification of novel small molecule TGF-ß antagonists using structure-based drug design.

Aberrant transforming growth factor-ß (TGF-ß) signalling has been associated with a number of disease pathologies, such as the development of fibrosis in the heart, lung and liver, cardiovascular disease and cancer, hence the TGF-ß pathway represents a promising target for a variety of diseases. However, highly specific ways to inhibit TGF-ß signalling need to be developed to prevent cross-talk with related receptors and minimise unwanted side effects. We have used used virtual screening and molecular docking to identify small molecule inhibitors of TGF-ß binding to TßRII. The crystal structure of TGF-ß3 in complex with the extracellular domain of the type II TGF-ß receptor was taken as a starting point for molecular docking and we developed a structure-based pharmacophore model to identify compounds that competitively inhibit the binding of TGF-ß to TßRII and antogonize TGF-ß signalling. We have experimentally tested 67 molecules suggested by in silico screening and similarity searching for their ability to inhibit TGF-ß signalling in TGF-ß-dependent luciferase assays in vitro and the molecule with the strongest inhibition had an IC50 of 18 µM. These compounds were selected to bind to the SS1 subsite (composed of F30, C31, D32, I50, T51 S52, I53, C54 and E55) of TßRII and all share the general property of being aromatic and fairly flat. Molecular dynamics simulations confirmed that this was the most likely binding mode. The computational methods used and the hits identified in this study provide an excellent guide to medicinal chemistry efforts to design tighter binding molecules to disrupt the TGF-ß/TßRII interaction.

SARS-CoV-2 spike variants differ in their allosteric responses to linoleic acid.

The SARS-CoV-2 spike protein contains a functionally important fatty acid (FA) binding site, which is also found in some other coronaviruses, e.g. SARS-CoV and MERS-CoV. The occupancy of the FA site by linoleic acid (LA) reduces infectivity by 'locking' the spike in a less infectious conformation. Here, we use dynamical-nonequilibrium molecular dynamics (D-NEMD) simulations to compare the allosteric responses of spike variants to LA removal. D-NEMD simulations show that the FA site is coupled to other functional regions of the protein, e.g. the receptor-binding motif (RBM), N-terminal domain (NTD), furin cleavage site, and regions surrounding the fusion peptide. D-NEMD simulations also identify the allosteric networks connecting the FA site to these functional regions. The comparison between the wild-type spike and four variants (Alpha, Delta, Delta plus, and Omicron BA.1) shows that the variants differ significantly in their responses to LA removal. The allosteric connections to the FA site on Alpha are generally similar to those on the wild-type protein, with the exception of the RBM and the S71-R78 region, which show a weaker link to the FA site. In contrast, Omicron is the most different variant, exhibiting significant differences in the RBM, NTD, V622-L629, and furin cleavage site. These differences in the allosteric modulation may be of functional relevance, potentially affecting transmissibility and virulence. Experimental comparison of the effects of LA on SARS-CoV-2 variants, including emerging variants, is warranted.

ATG8-dependent LMX1B-autophagy crosstalk shapes human midbrain dopaminergic neuronal resilience.

The LIM homeodomain transcription factors LMX1A and LMX1B are essential mediators of midbrain dopaminergic neuronal (mDAN) differentiation and survival. Here we show that LMX1A and LMX1B are autophagy transcription factors that provide cellular stress protection. Their suppression dampens the autophagy response, lowers mitochondrial respiration, and elevates mitochondrial ROS, and their inducible overexpression protects against rotenone toxicity in human iPSC-derived mDANs in vitro. Significantly, we show that LMX1A and LMX1B stability is in part regulated by autophagy, and that these transcription factors bind to multiple ATG8 proteins. Binding is dependent on subcellular localization and nutrient status, with LMX1B interacting with LC3B in the nucleus under basal conditions and associating with both cytosolic and nuclear LC3B during nutrient starvation. Crucially, ATG8 binding stimulates LMX1B-mediated transcription for efficient autophagy and cell stress protection, thereby establishing a novel LMX1B-autophagy regulatory axis that contributes to mDAN maintenance and survival in the adult brain.

Ca 2+ Regulates Dimerization of the BAR Domain Protein PICK1 and Consequent Membrane Curvature.

Bin-Amphiphysin-Rvs (BAR) domain proteins are critical regulators of membrane geometry. They induce and stabilize membrane curvature for processes, such as clathrin-coated pit formation and endosomal membrane tubulation. BAR domains form their characteristic crescent-shaped structure in the dimeric form, indicating that the formation of the dimer is critical to their function of inducing membrane curvature and suggesting that a dynamic monomer-dimer equilibrium regulated by cellular signaling would be a powerful mechanism for controlling BAR domain protein function. However, to the best of our knowledge, cellular mechanisms for regulating BAR domain dimerization remain unexplored. PICK1 is a Ca2+-binding BAR domain protein involved in the endocytosis and endosomal recycling of neuronal AMPA receptors and other transmembrane proteins. In this study, we demonstrated that PICK1 dimerization is regulated by a direct effect of Ca2+ ions via acidic regions in the BAR domain and at the N-terminus. While the cellular membrane tubulating activity of PICK1 is absent under basal conditions, Ca2+ influx causes the generation of membrane tubules that originate from the cell surface. Furthermore, in neurons, PICK1 dimerization increases transiently following NMDA receptor stimulation. We believe that this novel mechanism for regulating BAR domain dimerization and function represents a significant conceptual advance in our knowledge about the regulation of cellular membrane curvature.

Identification and validation of novel microtubule suppressors with an imidazopyridine scaffold through structure-based virtual screening and docking.

Targeting the colchicine binding site of α/ß tubulin microtubules can lead to suppression of microtubule dynamics, cell cycle arrest and apoptosis. Therefore, the development of microtubule (MT) inhibitors is considered a promising route to anticancer agents. Our approach to identify novel scaffolds as MT inhibitors depends on a 3D-structure-based pharmacophore approach and docking using three programs MOE, Autodock and BUDE (Bristol University Docking Engine) to screen a library of virtual compounds. From this work we identified the compound 7-(3-hydroxy-4-methoxy-phenyl)-3-(3-trifluoromethyl-phenyl)-6,7-dihydro-3H-imidazo[4,5-b]pyridin-5-ol (6) as a novel inhibitor scaffold. This compound inhibited several types of cancer cell proliferation at low micromolar concentrations with low toxicity. Compound 6 caused cell cycle arrest in the G2/M phase and blocked tubulin polymerization at low micromolar concentration (IC50 = 6.1 ±0.1 µM), inducing apoptosis via activation of caspase 9, increasing the level of the pro-apoptotic protein Bax and decreasing the level of the anti-apoptotic protein Bcl2. In summary, our approach identified a lead compound with potential antimitotic and antiproliferative activity.

The fatty acid site is coupled to functional motifs in the SARS-CoV-2 spike protein and modulates spike allosteric behaviour.

The SARS-CoV-2 spike protein is the first contact point between the SARS-CoV-2 virus and host cells and mediates membrane fusion. Recently, a fatty acid binding site was identified in the spike (Toelzer et al. Science 2020). The presence of linoleic acid at this site modulates binding of the spike to the human ACE2 receptor, stabilizing a locked conformation of the protein. Here, dynamical-nonequilibrium molecular dynamics simulations reveal that this fatty acid site is coupled to functionally relevant regions of the spike, some of them far from the fatty acid binding pocket. Removal of a ligand from the fatty acid binding site significantly affects the dynamics of distant, functionally important regions of the spike, including the receptor-binding motif, furin cleavage site and fusion-peptide-adjacent regions. Simulations of the D614G mutant show differences in behaviour between these clinical variants of the spike: the D614G mutant shows a significantly different conformational response for some structural motifs relevant for binding and fusion. The simulations identify structural networks through which changes at the fatty acid binding site are transmitted within the protein. These communication networks significantly involve positions that are prone to mutation, indicating that observed genetic variation in the spike may alter its response to linoleate binding and associated allosteric communication.

Structural insights in cell-type specific evolution of intra-host diversity by SARS-CoV-2.

As the global burden of SARS-CoV-2 infections escalates, so does the evolution of viral variants with increased transmissibility and pathology. In addition to this entrenched diversity, RNA viruses can also display genetic diversity within single infected hosts with co-existing viral variants evolving differently in distinct cell types. The BriSΔ variant, originally identified as a viral subpopulation from SARS-CoV-2 isolate hCoV-19/England/02/2020, comprises in the spike an eight amino-acid deletion encompassing a furin recognition motif and S1/S2 cleavage site. We elucidate the structure, function and molecular dynamics of this spike providing mechanistic insight into how the deletion correlates to viral cell tropism, ACE2 receptor binding and infectivity of this SARS-CoV-2 variant. Our results reveal long-range allosteric communication between functional domains that differ in the wild-type and the deletion variant and support a view of SARS-CoV-2 probing multiple evolutionary trajectories in distinct cell types within the same infected host.

Identification of ß-strand mediated protein-protein interaction inhibitors using ligand-directed fragment ligation.

ß-Strand mediated protein-protein interactions (PPIs) represent underexploited targets for chemical probe development despite representing a significant proportion of known and therapeutically relevant PPI targets. ß-Strand mimicry is challenging given that both amino acid side-chains and backbone hydrogen-bonds are typically required for molecular recognition, yet these are oriented along perpendicular vectors. This paper describes an alternative approach, using GKAP/SHANK1 PDZ as a model and dynamic ligation screening to identify small-molecule replacements for tranches of peptide sequence. A peptide truncation of GKAP functionalized at the N- and C-termini with acylhydrazone groups was used as an anchor. Reversible acylhydrazone bond exchange with a library of aldehyde fragments in the presence of the protein as template and in situ screening using a fluorescence anisotropy (FA) assay identified peptide hybrid hits with comparable affinity to the GKAP peptide binding sequence. Identified hits were validated using FA, ITC, NMR and X-ray crystallography to confirm selective inhibition of the target PDZ-mediated PPI and mode of binding. These analyses together with molecular dynamics simulations demonstrated the ligands make transient interactions with an unoccupied basic patch through electrostatic interactions, establishing proof-of-concept that this unbiased approach to ligand discovery represents a powerful addition to the armory of tools that can be used to identify PPI modulators.

A conserved arginine with non-conserved function is a key determinant of agonist selectivity in α7 nicotinic ACh receptors.

BACKGROUND AND PURPOSE: The α7 and α4ß2* ("*" denotes possibly assembly with another subunit) nicotinic acetylcholine receptors (nAChRs) are the most abundant nAChRs in the mammalian brain. These receptors are the most targeted nAChRs in drug discovery programmes for brain disorders. However, the development of subtype-specific agonists remains challenging due to the high degree of sequence homology and conservation of function in nAChRs. We have developed C(10) variants of cytisine, a partial agonist of α4ß2 nAChR that has been used for smoking cessation. The C(10) methyl analogue used in this study displays negligible affinity for α7 nAChR, while retaining high affinity for α4ß2 nAChR. EXPERIMENTAL APPROACH: The structural underpinning of the selectivity of 10-methylcytisine for α7 and α4ß2 nAChRs was investigated using molecular dynamic simulations, mutagenesis and whole-cell and single-channel current recordings. KEY RESULTS: We identified a conserved arginine in the ß3 strand that exhibits a non-conserved function in nAChRs. In α4ß2 nAChR, the arginine forms a salt bridge with an aspartate residue in loop B that is necessary for receptor expression, whereas in α7 nAChR, this residue is not stabilised by electrostatic interactions, making its side chain highly mobile. This lack of constrain produces steric clashes with agonists and affects the dynamics of residues involved in agonist binding and the coupling network. CONCLUSION AND IMPLICATIONS: We conclude that the high mobility of the ß3-strand arginine in the α7 nAChR influences agonist binding and possibly gating network and desensitisation. The findings have implications for rational design of subtype-selective nAChR agents.

Query-guided protein-protein interaction inhibitor discovery.

Protein-protein interactions (PPIs) are central to biological mechanisms, and can serve as compelling targets for drug discovery. Yet, the discovery of small molecule inhibitors of PPIs remains challenging given the large and typically shallow topography of the interacting protein surfaces. Here, we describe a general approach to the discovery of orthosteric PPI inhibitors that mimic specific secondary protein structures. Initially, hot residues at protein-protein interfaces are identified in silico or from experimental data, and incorporated into secondary structure-based queries. Virtual libraries of small molecules are then shape-matched against the queries, and promising ligands docked to target proteins. The approach is exemplified experimentally using two unrelated PPIs that are mediated by an α-helix (p53/hDM2) and a ß-strand (GKAP/SHANK1-PDZ). In each case, selective PPI inhibitors are discovered with low µM activity as determined by a combination of fluorescence anisotropy and 1H-15N HSQC experiments. In addition, hit expansion yields a series of PPI inhibitors with defined structure-activity relationships. It is envisaged that the generality of the approach will enable discovery of inhibitors of a wide range of unrelated secondary structure-mediated PPIs.