Antidepressant drugs act by directly binding to TRKB neurotrophin receptors.

It is unclear how binding of antidepressant drugs to their targets gives rise to the clinical antidepressant effect. We discovered that the transmembrane domain of tyrosine kinase receptor 2 (TRKB), the brain-derived neurotrophic factor (BDNF) receptor that promotes neuronal plasticity and antidepressant responses, has a cholesterol-sensing function that mediates synaptic effects of cholesterol. We then found that both typical and fast-acting antidepressants directly bind to TRKB, thereby facilitating synaptic localization of TRKB and its activation by BDNF. Extensive computational approaches including atomistic molecular dynamics simulations revealed a binding site at the transmembrane region of TRKB dimers. Mutation of the TRKB antidepressant-binding motif impaired cellular, behavioral, and plasticity-promoting responses to antidepressants in vitro and in vivo. We suggest that binding to TRKB and allosteric facilitation of BDNF signaling is the common mechanism for antidepressant action, which may explain why typical antidepressants act slowly and how molecular effects of antidepressants are translated into clinical mood recovery.

TrkB transmembrane domain: bridging structural understanding with therapeutic strategy.

TrkB (neuronal receptor tyrosine kinase-2, NTRK2) is the receptor for brain-derived neurotrophic factor (BDNF) and is a critical regulator of activity-dependent neuronal plasticity. The past few years have witnessed an increasing understanding of the structure and function of TrkB, including its transmembrane domain (TMD). TrkB interacts with membrane cholesterol, which bidirectionally regulates TrkB signaling. Additionally, TrkB has recently been recognized as a binding target of antidepressant drugs. A variety of different antidepressants, including typical and rapid-acting antidepressants, as well as psychedelic compounds, act as allosteric potentiators of BDNF signaling through TrkB. This suggests that TrkB is the common target of different antidepressant compounds. Although more research is needed, current knowledge suggests that TrkB is a promising target for further drug development.

High-content imaging and structure-based predictions reveal functional differences between Niemann-Pick C1 variants.

The human Niemann-Pick C1 (NPC1) gene encoding a 1278 amino acid protein is very heterogeneous. While some variants represent benign polymorphisms, NPC disease carriers and patients may possess rare variants, whose functional importance remains unknown. An NPC1 cDNA construct known as NPC1 wild-type variant (WT-V), distributed between laboratories and used as a WT control in several studies, also contains changes regarding specific amino acids compared to the NPC1 Genbank reference sequence. To improve the dissection of subtle functional differences, we generated human cells stably expressing NPC1 variants from the AAVS1 safe-harbor locus on an NPC1-null background engineered by CRISPR/Cas9 editing. We then employed high-content imaging with automated image analysis to quantitatively assess LDL-induced, time-dependent changes in lysosomal cholesterol content and lipid droplet formation. Our results indicate that the L472P change present in NPC1 WT-V compromises NPC1 functionality in lysosomal cholesterol export. All-atom molecular dynamics simulations suggest that the L472P change alters the relative position of the NPC1 middle and the C-terminal luminal domains, disrupting the recently characterized cholesterol efflux tunnel. These results reveal functional defects in NPC1 WT-V and highlight the strength of simulations and quantitative imaging upon stable protein expression in elucidating subtle differences in protein function.

Cholesterol-recognition motifs in the transmembrane domain of the tyrosine kinase receptor family: The case of TRKB.

Cholesterol is an essential constituent of cell membranes. The discovery of cholesterol-recognition amino acid consensus (CRAC) motif in proteins indicated a putative direct, non-covalent interaction between cholesterol and proteins. In the present study, we evaluated the presence of a CRAC motif and its inverted version (CARC) in the transmembrane region (TMR) of the tyrosine kinase receptor family (RTK) in several species using in silico methods. CRAC motifs were found across all species analyzed, while CARC was found only in vertebrates. The tropomyosin-related kinase B (TRKB), a member of the RTK family, through interaction with its endogenous ligand brain-derived neurotrophic factor (BDNF) is a core participant in the neuronal plasticity process and exhibits a CARC motif in its TMR. Upon identifying the conserved CARC motif in the TRKB, we performed molecular dynamics simulations of the mouse TRKB.TMR. The simulations indicated that cholesterol interaction with the TRKB CARC motif occurs mainly at the central Y433 residue. Our binding assay suggested a bell-shaped effect of cholesterol on BDNF interaction with TRKB receptors, and our results suggest that CARC/CRAC motifs may play a role in the function of the RTK family TMR.

