Acid activation mechanism of the influenza A M2 proton channel

Conformational heterogeneity and structural features for function of the prototype viroporin influenza AM2

The 97-residue influenza A matrix 2 (ΑM2) protein, a prototype for viroporins, transports protons through water molecules and His37. We discuss structural biology and molecular biophysics experiments and some functional assays that have transformed over 40 years our understanding of the structure and function of AM2. The structural studies on ΑM2 have been performed with different conditions (pH, temperature, lipid, constructs) and using various protein constructs, e.g., AM2 transmembrane (AM2TM) domain, AM2 conductance domain (AM2CD), ectodomain-containing or ectodomain-truncated, AM2 full length (AM2FL) and aimed to describe the different conformations and structural details that are necessary for the stability and function of AM2. However, the conclusions from these experiments appeared sometimes ambiguous and caused exciting debates. This was not due to inaccurate measurements, but instead because of the different membrane mimetic environment used, e.g., detergent, micelles or phospholipid bilayer, the method (e.g., X-ray crystallography, solid state NMR, solution NMR, native mass spectrometry), the used protein construct (e.g., AM2TM or AM2CD), or the amino acids residues to follow observables (e.g., NMR chemical shifts). We present these results according to the different used biophysical methods, the research groups and often by keeping a chronological order for presenting the progress in the research. We discuss ideas for additional research on structural details of AM2 and how the present findings can be useful to explore new routes of influenza A inhibition. The AM2 research can provide inspiration to study other viroporins as drug targets.

Quantitative insights into the mechanism of proton conduction and selectivity for the human voltage-gated proton channel Hv1

Human voltage-gated proton (hHv1) channels are crucial for regulating essential biological processes such as immune cell respiratory burst, sperm capacitation, and cancer cell migration. Despite the significant concentration difference between protons and other ions in physiological conditions, hHv1 demonstrates remarkable proton selectivity. Our calculations of single-proton, cation, and anion permeation free energy profiles quantitatively demonstrate that the proton selectivity of the wild-type channel originates from its strong proton affinity via the titration of the key residues D112 and D174, although the channel imposes similar kinetic blocking effects for protons compared to other ions. A two-proton knock-on model is proposed to mathematically explain the electrophysiological measurements of the pH-dependent proton conductance in the conductive state. Moreover, it is shown that the anion selectivity of the D112N mutant channel is tied to impaired proton transport and substantial anion leakage.

Structure and dynamics of the proton-selective histidine and the gating tryptophan in an inward rectifying hybrid influenza B and A virus M2 proton channel

The M2 proteins of influenza A and B viruses form acid-activated proton channels that are essential for the virus lifecycle. Proton selectivity is achieved by a transmembrane (TM) histidine whereas gating is achieved by a tryptophan residue. Although this functional apparatus is conserved between AM2 and BM2 channels, AM2 conducts protons exclusively inward whereas BM2 conducts protons in either direction depending on the pH gradient. Previous studies showed that in AM2, mutations of D44 abolished inward rectification of AM2, suggesting that the tryptophan gate is destabilized. To elucidate how charged residues C-terminal to the tryptophan regulates channel gating, here we investigate the structure and dynamics of H19 and W23 in a BM2 mutant, GDR-BM2, in which three BM2 residues are mutated to the corresponding AM2 residues, S16G, G26D and H27R. Whole-cell electrophysiological data show that GDR-BM2 conducts protons with inward rectification, identical to wild-type (WT) AM2 but different from WT-BM2. Solid-state NMR 15N and 13C spectra of H19 indicate that the mutant BM2 channel contains higher populations of cationic histidine and neutral τ tautomers compared to WT-BM2 at acidic pH. Moreover, 19F NMR spectra of 5-19F-labeled W23 resolve three peaks at acidic pH, suggesting three tryptophan sidechain conformations. Comparison of these spectra with the tryptophan spectra of other M2 peptides suggests that these indole sidechain conformations arise from interactions with the C-terminal charged residues and with the N-terminal cationic histidine. Taken together, these solid-state NMR data show that inward rectification in M2 proton channels is accomplished by tryptophan interactions with charged residues on both its C-terminal and N-terminal sides. Gating of these M2 proton channels is thus accomplished by a multi-residue complex with finely tuned electrostatic and aromatic interactions.

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.

Comparing Multifunctional Viral and Eukaryotic Proteins for Generating Scission Necks in Membranes

Deterministic formation of membrane scission necks by protein machinery with multiplexed functions is critical in biology. A microbial example is M2 viroporin, a proton pump from the influenza A virus that is multiplexed with membrane remodeling activity to induce budding and scission in the host membrane during viral maturation. In comparison, the dynamin family constitutes a class of eukaryotic proteins implicated in mitochondrial fission, as well as various budding and endocytosis pathways. In the case of Dnm1, the mitochondrial fission protein in yeast, the membrane remodeling activity is multiplexed with mechanoenzyme activity to create fission necks. It is not clear why these functions are combined in these scission processes, which occur in drastically different compositions and solution conditions. In general, direct experimental access to changing neck sizes induced by individual proteins or peptide fragments is challenging due to the nanoscale dimensions and influence of thermal fluctuations. Here, we use a mechanical model to estimate the size of scission necks by leveraging small-angle X-ray scattering structural data of protein-lipid systems under different conditions. The influence of interfacial tension, lipid composition, and membrane budding morphology on the size of the induced scission necks is systematically investigated using our data and molecular dynamic simulations. We find that the M2 budding protein from the influenza A virus has robust pH-dependent membrane activity that induces nanoscopic necks within the range of spontaneous hemifission for a broad range of lipid compositions. In contrast, the sizes of scission necks generated by mitochondrial fission proteins strongly depend on lipid composition, which suggests a role for mechanical constriction.

Molecular Dynamics Simulation of Complex Reactivity with the Rapid Approach for Proton Transport and Other Reactions (RAPTOR) Software Package

