Large Membrane Protein Research Articles

MiLoPYP: self-supervised molecular pattern mining and particle localization in situ

Cryo-electron tomography allows the routine visualization of cellular landscapes in three dimensions at nanometer-range resolutions. When combined with single-particle tomography, it is possible to obtain near-atomic resolution structures of frequently occurring macromolecules within their native environment. Two outstanding challenges associated with cryo-electron tomography/single-particle tomography are the automatic identification and localization of proteins, tasks that are hindered by the molecular crowding inside cells, imaging distortions characteristic of cryo-electron tomography tomograms and the sheer size of tomographic datasets. Current methods suffer from low accuracy, demand extensive and time-consuming manual labeling or are limited to the detection of specific types of proteins. Here, we present MiLoPYP, a two-step dataset-specific contrastive learning-based framework that enables fast molecular pattern mining followed by accurate protein localization. MiLoPYP’s ability to effectively detect and localize a wide range of targets including globular and tubular complexes as well as large membrane proteins, will contribute to streamline and broaden the applicability of high-resolution workflows for in situ structure determination.

Druglike Molecules Binding to Large Membrane Proteins: Absolute Binding Free Energy Computation.

In this research, we employed the alchemical double-decoupling method alongside restraining potentials, coupled with the FEPMD method, to ascertain the standard binding free energy of a drug-like molecule termed BHQ and three analogous compounds engineered with progressive addition of bulky para-alkyl groups binding to SERCA (Ca2+-ATPase of skeletal muscle sarcoplasmic reticulum). Integral transmembrane proteins represent crucial drug targets in numerous therapeutic interventions, presenting computational challenges due to their considerable system sizes. Our approach integrated the generalized born potential method and the spherical solvent boundary potential method, allowing us to explicitly focus on the active binding site while treating the remainder of the system implicitly. We evaluated contributions to the standard binding free energy from distinct interaction potentials: electrostatic, repulsive, dispersive, and restraining potentials, computed separately. The resulting absolute binding free energy of BHQ (11.63 kcal/mol) closely aligns with experimental measurements (10.56 kcal/mol). Notably, an accurate estimation of the absolute binding free energy was achieved for the simplest analog, created with the addition of a single para-methyl group. However, the analog with two para-methyl groups exhibited the highest binding free energy, which disagreed with experimental results. Determining the binding free energy of the BHQ analog engineered with three para-methyl groups presented challenges in convergence and resulted in the lowest free energy among the three.

In situ serial crystallography facilitates 96-well plate structural analysis at low symmetry.

The advent of serial crystallography has rejuvenated and popularized room-temperature X-ray crystal structure determination. Structures determined at physiological temperature reveal protein flexibility and dynamics. In addition, challenging samples (e.g. large complexes, membrane proteins and viruses) form fragile crystals that are often difficult to harvest for cryo-crystallography. Moreover, a typical serial crystallography experiment requires a large number of microcrystals, mainly achievable through batch crystallization. Many medically relevant samples are expressed in mammalian cell lines, producing a meager quantity of protein that is incompatible with batch crystallization. This can limit the scope of serial crystallography approaches. Direct in situ data collection from a 96-well crystallization plate enables not only the identification of the best diffracting crystallization condition but also the possibility for structure determination under ambient conditions. Here, we describe an in situ serial crystallography (iSX) approach, facilitating direct measurement from crystallization plates mounted on a rapidly exchangeable universal plate holder deployed at a microfocus beamline, ID23-2, at the European Synchrotron Radiation Facility. We applied our iSX approach on a challenging project, autotaxin, a therapeutic target expressed in a stable human cell line, to determine the structure in the lowest-symmetry P1 space group at 3.0 Å resolution. Our in situ data collection strategy provided a complete dataset for structure determination while screening various crystallization conditions. Our data analysis reveals that the iSX approach is highly efficient at a microfocus beamline, improving throughput and demonstrating how crystallization plates can be routinely used as an alternative method of presenting samples for serial crystallography experiments at synchrotrons.

Pleiotropic physiological functions of Piezo1 in human body and its effect on malignant behavior of tumors.