Multiscale Simulations of Biological Membranes: The Challenge To Understand Biological Phenomena in a Living Substance.

Biological membranes are tricky to investigate. They are complex in terms of molecular composition and structure, functional over a wide range of time scales, and characterized by nonequilibrium conditions. Because of all of these features, simulations are a great technique to study biomembrane behavior. A significant part of the functional processes in biological membranes takes place at the molecular level; thus computer simulations are the method of choice to explore how their properties emerge from specific molecular features and how the interplay among the numerous molecules gives rise to function over spatial and time scales larger than the molecular ones. In this review, we focus on this broad theme. We discuss the current state-of-the-art of biomembrane simulations that, until now, have largely focused on a rather narrow picture of the complexity of the membranes. Given this, we also discuss the challenges that we should unravel in the foreseeable future. Numerous features such as the actin-cytoskeleton network, the glycocalyx network, and nonequilibrium transport under ATP-driven conditions have so far received very little attention; however, the potential of simulations to solve them would be exceptionally high. A major milestone for this research would be that one day we could say that computer simulations genuinely research biological membranes, not just lipid bilayers.

Reduced level of docosahexaenoic acid shifts GPCR neuroreceptors to less ordered membrane regions.

G protein-coupled receptors (GPCRs) control cellular signaling and responses. Many of these GPCRs are modulated by cholesterol and polyunsaturated fatty acids (PUFAs) which have been shown to co-exist with saturated lipids in ordered membrane domains. However, the lipid compositions of such domains extracted from the brain cortex tissue of individuals suffering from GPCR-associated neurological disorders show drastically lowered levels of PUFAs. Here, using free energy techniques and multiscale simulations of numerous membrane proteins, we show that the presence of the PUFA DHA helps helical multi-pass proteins such as GPCRs partition into ordered membrane domains. The mechanism is based on hybrid lipids, whose PUFA chains coat the rough protein surface, while the saturated chains face the raft environment, thus minimizing perturbations therein. Our findings suggest that the reduction of GPCR partitioning to their native ordered environments due to PUFA depletion might affect the function of these receptors in numerous neurodegenerative diseases, where the membrane PUFA levels in the brain are decreased. We hope that this work inspires experimental studies on the connection between membrane PUFA levels and GPCR signaling.

Molecular mechanism for inhibition of twinfilin by phosphoinositides.

Membrane phosphoinositides control organization and dynamics of the actin cytoskeleton by regulating the activities of several key actin-binding proteins. Twinfilin is an evolutionarily conserved protein that contributes to cytoskeletal dynamics by interacting with actin monomers, filaments, and the heterodimeric capping protein. Twinfilin also binds phosphoinositides, which inhibit its interactions with actin, but the underlying mechanism has remained unknown. Here, we show that the high-affinity binding site of twinfilin for phosphoinositides is located at the C-terminal tail region, whereas the two actin-depolymerizing factor (ADF)/cofilin-like ADF homology domains of twinfilin bind phosphoinositides only with low affinity. Mutagenesis and biochemical experiments combined with atomistic molecular dynamics simulations reveal that the C-terminal tail of twinfilin interacts with membranes through a multivalent electrostatic interaction with a preference toward phosphatidylinositol 3,5-bisphosphate (PI(3,5)P2), PI(4,5)P2, and PI(3,4,5)P3 This initial interaction places the actin-binding ADF homology domains of twinfilin in close proximity to the membrane and subsequently promotes their association with the membrane, thus leading to inhibition of the actin interactions. In support of this model, a twinfilin mutant lacking the C-terminal tail inhibits actin filament assembly in a phosphoinositide-insensitive manner. Our mutagenesis data also reveal that the phosphoinositide- and capping protein-binding sites overlap in the C-terminal tail of twinfilin, suggesting that phosphoinositide binding additionally inhibits the interactions of twinfilin with the heterodimeric capping protein. The results demonstrate that the conserved C-terminal tail of twinfilin is a multifunctional binding motif, which is crucial for interaction with the heterodimeric capping protein and for tethering twinfilin to phosphoinositide-rich membranes.

