High-Resolution NMR Determination of the Dynamic Structure of Membrane Proteins

Pitfalls in measurements of R1 relaxation rates of protein backbone 15N nuclei

The dynamics of the backbone and side-chains of protein are routinely studied by interpreting experimentally determined 15N spin relaxation rates. R1(15N), the longitudinal relaxation rate, reports on fast motions and encodes, together with the transverse relaxation R2, structural information about the shape of the molecule and the orientation of the amide bond vectors in the internal diffusion frame. Determining error-free 15N longitudinal relaxation rates remains a challenge for small, disordered, and medium-sized proteins. Here, we show that mono-exponential fitting is sufficient, with no statistical preference for bi-exponential fitting up to 800 MHz. A detailed comparison of the TROSY and HSQC techniques at medium and high fields showed no statistically significant differences. The least error-prone DD/CSA interference removal technique is the selective inversion of amide signals while avoiding water resonance. The exchange of amide with solvent deuterons appears to affect the rate R1 of solvent-exposed amides in all fields tested and in each DD/CSA interference removal technique in a statistically significant manner. In summary, the most accurate R1(15N) rates in proteins are achieved by selective amide inversion, without the addition of D2O. Importantly, at high magnetic fields stronger than 800 MHz, when non-mono-exponential decay is involved, it is advisable to consider elimination of the shortest delays (typically up to 0.32 s) or bi-exponential fitting.

New insights into the structure and dynamics of the TOM complex in mitochondria

To date, there is no general physical model of the mechanism by which unfolded polypeptide chains with different properties are imported into the mitochondria. At the molecular level, it is still unclear how transit polypeptides approach, are captured by the protein translocation machinery in the outer mitochondrial membrane, and how they subsequently cross the entropic barrier of a protein translocation pore to enter the intermembrane space. This deficiency has been due to the lack of detailed structural and dynamic information about the membrane pores. In this review, we focus on the recently determined sub-nanometer cryo-EM structures and our current knowledge of the dynamics of the mitochondrial two-pore outer membrane protein translocation machinery (TOM core complex), which provide a starting point for addressing the above questions. Of particular interest are recent discoveries showing that the TOM core complex can act as a mechanosensor, where the pores close as a result of interaction with membrane-proximal structures. We highlight unusual and new correlations between the structural elements of the TOM complexes and their dynamic behavior in the membrane environment.

Sparse pseudocontact shift NMR data obtained from a non-canonical amino acid-linked lanthanide tag improves integral membrane protein structure prediction

A single experimental method alone often fails to provide the resolution, accuracy, and coverage needed to model integral membrane proteins (IMPs). Integrating computation with experimental data is a powerful approach to supplement missing structural information with atomic detail. We combine RosettaNMR with experimentally-derived paramagnetic NMR restraints to guide membrane protein structure prediction. We demonstrate this approach using the disulfide bond formation protein B (DsbB), an α-helical IMP. Here, we attached a cyclen-based paramagnetic lanthanide tag to an engineered non-canonical amino acid (ncAA) using a copper-catalyzed azide-alkyne cycloaddition (CuAAC) click chemistry reaction. Using this tagging strategy, we collected 203 backbone HN pseudocontact shifts (PCSs) for three different labeling sites and used these as input to guide de novo membrane protein structure prediction protocols in Rosetta. We find that this sparse PCS dataset combined with 44 long-range NOEs as restraints in our calculations improves structure prediction of DsbB by enhancements in model accuracy, sampling, and scoring. The inclusion of this PCS dataset improved the Cα-RMSD transmembrane segment values of the best-scoring and best-RMSD models from 9.57 Å and 3.06 Å (no NMR data) to 5.73 Å and 2.18 Å, respectively.

E as in Enigma: The Mysterious Role of the Voltage-Dependent Anion Channel Glutamate E73