Mechanosensitive ion channel protein 1 (Piezo1) is a large homotrimeric membrane protein. Piezo1 has various effects and plays an important and irreplaceable role in the maintenance of human life activities and homeostasis of the internal environment. In addition, recent studies have shown that Piezo1 plays a vital role in tumorigenesis, progression, malignancy and clinical prognosis. Piezo1 is involved in regulating the malignant behaviors of a variety of tumors, including cellular metabolic reprogramming, unlimited proliferation, inhibition of apoptosis, maintenance of stemness, angiogenesis, invasion and metastasis. Moreover, Piezo1 regulates tumor progression by affecting the recruitment, activation, and differentiation of multiple immune cells. Therefore, Piezo1 has excellent potential as an anti-tumor target. The article reviews the diverse physiological functions of Piezo1 in the human body and its major cellular pathways during disease development, and describes in detail the specific mechanisms by which Piezo1 affects the malignant behavior of tumors and its recent progress as a new target for tumor therapy, providing new perspectives for exploring more potential effects on physiological functions and its application in tumor therapy.

Role of Protonation States in the Stability of Molecular Dynamics Simulations of High-Resolution Membrane Protein Structures.

Classical molecular dynamics (MD) simulations provide unmatched spatial and time resolution of protein structure and function. However, the accuracy of MD simulations often depends on the quality of force field parameters and the time scale of sampling. Another limitation of conventional MD simulations is that the protonation states of titratable amino acid residues remain fixed during simulations, even though protonation state changes coupled to conformational dynamics are central to protein function. Due to the uncertainty in selecting protonation states, classical MD simulations are sometimes performed with all amino acids modeled in their standard charged states at pH 7. Here, we performed and analyzed classical MD simulations on high-resolution cryo-EM structures of two large membrane proteins that transfer protons by catalyzing protonation/deprotonation reactions. In simulations performed with titratable amino acids modeled in their standard protonation (charged) states, the structure diverges far from its starting conformation. In comparison, MD simulations performed with predetermined protonation states of amino acid residues reproduce the structural conformation, protein hydration, and protein-water and protein-protein interactions of the structure much better. The results support the notion that it is crucial to perform basic protonation state calculations, especially on structures where protonation changes play an important functional role, prior to the launch of any conventional MD simulations. Furthermore, the combined approach of fast protonation state prediction and MD simulations can provide valuable information about the charge states of amino acids in the cryo-EM sample. Even though accurate prediction of protonation states in proteinaceous environments currently remains a challenge, we introduce an approach of combining pKa prediction with cryo-EM density map analysis that helps in improving not only the protonation state predictions but also the atomic modeling of density data.

Identification of Ligands for Ion Channels: TRPM2.

Transient receptor potential melastatin 2 (TRPM2) is a calcium-permeable, nonselective cation channel with a widespread distribution throughout the body. It is involved in many pathological and physiological processes, making it a potential therapeutic target for various diseases, including Alzheimer's disease, Parkinson's disease, and cancers. New analytical techniques are beneficial for gaining a deeper understanding of its involvement in disease pathogenesis and for advancing the drug discovery for TRPM2-related diseases. In this work, we present the application of collision-induced affinity selection mass spectrometry (CIAS-MS) for the direct identification of ligands binding to TRPM2. CIAS-MS circumvents the need for high mass detection typically associated with mass spectrometry of large membrane proteins. Instead, it focuses on the detection of small molecules dissociated from the ligand-protein-detergent complexes. This affinity selection approach consolidates all affinity selection steps within the mass spectrometer, resulting in a streamlined process. We showed the direct identification of a known TRPM2 ligand dissociated from the protein-ligand complex. We demonstrated that CIAS-MS can identify binding ligands from complex mixtures of compounds and screened a compound library against TRPM2. We investigated the impact of voltage increments and ligand concentrations on the dissociation behavior of the binding ligand, revealing a dose-dependent relationship.

Transmembrane Helices 7 and 8 Confer Aggregation Sensitivity to the Cystic Fibrosis Transmembrane Conductance Regulator.

