MemProtMD: Automated Insertion of Membrane Protein Structures into Explicit Lipid Membranes

Three-Dimensional Interaction Homology: Deconstructing Residue-Residue and Residue-Lipid Interactions in Membrane Proteins

A method is described to deconstruct the network of hydropathic interactions within and between a protein's sidechain and its environment into residue-based three-dimensional maps. These maps encode favorable and unfavorable hydrophobic and polar interactions, in terms of spatial positions for optimal interactions, relative interaction strength, as well as character. In addition, these maps are backbone angle-dependent. After map calculation and clustering, a finite number of unique residue sidechain interaction maps exist for each backbone conformation, with the number related to the residue's size and interaction complexity. Structures for soluble proteins (~749,000 residues) and membrane proteins (~387,000 residues) were analyzed, with the latter group being subdivided into three subsets related to the residue's position in the membrane protein: soluble domain, core-facing transmembrane domain, and lipid-facing transmembrane domain. This work suggests that maps representing residue types and their backbone conformation can be reassembled to optimize the medium-to-high resolution details of a protein structure. In particular, the information encoded in maps constructed from the lipid-facing transmembrane residues appears to paint a clear picture of the protein-lipid interactions that are difficult to obtain experimentally.

MDverse: Shedding Light on the Dark Matter of Molecular Dynamics Simulations

The rise of open science and the absence of a global dedicated data repository for molecular dynamics (MD) simulations has led to the accumulation of MD files in generalist data repositories, constituting the dark matter of MD - data that is technically accessible, but neither indexed, curated, or easily searchable. Leveraging an original search strategy, we found and indexed about 250,000 files and 2,000 datasets from Zenodo, Figshare and Open Science Framework. With a focus on files produced by the Gromacs MD software, we illustrate the potential offered by the mining of publicly available MD data. We identified systems with specific molecular composition and were able to characterize essential parameters of MD simulation such as temperature and simulation length, and could identify model resolution, such as all-atom and coarse-grain. Based on this analysis, we inferred metadata to propose a search engine prototype to explore the MD data. To continue in this direction, we call on the community to pursue the effort of sharing MD data, and to report and standardize metadata to reuse this valuable matter.

From complex data to clear insights: visualizing molecular dynamics trajectories

Advances in simulations, combined with technological developments in high-performance computing, have made it possible to produce a physically accurate dynamic representation of complex biological systems involving millions to billions of atoms over increasingly long simulation times. The analysis of these computed simulations is crucial, involving the interpretation of structural and dynamic data to gain insights into the underlying biological processes. However, this analysis becomes increasingly challenging due to the complexity of the generated systems with a large number of individual runs, ranging from hundreds to thousands of trajectories. This massive increase in raw simulation data creates additional processing and visualization challenges. Effective visualization techniques play a vital role in facilitating the analysis and interpretation of molecular dynamics simulations. In this paper, we focus mainly on the techniques and tools that can be used for visualization of molecular dynamics simulations, among which we highlight the few approaches used specifically for this purpose, discussing their advantages and limitations, and addressing the future challenges of molecular dynamics visualization.

KCNQ potassium channels modulate Wnt activity in gastro-oesophageal adenocarcinomas

Voltage-sensitive potassium channels play an important role in controlling membrane potential and ionic homeostasis in the gut and have been implicated in gastrointestinal (GI) cancers. Through large-scale analysis of 897 patients with gastro-oesophageal adenocarcinomas (GOAs) coupled with in vitro models, we find KCNQ family genes are mutated in ∼30% of patients, and play therapeutically targetable roles in GOA cancer growth. KCNQ1 and KCNQ3 mediate the WNT pathway and MYC to increase proliferation through resultant effects on cadherin junctions. This also highlights novel roles of KCNQ3 in non-excitable tissues. We also discover that activity of KCNQ3 sensitises cancer cells to existing potassium channel inhibitors and that inhibition of KCNQ activity reduces proliferation of GOA cancer cells. These findings reveal a novel and exploitable role of potassium channels in the advancement of human cancer, and highlight that supplemental treatments for GOAs may exist through KCNQ inhibitors.

Structural basis for the activity regulation of Salt Overly Sensitive 1 in Arabidopsis salt tolerance

The plasma membrane Na+/H+ exchanger Salt Overly Sensitive 1 (SOS1) is crucial for plant salt tolerance. Unlike typical sodium/proton exchangers, SOS1 contains a large cytoplasmic domain (CPD) that regulates Na+/H+ exchange activity. However, the underlying modulation mechanism remains unclear. Here we report the structures of SOS1 from Arabidopsis thaliana in two conformations, primarily differing in CPD flexibility. The CPD comprises an interfacial domain, a cyclic nucleotide-binding domain-like domain (CNBD-like domain) and an autoinhibition domain. Through yeast cell-based Na+ tolerance test, we reveal the regulatory role of the interfacial domain and the activation role of the CNBD-like domain. The CPD forms a negatively charged cavity that is connected to the ion binding site. The transport of Na+ may be coupled with the conformational change of CPD. These findings provide structural and functional insight into SOS1 activity regulation.

AQcalc: A web server that identifies weak molecular interactions in protein structures

