Post-translational modification sites are present in hydrophilic cavities of alpha-synuclein, tau, FUS, and TDP-43 fibrils: A molecular dynamics study

Hydration plays a crucial role in the refolding of intrinsically disordered proteins into amyloid fibrils; however, the specific interactions between water and protein that may contribute to this process are still unknown. In our previous studies of alpha-synuclein (aSyn), we have shown that waters confined in fibril cavities are stabilizing features of this pathological fold; and that amino acids that hydrogen bond with these confined waters modulate primary and seeded aggregation. Here, we extend our aSyn molecular dynamics (MD) simulations with three new polymorphs and correlate MD trajectory information with known post-translational modifications (PTMs) and experimental data. We show that cavity residues are more evolutionarily conserved than non-cavity residues and are enriched with PTM sites. As expected, the confinement within hydrophilic cavities results in more stably hydrated amino acids. Interestingly, cavity PTM sites display the longest protein-water hydrogen bond lifetimes, three-fold greater than non-PTM cavity sites. Utilizing the deep mutational screen dataset by Newberry et al. and the Thioflavin T aggregation review by Pancoe et al. parsed using a fibril cavity/non-cavity definition, we show that hydrophobic changes to amino acids in cavities have a larger effect on fitness and aggregation rate than residues outside cavities, supporting our hypothesis that these sites are involved in the inhibition of aSyn toxic fibrillization. Finally, we expand our study to include analysis of fibril structures of tau, FUS, TDP-43, prion, and hnRNPA1; all of which contained hydrated cavities, with tau, FUS, and TDP-43 recapitulating our PTM results in aSyn fibril cavities.

Nimesulide, a COX-2 inhibitor, sensitizes pancreatic cancer cells to TRAIL-induced apoptosis by promoting DR5 clustering †

Nimesulide is a nonsteroidal anti-inflammatory drug and a COX-2 inhibitor with antitumor and antiproliferative activities that induces apoptosis in oral, esophagus, breast, and pancreatic cancer cells. Despite being removed from the market due to hepatotoxicity, nimesulide is still an important research tool being used to develop new anticancer drugs. Multiple studies have been done to modify the nimesulide skeleton to develop more potent anticancer agents and related compounds are promising scaffolds for future development. As such, establishing a mechanism of action for nimesulide remains an important part of realizing its potential. Here, we show that nimesulide enhances TRAIL-induced apoptosis in resistant pancreatic cancer cells by promoting clustering of DR5 in the plasma membrane. In this way, nimesulide acts like a related compound, DuP-697, which sensitizes TRAIL-resistant colon cancer cells in a similar manner. Our approach applies a time-resolved FRET-based biosensor that monitors DR5 clustering and conformational states in the plasma membrane. We show that this tool can be used for future high-throughput screens to identify novel, nontoxic small molecule scaffolds to overcome TRAIL resistance in cancer cells.

Proteasomal Stimulation by MK886 and Its Derivatives Can Rescue Tau-Induced Neurite Pathology

Proteasomal degradation of intrinsically disordered proteins, such as tau, is a critical component of proteostasis in both aging and neurodegenerative diseases. In this study, we investigated proteasomal activation by MK886 (MK). We previously identified MK as a lead compound capable of modulating tau oligomerization in a cellular FRET assay and rescuing P301L tau-induced cytotoxicity. We first confirmed robust proteasomal activation by MK using 20S proteasomal assays and a cellular proteasomal tau-GFP cleavage assay. We then show that MK treatment can significantly rescue tau-induced neurite pathology in differentiated SHSY5Y neurospheres. Due to this compelling result, we designed a series of seven MK analogs to determine if proteasomal activity is sensitive to structural permutations. Using the proteasome as the primary MOA, we examined tau aggregation, neurite outgrowth, inflammation, and autophagy assays to identify two essential substituents of MK that are required for compound activity: (1) removal of the N-chlorobenzyl group from MK negated both proteasomal and autophagic activity and reduced neurite outgrowth; and (2) removal of the indole-5-isopropyl group significantly improved neurite outgrowth and autophagy activity but reduced its anti-inflammatory capacity. Overall, our results suggest that the combination of proteasomal/autophagic stimulation and anti-inflammatory properties of MK and its derivatives can decrease tau-tau interactions and help rebalance dysfunctional proteostasis. Further development of MK to optimize its proteasomal, autophagic, and anti-inflammatory targets may lead to a novel therapeutic that would be beneficial in aging and neurodegenerative diseases.

Zafirlukast Is a Promising Scaffold for Selectively Inhibiting TNFR1 Signaling

Tumor necrosis factor (TNF) plays an important role in the pathogenesis of inflammatory and autoimmune diseases such as rheumatoid arthritis and Crohn's disease. The biological effects of TNF are mediated by binding to TNF receptors, TNF receptor 1 (TNFR1), or TNF receptor 2 (TNFR2), and this coupling makes TNFR1-specific inhibition by small-molecule therapies essential to avoid deleterious side effects. Recently, we engineered a time-resolved fluorescence resonance energy transfer biosensor for high-throughput screening of small molecules that modulate TNFR1 conformational states and identified zafirlukast as a compound that inhibits receptor activation, albeit at low potency. Here, we synthesized 16 analogues of zafirlukast and tested their potency and specificity for TNFR1 signaling. Using cell-based functional assays, we identified three analogues with significantly improved efficacy and potency, each of which induces a conformational change in the receptor (as measured by fluorescence resonance energy transfer (FRET) in cells). The best analogue decreased NF-κB activation by 2.2-fold, IκBα efficiency by 3.3-fold, and relative potency by two orders of magnitude. Importantly, we showed that the analogues do not block TNF binding to TNFR1 and that binding to the receptor's extracellular domain is strongly cooperative. Despite these improvements, the best candidate's maximum inhibition of NF-κB is only 63%, leaving room for further improvements to the zafirlukast scaffold to achieve full inhibition and prove its potential as a therapeutic lead. Interestingly, while we find that the analogues also bind to TNFR2 in vitro, they do not inhibit TNFR2 function in cells or cause any conformational changes upon binding. Thus, these lead compounds should also be used as reagents to study conformational-dependent activation of TNF receptors.