The Cystic Fibrosis Transmembrane Conductance Regulator (CFTR) is a large multi-spanning membrane protein that is susceptible to misfolding and aggregation. We have identified here the region responsible for this instability. Temperature-induced aggregation of C-terminally truncated versions of CFTR demonstrated that all truncations up to the second transmembrane domain (TMD2), including the R region, largely resisted aggregation. Limited proteolysis identified a folded structure that was prone to aggregation and consisted of TMD2 and at least part of the Regulatory Region R. Only when both TM7 (TransMembrane helix 7) and TM8 were present, TMD2 fragments became as aggregation-sensitive as wild-type CFTR, in line with increased thermo-instability of late CFTR nascent chains and in silico prediction of aggregation propensity. In accord, isolated TMD2 was degraded faster in cells than isolated TMD1. We conclude that TMD2 extended at its N-terminus with part of the R region forms a protease-resistant structure that induces heat instability in CFTR and may be responsible for its limited intracellular stability.

The Dysferlin C2A Domain Binds PI(4,5)P2 and Penetrates Membranes

Dysferlin is a large membrane protein found most prominently in striated muscle. Loss of dysferlin activity is associated with reduced exocytosis, abnormal intracellular Ca2+ and the muscle diseases limb-girdle muscular dystrophy and Miyoshi myopathy. The cytosolic region of dysferlin consists of seven C2 domains with mutations in the C2A domain at the N-terminus resulting in pathology. Despite the importance of Ca2+ and membrane binding activities of the C2A domain for dysferlin function, the mechanism of the domain remains poorly characterized. In this study we find that the C2A domain preferentially binds membranes containing PI(4,5)P2 through an interaction mediated by residues Y23, K32, K33, and R77 on the concave face of the domain. We also found that subsequent to membrane binding, the C2A domain inserts residues on the Ca2+ binding loops into the membrane. Analysis of solution NMR measurements indicate that the domain inhabits two distinct structural states, with Ca2+ shifting the population between states towards a more rigid structure with greater affinity for PI(4,5)P2. Based on our results, we propose a mechanism where Ca2+ converts C2A from a structurally dynamic, low PI(4,5)P2 affinity state to a high affinity state that targets dysferlin to PI(4,5)P2 enriched membranes through interaction with Tyr23, K32, K33, and R77. Binding also involves changes in lipid packing and insertion by the third Ca2+ binding loop of the C2 domain into the membrane, which would contribute to dysferlin function in exocytosis and Ca2+ regulation.

MIC26 and MIC27 are bona fide subunits of the MICOS complex in mitochondria and do not exist as glycosylated apolipoproteins.

Impairments of mitochondrial functions are linked to human ageing and pathologies such as cancer, cardiomyopathy, neurodegeneration and diabetes. Specifically, aberrations in ultrastructure of mitochondrial inner membrane (IM) and factors regulating them are linked to diabetes. The development of diabetes is connected to the 'Mitochondrial Contact Site and Cristae Organising System' (MICOS) complex which is a large membrane protein complex defining the IM architecture. MIC26 and MIC27 are homologous apolipoproteins of the MICOS complex. MIC26 has been reported as a 22 kDa mitochondrial and a 55 kDa glycosylated and secreted protein. The molecular and functional relationship between these MIC26 isoforms has not been investigated. In order to understand their molecular roles, we depleted MIC26 using siRNA and further generated MIC26 and MIC27 knockouts (KOs) in four different human cell lines. In these KOs, we used four anti-MIC26 antibodies and consistently detected the loss of mitochondrial MIC26 (22 kDa) and MIC27 (30 kDa) but not the loss of intracellular or secreted 55 kDa protein. Thus, the protein assigned earlier as 55 kDa MIC26 is nonspecific. We further excluded the presence of a glycosylated, high-molecular weight MIC27 protein. Next, we probed GFP- and myc-tagged variants of MIC26 with antibodies against GFP and myc respectively. Again, only the mitochondrial versions of these tagged proteins were detected but not the corresponding high-molecular weight MIC26, suggesting that MIC26 is indeed not post-translationally modified. Mutagenesis of predicted glycosylation sites in MIC26 also did not affect the detection of the 55 kDa protein band. Mass spectrometry of a band excised from an SDS gel around 55 kDa could not confirm the presence of any peptides derived from MIC26. Taken together, we conclude that both MIC26 and MIC27 are exclusively localized in mitochondria and that the observed phenotypes reported previously are exclusively due to their mitochondrial function.

Unraveling the mechanisms of cardiolipin function: The role of oxidative polymerization of unsaturated acyl chains