The voltage-dependent anion channel (VDAC) is the main passageway for ions and metabolites over the outer mitochondrial membrane. It was associated with many physiological processes, including apoptosis and modulation of intracellular Ca²⁺ signaling. The protein is formed by a barrel of 19 beta-sheets with an N-terminal helix lining the inner pore. Despite its large diameter, the channel can change its selectivity for ions and metabolites based on its open state to regulate transport into and out of mitochondria. VDAC was shown to be regulated by a variety of cellular factors and molecular partners including proteins, lipids and ions. Although the physiological importance of many of these modulatory effects are well described, the binding sites for molecular partners are still largely unknown. The highly symmetrical and sleek structure of the channel makes predictions of functional moieties difficult. However, one residue repeatedly sticks out when reviewing VDAC literature. A glutamate at position 73 (E73) located on the outside of the channel facing the hydrophobic membrane environment was repeatedly proposed to be involved in channel regulation on multiple levels. Here, we review the distinct hypothesized roles of E73 and summarize the open questions around this mysterious residue.

Increased slow dynamics defines ligandability of BTB domains

Efficient determination of protein ligandability, or the propensity to bind small-molecules, would greatly facilitate drug development for novel targets. Ligandability is currently assessed using computational methods that typically consider the static structural properties of putative binding sites or by experimental fragment screening. Here, we evaluate ligandability of conserved BTB domains from the cancer-relevant proteins LRF, KAISO, and MIZ1. Using fragment screening, we discover that MIZ1 binds multiple ligands. However, no ligands are uncovered for the structurally related KAISO or LRF. To understand the principles governing ligand-binding by BTB domains, we perform comprehensive NMR-based dynamics studies and find that only the MIZ1 BTB domain exhibits backbone µs-ms time scale motions. Interestingly, residues with elevated dynamics correspond to the binding site of fragment hits and recently defined HUWE1 interaction site. Our data argue that examining protein dynamics using NMR can contribute to identification of cryptic binding sites, and may support prediction of the ligandability of novel challenging targets.

Deactivation of catalysts in simultaneous reversible and irreversible parahydrogen NMR signal enhancement, and the role of co-ligands in the stabilization of the reversible method

Signal Amplification by Reversible Exchange (SABRE) and hydrogeneable Parahydrogen Induced Polarization (hPHIP) can enhance weak NMR signals, and thus increase the range of NMR applications. Here, using an N-heterocyclic carbene Ir-based catalyst, simultaneous SABRE and hPHIP was achieved for the compound with an N-donor site and an acetylene triple bond. It was demonstrated that the interplay between SABRE and hPHIP can be manipulated. Specifically, it was found that the hPHIP effect could be almost completely suppressed, while stable SABRE was observed in subsequent consecutive experiments. The presented results have the potential to increase the numbers of parahydrogen hyperpolarizable molecules.

Structure and Gating Behavior of the Human Integral Membrane Protein VDAC1 in a Lipid Bilayer

The voltage-dependent anion channel (VDAC), the most abundant protein in the outer mitochondrial membrane, is responsible for the transport of all ions and metabolites into and out of mitochondria. Larger than any of the β-barrel structures determined to date by magic-angle spinning (MAS) NMR, but smaller than the size limit of cryo-electron microscopy (cryo-EM), VDAC1’s 31 kDa size has long been a bottleneck in determining its structure in a near-native lipid bilayer environment. Using a single two-dimensional (2D) crystalline sample of human VDAC1 in lipids, we applied proton-detected fast magic-angle spinning NMR spectroscopy to determine the arrangement of β strands. Combining these data with long-range restraints from a spin-labeled sample, chemical shift-based secondary structure prediction, and previous MAS NMR and atomic force microscopy (AFM) data, we determined the channel’s structure at a 2.2 Å root-mean-square deviation (RMSD). The structure, a 19-stranded β-barrel, with an N-terminal α-helix in the pore is in agreement with previous data in detergent, which was questioned due to the potential for the detergent to perturb the protein’s functional structure. Using a quintuple mutant implementing the channel’s closed state, we found that dynamics are a key element in the protein’s gating behavior, as channel closure leads to the destabilization of not only the C-terminal barrel residues but also the α2 helix. We showed that cholesterol, previously shown to reduce the frequency of channel closure, stabilizes the barrel relative to the N-terminal helix. Furthermore, we observed channel closure through steric blockage by a drug shown to selectively bind to the channel, the Bcl2-antisense oligonucleotide G3139.