Discovery of a Non-competitive TNFR1 Antagonist Affibody with Picomolar Monovalent Potency That Does Not Affect TNFR2 Function

Tumor necrosis factor (TNF) is a key regulator of immune responses and plays a significant role in the initiation and maintenance of inflammation. Upregulation of TNF expression leads to several inflammatory diseases, such as Crohn's, ulcerative colitis, and rheumatoid arthritis. Despite the clinical success of anti-TNF treatments, the use of these therapies is limited because they can induce adverse side effects through inhibition of TNF biological activity, including blockade of TNF-induced immunosuppressive function of TNFR2. Using yeast display, we identified a synthetic affibody ligand (ABYTNFR1-1) with high binding affinity and specificity for TNFR1. Functional assays showed that the lead affibody potently inhibits TNF-induced NF-κB activation (IC50 of 0.23 nM) and, crucially, does not block the TNFR2 function. Additionally, ABYTNFR1-1 acts non-competitively─it does not block TNF binding or inhibit receptor-receptor interactions in pre-ligand-assembled dimers─thereby enhancing inhibitory robustness. The mechanism, monovalent potency, and affibody scaffold give this lead molecule uniquely strong potential as a therapeutic candidate for inflammatory diseases.

pH Effects Can Dominate Chemical Shift Perturbations in 1H,15N-HSQC NMR Spectroscopy for Studies of Small Molecule/α-Synuclein Interactions

¹H,¹⁵N-Heteronuclear Single Quantum Coherence (HSQC) NMR is a powerful technique that has been employed to characterize small-molecule interactions with intrinsically disordered monomeric α-Synuclein (aSyn). We report how solution pH can impact the interpretation of aSyn HSQC NMR spectra and demonstrate that small-molecule formulations (e.g., complexation with acidic salts) can lower sample pH and confound interpretation of drug binding and concomitant protein structural changes. Through stringent pH control, we confirm that several previously identified compounds (EGCG, Baicalin, and Dopamine (DOPA)) as well as a series of potent small-molecule inhibitors of aSyn pathology (Demeclocycline, Ro90-7501, and (±)-Bay K 8644) are capable of direct target engagement of aSyn. Previously, DOPA-aSyn interactions have been shown to elicit a dramatic chemical shift perturbation (CSP) localized to aSyn's H50 at low DOPA concentrations then expanding to aSyn's acidic C-terminal residues at increasing DOPA levels. Interestingly, this CSP profile mirrors our pH titration, where a small reduction in pH affects H50 CSP, and large pH changes induce robust C-terminal CSP. In contrast, under tightly controlled pH 5.0, DOPA induces significant CSPs observed at both ionizable and nonionizable residues. These results suggest that previous interpretations of DOPA-aSyn interactions were conflated with pH-induced CSP, highlighting the need for stringent pH control to minimize potential false-positive interpretations of ligand interactions in HSQC NMR experiments. Furthermore, DOPA's preferential interaction with aSyn under acidic pH represents a novel understanding of DOPA-aSyn interactions that may provide insight into the potential gain of toxic function of aSyn misfolding in α-synucleinopathies.

Indirubin Inhibits TRAIL-Induced Activation of Death Receptor 5 in Jurkat Cells

Death receptor 5 (DR5) is an apoptosis-inducing membrane receptor that mediates cell death in several life-threatening conditions. There is a crucial need for the discovery of DR5 antagonists for the therapeutic intervention of conditions in which the overactivation of DR5 underlies the pathophysiology. DR5 activation mediates cell death in non-alcoholic fatty liver disease (NAFLD) and neurodegenerative processes including amyloid-beta (Aβ) accumulation, spinal cord injury (SCI), and brain ischemia. In the current work, we used fluorescence resonance energy transfer (FRET) to monitor the conformational dynamics of DR5 that mediate death signaling. We used a time-resolved FRET screening platform to screen the Selleck library of 2863 U.S. Food and Drug Administration (FDA)-approved compounds. The high-throughput screen (HTS) identified 13 compounds that modulated the FRET between DR5 monomers beyond 5 median absolute deviations (MADs) from the DMSO controls. Of these 13 compounds, indirubin was identified to specifically inhibit tumor necrosis factor-related apoptosis-inducing ligand (TRAIL)-induced caspase-8 activity without modulating DR5 surface expression or TRAIL binding. Indirubin inhibited Fas-associated death domain (FADD) oligomerization and increased cellular FLICE-inhibitory protein (c-FLIP) expression; both are molecular mechanisms involved in inhibiting the DR5 signaling cascade. This study has elucidated previously unknown properties of indirubin that make it a promising candidate for therapeutic investigation of diseases in which overactivation of DR5 underlies pathology.

Advancements in a FRET Biosensor for Live-Cell Fluorescence-Lifetime High-Throughput Screening of Alpha-Synuclein

There is a critical need for small molecules capable of rescuing pathophysiological phenotypes induced by alpha-synuclein (aSyn) misfolding and oligomerization. Building upon our previous aSyn cellular fluorescence lifetime (FLT)-Förster resonance energy transfer (FRET) biosensors, we have developed an inducible cell model incorporating the red-shifted mCyRFP1/mMaroon1 (OFP/MFP) FRET pair. This new aSyn FRET biosensor improves the signal-to-noise ratio, reduces nonspecific background FRET, and results in a 4-fold increase (transient transfection) and 2-fold increase (stable, inducible cell lines) in FRET signal relative to our previous GFP/RFP aSyn biosensors. The inducible system institutes greater temporal control and scalability, allowing for fine-tuning of biosensor expression and minimizes cellular cytotoxicity due to overexpression of aSyn. Using these inducible aSyn-OFP/MFP biosensors, we screened the Selleck library of 2684 commercially available, FDA-approved compounds and identified proanthocyanidins and casanthranol as novel hits. Secondary assays validated the ability of these compounds to modulate aSyn FLT-FRET. Functional assays probing cellular cytotoxicity and aSyn fibrillization demonstrated their capability to inhibit seeded aSyn fibrillization. Proanthocyanidins completely rescued aSyn fibril-induced cellular toxicity with EC50 of 200 nM and casanthranol supported a 85.5% rescue with a projected EC50 of 34.2 μM. Furthermore, proanthocyanidins provide a valuable tool compound to validate our aSyn biosensor performance in future high-throughput screening campaigns of industrial-scale (million-compound) chemical libraries.