Cardiolipin is a unique phospholipid of the inner mitochondrial membrane (IMM) as well as in bacteria. It performs several vital functions such as resisting osmotic rupture and stabilizing the supramolecular structure of large membrane proteins, like ATP synthases and respirasomes. The process of cardiolipin biosynthesis results in the production of immature cardiolipin. A subsequent step is required for its maturation when its acyl groups are replaced with unsaturated acyl chains, primarily linoleic acid. Linoleic acid is the major fatty acid of cardiolipin across all organs and tissues, except for the brain. Linoleic acid is not synthesized by mammalian cells. It has the unique ability to undergo oxidative polymerization at a moderately accelerated rate compared to other unsaturated fatty acids. This property can enable cardiolipin to form covalently bonded net-like structures essential for maintaining the complex geometry of the IMM and gluing the quaternary structure of large IMM protein complexes. Unlike triglycerides, phospholipids possess only two covalently linked acyl chains, which constrain their capacity to develop robust and complicated structures through oxidative polymerization of unsaturated acyl chains. Cardiolipin, on the other hand, has four fatty acids at its disposal to form covalently bonded polymer structures. Despite its significance, the oxidative polymerization of cardiolipin has been overlooked due to the negative perception surrounding biological oxidation and methodological difficulties. Here, we discuss an intriguing hypothesis that oxidative polymerization of cardiolipin is essential for the structure and function of cardiolipin in the IMM in physiological conditions. In addition, we highlight current challenges associated with the identification and characterization of oxidative polymerization of cardiolipin in vivo. Altogether, the study provides a better understanding of the structural and functional role of cardiolipin in mitochondria.

Environmentally Ultrasensitive Fluorine Probe to Resolve Protein Conformational Ensembles by 19F NMR and Cryo-EM.

Limited chemical shift dispersion represents a significant barrier to studying multistate equilibria of large membrane proteins by 19F NMR. We describe a novel monofluoroethyl 19F probe that dramatically increases the chemical shift dispersion. The improved conformational sensitivity and line shape enable the detection of previously unresolved states in one-dimensional (1D) 19F NMR spectra of a 134 kDa membrane transporter. Changes in the populations of these states in response to ligand binding, mutations, and temperature correlate with population changes of distinct conformations in structural ensembles determined by single-particle cryo-electron microscopy (cryo-EM). Thus, 19F NMR can guide sample preparation to discover and visualize novel conformational states and facilitate image analysis and three-dimensional (3D) classification.

DNA-nanostructure-templated assembly of planar and curved lipid-bilayer membranes.

Lipid-bilayer nanodiscs and liposomes have been developed to stabilize membrane proteins in order to study their structures and functions. Nanodiscs are detergent-free, water-soluble, and size-controlled planar phospholipid-bilayer platforms. On the other hand, liposomes are curved phospholipid-bilayer spheres with an aqueous core used as drug delivery systems and model membrane platforms for studying cellular activities. A long-standing challenge is the generation of a homogenous and monodispersed lipid-bilayer system with a very wide range of dimensions and curvatures (elongation, bending, and twisting). A DNA-origami template provides a way to control the shapes, sizes, and arrangements of lipid bilayers via enforcing the assembly of lipid bilayers within the cavities created by DNA nanostructures. Here, we provide a concise overview and discuss how to design planar and curved lipid-bilayer membranes by using DNA-origami nanostructures as templates. Finally, we will discuss the potential applications of DNA-origami nanostructures in the structural and functional studies of large membrane proteins and their complexes.

Purification of Active Photosystem I-Light Harvesting Complex I from Plant Tissues.

This method is used to isolate Photosystem I (PSI) together with the Light Harvesting Complex I (LHCI), its native antenna, from plants. PSI-LHCI is a large membrane protein complex coordinating hundreds of light harvesting and electron transport factors and is the most efficient light harvesting system found in nature. Photons absorbed by the four LHCA antenna proteins that make up LHCI are transferred through excitonic interaction to the PSI core reaction center and are used to facilitate light-driven charge separation across the thylakoid membrane, providing reducing power and energy for carbon fixation in photoautotrophic organisms. The high quantum efficiency of PSI makes this complex an excellent model to study light-driven energy transfer. In this protocol, plant tissue is mechanically homogenized, and the chloroplasts are separated from the bulk cellular debris by filtration and centrifugation. The isolated chloroplasts are then osmotically lysed, and the thylakoid membranes are recovered via centrifugation and solubilized using the detergent n-dodecyl-beta-maltoside. The solubilized material is loaded onto an anion exchange column to collect most of the chlorophyll-containing complexes. Larger complexes are precipitated from the solution, resuspended in a small volume, and loaded on sucrose gradients to separate the major chlorophyll-containing complexes. The resulting sucrose gradient fractions are characterized to identify the band of interest containing PSI-LHCI. This protocol is highly similar to the protocol used in the crystallization of plant PSI-LHCI with some simplifications and relies on methods developed over the years in the lab of Nathan Nelson.