Regulation of Mitochondrial Respiration by VDAC Is Enhanced by Membrane-Bound Inhibitors with Disordered Polyanionic C-Terminal Domains

The voltage-dependent anion channel (VDAC) is the primary regulating pathway of water-soluble metabolites and ions across the mitochondrial outer membrane. When reconstituted into lipid membranes, VDAC responds to sufficiently large transmembrane potentials by transitioning to gated states in which ATP/ADP flux is reduced and calcium flux is increased. Two otherwise unrelated cytosolic proteins, tubulin, and α-synuclein (αSyn), dock with VDAC by a novel mechanism in which the transmembrane potential draws their disordered, polyanionic C-terminal domains into and through the VDAC channel, thus physically blocking the pore. For both tubulin and αSyn, the blocked state is observed at much lower transmembrane potentials than VDAC gated states, such that in the presence of these cytosolic docking proteins, VDAC's sensitivity to transmembrane potential is dramatically increased. Remarkably, the features of the VDAC gated states relevant for bioenergetics-reduced metabolite flux and increased calcium flux-are preserved in the blocked state induced by either docking protein. The ability of tubulin and αSyn to modulate mitochondrial potential and ATP production in vivo is now supported by many studies. The common physical origin of the interactions of both tubulin and αSyn with VDAC leads to a general model of a VDAC inhibitor, facilitates predictions of the effect of post-translational modifications of known inhibitors, and points the way toward the development of novel therapeutics targeting VDAC.

VDAC regulation of mitochondrial calcium flux: From channel biophysics to disease

Voltage-dependent anion channel (VDAC), the most abundant mitochondrial outer membrane protein, is important for a variety of mitochondrial functions including metabolite exchange, calcium transport, and apoptosis. While VDAC's role in shuttling metabolites between the cytosol and mitochondria is well established, there is a growing interest in understanding the mechanisms of its regulation of mitochondrial calcium transport. Here we review the current literature on VDAC's role in calcium signaling, its biophysical properties, physiological function, and pathology focusing on its importance in cardiac diseases. We discuss the specific biophysical properties of the three VDAC isoforms in mammalian cells-VDAC 1, 2, and 3-in relationship to calcium transport and their distinct roles in cell physiology and disease. Highlighting the emerging evidence that cytosolic proteins interact with VDAC and regulate its calcium permeability, we advocate for continued investigation into the VDAC interactome at the contact sites between mitochondria and organelles and its role in mitochondrial calcium transport.

Structure, gating and interactions of the voltage-dependent anion channel

The voltage-dependent anion channel (VDAC) is one of the most highly abundant proteins found in the outer mitochondrial membrane, and was one of the earliest discovered. Here we review progress in understanding VDAC function with a focus on its structure, discussing various models proposed for voltage gating as well as potential drug targets to modulate the channel's function. In addition, we explore the sensitivity of VDAC structure to variations in the membrane environment, comparing DMPC-only, DMPC with cholesterol, and near-native lipid compositions, and use magic-angle spinning NMR spectroscopy to locate cholesterol on the outside of the β-barrel. We find that the VDAC protein structure remains unchanged in different membrane compositions, including conditions with cholesterol.

The Structural Basis for Low Conductance in the Membrane Protein VDAC upon β-NADH Binding and Voltage Gating

The voltage-dependent anion channel (VDAC) forms the primary diffusion pore of the outer mitochondrial membrane. In its apo form, VDAC adopts an open conformation with high conductance. States of lower conductance can be induced by ligand binding or the application of voltage. Here, we clarify at the atomic level how β-NADH binding leads to a low-conductance state and characterize the role of the VDAC N-terminal helix in voltage gating. A high-resolution NMR structure of human VDAC-1 with bound NADH, combined with molecular dynamics simulation show that β-NADH binding reduces the pore conductance sterically without triggering a structural change. Electrophysiology recordings of crosslinked protein variants and NMR relaxation experiments probing different time scales show that increased helix dynamics is present in the open state and that motions of the N-terminal helices are involved in the VDAC voltage gating mechanism.

Towards the functional high-resolution coordination chemistry of blood plasma human serum albumin

