Insight into the mechanism of the influenza A proton channel from a structure in a lipid bilayer

The impact of transmembrane peptides on lipid bilayer structure and mechanics: A study of the transmembrane domain of the influenza A virus M2 protein

Transmembrane peptides play important roles in many biological processes by interacting with lipid membranes. This study investigates how the transmembrane domain of the influenza A virus M2 protein, M2TM, affects the structure and mechanics of model lipid bilayers. Atomic force microscopy (AFM) imaging revealed small decreases in bilayer thickness with increasing peptide concentrations. AFM-based force spectroscopy experiments complemented by theoretical model analysis demonstrated significant decreases in bilayer's Young's modulus (E) and lateral area compressibility modulus (KA). This suggests that M2TM disrupts the cohesive interactions between neighboring lipid molecules, leading to a decrease in both the bilayer's resistance to indentation (E) and its ability to resist lateral compression/expansion (KA). The large decreases in bilayer elastic parameters (i.e., E and KA) contrast with small changes in bilayer thickness, implying that bilayer mechanics are not solely dictated by bilayer thickness in the presence of transmembrane peptides. The observed significant reduction in bilayer mechanical properties suggests a softening effect on the bilayer, potentially facilitating membrane curvature generation, a crucial step for M2-mediated viral budding. In parallel, our Raman spectroscopy revealed small but statistically significant changes in hydrocarbon chain vibrational dynamics, indicative of minor disordering in lipid chain conformation. Our findings provide useful insights into the complex interplay between transmembrane peptides and lipid bilayers, highlighting the significance of peptide-lipid interactions in modulating membrane structure, mechanics, and molecular dynamics.

Advancing the field of viroporins-Structure, function and pharmacology: IUPHAR Review X

Viroporins possess important potential as antiviral targets due to their critical roles during virus life cycles, spanning from virus entry to egress. Although the antiviral amantadine targets the M2 viroporin of influenza A virus, successful progression of other viroporin inhibitors into clinical use remains challenging. These challenges relate in varying proportions to a lack of reliable full-length 3D-structures, difficulties in functionally characterising individual viroporins, and absence of verifiable direct binding between inhibitor and viroporin. This review offers perspectives to help overcome these challenges. We provide a comprehensive overview of the viroporin family, including their structural and functional features, highlighting the moldability of their energy landscapes and actions. To advance the field, we suggest a list of best practices to aspire towards unambiguous viroporin identification and characterisation, along with considerations of potential pitfalls. Finally, we present current and future scenarios of, and prospects for, viroporin targeting drugs.

Conformations of influenza A M2 protein in DOPC/DOPS and E. coli native lipids and proteins

We compared the conformations of the transmembrane domain (TMD) of influenza A M2 (IM2) protein reconstituted in 1,2-dioleoyl-sn-glycero-3-phosphocholine/1,2-dioleoyl-sn-glycero-3-phospho-L-serine (DOPC/DOPS) bilayers to those in isolated Escherichia coli (E. coli) membranes, having preserved its native proteins and lipids. IM2 is a single-pass transmembrane protein known to assemble into a homo-tetrameric proton channel. To represent this channel, we made a construct containing the IM2's TMD region flanked by the juxtamembrane residues. The single cysteine substitution, L43C, of leucine located in the bilayer polar region was paramagnetically tagged with a methanethiosulfonate nitroxide label for the electron spin resonance (ESR) study. For this particular residue, we probed the conformations of the spin-labeled IM2 reconstituted in DOPC/DOPS and isolated E. coli membranes using continuous-wave ESR and double electron-electron resonance (DEER) spectroscopy. The total protein-to-lipid molar ratio spanned the range from 1:230 to 1:10,400. The continuous-wave ESR spectra corresponded to very slow spin-label motion in both environments. In all cases, the DEER data were reconstructed into distance distributions with well-resolved peaks at 1.68 and 2.37 nm in distance and amplitude ratios of 1.41 ± 0.2 and 2:1, respectively. This suggests four nitroxide spin labels located at the corners of a square, indicative of an axially symmetric tetramer. The distance modeling of DEER data with molecular modeling software applied to the NMR molecular structures (PDB: 2L0J) confirmed the symmetry and closed state of the C-terminal exit pore of the IM2 TMD tetramer in agreement with the model. Thus, we can conclude that, under conditions of pH 7.4 used in this study, IM2 TMD has similar conformations in model lipid bilayers and membranes made of native E. coli lipids and proteins of comparable thickness and fluidity, notwithstanding the complexity of the E. coli membranes caused by their lipid diversity and the abundance of integral and peripheral membrane proteins.

Membrane curvature sensing and symmetry breaking of the M2 proton channel from Influenza A

The M2 proton channel aids in the exit of mature influenza viral particles from the host plasma membrane through its ability to stabilize regions of high negative Gaussian curvature (NGC) that occur at the neck of budding virions. The channels are homo-tetramers that contain a cytoplasm-facing amphipathic helix (AH) that is necessary and sufficient for NGC generation; however, constructs containing the transmembrane spanning helix, which facilitates tetramerization, exhibit enhanced curvature generation. Here, we used all-atom molecular dynamics (MD) simulations to explore the conformational dynamics of M2 channels in lipid bilayers revealing that the AH is dynamic, quickly breaking the fourfold symmetry observed in most structures. Next, we carried out MD simulations with the protein restrained in four- and twofold symmetric conformations to determine the impact on the membrane shape. While each pattern was distinct, all configurations induced pronounced curvature in the outer leaflet, while conversely, the inner leaflets showed minimal curvature and significant lipid tilt around the AHs. The MD-generated profiles at the protein-membrane interface were then extracted and used as boundary conditions in a continuum elastic membrane model to calculate the membrane-bending energy of each conformation embedded in different membrane surfaces characteristic of a budding virus. The calculations show that all three M2 conformations are stabilized in inward-budding, concave spherical caps and destabilized in outward-budding, convex spherical caps, the latter reminiscent of a budding virus. One of the C2-broken symmetry conformations is stabilized by 4 kT in NGC surfaces with the minimum energy conformation occurring at a curvature corresponding to 33 nm radii. In total, our work provides atomistic insight into the curvature sensing capabilities of M2 channels and how enrichment in the nascent viral particle depends on protein shape and membrane geometry.

Activation of the influenza B M2 proton channel (BM2)

Influenza B viruses have co-circulated during most seasonal flu epidemics and can cause significant human morbidity and mortality due to their rapid mutation, emerging drug resistance, and severe impact on vulnerable populations. The influenza B M2 proton channel (BM2) plays an essential role in viral replication, but the mechanisms behind its symmetric proton conductance and the involvement of a second histidine (His27) cluster remain unclear. Here we perform the membrane-enabled continuous constant-pH molecular dynamics simulations on wildtype BM2 and a key H27A mutant to explore its pH-dependent conformational switch. Simulations capture the activation as the first histidine (His19) protonates and reveal the transition at lower pH values compared to AM2 is a result of electrostatic repulsions between His19 and pre-protonated His27. Crucially, we provide an atomic-level understanding of the symmetric proton conduction by identifying pre-activating channel hydration in the C-terminal portion. This research advances our understanding of the function of BM2 function and lays the groundwork for further chemically reactive modeling of the explicit proton transport process as well as possible anti-flu drug design efforts.

Solid state NMR spectral editing of histidine, arginine and lysine using Hadamard encoding