Concerted regulation of npc2 binding to endosomal/lysosomal membranes by bis(monoacylglycero)phosphate and sphingomyelin.

Niemann-Pick Protein C2 (npc2) is a small soluble protein critical for cholesterol transport within and from the lysosome and the late endosome. Intriguingly, npc2-mediated cholesterol transport has been shown to be modulated by lipids, yet the molecular mechanism of npc2-membrane interactions has remained elusive. Here, based on an extensive set of atomistic simulations and free energy calculations, we clarify the mechanism and energetics of npc2-membrane binding and characterize the roles of physiologically relevant key lipids associated with the binding process. Our results capture in atomistic detail two competitively favorable membrane binding orientations of npc2 with a low interconversion barrier. The first binding mode (Prone) places the cholesterol binding pocket in direct contact with the membrane and is characterized by membrane insertion of a loop (V59-M60-G61-I62-P63-V64-P65). This mode is associated with cholesterol uptake and release. On the other hand, the second mode (Supine) places the cholesterol binding pocket away from the membrane surface, but has overall higher membrane binding affinity. We determined that bis(monoacylglycero)phosphate (bmp) is specifically required for strong membrane binding in Prone mode, and that it cannot be substituted by other anionic lipids. Meanwhile, sphingomyelin counteracts bmp by hindering Prone mode without affecting Supine mode. Our results provide concrete evidence that lipids modulate npc2-mediated cholesterol transport either by favoring or disfavoring Prone mode and that they impose this by modulating the accessibility of bmp for interacting with npc2. Overall, we provide a mechanism by which npc2-mediated cholesterol transport is controlled by the membrane composition and how npc2-lipid interactions can regulate the transport rate.

Proton-coupled electron transfer and the role of water molecules in proton pumping by cytochrome c oxidase.

Molecular oxygen acts as the terminal electron sink in the respiratory chains of aerobic organisms. Cytochrome c oxidase in the inner membrane of mitochondria and the plasma membrane of bacteria catalyzes the reduction of oxygen to water, and couples the free energy of the reaction to proton pumping across the membrane. The proton-pumping activity contributes to the proton electrochemical gradient, which drives the synthesis of ATP. Based on kinetic experiments on the O-O bond splitting transition of the catalytic cycle (A → P(R)), it has been proposed that the electron transfer to the binuclear iron-copper center of O2 reduction initiates the proton pump mechanism. This key electron transfer event is coupled to an internal proton transfer from a conserved glutamic acid to the proton-loading site of the pump. However, the proton may instead be transferred to the binuclear center to complete the oxygen reduction chemistry, which would constitute a short-circuit. Based on atomistic molecular dynamics simulations of cytochrome c oxidase in an explicit membrane-solvent environment, complemented by related free-energy calculations, we propose that this short-circuit is effectively prevented by a redox-state-dependent organization of water molecules within the protein structure that gates the proton transfer pathway.

Transient formation of water-conducting states in membrane transporters.

Membrane transporters rely on highly coordinated structural transitions between major conformational states for their function, to prevent simultaneous access of the substrate binding site to both sides of the membrane--a mode of operation known as the alternating access model. Although this mechanism successfully accounts for the efficient exchange of the primary substrate across the membrane, accruing evidence on significant water transport and even uncoupled ion transport mediated by transporters has challenged the concept of perfect mechanical coupling and coordination of the gating mechanism in transporters, which might be expected from the alternating access model. Here, we present a large set of extended equilibrium molecular dynamics simulations performed on several classes of membrane transporters in different conformational states, to test the presence of the phenomenon in diverse transporter classes and to investigate the underlying molecular mechanism of water transport through membrane transporters. The simulations reveal spontaneous formation of transient water-conducting (channel-like) states allowing passive water diffusion through the lumen of the transporters. These channel-like states are permeable to water but occluded to substrate, thereby not hindering the uphill transport of the primary substrate, i.e., the alternating access model remains applicable to the substrate. The rise of such water-conducting states during the large-scale structural transitions of the transporter protein is indicative of imperfections in the coordinated closing and opening motions of the cytoplasmic and extracellular gates. We propose that the observed water-conducting states likely represent a universal phenomenon in membrane transporters, which is consistent with their reliance on large-scale motion for function.