Broad Tricyclic Ring Inhibitors Block SARS-CoV-2 Spike Function Required for Viral Entry

The entry of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) into host cells requires binding of the viral spike glycoprotein to the angiotensin-converting enzyme 2 (ACE2) receptor, which triggers subsequent conformational changes to facilitate viral and cellular fusion at the plasma membrane or following endocytosis. Here, we experimentally identified selective and broad inhibitors of SARS-CoV-2 entry that share a tricyclic ring (or similar) structure. The inhibitory effect was restricted to early steps during infection and the entry inhibitors interacted with the receptor binding domain of the SARS-CoV-2 spike but did not significantly interfere with receptor (ACE2) binding. Instead, some of these compounds induced conformational changes or affected spike assembly and blocked SARS-CoV-2 spike cell-cell fusion activity. The broad inhibitors define a highly conserved binding pocket that is present on the spikes of SARS-CoV-1, SARS-CoV-2, and all circulating SARS-CoV-2 variants tested and block SARS-CoV spike activity required for mediating viral entry. These compounds provide new insights into the SARS-CoV-2 spike topography, as well as into critical steps on the entry pathway, and can serve as lead candidates for the development of broad-range entry inhibitors against SARS-CoVs.

Fluorescence Lifetime Measurement of Prefibrillar Sickle Hemoglobin Oligomers as a Platform for Drug Discovery in Sickle Cell Disease

The molecular origin of sickle cell disease (SCD) has been known since 1949, but treatments remain limited. We present the first high-throughput screening (HTS) platform for discovering small molecules that directly inhibit sickle hemoglobin (HbS) oligomerization and improve blood flow, potentially overcoming a long-standing bottleneck in SCD drug discovery. We show that at concentrations far below the threshold for nucleation and rapid polymerization, deoxygenated HbS forms small assemblies of multiple α2β2 tetramers. Our HTS platform leverages high-sensitivity fluorescence lifetime measurements that monitor these temporally stable prefibrillar HbS oligomers. We show that this approach is sensitive to compounds that inhibit HbS polymerization with or without modulating hemoglobin oxygen binding affinity. We also report the results of a pilot small-molecule screen in which we discovered and validated several novel inhibitors of HbS oligomerization.

Threonine Cavities Are Targetable Motifs That Control Alpha-Synuclein Fibril Growth

Recent high-resolution structures of alpha-synuclein (aSyn) fibrils offer promise for rational approaches to drug discovery for Parkinson's disease and Lewy body dementia. Harnessing the first such structures, we previously used molecular dynamics and free energy calculations to suggest that threonines 72 and 75─which line water-filled cavities within the fibril stacks─may be of central importance in stabilizing fibrils. Here, we used experimental mutagenesis of both wild-type and A53T aSyn to show that both threonine residues play important but surprisingly disparate roles in fibril nucleation and elongation. The T72A mutant, but not T75A, resulted in a large increase in the extent of fibrillization during primary nucleation, leading us to posit that T72 acts as a "brake" on run-away aggregation. An expanded set of simulations of five recent high-resolution fibril structures suggests that confinement of cavity waters around T72 correlates with this finding. In contrast, the T75A mutation led to a modest decrease in the extent of fibrillization. Furthermore, both T72A and T75A completely blocked the initial fibril elongation in seeded fibrillization. To test whether these threonine-lined cavities are druggable targets, we used computational docking to identify potential small-molecule binders. We show that the top-scoring hit, aprepitant, strongly promotes fibril growth while specifically interacting with aSyn fibrils and not monomer, and we offer speculation as to how such compounds could be used therapeutically.

The role of wild-type tau in Alzheimer's disease and related tauopathies

Tau oligomers have recently emerged as the principal toxic species in Alzheimer's disease (AD) and tauopathies. Tau oligomers are spontaneously self-assembled soluble tau proteins that are formed prior to fibrils, and they have been shown to play a central role in neuronal cell death and in the induction of neurodegeneration in animal models. As the therapeutic paradigm shifts to targeting toxic tau oligomers, this suggests the focus to study tau oligomerization in species that are less susceptible to fibrillization. While truncated and mutation containing tau as well as the isolated repeat domains are particularly prone to fibrillization, the wild-type (WT) tau proteins have been shown to be resistant to fibril formation in the absence of aggregation inducers. In this review, we will summarize and discuss the toxicity of WT tau both in vitro and in vivo, as well as its involvement in tau oligomerization and cell-to-cell propagation of pathology. Understanding the role of WT tau will enable more effective biomarker development and therapeutic discovery for treatment of AD and tauopathies.

Noncompetitive Allosteric Antagonism of Death Receptor 5 by a Synthetic Affibody Ligand

Fatty acid-induced upregulation of death receptor 5 (DR5) and its cognate ligand, tumor necrosis factor-related apoptosis-inducing ligand (TRAIL), promotes hepatocyte lipoapoptosis, which is a key mechanism in the progression of fatty liver disease. Accordingly, inhibition of DR5 signaling represents an attractive strategy for treating fatty liver disease. Ligand competition strategies are prevalent in tumor necrosis factor receptor antagonism, but recent studies have suggested that noncompetitive inhibition through perturbation of the receptor conformation may be a compelling alternative. To this end, we used yeast display and a designed combinatorial library to identify a synthetic 58-amino acid affibody ligand that specifically binds DR5. Biophysical and biochemical studies show that the affibody neither blocks TRAIL binding nor prevents the receptor-receptor interaction. Live-cell fluorescence lifetime measurements indicate that the affibody induces a conformational change in transmembrane dimers of DR5 and favors an inactive state of the receptor. The affibody inhibits apoptosis in TRAIL-treated Huh-7 cells, an in vitro model of fatty liver disease. Thus, this lead affibody serves as a potential drug candidate, with a unique mechanism of action, for fatty liver disease.

Discovery of Small Molecule Inhibitors of Huntingtin Exon 1 Aggregation by FRET-Based High-Throughput Screening in Living Cells