The NMR signals from protein sidechains are rich in information about intra- and inter-molecular interactions, but their detection can be complicated due to spectral overlap as well as conformational and hydrogen exchange. In this work, we demonstrate a protocol for multi-dimensional solid-state NMR spectral editing of signals from basic sidechains based on Hadamard matrix encoding. The Hadamard method acquires multi-dimensional experiments in such a way that both the backbone and under-sampled sidechain signals can be decoded for unambiguous editing in the 15 N spectral frequency dimension. All multi-dimensional 15 N-edited solid-state NMR experiments can be acquired using this strategy, thereby accelerating the acquisition of spectra spanning broad frequency bandwidth. Application of these methods to the ferritin nanocage, reveals signals from N atoms from His, Arg, Lys and Trp sidechains, as well as their tightly bound, ordered water molecules. The Hadamard approach adds to the arsenal of spectroscopic approaches for protein NMR signal detection.

Simulating asymmetric membranes using P21 periodic boundary conditions

Molecular dynamics (MD) simulations of symmetric lipid bilayers are now well established, while those of asymmetric ones are considerably less developed. This disjunction arises in part because the surface tensions of leaflets in asymmetric bilayers can differ (unlike those of symmetric ones), and there is no simple way to determine them without assumptions. This chapter describes the use of P21 periodic boundary conditions (PBC), which allow lipids to switch leaflets, to generate asymmetric bilayers under the assumption of equal chemical potentials of lipids in opposing leaflets. A series of examples, ranging from bilayers with one lipid type to those with peptides and proteins, provides a guide for the use of P21 PBC. Critical properties of asymmetric membranes, such as spontaneous curvature, are highly sensitive to differences in the leaflet surface tensions (or differential stress), and equilibration with P21 PBC substantially reduces differential stress of asymmetric bilayers assembled with surface area-based methods. Limitations of the method are discussed. Technically, the nonstandard unit cell is difficult to parallelize and to incorporate restraints. Inherently, the assumption of equal chemical potentials, and therefore the method itself, is not applicable to all target systems. Despite these limitations, it is argued that P21 simulations should be considered when designing equilibration protocols for MD studies of most asymmetric membranes.

Tackling Influenza A virus by M2 ion channel blockers: Latest progress and limitations

Influenza outbreaks cause pandemics in millions of people. The treatment of influenza remains a challenge due to significant genetic polymorphism in the influenza virus. Also, developing vaccines to protect against seasonal and pandemic influenza infections is constantly impeded. Thus, antibiotics are the only first line of defense against antigenically distinct strains or new subtypes of influenza viruses. Among several anti-influenza targets, the M2 protein of the influenza virus performs several activities. M2 protein is an ion channel that permits proton conductance through the virion envelope and the deacidification of the Golgi apparatus. Both these functions are critical for viral replication. Thus, targeting the M2 protein of the influenza virus is an essential target. Rimantadine and amantadine are two well-known drugs that act on the M2 protein. However, these drugs acquired resistance to influenza and thus are not recommended to treat influenza infections. This review discusses an overview of anti-influenza therapy, M2 ion channel functions, and its working principle. It also discusses the M2 structure and its role, and the change in the structure leads to mutant variants of influenza A virus. We also shed light on the recently identified compounds acting against wild-type and mutated M2 proteins of influenza virus A. These scaffolds could be an alternative to M2 inhibitors and be developed as antibiotics for treating influenza infections.

The Influenza A Virus Replication Cycle: A Comprehensive Review

Influenza A virus (IAV) is the primary causative agent of influenza, colloquially called the flu. Each year, it infects up to a billion people, resulting in hundreds of thousands of human deaths, and causes devastating avian outbreaks with worldwide losses worth billions of dollars. Always present is the possibility that a highly pathogenic novel subtype capable of direct human-to-human transmission will spill over into humans, causing a pandemic as devastating if not more so than the 1918 influenza pandemic. While antiviral drugs for influenza do exist, they target very few aspects of IAV replication and risk becoming obsolete due to antiviral resistance. Antivirals targeting other areas of IAV replication are needed to overcome this resistance and combat the yearly epidemics, which exact a serious toll worldwide. This review aims to summarise the key steps in the IAV replication cycle, along with highlighting areas of research that need more focus.

A comparative study of influenza A M2 protein conformations in DOPC/DOPS liposomes and in native E. coli membranes

We compared the conformations of the transmembrane domain (TMD) of influenza A M2 (IAM2) protein reconstituted at pH 7.4 in DOPC/DOPS bilayers to those in isolated E. coli membranes, having preserved its native proteins and lipids. IAM2 is a single-pass transmembrane protein known to assemble into homo-tetrameric proton channel. To represent this channel, we made a construct containing the IAM2's TMD region flanked by the juxtamembrane residues. The single cysteine substitute, L43C, of leucine located in the bilayer polar region was paramagnetically tagged with a methanethiosulfonate nitroxide label for the ESR (electron spin resonance) study. We compared the conformations of the spin-labeled IAM2 residing in DOPC/DOPS and native E. coli membranes using continuous-wave (CW) ESR and double electron-electron resonance (DEER) spectroscopy. The total protein-to-lipid molar ratio spanned the range from 1:230 to 1:10,400⩦ The CW ESR spectra corresponded to a nearly rigid limit spin label dynamics in both environments. In all cases, the DEER data were reconstructed into the distance distributions showing well-resolved peaks at 1.68 nm and 2.37 nm. The peak distance ratio was 1.41±0.2 and the amplitude ratio was 2:1. This is what one expects from four nitroxide spin-labels located at the corners of a square, indicative of an axially symmetric tetramer. Distance modeling of DEER data with molecular modeling software applied to the NMR molecular structures (PDB: 2L0J) confirmed the symmetry and closed state of the C-terminal exit pore of the IAM2 tetramer in agreement with the NMR model. Thus, we can conclude that IAM2 TMD has similar conformations in model and native E. coli membranes of comparable thickness and fluidity, notwithstanding the complexity of the E. coli membranes caused by their lipid diversity and the abundance of integral and peripheral membrane proteins.

Real-time tracking of drug binding to influenza A M2 reveals a high energy barrier

The drug Rimantadine binds to two different sites in the M2 protein from influenza A, a peripheral site and a pore site that is the primary site of efficacy. It remained enigmatic that pore binding did not occur in certain detergent micelles, and in particular incomplete binding was observed in a mixture of lipids selected to match the viral membrane. Here we show that two effects are responsible, namely changes in the protein upon pore binding that prevented detergent solubilization, and slow binding kinetics in the lipid samples. Using 55-100 kHz magic-angle spinning NMR, we characterize kinetics of drug binding in three different lipid environments: DPhPC, DPhPC with cholesterol and viral mimetic membrane lipid bilayers. Slow pharmacological binding kinetics allowed the characterization of spectral changes associated with non-specific binding to the protein periphery in the kinetically trapped pore-apo state. Resonance assignments were determined from a set of proton-detected 3D spectra. Chemical shift changes associated with functional binding in the pore of M2 were tracked in real time in order to estimate the activation energy. The binding kinetics are affected by pH and the lipid environment and in particular cholesterol. We found that the imidazole-imidazole hydrogen bond at residue histidine 37 is a stable feature of the protein across several lipid compositions. Pore binding breaks the imidazole-imidazole hydrogen bond and limits solubilization in DHPC detergent.

In Situ Spectroscopic Detection of Large-Scale Reorientations of Transmembrane Helices During Influenza A M2 Channel Opening