Structure and permeation mechanism of a mammalian urea transporter.

As an adaptation to infrequent access to water, terrestrial mammals produce urine that is hyperosmotic to plasma. To prevent osmotic diuresis by the large quantity of urea generated by protein catabolism, the kidney epithelia contain facilitative urea transporters (UTs) that allow rapid equilibration between the urinary space and the hyperosmotic interstitium. Here we report the first X-ray crystal structure of a mammalian UT, UT-B, at a resolution of 2.36 Å. UT-B is a homotrimer and each protomer contains a urea conduction pore with a narrow selectivity filter. Structural analyses and molecular dynamics simulations showed that the selectivity filter has two urea binding sites separated by an approximately 5.0 kcal/mol energy barrier. Functional studies showed that the rate of urea conduction in UT-B is increased by hypoosmotic stress, and that the site of osmoregulation coincides with the location of the energy barrier.

Martini 3 OliGoÌ mers: A Scalable Approach for Multimers and Fibrils in GROMACS.

Martini 3 is a widely used coarse-grained simulation method for large-scale biomolecular simulations. It can be combined with a Go̅ model to realistically describe higher-order protein structures while allowing the folding and unfolding events. However, as of today, this method has largely been used only for individual monomers. In this article, we describe how the Go̅ model can be implemented within the framework of Martini 3 for a multimer system, taking into account both intramolecular and intermolecular interactions in an oligomeric protein system. We demonstrate the method by showing how it can be applied to both structural stability maintenance and assembly/disassembly of protein oligomers, using aquaporin tetramer, insulin dimer, and amyloid-ß fibril as examples. We find that addition of intermolecular Go̅ potentials stabilizes the quaternary structure of proteins. The strength of the Go̅ potentials can be tuned so that the internal fluctuations of proteins match the behavior of atomistic simulation models, however, the results also show that the use of too strong intermolecular Go̅ potentials weakens the chemical specificity of oligomerization. The Martini-Go̅ model presented here enables the use of Go̅ potentials in oligomeric molecular systems in a computationally efficient and parallelizable manner, especially in the case of homopolymers, where the number of identical protein monomers is high. This paves the way for coarse-grained simulations of large protein complexes, such as viral protein capsids and prion fibrils, in complex biological environments.

Visualizing functional motions of membrane transporters with molecular dynamics simulations.

Computational modeling and molecular simulation techniques have become an integral part of modern molecular research. Various areas of molecular sciences continue to benefit from, indeed rely on, the unparalleled spatial and temporal resolutions offered by these technologies, to provide a more complete picture of the molecular problems at hand. Because of the continuous development of more efficient algorithms harvesting ever-expanding computational resources, and the emergence of more advanced and novel theories and methodologies, the scope of computational studies has expanded significantly over the past decade, now including much larger molecular systems and far more complex molecular phenomena. Among the various computer modeling techniques, the application of molecular dynamics (MD) simulation and related techniques has particularly drawn attention in biomolecular research, because of the ability of the method to describe the dynamical nature of the molecular systems and thereby to provide a more realistic representation, which is often needed for understanding fundamental molecular properties. The method has proven to be remarkably successful in capturing molecular events and structural transitions highly relevant to the function and/or physicochemical properties of biomolecular systems. Herein, after a brief introduction to the method of MD, we use a number of membrane transport proteins studied in our laboratory as examples to showcase the scope and applicability of the method and its power in characterizing molecular motions of various magnitudes and time scales that are involved in the function of this important class of membrane proteins.

Movement of NH3 through the human urea transporter B: a new gas channel.

