Myr-Arf1 conformational flexibility at the membrane surface sheds light on the interactions with ArfGAP ASAP1.

ADP-ribosylation factor 1 (Arf1) interacts with multiple cellular partners and membranes to regulate intracellular traffic, organelle structure and actin dynamics. Defining the dynamic conformational landscape of Arf1 in its active form, when bound to the membrane, is of high functional relevance and key to understanding how Arf1 can alter diverse cellular processes. Through concerted application of nuclear magnetic resonance (NMR), neutron reflectometry (NR) and molecular dynamics (MD) simulations, we show that, while Arf1 is anchored to the membrane through its N-terminal myristoylated amphipathic helix, the G domain explores a large conformational space, existing in a dynamic equilibrium between membrane-associated and membrane-distal conformations. These configurational dynamics expose different interfaces for interaction with effectors. Interaction with the Pleckstrin homology domain of ASAP1, an Arf-GTPase activating protein (ArfGAP), restricts motions of the G domain to lock it in what seems to be a conformation exposing functionally relevant regions.

Targeting lipid-protein interaction to treat Syk-mediated acute myeloid leukemia.

Membrane lipids control the cellular activity of kinases containing the Src homology 2 (SH2) domain through direct lipid-SH2 domain interactions. Here we report development of new nonlipidic small molecule inhibitors of the lipid-SH2 domain interaction that block the cellular activity of their host proteins. As a pilot study, we evaluated the efficacy of lipid-SH2 domain interaction inhibitors for spleen tyrosine kinase (Syk), which is implicated in hematopoietic malignancies, including acute myeloid leukemia (AML). An optimized inhibitor (WC36) specifically and potently suppressed oncogenic activities of Syk in AML cell lines and patient-derived AML cells. Unlike ATP-competitive Syk inhibitors, WC36 was refractory to de novo and acquired drug resistance due to its ability to block not only the Syk kinase activity, but also its noncatalytic scaffolding function that is linked to drug resistance. Collectively, our study shows that targeting lipid-protein interaction is a powerful approach to developing new small molecule drugs.

Cryo-EM structures of human ABCA7 provide insights into its phospholipid translocation mechanisms.

Phospholipid extrusion by ABC subfamily A (ABCA) exporters is central to cellular physiology, although the specifics of the underlying substrate interactions and transport mechanisms remain poorly resolved at the molecular level. Here we report cryo-EM structures of lipid-embedded human ABCA7 in an open state and in a nucleotide-bound, closed state at resolutions between 3.6 and 4.0 Å. The former reveals an ordered patch of bilayer lipids traversing the transmembrane domain (TMD), while the latter reveals a lipid-free, closed TMD with a small extracellular opening. These structures offer a structural framework for both substrate entry and exit from the ABCA7 TMD and highlight conserved rigid-body motions that underlie the associated conformational transitions. Combined with functional analysis and molecular dynamics (MD) simulations, our data also shed light on lipid partitioning into the ABCA7 TMD and localized membrane perturbations that underlie ABCA7 function and have broader implications for other ABCA family transporters.

Accelerating All-Atom Simulations and Gaining Mechanistic Understanding of Biophysical Systems through State Predictive Information Bottleneck.

An effective implementation of enhanced sampling algorithms for molecular dynamics simulations requires a priori knowledge of the approximate reaction coordinate describing the relevant mechanisms in the system. In this work, we focus on the recently developed artificial intelligence-based State Predictive Information Bottleneck (SPIB) approach and demonstrate how SPIB can learn such a reaction coordinate as a deep neural network even from undersampled trajectories. We exemplify its usefulness by achieving more than 40 times acceleration in simulating two model biophysical systems through well-tempered metadynamics performed by biasing along the SPIB-learned reaction coordinate. These include left- to right-handed chirality transitions in a synthetic helical peptide (Aib)9 and permeation of a small benzoic acid molecule through a synthetic, symmetric phospholipid bilayer. In addition to significantly accelerating the dynamics and achieving back and forth movement between different metastable states, the SPIB-based reaction coordinate gives mechanistic insights into the processes driving these two important problems.

