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

Exploring the structural and functional dynamics of trimeric and tetrameric states of influenza encoded PB1-F2 viroporin through molecular dynamics simulations

Influenza Viruses have always been a major health concern due to their highly contagious nature. The PB1-F2 viroporin encoded by the influenza A virus is known to be a pro-apoptotic protein involved in cell death induction of the host immune cells. The structural arrangement and the mode of action of PB1-F2 viroporin have not been fully understood yet. Nonetheless, there is limited information on the oligomeric state of PB1-F2 and its possible role in the pore formation which could act as a channel for ion transport. The probable oligomeric structural existences of the viroporin and their channel-like behavior need to be explored in light of experimental reports cited in the literature. In our study, we report on the structural and dynamical properties of the trimeric and tetrameric state of PB1-F2, investigated by molecular dynamics simulations with improved sampling of conformational states as the initial focus of the study is to establish a rationale for their existence in a lipid environment. The simulation study provides detailed information on the mitochondrial membrane permeation pathway which causes the leakage of mitochondrial contents like cytochrome C and induces apoptosis. By focusing on low-order oligomers, trimer, and tetramer, we have identified key pore-forming characteristics that serve as a foundation for understanding the pro-apoptotic activity of PB1-F2. The structural and dynamical properties of these states were evaluated in the light of experimental reports, which reveal the tetrameric form to be the preferable state in the lipid environment, demonstrating superior structural stability, effective channel symmetry, and ion permeation compared to the higher-order oligomers besides trimer including pentameric and hexameric assemblies. The simulation results also explore the typical ion transportation criteria based on finding a less energetic barrier for ions/water molecules crossing the membrane.

Reaction Mechanism Path Sampling Based on Parallel Cascade Selection QM/MM Molecular Dynamics Simulation: PaCS-Q

Quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) simulations are essential for elucidating complex biochemical reaction mechanisms. However, conventional enhanced sampling methods, such as umbrella sampling and metadynamics, often face limitations in computational cost, sampling completeness, and reliance on predefined reaction coordinates. To address these challenges, we developed Parallel Cascade Selection QM/MM MD (PaCS-Q) simulation, a novel strategy that efficiently explores reaction pathways by iteratively identifying high-potential structures for configurational transitions without predefined biases or external constraints. PaCS-Q directly tracks changes in bond distances over time, enabling the identification of transition states and intermediates. Validation of the Claisen rearrangement in chorismate mutase and the peptidyl aldehyde reaction in the Zika virus NS2B/NS3 serine protease demonstrated accurate pathway capture, reduced computational costs, and efficient sampling. With its user-friendly workflow, PaCS-Q broadens accessibility for computational and experimental researchers, offering a robust tool for studying enzymatic mechanisms with high accuracy and efficiency.

Activation of the Influenza B M2 Proton Channel (BM2)

Influenza B viruses have cocirculated 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 performed membrane-enabled continuous constant-pH molecular dynamics simulations on wildtype BM2 and a key H27A mutant channel to explore its pH-dependent conformational switch. Simulations captured the activation as the first histidine (His19) protonates and revealed the transition at lower pH values compared to AM2 is a result of electrostatic repulsions between His19 and preprotonated His27. Crucially, we provided an atomic-level understanding of the symmetric proton conduction by identifying preactivating 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 antiflu drug design efforts.

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.

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.

Ionic Liquid Gating Induces Anomalous Permeation through Membrane Channel Proteins

Membrane channel proteins (MCPs) play key roles in matter transport through cell membranes and act as major targets for vaccines and drugs. For emerging ionic liquid (IL) drugs, a rational understanding of how ILs affect the structure and transport function of MCP is crucial to their design. In this work, GPU-accelerated microsecond-long molecular dynamics simulations were employed to investigate the modulating mechanism of ILs on MCP. Interestingly, ILs prefer to insert into the lipid bilayer and channel of aquaporin-2 (AQP2) but adsorb on the entrance of voltage-gated sodium channels (Nav). Molecular trajectory and free energy analysis reflect that ILs have a minimal impact on the structure of MCPs but significantly influence MCP functions. It demonstrates that ILs can decrease the overall energy barrier for water through AQP2 by 1.88 kcal/mol, whereas that for Na+ through Nav is increased by 1.70 kcal/mol. Consequently, the permeation rates of water and Na+ can be enhanced and reduced by at least 1 order of magnitude, respectively. Furthermore, an abnormal IL gating mechanism was proposed by combining the hydrophobic nature of MCP and confined water/ion coordination effects. More importantly, we performed experiments to confirm the influence of ILs on AQP2 in human cells and found that treatment with ILs significantly accelerated the changes in cell volume in response to altered external osmotic pressure. Overall, these quantitative results will not only deepen the understanding of IL-cell interactions but may also shed light on the rational design of drugs and disease diagnosis.

Water, Protons, and the Gating of Voltage-Gated Potassium Channels

Ion channels are ubiquitous throughout all forms of life. Potassium channels are even found in viruses. Every cell must communicate with its surroundings, so all cells have them, and excitable cells, in particular, especially nerve cells, depend on the behavior of these channels. Every channel must be open at the appropriate time, and only then, so that each channel opens in response to the stimulus that tells that channel to open. One set of channels, including those in nerve cells, responds to voltage. There is a standard model for the gating of these channels that has a section of the protein moving in response to the voltage. However, there is evidence that protons are moving, rather than protein. Water is critical as part of the gating process, although it is hard to see how this works in the standard model. Here, we review the extensive evidence of the importance of the role of water and protons in gating these channels. Our principal example, but by no means the only example, will be the Kv1.2 channel. Evidence comes from the effects of D2O, from mutations in the voltage sensing domain, as well as in the linker between that domain and the gate, and at the gate itself. There is additional evidence from computations, especially quantum calculations. Structural evidence comes from X-ray studies. The hydration of ions is critical in the transfer of ions in constricted spaces, such as the gate region and the pore of a channel; we will see how the structure of the hydrated ion fits with the structure of the channel. In addition, there is macroscopic evidence from osmotic experiments and streaming current measurements. The combined evidence is discussed in the context of a model that emphasizes the role of protons and water in gating these channels.

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.

Functions of Viroporins in the Viral Life Cycle and Their Regulation of Host Cell Responses

Viroporins are virally encoded transmembrane proteins that are essential for viral pathogenicity and can participate in various stages of the viral life cycle, thereby promoting viral proliferation. Viroporins have multifaceted effects on host cell biological functions, including altering cell membrane permeability, triggering inflammasome formation, inducing apoptosis and autophagy, and evading immune responses, thereby ensuring that the virus completes its life cycle. Viroporins are also virulence factors, and their complete or partial deletion often reduces virion release and reduces viral pathogenicity, highlighting the important role of these proteins in the viral life cycle. Thus, viroporins represent a common drug-protein target for inhibiting drugs and the development of antiviral therapies. This article reviews current studies on the functions of viroporins in the viral life cycle and their regulation of host cell responses, with the aim of improving the understanding of this growing family of viral proteins.