Biological function relies on the fact that biomolecules can switch between different conformations and aggregation states. Such transitions involve a rearrangement of parts of the biomolecules involved that act as dynamic domains. The reliable identification of such domains is thus a key problem in biophysics. In this work we present a method to identify semi-rigid domains based on dynamical data that can be obtained from molecular dynamics simulations or experiments. To this end the average inter-atomic distance-deviations are computed. The resulting matrix is then clustered by a constrained quadratic optimization problem. The reliability and performance of the method are demonstrated for two artificial peptides. Furthermore we correlate the mechanical properties with biological malfunction in three variants of amyloidogenic transthyretin protein, where the method reveals that a pathological mutation destabilizes the natural dimer structure of the protein. Finally the method is used to identify functional domains of the GroEL-GroES chaperone, thus illustrating the efficiency of the method for large biomolecular machines.

Synaptic vesicle fusion is mediated by SNARE proteins forming in between synaptic vesicle (v-SNARE) and plasma membrane (t-SNARE), one of which is Syntaxin-1A. Although exocytosis mainly occurs at active zones, Syntaxin-1A appears to cover the entire neuronal membrane. By using STED super-resolution light microscopy and image analysis of Drosophila neuro-muscular junctions, we show that Syntaxin-1A clusters are more abundant and have an increased size at active zones. A computational particle-based model of syntaxin cluster formation and dynamics is developed. The model is parametrized to reproduce Syntaxin cluster-size distributions found by STED analysis, and successfully reproduces existing FRAP results. The model shows that the neuronal membrane is adjusted in a way to strike a balance between having most syntaxins stored in large clusters, while still keeping a mobile fraction of syntaxins free or in small clusters that can efficiently search the membrane or be traded between clusters. This balance is subtle and can be shifted toward almost no clustering and almost complete clustering by modifying the syntaxin interaction energy on the order of only 1 kBT. This capability appears to be exploited at active zones. The larger active-zone syntaxin clusters are more stable and provide regions of high docking and fusion capability, whereas the smaller clusters outside may serve as flexible reserve pool or sites of spontaneous ectopic release.

Studying atherosclerotic calcification in vivo requires mouse models with genetic modifications. Previous studies showed that injection of recombinant adeno-associated virus vector (AAV) encoding a gain-of-function mutant PCSK9 into mice promotes atherosclerosis. We aimed to study cardiovascular calcification induced by PCSK9 AAV in C57BL/6J mice.

Clathrin-mediated endocytosis (CME) involves membrane-associated scaffolds of the bin-amphiphysin-rvs (BAR) domain protein family as well as the GTPase dynamin, and is accompanied and perhaps triggered by changes in local lipid composition. How protein recruitment, scaffold assembly and membrane deformation is spatiotemporally controlled and coupled to fission is poorly understood. We show by computational modelling and super-resolution imaging that phosphatidylinositol 3,4-bisphosphate [PI(3,4)P2] synthesis within the clathrin-coated area of endocytic intermediates triggers selective recruitment of the PX-BAR domain protein SNX9, as a result of complex interactions of endocytic proteins competing for phospholipids. The specific architecture induces positioning of SNX9 at the invagination neck where its self-assembly regulates membrane constriction, thereby providing a template for dynamin fission. These data explain how lipid conversion at endocytic pits couples local membrane constriction to fission. Our work demonstrates how computational modelling and super-resolution imaging can be combined to unravel function and mechanisms of complex cellular processes.

Recombinant adeno-associated virus (rAAV) vectors have emerged as the safe vehicles of choice for long-term gene transfer in mammalian nervous system. Recombinant adeno-associated virus-mediated localized gene transfer in adult nervous system following direct inoculation, that is, intracerebral or intrathecal, is well documented. However, recombinant adeno-associated virus delivery in defined neuronal populations in adult animals using less-invasive methods as well as avoiding ectopic gene expression following systemic inoculation remain challenging. Harnessing the capability of some recombinant adeno-associated virus serotypes for retrograde transduction may potentially address such limitations (Note: The term retrograde transduction in this manuscript refers to the uptake of injected recombinant adeno-associated virus particles at nerve terminals, retrograde transport, and subsequent transduction of nerve cell soma). In some studies, recombinant adeno-associated virus serotypes 2/6, 2/8, and 2/9 have been shown to exhibit transduction of connected neuroanatomical tracts in adult animals following lower limb intramuscular recombinant adeno-associated virus delivery in a pattern suggestive of retrograde transduction. However, an extensive side-by-side comparison of these serotypes following intramuscular delivery regarding tissue viral load, and the effect of promoter on transgene expression, has not been performed. Hence, we delivered recombinant adeno-associated virus serotypes 2/6, 2/8, or 2/9 encoding enhanced green fluorescent protein (eGFP), under the control of either cytomegalovirus (CMV) or human synapsin (hSyn) promoter, via a single unilateral hindlimb intramuscular injection in the bicep femoris of adult C57BL/6J mice. Four weeks post injection, we quantified viral load and transgene (enhanced green fluorescent protein) expression in muscle and related nervous tissues. Our data show that the select recombinant adeno-associated virus serotypes transduce sciatic nerve and groups of neurons in the dorsal root ganglia on the injected side, indicating that the intramuscular recombinant adeno-associated virus delivery is useful for achieving gene transfer in local neuroanatomical tracts. We also observed sparse recombinant adeno-associated virus viral delivery or eGFP transduction in lumbar spinal cord and a noticeable lack thereof in brain. Therefore, further improvements in recombinant adeno-associated virus design are warranted to achieve efficient widespread retrograde transduction following intramuscular and possibly other peripheral routes of delivery.