High-level expression of voltage-gated sodium channels in thymidine-arrested human embryonic kidney cells.

The benefit of using heterologous mammalian cell lines such as human embryonic kidney (HEK293) cells for expressing voltage-gated ion channels includes efficient translation and processing of such large membrane proteins. Still, the biophysical assessment of these channels is hindered by their low expression and reduced localization on the surface of a cell. Having a higher expression level of ion channels is crucial for investigating mutations that cause a decrease in the current amplitude as well as for measuring tiny currents such as gating currents and gating pore currents that are in the range of 0.1 to 1% of the peak current. We aimed at improving the efficiency of ion channel expression and enhancing their localization in cell membranes by modifying the protocol of cell cycle arrest at the G1/S boundary using thymidine to direct the cells into increasing the expression machinery during viral transduction of mammalian voltage-gated sodium channels. The combination of cell arrest at the G1/S boundary and mammalian baculovirus system (BacMam) resulted in a 5-fold increase in the current density of different mammalian voltage-gated sodium channels in HEK293 cells compared to the current literature. By applying this modified protocol for expression, we succeeded to measure a pathogenic gating pore current in HEK293 cells of ∼0.5% of the central pore current induced by an autism-related mutation in the R2 gating charge in Domain II (R853Q) of Nav1.2, which is blocked in a voltage-dependent manner when the voltage sensors activate. These results show that our protocol is a universal one that can be applied on different voltage-gated ion channels and possibly other membrane proteins. This ultimately will pave the way to decipher the structure and function of ion channels and their association with ion channelopathies and other diseases.

Acting on the CFTR Membrane-Spanning Domains Interface Rescues Some Misfolded Mutants.

ABC transporters are large membrane proteins sharing a complex architecture, which comprises two nucleotide-binding domains (NBDs) and two membrane-spanning domains (MSDs). These domains are susceptible to mutations affecting their folding and assembly. In the CFTR (ABCC7) protein, a groove has been highlighted in the MSD1 at the level of the membrane inner leaflet, containing both multiple mutations affecting folding and a binding site for pharmaco-chaperones that stabilize this region. This groove is also present in ABCB proteins, however it is covered by a short elbow helix, while in ABCC proteins it remains unprotected, due to a lower position of the elbow helix in the presence of the ABCC-specific lasso motif. Here, we identified a MSD1 second-site mutation located in the vicinity of the CFTR MSD1 groove that partially rescued the folding defect of cystic fibrosis causing mutations located within MSD1, while having no effect on the most frequent mutation, F508del, located within NBD1. A model of the mutated protein 3D structure suggests additional interaction between MSD1 and MSD2, strengthening the assembly at the level of the MSD intracellular loops. Altogether, these results provide insightful information in understanding key features of the folding and function of the CFTR protein in particular, and more generally, of type IV ABC transporters.

A Method to Monitor Activity of SARS-CoV-2 Nsp3 from Cells.

SARS-CoV-2 protease Nsp3 is a therapeutic target for developing anti-SARS-CoV-2 drugs. Nsp3 is a large multi-spanning membrane protein, and its characterization in vitro has been challenging. Here we describe an in vitro assay to characterize the biochemical activity of full-length Nsp3 isolated from cells. The assay can be used to evaluate Nsp3 inhibitors.

3,4-Bis(hydroxymethyl)hexane-1,6-diol-based Maltosides (HDMs) for Membrane-Protein Study: Importance of Detergent Rigidity-Flexibility Balance in Protein Stability.