Human serum albumin (HSA) is a monomeric, globular, multi-carrier and the most abundant protein in the blood. HSA displays multiple ligand binding sites with extraordinary binding capacity for a wide range of ions and molecules. For decades, HSA's ability to bind to various ligands has led many scientists to study its physiological properties and protein structure; indeed, a better understanding of HSA-ligand interactions in human blood, at the atomic level, will likely foster the development of more potent, and overall more performant, diagnostic and therapeutic tools against serious human disorders such as diabetes, cardiovascular disorders, and cancer. Here, we present a concise overview of the current knowledge of HSA's structural characteristics, and its coordination chemistry with transition metal ions, within the scope and limitations of current techniques and biophysical methods to reach atomic resolution in solution and in blood serum. We also highlight the overwhelming need of a detailed atomistic understanding of HSA dynamic structures and interactions that are transient, weak, multi-site and multi-step, and allosterically affected by each other. Considering the fact that HSA is a current clinical tool for drug delivery systems and a potential contender as molecular cargo and nano-vehicle used in biophysical, clinical and industrial fields, we underline the emerging need for novel approaches to target the dynamic functional coordination chemistry of the human blood serum albumin in solution, at the atomic level.

Assessing the role of residue E73 and lipid headgroup charge in VDAC1 voltage gating

The voltage-dependent anion channel (VDAC) is the most abundant protein of the mitochondrial outer membrane (MOM) where it regulates transport of ions and metabolites in and out of the organelle. VDAC function is extensively studied in a lipid bilayer system that allows conductance monitoring of reconstituted channels under applied voltage. The process of switching from a high-conductance state, open to metabolites, to a variety of low-conducting states, which excludes metabolite transport, is termed voltage gating and the mechanism remains poorly understood. Recent studies have implicated the involvement of the membrane-solvated residue E73 in the gating process through β-barrel destabilization. However, there has been no direct experimental evidence of E73 involvement in VDAC1 voltage gating. Here, using electrophysiology measurements, we exclude the involvement of E73 in murine VDAC1 (mVDAC1) voltage gating process. With an established protocol of assessing voltage gating of VDACs reconstituted into planar lipid membranes, we definitively show that mVDAC1 gating properties do not change when E73 is replaced by either a glutamine or an alanine. We further demonstrate that cholesterol has no effect on mVDAC1 gating characteristics, though it was shown that E73 is coordinating residue in the cholesterol binding site. In contrast, we found a pronounced gating effect based on the charge of the phospholipid headgroup, where the positive charge stimulates and negative charge suppresses gating. These findings call for critical evaluation of the existing models of VDAC gating and contribute to our understanding of VDAC's role in control of MOM permeability and regulation of mitochondrial respiration and metabolism.

Probing Membrane Protein Insertion into Lipid Bilayers by Solid-State NMR

Determination of the environment surrounding a protein is often key to understanding its function and can also be used to infer the structural properties of the protein. By using proton-detected solid-state NMR, we show that reduced spin diffusion within the protein under conditions of fast magic-angle spinning, high magnetic field, and sample deuteration allows the efficient measurement of site-specific exposure to mobile water and lipids. We demonstrate this site specificity on two membrane proteins, the human voltage dependent anion channel, and the alkane transporter AlkL from Pseudomonas putida. Transfer from lipids is observed selectively in the membrane spanning region, and an average lipid-protein transfer rate of 6 s^-1 was determined for residues protected from exchange. Transfer within the protein, as tracked in the ¹⁵ N-¹ H 2D plane, was estimated from initial rates and found to be in a similar range of about 8 to 15 s^-1 for several resolved residues, explaining the site specificity.

Lipid composition and salt concentration as regulatory factors of the anion selectivity of VDAC studied by coarse-grained molecular dynamics simulations