Microscopic Characterization of the Chloride Permeation Pathway in the Human Excitatory Amino Acid Transporter 1 (EAAT1).

Excitatory amino acid transporters (EAATs) are glutamate transporters that belong to the solute carrier 1A (SLC1A) family. They couple glutamate transport to the cotransport of three sodium (Na+) ions and one proton (H+) and the counter-transport of one potassium (K+) ion. In addition to this coupled transport, binding of cotransported species to EAATs activates a thermodynamically uncoupled chloride (Cl-) conductance. Structures of SLC1A family members have revealed that these transporters use a twisting elevator mechanism of transport, where a mobile transport domain carries substrate and coupled ions across the membrane, while a static scaffold domain anchors the transporter in the membrane. We recently demonstrated that the uncoupled Cl- conductance is activated by the formation of an aqueous pore at the domain interface during the transport cycle in archaeal GltPh. However, a pathway for the uncoupled Cl- conductance has not been reported for the EAATs, and it is unclear if such a pathway is conserved. Here, we employ all-atom molecular dynamics (MD) simulations combined with enhanced sampling, free-energy calculations, and experimental mutagenesis to approximate large-scale conformational changes during the transport process and identified a Cl--conducting conformation in human EAAT1 (hEAAT1). Sampling the large-scale structural transitions in hEAAT1 allowed us to capture an intermediate conformation formed during the transport cycle with a continuous aqueous pore at the domain interface. The free-energy calculations performed for the conduction of Cl- and Na+ ions through the captured conformation highlight the presence of two hydrophobic gates that control low-barrier movement of Cl- through the aqueous pathway. Overall, our findings provide insights into the mechanism by which a human neurotransmitter transporter supports functional duality of active transport and passive Cl- permeation and confirm the commonality of this mechanism in different members of the SLC1A family.

A Companion Guide to the String Method with Swarms of Trajectories: Characterization, Performance, and Pitfalls.

The string method with swarms of trajectories (SMwST) is an algorithm that identifies a physically meaningful transition pathway─a one-dimensional curve, embedded within a high-dimensional space of selected collective variables. The SMwST algorithm leans on a series of short, unbiased molecular dynamics simulations spawned at different locations of the discretized path, from whence an average dynamic drift is determined to evolve the string toward an optimal pathway. However conceptually simple in both its theoretical formulation and practical implementation, the SMwST algorithm is computationally intensive and requires a careful choice of parameters for optimal cost-effectiveness in applications to challenging problems in chemistry and biology. In this contribution, the SMwST algorithm is presented in a self-contained manner, discussing with a critical eye its theoretical underpinnings, applicability, inherent limitations, and use in the context of path-following free-energy calculations and their possible extension to kinetics modeling. Through multiple simulations of a prototypical polypeptide, combining the search of the transition pathway and the computation of the potential of mean force along it, several practical aspects of the methodology are examined with the objective of optimizing the computational effort, yet without sacrificing accuracy. In light of the results reported here, we propose some general guidelines aimed at improving the efficiency and reliability of the computed pathways and free-energy profiles underlying the conformational transitions at hand.

Identification of structural transitions in bacterial fatty acid binding proteins that permit ligand entry and exit at membranes.

Fatty acid (FA) transfer proteins extract FA from membranes and sequester them to facilitate their movement through the cytosol. Detailed structural information is available for these soluble protein-FA complexes, but the structure of the protein conformation responsible for FA exchange at the membrane is unknown. Staphylococcus aureus FakB1 is a prototypical bacterial FA transfer protein that binds palmitate within a narrow, buried tunnel. Here, we define the conformational change from a "closed" FakB1 state to an "open" state that associates with the membrane and provides a path for entry and egress of the FA. Using NMR spectroscopy, we identified a conformationally flexible dynamic region in FakB1, and X-ray crystallography of FakB1 mutants captured the conformation of the open state. In addition, molecular dynamics simulations show that the new amphipathic α-helix formed in the open state inserts below the phosphate plane of the bilayer to create a diffusion channel for the hydrophobic FA tail to access the hydrocarbon core and place the carboxyl group at the phosphate layer. The membrane binding and catalytic properties of site-directed mutants were consistent with the proposed membrane docked structure predicted by our molecular dynamics simulations. Finally, the structure of the bilayer-associated conformation of FakB1 has local similarities with mammalian FA binding proteins and provides a conceptual framework for how these proteins interact with the membrane to create a diffusion channel from the FA location in the bilayer to the protein interior.