Simulating chemically reactive phenomena such as proton transport on nanosecond to microsecond and beyond time scales is a challenging task. Ab initio methods are unable to currently access these time scales routinely, and traditional molecular dynamics methods feature fixed bonding arrangements that cannot account for changes in the system's bonding topology. The Multiscale Reactive Molecular Dynamics (MS-RMD) method, as implemented in the Rapid Approach for Proton Transport and Other Reactions (RAPTOR) software package for the LAMMPS molecular dynamics code, offers a method to routinely sample longer time scale reactive simulation data with statistical precision. RAPTOR may also be interfaced with enhanced sampling methods to drive simulations toward the analysis of reactive rare events, and a number of collective variables (CVs) have been developed to facilitate this. Key advances to this methodology, including GPU acceleration efforts and novel CVs to model water wire formation are reviewed, along with recent applications of the method which demonstrate its versatility and robustness.

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.

Transcendent Aspects of Proton Channels

A handful of biological proton-selective ion channels exist. Some open at positive or negative membrane potentials, others open at low or high pH, and some are light activated. This review focuses on common features that result from the unique properties of protons. Proton conduction through water or proteins differs qualitatively from that of all other ions. Extraordinary proton selectivity is needed to ensure that protons permeate and other ions do not. Proton selectivity arises from a proton pathway comprising a hydrogen-bonded chain that typically includes at least one titratable amino acid side chain. The enormously diverse functions of proton channels in disparate regions of the phylogenetic tree can be summarized by considering the chemical and electrical consequences of proton flux across membranes. This review discusses examples of cells in which proton efflux serves to increase pHi, decrease pHo, control the membrane potential, generate action potentials, or compensate transmembrane movement of electrical charge.

Comparing multifunctional viral and eukaryotic proteins for generating scission necks in membranes

Deterministic formation of membrane scission necks by protein machinery with multiplexed functions is critical in biology. A microbial example is the M2 viroporin, a proton pump from the influenza A virus which is multiplexed with membrane remodeling activity to induce budding and scission in the host membrane during viral maturation. In comparison, the dynamin family constitutes a class of eukaryotic proteins implicated in mitochondrial fission, as well as various budding and endocytosis pathways. In the case of Dnm1, the mitochondrial fission protein in yeast, the membrane remodeling activity is multiplexed with mechanoenzyme activity to create fission necks. It is not clear why these functions are combined in these scission processes, which occur in drastically different compositions and solution conditions. In general, direct experimental access to changing neck sizes induced by individual proteins or peptide fragments is challenging due to the nanoscale dimensions and influence of thermal fluctuations. Here, we use a mechanical model to estimate the size of scission necks by leveraging Small-Angle X-ray Scattering (SAXS) structural data of protein-lipid systems under different conditions. The influence of interfacial tension, lipid composition, and membrane budding morphology on the size of the induced scission necks is systematically investigated using our data and molecular dynamic simulations. We find that the M2 budding protein from the influenza A virus has robust pH-dependent membrane activity that induces nanoscopic necks within the range of spontaneous hemi-fission for a broad range of lipid compositions. In contrast, the sizes of scission necks generated by mitochondrial fission proteins strongly depend on lipid composition, which suggests a role for mechanical constriction.

Computational Characterization of the Reactivity of Compound I in Unspecific Peroxygenases

Unspecific peroxygenases (UPOs) are emerging as promising biocatalysts for selective oxyfunctionalization of unactivated C-H bonds. However, their potential in large-scale synthesis is currently constrained by suboptimal chemical selectivity. Improving the selectivity of UPOs requires a deep understanding of the molecular basis of their catalysis. Recent molecular simulations have sought to unravel UPO's selectivity and inform their design principles. However, most of these studies focused on substrate-binding poses. Few researchers have investigated how the reactivity of CpdI, the principal oxidizing intermediate in the catalytic cycle, influences selectivity in a realistic protein environment. Moreover, the influence of protein electrostatics on the reaction kinetics of CpdI has also been largely overlooked. To bridge this gap, we used multiscale simulations to interpret the regio- and enantioselective hydroxylation of the n-heptane substrate catalyzed by Agrocybe aegerita UPO (AaeUPO). We comprehensively characterized the energetics and kinetics of the hydrogen atom-transfer (HAT) step, initiated by CpdI, and the subsequent oxygen rebound step forming the product. Notably, our approach involved both free energy and potential energy evaluations in a quantum mechanics/molecular mechanics (QM/MM) setting, mitigating the dependence of results on the choice of initial conditions. These calculations illuminate the thermodynamics and kinetics of the HAT and oxygen rebound steps. Our findings highlight that both the conformational selection and the distinct chemical reactivity of different substrate hydrogen atoms together dictate the regio- and enantio-selectivity. Building on our previous study of CpdI's formation in AaeUPO, our results indicate that the HAT step is the rate-limiting step in the overall catalytic cycle. The subsequent oxygen rebound step is swift and retains the selectivity determined by the HAT step. We also pinpointed several polar and charged amino acid residues whose electrostatic potentials considerably influence the reaction barrier of the HAT step. Notably, the Glu196 residue is pivotal for both the CpdI's formation and participation in the HAT step. Our research offers in-depth insights into the catalytic cycle of AaeUPO, which will be instrumental in the rational design of UPOs with enhanced properties.

Understanding the Multifaceted Mechanism of Compound I Formation in Unspecific Peroxygenases through Multiscale Simulations