Huntington's disease (HD) is the most common inherited neurodegenerative disorder and one of the nine polyglutamine (polyQ) diseases. HD is characterized by the pathological aggregation of the misfolded huntingtin exon 1 protein (Httex1) with abnormally long polyQ expansion due to genetic mutation. While there is currently no effective treatment for HD, inhibition of aggregate formation represents a direct approach in mediating the toxicity associated with Httex1 misfolding. To exploit this therapeutic window, we engineered two fluorescence resonance energy transfer (FRET) based biosensors that monitor the aggregation of Httex1 with different expanded Q-lengths (Q39 and Q72) in living cells. These FRET biosensors, together with a high-precision fluorescence lifetime detection platform, enable high-throughput screening of small molecules that target Httex1 aggregation. We found six small molecules that decreased the FRET of the biosensors and reduced Httex1-Q72-induced neuronal cytotoxicity in N2a cells with nanomolar potency. Using advanced SPR and EPR techniques, we confirmed that the compounds directly bind to Httex1 fibrils and inhibit aggregate formation. This strategy in targeting the Httex1 aggregates can be applicable to other proteins involved in polyQ related diseases.

How Do Ethanolamine Plasmalogens Contribute to Order and Structure of Neurological Membranes?

Ethanolamine plasmalogen (EtnPLA) is a conical-shaped ether lipid and an essential component of neurological membranes. Low stability against oxidation limits its study in experiments. The concentration of EtnPLA in the bilayer varies depending on cell type and disease progression. Here we report on mixed bilayers of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) and 1-(1Z-octadecenyl)-2-oleoyl-sn-glycero-3-phosphoethanolamine (C18(Plasm)-18:1PE, PLAPE), an EtnPLA lipid subtype, at mole ratios of 2:1, 1:1, and 1:2. We present X-ray diffuse scattering (XDS) form factors F(qz) from oriented stacks of bilayers, related electron-density profiles, and hydrocarbon chain NMR order parameters. To aid future research on EtnPLA lipids and associated proteins, we have also extended the CHARMM36 all-atom force field to include the PLAPE lipid. The ability of the new force-field parameters to reproduce both X-ray and NMR structural properties of the mixed bilayer is remarkable. Our results indicate a thickening of the bilayer upon incorporation of increasing amounts of PLAPE into mixed bilayers, a reduction of lateral area per molecule, and an increase in lipid tail-ordering. The lateral compressibility modulus (KA) calculated from simulations yielded values for PLAPE similar to POPC.

Conformational states of TNFR1 as a molecular switch for receptor function

Tumor necrosis factor receptor 1 (TNFR1) is a transmembrane receptor that plays a key role in the regulation of the inflammatory pathway. While inhibition of TNFR1 has been the focus of many studies for the treatment of autoimmune diseases such as rheumatoid arthritis, activation of the receptor is important for the treatment of immunodeficiency diseases such as HIV and neurodegenerative diseases such as Alzheimer's disease where a boost in immune signaling is required. In addition, activation of other TNF receptors such as death receptor 5 or FAS receptor is important for cancer therapy. Here, we used a previously established TNFR1 fluorescence resonance energy transfer (FRET) biosensor together with a fluorescence lifetime technology as a high-throughput screening platform to identify a novel small molecule that activates TNFR1 by increasing inter-monomeric spacing in a ligand-independent manner. This shows that the conformational rearrangement of pre-ligand assembled receptor dimers can determine the activity of the receptor. By probing the interaction between the receptor and its downstream signaling molecule (TRADD) our findings support a new model of TNFR1 activation in which varying conformational states of the receptor act as a molecular switch in determining receptor function.

Targeting the ensemble of heterogeneous tau oligomers in cells: A novel small molecule screening platform for tauopathies

OBJECTIVE

Understanding the heterogeneous pathology in Alzheimer's disease and related tauopathies is one of the most urgent and fundamental challenges facing the discovery of novel disease-modifying therapies. Through monitoring ensembles of toxic and nontoxic tau oligomers spontaneously formed in cells, our biosensor technology can identify tool compounds that modulate tau oligomer structure and toxicity, providing much needed insight into the nature and properties of toxic tau oligomers.

BACKGROUND

Tauopathies are a group of neurodegenerative disorders characterized by pathologic aggregation of the microtubule binding protein tau. Recent studies suggest that tau oligomers are the primary toxic species in tauopathies.

NEW/UPDATED HYPOTHESIS

We hypothesize that tau biosensors capable of monitoring tau oligomer conformation are able to identify tool compounds that modulate the structure and conformation of these tau assemblies, providing key insight into the unique structural fingerprints of toxic tau oligomers. These fingerprints will provide gravely needed biomarker profiles to improve staging of early tauopathy pathology and generate lead compounds for potential new therapeutics. Our time-resolved fluorescence resonance energy transfer biosensors provide us an exquisitely sensitive technique to monitor minute structural changes in monomer and oligomer conformation. In this proof-of-concept study, we identified a novel tool compound, MK-886, which directly binds tau, perturbs the conformation of toxic tau oligomers, and rescues tau-induced cytotoxicity. Furthermore, we show that MK-886 alters the conformation of tau monomer at the proline-rich and microtubule binding regions, stabilizing an on-pathway oligomer.

MAJOR CHALLENGES FOR THE HYPOTHESIS

Our approach monitors changes in the ensemble of assemblies that are spontaneously formed in cells but does not specifically isolate or enrich unique toxic tau species. However, time-resolved fluorescence resonance energy transfer does not provide high-resolution, atomic scale information, requiring additional experimental techniques to resolve the structural features stabilized by different tool compounds.

LINKAGE TO OTHER MAJOR THEORIES

Our biosensor technology is broadly applicable to other areas of tauopathy therapeutic development. These biosensors can be readily modified for different isoforms of tau, specific post-translational modifications, and familial Alzheimer's disease-associated mutations. We are eager to explore tau interactions with chaperone proteins, monitor cross-reactivity with other intrinsically disordered proteins, and target seeded oligomer pathology.

Multiscale Computational Modeling of Tubulin-Tubulin Lateral Interaction