Antivirals are needed to combat the COVID-19 pandemic, which is caused by SARS-CoV-2. The clinically-proven protease inhibitor Camostat mesylate inhibits SARS-CoV-2 infection by blocking the virus-activating host cell protease TMPRSS2. However, antiviral activity of Camostat mesylate metabolites and potential viral resistance have not been analyzed. Moreover, antiviral activity of Camostat mesylate in human lung tissue remains to be demonstrated.

Host proteases have been suggested to be crucial for dissemination of MERS, SARS-CoV, and SARS-CoV-2 coronaviruses, but the relative contribution of membrane versus intracellular proteases remains controversial. Transmembrane serine protease 2 (TMPRSS2) is regarded as one of the main proteases implicated in the coronavirus S protein priming, an important step for binding of the S protein to the angiotensin-converting enzyme 2 (ACE2) receptor before cell entry.

SORL1 is strongly associated with both sporadic and familial forms of Alzheimer's disease (AD), but a lack of information about alternatively spliced transcripts currently limits our understanding of the role of SORL1 in AD. Here, we describe a SORL1 transcript (SORL1-38b) characterized by inclusion of a novel exon (E38b) that encodes a truncated protein. We identified E38b-containing transcripts in several brain regions, with the highest expression in the cerebellum and showed that SORL1-38b is largely located in neuronal dendrites, which is in contrast to the somatic distribution of transcripts encoding the full-length SORLA protein (SORL1-fl). SORL1-38b transcript levels were significantly reduced in AD cerebellum in three independent cohorts of postmortem brains, whereas no changes were observed for SORL1-fl. A trend of lower 38b transcript level in cerebellum was found for individuals carrying the risk variant at rs2282649 (known as SNP24), although not reaching statistical significance. These findings suggest synaptic functions for SORL1-38b in the brain, uncovering novel aspects of SORL1 that can be further explored in AD research.

We introduce the software package ReaDDy for simulation of detailed spatiotemporal mechanisms of dynamical processes in the cell, based on reaction-diffusion dynamics with particle resolution. In contrast to other particle-based reaction kinetics programs, ReaDDy supports particle interaction potentials. This permits effects such as space exclusion, molecular crowding and aggregation to be modeled. The biomolecules simulated can be represented as a sphere, or as a more complex geometry such as a domain structure or polymer chain. ReaDDy bridges the gap between small-scale but highly detailed molecular dynamics or Brownian dynamics simulations and large-scale but little-detailed reaction kinetics simulations. ReaDDy has a modular design that enables the exchange of the computing core by efficient platform-specific implementations or dynamical models that are different from Brownian dynamics.

The ADAMTS9 rs4607103 C allele is one of the few gene variants proposed to increase the risk of type 2 diabetes through an impairment of insulin sensitivity. We show that the variant is associated with increased expression of the secreted ADAMTS9 and decreased insulin sensitivity and signaling in human skeletal muscle. In line with this, mice lacking Adamts9 selectively in skeletal muscle have improved insulin sensitivity. The molecular link between ADAMTS9 and insulin signaling was characterized further in a model where ADAMTS9 was overexpressed in skeletal muscle. This selective overexpression resulted in decreased insulin signaling presumably mediated through alterations of the integrin β1 signaling pathway and disruption of the intracellular cytoskeletal organization. Furthermore, this led to impaired mitochondrial function in mouse muscle-an observation found to be of translational character because humans carrying the ADAMTS9 risk allele have decreased expression of mitochondrial markers. Finally, we found that the link between ADAMTS9 overexpression and impaired insulin signaling could be due to accumulation of harmful lipid intermediates. Our findings contribute to the understanding of the molecular mechanisms underlying insulin resistance and type 2 diabetes and point to inhibition of ADAMTS9 as a potential novel mode of treating insulin resistance.