Unspecific peroxygenases (UPOs) can selectively oxyfunctionalize unactivated hydrocarbons by using peroxides under mild conditions. They circumvent the oxygen dilemma faced by cytochrome P450s and exhibit greater stability than the latter. As such, they hold great potential for industrial applications. A thorough understanding of their catalysis is needed to improve their catalytic performance. However, it remains elusive how UPOs effectively convert peroxide to Compound I (CpdI), the principal oxidizing intermediate in the catalytic cycle. Previous computational studies of this process primarily focused on heme peroxidases and P450s, which have significant differences in the active site from UPOs. Additionally, the roles of peroxide unbinding in the kinetics of CpdI formation, which is essential for interpreting existing experiments, have been understudied. Moreover, there has been a lack of free energy characterizations with explicit sampling of protein and hydration dynamics, which is critical for understanding the thermodynamics of the proton transport (PT) events involved in CpdI formation. To bridge these gaps, we employed multiscale simulations to comprehensively characterize the CpdI formation in wild-type UPO from Agrocybe aegerita (AaeUPO). Extensive free energy and potential energy calculations were performed in a quantum mechanics/molecular mechanics setting. Our results indicate that substrate-binding dehydrates the active site, impeding the PT from H2O2 to a nearby catalytic base (Glu196). Furthermore, the PT is coupled with considerable hydrogen bond network rearrangements near the active site, facilitating subsequent O-O bond cleavage. Finally, large unbinding free energy barriers kinetically stabilize H2O2 at the active site. These findings reveal a delicate balance among PT, hydration dynamics, hydrogen bond rearrangement, and cosubstrate unbinding, which collectively enable efficient CpdI formation. Our simulation results are consistent with kinetic measurements and offer new insights into the CpdI formation mechanism at atomic-level details, which can potentially aid the design of next-generation biocatalysts for sustainable chemical transformations of feedstocks.

The M2 proteins of bat influenza A viruses reveal atypical features compared to conventional M2 proteins

The influenza A virus (IAV) M2 protein has proton channel activity, which plays a role in virus uncoating and may help to preserve the metastable conformation of the IAV hemagglutinin (HA). In contrast to the highly conserved M2 proteins of conventional IAV, the primary sequences of bat IAV H17N10 and H18N11 M2 proteins show remarkable divergence, suggesting that these proteins may differ in their biological function. We, therefore, assessed the proton channel activity of bat IAV M2 proteins and investigated its role in virus replication. Here, we show that the M2 proteins of bat IAV did not fully protect acid-sensitive HA of classical IAV from low pH-induced conformational change, indicating low proton channel activity. Interestingly, the N31S substitution not only rendered bat IAV M2 proteins sensitive to inhibition by amantadine but also preserved the metastable conformation of acid-sensitive HA to a greater extent. In contrast, the acid-stable HA of H18N11 did not rely on such support by M2 protein. When mutant M2(N31S) protein was expressed in the context of chimeric H18N11/H5N1(6:2) encoding HA and NA of avian IAV H5N1, amantadine significantly inhibited virus entry, suggesting that ion channel activity supported virus uncoating. Finally, the cytoplasmic domain of the H18N11 M2 protein mediated rapid internalization of the protein from the plasma membrane leading to low-level expression at the cell surface. However, cell surface levels of H18N11 M2 protein were significantly enhanced in cells infected with the chimeric H18N11/H5N1(6:2) virus. The potential role of the N1 sialidase in arresting M2 internalization is discussed. IMPORTANCE Bat IAV M2 proteins not only differ from the homologous proteins of classical IAV by their divergent primary sequence but are also unable to preserve the metastable conformation of acid-sensitive HA, indicating low proton channel activity. This unusual feature may help to avoid M2-mediated cytotoxic effects and inflammation in bats infected with H17N10 or H18N11. Unlike classical M2 proteins, bat IAV M2 proteins with the N31S substitution mediated increased protection of HA from acid-induced conformational change. This remarkable gain of function may help to understand how single point mutations can modulate proton channel activity. In addition, the cytoplasmic domain was found to be responsible for the low cell surface expression level of bat IAV M2 proteins. Given that the M2 cytoplasmic domain of conventional IAV is well known to participate in virus assembly at the plasma membrane, this atypical feature might have consequences for bat IAV budding and egress.

Antiviral Approaches against Influenza Virus

Preventing and controlling influenza virus infection remains a global public health challenge, as it causes seasonal epidemics to unexpected pandemics. These infections are responsible for high morbidity, mortality, and substantial economic impact. Vaccines are the prophylaxis mainstay in the fight against influenza. However, vaccination fails to confer complete protection due to inadequate vaccination coverages, vaccine shortages, and mismatches with circulating strains. Antivirals represent an important prophylactic and therapeutic measure to reduce influenza-associated morbidity and mortality, particularly in high-risk populations. Here, we review current FDA-approved influenza antivirals with their mechanisms of action, and different viral- and host-directed influenza antiviral approaches, including immunomodulatory interventions in clinical development. Furthermore, we also illustrate the potential utility of machine learning in developing next-generation antivirals against influenza.

Tracking the Delocalized Proton in Concerted Proton Transfer in Bulk Water

A solvated proton in water is often characterized as a charge or structural defect, and it is important to track its evolution on-the-fly in certain dynamics simulations. Previously, we introduced the proton indicator, a pseudo-atom, whose position approximates the location of the excess proton modeled as a structural defect. The proton indicator generally yields a smooth trajectory of a hydrated proton diffusing in aqueous solutions, including in the events of stepwise proton transfer via the Grotthuss mechanism; however, the proton indicator did not perform well in the events of concerted proton transfer, for which it occasionally yielded large position displacements between two successive time steps. To overcome this hurdle, we develop a new algorithm of a proton indicator with greatly enhanced performance for concerted proton transfer in bulk water. A protocol is proposed to exhaustively explore the hydrogen-bonding network of the water wires over which the excess proton is delocalized and to properly account for the contributions of the water molecules in this network as the geometry evolves. The new proton indicator (called Indicator 2.0) is assessed in dynamics simulations of an excess proton in bulk water and in specially constructed model systems of more complex architectures. The results demonstrate that the new indicator yields a smooth trajectory in both stepwise and concerted proton transfers.