Microtubules are multistranded polymers in eukaryotic cells that support key cellular functions such as chromosome segregation, motor-based cargo transport, and maintenance of cell polarity. Microtubules self-assemble via "dynamic instability," in which the dynamic plus ends switch stochastically between alternating phases of polymerization and depolymerization. A key question in the field is what are the atomistic origins of this switching, i.e., what is different between the GTP- and GDP-tubulin states that enables microtubule growth and shortening, respectively? More generally, a major challenge in biology is how to connect theoretical frameworks across length- and timescales, from atoms to cellular behavior. In this study, we describe a multiscale model by linking atomistic molecular dynamics (MD), molecular Brownian dynamics (BD), and cellular-level thermokinetic modeling of microtubules. Here, we investigated the underlying interaction energy when tubulin dimers associate laterally by performing all-atom MD simulations. We found that the lateral potential energy is not significantly different among three nucleotide states of tubulin, GTP, GDP, and GMPCPP and is estimated to be ≅ -11 k_BT. Furthermore, using MD potential energy in our BD simulations of tubulin dimers confirms that the lateral bond is weak on its own, with a mean lifetime of ∼0.1 μs, implying that the longitudinal bond is required for microtubule assembly. We conclude that nucleotide-dependent lateral-bond strength is not the key mediator microtubule dynamic instability, implying that GTP acts elsewhere to exert its stabilizing influence on microtubule polymer. Furthermore, the estimated lateral-bond strength (ΔG_lat⁰≅ -5 k_BT) is well-aligned with earlier estimates based on thermokinetic modeling and light microscopy measurements. Thus, we have computationally connected atomistic-level structural information, obtained by cryo-electron microscopy, to cellular-scale microtubule assembly dynamics using a combination of MD, BD, and thermokinetic models to bridge from Ångstroms to micrometers and from femtoseconds to minutes.

Noncompetitive inhibitors of TNFR1 probe conformational activation states

Tumor necrosis factor receptor 1 (TNFR1) is a central mediator of the inflammatory pathway and is associated with several autoimmune diseases such as rheumatoid arthritis. A revision to the canonical model of TNFR1 activation suggests that activation involves conformational rearrangements of preassembled receptor dimers. Here, we identified small-molecule allosteric inhibitors of TNFR1 activation and probed receptor dimerization and function. Specifically, we used a fluorescence lifetime-based high-throughput screen and biochemical, biophysical, and cellular assays to identify small molecules that noncompetitively inhibited the receptor without reducing ligand affinity or disrupting receptor dimerization. We also found that residues in the ligand-binding loop that are critical to the dynamic coupling between the extracellular and the transmembrane domains played a key gatekeeper role in the conformational dynamics associated with signal propagation. Last, using a simple structure-activity relationship analysis, we demonstrated that these newly found molecules could be further optimized for improved potency and specificity. Together, these data solidify and deepen the new model for TNFR1 activation.

Minimal Nucleation State of α-Synuclein Is Stabilized by Dynamic Threonine-Water Networks

The first structures of α-synuclein (αSyn) fibrils have recently been solved. Here, we use a unique combination of molecular dynamics simulation strategies to address the minimal nucleation size of the 11-amino acid NAC protofibril solved by X-ray and to interrogate the dynamic behavior of unexpected crystal waters in the steric zipper. We found that protofibrils of >8 chains are thermodynamically stabilized due to protection of the fibril core from solvent influx and ordering of the end strands by the fibril core. In these stable oligomers, water molecules resolved in the crystal structure freely exchange with bulk solvent but are, on average, stably coordinated along the β-sheet by inward-facing Thr72 and Thr75. We confirm the persistence of this water coordination via simulations of the full-length Greek-key structure solved by NMR and speculate that these Thr-water networks are important in the context of enhanced fibril nucleation in the familial A53T mutation.

Soluble Extracellular Domain of Death Receptor 5 Inhibits TRAIL-Induced Apoptosis by Disrupting Receptor-Receptor Interactions

Dysregulation of tumor necrosis factor (TNF) receptor signaling is a key feature of various inflammatory disorders. Current treatments for TNF-related diseases function either by sequestering ligand or blocking ligand-receptor interactions, which can cause dangerous side effects by inhibiting the receptors that are not involved in the disease condition. Thus, alternate strategies that target receptor-receptor interactions are needed. We hypothesized that the soluble extracellular domain (ECD) of long isoform of death receptor 5 (DR5) could block endogenous receptor assembly, mimicking the biological effect of decoy receptors that lack the death domain to trigger apoptosis. Using live-cell fluorescence resonance energy transfer studies, we demonstrated that soluble ECD disrupts endogenous DR5-DR5 interactions. Cell viability assays were used to demonstrate the complete inhibition of TNF-related apoptosis-inducing ligand (TRAIL)-induced apoptosis by the ECD, although TRAIL is still able to bind to the receptor. Importantly, we used mutagenesis to prove that the inhibition of TRAIL-induced apoptosis by the ECD predominantly comes from the disruption of DR5 oligomerization and not ligand sequestration. Inhibition of death receptor activation should have important therapeutic applications in diseases such as nonalcoholic fatty liver disease. More generally, this approach should be generalized to enable the inhibition of other TNF receptor signaling mechanisms that are associated in a wide range of clinical conditions.

An Innovative High-Throughput Screening Approach for Discovery of Small Molecules That Inhibit TNF Receptors

Tumor necrosis factor receptor 1 (TNFR1) is a transmembrane receptor that binds tumor necrosis factor or lymphotoxin-alpha and plays a critical role in regulating the inflammatory response. Upregulation of these ligands is associated with inflammatory and autoimmune diseases. Current treatments reduce symptoms by sequestering free ligands, but this can cause adverse side effects by unintentionally inhibiting ligand binding to off-target receptors. Hence, there is a need for new small molecules that specifically target the receptors, rather than the ligands. Here, we developed a TNFR1 FRET biosensor expressed in living cells to screen compounds from the NIH Clinical Collection. We used an innovative high-throughput fluorescence lifetime screening platform that has exquisite spatial and temporal resolution to identify two small-molecule compounds, zafirlukast and triclabendazole, that inhibit the TNFR1-induced IκBα degradation and NF-κB activation. Biochemical and computational docking methods were used to show that zafirlukast disrupts the interactions between TNFR1 pre-ligand assembly domain (PLAD), whereas triclabendazole acts allosterically. Importantly, neither compound inhibits ligand binding, proving for the first time that it is possible to inhibit receptor activation by targeting TNF receptor-receptor interactions. This strategy should be generally applicable to other members of the TNFR superfamily, as well as to oligomeric receptors in general.