Viroporins are small ion channels in membranes of enveloped viruses that play key roles during viral life cycles. To use viroporins as drug targets against viral infection requires in-depth mechanistic understanding and, with that, methods that enable investigations under in situ conditions. Here, we apply surface-enhanced infrared absorption (SEIRA) spectroscopy to Influenza A M2 reconstituted within a solid-supported membrane, to shed light on the mechanics of its viroporin function. M2 is a paradigm of pH-activated proton channels and controls the proton flux into the viral interior during viral infection. We use SEIRA to track the large-scale reorientation of M2's transmembrane α-helices in situ during pH-activated channel opening. We quantify this event as a helical tilt from 26° to 40° by correlating the experimental results with solid-state nuclear magnetic resonance-informed computational spectroscopy. This mechanical motion is impeded upon addition of the inhibitor rimantadine, giving a direct spectroscopic marker to test antiviral activity. The presented approach provides a spectroscopic tool to quantify large-scale structural changes and to track the function and inhibition of the growing number of viroporins from pathogenic viruses in future studies.

The role of an amphiphilic helix and transmembrane region in the efficient acylation of the M2 protein from influenza virus

Protein palmitoylation, a cellular process occurring at the membrane-cytosol interface, is orchestrated by members of the DHHC enzyme family and plays a pivotal role in regulating various cellular functions. The M2 protein of the influenza virus, which is acylated at a membrane-near amphiphilic helix serves as a model for studying the intricate signals governing acylation and its interaction with the cognate enzyme, DHHC20. We investigate it here using both experimental and computational assays. We report that altering the biophysical properties of the amphiphilic helix, particularly by shortening or disrupting it, results in a substantial reduction in M2 palmitoylation, but does not entirely abolish the process. Intriguingly, DHHC20 exhibits an augmented affinity for some M2 mutants compared to the wildtype M2. Molecular dynamics simulations unveil interactions between amino acids of the helix and the catalytically significant DHHC and TTXE motifs of DHHC20. Our findings suggest that the binding of M2 to DHHC20, while not highly specific, is mediated by requisite contacts, possibly instigating the transfer of fatty acids. A comprehensive comprehension of protein palmitoylation mechanisms is imperative for the development of DHHC-specific inhibitors, holding promise for the treatment of diverse human diseases.

Dimeric Transmembrane Structure of the SARS-CoV-2 E Protein

The SARS-CoV-2 E protein is a transmembrane (TM) protein with its N-terminus exposed on the external surface of the virus. At debate is its oligomeric state, let alone its function. Here, the TM structure of the E protein is characterized by oriented sample and magic angle spinning solid-state NMR in lipid bilayers and refined by molecular dynamics simulations. This protein was previously found to be a pentamer, with a hydrophobic pore that appears to function as an ion channel. We identify only a front-to-front, symmetric helix-helix interface, leading to a dimeric structure that does not support channel activity. The two helices have a tilt angle of only 6°, resulting in an extended interface dominated by Leu and Val sidechains. While residues Val14-Thr35 are almost all buried in the hydrophobic region of the membrane, Asn15 lines a water-filled pocket that potentially serves as a drug-binding site. The E and other viral proteins may adopt different oligomeric states to help perform multiple functions.

Rapid Determination of the Topology of Oligomeric α-Helical Membrane Proteins by Water- and Lipid-Edited Methyl NMR

Single-span oligomeric α-helical transmembrane proteins are common in virus ion channels, which are targets of antiviral drugs. Knowledge about the high-resolution structures of these oligomeric α-helical bundles is so far scarce. Structure determination of these membrane proteins by solid-state NMR traditionally requires resolving and assigning protein chemical shifts and measuring many interhelical distances, which are time-consuming. To accelerate experimental structure determination, here we introduce a simple solid-state NMR approach that uses magnetization transfer from water and lipid protons to the protein. By detecting the water- and lipid-transferred intensities of the high-sensitivity methyl 13C signals of Leu, Val, and Ile residues, which are highly enriched in these membrane proteins, we can derive models of the topology of these homo-oligomeric helical bundles. The topology is specified by the positions of amino acid residues in heptad repeats and the orientations of residues relative to the channel pore, lipids, and the helical interface. We demonstrate this water- and lipid-edited methyl NMR approach on the envelope (E) protein of SARS-CoV-2, the causative agent of the COVID-19 pandemic. We show that water-edited and lipid-edited 2D 13C-13C correlation spectra can be measured with sufficient sensitivity. Even without resolving multiple residues of the same type in the NMR spectra, we can obtain the helical bundle topology. We apply these experiments to the structurally unknown E proteins of the MERS coronavirus and the human coronavirus NL63. The resulting structural topologies show interesting differences in the positions of the aromatic residues in these three E proteins, suggesting that these viroporins may have different mechanisms of activation and ion conduction.

A Study of the Activity of Adamantyl Amines against Mutant Influenza A M2 Channels Identified a Polycyclic Cage Amine Triple Blocker, Explored by Molecular Dynamics Simulations and Solid-State NMR

We compared the anti-influenza potencies of 57 adamantyl amines and analogs against influenza A virus with serine-31 M2 proton channel, usually termed as WT M2 channel, which is amantadine sensitive. We also tested a subset of these compounds against viruses with the amantadine-resistant L26F, V27A, A30T, G34E M2 mutant channels. Four compounds inhibited WT M2 virus in vitro with mid-nanomolar potency, with 27 compounds showing sub-micromolar to low micromolar potency. Several compounds inhibited L26F M2 virus in vitro with sub-micromolar to low micromolar potency, but only three compounds blocked L26F M2-mediated proton current as determined by electrophysiology (EP). One compound was found to be a triple blocker of WT, L26F, V27A M2 channels by EP assays, but did not inhibit V27A M2 virus in vitro, and one compound inhibited WT, L26F, V27A M2 in vitro without blocking V27A M2 channel. One compound blocked only L26F M2 channel by EP, but did not inhibit virus replication. The triple blocker compound is as long as rimantadine, but could bind and block V27A M2 channel due to its larger girth as revealed by molecular dynamics simulations, while MAS NMR informed on the interaction of the compound with M2(18-60) WT or L26F or V27A.

Dimeric Transmembrane Structure of the SARS-CoV-2 E Protein

The SARS-CoV-2 E protein is a transmembrane (TM) protein with its N-terminus exposed on the external surface of the virus. Here, the TM structure of the E protein is characterized by oriented sample and magic angle spinning solid-state NMR in lipid bilayers and refined by molecular dynamics simulations. This protein has been found to be a pentamer, with a hydrophobic pore that appears to function as an ion channel. We identified only a symmetric helix-helix interface, leading to a dimeric structure that does not support channel activity. The two helices have a tilt angle of only 6°, resulting in an extended interface dominated by Leu and Val sidechains. While residues Val14-Thr35 are almost all buried in the hydrophobic region of the membrane, Asn15 lines a water-filled pocket that potentially serves as a drug-binding site. The E and other viral proteins may adopt different oligomeric states to help perform multiple functions.

Activation of mouse Otop3 proton channels by Zn2+