FGF2 is exported from cells by an unconventional secretory mechanism. Here, we directly visualized individual FGF2 membrane translocation events at the plasma membrane using live cell TIRF microscopy. This process was dependent on both PI(4,5)P2-mediated recruitment of FGF2 at the inner leaflet and heparan sulfates capturing FGF2 at the outer plasma membrane leaflet. By simultaneous imaging of both FGF2 membrane recruitment and the appearance of FGF2 at the cell surface, we revealed the kinetics of FGF2 membrane translocation in living cells with an average duration of ∼200 ms. Furthermore, we directly demonstrated FGF2 oligomers at the inner leaflet of living cells with a FGF2 dimer being the most prominent species. We propose this dimer to represent a key intermediate in the formation of higher FGF2 oligomers that form membrane pores and put forward a kinetic model explaining the mechanism by which membrane-inserted FGF2 oligomers serve as dynamic translocation intermediates during unconventional secretion of FGF2.

The visual pigment rhodopsin belongs to the family of G protein-coupled receptors that can form higher oligomers. It is controversial whether rhodopsin forms oligomers and whether oligomers are functionally relevant. Here, we study rhodopsin organization in cryosections of dark-adapted mouse rod photoreceptors by cryoelectron tomography. We identify four hierarchical levels of organization. Rhodopsin forms dimers; at least ten dimers form a row. Rows form pairs (tracks) that are aligned parallel to the disk incisures. Particle-based simulation shows that the combination of tracks with fast precomplex formation, i.e. rapid association and dissociation between inactive rhodopsin and the G protein transducin, leads to kinetic trapping: rhodopsin first activates transducin from its own track, whereas recruitment of transducin from other tracks proceeds more slowly. The trap mechanism could produce uniform single-photon responses independent of rhodopsin lifetime. In general, tracks might provide a platform that coordinates the spatiotemporal interaction of signaling molecules.

Individuals with type 1 diabetes mellitus are at high risk of developing atherosclerotic cardiovascular disease, but the underlying mechanisms by which type 1 diabetes accelerates atherosclerosis remain unknown. Increased retention of low-density lipoprotein (LDL) in atherosclerosis-prone sites of the diabetic vascular wall has been suggested, but direct evidence is lacking. In the present study, we investigated whether retention of LDL is increased in atherosclerotic-prone areas using a murine model of type 1 diabetes.

Among cyclic nucleotide phosphodiesterases (PDEs), PDE6 is unique in serving as an effector enzyme in G protein-coupled signal transduction. In retinal rods and cones, PDE6 is membrane-bound and activated to hydrolyse its substrate, cGMP, by binding of two active G protein α-subunits (Gα*). To investigate the activation mechanism of mammalian rod PDE6, we have collected functional and structural data, and analysed them by reaction-diffusion simulations. Gα* titration of membrane-bound PDE6 reveals a strong functional asymmetry of the enzyme with respect to the affinity of Gα* for its two binding sites on membrane-bound PDE6 and the enzymatic activity of the intermediary 1 : 1 Gα* · PDE6 complex. Employing cGMP and its 8-bromo analogue as substrates, we find that Gα* · PDE6 forms with high affinity but has virtually no cGMP hydrolytic activity. To fully activate PDE6, it takes a second copy of Gα* which binds with lower affinity, forming Gα* · PDE6 · Gα*. Reaction-diffusion simulations show that the functional asymmetry of membrane-bound PDE6 constitutes a coincidence switch and explains the lack of G protein-related noise in visual signal transduction. The high local concentration of Gα* generated by a light-activated rhodopsin molecule efficiently activates PDE6, whereas the low density of spontaneously activated Gα* fails to activate the effector enzyme.

Interleukin-1beta (IL-1beta) and tumor necrosis factor-alpha (TNF-alpha) are expressed by microglia and infiltrating macrophages following ischemic stroke. Whereas IL-1beta is primarily neurotoxic in ischemic stroke, TNF-alpha may have neurotoxic and/or neuroprotective effects. We investigated whether IL-1beta and TNF-alpha are synthesized by overlapping or segregated populations of cells after ischemic stroke in mice.

