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

2021 ◽  
Author(s):  
Laura C. Watkins ◽  
William F. DeGrado ◽  
Gregory A. Voth

ABSTRACTThe influenza A M2 channel, a prototype for the viroporin class of viral channels, is an acid-activated viroporin that conducts protons across the viral membrane, a critical step in the viral life cycle. As the protons enter from the viral exterior, four central His37 residues control the channel activation by binding subsequent protons, which opens the Trp41 gate and allows proton flux to the viral interior. Asp44 is essential for maintaining the Trp41 gate in a closed state at high pH, which results in asymmetric conduction. The prevalent D44N mutant disrupts this gate and opens the C-terminal end of the channel, resulting in overall increased conduction in the physiologically relevant pH range and a loss of this asymmetric conduction. Here, we use extensive Multiscale Reactive Molecular Dynamics (MS-RMD) and Quantum Mechanics/Molecular mechanics (QM/MM) simulations with an explicit, reactive excess proton to calculate the free energy of proton transport in the 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 up 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 our understanding of the role of emergent mutations in viroporins.

2020 ◽  
Author(s):  
Chenghan Li ◽  
Zhi Yue ◽  
L. Michel Espinoza-Fonseca ◽  
Gregory A. Voth

ABSTRACTThe sarcoplasmic reticulum Ca2+-ATPase (SERCA) transports two Ca2+ ions from the cytoplasm to the reticulum lumen at the expense of ATP hydrolysis. In addition to transporting Ca2+, 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 Ca2+. Previous studies have shown the existence of a transient water-filled pore in SERCA that connects the Ca2+-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 (MS-RMD) 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 Ca2+ transport. We propose that this proton transport mechanism is operational and serves as a functional conduit for passive proton transport across the sarcoplasmic reticulum.SIGNIFICANCEMultiscale reactive molecular dynamics combined with free energy sampling was applied to study proton transport through a transient water pore connecting the Ca2+-binding site to the lumen in SERCA. This is the first computational study of this large biomolecular system that treats the hydrated excess proton and its transport through water structures and amino acids explicitly. When also correctly accounting for the hydration fluctuations of the pore, it is found that a transiently hydrated channel can transport protons on a microsecond timescale. These results quantitatively support the hypothesis of the proton intake into the sarcoplasm via SERCA, in addition to the well-known proton pumping by SERCA to the cytoplasm along with Ca2+ transport.


2009 ◽  
Vol 10 (2) ◽  
pp. R18 ◽  
Author(s):  
Rachel Brower-Sinning ◽  
Donald M Carter ◽  
Corey J Crevar ◽  
Elodie Ghedin ◽  
Ted M Ross ◽  
...  

2016 ◽  
Vol 113 (45) ◽  
pp. E6955-E6964 ◽  
Author(s):  
Ruibin Liang ◽  
Jessica M. J. Swanson ◽  
Jesper J. Madsen ◽  
Mei Hong ◽  
William F. DeGrado ◽  
...  

The homotetrameric influenza A M2 channel (AM2) is an acid-activated proton channel responsible for the acidification of the influenza virus interior, an important step in the viral lifecycle. Four histidine residues (His37) in the center of the channel act as a pH sensor and proton selectivity filter. Despite intense study, the pH-dependent activation mechanism of the AM2 channel has to date not been completely understood at a molecular level. Herein we have used multiscale computer simulations to characterize (with explicit proton transport free energy profiles and their associated calculated conductances) the activation mechanism of AM2. All proton transfer steps involved in proton diffusion through the channel, including the protonation/deprotonation of His37, are explicitly considered using classical, quantum, and reactive molecular dynamics methods. The asymmetry of the proton transport free energy profile under high-pH conditions qualitatively explains the rectification behavior of AM2 (i.e., why the inward proton flux is allowed when the pH is low in viral exterior and high in viral interior, but outward proton flux is prohibited when the pH gradient is reversed). Also, in agreement with electrophysiological results, our simulations indicate that the C-terminal amphipathic helix does not significantly change the proton conduction mechanism in the AM2 transmembrane domain; the four transmembrane helices flanking the channel lumen alone seem to determine the proton conduction mechanism.


2016 ◽  
Author(s):  
Sangyun Lee ◽  
Heather B. Mayes ◽  
Jessica M. J. Swanson ◽  
Gregory A. Voth

AbstractThe ClC family of transmembrane proteins functions throughout nature to control the transport of Cl− ions across biological membranes. ClC-ec1 from Escherichia coli is an antiporter, coupling the transport of Cl− and H+ ions in opposite directions and driven by the concentration gradients of the ions. Despite keen interest in this protein, the molecular mechanism of the Cl−/H+ coupling has not been fully elucidated. Here, we have used multiscale simulation to help identify the essential mechanism of the Cl−/H+ coupling. We find that the highest barrier for proton transport (PT) from the intra- to extracellular solution is attributable to a chemical reaction—the deprotonation of glutamic acid 148 (E148). This barrier is significantly reduced by the binding of Cl− in the “central” site (Cl−cen), which displaces E148 and thereby facilitates its deprotonation. Conversely, in the absence of Cl−cen E148 favors the “down” conformation, which results in a much higher cumulative rotation and deprotonation barrier that effectively blocks PT to the extracellular solution. Thus, the rotation of E148 plays a critical role in defining the Cl−/H+ coupling. As a control, we have also simulated PT in the ClC-ec1 E148A mutant to further understand the role of this residue. Replacement with a non-protonatable residue greatly increases the free energy barrier for PT from E203 to the extracellular solution, explaining the experimental result that PT in E148A is blocked whether or not Cl−cen is present. The results presented here suggest both how a chemical reaction can control the rate of PT and also how it can provide a mechanism for a coupling of the two ion transport processes.