Otopetrins (Otop1-Otop3) belong to a newly identified family of proton (H⁺) channels activated by extracellular acidification. Here, we found that Zn²⁺ activates the mouse Otop3 (mOtop3) proton channels by using electrophysiological patch-clamp techniques. In mOtop3-expressing human embryonic kidney HEK293T cells, a biphasic inward mOtop3 H⁺ current comprising a fast transient current followed by a sustained current was observed upon extracellular acidification at pH 5.0. No significant activation of the mOtop3 channel was observed at pH 6.5 and 7.4, but interestingly, Zn²⁺ dose-dependently induced a sustained activation of mOtop3 under these pH conditions. Increasing the Zn²⁺ concentration had no effect on the reversal potential of the channel currents, suggesting that Zn²⁺ does not permeate through the mOtop3. The activation of the mOtop3 channel was specific to Zn²⁺ among divalent metal cations. Our findings reveal a novel modulatory mechanism of mOtop3 proton channels by Zn²⁺.

The Cytoplasmic Domain of the SARS-CoV-2 Envelope Protein Assembles into a β-Sheet Bundle in Lipid Bilayers

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) envelope (E) protein forms a pentameric ion channel in the lipid membrane of the endoplasmic reticulum Golgi intermediate compartment (ERGIC) of the infected cell. The cytoplasmic domain of E interacts with host proteins to cause virus pathogenicity and may also mediate virus assembly and budding. To understand the structural basis of these functions, here we investigate the conformation and dynamics of an E protein construct (residues 8-65) that encompasses the transmembrane domain and the majority of the cytoplasmic domain using solid-state NMR. ¹³C and ¹⁵N chemical shifts indicate that the cytoplasmic domain adopts a β-sheet-rich conformation that contains three β-strands separated by turns. The five subunits associate into an umbrella-shaped bundle that is attached to the transmembrane helices by a disordered loop. Water-edited NMR spectra indicate that the third β-strand at the C terminus of the protein is well hydrated, indicating that it is at the surface of the β-bundle. The structure of the cytoplasmic domain cannot be uniquely determined from the inter-residue correlations obtained here due to ambiguities in distinguishing intermolecular and intramolecular contacts for a compact pentameric assembly of this small domain. Instead, we present four structural topologies that are consistent with the measured inter-residue contacts. These data indicate that the cytoplasmic domain of the SARS-CoV-2 E protein has a strong propensity to adopt β-sheet conformations when the protein is present at high concentrations in lipid bilayers. The equilibrium between the β-strand conformation and the previously reported α-helical conformation may underlie the multiple functions of E in the host cell and in the virion.

Molecular origins of asymmetric proton conduction in the influenza M2 channel

The M2 proton channel of influenza A is embedded into the viral envelope and allows acidification of the virion when the external pH is lowered. In contrast, no outward proton conductance is observed when the internal pH is lowered, although outward current is observed at positive voltage. Residues Trp41 and Asp44 are known to play a role in preventing pH-driven outward conductance, but the mechanism for this is unclear. We investigate this issue using classical molecular dynamics simulations with periodic proton hops. When all key His37 residues are neutral, inward proton movement is much more facile than outward movement if the His are allowed to shuttle the proton. The preference for inward movement increases further as the charge on the His37 increases. Analysis of the trajectories reveals three factors accounting for this asymmetry. First, in the outward direction, Asp44 traps the hydronium by strong electrostatic interactions. Secondly, Asp44 and Trp41 orient the hydronium with the protons pointing inward, hampering outward Grotthus hopping. As a result, the effective barrier is lower in the inward direction. Trp41 adds to the barrier by weakly H-bonding to potential H+ acceptors. Finally, for charged His, the H3O+ in the inner vestibule tends to get trapped at lipid-lined fenestrations of the cone-shaped channel. Simulations qualitatively reproduce the experimentally observed higher outward conductance of mutants. The ability of positive voltage, unlike proton gradient, to induce an outward current appears to arise from its ability to bias H3O+ and the waters around it toward more H-outward orientations.

Cholesterol and M2 Rendezvous in Budding and Scission of Influenza A Virus

The cholesterol of the host cell plasma membrane and viral M2 protein plays a crucial role in multiple stages of infection and replication of the influenza A virus. Cholesterol is required for the formation of heterogeneous membrane microdomains (or rafts) in the budozone of the host cell that serves as assembly sites for the viral components. The raft microstructures act as scaffolds for several proteins. Cholesterol may further contribute to the mechanical forces necessary for membrane scission in the last stage of budding and help to maintain the stability of the virus envelope. The M2 protein has been shown to cause membrane scission in model systems by promoting the formation of curved lipid bilayer structures that, in turn, can lead to membrane vesicles budding off or scission intermediates. Membrane remodeling by M2 is intimately linked with cholesterol as it affects local lipid composition, fluidity, and stability of the membrane. Thus, both cholesterol and M2 protein contribute to the efficient and proper release of newly formed influenza viruses from the virus-infected cells.

A proof-of-concept study of the secondary structure of influenza A, B M2 and MERS- and SARS-CoV E transmembrane peptides using folding molecular dynamics simulations in a membrane mimetic solvent

Here, we have carried out a proof-of-concept molecular dynamics (MD) simulation with adaptive tempering in a membrane mimetic environment to study the folding of single-pass membrane peptides. We tested the influenza A M2 viroporin, influenza B M2 viroporin, and protein E from coronaviruses MERS-Cov-2 and SARS-CoV-2 peptides with known experimental secondary structures in membrane bilayers. The two influenza-derived peptides are significantly different in the peptide sequence and secondary structure and more polar than the two coronavirus-derived peptides. Through a total of more than 50 μs of simulation time that could be accomplished in trifluoroethanol (TFE), as a membrane model, we characterized comparatively the folding behavior, helical stability, and helical propensity of these transmembrane peptides that match perfectly their experimental secondary structures, and we identified common motifs that reflect their quaternary organization and known (or not) biochemical function. We showed that BM2 is organized into two structurally distinct parts: a significantly more stable N-terminal half, and a fast-converting C-terminal half that continuously folds and unfolds between α-helical structures and non-canonical structures, which are mostly turns. In AM2, both the N-terminal half and C-terminal half are very flexible. In contrast, the two coronavirus-derived transmembrane peptides are much more stable and fast helix-formers when compared with the influenza ones. In particular, the SARS-derived peptide E appears to be the fastest and most stable helix-former of all the four viral peptides studied, with a helical structure that persists almost without disruption for the whole of its 10 μs simulation. By comparing the results with experimental observations, we benchmarked TFE in studying the conformation of membrane and hydrophobic peptides. This work provided accurate results suggesting a methodology to run long MD simulations and predict structural properties of biologically important membrane peptides.

Divide and Conquer: A Tailored Solid‐state NMR Approach to Study Large Membrane Protein Complexes

Membrane proteins are known to exert many essential biological functions by forming complexes in cell membranes. An example refers to the β‐barrel assembly machinery (BAM), a 200 kDa pentameric complex containing BAM proteins A–E that catalyzes the essential process of protein insertion into the outer membrane of gram‐negative bacteria. While progress has been made in capturing three‐dimensional structural snapshots of the BAM complex, the role of the lipoprotein BamC in the complex assembly in functional lipid bilayers has remained unclear. We have devised a component‐selective preparation scheme to directly study BamC as part of the entire BAM complex in lipid bilayers. Combination with proton‐detected solid‐state NMR methods allowed us to probe the structure, dynamics, and supramolecular topology of full‐length BamC embedded in the entire complex in lipid bilayers. Our approach may help decipher how individual proteins contribute to the dynamic formation and functioning of membrane protein complexes in membranes.

Breathing and tilting: mesoscale simulations illuminate influenza glycoprotein vulnerabilities