α-Synuclein's Uniquely Long Amphipathic Helix Enhances its Membrane Binding and Remodeling Capacity

α-Synuclein is the primary protein found in Lewy bodies, the protein and lipid aggregates associated with Parkinson's disease and Lewy body dementia. The protein folds into a uniquely long amphipathic α-helix (AH) when bound to a membrane, and at high enough concentrations, it induces large-scale remodeling of membranes (tubulation and vesiculation). By engineering a less hydrophobic variant of α-Synuclein, we previously showed that the energy associated with binding of α-Synuclein's AH correlates with the extent of membrane remodeling (Braun et al. in J Am Chem Soc 136:9962-9972, 2014). In this study, we combine fluorescence correlation spectroscopy, electron microscopy, and vesicle clearance assays with coarse-grained molecular dynamics simulations to test the impact of decreasing the length of the amphipathic helix on membrane binding energy and tubulation. We show that truncation of α-Synuclein's AH length by approximately 15% reduces both its membrane binding affinity (by fivefold) and membrane remodeling capacity (by nearly 50% on per mole of bound protein basis). Results from simulations correlate well with the experiments and lend support to the idea that at high protein density there is a stabilization of individual, protein-induced membrane curvature fields. The extent to which these curvature fields are stabilized, a function of binding energy, dictates the extent of tubulation. Somewhat surprisingly, we find that this stabilization does not correlate directly with the geometric distribution of the proteins on the membrane surface.

Piecing it together: Unraveling the elusive structure-function relationship in single-pass membrane receptors

The challenge of crystallizing single-pass plasma membrane receptors has remained an obstacle to understanding the structural mechanisms that connect extracellular ligand binding to cytosolic activation. For example, the complex interplay between receptor oligomerization and conformational dynamics has been, historically, only inferred from static structures of isolated receptor domains. A fundamental challenge in the field of membrane receptor biology, then, has been to integrate experimentally observable dynamics of full-length receptors (e.g. diffusion and conformational flexibility) into static structural models of the disparate domains. In certain receptor families, e.g. the ErbB receptors, structures have led somewhat linearly to a putative model of activation. In other families, e.g. the tumor necrosis factor (TNF) receptors, structures have produced divergent hypothetical mechanisms of activation and transduction. Here, we discuss in detail these and other related receptors, with the goal of illuminating the current challenges and opportunities in building comprehensive models of single-pass receptor activation. The deepening understanding of these receptors has recently been accelerated by new experimental and computational tools that offer orthogonal perspectives on both structure and dynamics. As such, this review aims to contextualize those technological developments as we highlight the elegant and complex conformational communication between receptor domains. This article is part of a Special Issue entitled: Interactions between membrane receptors in cellular membranes edited by Kalina Hristova.

The IM30/Vipp1 C-terminus associates with the lipid bilayer and modulates membrane fusion

IM30/Vipp1 proteins are crucial for thylakoid membrane biogenesis in chloroplasts and cyanobacteria. A characteristic C-terminal extension distinguishes these proteins from the homologous bacterial PspA proteins, and this extension has been discussed to be key for the IM30/Vipp1 activity. Here we report that the extension of the Synechocystis IM30 protein is indispensable, and argue that both, the N-terminal PspA-domain as well as the C-terminal extension are needed in order for the IM30 protein to conduct its in vivo function. In vitro, we show that the PspA-domain of IM30 is vital for stability/folding and oligomer formation of IM30 as well as for IM30-triggered membrane fusion. In contrast, the IM30 C-terminal domain is involved in and necessary to stabilize defined contacts to negatively charged membrane surfaces, and to modulate the IM30-induced membrane fusion activity. Although the two IM30 protein domains have distinct functional roles, only together they enable IM30 to work properly.

Polyunsaturated chains in asymmetric lipids disorder raft mixtures and preferentially associate with α-Synuclein

Using molecular dynamics simulations, we have explored the effect of asymmetric lipids-specifically those that contain one polyunsaturated (PUFA) and one saturated fatty acid chain-on phase separation in heterogeneous membranes. These lipids are prevalent in neuronal membranes, particularly in synaptic membranes, where the Parkinson's Disease protein α-Synuclein (αS) is found. We have therefore explored the relationship between asymmetric, PUFA-containing lipids, and αS. The simulations show that asymmetric lipids partition to the liquid disordered (Ld) phase of canonical raft mixtures because of the highly disordered PUFA chain. In the case of a membrane built to mimic the lipid composition of a synaptic vesicle, the PUFA-containing asymmetric lipids completely disrupt phase separation. Because αS is positively charged, we show that it partitions with negatively charged lipids, regardless of the saturation state of the chains. Additionally, αS preferentially associates with the polyunsaturated fatty acid tails of both charged and neutral lipids. This is a consequence of those chains' ability to accommodate the void beneath the amphipathic helix. This article is part of a Special Issue entitled: Lipid order/lipid defects and lipid-control of protein activity edited by Dirk Schneider.

Death Receptor 5 Networks Require Membrane Cholesterol for Proper Structure and Function

Death receptor 5 (DR5) is an apoptosis-inducing member of the tumor necrosis factor receptor superfamily, whose activity has been linked to membrane cholesterol content. Upon ligand binding, DR5 forms large clusters within the plasma membrane that have often been assumed to be manifestations of receptor co-localization in cholesterol-rich membrane domains. However, we have recently shown that DR5 clusters are more than just randomly aggregated receptors. Instead, these are highly structured networks held together by receptor dimers. These dimers are stabilized by specific transmembrane helix-helix interactions, including a disulfide bond in the long isoform of the receptor. The complex relationships among DR5 network formation, transmembrane helix dimerization, membrane cholesterol, and receptor activity has not been established. It is unknown whether the membrane itself plays an active role in driving DR5 transmembrane helix interactions or in the formation of the networks. We show that cholesterol depletion in cells does not inhibit the formation of DR5 networks. However, the networks that form in cholesterol-depleted cells fail to induce caspase cleavage. These results suggest a potential structural difference between active and inactive networks. As evidence, we show that cholesterol is necessary for the covalent dimerization of DR5 transmembrane domains. Molecular simulations and experiments in synthetic vesicles on the DR5 transmembrane dimer suggest that dimerization is facilitated by increased helicity in a thicker bilayer.