Ataxia-linked SLC1A3 mutations alter EAAT1 chloride channel activity and glial regulation of CNS function.

Glutamate is the predominant excitatory neurotransmitter in the mammalian central nervous system (CNS). Excitatory amino acid transporters (EAATs) regulate extracellular glutamate by transporting it into cells, mostly glia, to terminate neurotransmission and to avoid neurotoxicity. EAATs are also chloride (Cl-) channels, but the physiological role of Cl- conductance through EAATs is poorly understood. Mutations of human EAAT1 (hEAAT1) have been identified in patients with episodic ataxia type 6 (EA6). One mutation showed increased Cl- channel activity and decreased glutamate transport, but the relative contributions of each function of hEAAT1 to mechanisms underlying the pathology of EA6 remain unclear. Here we investigated the effects of 5 additional EA6-related mutations on hEAAT1 function in Xenopus laevis oocytes, and on CNS function in a Drosophila melanogaster model of locomotor behavior. Our results indicate that mutations resulting in decreased hEAAT1 Cl- channel activity but with functional glutamate transport can also contribute to the pathology of EA6, highlighting the importance of Cl- homeostasis in glial cells for proper CNS function. We also identified what we believe is a novel mechanism involving an ectopic sodium (Na+) leak conductance in glial cells. Together, these results strongly support the idea that EA6 is primarily an ion channelopathy of CNS glia.

PIP2-dependent coupling of voltage sensor and pore domains in Kv7.2 channel.

Phosphatidylinositol-4,5-bisphosphate (PIP2) is a signaling lipid which regulates voltage-gated Kv7/KCNQ potassium channels. Altered PIP2 sensitivity of neuronal Kv7.2 channel is involved in KCNQ2 epileptic encephalopathy. However, the molecular action of PIP2 on Kv7.2 gating remains largely elusive. Here, we use molecular dynamics simulations and electrophysiology to characterize PIP2 binding sites in a human Kv7.2 channel. In the closed state, PIP2 localizes to the periphery of the voltage-sensing domain (VSD). In the open state, PIP2 binds to 4 distinct interfaces formed by the cytoplasmic ends of the VSD, the gate, intracellular helices A and B and their linkers. PIP2 binding induces bilayer-interacting conformation of helices A and B and the correlated motion of the VSD and the pore domain, whereas charge-neutralizing mutations block this coupling and reduce PIP2 sensitivity of Kv7.2 channels by disrupting PIP2 binding. These findings reveal the allosteric role of PIP2 in Kv7.2 channel activation.

Active participation of membrane lipids in inhibition of multidrug transporter P-glycoprotein.

P-glycoprotein (Pgp) is a major efflux pump in humans, overexpressed in a variety of cancers and associated with the development of multi-drug resistance. Allosteric modulation by various ligands (e.g., transport substrates, inhibitors, and ATP) has been biochemically shown to directly influence structural dynamics, and thereby, the function of Pgp. However, the molecular details of such effects, particularly with respect to the role and involvement of the surrounding lipids, are not well established. Here, we employ all-atom molecular dynamics (MD) simulations to study the conformational landscape of Pgp in the presence of a high-affinity, third-generation inhibitor, tariquidar, in comparison to the nucleotide-free (APO) and the ATP-bound states, in order to characterize the mechanical effects of the inhibitor that might be of relevance to its blocking mechanism of Pgp. Simulations in a multi-component lipid bilayer show a dynamic equilibrium between open(er) and more closed inward-facing (IF) conformations in the APO state, with binding of ATP shifting the equilibrium towards conformations more prone to ATP hydrolysis and subsequent events in the transport cycle. In the presence of the inhibitor bound to the drug-binding pocket within the transmembrane domain (TMD), Pgp samples more open IF conformations, and the nucleotide binding domains (NBDs) become highly dynamic. Interestingly, and reproduced in multiple independent simulations, the inhibitor is observed to facilitate recruitment of lipid molecules into the Pgp lumen through the two proposed drug-entry portals, where the lipid head groups from the cytoplasmic leaflet penetrate into and, in some cases, translocate inside the TMD, while the lipid tails remain extended into the bulk lipid environment. These "wedge" lipids likely enhance the inhibitor-induced conformational restriction of the TMD leading to the differential modulation of coupling pathways observed with the NBDs downstream. We suggest a novel inhibitory mechanism for tariquidar, and potentially for related third-generation Pgp inhibitors, where lipids are seen to enhance the inhibitory role in the catalytic cycle of membrane transporters.