Influenza virus has resurfaced recently from inactivity during the early stages of the COVID-19 pandemic, raising serious concerns about the nature and magnitude of future epidemics. The main antigenic targets of influenza virus are two surface glycoproteins, hemagglutinin (HA) and neuraminidase (NA). Whereas the structural and dynamical properties of both glycoproteins have been studied previously, the understanding of their plasticity in the whole-virion context is fragmented. Here, we investigate the dynamics of influenza glycoproteins in a crowded protein environment through mesoscale all-atom molecular dynamics simulations of two evolutionary-linked glycosylated influenza A whole-virion models. Our simulations reveal and kinetically characterize three main molecular motions of influenza glycoproteins: NA head tilting, HA ectodomain tilting, and HA head breathing. The flexibility of HA and NA highlights antigenically relevant conformational states, as well as facilitates the characterization of a novel monoclonal antibody, derived from human convalescent plasma, that binds to the underside of the NA head. Our work provides previously unappreciated views on the dynamics of HA and NA, advancing the understanding of their interplay and suggesting possible strategies for the design of future vaccines and antivirals against influenza.

Clustering of tetrameric influenza M2 peptides in lipid bilayers investigated by 19F solid-state NMR

The influenza M2 protein forms a drug-targeted tetrameric proton channel to mediate virus uncoating, and carries out membrane scission to enable virus release. While the proton channel function of M2 has been extensively studied, the mechanism by which M2 catalyzes membrane scission is still not well understood. Previous fluorescence and electron microscopy studies indicated that M2 tetramers concentrate at the neck of the budding virus in the host plasma membrane. However, molecular evidence for this clustering is scarce. Here, we use ¹⁹F solid-state NMR to investigate M2 clustering in phospholipid bilayers. By mixing equimolar amounts of 4F-Phe47 labeled M2 peptide and CF₃-Phe47 labeled M2 peptide and measuring F-CF₃ cross peaks in 2D ¹⁹F¹⁹F correlation spectra, we show that M2 tetramers form nanometer-scale clusters in lipid bilayers. This clustering is stronger in cholesterol-containing membranes and phosphatidylethanolamine (PE) membranes than in cholesterol-free phosphatidylcholine and phosphatidylglycerol membranes. The observed correlation peaks indicate that Phe47 sidechains from different tetramers are less than ~2 nm apart. ¹H¹⁹F correlation peaks between lipid chain protons and fluorinated Phe47 indicate that Phe47 is more deeply inserted into the lipid bilayer in the presence of cholesterol than in its absence, suggesting that Phe47 preferentially interacts with cholesterol. Static ³¹P NMR spectra indicate that M2 induces negative Gaussian curvature in the PE membrane. These results suggest that M2 tetramers cluster at cholesterol- and PE-rich regions of cell membranes to cause membrane curvature, which in turn can facilitate membrane scission in the last step of virus budding and release.

Gating Mechanism of the Voltage-Gated Proton Channel Studied by Molecular Dynamics Simulations

The voltage-gated proton channel Hv1 has important roles in proton extrusion, pH homeostasis, sperm motility, and cancer progression. The Hv1 channel has also been found to be highly expressed in cell lines and tissue samples from patients with breast cancer. A high-resolution closed-state structure has been reported for the mouse Hv1 chimera channel (mHv1cc), solved by X-ray crystallography, but the open-state structure of Hv1 has not been solved. Since Hv1 is a promising drug target, various groups have proposed open conformations by molecular modeling and simulation studies. However, the gating mechanism and the open-state conformation under the membrane potential are still debate. Here, we present a molecular dynamics study considering membrane potential and pH conditions. The closed-state structure of mHv1cc was used to run molecular dynamics (MD) simulations with respect to electric field and pH conditions in order to investigate the mechanism of proton transfer. We observed a continuous hydrogen bond chain of water molecules called a water-wire to be formed through the channel pore in the channel opening, triggered by downward displacement of the S2 helix and upward movement of the S4 helix relative to other helices. Due to the movement of the S2 and S4 helices, the internal salt bridge network was rearranged, and the hydrophobic gating layers were destroyed. In line with previous experimental and simulation observations, our simulation results led us to propose a new gating mechanism for the Hv1 proton channel, and may provide valuable information for novel drug discovery.

Rotational Dynamics of The Transmembrane Domains Play an Important Role in Peptide Dynamics of Viral Fusion and Ion Channel Forming Proteins—A Molecular Dynamics Simulation Study

Focusing on the transmembrane domains (TMDs) of viral fusion and channel-forming proteins (VCPs), experimentally available and newly generated peptides in an ideal conformation of the S and E proteins of severe acute respiratory syndrome coronavirus type 2 (SARS-CoV-2) and SARS-CoV, gp41 and Vpu, both of human immunodeficiency virus type 1 (HIV-1), haemagglutinin and M2 of influenza A, as well as gB of herpes simplex virus (HSV), are embedded in a fully hydrated lipid bilayer and used in multi-nanosecond molecular dynamics simulations. It is aimed to identify differences in the dynamics of the individual TMDs of the two types of viral membrane proteins. The assumption is made that the dynamics of the individual TMDs are decoupled from their extra-membrane domains, and that the mechanics of the TMDs are distinct from each other due to the different mechanism of function of the two types of proteins. The diffusivity coefficient (DC) of the translational and rotational diffusion is decreased in the oligomeric state of the TMDs compared to those values when calculated from simulations in their monomeric state. When comparing the calculations for two different lengths of the TMD, a longer full peptide and a shorter purely TMD stretch, (i) the difference of the calculated DCs begins to level out when the difference exceeds approximately 15 amino acids per peptide chain, and (ii) the channel protein rotational DC is the most affected diffusion parameter. The rotational dynamics of the individual amino acids within the middle section of the TMDs of the fusion peptides remain high upon oligomerization, but decrease for the channel peptides, with an increasing number of monomers forming the oligomeric state, suggesting an entropic penalty on oligomerization for the latter.

Structural and functional analysis of the roles of influenza C virus membrane proteins in assembly and budding

Assembly and budding of the influenza C virus is mediated by three membrane proteins: the hemagglutinin-esterase-fusion glycoprotein (HEF), the matrix protein (CM1), and the ion channel (CM2). Here we investigated whether the formation of the hexagonal HEF arrangement, a distinctive feature of influenza C virions is important for virus budding. We used super resolution microscopy and found 250-nm sized HEF clusters at the plasma membrane of transfected cells, which were insensitive to cholesterol extraction and cytochalasin treatment. Overexpression of either CM1, CM2, or HEF caused the release of membrane-enveloped particles. Cryo-electron microscopy of the latter revealed spherical vesicles exhibiting the hexagonal HEF clusters. We subsequently used reverse genetics to identify elements in HEF required for this clustering. We found that deletion of the short cytoplasmic tail of HEF reduced virus titer and hexagonal HEF arrays, suggesting that an interaction with CM1 stabilizes the HEF clusters. In addition, we substituted amino acids at the surface of the closed HEF conformation and identified specific mutations that prevented virus rescue, others reduced virus titers and the number of HEF clusters in virions. Finally, mutation of two regions that mediate contacts between trimers in the in-situ structure of HEF was shown to prevent rescue of infectious virus particles. Mutations at residues thought to mediate lateral interactions were revealed to promote intracellular trafficking defects. Taken together, we propose that lateral interactions between the ectodomains of HEF trimers are a driving force for virus budding, although CM2 and CM1 also play important roles in this process.

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.

Emulating Membrane Protein Environments─How Much Lipid Is Required for a Native Structure: Influenza S31N M2