Antiviral therapy is urgently needed to combat the coronavirus disease 2019 (COVID-19) pandemic, which is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The protease inhibitor camostat mesylate inhibits SARS-CoV-2 infection of lung cells by blocking the virus-activating host cell protease TMPRSS2. Camostat mesylate has been approved for treatment of pancreatitis in Japan and is currently being repurposed for COVID-19 treatment. However, potential mechanisms of viral resistance as well as camostat mesylate metabolization and antiviral activity of metabolites are unclear. Here, we show that SARS-CoV-2 can employ TMPRSS2-related host cell proteases for activation and that several of them are expressed in viral target cells. However, entry mediated by these proteases was blocked by camostat mesylate. The camostat metabolite GBPA inhibited the activity of recombinant TMPRSS2 with reduced efficiency as compared to camostat mesylate and was rapidly generated in the presence of serum. Importantly, the infection experiments in which camostat mesylate was identified as a SARS-CoV-2 inhibitor involved preincubation of target cells with camostat mesylate in the presence of serum for 2 h and thus allowed conversion of camostat mesylate into GBPA. Indeed, when the antiviral activities of GBPA and camostat mesylate were compared in this setting, no major differences were identified. Our results indicate that use of TMPRSS2-related proteases for entry into target cells will not render SARS-CoV-2 camostat mesylate resistant. Moreover, the present and previous findings suggest that the peak concentrations of GBPA established after the clinically approved camostat mesylate dose (600 mg/day) will result in antiviral activity.

One of the main challenges in small molecule drug discovery is finding novel chemical compounds with desirable properties. In this work, we propose a novel method that combines in silico prediction of molecular properties such as biological activity or pharmacokinetics with an in silico optimization algorithm, namely Particle Swarm Optimization. Our method takes a starting compound as input and proposes new molecules with more desirable (predicted) properties. It navigates a machine-learned continuous representation of a drug-like chemical space guided by a defined objective function. The objective function combines multiple in silico prediction models, defined desirability ranges and substructure constraints. We demonstrate that our proposed method is able to consistently find more desirable molecules for the studied tasks in relatively short time. We hope that our method can support medicinal chemists in accelerating and improving the lead optimization process.

Host proteases have been suggested to be crucial for dissemination of MERS, SARS-CoV, and SARS-CoV-2 coronaviruses, but the relative contribution of membrane versus intracellular proteases remains controversial. Transmembrane serine protease 2 (TMPRSS2) is regarded as one of the main proteases implicated in the coronavirus S protein priming, an important step for binding of the S protein to the angiotensin-converting enzyme 2 (ACE2) receptor before cell entry. The main cellular location where the SARS-CoV-2 S protein priming occurs remains debatable, therefore hampering the development of targeted treatments. Herein, we identified the human extracellular serine protease inhibitor (serpin) alpha 1 antitrypsin (A1AT) as a novel TMPRSS2 inhibitor. Structural modeling revealed that A1AT docked to an extracellular domain of TMPRSS2 in a conformation that is suitable for catalysis, resembling similar serine protease-inhibitor complexes. Inhibitory activity of A1AT was established in a SARS-CoV-2 viral load system. Notably, plasma A1AT levels were associated with COVID-19 disease severity. Our data support the key role of extracellular serine proteases in SARS-CoV-2 infections and indicate that treatment with serpins, particularly the FDA-approved drug A1AT, may be effective in limiting SARS-CoV-2 dissemination by affecting the surface of the host cells.

The axon initial segment (AIS) is enriched in specific adaptor, cytoskeletal, and transmembrane molecules. During AIS establishment, a membrane diffusion barrier is formed between the axonal and somatodendritic domains. Recently, an axonal periodic pattern of actin, spectrin, and ankyrin forming 190-nm-spaced, ring-like structures has been discovered. However, whether this structure is related to the diffusion barrier function is unclear. Here, we performed single-particle tracking time-course experiments on hippocampal neurons during AIS development. We analyzed the mobility of lipid-anchored molecules by high-speed single-particle tracking and correlated positions of membrane molecules with the nanoscopic organization of the AIS cytoskeleton. We observe a strong reduction in mobility early in AIS development. Membrane protein motion in the AIS plasma membrane is confined to a repetitive pattern of ∼190-nm-spaced segments along the AIS axis as early as day in vitro 4, and this pattern alternates with actin rings. Mathematical modeling shows that diffusion barriers between the segments significantly reduce lateral diffusion along the axon.

The trans-membrane protease serine 2 (TMPRSS2) is essential for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) cell entry and infection. Efficacy and safety of TMPRSS2 inhibitors in patients with coronavirus disease 2019 (Covid-19) have not been evaluated in randomized trials.