Evaluation and Selection of Video Stabilization Techniques for UAV-Based Active Infrared Thermography Application.

Unmanned Aerial Vehicles (UAVs) that can fly around an aircraft carrying several sensors, e.g., thermal and optical cameras, to inspect the parts of interest without removing them can have significant impact in reducing inspection time and cost. One of the main challenges in the UAV based active InfraRed Thermography (IRT) inspection is the UAV's unexpected motions. Since active thermography is mainly concerned with the analysis of thermal sequences, unexpected motions can disturb the thermal profiling and cause data misinterpretation especially for providing an automated process pipeline of such inspections. Additionally, in the scenarios where post-analysis is intended to be applied by an inspector, the UAV's unexpected motions can increase the risk of human error, data misinterpretation, and incorrect characterization of possible defects. Therefore, post-processing is required to minimize/eliminate such undesired motions using digital video stabilization techniques. There are number of video stabilization algorithms that are readily available; however, selecting the best suited one is also challenging. Therefore, this paper evaluates video stabilization algorithms to minimize/mitigate undesired UAV motion and proposes a simple method to find the best suited stabilization algorithm as a fundamental first step towards a fully operational UAV-IRT inspection system.

Glutamate transporters have a chloride channel with two hydrophobic gates.

Glutamate is the most abundant excitatory neurotransmitter in the central nervous system, and its precise control is vital to maintain normal brain function and to prevent excitotoxicity1. The removal of extracellular glutamate is achieved by plasma-membrane-bound transporters, which couple glutamate transport to sodium, potassium and pH gradients using an elevator mechanism2-5. Glutamate transporters also conduct chloride ions by means of a channel-like process that is thermodynamically uncoupled from transport6-8. However, the molecular mechanisms that enable these dual-function transporters to carry out two seemingly contradictory roles are unknown. Here we report the cryo-electron microscopy structure of a glutamate transporter homologue in an open-channel state, which reveals an aqueous cavity that is formed during the glutamate transport cycle. The functional properties of this cavity, combined with molecular dynamics simulations, reveal it to be an aqueous-accessible chloride permeation pathway that is gated by two hydrophobic regions and is conserved across mammalian and archaeal glutamate transporters. Our findings provide insight into the mechanism by which glutamate transporters support their dual function, and add information that will assist in mapping the complete transport cycle shared by the solute carrier 1A transporter family.

Confronting pitfalls of AI-augmented molecular dynamics using statistical physics.

Artificial intelligence (AI)-based approaches have had indubitable impact across the sciences through the ability to extract relevant information from raw data. Recently, AI has also found use in enhancing the efficiency of molecular simulations, wherein AI derived slow modes are used to accelerate the simulation in targeted ways. However, while typical fields where AI is used are characterized by a plethora of data, molecular simulations, per construction, suffer from limited sampling and thus limited data. As such, the use of AI in molecular simulations can suffer from a dangerous situation where the AI-optimization could get stuck in spurious regimes, leading to incorrect characterization of the reaction coordinate (RC) for the problem at hand. When such an incorrect RC is then used to perform additional simulations, one could start to deviate progressively from the ground truth. To deal with this problem of spurious AI-solutions, here, we report a novel and automated algorithm using ideas from statistical mechanics. It is based on the notion that a more reliable AI-solution will be one that maximizes the timescale separation between slow and fast processes. To learn this timescale separation even from limited data, we use a maximum caliber-based framework. We show the applicability of this automatic protocol for three classic benchmark problems, namely, the conformational dynamics of a model peptide, ligand-unbinding from a protein, and folding/unfolding energy landscape of the C-terminal domain of protein G. We believe that our work will lead to increased and robust use of trustworthy AI in molecular simulations of complex systems.