2020 ◽  
Vol 117 (7) ◽  
pp. 3583-3591 ◽  
Author(s):  
Hedieh Torabifard ◽  
Afra Panahi ◽  
Charles L. Brooks

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 (CpHMDMSλ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 CpHMDMSλD provides a realistic pH-dependent model for membrane proteins.


2016 ◽  
Vol 113 (27) ◽  
pp. 7420-7425 ◽  
Author(s):  
Ruibin Liang ◽  
Jessica M. J. Swanson ◽  
Yuxing Peng ◽  
Mårten Wikström ◽  
Gregory A. Voth

Cytochrome c oxidase (CcO) reduces oxygen to water and uses the released free energy to pump protons across the membrane. We have used multiscale reactive molecular dynamics simulations to explicitly characterize (with free-energy profiles and calculated rates) the internal proton transport events that enable proton pumping during first steps of oxidation of the fully reduced enzyme. Our results show that proton transport from amino acid residue E286 to both the pump loading site (PLS) and to the binuclear center (BNC) are thermodynamically driven by electron transfer from heme a to the BNC, but that the former (i.e., pumping) is kinetically favored whereas the latter (i.e., transfer of the chemical proton) is rate-limiting. The calculated rates agree with experimental measurements. The backflow of the pumped proton from the PLS to E286 and from E286 to the inside of the membrane is prevented by large free-energy barriers for the backflow reactions. Proton transport from E286 to the PLS through the hydrophobic cavity and from D132 to E286 through the D-channel are found to be strongly coupled to dynamical hydration changes in the corresponding pathways and, importantly, vice versa.


2021 ◽  
Vol 22 (5) ◽  
pp. 2409
Author(s):  
Anastasia A. Bizyaeva ◽  
Dmitry A. Bunin ◽  
Valeria L. Moiseenko ◽  
Alexandra S. Gambaryan ◽  
Sonja Balk ◽  
...  

Nucleic acid aptamers are generally accepted as promising elements for the specific and high-affinity binding of various biomolecules. It has been shown for a number of aptamers that the complexes with several related proteins may possess a similar affinity. An outstanding example is the G-quadruplex DNA aptamer RHA0385, which binds to the hemagglutinins of various influenza A virus strains. These hemagglutinins have homologous tertiary structures but moderate-to-low amino acid sequence identities. Here, the experiment was inverted, targeting the same protein using a set of related, parallel G-quadruplexes. The 5′- and 3′-flanking sequences of RHA0385 were truncated to yield parallel G-quadruplex with three propeller loops that were 7, 1, and 1 nucleotides in length. Next, a set of minimal, parallel G-quadruplexes with three single-nucleotide loops was tested. These G-quadruplexes were characterized both structurally and functionally. All parallel G-quadruplexes had affinities for both recombinant hemagglutinin and influenza virions. In summary, the parallel G-quadruplex represents a minimal core structure with functional activity that binds influenza A hemagglutinin. The flanking sequences and loops represent additional features that can be used to modulate the affinity. Thus, the RHA0385–hemagglutinin complex serves as an excellent example of the hypothesis of a core structure that is decorated with additional recognizing elements capable of improving the binding properties of the aptamer.


2021 ◽  
Vol 22 (12) ◽  
pp. 6458
Author(s):  
Aleksandra Pieniężna ◽  
Weronika Witak ◽  
Aneta Szymańska ◽  
Justyna Brasuń

In this paper, we present studies on the influence of the disulfide bridge on the copper (II) ions’ binding abilities by the cyclic His4-peptide. The studied ligand HKHPHRHC-S-S-C consists of nine amino acids. The cyclic structure was obtained through a disulfide bridge between two cysteinyl groups. Moreover, this peptide is characterized by the presence of four His residues in the sequence, which makes it an interesting ligand for transition metal ions. The potentiometric and spectroscopic (UV-Vis spectroscopy and circular dichroism spectroscopy (CD)) studies were carried out in various molar ligand to metal ratios: 2:1, 1:1, and 1:2, in the pH range of 2.5–11 at 25 °C. The results showed that the cyclic His4-peptide promotes dinuclear complexes in each of these systems and forms the final dinuclear species with the {NIm, 3N-amide}{NIm, 3N-amide} coordination mode. The obtained data shows that cyclization by the formation of the disulfide bond has an impact on the peptide chain flexibility and appearance of additional potential donors for metal ions and influences the copper (II) ions’ coordination.


Author(s):  
Nayara Dantas Coutinho ◽  
Hugo Gontijo Machado ◽  
Valter Henrique Carvalho-Silva ◽  
Wender A. Silva

Recent studies have assigned hydroxide elimination and C=C bond formation step in base-promoted aldol condensation the role of having a strong influence in the overall rate reaction, in contrast to...


Sign in / Sign up

Export Citation Format

Share Document