Membrane remodeling and mechanics: Experiments and simulations of α-Synuclein

We review experimental and simulation approaches that have been used to determine curvature generation and remodeling of lipid bilayers by membrane-bending proteins. Particular emphasis is placed on the complementary approaches used to study α-Synuclein (αSyn), a major protein involved in Parkinson's disease (PD). Recent cellular and biophysical experiments have shown that the protein 1) deforms the native structure of mitochondrial and model membranes; and 2) inhibits vesicular fusion. Today's advanced experimental and computational technology has made it possible to quantify these protein-induced changes in membrane shape and material properties. Collectively, experiments, theory and multi-scale simulation techniques have established the key physical determinants of membrane remodeling and rigidity: protein binding energy, protein partition depth, protein density, and membrane tension. Despite the exciting and significant progress made in recent years in these areas, challenges remain in connecting biophysical insights to the cellular processes that lead to disease. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

HuR Contributes to TRAIL Resistance by Restricting Death Receptor 4 Expression in Pancreatic Cancer Cells

Pancreatic ductal adenocarcinoma (PDA) is one of the most lethal cancers, in part, due to resistance to both conventional and targeted therapeutics. TRAIL directly induces apoptosis through engagement of cell surface Death Receptors (DR4 and DR5), and has been explored as a molecular target for cancer treatment. Clinical trials with recombinant TRAIL and DR-targeting agents, however, have failed to show overall positive outcomes. Herein, we identify a novel TRAIL resistance mechanism governed by Hu antigen R (HuR, ELAV1), a stress-response protein abundant and functional in PDA cells. Exogenous HuR overexpression in TRAIL-sensitive PDA cell lines increases TRAIL resistance whereas silencing HuR in TRAIL-resistant PDA cells, by siRNA oligo-transfection, decreases TRAIL resistance. PDA cell exposure to soluble TRAIL induces HuR translocation from the nucleus to the cytoplasm. Furthermore, it is demonstrated that HuR interacts with the 3'-untranslated region (UTR) of DR4 mRNA. Pre-treatment of PDA cells with MS-444 (Novartis), an established small molecule inhibitor of HuR, substantially increased DR4 and DR5 cell surface levels and enhanced TRAIL sensitivity, further validating HuR's role in affecting TRAIL apoptotic resistance. NanoString analyses on the transcriptome of TRAIL-exposed PDA cells identified global HuR-mediated increases in antiapoptotic processes. Taken together, these data extend HuR's role as a key regulator of TRAIL-induced apoptosis.

IMPLICATIONS

Discovery of an important new HuR-mediated TRAIL resistance mechanism suggests that tumor-targeted HuR inhibition increases sensitivity to TRAIL-based therapeutics and supports their re-evaluation as an effective treatment for PDA patients. Mol Cancer Res; 14(7); 599-611. ©2016 AACR.

α-Synuclein Reduces Tension and Increases Undulations in Simulations of Small Unilamellar Vesicles

Using coarse-grained molecular dynamics simulations we have explored the effect of α-Synuclein (αSyn) on the structural and mechanical properties of small unilamellar vesicles in the fluid-phase. The study is motivated by observations that a high density of membrane-bound αSyn inhibits the fusion of synthetic small unilamellar vesicles. By combining three-dimensional pressure tensor calculations with our recently developed spherical harmonics fluctuation analysis approach, we show a reduction in membrane surface tension and increased membrane undulations when αSyn is bound to the vesicle's outer leaflet at a 200:1 L/P. The protein effects these changes by decreasing the negative pressure in the headgroup region of the outer leaflet and increasing the positive pressure throughout the hydrocarbon core.

Biophysics of α-synuclein induced membrane remodelling

α-Synuclein is an intrinsically disordered protein whose aggregation is a hallmark of Parkinson's disease. In neurons, α-synuclein is thought to play important roles in mediating both endo- and exocytosis of synaptic vesicles through interactions with either the lipid bilayer or other proteins. Upon membrane binding, the N-terminus of α-synuclein forms a helical structure and inserts into the hydrophobic region of the outer membrane leaflet. However, membrane structural changes induced by α-synuclein are still largely unclear. Here we report a substantial membrane area expansion induced by the binding of α-synuclein monomers. This measurement is accomplished by observing the increase of membrane area during the binding of α-synuclein to pipette-aspirated giant vesicles. The extent of membrane area expansion correlates linearly with the density of α-synuclein on the membrane, revealing a constant area increase induced by the binding per α-synuclein molecule. The area expansion per synuclein is found to exhibit a strong dependence on lipid composition, but is independent of membrane tension and vesicle size. Fragmentation or tubulation of the membrane follows the membrane expansion process. However, contrary to BAR domain proteins, no distinct tubulation-transition density can apparently be identified for α-synuclein, suggesting a more complex membrane curvature generation mechanism. Consideration of α-synuclein's membrane binding free energy and biophysical properties of the lipid bilayer leads us to conclude that membrane expansion by α-synuclein results in thinning of the bilayer. These membrane thinning and tubulation effects may underlie α-synuclein's role in mediating cell trafficking processes such as endo- and exocytosis.

Determining Structural and Mechanical Properties from Molecular Dynamics Simulations of Lipid Vesicles