The voltage-dependent anion channel (VDAC) is a mitochondrial outer membrane protein whose fundamental function is to facilitate and regulate the flow of metabolites between the cytosol and the mitochondrial intermembrane space. Using coarse-grained molecular dynamics simulations, we investigated the dependence of VDAC selectivity towards small inorganic anions on two factors: the ionic strength and the lipid composition. In agreement with experimental data we found that VDAC becomes less anion selective with increasing salt concentration due to the screening of a few basic residues that point into the pore lumen. The molecular dynamics simulations provide insight into the regulation mechanism of VDAC selectivity by the composition in the lipid membrane and suggest that the ion distribution is differently modulated by POPE compared to the POPC bilayer. This occurs through the more persistent interactions of acidic residues located at both rims of the β-barrel with head groups of POPE which in turn impact the electrostatic potential and thereby the selectivity of the pore. This mechanism occurs not only in POPE single component membranes but also in a mixed POPE/POPC bilayer by an enrichment of POPE over POPC lipids on the surface of VDAC. Thus we show here that computationally-inexpensive coarse-grained simulations are able to capture, in a semi-quantitative way, essential features of VDAC anion selectivity and could pave the way toward a molecular level understanding of metabolite transport in natural membranes.

Mitochondrial Proteolipid Complexes of Creatine Kinase

Isoforms of creatine kinase (CK) generate and use phosphocreatine, a concentrated and highly diffusible cellular "high energy" intermediate, for the main purpose of energy buffering and transfer in order to maintain cellular energy homeostasis. The mitochondrial CK isoform (mtCK) localizes to the mitochondrial intermembrane and cristae space, where it assembles into peripherally membrane-bound, large cuboidal homooctamers. These are part of proteolipid complexes wherein mtCK directly interacts with cardiolipin and other anionic phospholipids, as well as with the VDAC channel in the outer membrane. This leads to a stabilization and cross-linking of inner and outer mitochondrial membrane, forming so-called contact sites. Also the adenine nucleotide translocator of the inner membrane can be recruited into these proteolipid complexes, probably mediated by cardiolipin. The complexes have functions mainly in energy transfer to the cytosol and stimulation of oxidative phosphorylation, but also in restraining formation of reactive oxygen species and apoptosis. In vitro evidence indicates a putative role of mtCK in mitochondrial phospholipid distribution, and most recently a role in thermogenesis has been proposed. This review summarizes the essential structural and functional data of these mtCK complexes and describes in more detail the more recent advances in phospholipid interaction, thermogenesis, cancer and evolution of mtCK.

Ten Years of High Resolution Structural Research on the Voltage Dependent Anion Channel (VDAC)-Recent Developments and Future Directions

Mitochondria are evolutionarily related to Gram-negative bacteria and both comprise two membrane systems with strongly differing protein composition. The major protein in the outer membrane of mitochondria is the voltage-dependent anion channel (VDAC), which mediates signal transmission across the outer membrane but also the exchange of metabolites, most importantly ADP and ATP. More than 30 years after its discovery three identical high-resolution structures were determined in 2008. These structures show a 19-stranded anti-parallel beta-barrel with an N-terminal helix located inside. An odd number of beta-strands is also shared by Tom40, another member of the VDAC superfamily. This indicates that this superfamily is evolutionarily relatively young and that it has emerged in the context of mitochondrial evolution. New structural information obtained during the last decade on Tom40 can be used to cross-validate the structure of VDAC and vice versa. Interpretation of biochemical and biophysical studies on both protein channels now rests on a solid basis of structural data. Over the past 10 years, complementary structural and functional information on proteins of the VDAC superfamily has been collected from in-organello, in-vitro, and in silico studies. Most of these findings have confirmed the validity of the original structures. This short article briefly reviews the most important advances on the structure and function of VDAC superfamily members collected during the last decade and summarizes how they enhanced our understanding of the channel.

Mitochondrial VDAC1: A Key Gatekeeper as Potential Therapeutic Target

Mitochondria are the key source of ATP that fuels cellular functions, and they are also central in cellular signaling, cell division and apoptosis. Dysfunction of mitochondria has been implicated in a wide range of diseases, including neurodegenerative and cardiac diseases, and various types of cancer. One of the key proteins that regulate mitochondrial function is the voltage-dependent anion channel 1 (VDAC1), the most abundant protein on the outer membrane of mitochondria. VDAC1 is the gatekeeper for the passages of metabolites, nucleotides, and ions; it plays a crucial role in regulating apoptosis due to its interaction with apoptotic and anti-apoptotic proteins, namely members of the Bcl-2 family of proteins and hexokinase. Therefore, regulation of VDAC1 is crucial not only for metabolic functions of mitochondria, but also for cell survival. In fact, multiple lines of evidence have confirmed the involvement of VDAC1 in several diseases. Consequently, modulation or dysregulation of VDAC1 function can potentially attenuate or exacerbate pathophysiological conditions. Understanding the role of VDAC1 in health and disease could lead to selective protection of cells in different tissues and diverse diseases. The purpose of this review is to discuss the role of VDAC1 in the pathogenesis of diseases and as a potentially effective target for therapeutic management of various pathologies.