Weak molecular interactions play an important role in protein structure and function. Computational tools that identify weak molecular interactions are, therefore, valuable for the study of proteins. Here, we present AQcalc, a web server (https://aqcalcbiocomputing.com/) that can be used to identify anion-quadrupole (AQ) interactions, which are weak interactions involving aromatic residue (Trp, Tyr, and Phe) ring edges and anions (Asp, Glu, and phosphate ion) both within proteins and at their interfaces (protein-protein, protein-nucleic acids, and protein-lipid bilayer). AQcalc identifies AQ interactions as well as clusters involving AQ, cation-π, and salt bridges, among others. Utilizing AQcalc we analyzed weak interactions in protein models, even in the absence of experimental structures, to understand the contributions of weak interactions to deleterious structural changes, including those associated with oncogenic and germline disease variants. We identified several deleterious variants with disrupted AQ interactions (comparable in frequency to cation-π disruptions). Amyloid fibrils utilize AQ to bury anions at frequencies that far exceed those observed for globular proteins. AQ interactions were detected three and five times more frequently than the hydrogen-bonded AQ (HBAQ) in fibril structures and protein-lipid bilayer interfaces, respectively. By contrast, AQ and HBAQ interactions were detected with similar frequencies in globular proteins. Collectively, these findings suggest AQcalc will be effective in facilitating fine structural analysis. As other web utilities designed to identify protein residue interaction networks do not report AQ interactions, wide use of AQcalc will enrich our understanding of residue interaction networks and facilitate hypothesis testing by identifying and experimentally characterizing these comparably weak but important interactions.

The energetics and ion coupling of cholesterol transport through Patched1

Patched1 (PTCH1) is a tumor suppressor protein of the mammalian Hedgehog (HH) signaling pathway, implicated in embryogenesis and tissue homeostasis. PTCH1 inhibits the G protein-coupled receptor Smoothened (SMO) via a debated mechanism involving modulating ciliary cholesterol accessibility. Using extensive molecular dynamics simulations and free energy calculations to evaluate cholesterol transport through PTCH1, we find an energetic barrier of ~15 to 20 kilojoule per mole for cholesterol export. In silico data are coupled to in vivo biochemical assays of PTCH1 mutants to probe coupling between cation binding sites, transmembrane motions, and PTCH1 activity. Using complementary simulations of Dispatched1, we find that transition between "inward-open" and solvent "occluded" states is accompanied by Na+-induced pinching of intracellular helical segments. Thus, our findings illuminate the energetics and ion coupling stoichiometries of PTCH1 transport mechanisms, whereby one to three Na+ or two to three K+ couple to cholesterol export, and provide the first molecular description of transitions between distinct transport states.

Expulsion mechanism of the substrate-translocating subunit in ECF transporters

Energy-coupling factor (ECF)-type transporters mediate the uptake of micronutrients in many bacteria. They consist of a substrate-translocating subunit (S-component) and an ATP-hydrolysing motor (ECF module) Previous data indicate that the S-component topples within the membrane to alternately expose the binding site to either side of the membrane. In many ECF transporters, the substrate-free S-component can be expelled from the ECF module. Here we study this enigmatic expulsion step by cryogenic electron microscopy and reveal that ATP induces a concave-to-convex shape change of two long helices in the motor, thereby destroying the S-component's docking site and allowing for its dissociation. We show that adaptation of the membrane morphology to the conformational state of the motor may favour expulsion of the substrate-free S-component when ATP is bound and docking of the substrate-loaded S-component after hydrolysis. Our work provides a picture of bilayer-assisted chemo-mechanical coupling in the transport cycle of ECF transporters.

Physicochemical Nature of SARS-CoV-2 Spike Protein Binding to Human Vimentin

Vimentin, a protein that builds part of the cytoskeleton and is involved in many aspects of cellular function, was recently identified as a cell surface attachment site for the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The present study investigated the physicochemical nature of the binding between the SARS-CoV-2 S1 glycoprotein receptor binding domain (S1 RBD) and human vimentin using atomic force microscopy and a quartz crystal microbalance. The molecular interactions of S1 RBD and vimentin proteins were quantified using vimentin monolayers attached to the cleaved mica or a gold microbalance sensor as well as in its native extracellular form present on the live cell surface. The presence of specific interactions between vimentin and S1 RBD was also confirmed using in silico studies. This work provides new evidence that cell-surface vimentin (CSV) functions as a site for SARS-CoV-2 virus attachment and is involved in the pathogenesis of Covid-19, providing a potential target for therapeutic countermeasures.

Class B1 GPCR activation by an intracellular agonist

G protein-coupled receptors (GPCRs) generally accommodate specific ligands in the orthosteric-binding pockets. Ligand binding triggers a receptor allosteric conformational change that leads to the activation of intracellular transducers, G proteins and β-arrestins. Because these signals often induce adverse effects, the selective activation mechanism for each transducer must be elucidated. Thus, many orthosteric-biased agonists have been developed, and intracellular-biased agonists have recently attracted broad interest. These agonists bind within the receptor intracellular cavity and preferentially tune the specific signalling pathway over other signalling pathways, without allosteric rearrangement of the receptor from the extracellular side1-3. However, only antagonist-bound structures are currently available1,4-6, and there is no evidence to support that biased agonist binding occurs within the intracellular cavity. This limits the comprehension of intracellular-biased agonism and potential drug development. Here we report the cryogenic electron microscopy structure of a complex of Gs and the human parathyroid hormone type 1 receptor (PTH1R) bound to a PTH1R agonist, PCO371. PCO371 binds within an intracellular pocket of PTH1R and directly interacts with Gs. The PCO371-binding mode rearranges the intracellular region towards the active conformation without extracellularly induced allosteric signal propagation. PCO371 stabilizes the significantly outward-bent conformation of transmembrane helix 6, which facilitates binding to G proteins rather than β-arrestins. Furthermore, PCO371 binds within the highly conserved intracellular pocket, activating 7 out of the 15 class B1 GPCRs. Our study identifies a new and conserved intracellular agonist-binding pocket and provides evidence of a biased signalling mechanism that targets the receptor-transducer interface.

3D interaction homology: The hydrophobic residues alanine, isoleucine, leucine, proline and valine play different structural roles in soluble and membrane proteins

The aliphatic hydrophobic amino acid residues—alanine, isoleucine, leucine, proline and valine—are among the most common found in proteins. Their structural role in proteins is seemingly obvious: engage in hydrophobic interactions to stabilize secondary, and to a lesser extent, tertiary and quaternary structure. However, favorable hydrophobic interactions involving the sidechains of these residue types are generally less significant than the unfavorable set arising from interactions with polar atoms. Importantly, the constellation of interactions between residue sidechains and their environments can be recorded as three-dimensional maps that, in turn, can be clustered. The clustered average map sets compose a library of interaction profiles encoding interaction strengths, interaction types and the optimal 3D position for the interacting partners. This library is backbone angle-dependent and suggests solvent and lipid accessibility for each unique interaction profile. In this work, in addition to analysis of soluble proteins, a large set of membrane proteins that contained optimized artificial lipids were evaluated by parsing the structures into three distinct components: soluble extramembrane domain, lipid facing transmembrane domain, core transmembrane domain. The aliphatic residues were extracted from each of these sets and passed through our calculation protocol. Notable observations include: the roles of aliphatic residues in soluble proteins and in the membrane protein’s soluble domains are nearly identical, although the latter are slightly more solvent accessible; by comparing maps calculated with sidechain-lipid interactions to maps ignoring those interactions, the potential extent of residue-lipid and residue-interactions can be assessed and likely exploited in structure prediction and modeling; amongst these residue types, the levels of lipid engagement show isoleucine as the most engaged, while the other residues are largely interacting with neighboring helical residues.

UNILIPID, a Methodology for Energetically Accurate Prediction of Protein Insertion into Implicit Membranes of Arbitrary Shape

The insertion of proteins into membranes is crucial for understanding their function in many biological processes. In this work, we present UNILIPID, a universal implicit lipid-protein description as a methodology for dealing with implicit membranes. UNILIPID is independent of the scale of representation and can be applied at the level of all atoms, coarse-grained particles down to the level of a single bead per amino acid. We provide example implementations for these scales and demonstrate the versatility of our approach by accurately reflecting the free energy of transfer for each amino acid. In addition to single membranes, we describe the analytical implementation of double membranes and show that UNILIPID is well suited for modeling at multiple scales. We generalize to membranes of arbitrary shape. With UNILIPID, we provide a methodological framework for a simple and general parameterization tuned to reproduce a selected reference hydrophobicity scale. The software we provide along with the methodological description is optimized for specific user features such as real-time response, live visual analysis, and virtual reality experiences.

The Energetics and Ion Coupling of Cholesterol Transport Through Patched1

Patched1 (PTCH1) is the principal tumour suppressor protein of the mammalian Hedgehog (HH) signalling pathway, implicated in embryogenesis and tissue homeostasis. PTCH1 inhibits the Class F G protein-coupled receptor Smoothened (SMO) via a debated mechanism involving modulating accessible cholesterol levels within ciliary membranes. Using extensive molecular dynamics (MD) simulations and free energy calculations to evaluate cholesterol transport through PTCH1, we find an energetic barrier of ~15-20 kJ mol -1 for cholesterol export. In simulations we identify cation binding sites within the PTCH1 transmembrane domain (TMD) which may provide the energetic impetus for cholesterol transport. In silico data are coupled to in vivo biochemical assays of PTCH1 mutants to probe coupling between transmembrane motions and PTCH1 activity. Using complementary simulations of Dispatched1 (DISP1) we find that transition between 'inward-open' and solvent 'occluded' states is accompanied by Na + induced pinching of intracellular helical segments. Thus, our findings illuminate the energetics and ion-coupling stoichiometries of PTCH1 transport mechanisms, whereby 1-3 Na + or 2-3 K + couple to cholesterol export, and provide the first molecular description of transitions between distinct transport states.

ProtRAP: Predicting Lipid Accessibility Together with Solvent Accessibility of Proteins in One Run

Solvent accessibility has been extensively used to characterize and predict the chemical properties of the surface residues of soluble proteins. However, there is not yet a widely accepted quantity of the same dimension for the study of lipid-accessible residues of membrane proteins. In this study, we propose that lipid accessibility, defined in a similar way to solvent accessibility, can be used to characterize the lipid-accessible residues of membrane proteins. Moreover, we developed a deep learning-based method, ProtRAP (Protein Relative Accessibility Predictor), to predict the relative lipid accessibility and relative solvent accessibility of residues from a given protein sequence, which can infer which residues are likely accessible to lipids, accessible to solvent, or buried in the protein interior in one run.

Machine learning in computational modelling of membrane protein sequences and structures: From methodologies to applications

Membrane proteins mediate a wide spectrum of biological processes, such as signal transduction and cell communication. Due to the arduous and costly nature inherent to the experimental process, membrane proteins have long been devoid of well-resolved atomic-level tertiary structures and, consequently, the understanding of their functional roles underlying a multitude of life activities has been hampered. Currently, computational tools dedicated to furthering the structure-function understanding are primarily focused on utilizing intelligent algorithms to address a variety of site-wise prediction problems (e.g., topology and interaction sites), but are scattered across different computing sources. Moreover, the recent advent of deep learning techniques has immensely expedited the development of computational tools for membrane protein-related prediction problems. Given the growing number of applications optimized particularly by manifold deep neural networks, we herein provide a review on the current status of computational strategies mainly in membrane protein type classification, topology identification, interaction site detection, and pathogenic effect prediction. Meanwhile, we provide an overview of how the entire prediction process proceeds, including database collection, data pre-processing, feature extraction, and method selection. This review is expected to be useful for developing more extendable computational tools specific to membrane proteins.

Molecular Insights into Substrate Binding of the Outer Membrane Enzyme OmpT

The enzyme OmpT of the outer membrane of Escherichia coli shows proteolytic activity and cleaves peptides and proteins. Using molecular dynamics simulations in a fully hydrated lipid bilayer on a time scale of hundreds of nanoseconds, we draw a detailed atomic picture of substrate recognition in the OmpT-holo enzyme complex. Hydrogen bonds and salt bridges are essential for maintaining the integrity of the active site and play a central role for OmpT in recognizing its substrate. Electrostatic interactions are critical at all stages from approaching the substrate to docking at the active site. Computational alanine scanning based on the Molecular Mechanics Generalized Born Surface Area (MM-GBSA) approach confirms the importance of multiple residues in the active site that form salt bridges. The substrate fluctuates along the axis of the β-barrel, which is associated with oscillations of the binding cleft formed by the residue pairs D210-H212 and D83-D85. Principal component analysis suggests that substrate and protein movements are correlated. We observe the transient presence of putative catalytic water molecules near the active site, which may be involved in the nucleophilic attack on the cleavable peptide bond of the substrate.

Strategies for the Purification of Membrane Proteins

Membrane proteins account for approximately 30% of the coding regions of all sequenced genomes, and they play crucial roles in many fundamental cell processes. However, there are relatively few membrane proteins with known three-dimensional structures. This is likely due to technical challenges associated with membrane protein extraction, solubilization, and purification. Membrane proteins are classified based on the level of interaction with membrane lipid bilayers, with peripheral membrane proteins associating non-covalently with the membrane, and integral membrane proteins associating more strongly by means of hydrophobic interactions. Generally speaking, peripheral membrane proteins can be purified by milder techniques than integral membrane proteins, with the latter's extraction requiring phospholipid bilayer disruption using detergents or organic solvents. In this chapter, important considerations for membrane protein purification are addressed, with a focus on the initial stages of membrane protein solubilization, where problems are most frequently encountered. Protocols are outlined for the extraction of peripheral membrane proteins, solubilization of integral membrane proteins, and sample clean-up and concentration.

The NALCN channel regulates metastasis and nonmalignant cell dissemination

We identify the sodium leak channel non-selective protein (NALCN) as a key regulator of cancer metastasis and nonmalignant cell dissemination. Among 10,022 human cancers, NALCN loss-of-function mutations were enriched in gastric and colorectal cancers. Deletion of Nalcn from gastric, intestinal or pancreatic adenocarcinomas in mice did not alter tumor incidence, but markedly increased the number of circulating tumor cells (CTCs) and metastases. Treatment of these mice with gadolinium-a NALCN channel blocker-similarly increased CTCs and metastases. Deletion of Nalcn from mice that lacked oncogenic mutations and never developed cancer caused shedding of epithelial cells into the blood at levels equivalent to those seen in tumor-bearing animals. These cells trafficked to distant organs to form normal structures including lung epithelium, and kidney glomeruli and tubules. Thus, NALCN regulates cell shedding from solid tissues independent of cancer, divorcing this process from tumorigenesis and unmasking a potential new target for antimetastatic therapies.

Predicted Structure and Functions of the Prototypic Alphaherpesvirus Herpes Simplex Virus Type-1 UL37 Tegument Protein

The alphaherpesvirus UL37 tegument protein is a highly conserved, multi-functional protein. Mutagenesis analysis delineated the UL37 domains necessary for retrograde transport and viral replication. Specifically, the amino-terminal 480 amino acids are dispensable for virus replication in epithelial cell culture, but it is unknown whether this amino-terminal deletion affects UL37 structure and intracellular transport in epithelial cells and neurons. To investigate the structure and function of UL37, we utilized multiple computational approaches to predict and characterize the secondary and tertiary structure and other functional features. The structure of HSV-1 UL37 and Δ481N were deduced using publicly available predictive algorithms. The predicted model of HSV-1 UL37 is a stable, multi-functional, globular monomer, rich in alpha helices, with unfolded regions within the linker and the C-tail domains. The highly flexible C-tail contains predicted binding sites to the dynein intermediate chain, as well as DNA and RNA. Predicted interactions with the cytoplasmic surface of the lipid membrane suggest UL37 is a peripheral membrane protein. The Δ481N truncation did not alter the predicted structure of the UL37 C-terminus protein and its predicted interaction with dynein. We validated these models by examining the replication kinetics and transport of the Δ481N virus toward the nuclei of infected epithelial and neuronal cells. The Δ481N virus had substantial defects in virus spread; however, it exhibited no apparent defects in virus entry and intracellular transport. Using computational analyses, we identified several key features of UL37, particularly the flexible unstructured tail; we then demonstrated that the UL37 C-terminus alone is sufficient to effectively transport the virus towards the nucleus of infected epithelial and neuronal cells.

The allosteric mechanism leading to an open-groove lipid conductive state of the TMEM16F scramblase

TMEM16F is a Ca²⁺-activated phospholipid scramblase in the TMEM16 family of membrane proteins. Unlike other TMEM16s exhibiting a membrane-exposed hydrophilic groove that serves as a translocation pathway for lipids, the experimentally determined structures of TMEM16F shows the groove in a closed conformation even under conditions of maximal scramblase activity. It is currently unknown if/how TMEM16F groove can open for lipid scrambling. Here we describe the analysis of ~400 µs all-atom molecular dynamics (MD) simulations of the TMEM16F revealing an allosteric mechanism leading to an open-groove, lipid scrambling competent state of the protein. The groove opens into a continuous hydrophilic conduit that is highly similar in structure to that seen in other activated scramblases. The allosteric pathway connects this opening to an observed destabilization of the Ca²⁺ ion bound at the distal site near the dimer interface, to the dynamics of specific protein regions that produces the open-groove state to scramble phospholipids.

Molecular dynamics simulations reveal the allosteric mechanism leading to an open, lipid scrambling competent state of a mammalian TMEM16F.

Endogenous ligand recognition and structural transition of a human PTH receptor

Endogenous parathyroid hormone (PTH) and PTH-related peptide (PTHrP) bind to the parathyroid hormone receptor 1 (PTH1R) and activate the stimulatory G-protein (Gs) signaling pathway. Intriguingly, the two ligands have distinct signaling and physiological properties: PTH evokes prolonged Gs activation, whereas PTHrP evokes transient Gs activation with reduced bone-resorption effects. The distinct molecular actions are ascribed to the differences in ligand recognition and dissociation kinetics. Here, we report cryoelectron microscopic structures of six forms of the human PTH1R-Gs complex in the presence of PTH or PTHrP at resolutions of 2.8 -4.1 Å. A comparison of the PTH-bound and PTHrP-bound structures reveals distinct ligand-receptor interactions underlying the ligand affinity and selectivity. Furthermore, five distinct PTH-bound structures, combined with computational analyses, provide insights into the unique and complex process of ligand dissociation from the receptor and shed light on the distinct durations of signaling induced by PTH and PTHrP.

Self-Assembly of Lipid Mixtures in Solutions: Structures, Dynamics Processes and Mechanical Properties

The self-assembly of lipid mixtures in aqueous solution was investigated by dissipative particle dynamics simulation. Two types of lipid molecules were modelled, where three mixed structures, i.e., the membrane, perforated membrane and vesicle, were determined in the self-assembly processes. Phase behaviour was investigated by using the phase diagrams based on the tail chain lengths for the two types of lipids. Several parameters, such as chain number and average radius of gyration, were employed to explore the structural formations of the membrane and perforated membrane in the dynamic processes. Interface tension was used to demonstrate the mechanical properties of the membrane and perforated membrane in the equilibrium state and dynamics processes. Results help us to understand the self-assembly mechanism of the biomolecule mixtures, which has a potential application for designing the lipid molecule-based bio-membranes in solutions.

TMEM16 scramblases thin the membrane to enable lipid scrambling

TMEM16 scramblases dissipate the plasma membrane lipid asymmetry to activate multiple eukaryotic cellular pathways. Scrambling was proposed to occur with lipid headgroups moving between leaflets through a membrane-spanning hydrophilic groove. Direct information on lipid-groove interactions is lacking. We report the 2.3 Å resolution cryogenic electron microscopy structure of the nanodisc-reconstituted Ca2+-bound afTMEM16 scramblase showing how rearrangement of individual lipids at the open pathway results in pronounced membrane thinning. Only the groove's intracellular vestibule contacts lipids, and mutagenesis suggests scrambling does not require specific protein-lipid interactions with the extracellular vestibule. We find scrambling can occur outside a closed groove in thinner membranes and is inhibited in thicker membranes, despite an open pathway. Our results show afTMEM16 thins the membrane to enable scrambling and that an open hydrophilic pathway is not a structural requirement to allow rapid transbilayer movement of lipids. This mechanism could be extended to other scramblases lacking a hydrophilic groove.

Membrane contact probability: An essential and predictive character for the structural and functional studies of membrane proteins

One of the unique traits of membrane proteins is that a significant fraction of their hydrophobic amino acids is exposed to the hydrophobic core of lipid bilayers rather than being embedded in the protein interior, which is often not explicitly considered in the protein structure and function predictions. Here, we propose a characteristic and predictive quantity, the membrane contact probability (MCP), to describe the likelihood of the amino acids of a given sequence being in direct contact with the acyl chains of lipid molecules. We show that MCP is complementary to solvent accessibility in characterizing the outer surface of membrane proteins, and it can be predicted for any given sequence with a machine learning-based method by utilizing a training dataset extracted from MemProtMD, a database generated from molecular dynamics simulations for the membrane proteins with a known structure. As the first of many potential applications, we demonstrate that MCP can be used to systematically improve the prediction precision of the protein contact maps and structures.

The distribution of residues on protein surfaces is largely determined by the surrounding environment. For soluble proteins, most of the residues on the outer surface are hydrophilic, and people use the quantity “solvent accessibility” to describe and predict these surface residues. In contrast, for membrane proteins that are embedded in a lipid bilayer, many of their surface residues are hydrophobic and membrane-contacting, but there is yet a widely-accepted quantity for the description or prediction of this characteristic property. Here, we propose a new quantity termed “membrane contact probability (MCP)”, which can be used to describe and predict the membrane-contacting surface residues of proteins. We also propose a machine learning-based method to predict MCP from protein sequences, utilizing the dataset generated by physics-based computer simulations. We demonstrate that a quantity such as MCP is helpful for protein structure prediction, and we believe that it will find broad applications in the structure and function studies of membrane proteins.

Assembly and Analysis of Cell-Scale Membrane Envelopes

The march toward exascale computing will enable routine molecular simulation of larger and more complex systems, for example, simulation of entire viral particles, on the scale of approximately billions of atoms─a simulation size commensurate with a small bacterial cell. Anticipating the future hardware capabilities that will enable this type of research and paralleling advances in experimental structural biology, efforts are currently underway to develop software tools, procedures, and workflows for constructing cell-scale structures. Herein, we describe our efforts in developing and implementing an efficient and robust workflow for construction of cell-scale membrane envelopes and embedding membrane proteins into them. A new approach for construction of massive membrane structures that are stable during the simulations is built on implementing a subtractive assembly technique coupled with the development of a structure concatenation tool (fastmerge), which eliminates overlapping elements based on volumetric criteria rather than adding successive molecules to the simulation system. Using this approach, we have constructed two "protocells" consisting of MARTINI coarse-grained beads to represent cellular membranes, one the size of a cellular organelle and another the size of a small bacterial cell. The membrane envelopes constructed here remain whole during the molecular dynamics simulations performed and exhibit water flux only through specific proteins, demonstrating the success of our methodology in creating tight cell-like membrane compartments. Extended simulations of these cell-scale structures highlight the propensity for nonspecific interactions between adjacent membrane proteins leading to the formation of protein microclusters on the cell surface, an insight uniquely enabled by the scale of the simulations. We anticipate that the experiences and best practices presented here will form the basis for the next generation of cell-scale models, which will begin to address the addition of soluble proteins, nucleic acids, and small molecules essential to the function of a cell.

Lateral access mechanism of LPA receptor probed by molecular dynamics simulation

G-protein-coupled receptors (GPCR) are a family of membrane receptors that play important roles in the regulation of various physiological phenomena. LPA receptors (LPA_1-6) are members of the class A GPCRs, which transduce a lysophosphatidic acid (LPA) signal across the cell membrane and evoke various responses, including cellular survival, proliferation, differentiation, and migration. The crystal structure of LPA₆ revealed a gap between its transmembrane helices (TMs), which is opened toward the membrane side. This led to the proposal of the “lateral access model,” in which its lipophilic ligand directly enters the binding pocket through the gap structure at the membrane. In this study, we performed molecular dynamics (MD) simulations and Markov state model (MSM) analyses of LPA₆ and LPA, to elucidate the long timescale dynamics of the ligand binding process. The results from the 71.4-μs MD simulation suggested that the flexibility of the TMs constituting the gap structure enables the lateral entrance of the ligand, and the key interactions between the receptor and ligand facilitate the transition state of the ligand binding process.

Polymodal Control of TMEM16x Channels and Scramblases

The TMEM16A/anoctamin-1 calcium-activated chloride channel (CaCC) contributes to a range of vital functions, such as the control of vascular tone and epithelial ion transport. The channel is a founding member of a family of 10 proteins (TMEM16x) with varied functions; some members (i.e., TMEM16A and TMEM16B) serve as CaCCs, while others are lipid scramblases, combine channel and scramblase function, or perform additional cellular roles. TMEM16x proteins are typically activated by agonist-induced Ca2+ release evoked by Gq-protein-coupled receptor (GqPCR) activation; thus, TMEM16x proteins link Ca2+-signalling with cell electrical activity and/or lipid transport. Recent studies demonstrate that a range of other cellular factors—including plasmalemmal lipids, pH, hypoxia, ATP and auxiliary proteins—also control the activity of the TMEM16A channel and its paralogues, suggesting that the TMEM16x proteins are effectively polymodal sensors of cellular homeostasis. Here, we review the molecular pathophysiology, structural biology, and mechanisms of regulation of TMEM16x proteins by multiple cellular factors.

Prediction of molecular interactions and physicochemical properties relevant for vasopressin V2 receptor antagonism

We have developed two ligand- and receptor-based computational approaches to study the physicochemical properties relevant to the biological activity of vasopressin V2 receptor (V2R) antagonist and eventually to predict the expected binding mode to V2R. The obtained quantitative structure activity relationship (QSAR) model showed a correlation of the antagonist activity with the hydration energy (E_H2O), the polarizability (P), and the calculated partial charge on atom N7 (q6) of the common substructure. The first two descriptors showed a positive contribution to antagonist activity, while the third one had a negative contribution. V2R was modeled and further relaxed on a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocoline (POPC) membrane by molecular dynamics simulations. The receptor antagonist complexes were guessed by molecular docking, and the stability of the most relevant structures was also evaluated by molecular dynamics simulations. As a result, amino acid residues Q96, W99, F105, K116, F178, A194, F307, and M311 were identified with the probably most relevant antagonist-receptor interactions on the studied complexes. The proposed QSAR model could explain the molecular properties relevant to the antagonist activity. The contributions to the antagonist-receptor interaction appeared also in agreement with the binding mode of the complexes obtained by molecular docking and molecular dynamics. These models will be used in further studies to look for new V2R potential antagonist molecules.

Spatial arrangement of proteins in planar and curved membranes by PPM 3.0

Cellular protrusions, invaginations, and many intracellular organelles have strongly curved membrane regions. Transmembrane and peripheral membrane proteins that induce, sense, or stabilize such regions cannot be properly fitted into a single flat bilayer. To treat such proteins, we developed a new method and a web tool, PPM 3.0, for positioning proteins in curved or planar, single or multiple membranes. This method determines the energetically optimal spatial position, the hydrophobic thickness, and the radius of intrinsic curvature of a membrane-deforming protein structure by arranging it in a single or several sphere-shaped or planar membrane sections. In addition, it can define the lipid-embedded regions of a protein that simultaneously spans several membranes or determine the optimal position of a peptide in a spherical micelle. The PPM 3.0 web server operates with 17 types of biological membranes and 4 types of artificial bilayers. It is publicly available at https://opm.phar.umich.edu/ppm_server3. PPM 3.0 was applied to identify and characterize arrangements in membranes of 128 proteins with a significant intrinsic curvature, such as BAR domains, annexins, Piezo, and MscS mechanosensitive channels, cation-chloride cotransporters, as well as mitochondrial ATP synthases, calcium uniporters, and TOM complexes. These proteins form large complexes that are mainly localized in mitochondria, plasma membranes, and endosomes. Structures of bacterial drug efflux pumps, AcrAB-TolC, MexAB-OrpM, and MacAB-TolC, were positioned in both membranes of the bacterial cell envelop, while structures of multimeric gap-junction channels were arranged in two opposed cellular membranes.

Water Nanoconfined in a Hydrophobic Pore: Molecular Dynamics Simulations of Transmembrane Protein 175 and the Influence of Water Models

Water molecules within biological ion channels are in a nanoconfined environment and therefore exhibit behaviors which differ from that of bulk water. Here, we investigate the phenomenon of hydrophobic gating, the process by which a nanopore may spontaneously dewet to form a "vapor lock" if the pore is sufficiently hydrophobic and/or narrow. This occurs without steric occlusion of the pore. Using molecular dynamics simulations with both rigid fixed-charge and polarizable (AMOEBA) force fields, we investigate this wetting/dewetting behavior in the transmembrane protein 175 ion channel. We examine how a range of rigid fixed-charge and polarizable water models affect wetting/dewetting in both the wild-type structure and in mutants chosen to cover a range of nanopore radii and pore-lining hydrophobicities. Crucially, we find that the rigid fixed-charge water models lead to similar wetting/dewetting behaviors, but that the polarizable water model resulted in an increased wettability of the hydrophobic gating region of the pore. This has significant implications for molecular simulations of nanoconfined water, as it implies that polarizability may need to be included if we are to gain detailed mechanistic insights into wetting/dewetting processes. These findings are of importance for the design of functionalized biomimetic nanopores (e.g., sensing or desalination) as well as for furthering our understanding of the mechanistic processes underlying biological ion channel function.

Cryo-EM reveals mechanistic insights into lipid-facilitated polyamine export by human ATP13A2

The cytoplasmic polyamine maintains cellular homeostasis by chelating toxic metal cations, regulating transcriptional activity, and protecting DNA. ATP13A2 was identified as a lysosomal polyamine exporter responsible for polyamine release into the cytosol, and its dysfunction is associated with Alzheimer's disease and other neural degradation diseases. ATP13A2 belongs to the P5 subfamily of the P-type ATPase family, but its mechanisms remain unknown. Here, we report the cryoelectron microscopy (cryo-EM) structures of human ATP13A2 under four different conditions, revealing the structural coupling between the polyamine binding and the dephosphorylation. Polyamine is bound at the luminal tunnel and recognized through numerous electrostatic and π-cation interactions, explaining its broad specificity. The unique N-terminal domain is anchored to the lipid membrane to stabilize the E2P conformation, thereby accelerating the E1P-to-E2P transition. These findings reveal the distinct mechanism of P5B ATPases, thereby paving the way for neuroprotective therapy by activating ATP13A2.

Relative Affinities of Protein-Cholesterol Interactions from Equilibrium Molecular Dynamics Simulations

Specific interactions of lipids with membrane proteins contribute to protein stability and function. Multiple lipid interactions surrounding a membrane protein are often identified in molecular dynamics (MD) simulations and are, increasingly, resolved in cryo-electron microscopy (cryo-EM) densities. Determining the relative importance of specific interaction sites is aided by determination of lipid binding affinities using experimental or simulation methods. Here, we develop a method for determining protein-lipid binding affinities from equilibrium coarse-grained MD simulations using binding saturation curves, designed to mimic experimental protocols. We apply this method to directly obtain affinities for cholesterol binding to multiple sites on a range of membrane proteins and compare our results with free energies obtained from density-based equilibrium methods and with potential of mean force calculations, getting good agreement with respect to the ranking of affinities for different sites. Thus, our binding saturation method provides a robust, high-throughput alternative for determining the relative consequence of individual sites seen in, e.g., cryo-EM derived membrane protein structures surrounded by an array of ancillary lipid densities.

3D interaction homology: Hydropathic interaction environments of serine and cysteine are strikingly different and their roles adapt in membrane proteins

Atomic-resolution protein structural models are prerequisites for many downstream activities like structure-function studies or structure-based drug discovery. Unfortunately, this data is often unavailable for some of the most interesting and therapeutically important proteins. Thus, computational tools for building native-like structural models from less-than-ideal experimental data are needed. To this end, interaction homology exploits the character, strength and loci of the sets of interactions that define a structure. Each residue type has its own limited set of backbone angle-dependent interaction motifs, as defined by their environments. In this work, we characterize the interactions of serine, cysteine and S-bridged cysteine in terms of 3D hydropathic environment maps. As a result, we explore several intriguing questions. Are the environments different between the isosteric serine and cysteine residues? Do some environments promote the formation of cystine S-S bonds? With the increasing availability of structural data for water-insoluble membrane proteins, are there environmental differences for these residues between soluble and membrane proteins? The environments surrounding serine and cysteine residues are dramatically different: serine residues are about 50% solvent exposed, while cysteines are only 10% exposed; the latter are more involved in hydrophobic interactions although there are backbone angle-dependent differences. Our analysis suggests that one driving force for -S-S- bond formation is a rather substantial increase in burial and hydrophobic interactions in cystines. Serine and cysteine become less and more, respectively, solvent-exposed in membrane proteins. 3D hydropathic environment maps are an evolving structure analysis tool showing promise as elements in a new protein structure prediction paradigm.

The Role of the Membrane in Transporter Folding and Activity

The synthesis, folding, and function of membrane transport proteins are critical factors for defining cellular physiology. Since the stability of these proteins evolved amidst the lipid bilayer, it is no surprise that we are finding that many of these membrane proteins demonstrate coupling of their structure or activity in some way to the membrane. More and more transporter structures are being determined with some information about the surrounding membrane, and computational modeling is providing further molecular details about these solvation structures. Thus, the field is moving towards identifying which molecular mechanisms - lipid interactions, membrane perturbations, differential solvation, and bulk membrane effects - are involved in linking membrane energetics to transporter stability and function. In this review, we present an overview of these mechanisms and the growing evidence that the lipid bilayer is a major determinant of the fold, form, and function of membrane transport proteins in membranes.

Identification and assessment of cardiolipin interactions with E. coli inner membrane proteins

Integral membrane proteins are localized and/or regulated by lipids present in the surrounding bilayer. While bacteria have relatively simple membranes, there is ample evidence that many bacterial proteins bind to specific lipids, especially the anionic lipid cardiolipin. Here, we apply molecular dynamics simulations to assess lipid binding to 42 different Escherichia coli inner membrane proteins. Our data reveal an asymmetry between the membrane leaflets, with increased anionic lipid binding to the inner leaflet regions of the proteins, particularly for cardiolipin. From our simulations, we identify >700 independent cardiolipin binding sites, allowing us to identify the molecular basis of a prototypical cardiolipin binding site, which we validate against structures of bacterial proteins bound to cardiolipin. This allows us to construct a set of metrics for defining a high-affinity cardiolipin binding site on bacterial membrane proteins, paving the way for a heuristic approach to defining other protein-lipid interactions.

Cryo-EM structure of PepT2 reveals structural basis for proton-coupled peptide and prodrug transport in mammals

The SLC15 family of proton-coupled solute carriers PepT1 and PepT2 play a central role in human physiology as the principal route for acquiring and retaining dietary nitrogen. A remarkable feature of the SLC15 family is their extreme substrate promiscuity, which has enabled the targeting of these transporters for the improvement of oral bioavailability for several prodrug molecules. Although recent structural and biochemical studies on bacterial homologs have identified conserved sites of proton and peptide binding, the mechanism of peptide capture and ligand promiscuity remains unclear for mammalian family members. Here, we present the cryo-electron microscopy structure of the outward open conformation of the rat peptide transporter PepT2 in complex with an inhibitory nanobody. Our structure, combined with molecular dynamics simulations and biochemical and cell-based assays, establishes a framework for understanding peptide and prodrug recognition within this pharmaceutically important transporter family.

Cryo-EM structure of the β3-adrenergic receptor reveals the molecular basis of subtype selectivity

The β₃-adrenergic receptor (β₃AR) is predominantly expressed in adipose tissue and urinary bladder and has emerged as an attractive drug target for the treatment of type 2 diabetes, obesity, and overactive bladder (OAB). Here, we report the cryogenic electron microscopy structure of the β₃AR-G_s signaling complex with the selective agonist mirabegron, a first-in-class drug for OAB. Comparison of this structure with the previously reported β₁AR and β₂AR structures reveals a receptor activation mechanism upon mirabegron binding to the orthosteric site. Notably, the narrower exosite in β₃AR creates a perpendicular pocket for mirabegron. Mutational analyses suggest that a combination of both the exosite shape and the amino-acid-residue substitutions defines the drug selectivity of the βAR agonists. Our findings provide a molecular basis for βAR subtype selectivity, allowing the design of more-selective agents with fewer adverse effects.

Cryo-EM structures of the caspase-activated protein XKR9 involved in apoptotic lipid scrambling

The exposure of the negatively charged lipid phosphatidylserine on the cell surface, catalyzed by lipid scramblases, is an important signal for the clearance of apoptotic cells by macrophages. The protein XKR9 is a member of a conserved family that has been associated with apoptotic lipid scrambling. Here, we describe structures of full-length and caspase-treated XKR9 from Rattus norvegicus in complex with a synthetic nanobody determined by cryo-electron microscopy. The 43 kDa monomeric membrane protein can be divided into two structurally related repeats, each containing four membrane-spanning segments and a helix that is partly inserted into the lipid bilayer. In the full-length protein, the C-terminus interacts with a hydrophobic pocket located at the intracellular side acting as an inhibitor of protein function. Cleavage by caspase-3 at a specific site releases 16 residues of the C-terminus, thus making the pocket accessible to the cytoplasm. Collectively, the work has revealed the unknown architecture of the XKR family and has provided initial insight into its activation by caspases.

Maintenance of organellar protein homeostasis by ER-associated degradation and related mechanisms

Protein homeostasis mechanisms are fundamentally important to match cellular needs and to counteract stress conditions. A fundamental challenge is to understand how defective proteins are recognized and extracted from cellular organelles to be degraded in the cytoplasm. The endoplasmic reticulum (ER)-associated degradation (ERAD) pathway is the best-understood organellar protein quality control system. Here, we review new insights into the mechanism of recognition and retrotranslocation of client proteins in ERAD. In addition to the membrane-integral ERAD E3 ubiquitin ligases, we highlight one protein family that is remarkably often involved in various aspects of membrane protein quality control and protein dislocation: the rhomboid superfamily, which includes derlins and intramembrane serine proteases. Rhomboid-like proteins have been found to control protein homeostasis in the ER, but also in other eukaryotic organelles and in bacteria, pointing toward conserved principles of membrane protein quality control across organelles and evolution.

Structural dynamics bridge the gap between the genetic and functional levels of GPCRs

G protein-coupled receptors (GPCRs) are implicated in nearly all physiological processes in the human body and represent an important drug targeting class. The genes encoding the different GPCR (sub)types determine their specific functionality, which can be altered by natural genetic variants and isoforms. Deciphering the molecular link between sequence diversity and its functional consequences is a current challenge and critical for the comprehension of the physiological response of GPCRs. It requires a global understanding of how protein sequence translates into protein structure, how this impacts the structural motions of the protein, and, finally, how all these factors determine the receptor functionality. Here, we discuss available resources and state-of-the-art computational approaches to address this question.

Factors Governing Selectivity of Dopamine Receptor Binding Compounds for D2R and D3R Subtypes

Targeting the D3 dopamine receptor (D3R) is a promising pharmacotherapeutic strategy for the treatment of many disorders. The structure of the D3R is similar to the D2 dopamine receptor (D2R), especially in the transmembrane spanning regions that form the orthosteric binding site, making it difficult to identify D3R selective pharmacotherapeutic agents. Here, we examine the molecular basis for the high affinity D3R binding and D3R vs D2R binding selectivity of substituted phenylpiperazine thiopheneamides. We show that removing the thiophenearylamide portion of the ligand consistently decreases the affinity of these ligands at D3R, while not affecting their affinity at the D2R. Our long (>10 μs) molecular dynamics simulations demonstrated that both dopamine receptor subtypes adopt two major conformations that we refer to as closed or open conformations, with D3R sampling the open conformation more frequently than D2R. The binding of ligands with conjoined orthosteric-allosteric binding moieties causes the closed conformation to populate more often in the trajectories. Also, significant differences were observed in the extracellular loops (ECL) of these two receptor subtypes leading to the identification of several residues that contribute differently to the ligand binding for the two receptors that could potentially contribute to ligand binding selectivity. Our observations also suggest that the displacement of ordered water in the binding pocket of D3R contributes to the affinity of the compounds containing an allosteric binding motif. These studies provide a better understanding of how a bitopic mode of engagement can determine ligands that bind selectively to D2 and D3 dopamine receptor subtypes.

Structural basis for allosteric control of the SERCA-Phospholamban membrane complex by Ca2+ and phosphorylation

Phospholamban (PLN) is a mini-membrane protein that directly controls the cardiac Ca²⁺-transport response to β-adrenergic stimulation, thus modulating cardiac output during the fight-or-flight response. In the sarcoplasmic reticulum membrane, PLN binds to the sarco(endo)plasmic reticulum Ca²⁺-ATPase (SERCA), keeping this enzyme's function within a narrow physiological window. PLN phosphorylation by cAMP-dependent protein kinase A or increase in Ca²⁺ concentration reverses the inhibitory effects through an unknown mechanism. Using oriented-sample solid-state NMR spectroscopy and replica-averaged NMR-restrained structural refinement, we reveal that phosphorylation of PLN's cytoplasmic regulatory domain signals the disruption of several inhibitory contacts at the transmembrane binding interface of the SERCA-PLN complex that are propagated to the enzyme's active site, augmenting Ca²⁺ transport. Our findings address long-standing questions about SERCA regulation, epitomizing a signal transduction mechanism operated by posttranslationally modified bitopic membrane proteins.

The Groovy TMEM16 Family: Molecular Mechanisms of Lipid Scrambling and Ion Conduction

The TMEM16 family of membrane proteins displays a remarkable functional dichotomy - while some family members function as Ca²⁺-activated anion channels, the majority of characterized TMEM16 homologs are Ca²⁺-activated lipid scramblases, which catalyze the exchange of phospholipids between the two membrane leaflets. Furthermore, some TMEM16 scramblases can also function as channels. Due to their involvement in important physiological processes, the family has been actively studied ever since their molecular identity was unraveled. In this review, we will summarize the recent advances in the field and how they influenced our view of TMEM16 family function and evolution. Structural, functional and computational studies reveal how relatively small rearrangements in the permeation pathway are responsible for the observed functional duality: while TMEM16 scramblases can adopt both ion- and lipid conductive conformations, TMEM16 channels can only populate the former. Recent data further provides the molecular details of a stepwise activation mechanism, which is initiated by Ca²⁺ binding and modulated by various cellular factors, including lipids. TMEM16 function and the surrounding membrane properties are inextricably intertwined, with the protein inducing bilayer deformations associated with scrambling, while the surrounding lipids modulate TMEM16 conformation and activity.

Visualizing protein structures - tools and trends

Molecular visualization is fundamental in the current scientific literature, textbooks and dissemination materials. It provides an essential support for presenting results, reasoning on and formulating hypotheses related to molecular structure. Tools for visual exploration of structural data have become easily accessible on a broad variety of platforms thanks to advanced software tools that render a great service to the scientific community. These tools are often developed across disciplines bridging computer science, biology and chemistry. This mini-review was written as a short and compact overview for scientists who need to visualize protein structures and want to make an informed decision which tool they should use. Here, we first describe a few 'Swiss Army knives' geared towards protein visualization for everyday use with an existing large user base, then focus on more specialized tools for peculiar needs that are not yet as broadly known. Our selection is by no means exhaustive, but reflects a diverse snapshot of scenarios that we consider informative for the reader. We end with an account of future trends and perspectives.

Guardians of the Cell: State-of-the-Art of Membrane Proteins from a Computational Point-of-View

Membrane proteins (MPs) encompass a large family of proteins with distinct cellular functions, and although representing over 50% of existing pharmaceutical drug targets, their structural and functional information is still very scarce. Over the last years, in silico analysis and algorithm development were essential to characterize MPs and overcome some limitations of experimental approaches. The optimization and improvement of these methods remain an ongoing process, with key advances in MPs' structure, folding, and interface prediction being continuously tackled. Herein, we discuss the latest trends in computational methods toward a deeper understanding of the atomistic and mechanistic details of MPs.

Structure-Function of TMEM16 Ion Channels and Lipid Scramblases

The TMEM16 protein family comprises two novel classes of structurally conserved but functionally distinct membrane transporters that function as Ca2+-dependent Cl- channels (CaCCs) or dual functional Ca2+-dependent ion channels and phospholipid scramblases. Extensive functional and structural studies have advanced our understanding of TMEM16 molecular mechanisms and physiological functions. TMEM16A and TMEM16B CaCCs control transepithelial fluid transport, smooth muscle contraction, and neuronal excitability, whereas TMEM16 phospholipid scramblases mediate the flip-flop of phospholipids across the membrane to allow phosphatidylserine externalization, which is essential in a plethora of important processes such as blood coagulation, bone development, and viral and cell fusion. In this chapter, we summarize the major methods in studying TMEM16 ion channels and scramblases and then focus on the current mechanistic understanding of TMEM16 Ca2+- and voltage-dependent channel gating as well as their ion and phospholipid permeation.

The MscS-like channel YnaI has a gating mechanism based on flexible pore helices

The mechanosensitive channel of small conductance (MscS) is the prototype of an evolutionarily diversified large family that fine-tunes osmoregulation but is likely to fulfill additional functions. Escherichia coli has six osmoprotective paralogs with different numbers of transmembrane helices. These helices are important for gating and sensing in MscS but the role of the additional helices in the paralogs is not understood. The medium-sized channel YnaI was extracted and delivered in native nanodiscs in closed-like and open-like conformations using the copolymer diisobutylene/maleic acid (DIBMA) for structural studies. Here we show by electron cryomicroscopy that YnaI has an extended sensor paddle that during gating relocates relative to the pore concomitant with bending of a GGxGG motif in the pore helices. YnaI is the only one of the six paralogs that has this GGxGG motif allowing the sensor paddle to move outward. Access to the pore is through a vestibule on the cytosolic side that is fenestrated by side portals. In YnaI, these portals are obstructed by aromatic side chains but are still fully hydrated and thus support conductance. For comparison with large-sized channels, we determined the structure of YbiO, which showed larger portals and a wider pore with no GGxGG motif. Further in silico comparison of MscS, YnaI, and YbiO highlighted differences in the hydrophobicity and wettability of their pores and vestibule interiors. Thus, MscS-like channels of different sizes have a common core architecture but show different gating mechanisms and fine-tuned conductive properties.

Evaluating inositol phospholipid interactions with inward rectifier potassium channels and characterising their role in disease

Membrane proteins are frequently modulated by specific protein-lipid interactions. The activation of human inward rectifying potassium (hKir) channels by phosphoinositides (PI) has been well characterised. Here, we apply a coarse-grained molecular dynamics free-energy perturbation (CG-FEP) protocol to capture the energetics of binding of PI lipids to hKir channels. By using either a single- or multi-step approach, we establish a consistent value for the binding of PIP2 to hKir channels, relative to the binding of the bulk phosphatidylcholine phospholipid. Furthermore, by perturbing amino acid side chains on hKir6.2, we show that the neonatal diabetes mutation E179K increases PIP2 affinity, while the congenital hyperinsulinism mutation K67N results in a reduced affinity. We show good agreement with electrophysiological data where E179K exhibits a reduction in neomycin sensitivity, implying that PIP2 binds more tightly E179K channels. This illustrates the application of CG-FEP to compare affinities between lipid species, and for annotating amino acid residues.