We have developed an algorithm to determine membrane structure, area per lipid, and bending rigidity from molecular dynamics simulations of lipid vesicles. Current methods to extract structure from vesicle simulations define densities relative to the global center of mass of the vesicle. This approach ignores the long-wavelength fluctuations (undulations) that develop across the sphere and broaden the underlying structure. Our method establishes a local reference frame by defining a radially undulating reference surface (URS) and thereby removes the broadening effect of the undulations. Using an arc-length low-pass filter, we render the URS by defining the bilayer midplane on an equi-angular θ, ϕ-grid (colatitude, longitude). This surface is then expanded onto a truncated series of spherical harmonics. The spherical harmonic coefficients characterize the long-wavelength fluctuations that define both the local reference frame-used to determine the bilayer's structure-and the area per lipid (AL) along the undulating surface. Additionally, the resulting power spectrum of spherical harmonic coefficients can be fit to a Helfrich continuum model for membrane bending in spherical geometry to extract bending rigidity (kc). kc values determined for both DMPC and DMPC + cholesterol (30 mol %) vesicles are consistent with values from corresponding flat-patch systems determined using an independent, previously published spectral method. These new tools to accurately extract structure, AL, and kc should prove invaluable in evaluating the construction and equilibration of lipid vesicle simulations.

α-Synuclein-induced membrane remodeling is driven by binding affinity, partition depth, and interleaflet order asymmetry

We have investigated the membrane remodeling capacity of the N-terminal membrane-binding domain of α-synuclein (α-Syn₁₀₀). Using fluorescence correlation spectroscopy and vesicle clearance assays, we show that α-Syn₁₀₀ fully tubulates POPG vesicles, the first demonstration that the amphipathic helix on its own is capable of this effect. We also show that at equal density of membrane-bound protein, α-Syn has dramatically reduced affinity for, and does not tubulate, vesicles composed of a 1:1 POPG:POPC mixture. Coarse-grained molecular dynamics simulations suggested that the difference between the pure POPG and mixture results may be attributed to differences in the protein’s partition depth, the membrane’s hydrophobic thickness, and disruption of acyl chain order. To explore the importance of these attributes compared with the role of the reduced binding energy, we created an α-Syn₁₀₀ variant in which we removed the hydrophobic core of the non-amyloid component (NAC) domain and tested its impact on pure POPG vesicles. We observed a substantial reduction in binding affinity and tubulation, and simulations of the NAC-null protein suggested that the reduced binding energy increases the protein mobility on the bilayer surface, likely impacting the protein’s ability to assemble into organized pretubule structures. We also used simulations to explore a potential role for interleaflet coupling as an additional driving force for tubulation. We conclude that symmetry across the leaflets in the tubulated state maximizes the interaction energy of the two leaflets and relieves the strain induced by the hydrophobic void beneath the amphipathic helix.

Comparing simulations of lipid bilayers to scattering data: the GROMOS 43A1-S3 force field

Simulations of DOPC at T = 303 K were performed using the united atom force field 43A1-S3 at six fixed projected areas, A(P) = 62, 64, 66, 68, 70, and 72 Å(2), as well as a tensionless simulation that produced an average A(NPT) = 65.8 Å(2). After a small undulation correction for the system size consisting of 288 lipids, results were compared to experimental data. The best, and excellent, fit to neutron scattering data occurs at an interpolated A(N) = 66.6 Å(2) and the best, but not as good, fit to the more extensive X-ray scattering data occurs at A(X) = 68.7 Å(2). The distance ΔDB-H between the Gibbs dividing surface for water and the peak in the electron density profile agrees with scattering experiments. The calculated area compressibility K(A) = 277 ± 10 mN/m is in excellent agreement with the micromechanical experiment. The volume per lipid V(L) is smaller than volume experiments which suggests a workaround that raises all the areas by about 1.5%. Although A(X) ≠ A(N) ≠ A(NPT), this force field obtains acceptable agreement with experiment for A(L) = 67.5 Å(2) (68.5 Å(2) in the workaround), which we suggest is a better DOPC result from 43A1-S3 simulations than its value from the tensionless NPT simulation. However, nonsimulation modeling obtains better simultaneous fits to both kinds of scattering data, which suggests that the force fields can still be improved.

HuR's post-transcriptional regulation of Death Receptor 5 in pancreatic cancer cells

Apoptosis is one of the core signaling pathways disrupted in pancreatic ductal adenocarcinoma (PDA). Death receptor 5 (DR5) is a member of the tumor necrosis factor (TNF)-receptor superfamily that is expressed in cancer cells. Binding of TNF-related apoptosis-inducing ligand (TRAIL) to DR5 is a potent trigger of the extrinsic apoptotic pathway, and numerous clinical trials are based on DR5-targeted therapies for cancer, including PDA. Human antigen R (HuR), an RNA-binding protein, regulates a select number of transcripts under stress conditions. Here we report that HuR translocates from the nucleus to the cytoplasm of PDA cells upon treatment with a DR5 agonist. High doses of DR5 agonist induce cleavage of both HuR and caspase 8. HuR binds to DR5 mRNA at the 5'-untranslated region (UTR) in PDA cells in response to different cancer-associated stressors and subsequently represses DR5 protein expression; silencing HuR augments DR5 protein production by enabling its translation and thus enhances apoptosis. In PDA specimens (n = 53), negative HuR cytoplasmic expression correlated with elevated DR5 expression (odds ratio 16.1, p < 0.0001). Together, these data demonstrate a feedback mechanism elicited by HuR-mediated repression of the key apoptotic membrane protein DR5.

The methionine-aromatic motif plays a unique role in stabilizing protein structure

Of the 20 amino acids, the precise function of methionine (Met) remains among the least well understood. To establish a determining characteristic of methionine that fundamentally differentiates it from purely hydrophobic residues, we have used in vitro cellular experiments, molecular simulations, quantum calculations, and a bioinformatics screen of the Protein Data Bank. We show that approximately one-third of all known protein structures contain an energetically stabilizing Met-aromatic motif and, remarkably, that greater than 10,000 structures contain this motif more than 10 times. Critically, we show that as compared with a purely hydrophobic interaction, the Met-aromatic motif yields an additional stabilization of 1-1.5 kcal/mol. To highlight its importance and to dissect the energetic underpinnings of this motif, we have studied two clinically relevant TNF ligand-receptor complexes, namely TRAIL-DR5 and LTα-TNFR1. In both cases, we show that the motif is necessary for high affinity ligand binding as well as function. Additionally, we highlight previously overlooked instances of the motif in several disease-related Met mutations. Our results strongly suggest that the Met-aromatic motif should be exploited in the rational design of therapeutics targeting a range of proteins.