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.

Martini 3: a general purpose force field for coarse-grained molecular dynamics.

The coarse-grained Martini force field is widely used in biomolecular simulations. Here we present the refined model, Martini 3 ( http://cgmartini.nl ), with an improved interaction balance, new bead types and expanded ability to include specific interactions representing, for example, hydrogen bonding and electronic polarizability. The updated model allows more accurate predictions of molecular packing and interactions in general, which is exemplified with a vast and diverse set of applications, ranging from oil/water partitioning and miscibility data to complex molecular systems, involving protein-protein and protein-lipid interactions and material science applications as ionic liquids and aedamers.

Seipin traps triacylglycerols to facilitate their nanoscale clustering in the endoplasmic reticulum membrane.

Seipin is a disk-like oligomeric endoplasmic reticulum (ER) protein important for lipid droplet (LD) biogenesis and triacylglycerol (TAG) delivery to growing LDs. Here we show through biomolecular simulations bridged to experiments that seipin can trap TAGs in the ER bilayer via the luminal hydrophobic helices of the protomers delineating the inner opening of the seipin disk. This promotes the nanoscale sequestration of TAGs at a concentration that by itself is insufficient to induce TAG clustering in a lipid membrane. We identify Ser166 in the α3 helix as a favored TAG occupancy site and show that mutating it compromises the ability of seipin complexes to sequester TAG in silico and to promote TAG transfer to LDs in cells. While the S166D-seipin mutant colocalizes poorly with promethin, the association of nascent wild-type seipin complexes with promethin is promoted by TAGs. Together, these results suggest that seipin traps TAGs via its luminal hydrophobic helices, serving as a catalyst for seeding the TAG cluster from dissolved monomers inside the seipin ring, thereby generating a favorable promethin binding interface.

Exposure to Aldehyde Cherry e-Liquid Flavoring and Its Vaping Byproduct Disrupt Pulmonary Surfactant Biophysical Function.

Over the past decade, there has been a significant rise in the use of vaping devices, particularly among adolescents, raising concerns for effects on respiratory health. Pressingly, many recent vaping-related lung injuries are unexplained by current knowledge, and the overall implications of vaping for respiratory health are poorly understood. This study investigates the effect of hydrophobic vaping liquid chemicals on the pulmonary surfactant biophysical function. We focus on the commonly used flavoring benzaldehyde and its vaping byproduct, benzaldehyde propylene glycol acetal. The study involves rigorous testing of the surfactant biophysical function in Langmuir trough and constrained sessile drop surfactometer experiments with both protein-free synthetic surfactant and hydrophobic protein-containing clinical surfactant models. The study reveals that exposure to these vaping chemicals significantly interferes with the synthetic and clinical surfactant biophysical function. Further atomistic simulations reveal preferential interactions with SP-B and SP-C surfactant proteins. Additionally, data show surfactant lipid-vaping chemical interactions and suggest significant transfer of vaping chemicals to the experimental subphase, indicating a toxicological mechanism for the alveolar epithelium. Our study, therefore, reveals novel mechanisms for the inhalational toxicity of vaping. This highlights the need to reassess the safety of vaping liquids for respiratory health, particularly the use of aldehyde chemicals as vaping flavorings.

Surfactant Proteins SP-B and SP-C in Pulmonary Surfactant Monolayers: Physical Properties Controlled by Specific Protein-Lipid Interactions.

The lining of the alveoli is covered by pulmonary surfactant, a complex mixture of surface-active lipids and proteins that enables efficient gas exchange between inhaled air and the circulation. Despite decades of advancements in the study of the pulmonary surfactant, the molecular scale behavior of the surfactant and the inherent role of the number of different lipids and proteins in surfactant behavior are not fully understood. The most important proteins in this complex system are the surfactant proteins SP-B and SP-C. Given this, in this work we performed nonequilibrium all-atom molecular dynamics simulations to study the interplay of SP-B and SP-C with multicomponent lipid monolayers mimicking the pulmonary surfactant in composition. The simulations were complemented by z-scan fluorescence correlation spectroscopy and atomic force microscopy measurements. Our state-of-the-art simulation model reproduces experimental pressure-area isotherms and lateral diffusion coefficients. In agreement with previous research, the inclusion of either SP-B and SP-C increases surface pressure, and our simulations provide a molecular scale explanation for this effect: The proteins display preferential lipid interactions with phosphatidylglycerol, they reside predominantly in the lipid acyl chain region, and they partition into the liquid expanded phase or even induce it in an otherwise packed monolayer. The latter effect is also visible in our atomic force microscopy images. The research done contributes to a better understanding of the roles of specific lipids and proteins in surfactant function, thus helping to develop better synthetic products for surfactant replacement therapy used in the treatment of many fatal lung-related injuries and diseases.

Unbreaking Assemblies in Molecular Simulations with Periodic Boundaries.

In molecular simulations, periodic boundary conditions are typically used to avoid surface effects occurring at the boundaries of the simulation box. A consequence of this is that molecules and assemblies may appear split over the boundaries. Broken molecular assemblies make it difficult to interpret, analyze, and visualize molecular simulation data. We present a general and fast algorithm that repairs molecular assemblies that are broken due to periodic boundary conditions. The open source method presented here, MDVWhole, works for all translation-only crystallographic periodic boundary conditions. The method consumes little memory and can fix the visualization of the assembly of millions of particles in a few seconds. Thus, it is suitable for processing both single simulation frames and long trajectories with millions of points.

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.

The F1 loop of the talin head domain acts as a gatekeeper in integrin activation and clustering.

Integrin activation and clustering by talin are early steps of cell adhesion. Membrane-bound talin head domain and kindlin bind to the ß integrin cytoplasmic tail, cooperating to activate the heterodimeric integrin, and the talin head domain induces integrin clustering in the presence of Mn2+ Here we show that kindlin-1 can replace Mn2+ to mediate ß3 integrin clustering induced by the talin head, but not that induced by the F2-F3 fragment of talin. Integrin clustering mediated by kindlin-1 and the talin head was lost upon deletion of the flexible loop within the talin head F1 subdomain. Further mutagenesis identified hydrophobic and acidic motifs in the F1 loop responsible for ß3 integrin clustering. Modeling, computational and cysteine crosslinking studies showed direct and catalytic interactions of the acidic F1 loop motif with the juxtamembrane domains of α- and ß3-integrins, in order to activate the ß3 integrin heterodimer, further detailing the mechanism by which the talin-kindlin complex activates and clusters integrins. Moreover, the F1 loop interaction with the ß3 integrin tail required the newly identified compact FERM fold of the talin head, which positions the F1 loop next to the inner membrane clasp of the talin-bound integrin heterodimer.This article has an associated First Person interview with the first author of the paper.

Author Correction: An efficient auxin-inducible degron system with low basal degradation in human cells.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

An efficient auxin-inducible degron system with low basal degradation in human cells.

Auxin-inducible degron technology allows rapid and controlled protein depletion. However, basal degradation without auxin and inefficient auxin-inducible depletion have limited its utility. We have identified a potent auxin-inducible degron system composed of auxin receptor F-box protein AtAFB2 and short degron miniIAA7. The system showed minimal basal degradation and enabled rapid auxin-inducible depletion of endogenous human transmembrane, cytoplasmic and nuclear proteins in 1 h with robust functional phenotypes.

Cryo-EM, X-ray diffraction, and atomistic simulations reveal determinants for the formation of a supramolecular myelin-like proteolipid lattice.

Myelin protein P2 is a peripheral membrane protein of the fatty acid-binding protein family that functions in the formation and maintenance of the peripheral nerve myelin sheath. Several P2 gene mutations cause human Charcot-Marie-Tooth neuropathy, but the mature myelin sheath assembly mechanism is unclear. Here, cryo-EM of myelin-like proteolipid multilayers revealed an ordered three-dimensional (3D) lattice of P2 molecules between stacked lipid bilayers, visualizing supramolecular assembly at the myelin major dense line. The data disclosed that a single P2 layer is inserted between two bilayers in a tight intermembrane space of ∼3 nm, implying direct interactions between P2 and two membrane surfaces. X-ray diffraction from P2-stacked bicelle multilayers revealed lateral protein organization, and surface mutagenesis of P2 coupled with structure-function experiments revealed a role for both the portal region of P2 and its opposite face in membrane interactions. Atomistic molecular dynamics simulations of P2 on model membrane surfaces suggested that Arg-88 is critical for P2-membrane interactions, in addition to the helical lid domain. Negatively charged lipid headgroups stably anchored P2 on the myelin-like bilayer surface. Membrane binding may be accompanied by opening of the P2 ß-barrel structure and ligand exchange with the apposing bilayer. Our results provide an unprecedented view into an ordered, multilayered biomolecular membrane system induced by the presence of a peripheral membrane protein from human myelin. This is an important step toward deciphering the 3D assembly of a mature myelin sheath at the molecular level.

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.

Redox-coupled quinone dynamics in the respiratory complex I.

Complex I couples the free energy released from quinone (Q) reduction to pump protons across the biological membrane in the respiratory chains of mitochondria and many bacteria. The Q reduction site is separated by a large distance from the proton-pumping membrane domain. To address the molecular mechanism of this long-range proton-electron coupling, we perform here full atomistic molecular dynamics simulations, free energy calculations, and continuum electrostatics calculations on complex I from Thermus thermophilus We show that the dynamics of Q is redox-state-dependent, and that quinol, QH2, moves out of its reduction site and into a site in the Q tunnel that is occupied by a Q analog in a crystal structure of Yarrowia lipolytica We also identify a second Q-binding site near the opening of the Q tunnel in the membrane domain, where the Q headgroup forms strong interactions with a cluster of aromatic and charged residues, while the Q tail resides in the lipid membrane. We estimate the effective diffusion coefficient of Q in the tunnel, and in turn the characteristic time for Q to reach the active site and for QH2 to escape to the membrane. Our simulations show that Q moves along the Q tunnel in a redox-state-dependent manner, with distinct binding sites formed by conserved residue clusters. The motion of Q to these binding sites is proposed to be coupled to the proton-pumping machinery in complex I.

Tail-Oxidized Cholesterol Enhances Membrane Permeability for Small Solutes.

Cholesterol renders mammalian cell membranes more compact by reducing the amount of voids in the membrane structure. Because of this, cholesterol is known to regulate the ability of cell membranes to prevent the permeation of water and water-soluble molecules through the membranes. Meanwhile, it is also known that even seemingly tiny modifications in the chemical structure of cholesterol can lead to notable changes in membrane properties. The question is, how significantly do these small changes in cholesterol structure affect the permeability barrier function of cell membranes? In this work, we applied fluorescence methods as well as atomistic molecular dynamics simulations to characterize changes in lipid membrane permeability induced by cholesterol oxidation. The studied 7ß-hydroxycholesterol (7ß-OH-chol) and 27-hydroxycholesterol (27-OH-chol) represent two distinct groups of oxysterols, namely, ring- and tail-oxidized cholesterols, respectively. Our previous research showed that the oxidation of the cholesterol tail has only a marginal effect on the structure of a lipid bilayer; however, oxidation was found to disturb membrane dynamics by introducing a mechanism that allows sterol molecules to move rapidly back and forth across the membrane-bobbing. Herein, we show that bobbing of 27-OH-chol accelerates fluorescence quenching of NBD-lipid probes in the inner leaflet of liposomes by dithionite added to the liposomal suspension. Systematic experiments using fluorescence quenching spectroscopy and microscopy led to the conclusion that the presence of 27-OH-chol increases membrane permeability to the dithionite anion. Atomistic molecular dynamics simulations demonstrated that 27-OH-chol also facilitates water transport across the membrane. The results support the view that oxysterol bobbing gives rise to successive perturbations to the hydrophobic core of the membrane, and these perturbations promote the permeation of water and small water-soluble molecules through a lipid bilayer. The observed impairment of permeability can have important consequences for eukaryotic organisms. The effects described for 27-OH-chol were not observed for 7ß-OH-chol which represents ring-oxidized sterols.

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.

Mechanistic principles underlying regulation of the actin cytoskeleton by phosphoinositides.

The actin cytoskeleton powers membrane deformation during many cellular processes, such as migration, morphogenesis, and endocytosis. Membrane phosphoinositides, especially phosphatidylinositol 4,5-bisphosphate [PI(4,5)P2], regulate the activities of many actin-binding proteins (ABPs), including profilin, cofilin, Dia2, N-WASP, ezrin, and moesin, but the underlying molecular mechanisms have remained elusive. Moreover, because of a lack of available methodology, the dynamics of membrane interactions have not been experimentally determined for any ABP. Here, we applied a combination of biochemical assays, photobleaching/activation approaches, and atomistic molecular dynamics simulations to uncover the molecular principles by which ABPs interact with phosphoinositide-rich membranes. We show that, despite using different domains for lipid binding, these proteins associate with membranes through similar multivalent electrostatic interactions, without specific binding pockets or penetration into the lipid bilayer. Strikingly, our experiments reveal that these proteins display enormous differences in the dynamics of membrane interactions and in the ranges of phosphoinositide densities that they sense. Profilin and cofilin display transient, low-affinity interactions with phosphoinositide-rich membranes, whereas F-actin assembly factors Dia2 and N-WASP reside on phosphoinositide-rich membranes for longer periods to perform their functions. Ezrin and moesin, which link the actin cytoskeleton to the plasma membrane, bind membranes with very high affinity and slow dissociation dynamics. Unlike profilin, cofilin, Dia2, and N-WASP, they do not require high "stimulus-responsive" phosphoinositide density for membrane binding. Moreover, ezrin can limit the lateral diffusion of PI(4,5)P2 along the lipid bilayer. Together, these findings demonstrate that membrane-interaction mechanisms of ABPs evolved to precisely fulfill their specific functions in cytoskeletal dynamics.

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.

The Role of Temperature and Lipid Charge on Intake/Uptake of Cationic Gold Nanoparticles into Lipid Bilayers.

Understanding the molecular mechanisms governing nanoparticle-membrane interactions is of prime importance for drug delivery and biomedical applications. Neutron reflectometry (NR) experiments are combined with atomistic and coarse-grained molecular dynamics (MD) simulations to study the interaction between cationic gold nanoparticles (AuNPs) and model lipid membranes composed of a mixture of zwitterionic di-stearoyl-phosphatidylcholine (DSPC) and anionic di-stearoyl-phosphatidylglycerol (DSPG). MD simulations show that the interaction between AuNPs and a pure DSPC lipid bilayer is modulated by a free energy barrier. This can be overcome by increasing temperature, which promotes an irreversible AuNP incorporation into the lipid bilayer. NR experiments confirm the encapsulation of the AuNPs within the lipid bilayer at temperatures around 55 °C. In contrast, the AuNP adsorption is weak and impaired by heating for a DSPC-DSPG (3:1) lipid bilayer. These results demonstrate that both the lipid charge and the temperature play pivotal roles in AuNP-membrane interactions. Furthermore, NR experiments indicate that the (negative) DSPG lipids are associated with lipid extraction upon AuNP adsorption, which is confirmed by coarse-grained MD simulations as a lipid-crawling effect driving further AuNP aggregation. Overall, the obtained detailed molecular view of the interaction mechanisms sheds light on AuNP incorporation and membrane destabilization.