This report investigates the homotetrameric membrane protein structure of the S31N M2 protein from Influenza A virus in the presence of a high molar ratio of lipid. The structured regions of this protein include a single transmembrane helix and an amphipathic helix. Two structures of the S31N M2 conductance domain from Influenza A virus have been deposited in the Protein Data Bank (PDB). These structures present different symmetries about the channel main axis. We present new magic angle spinning and oriented sample solid-state NMR spectroscopic data for S31N M2 in liquid crystalline lipid bilayers using protein tetramer:lipid molar ratios ranging from 1:120 to 1:240. The data is consistent with an essentially 4-fold-symmetric structure very similar to the M2 WT structure that also has a single conformation for the four monomers, except at the His37 and Trp41 functional sites when characterized in samples with a high molar ratio of lipid. While detergent solubilization is well recognized today as a nonideal environment for small membrane proteins, here we discuss the influence of a high lipid to protein ratio for samples of the S31N M2 protein to stabilize an essentially 4-fold-symmetric conformation of the M2 membrane protein. While it is generally accepted that the chemical and physical properties of the native environment of membrane proteins needs to be reproduced judiciously to achieve the native protein structure, here we show that not only the character of the emulated membrane environment is important but also the abundance of the environment is important for achieving the native structure. This is a critical finding as a membrane protein spectroscopist's goal is always to generate a sample with the highest possible protein sensitivity while obtaining spectra of the native-like structure.

Multiscale Simulation of an Influenza A M2 Channel Mutant Reveals Key Features of Its Markedly Different Proton Transport Behavior

The influenza A M2 channel, a prototype for viroporins, is an acid-activated viroporin that conducts protons across the viral membrane, a critical step in the viral life cycle. Four central His37 residues control channel activation by binding subsequent protons from the viral exterior, which opens the Trp41 gate and allows proton flux to the interior. Asp44 is essential for maintaining the Trp41 gate in a closed state at high pH, resulting in asymmetric conduction. The prevalent D44N mutant disrupts this gate and opens the C-terminal end of the channel, resulting in increased conduction and a loss of this asymmetric conduction. Here, we use extensive Multiscale Reactive Molecular Dynamics (MS-RMD) and quantum mechanics/molecular mechanics (QM/MM) molecular dynamics simulations with an explicit, reactive excess proton to calculate the free energy of proton transport in this M2 mutant and to study the dynamic molecular-level behavior of D44N M2. We find that this mutation significantly lowers the barrier of His37 deprotonation in the activated state and shifts the barrier for entry to the Val27 tetrad. These free energy changes are reflected in structural shifts. Additionally, we show that the increased hydration around the His37 tetrad diminishes the effect of the His37 charge on the channel's water structure, facilitating proton transport and enabling activation from the viral interior. Altogether, this work provides key insight into the fundamental characteristics of PT in WT M2 and how the D44N mutation alters this PT mechanism, and it expands understanding of the role of emergent mutations in viroporins.

From Acid Activation Mechanisms of Proton Conduction to Design of Inhibitors of the M2 Proton Channel of Influenza A Virus

With an estimated 1 billion people affected across the globe, influenza is one of the most serious health concerns worldwide. Therapeutic treatments have encompassed a number of key functional viral proteins, mainly focused on the M2 proton channel and neuraminidase. This review highlights the efforts spent in targeting the M2 proton channel, which mediates the proton transport toward the interior of the viral particle as a preliminary step leading to the release of the fusion peptide in hemagglutinin and the fusion of the viral and endosomal membranes. Besides the structural and mechanistic aspects of the M2 proton channel, attention is paid to the challenges posed by the development of efficient small molecule inhibitors and the evolution toward novel ligands and scaffolds motivated by the emergence of resistant strains.

Backbone Torsion Angle Determination Using Proton Detected Magic-Angle Spinning Nuclear Magnetic Resonance

Protein torsion angles define the backbone secondary structure of proteins. Magic-angle spinning (MAS) NMR methods using carbon detection have been developed to measure torsion angles by determining the relative orientation between two anisotropic interactions─dipolar coupling or chemical shift anisotropy. Here we report a new proton-detection based method to determine the backbone torsion angle by recoupling NH and CH dipolar couplings within the HCANH pulse sequence, for protonated or partly deuterated samples. We demonstrate the efficiency and precision of the method with microcrystalline chicken α spectrin SH3 protein and the influenza A matrix 2 (M2) membrane protein, using 55 or 90 kHz MAS. For M2, pseudo-4D data detect a turn between transmembrane and amphipathic helices.

Spiers Memorial Lecture: Analysis and de novo design of membrane-interactive peptides

Membrane-peptide interactions play critical roles in many cellular and organismic functions, including protection from infection, remodeling of membranes, signaling, and ion transport. Peptides interact with membranes in a variety of ways: some associate with membrane surfaces in either intrinsically disordered conformations or well-defined secondary structures. Peptides with sufficient hydrophobicity can also insert vertically as transmembrane monomers, and many associate further into membrane-spanning helical bundles. Indeed, some peptides progress through each of these stages in the process of forming oligomeric bundles. In each case, the structure of the peptide and the membrane represent a delicate balance between peptide-membrane and peptide-peptide interactions. We will review this literature from the perspective of several biologically important systems, including antimicrobial peptides and their mimics, α-synuclein, receptor tyrosine kinases, and ion channels. We also discuss the use of de novo design to construct models to test our understanding of the underlying principles and to provide useful leads for pharmaceutical intervention of diseases.

Inside and Out of the Pore: Comparing Interactions and Molecular Dynamics of Influenza A M2 Viroporin Complexes in Standard Lipid Bilayers

Ion channels located at viral envelopes (viroporins) have a critical function for the replication of infectious viruses and are important drug targets. Over the last decade, the number and duration of molecular dynamics (MD) simulations of the influenza A M2 ion channel owing to the increased computational efficiency. Here, we aimed to define the system setup and simulation conditions for the correct description of the protein-pore and the protein-lipid interactions for influenza A M2 in comparison with experimental data. We performed numerous MD simulations of the influenza A M2 protein in complex with adamantane blockers in standard lipid bilayers using OPLS2005 and CHARMM36 (C36) force fields. We explored the effect of varying the M2 construct (M2(22-46) and M2(22-62)), the lipid buffer size and type (stiffer DMPC or softer POPC with or without 20% cholesterol), the simulation time, the H37 protonation site (N_δ or Ν_ε), the conformational state of the W41 channel gate, and M2's cholesterol binding sites (BSs). We report that the 200 ns MD with M2(22-62) (having N_ε Η37) in the 20 Å lipid buffer with the C36 force field accurately describe: (a) the M2 pore structure and interactions inside the pore, that is, adamantane channel blocker location, water clathrate structure, and water or chloride anion blockage/passage from the M2 pore in the presence of a channel blocker and (b) interactions between M2 and the membrane environment as reflected by the calculation of the M2 bundle tilt, folding of amphipathic helices, and cholesterol BSs. Strikingly, we also observed that the C36 1 μs MD simulations using M2(22-62) embedded in a 20 Å POPC:cholesterol (5:1) scrambled membrane produced frequent interactions with cholesterol, which when combined with computational kinetic analysis, revealed the experimentally observed BSs of cholesterol and suggested three similarly long-interacting positions in the top leaflet that have previously not been observed experimentally. These findings promise to be useful for other viroporin systems.

Pore-Bound Water at the Key Residue Histidine 37 in Influenza A M2