Influenza A Virus: Cellular Entry

The frequent emergence of pathogenic viruses with pandemic potential has posed a significant threat to human health and economy, despite enormous advances in our understanding of infection mechanisms and devising countermeasures through developing various prophylactic and therapeutic strategies. The recent coronavirus disease (COVID-19) pandemic has re-emphasised the importance of rigorous research on virus infection mechanisms and highlighted the need for our preparedness for potential pandemics. Although viruses cannot self-replicate, they tap into host cell factors and processes for their entry, propagation and dissemination. Upon entering the host cells, viruses ingeniously utilise the innate biological functions of the host cell to replicate themselves and maintain their existence in the hosts. Influenza A virus (IAV), which has a negative-sense, single-stranded RNA as its genome, is no exception. IAVs are enveloped viruses with a lipid bilayer derived from the host cell membrane and have a surface covered with the spike glycoprotein haemagglutinin (HA) and neuraminidase (NA). Viral genome is surrounded by an M1 shell, forming a "capsid" in the virus particle. IAV particles use HA to recognise sialic acids on the cell surface of lung epithelial cells for their attachment. After attachment to the cell surface, IAV particles are endocytosed and sorted into the early endosomes. Subsequently, as the early endosomes mature into late endosomes, the endosomal lumen becomes acidified, and the low pH of the late endosomes induces conformational reaggangements in the HA to initiate fusion between the endosomal and viral membranes. Upon fusion, the viral capsid disintegrates and the viral ribonucleoprotein (vRNP) complexes containing the viral genome are released into the cytosol. The process of viral capsid disintegration is called "uncoating". After successful uncoating, the vRNPs are imported into the nucleus by importin α/β (IMP α/β), where viral replication and transcription take place and the new vRNPs are assembled. Recently, we have biochemically elucidated the molecular mechanisms of the processes of viral capsid uncoating subsequent viral genome dissociation. In this chapter, we present the molecular details of the viral uncoating process.

Generalized Transition State Theory Treatment of Water-Assisted Proton Transport Processes in Proteins

Transition state theory (TST) is widely employed for estimating the transition rate of a reaction when combined with free energy sampling techniques. A derivation of the transition theory rate expression for a general n-dimensional case is presented in this work which specifically focuses on water-assisted proton transfer/transport reactions, especially for protein systems. Our work evaluates the TST prefactor calculated at the transition state dividing surface compared to one sampled, as an approximation, in the reactant state in four case studies of water-assisted proton transport inside membrane proteins and highlights the significant impact of the prefactor position dependence in proton transport processes.

Pre-Training of Equivariant Graph Matching Networks with Conformation Flexibility for Drug Binding

The latest biological findings observe that the motionless "lock-and-key" theory is not generally applicable and that changes in atomic sites and binding pose can provide important information for understanding drug binding. However, the computational expenditure limits the growth of protein trajectory-related studies, thus hindering the possibility of supervised learning. A spatial-temporal pre-training method based on the modified equivariant graph matching networks, dubbed ProtMD which has two specially designed self-supervised learning tasks: atom-level prompt-based denoising generative task and conformation-level snapshot ordering task to seize the flexibility information inside molecular dynamics (MD) trajectories with very fine temporal resolutions is presented. The ProtMD can grant the encoder network the capacity to capture the time-dependent geometric mobility of conformations along MD trajectories. Two downstream tasks are chosen to verify the effectiveness of ProtMD through linear detection and task-specific fine-tuning. A huge improvement from current state-of-the-art methods, with a decrease of 4.3% in root mean square error for the binding affinity problem and an average increase of 13.8% in the area under receiver operating characteristic curve and the area under the precision-recall curve for the ligand efficacy problem is observed. The results demonstrate a strong correlation between the magnitude of conformation's motion in the 3D space and the strength with which the ligand binds with its receptor.

Biomimetic Artificial Proton Channels

One of the most common biochemical processes is the proton transfer through the cell membranes, having significant physiological functions in living organisms. The proton translocation mechanism has been extensively studied; however, mechanistic details of this transport are still needed. During the last decades, the field of artificial proton channels has been in continuous growth, and understanding the phenomena of how confined water and channel components mediate proton dynamics is very important. Thus, proton transfer continues to be an active area of experimental and theoretical investigations, and acquiring insights into the proton transfer mechanism is important as this enlightenment will provide direct applications in several fields. In this review, we present an overview of the development of various artificial proton channels, focusing mostly on their design, self-assembly behavior, proton transport activity performed on bilayer membranes, and comparison with protein proton channels. In the end, we discuss their potential applications as well as future development and perspectives.

Amino acid deprotonation rates from classical force fields