Aquaporins and Rh proteins can function as gas (CO2 and NH3) channels. The present study explores the urea, H2O, CO2, and NH3 permeability of the human urea transporter B (UT-B) (SLC14A1), expressed in Xenopus oocytes. We monitored urea uptake using [¹4C]urea and measured osmotic water permeability (Pf) using video microscopy. To obtain a semiquantitative measure of gas permeability, we used microelectrodes to record the maximum transient change in surface pH (ΔpHS) caused by exposing oocytes to 5% CO2/33 mM HCO3⁻ (pHS increase) or 0.5 mM NH3/NH4⁺ (pHS decrease). UT-B expression increased oocyte permeability to urea by >20-fold, and Pf by 8-fold vs. H2O-injected control oocytes. UT-B expression had no effect on the CO2-induced ΔpHS but doubled the NH3-induced ΔpHS. Phloretin reduced UT-B-dependent urea uptake (Jurea*) by 45%, Pf* by 50%, and (- ΔpHS*)NH3 by 70%. p-Chloromercuribenzene sulfonate reduced Jurea* by 25%, Pf* by 30%, and (ΔpHS*)NH3 by 100%. Molecular dynamics (MD) simulations of membrane-embedded models of UT-B identified the monomeric UT-B pores as the main conduction pathway for both H2O and NH3 and characterized the energetics associated with permeation of these species through the channel. Mutating each of two conserved threonines lining the monomeric urea pores reduced H2O and NH3 permeability. Our data confirm that UT-B has significant H2O permeability and for the first time demonstrate significant NH3 permeability. Thus the UTs become the third family of gas channels. Inhibitor and mutagenesis studies and results of MD simulations suggest that NH3 and H2O pass through the three monomeric urea channels in UT-B.

How Neuromembrane Lipids Modulate Membrane Proteins: Insights from G-Protein-Coupled Receptors (GPCRs) and Receptor Tyrosine Kinases (RTKs).

Lipids play a diverse and critical role in cellular processes in all tissues. The unique lipid composition of nerve membranes is particularly interesting because it contains, among other things, polyunsaturated lipids, such as docosahexaenoic acid, which the body only gets through the diet. The crucial role of lipids in neurological processes, especially in receptor-mediated cell signaling, is emphasized by the fact that in many neuropathological diseases there are significant deviations in the lipid composition of nerve membranes compared to healthy individuals. The lipid composition of neuromembranes can significantly affect the function of receptors by regulating the physical properties of the membrane or by affecting specific interactions between receptors and lipids. In addition, it is worth noting that the ligand-binding pocket of many receptors is located inside the cell membrane, due to which lipids can even modulate the binding of ligands to their receptors. These mechanisms highlight the importance of lipids in the regulation of membrane receptor activation and function. In this article, we focus on two major protein families: G-protein-coupled receptors (GPCRs) and receptor tyrosine kinases (RTKs) and discuss how lipids affect their function in neuronal membranes, elucidating the basic mechanisms underlying neuronal function and dysfunction.

Cholesterol esters form supercooled lipid droplets whose nucleation is facilitated by triacylglycerols.

Cellular cholesterol can be metabolized to its fatty acid esters, cholesteryl esters (CEs), to be stored in lipid droplets (LDs). With triacylglycerols (TGs), CEs represent the main neutral lipids in LDs. However, while TG melts at ~4 °C, CE melts at ~44 °C, raising the question of how CE-rich LDs form in cells. Here, we show that CE forms supercooled droplets when the CE concentration in LDs is above 20% to TG and, in particular, liquid-crystalline phases when the fraction of CEs is above 90% at 37 °C. In model bilayers, CEs condense and nucleate droplets when the CE/phospholipid ratio reaches over 10-15%. This concentration is reduced by TG pre-clusters in the membrane that thereby facilitate CE nucleation. Accordingly, blocking TG synthesis in cells is sufficient to strongly dampen CE LD nucleation. Finally, CE LDs emerged at seipins, which cluster and nucleate TG LDs in the ER. However, when TG synthesis is inhibited, similar numbers of LDs are generated in the presence and absence of seipin, suggesting that seipin controls CE LD formation via its TG clustering capacity. Our data point to a unique model whereby TG pre-clusters, favorable at seipins, catalyze the nucleation of CE LDs.