Hydrogen-deuterium exchange mass spectrometry captures distinct dynamics upon substrate and inhibitor binding to a transporter.

Proton-coupled transporters use transmembrane proton gradients to power active transport of nutrients inside the cell. High-resolution structures often fail to capture the coupling between proton and ligand binding, and conformational changes associated with transport. We combine HDX-MS with mutagenesis and MD simulations to dissect the molecular mechanism of the prototypical transporter XylE. We show that protonation of a conserved aspartate triggers conformational transition from outward-facing to inward-facing state. This transition only occurs in the presence of substrate xylose, while the inhibitor glucose locks the transporter in the outward-facing state. MD simulations corroborate the experiments by showing that only the combination of protonation and xylose binding, and not glucose, sets up the transporter for conformational switch. Overall, we demonstrate the unique ability of HDX-MS to distinguish between the conformational dynamics of inhibitor and substrate binding, and show that a specific allosteric coupling between substrate binding and protonation is a key step to initiate transport.

An Investigation of Rotary Drone HERM Line Spectrum under Manoeuvering Conditions.

Detecting and identifying drones is of great interest due to the proliferation of highly manoeuverable drones with on-board sensors of increasing sensing capabilities. In this paper, we investigate the use of radars for tackling this problem. In particular, we focus on the problem of detecting rotary drones and distinguishing between single-propeller and multi-propeller drones using a micro-Doppler analysis. Two different radars were used, an ultra wideband (UWB) continuous wave (CW) C-band radar and an automotive frequency modulated continuous wave (FMCW) W-band radar, to collect micro-Doppler signatures of the drones. By taking a closer look at HElicopter Rotor Modulation (HERM) lines, the spool and chopping lines are identified for the first time in the context of drones to determine the number of propeller blades. Furthermore, a new multi-frequency analysis method using HERM lines is developed, which allows the detection of propeller rotation rates (spool and chopping frequencies) of single and multi-propeller drones. Therefore, the presented method is a promising technique to aid in the classification of drones.

Membrane surface recognition by the ASAP1 PH domain and consequences for interactions with the small GTPase Arf1.

Adenosine diphosphate-ribosylation factor (Arf) guanosine triphosphatase-activating proteins (GAPs) are enzymes that need to bind to membranes to catalyze the hydrolysis of guanosine triphosphate (GTP) bound to the small GTP-binding protein Arf. Binding of the pleckstrin homology (PH) domain of the ArfGAP With SH3 domain, ankyrin repeat and PH domain 1 (ASAP1) to membranes containing phosphatidylinositol 4,5-bisphosphate [PI(4,5)P2] is key for maximum GTP hydrolysis but not fully understood. By combining nuclear magnetic resonance, neutron reflectometry, and molecular dynamics simulation, we show that binding of multiple PI(4,5)P2 molecules to the ASAP1 PH domain (i) triggers a functionally relevant allosteric conformational switch and (ii) maintains the PH domain in a well-defined orientation, allowing critical contacts with an Arf1 mimic to occur. Our model provides a framework to understand how binding of the ASAP1 PH domain to PI(4,5)P2 at the membrane may play a role in the regulation of ASAP1.

Identifying mutation hotspots reveals pathogenetic mechanisms of KCNQ2 epileptic encephalopathy.