Crystal structural characterization reveals novel oligomeric interactions of human voltage-dependent anion channel 1

Voltage-dependent anion channel 1 (VDAC1), which is located in the outer mitochondrial membrane, plays important roles in various cellular processes. For example, oligomerization of VDAC1 is involved in the release of cytochrome c to the cytoplasm, leading to apoptosis. However, it is unknown how VDAC1 oligomerization occurs in the membrane. In the present study, we determined high-resolution crystal structures of oligomeric human VDAC1 (hVDAC1) prepared by using an Escherichia coli cell-free protein synthesis system, which avoided the need for denaturation and refolding of the protein. Broad-range screening using a bicelle crystallization method produced crystals in space groups C222 and P22₁ 2₁ , which diffracted to a resolution of 3.10 and 3.15 Å, respectively. Each crystal contained two hVDAC1 protomers in the asymmetric unit. Dimer within the asymmetrical unit of the crystal in space group C222 were oriented parallel, whereas those of the crystal in space group P22₁ 2₁ were oriented anti-parallel. From a model of the crystal in space group C222, which we constructed by using crystal symmetry operators, a heptameric structure with eight patterns of interaction between protomers, including hydrophobic interactions with β-strands, hydrophilic interactions with loop regions, and protein-lipid interactions, was observed. It is possible that by having multiple patterns of interaction, VDAC1 can form homo- or hetero-oligomers not only with other VDAC1 protomers but also with other proteins such as VDAC2, VDAC3 and apoptosis-regulating proteins in the Bcl-2 family.

Applications of NMR to membrane proteins

Membrane proteins present a challenge for structural biology. In this article, we review some of the recent developments that advance the application of NMR to membrane proteins, with emphasis on structural studies in detergent-free, lipid bilayer samples that resemble the native environment. NMR spectroscopy is not only ideally suited for structure determination of membrane proteins in hydrated lipid bilayer membranes, but also highly complementary to the other principal techniques based on X-ray and electron diffraction. Recent advances in NMR instrumentation, spectroscopic methods, computational methods, and sample preparations are driving exciting new efforts in membrane protein structural biology.

Micelles, Bicelles, and Nanodiscs: Comparing the Impact of Membrane Mimetics on Membrane Protein Backbone Dynamics

Detergents are often used to investigate the structure and dynamics of membrane proteins. Whereas the structural integrity seems to be preserved in detergents for many membrane proteins, their functional activity is frequently compromised, but can be restored in a lipid environment. Herein we show with per‐residue resolution that while OmpX forms a stable β‐barrel in DPC detergent micelles, DHPC/DMPC bicelles, and DMPC nanodiscs, the pico‐ to nanosecond and micro‐ to millisecond motions differ substantially between the detergent and lipid environment. In particular for the β‐strands, there is pronounced dynamic variability in the lipid environment, which appears to be suppressed in micelles. This unexpected complex and membrane‐mimetic‐dependent dynamic behavior indicates that the frequent loss of membrane protein activity in detergents might be related to reduced internal dynamics and that membrane protein activity correlates with lipid flexibility.

High-Resolution NMR Determination of the Dynamic Structure of Membrane Proteins

¹⁵N spin‐relaxation rates are demonstrated to provide critical information about the long‐range structure and internal motions of membrane proteins. Combined with an improved calculation method, the relaxation‐rate‐derived structure of the 283‐residue human voltage‐dependent anion channel revealed an anisotropically shaped barrel with a rigidly attached N‐terminal helix. Our study thus establishes an NMR spectroscopic approach to determine the structure and dynamics of mammalian membrane proteins at high accuracy and resolution.