The GET insertase exhibits conformational plasticity and induces membrane thinning.

The eukaryotic guided entry of tail-anchored proteins (GET) pathway mediates the biogenesis of tail-anchored (TA) membrane proteins at the endoplasmic reticulum. In the cytosol, the Get3 chaperone captures the TA protein substrate and delivers it to the Get1/Get2 membrane protein complex (GET insertase), which then inserts the substrate via a membrane-embedded hydrophilic groove. Here, we present structures, atomistic simulations and functional data of human and Chaetomium thermophilum Get1/Get2/Get3. The core fold of the GET insertase is conserved throughout eukaryotes, whilst thinning of the lipid bilayer occurs in the vicinity of the hydrophilic groove to presumably lower the energetic barrier of membrane insertion. We show that the gating interaction between Get2 helix α3' and Get3 drives conformational changes in both Get3 and the Get1/Get2 membrane heterotetramer. Thus, we provide a framework to understand the conformational plasticity of the GET insertase and how it remodels its membrane environment to promote substrate insertion.

Zika virus prM protein contains cholesterol binding motifs required for virus entry and assembly.

For successful infection of host cells and virion production, enveloped viruses, including Zika virus (ZIKV), extensively rely on cellular lipids. However, how virus protein-lipid interactions contribute to the viral life cycle remains unclear. Here, we employ a chemo-proteomics approach with a bifunctional cholesterol probe and show that cholesterol is closely associated with the ZIKV structural protein prM. Bioinformatic analyses, reverse genetics alongside with photoaffinity labeling assays, and atomistic molecular dynamics simulations identified two functional cholesterol binding motifs within the prM transmembrane domain. Loss of prM-cholesterol association has a bipartite effect reducing ZIKV entry and leading to assembly defects. We propose a model in which membrane-resident M facilitates cholesterol-supported lipid exchange during endosomal entry and, together with cholesterol, creates a platform promoting virion assembly. In summary, we identify a bifunctional role of prM in the ZIKV life cycle by mediating viral entry and virus assembly in a cholesterol-dependent manner.

Psychedelics promote plasticity by directly binding to BDNF receptor TrkB.

Psychedelics produce fast and persistent antidepressant effects and induce neuroplasticity resembling the effects of clinically approved antidepressants. We recently reported that pharmacologically diverse antidepressants, including fluoxetine and ketamine, act by binding to TrkB, the receptor for BDNF. Here we show that lysergic acid diethylamide (LSD) and psilocin directly bind to TrkB with affinities 1,000-fold higher than those for other antidepressants, and that psychedelics and antidepressants bind to distinct but partially overlapping sites within the transmembrane domain of TrkB dimers. The effects of psychedelics on neurotrophic signaling, plasticity and antidepressant-like behavior in mice depend on TrkB binding and promotion of endogenous BDNF signaling but are independent of serotonin 2A receptor (5-HT2A) activation, whereas LSD-induced head twitching is dependent on 5-HT2A and independent of TrkB binding. Our data confirm TrkB as a common primary target for antidepressants and suggest that high-affinity TrkB positive allosteric modulators lacking 5-HT2A activity may retain the antidepressant potential of psychedelics without hallucinogenic effects.

Maturation of the SARS-CoV-2 virus is regulated by dimerization of its main protease.

SARS-CoV-2 main protease (Mpro) involved in COVID-19 is required for maturation of the virus and infection of host cells. The key question is how to block the activity of Mpro. By combining atomistic simulations with machine learning, we found that the enzyme regulates its own activity by a collective allosteric mechanism that involves dimerization and binding of a single substrate. At the core of the collective mechanism is the coupling between the catalytic site residues, H41 and C145, which direct the activity of Mpro dimer, and two salt bridges formed between R4 and E290 at the dimer interface. If these salt bridges are mutated, the activity of Mpro is blocked. The results suggest that dimerization of main proteases is a general mechanism to foster coronavirus proliferation, and propose a robust drug-based strategy that does not depend on the frequently mutating spike proteins at the viral envelope used to develop vaccines.