Kv7 channels are enriched at the axonal plasma membrane where their voltage-dependent potassium currents suppress neuronal excitability. Mutations in Kv7.2 and Kv7.3 subunits cause epileptic encephalopathy (EE), yet the underlying pathogenetic mechanism is unclear. Here, we used novel statistical algorithms and structural modeling to identify EE mutation hotspots in key functional domains of Kv7.2 including voltage sensing S4, the pore loop and S6 in the pore domain, and intracellular calmodulin-binding helix B and helix B-C linker. Characterization of selected EE mutations from these hotspots revealed that L203P at S4 induces a large depolarizing shift in voltage dependence of Kv7.2 channels and L268F at the pore decreases their current densities. While L268F severely reduces expression of heteromeric channels in hippocampal neurons without affecting internalization, K552T and R553L mutations at distal helix B decrease calmodulin-binding and axonal enrichment. Importantly, L268F, K552T, and R553L mutations disrupt current potentiation by increasing phosphatidylinositol 4,5-bisphosphate (PIP2), and our molecular dynamics simulation suggests PIP2 interaction with these residues. Together, these findings demonstrate that each EE variant causes a unique combination of defects in Kv7 channel function and neuronal expression, and suggest a critical need for both prediction algorithms and experimental interrogations to understand pathophysiology of Kv7-associated EE.

Computational Dissection of Membrane Transport at a Microscopic Level.

Membrane transporters are key gatekeeper proteins at cellular membranes that closely control the traffic of materials. Their function relies on structural rearrangements of varying degrees that facilitate substrate translocation across the membrane. Characterizing these functionally important molecular events at a microscopic level is key to our understanding of membrane transport, yet challenging to achieve experimentally. Recent advances in simulation technology and computing power have rendered molecular dynamics (MD) simulation a powerful biophysical tool to investigate a wide range of dynamical events spanning multiple spatial and temporal scales. Here, we review recent studies of diverse membrane transporters using computational methods, with an emphasis on highlighting the technical challenges, key lessons learned, and new opportunities to illuminate transporter structure and function.

Microscopic Characterization of GRP1 PH Domain Interaction with Anionic Membranes.

The pleckstrin homology (PH) domain of general receptor for phosphoionositides 1 (GRP1-PHD) binds specifically to phosphatidylinositol (3,4,5)-triphosphate (PIP3 ), and acts as a second messenger. Using an extensive array of molecular dynamics (MD) simulations employing highly mobile membrane mimetic (HMMM) model as well as complementary full membrane simulations, we capture differentiable binding and dynamics of GRP1-PHD in the presence of membranes containing PC, PS, and PIP3 lipids in varying compositions. While GRP1-PHD forms only transient interactions with pure PC membranes, incorporation of anionic lipids resulted in stable membrane-bound configurations. We report the first observation of two distinct PIP3 binding modes on GRP1-PHD, involving PIP3 interactions at a "canonical" and at an "alternate" site, suggesting the possibility of simultaneous binding of multiple anionic lipids. The full membrane simulations confirmed the stability of the membrane bound pose of GRP1-PHD as captured from our HMMM membrane binding simulations. By performing additional steered membrane unbinding simulations and calculating nonequilibrium work associated with the process, as well as metadynamics simulations, on the protein bound to full membranes, allowing for more quantitative examination of the binding strength of the GRP1-PHD to the membrane, we demonstrate that along with the bound PIP3 , surrounding anionic PS lipids increase the energetic cost of unbinding of GRP1-PHD from the canonical mode, causing them to dissociate more slowly than the alternate mode. Our results demonstrate that concurrent binding of multiple anionic lipids by GRP1-PHD contributes to its membrane affinity, which in turn control its signaling activity. © 2019 Wiley Periodicals, Inc.

Characterization of Lipid-Protein Interactions and Lipid-Mediated Modulation of Membrane Protein Function through Molecular Simulation.

The cellular membrane constitutes one of the most fundamental compartments of a living cell, where key processes such as selective transport of material and exchange of information between the cell and its environment are mediated by proteins that are closely associated with the membrane. The heterogeneity of lipid composition of biological membranes and the effect of lipid molecules on the structure, dynamics, and function of membrane proteins are now widely recognized. Characterization of these functionally important lipid-protein interactions with experimental techniques is however still prohibitively challenging. Molecular dynamics (MD) simulations offer a powerful complementary approach with sufficient temporal and spatial resolutions to gain atomic-level structural information and energetics on lipid-protein interactions. In this review, we aim to provide a broad survey of MD simulations focusing on exploring lipid-protein interactions and characterizing lipid-modulated protein structure and dynamics that have been successful in providing novel insight into the mechanism of membrane protein function.