Atomic details of structured water molecules are indispensable to understand the thermodynamics of important biological processes including the proton conduction mechanism of the M2 protein. Despite the expectation of structured water molecules based on crystal structures of Influenza A M2, only two water populations have been observed by NMR in reconstituted lipid bilayer samples. These are the bulk‐ and lipid‐associated water populations typically seen in membrane samples. Here, we detect a bound water molecule at a chemical shift of 11 ppm, located near the functional histidine 37 residue in the M2 conductance domain, which comprises residues 18 to 60. Combining 100 kHz magic‐angle spinning NMR, dynamic nuclear polarization and density functional theory calculations, we show that the bound water forms a hydrogen bond to the δ1 nitrogen of histidine 37.

Virus Structures and Dynamics by Magic-Angle Spinning NMR

Techniques for atomic-resolution structural biology have evolved during the past several decades. Breakthroughs in instrumentation, sample preparation, and data analysis that occurred in the past decade have enabled characterization of viruses with an unprecedented level of detail. Here we review the recent advances in magic-angle spinning (MAS) nuclear magnetic resonance (NMR) spectroscopy for structural analysis of viruses and viral assemblies. MAS NMR is a powerful method that yields information on 3D structures and dynamics in a broad range of experimental conditions. After a brief introduction, we discuss recent structural and functional studies of several viruses investigated with atomic resolution at various levels of structural organization, from individual domains of a membrane protein reconstituted into lipid bilayers to virus-like particles and intact viruses. We present examples of the unique information revealed by MAS NMR about drug binding, conduction mechanisms, interactions with cellular host factors, and DNA packaging in biologically relevant environments that are inaccessible by other methods.

Membrane Materials for Selective Ion Separations at the Water-Energy Nexus

Synthetic polymer membranes are enabling components in key technologies at the water-energy nexus, including desalination and energy conversion, because of their high water/salt selectivity or ionic conductivity. However, many applications at the water-energy nexus require ion selectivity, or separation of specific ionic species from other similar species. Here, the ion selectivity of conventional polymeric membrane materials is assessed and recent progress in enhancing selective transport via tailored free volume elements and ion-membrane interactions is described. In view of the limitations of polymeric membranes, three material classes-porous crystalline materials, 2D materials, and discrete biomimetic channels-are highlighted as possible candidates for ion-selective membranes owing to their molecular-level control over physical and chemical properties. Lastly, research directions and critical challenges for developing bioinspired membranes with molecular recognition are provided.

Effects of Cholesterol on the Partitioning of a Drug Molecule in Lipid Bilayers

Drug molecules either bind to membrane-bound targets or permeate through cell membranes to reach intracellular targets, and hence, their membrane partition and permeation are of great importance. Here, we studied the effects of cholesterol on the partition of amantadine, an antiflu drug molecule, into 1,2-dimyristoyl-sn-glycero-3-phosphocholine bilayers using molecular dynamics simulations. The membrane partition of amantadine is sensitive to the cholesterol mole fraction (x_chol). In the absence of cholesterol, amantadine is stably bound in membranes, but at x_chol = 32%, it can escape to the aqueous phase, in agreement with recent experiments. The reduced membrane partition of amantadine at a high cholesterol content is mainly due to the perturbation of the bilayer structure and dynamics. Surrounding lipids stabilize amantadine by having their tails wrapped around the drug molecule, and this ability is compromised when cholesterol is present to increase the order in lipid tails. The atomic details on interactions with lipids and perturbations by cholesterol revealed here provide insight into membrane partition and delivery of drug molecules to their targets.

Solid-State NMR of Membrane Proteins in Lipid Bilayers: To Spin or Not To Spin?

Membrane proteins mediate a plethora of cellular functions and represent important targets for drug development. Unlike soluble proteins, membrane proteins require native-like environments to fold correctly and be active. Therefore, modern structural biology techniques have aimed to determine the structure and dynamics of these membrane proteins at physiological temperature and in liquid crystalline lipid bilayers. With the flourishing of new NMR methodologies and improvements in sample preparations, magic angle spinning (MAS) and oriented sample solid-state NMR (OS-ssNMR) spectroscopy of membrane proteins is experiencing a new renaissance. Born as antagonistic approaches, these techniques nowadays offer complementary information on the structural topology and dynamics of membrane proteins reconstituted in lipid membranes. By spinning biosolid samples at the magic angle (θ = 54.7°), MAS NMR experiments remove the intrinsic anisotropy of the NMR interactions, increasing spectral resolution. Internuclear spin interactions (spin exchange) are reintroduced by RF pulses, providing distances and torsion angles to determine secondary, tertiary, and quaternary structures of membrane proteins. OS-ssNMR, on the other hand, directly detects anisotropic NMR parameters such as dipolar couplings (DC) and anisotropic chemical shifts (CS), providing orientational constraints to determine the architecture (i.e., topology) of membrane proteins relative to the lipid membrane. Defining the orientation of membrane proteins and their interactions with lipid membranes is of paramount importance since lipid-protein interactions can shape membrane protein conformations and ultimately define their functional states.In this Account, we report selected studies from our group integrating MAS and OS-ssNMR techniques to give a comprehensive view of the biological processes occurring at cellular membranes. We focus on the main experiments for both techniques, with an emphasis on new implementation to increase both sensitivity and spectral resolution. We also describe how the structural constraints derived from both isotropic and anisotropic NMR parameters are integrated into dynamic structural modeling using replica-averaged orientational-restrained molecular dynamics simulations (RAOR-MD). We showcase small membrane proteins that are involved in Ca²⁺ transport and regulate cardiac and skeletal muscle contractility: phospholamban (PLN, 6 kDa), sarcolipin (SLN, 4 kDa), and DWORF (4 kDa). We summarize our results for the structures of these polypeptides free and in complex with the sarcoplasmic reticulum (SR) Ca²⁺-ATPase (SERCA, 110 kDa). Additionally, we illustrate the progress toward the determination of the structural topology of a six transmembrane protein associated with succinate and acetate transport (SatP, hexamer 120 kDa). From these examples, the integrated MAS and OS-ssNMR approach, in combination with modern computational methods, emerges as a way to overcome the challenges posed by studying large membrane protein systems.

Influenza Viruses: Harnessing the Crucial Role of the M2 Ion-Channel and Neuraminidase toward Inhibitor Design

As a member of the Orthomyxoviridae family of viruses, influenza viruses (IVs) are known causative agents of respiratory infection in vertebrates. They remain a major global threat responsible for the most virulent diseases and global pandemics in humans. The virulence of IVs and the consequential high morbidity and mortality of IV infections are primarily attributed to the high mutation rates in the IVs' genome coupled with the numerous genomic segments, which give rise to antiviral resistant and vaccine evading strains. Current therapeutic options include vaccines and small molecule inhibitors, which therapeutically target various catalytic processes in IVs. However, the periodic emergence of new IV strains necessitates the continuous development of novel anti-influenza therapeutic options. The crux of this review highlights the recent studies on the biology of influenza viruses, focusing on the structure, function, and mechanism of action of the M2 channel and neuraminidase as therapeutic targets. We further provide an update on the development of new M2 channel and neuraminidase inhibitors as an alternative to existing anti-influenza therapy. We conclude by highlighting therapeutic strategies that could be explored further towards the design of novel anti-influenza inhibitors with the ability to inhibit resistant strains.

1H/13C/15N triple-resonance experiments for structure determinaton of membrane proteins by oriented-sample NMR