Acid ionization constants (pKa's) of titratable amino acid side chains have received a large amount of experimental and theoretical attention. In many situations, however, the rates of protonation and deprotonation, kon & koff, may also be important, for example in understanding the mechanism of action of proton channels or membrane proteins that couple proton transport to other processes. Protonation and deprotonation involve making and breaking of covalent bonds, which cannot be studied by classical force fields. However, environment effects on the rates should be captured by such methods. Here we present an approach for estimating deprotonation rates based on Warshel's extension of Marcus's theory of electron transfer, with input from molecular simulations. The missing bond dissociation energy is represented by a constant term determined by fitting the pKa value in solution. The statistics of the energy gap between protonated and deprotonated states is used to compute free energy curves of the two states and thus free energy barriers, from which the rate can be deduced. The method is applied to Glu, Asp, and His in bulk solution and select membrane proteins: the M2 proton channel, bacteriorhodopsin, and cytochrome c oxidase.

Multiscale simulation unravels the light-regulated reversible inhibition of dihydrofolate reductase by phototrexate

Molecular photoswitches are widely used in photopharmacology, where the biomolecular functions are photo-controlled reversibly with high spatiotemporal precision. Despite the success of this field, it remains elusive how the protein environment modulates the photochemical properties of photoswitches. Understanding this fundamental question is critical for designing more effective light-regulated drugs with mitigated side effects. In our recent work, we employed first-principles non-adiabatic dynamics simulations to probe the effects of protein on the trans to cis photoisomerization of phototrexate (PTX), a photochromic analog of the anticancer therapeutic methotrexate that inhibits the target enzyme dihydrofolate reductase (DHFR). Building upon this study, in this work, we employ multiscale simulations to unravel the full photocycle underlying the light-regulated reversible inhibition of DHFR by PTX, which remains elusive until now. First-principles non-adiabatic dynamics simulations reveal that the cis to trans photoisomerization quantum yield is hindered in the protein due to backward isomerization on the ground-state following non-adiabatic transition, which arises from the favorable binding of the cis isomer with the protein. However, free energy simulations indicate that cis to trans photoisomerization significantly decreases the binding affinity of the PTX. Thus, the cis to trans photoisomerization most likely precedes the ligand unbinding from the protein. We propose the most probable photocycle of the PTX-DHFR system. Our comprehensive simulations highlight the trade-offs among the binding affinity, photoisomerization quantum yield, and the thermal stability of the ligand's different isomeric forms. As such, our work reveals new design principles of light-regulated drugs in photopharmacology.

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.

Structural motifs for subtype-specific pH-sensitive gating of vertebrate otopetrin proton channels

Otopetrin (OTOP) channels are proton-selective ion channels conserved among vertebrates and invertebrates, with no structural similarity to other ion channels. There are three vertebrate OTOP channels (OTOP1, OTOP2, and OTOP3), of which one (OTOP1) functions as a sour taste receptor. Whether extracellular protons gate OTOP channels, in addition to permeating them, was not known. Here, we compare the functional properties of the three murine OTOP channels using patch-clamp recording and cytosolic pH microfluorimetry. We find that OTOP1 and OTOP3 are both steeply activated by extracellular protons, with thresholds of pH_o <6.0 and 5.5, respectively, and kinetics that are pH-dependent. In contrast, OTOP2 channels are broadly active over a large pH range (pH 5 pH 10) and carry outward currents in response to extracellular alkalinization (>pH 9.0). Strikingly, we could change the pH-sensitive gating of OTOP2 and OTOP3 channels by swapping extracellular linkers that connect transmembrane domains. Swaps of extracellular linkers in the N domain, comprising transmembrane domains 1-6, tended to change the relative conductance at alkaline pH of chimeric channels, while swaps within the C domain, containing transmembrane domains 7-12, tended to change the rates of OTOP3 current activation. We conclude that members of the OTOP channel family are proton-gated (acid-sensitive) proton channels and that the gating apparatus is distributed across multiple extracellular regions within both the N and C domains of the channels. In addition to the taste system, OTOP channels are expressed in the vertebrate vestibular and digestive systems. The distinct gating properties we describe may allow them to subserve varying cell-type specific functions in these and other biological systems.

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.

Molecular interactions of the M and E integral membrane proteins of SARS-CoV-2

Specific lipid-protein interactions are key for cellular processes, and even more so for the replication of pathogens. The COVID-19 pandemic has drastically changed our lives and caused the death of nearly four million people worldwide, as of this writing. SARS-CoV-2 is the virus that causes the disease and has been at the center of scientific research over the past year. Most of the research on the virus is focused on key players during its initial attack and entry into the cellular host; namely the S protein, its glycan shield, and its interactions with the ACE2 receptors of human cells. As cases continue to rise around the globe, and new mutants are identified, there is an urgent need to understand the mechanisms of this virus during different stages of its life cycle. Here, we consider two integral membrane proteins of SARS-CoV-2 known to be important for viral assembly and infectivity. We have used microsecond-long all-atom molecular dynamics to examine the lipid-protein and protein-protein interactions of the membrane (M) and envelope (E) structural proteins of SARS-CoV-2 in a complex membrane model. We contrast the two proposed protein complexes for each of these proteins, and quantify their effect on their local lipid environment. This ongoing work also aims to provide molecular-level understanding of the mechanisms of action of this virus to possibly aid in the design of novel treatments.

Studying the Effect of Amino Acid Substitutions in the M2 Ion Channel of the Influenza Virus on the Antiviral Activity of the Aminoadamantane Derivative In Vitro and In Silico

Purpose: The aminoadamantane derivative of L-histidyl-1-adamantayl ethylamine hydrochloride (HCl*H-His-Rim) has showed a high inhibition level against influenza A virus strains in vitro. The aim of this work is to search and establish evidence of the direct effect of the drug on influenza A virus proton channel M2.

Methods: The compound HCl*H-His-Rim was obtained by classical peptide synthesis methods. Influenza A virus mutants of A/PuertoRico/8/34(H1N1) strain were obtained by reverse genetics methods. The mutant samples of the virus were cultured on chicken embryos with a virus titer in the hemagglutination test. ELISA was carried out on Madin-Darby canine kidney (MDCK) monolayer cells when multiplying the virus 10^-4-10^-6. The binding stability of HCl*H-His-Rim was compared to those of M2 (S31N) and M2 (S31N_A30T) channels by molecular dynamic (MD) modeling. The calculation was performed taking into account the interaction with the model lipid bilayer (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) in the presence of water molecules in accordance with the three-center model.

Results: It was found that HCl*H-His-Rim is a direct action drug against influenza A. The most likely conformation of drug binding to target protein has been shown. It has been found that the A30T mutation reduces the binding energy of the drug, and the results obtained in vitro have confirmed the data calculated in silico.

Conclusion: The mechanism of action of HCl*H-His-Rim is directly related to the suppression of the function of the proton channel M2 of influenza A virus.

A quantitative paradigm for water-assisted proton transport through proteins and other confined spaces

Water-assisted proton transport through confined spaces influences many phenomena in biomolecular and nanomaterial systems. In such cases, the water molecules that fluctuate in the confined pathways provide the environment and the medium for the hydrated excess proton migration via Grotthuss shuttling. However, a definitive collective variable (CV) that accurately couples the hydration and the connectivity of the proton wire with the proton translocation has remained elusive. To address this important challenge-and thus to define a quantitative paradigm for facile proton transport in confined spaces-a CV is derived in this work from graph theory, which is verified to accurately describe water wire formation and breakage coupled to the proton translocation in carbon nanotubes and the Cl^-/H⁺ antiporter protein, ClC-ec1. Significant alterations in the conformations and thermodynamics of water wires are uncovered after introducing an excess proton into them. Large barriers in the proton translocation free-energy profiles are found when water wires are defined to be disconnected according to the new CV, even though the pertinent confined space is still reasonably well hydrated and-by the simple measure of the mere existence of a water structure-the proton transport would have been predicted to be facile via that oversimplified measure. In this paradigm, however, the simple presence of water is not sufficient for inferring proton translocation, since an excess proton itself is able to drive hydration, and additionally, the water molecules themselves must be adequately connected to facilitate any successful proton transport.

Influenza AM2 Channel Oligomerization Is Sensitive to Its Chemical Environment

Viroporins are small viral ion channels that play important roles in the viral infection cycle and are proven antiviral drug targets. Matrix protein 2 from influenza A (AM2) is the best-characterized viroporin, and the current paradigm is that AM2 forms monodisperse tetramers. Here, we used native mass spectrometry and other techniques to characterize the oligomeric state of both the full-length and transmembrane (TM) domain of AM2 in a variety of different pH and detergent conditions. Unexpectedly, we discovered that AM2 formed a range of different oligomeric complexes that were strongly influenced by the local chemical environment. Native mass spectrometry of AM2 in nanodiscs with different lipids showed that lipids also affected the oligomeric states of AM2. Finally, nanodiscs uniquely enabled the measurement of amantadine binding stoichiometries to AM2 in the intact lipid bilayer. These unexpected results reveal that AM2 can form a wider range of oligomeric states than previously thought possible, which may provide new potential mechanisms of influenza pathology and pharmacology.

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.

Improved Indicator Algorithms for Tracking a Hydrated Proton as A Local Structural Defect in Grotthuss Diffusion in Aqueous Solutions

Keeping track of a hydrated proton in dynamics simulations is important and nontrivial. Here, we report two revised algorithms for the proton indicator, a pseudo-atom whose position approximates the location of an excess proton diffusing via the Grotthuss mechanism in aqueous solution. The new methods describe the delocalized proton as a structural defect. Encouragingly, in test simulations of a hydrated proton in bulk water, the new algorithms substantially outperform the original scheme by significantly reducing large displacements in the indicator positions upon donor switch, yielding smoother trajectories that effectively track the movement of the solvated proton.

Effect of cholesterol on the membrane partitioning dynamics of hepatitis A virus-2B peptide

Understanding viral peptide detection and partitioning and the subsequent host membrane composition-based response is essential for gaining insights into the viral mechanism. Here, we probe the crucial role of the presence of membrane lipid packing defects, depending on the membrane composition, in allowing the viral peptide belonging to C-terminal Hepatitis A Virus-2B (HAV-2B) to detect, attach and subsequently partition into host cell membrane mimics. Using molecular dynamics simulations, we conclusively show that the hydrophobic residues in the viral peptide detect transiently present lipid packing defects, insert themselves into such defects, form anchor points and facilitate the partitioning of the peptide, thereby inducing membrane disruption. We also show that the presence of cholesterol significantly alters such lipid packing defects, both in size and in number, thus mitigating the partitioning of the membrane active viral peptide into cholesterol-rich membranes. Our results are in excellent agreement with previously published experimental data and further explain the role of lipid defects in understanding such data. These results show differential ways in which the presence and absence of cholesterol can alter the permeability of the host membranes to the membrane active peptide component of HAV-2B virus, via lipid packing defects, and can possibly be a part of the general membrane detection mechanism for viroporins.

Enhancing acid tolerance of Escherichia coli via viroporin-mediated export of protons and its application for efficient whole-cell biotransformation

Escherichia coli-based whole-cell biocatalysts are widely used for the sustainable production of value-added chemicals. However, weak acids present as substrates and/or products obstruct the growth and fermentation capability of E. coli. Here, we show that a viroporin consisting of the influenza A matrix-2 (M2) protein, is activated by low pH and has proton channel activity in E. coli. The heterologous expression of the M2 protein in E. coli resulted in a significant increase in the intracellular pH and cell viability in the presence of various weak acids with different lengths of carbon chains. In addition, the feasibility of developing a robust and efficient E. coli-based whole-cell biocatalyst via introduction of the proton-selective viroporin was explored by employing (Z)-11-(heptanolyoxy)undec-9-enoic acid (ester) and 2-fucosyllactose (2'-FL) as model products, whose production is hampered by cytosolic acidification. The engineered E. coli strains containing the proton-selective viroporin exhibited approximately 80% and 230% higher concentrations of the ester and 2'-FL, respectively, than the control strains without the M2 protein. The simple and powerful strategy developed in this study can be applied to produce other valuable chemicals whose production involves substrates and/or products that cause cytosolic acidification.

How Influenza Virus Uses Host Cell Pathways during Uncoating

Influenza is a zoonotic respiratory disease of major public health interest due to its pandemic potential, and a threat to animals and the human population. The influenza A virus genome consists of eight single-stranded RNA segments sequestered within a protein capsid and a lipid bilayer envelope. During host cell entry, cellular cues contribute to viral conformational changes that promote critical events such as fusion with late endosomes, capsid uncoating and viral genome release into the cytosol. In this focused review, we concisely describe the virus infection cycle and highlight the recent findings of host cell pathways and cytosolic proteins that assist influenza uncoating during host cell entry.

Molecular interactions of the M and E integral membrane proteins of SARS-CoV-2

Specific lipid-protein interactions are key for cellular processes, and even more so for the replication of pathogens. The COVID-19 pandemic has drastically changed our lives and cause the death of nearly three million people worldwide, as of this writing. SARS-CoV-2 is the virus that causes the disease and has been at the center of scientific research over the past year. Most of the research on the virus is focused on key players during its initial attack and entry into the cellular host; namely the S protein, its glycan shield, and its interactions with the ACE2 receptors of human cells. As cases continue to raise around the globe, and new mutants are identified, there is an urgent need to understand the mechanisms of this virus during different stages of its life cycle. Here, we consider two integral membrane proteins of SARS-CoV-2 known to be important for viral assembly and infectivity. We have used microsecond-long all-atom molecular dynamics to examine the lipid-protein and protein-protein interactions of the membrane (M) and envelope (E) structural proteins of SARS-CoV-2 in a complex membrane model. We contrast the two proposed protein complexes for each of these proteins, and quantify their effect on their local lipid environment. This ongoing work also aims to provide molecular-level understanding of the mechanisms of action of this virus to possibly aid in the design of novel treatments.

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.

Membrane Curvature Revisited-the Archetype of Rhodopsin Studied by Time-Resolved Electronic Spectroscopy

G-protein-coupled receptors (GPCRs) comprise the largest and most pharmacologically targeted membrane protein family. Here, we used the visual receptor rhodopsin as an archetype for understanding membrane lipid influences on conformational changes involved in GPCR activation. Visual rhodopsin was recombined with lipids varying in their degree of acyl chain unsaturation and polar headgroup size using 1-palmitoyl-2-oleoyl-sn-glycero- and 1,2-dioleoyl-sn-glycerophospholipids with phosphocholine (PC) or phosphoethanolamine (PE) substituents. The receptor activation profile after light excitation was measured using time-resolved ultraviolet-visible spectroscopy. We discovered that more saturated POPC lipids back shifted the equilibrium to the inactive state, whereas the small-headgroup, highly unsaturated DOPE lipids favored the active state. Increasing unsaturation and decreasing headgroup size have similar effects that combine to yield control of rhodopsin activation, and necessitate factors beyond proteolipid solvation energy and bilayer surface electrostatics. Hence, we consider a balance of curvature free energy with hydrophobic matching and demonstrate how our data support a flexible surface model (FSM) for the coupling between proteins and lipids. The FSM is based on the Helfrich formulation of membrane bending energy as we previously first applied to lipid-protein interactions. Membrane elasticity and curvature strain are induced by lateral pressure imbalances between the constituent lipids and drive key physiological processes at the membrane level. Spontaneous negative monolayer curvature toward water is mediated by unsaturated, small-headgroup lipids and couples directly to GPCR activation upon light absorption by rhodopsin. For the first time to our knowledge, we demonstrate this modulation in both the equilibrium and pre-equilibrium evolving states using a time-resolved approach.

Influenza A M2 Inhibitor Binding Understood through Mechanisms of Excess Proton Stabilization and Channel Dynamics

Prevalent resistance to inhibitors that target the influenza A M2 proton channel has necessitated a continued drug design effort, supported by a sustained study of the mechanism of channel function and inhibition. Recent high-resolution X-ray crystal structures present the first opportunity to see how the adamantyl amine class of inhibitors bind to M2 and disrupt and interact with the channel’s water network, providing insight into the critical properties that enable their effective inhibition in wild-type M2. In this work, we examine the hypothesis that these drugs act primarily as mechanism-based inhibitors by comparing hydrated excess proton stabilization during proton transport in M2 with the interactions revealed in the crystal structures, using the Multiscale Reactive Molecular Dynamics (MS-RMD) methodology. MS-RMD, unlike classical molecular dynamics, models the hydrated proton (hydronium-like cation) as a dynamic excess charge defect and allows bonds to break and form, capturing the intricate interactions between the hydrated excess proton, protein atoms, and water. Through this, we show that the ammonium group of the inhibitors is effectively positioned to take advantage of the channel’s natural ability to stabilize an excess protonic charge and act as a hydronium mimic. Additionally, we show that the channel is especially stable in the drug binding region, highlighting the importance of this property for binding the adamantane group. Finally, we characterize an additional hinge point near Val27, which dynamically responds to charge and inhibitor binding. Altogether, this work further illuminates a dynamic understanding of the mechanism of drug inhibition in M2, grounded in the fundamental properties that enable the channel to transport and stabilize excess protons, with critical implications for future drug design efforts.

Multiscale Simulation Reveals Passive Proton Transport Through SERCA on the Microsecond Timescale

The sarcoplasmic reticulum Ca²⁺-ATPase (SERCA) transports two Ca²⁺ ions from the cytoplasm to the reticulum lumen at the expense of ATP hydrolysis. In addition to transporting Ca²⁺, SERCA facilitates bidirectional proton transport across the sarcoplasmic reticulum to maintain the charge balance of the transport sites and to balance the charge deficit generated by the exchange of Ca²⁺. Previous studies have shown the existence of a transient water-filled pore in SERCA that connects the Ca²⁺ binding sites with the lumen, but the capacity of this pathway to sustain passive proton transport has remained unknown. In this study, we used the multiscale reactive molecular dynamics method and free energy sampling to quantify the free energy profile and timescale of the proton transport across this pathway while also explicitly accounting for the dynamically coupled hydration changes of the pore. We find that proton transport from the central binding site to the lumen has a microsecond timescale, revealing a novel passive cytoplasm-to-lumen proton flow beside the well-known inverse proton countertransport occurring in active Ca²⁺ transport. We propose that this proton transport mechanism is operational and serves as a functional conduit for passive proton transport across the sarcoplasmic reticulum.

Cholesterol Alters the Orientation and Activity of the Influenza Virus M2 Amphipathic Helix in the Membrane

The influenza virus M2 amphipathic helix (M2AH) alters membrane curvature in a cholesterol-dependent manner, mediating viral membrane scission during influenza virus budding. Here, we have investigated the biophysical effects of cholesterol on the ability of an M2AH peptide to manipulate membrane properties. We see that the ability of the M2AH to interact with membranes and form an α-helix is independent of membrane cholesterol concentration; however, cholesterol affects the angle of the M2AH peptide within the membrane. This change in membrane orientation affects the ability of the M2AH to alter lipid order. In low-cholesterol membranes, the M2AH is inserted near the level of the lipid head groups, increasing lipid order, which may contribute to generation of the membrane curvature. As the cholesterol content increases, the M2AH insertion becomes flatter and slightly deeper in the membrane below the lipid headgroups, where the polar face can continue to interact with the headgroups while the hydrophobic face binds cholesterol. This changed orientation minimizes lipid packing defects and lipid order changes, likely reducing the generation of membrane curvature. Thus, cholesterol regulates M2 membrane scission by precisely modulating M2AH positioning within the membrane. This has implications for the understanding of many of amphipathic-helix-driven cellular budding processes that occur in specific lipid environments.

In Silico Study of Membrane Lipid Composition Regulating Conformation and Hydration of Influenza Virus B M2 Channel

The proton conduction of transmembrane influenza virus B M2 (BM2) proton channel is possibly mediated by the membrane environment, but the detailed molecular mechanism is challenging to determine. In this work, how membrane lipid composition regulates the conformation and hydration of BM2 channel is elucidated in silico. The appearance of several important hydrogen-bond networks has been discovered, as the addition of negatively charged lipid palmitoyloleoyl phosphatidylglycerol (POPG) and cholesterol reduces membrane fluidity and augments membrane rigidity. A more rigid membrane environment is beneficial to expand the channel, allow more water to enter the channel, promote channel hydration, and then even affect the proton conduction facilitated by the hydrated channel. Thus, membrane environment could be identified as an important influence factor of conformation and hydration of BM2. These findings can provide a unique perspective for understanding the mechanism of membrane lipid composition regulating conformation and hydration of BM2 and have important significance to the further study of anti-influenza virus B drugs.

Understanding and Tracking the Excess Proton in Ab Initio Simulations; Insights from IR Spectra

Proton transport in aqueous media is ubiquitously important in chemical and biological processes. Although ab initio molecular dynamics (AIMD) simulations have made great progress in characterizing proton transport, there has been a long-standing challenge in defining and tracking the excess proton, or more properly, the center of excess charge (CEC) created when a hydrogen nucleus distorts the electron distributions of water molecules in a delocalized and highly dynamic nature. Yet, defining (and biasing) such a CEC is essential when combining AIMD with enhanced sampling methods to calculate the relevant macroscopic properties via free-energy landscapes, which is the standard practice for most processes of interest. Several CEC formulas have been proposed and used, but none have yet been systematically tested or rigorously derived. In this paper, we show that the CEC can be used as a computational tool to disentangle IR features of the solvated excess proton from its surrounding solvent, and in turn, how correlating the features in the excess charge spectrum with the behavior of CEC in simulations enables a systematic evaluation of various CEC definitions. We present a new definition of CEC and show how it overcomes the limitations of those currently available both from a spectroscopic point of view and from a practical perspective of performance in enhanced sampling simulations.

M2 amphipathic helices facilitate pH-dependent conformational transition in influenza A virus

The matrix-2 (M2) protein from influenza A virus is a tetrameric, integral transmembrane (TM) protein that plays a vital role in viral replication by proton flux into the virus. The His37 tetrad is a pH sensor in the center of the M2 TM helix that activates the channel in response to the low endosomal pH. M2 consists of different regions that are believed to be involved in membrane targeting, packaging, nucleocapsid binding, and proton transport. Although M2 has been the target of many experimental and theoretical studies that have led to significant insights into its structure and function under differing conditions, the main mechanism of proton transport, its conformational dynamics, and the role of the amphipathic helices (AHs) on proton conductance remain elusive. To this end, we have applied explicit solvent constant pH molecular dynamics using the multisite λ-dynamics approach (CpHMD^MSλD) to investigate the buried ionizable residues comprehensively and to elucidate their effect on the conformational transition. Our model recapitulates the pH-dependent conformational transition of M2 from closed to open state when the AH domain is included in the M2 construct, revealing the role of the amphipathic helices on this transition and shedding light on the proton-transport mechanism. This work demonstrates the importance of including the amphipathic helices in future experimental and theoretical studies of ion channels. Finally, our work shows that explicit solvent CpHMD^MSλD provides a realistic pH-dependent model for membrane proteins.

Atomic structures of closed and open influenza B M2 proton channel reveal the conduction mechanism

The influenza B M2 (BM2) proton channel is activated by acidic pH to mediate virus uncoating. Unlike influenza A M2 (AM2), which conducts protons with strong inward rectification, BM2 conducts protons both inward and outward. Here we report 1.4- and 1.5-Å solid-state NMR structures of the transmembrane domain of the closed and open BM2 channels in a phospholipid environment. Upon activation, the transmembrane helices increase the tilt angle by 6° and the average pore diameter enlarges by 2.1 Å. BM2 thus undergoes a scissor motion for activation, which differs from the alternating-access motion of AM2. These results indicate that asymmetric proton conduction requires a backbone hinge motion, whereas bidirectional conduction is achieved by a symmetric scissor motion. The proton-selective histidine and gating tryptophan in the open BM2 reorient on the microsecond timescale, similar to AM2, indicating that side chain dynamics are the essential driver of proton shuttling.