Detergents have been major contributors to membrane-protein structural study for decades. However, membrane proteins solubilized in conventional detergents tend to aggregate or denature over time. Stability of large eukaryotic membrane proteins with complex structures tends to be particularly poor, necessitating development of novel detergents with improved properties. Here, we prepared a novel class of detergents, designated 3,4-bis(hydroxymethyl)hexane-1,6-diol-based maltosides (HDMs). When tested on three membrane proteins, including two G-protein-coupled receptors (GPCRs), the new detergents displayed significantly better behaviors compared with DDM. Moreover, the HDMs were superior or comparable to LMNG, an amphiphile widely used for GPCR structural study. An optimal balance of detergent rigidity vs. flexibility of the HDMs is likely responsible for their favorable behaviors toward membrane-protein stability. Thus, the current study not only introduces the HDMs, with significant potential for membrane-protein structural study, but also suggests a useful guideline for designing novel detergents for membrane-protein research.

Importance of molecular dynamics equilibrium protocol on protein-lipid interaction near channel pore

Multiscale molecular dynamics simulations using Martini coarse-grained (CG) and all-atom (AA) force fields are commonly used in membrane protein studies. In particular, reverse mapping an equilibrated CG model to an AA model offers an efficient way for preparing large membrane protein systems with complex protein shapes and lipid compositions. Here, we report that this hybrid CG-equilibrium-AA-production protocol may artificially increase lipid density and decrease hydration in ion channel pores walled with transmembrane gaps. To understand the origin of this conundrum, we conducted replicas of CG, AA, and CG reverse-mapped AA simulations of the pore domain of the mechanosensitive Piezo1 channel in a nonconducting conformation. Lipid/water density analysis and free energy calculations reveal that the lack of initial pore hydration allows excessive lipids to enter the upper pore lumen through gaps between pore helices during CG simulation. Due to the mismatch between CG and AA lipid kinetics, these pore lipids remain trapped in the subsequent AA simulations, despite unfavorable binding free energy. We tested several CG equilibrium protocols and found that a protocol restraining the whole lipid produces pore hydration consistent with AA results, thus eliminating this artifact for further studies of lipid gating and protein-lipid interactions.

COGRIMEN: Coarse-Grained Method for Modeling of Membrane Proteins in Implicit Environments.

The presented methodology is based on coarse-grained representation of biomolecules in implicit environments and is designed for the molecular dynamics simulations of membrane proteins and their complexes. The membrane proteins are not only found in the cell membrane but also in all membranous compartments of the cell: Golgi apparatus, mitochondria, endosomes and lysosomes, and they usually form large complexes. To investigate such systems the methodology is proposed based on two independent approaches combining the coarse-grained MARTINI model for proteins and the effective energy function to mimic the water/membrane environments. The latter is based on the implicit environment developed for all-atom simulations in the IMM1 method. The force field solvation parameters for COGRIMEN were initially calculated from IMM1 all-atom parameters and then optimized using Genetic Algorithms. The new methodology was tested on membrane proteins, their complexes and oligomers. COGRIMEN method is implemented as a patch for NAMD program and can be useful for fast and brief studies of large membrane protein complexes.

Electrostatic interactions between single arginine and phospholipids modulate physiological properties of sarcoplasmic reticulum Ca2+-ATPase

Arg324 of sarcoplasmic reticulum Ca2+-ATPase forms electrostatic interactions with the phosphate moiety of phospholipids in most reaction states, and a hydrogen bond with Tyr122 in other states. Using site-directed mutagenesis, we explored the functional roles of Arg324 interactions, especially those with lipids, which at first glance might seem too weak to modulate the function of such a large membrane protein. The hydrogen bond forms transiently and facilitates Ca2+ binding from the cytoplasmic side. The contributions of the electrostatic interactions to the reaction steps were quantified using a rate vs activity coefficient plot. We found that the interaction between Arg324 and lipids decreases the affinity for luminal Ca2+. The transformation rate of the phosphoenzyme intermediate is facilitated by the electrostatic interactions, and the function of these interactions depends not only on the type but also on the composition of the phospholipids. The properties observed in microsomes could not be reproduced with any single phospholipid, but with a mixture of phospholipids that mimics the native membrane. These results suggest the importance of swapping of the lipid partners of different headgroups in the reaction step. This study shows that Arg324 plays a role in the reaction cycle via complex intra-protein and protein-lipid interactions.