The benefits of triple-resonance experiments for structure determination of macroscopically oriented membrane proteins by solid-state NMR are discussed. While double-resonance ¹H/¹⁵N experiments are effective for structure elucidation of alpha-helical domains, extension of the method of oriented samples to more complex topologies and assessing side-chain conformations necessitates further development of triple-resonance (¹H/¹³C/¹⁵N) NMR pulse sequences. Incorporating additional spectroscopic dimensions involving ¹³C spin-bearing nuclei, however, introduces essential complications arising from the wide frequency range of the ¹H-¹³C dipolar couplings and ¹³C CSA (>20 ​kHz), and the presence of the ¹³C-¹³C homonuclear dipole-dipole interactions. The recently reported ROULETTE-CAHA pulse sequence, in combination with the selective z-filtering, can be used to evolve the structurally informative ¹H-¹³C dipolar coupling arising from the aliphatic carbons while suppressing the signals from the carbonyl and methyl regions. Proton-mediated magnetization transfer under mismatched Hartman-Hahn conditions (MMHH) can be used to correlate ¹³C and ¹⁵N nuclei in such triple-resonance experiments for the subsequent ¹⁵N detection. The recently developed pulse sequences are illustrated for n-acetyl Leucine (NAL) single crystal and doubly labeled Pf1 coat protein reconstituted in magnetically aligned bicelles. An interesting observation is that in the case of ¹⁵N-labeled NAL measured at ¹³C natural abundance, the triple (¹H/¹³C/¹⁵N) MMHH scheme predominantly gives rise to long-range intermolecular magnetization transfers from ¹³C to ¹⁵N spins; whereas direct Hartmann-Hahn ¹³C/¹⁵N transfer is entirely intramolecular. The presented developments advance NMR of oriented samples for structure determination of membrane proteins and liquid crystals.

Native Chemical Ligation of Highly Hydrophobic Peptides in Ionic Liquid-Containing Media

The chemical synthesis of a highly hydrophobic membrane-associated peptide by native chemical ligation (NCL) in an ionic liquid (IL) [C₂mim][OAc]/buffer mixture was achieved by employing peptide concentrations up to 11 mM. NCL was studied at different pH and water content and compared to several “gold-standard” ligation protocols. The optimized reaction protocol for the NCL in IL required the addition of 40% water and pH adjustment to 7.0–7.5, resulting in ligation yields of up to 80–95% within 1 to 4 h. This new ligation protocol is generally applicable and outperforms current “gold-standard” NCL methods.

Solid-state NMR spectroscopy

Solid-state nuclear magnetic resonance (NMR) spectroscopy is an atomic-level method used to determine the chemical structure, three-dimensional structure, and dynamics of solids and semi-solids. This Primer summarizes the basic principles of NMR as applied to the wide range of solid systems. The fundamental nuclear spin interactions and the effects of magnetic fields and radiofrequency pulses on nuclear spins are the same as in liquid-state NMR. However, because of the anisotropy of the interactions in the solid state, the majority of high-resolution solid-state NMR spectra is measured under magic-angle spinning (MAS), which has profound effects on the types of radiofrequency pulse sequences required to extract structural and dynamical information. We describe the most common MAS NMR experiments and data analysis approaches for investigating biological macromolecules, organic materials, and inorganic solids. Continuing development of sensitivity-enhancement approaches, including 1H-detected fast MAS experiments, dynamic nuclear polarization, and experiments tailored to ultrahigh magnetic fields, is described. We highlight recent applications of solid-state NMR to biological and materials chemistry. The Primer ends with a discussion of current limitations of NMR to study solids, and points to future avenues of development to further enhance the capabilities of this sophisticated spectroscopy for new applications.

NMR-Based Methods for Protein Analysis

Nuclear magnetic resonance (NMR) spectroscopy is a well-established method for analyzing protein structure, interaction, and dynamics at atomic resolution and in various sample states including solution state, solid state, and membranous environment. Thanks to rapid NMR methodology development, the past decade has witnessed a growing number of protein NMR studies in complex systems ranging from membrane mimetics to living cells, which pushes the research frontier further toward physiological environments and offers unique insights in elucidating protein functional mechanisms. In particular, in-cell NMR has become a method of choice for bridging the huge gap between structural biology and cell biology. Herein, we review the recent developments and applications of NMR methods for protein analysis in close-to-physiological environments, with special emphasis on in-cell protein structural determination and the analysis of protein dynamics, both difficult to be accessed by traditional methods.

Highlighting membrane protein structure and function: A celebration of the Protein Data Bank

Biological membranes define the boundaries of cells and compartmentalize the chemical and physical processes required for life. Many biological processes are carried out by proteins embedded in or associated with such membranes. Determination of membrane protein (MP) structures at atomic or near-atomic resolution plays a vital role in elucidating their structural and functional impact in biology. This endeavor has determined 1198 unique MP structures as of early 2021. The value of these structures is expanded greatly by deposition of their three-dimensional (3D) coordinates into the Protein Data Bank (PDB) after the first atomic MP structure was elucidated in 1985. Since then, free access to MP structures facilitates broader and deeper understanding of MPs, which provides crucial new insights into their biological functions. Here we highlight the structural and functional biology of representative MPs and landmarks in the evolution of new technologies, with insights into key developments influenced by the PDB in magnifying their impact.

Anti-Influenza Drug Discovery and Development: Targeting the Virus and Its Host by All Possible Means

Infections by influenza virus constitute a major and recurrent threat for human health. Together with vaccines, antiviral drugs play a key role in the prevention and treatment of influenza virus infection and disease. Today, the number of antiviral molecules approved for the treatment of influenza is relatively limited, and their use is threatened by the emergence of viral strains with resistance mutations. There is therefore a real need to expand the prophylactic and therapeutic arsenal. This chapter summarizes the state of the art in drug discovery and development for the treatment of influenza virus infections, with a focus on both virus-targeting and host cell-targeting strategies. Novel antiviral strategies targeting other viral proteins or targeting the host cell, some of which are based on drug repurposing, may be used in combination to strengthen our therapeutic arsenal against this major pathogen.

Mechanism and Kinetics of Copper Complexes Binding to the Influenza A M2 S31N and S31N/G34E Channels

Copper(II) is known to bind in the influenza virus His37 cluster in the homotetrameric M2 proton channel and block the proton current needed for uncoating. Copper complexes based on iminodiacetate also block the M2 proton channel and show reduced cytotoxicity and zebrafish-embryo toxicity. In voltage-clamp oocyte studies using the ubiquitous amantadine-insensitive M2 S31N variant, the current block showed fast and slow phases, in contrast to the single phase found for amantadine block of wild-type M2. Here, we evaluate the mechanism of block by copper adamantyl iminodiacitate and copper cyclooctyl iminodiacitate complexes and address whether the complexes can coordinate with one or more of the His37 imidazoles. The current traces were fitted to parametrized master equations. The energetics of binding and the rate constants suggest that the first step is copper complex binding within the channel, and the slow step in the current block is the formation of a Cu-histidine coordination complex. Solution-phase isothermal titration calorimetry and density functional theory (DFT) calculations indicate that imidazole binds to the copper complexes. Structural optimization using DFT reveals that the complexes fit inside the channel and project the Cu(II) toward the His37 cluster, allowing one imidazole to form a coordination complex with Cu(II). Electrophysiology and DFT studies also show that the complexes block the G34E amantadine-resistant mutant despite some crowding in the binding site by the glutamates.