Exploring the proton conductance and drug resistance of BM2 channel through molecular dynamics simulations and free energy calculations at different pH conditions

Understanding the conformational changes in the influenza B M2 ion channel at various protonation states

The characterization of influenza (A/B M2) ion channels is very important as they are potential binding sites for the drugs. We report the all-atom molecular dynamics study of the influenza B M2 ion channel in the presence of explicit solvent and lipid bilayers using the high resolution solid-state NMR structures. The importance of the various protonation states of histidine in the activation of the ion channel is discussed. The conformational changes at the closed and the open structures clearly show that the increase in tilt angle is necessary for the activation of the ion channel. Additionally, the free energy surfaces of the eight systems show the importance of the protonation state of the histidine residues in the activation of the influenza B M2 ion channel. The protonation of the histidine residues increases the tilt angle and the intra-helix distance which is evident from the superimposition of the structures corresponding to the maxima and the minima in the free energy landscape. The findings imply differences in the singly protonated and double protonated conformational states of BM2 ion channel and provide insights to help further studies of these ion channels as the drug targets for the influenza virus.

Helix-coil transition and conformational deformity in Aβ42-monomer: a case study using the Zn2+ cation

The metal ions (like Fe2+, Zn²⁺, Cu2+) are known to influence the amyloid beta (Aβ) aggregation. In this study, we have examined the conformational and dynamical changes during the coordination of Aβ-monomer with the Zn2+ ion using all-atom molecular dynamics (MD) simulations using explicit solvent models. We have probed the unfolding of the full-length Aβ₄₂ monomer both inclusive and exclusive of the Zn2+ cation, with 1:1 ratio of the peptide and the Zn2+ cation. The inclusion of the Zn2+ cation shows differential intra-peptide interactions which has been probed using various analyses. The Helix - Coil transition of the wild type Aβ42 monomer is studied using the steered molecular dynamics simulations by taking the end-to-end C-α distance across the peptide. This gives an idea of the unequal intra - peptide and peptide - water interactions being found across the length of the Aβ monomer. The transition of an α-helix dominated wild-type (WT) Aβ structure to the unfolded coil structure gives significant evidence of the intra-peptide hydrogen bonding shifts in the presence of the Zn2+ cation. This accounts for the structural and the dynamical variations that take place in the Aβ monomer in the presence of the Zn2+ cation to mimic the conditions/environment at the onset of fibrillation.Communicated by Ramaswamy H. Sarma.

In Silico Analysis Revealed a Unique Binding but Ineffective Mode of Amantadine to Influenza Virus B M2 Channel

The M2 proton channel of influenza A (AM2) and B (BM2) have a highly conserved function motif, considered as the effective target. As yet, there is no effective drug against BM2. Research showed that AM2 channel blocker, amantadine (AMT), was able to bind to BM2 channel, but AMT lacked inhibition against BM2. Nevertheless, the study of the binding but ineffective mode of AMT to BM2 is challenging. To resolve the challenge and obtain more information for drug design of inhibitors targeting BM2, multiple molecular dynamics simulations were performed. We discovered AMT mainly adopted up binding mode in BM2, involved in a transition flipping from down mode to up mode. Furthermore, we discovered a new key factor to explain ineffective inhibition of AMT to BM2 because of the unmatched spatial geometry between AMT and BM2. Our work could enrich structural feature information on BM2 and provide a new perspective for rational drug design of anti-influenza B.

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.

What are the effects of the serine triad on proton conduction of an influenza B M2 channel? An investigation by molecular dynamics simulations

The tetrameric influenza B M2 channel (BM2), an acid activated proton channel, is important in the influenza virus B lifecycle. A conserved HxxxW motif is responsible for proton conduction and channel gating. In this study, to explore the effects of the serine triad (S9, S12 and S16) on proton conduction, we performed classical molecular dynamics (CMD) simulations and adaptive steered molecular dynamics (ASMD) simulations at different protonation states of the H19 tetrad. The results of the pore radius and the C-terminal tilt angle show that the electrostatic repulsion induced by protonated H19 is the key driving force for opening the BM2 channel. The open states could be stabilized by the hydrogen bonds between S16 and protonated H19. The solvent accessible surface area and water density indicate that the polar hydrophilic environment provided by the serine triad facilitates the formation of a water wire, and then exhibits favourable effects on proton conduction. The mutant research verifies and supports these views. Our work clarifies the effects of the serine triad on proton conduction in the BM2 channel, which would help us deeply understand the proton conduction mechanism in BM2 and provides a new perspective for antiviral drug design against BM2.

A unique activation-promotion mechanism of the influenza B M2 proton channel uncovered by multiscale simulations

The influenza B M2 protein (BM2) forms an acid-activated proton channel that is important for the virus's lifecycle. Despite extensive research efforts, the detailed activation mechanism of the BM2 proton channel is often elusive. Herein a pH-regulated mechanism of the BM2TM domain has been systematically characterized using multiscale computer simulations, including classical molecular dynamics, constant pH molecular dynamics (CpHMD) and quantum mechanics/molecular mechanics (QM/MM) approaches. Our simulations reveal a pH-dependent conformational switch from the C-terminal closed to the C-terminal open conformers, and provide the free energy of conformational activation coupled to the titration of the His19 tetrad. Importantly, our results confirm the coupling titration between the His19 tetrad and His27 tetrad, and identify that the full-cationic state (His2744+) dominates at the low pH (the His19 tetrad at +2, +3 and +4 charge states). Our QM/MM simulations indicate that the second titratable histidine, His27, could further promote the BM2 acid activation and speed up proton dissociation from the HxxxW motif, thus facilitating proton conduction by BM2. Taken together, a unique "activation-promotion mechanism" about the BM2 proton channel is proposed, and these results may be helpful for the understanding of other similar proton channels and the development of BM2 inhibitors.

Viral channel forming proteins--How to assemble and depolarize lipid membranes in silico

Viral channel forming proteins (VCPs) have been discovered in the late 70s and are found in many viruses to date. Usually they are small and have to assemble to form channels which depolarize the lipid membrane of the host cells. Structural information is just about to emerge for just some of them. Thus, computational methods play a pivotal role in generating plausible structures which can be used in the drug development process. In this review the accumulation of structural data is introduced from a historical perspective. Computational performances and their predictive power are reported guided by biological questions such as the assembly, mechanism of function and drug–protein interaction of VCPs. An outlook of how coarse grained simulations can contribute to yet unexplored issues of these proteins is given. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

•

Early references about the discovery of viral channel forming proteins.

•

Latest structural information about the class of proteins.

•

Identification of structural motifs, assembly mechanism of function and drug action using computational methods.

•

Outlook for the use of coarse grained techniques to address assembly and integration into cellular processes.

Molecular dynamics, monte carlo simulations, and langevin dynamics: a computational review

Macromolecular structures, such as neuraminidases, hemagglutinins, and monoclonal antibodies, are not rigid entities. Rather, they are characterised by their flexibility, which is the result of the interaction and collective motion of their constituent atoms. This conformational diversity has a significant impact on their physicochemical and biological properties. Among these are their structural stability, the transport of ions through the M2 channel, drug resistance, macromolecular docking, binding energy, and rational epitope design. To assess these properties and to calculate the associated thermodynamical observables, the conformational space must be efficiently sampled and the dynamic of the constituent atoms must be simulated. This paper presents algorithms and techniques that address the abovementioned issues. To this end, a computational review of molecular dynamics, Monte Carlo simulations, Langevin dynamics, and free energy calculation is presented. The exposition is made from first principles to promote a better understanding of the potentialities, limitations, applications, and interrelations of these computational methods.

Mechanistic insight into the functional transition of the enzyme guanylate kinase induced by a single mutation

Dramatic functional changes of enzyme usually require scores of alterations in amino acid sequence. However, in the case of guanylate kinase (GK), the functional novelty is induced by a single (S→P) mutation, leading to the functional transition of the enzyme from a phosphoryl transfer kinase into a phosphorprotein interaction domain. Here, by using molecular dynamic (MD) and metadynamics simulations, we provide a comprehensive description of the conformational transitions of the enzyme after mutating serine to proline. Our results suggest that the serine plays a crucial role in maintaining the closed conformation of wild-type GK and the GMP recognition. On the contrary, the S→P mutant exhibits a stable open conformation and loses the ability of ligand binding, which explains its functional transition from the GK enzyme to the GK domain. Furthermore, the free energy profiles (FEPs) obtained by metadymanics clearly demonstrate that the open-closed conformational transition in WT GK is positive correlated with the process of GMP binding, indicating the GMP-induced closing motion of GK enzyme, which is not observed in the mutant. In addition, the FEPs show that the S→P mutation can also leads to the mis-recognition of GMP, explaining the vanishing of catalytic activity of the mutant.

Interactions of lipids and detergents with a viral ion channel protein: molecular dynamics simulation studies

Structural studies of membrane proteins have highlighted the likely influence of membrane mimetic environments (i.e., lipid bilayers versus detergent micelles) on the conformation and dynamics of small α-helical membrane proteins. We have used molecular dynamics simulations to compare the conformational dynamics of BM2 (a small α-helical protein from the membrane of influenza B) in a model phospholipid bilayer environment with its behavior in protein–detergent complexes with either the zwitterionic detergent dihexanoylphosphatidylcholine (DHPC) or the nonionic detergent dodecylmaltoside (DDM). We find that DDM more closely resembles the lipid bilayer in terms of its interaction with the protein, while the short-tailed DHPC molecule forms “nonphysiological” interactions with the protein termini. We find that the intrinsic micelle properties of each detergent are conserved upon formation of the protein–detergent complex. This implies that simulations of detergent micelles may be used to help select optimal conditions for experimental studies of membrane proteins.

COMPARING FOLDING MECHANISMS OF DIFFERENT PRION PROTEINS BY Gō MODEL

Prions are associated with neurodegenerative diseases induced by transmissible spongiform encephalopathies. The infectious scrapie form is referred to as PrP Sc , which has conformational change from normal prion with predominant α-helical conformation to the abnormal PrP Sc that is rich in β-sheet content. Neurodegenerative diseases have been found from both human and bovine sources, but there are no reports about infected by transmissible spongiform encephalopathies from rabbit, canine and horse sources. Here we used coarse-grained Gō model to compare the difference among human, bovine, rabbit, canine, and horse normal (cellular) prion proteins. The denatured state of normal prion has relation with the conversion from normal to abnormal prion protein, so we used all-atom Gō model to investigate the folding pathway and energy landscape for human prion protein. Through using coarse-grained Gō model, the cooperativity of the five prion proteins was characterized in terms of calorimetric criterion, sigmoidal transition, and free-energy profile. The rabbit and horse prion proteins have higher folding free-energy barrier and cooperativity, and canine prion protein has slightly higher folding free-energy barrier comparing with human and bovine prion proteins. The results from all-atom Gō model confirmed the validity of C α-Gō model. The correlations of our results with previous experimental and theoretical researches were discussed.

UNDERSTANDING THE MOLECULAR MECHANISM OF BINDING MODES OF AURORA A INHIBITORS BY LONG TIME SCALE GPU DYNAMICS

Inhibition of Aurora A kinase interaction is considered to be a promising approach for the discovery of new molecularly targeted cancer therapeutics. In this study, the binding mechanisms of two different inhibitors with a contrasting binding affinity to Aurora A were investigated by long time scale GPU molecular dynamics (MD) simulations coupled with molecular mechanics-Poisson–Boltzmann/generalized Born surface area (MM-PB/GBSA) method. The results showed that the predicted binding free energies of these two complexes were consistent with the experimental data. Through analyzing the individual energy components of binding free energy, we found that the van der Waals contribution was the main force to drive the inhibitor–protein binding and the electrostatic contribution was also a crucial factor for the inhibitor–Aurora A binding. The structural analysis demonstrated that the inhibitor HPM could produce more hydrophobic interaction contacts with Aurora A than that of 2JZ, and the loss of key hydrogen bonds between the inhibitor and residue Arg137 in the hinge region of Aurora A was another important reason for the weaker binding affinity of 2JZ to Aurora A. This study sheds more light on the development of the efficient inhibitors targeting the Aurora A.

Variation in one residue associated with the metal ion-dependent adhesion site regulates αIIbβ3 integrin ligand binding affinity

The Asp of the RGD motif of the ligand coordinates with the β I domain metal ion dependent adhesion site (MIDAS) divalent cation, emphasizing the importance of the MIDAS in ligand binding. There appears to be two distinct groups of integrins that differ in their ligand binding affinity and adhesion ability. These differences may be due to a specific residue associated with the MIDAS, particularly the β3 residue Ala²⁵² and corresponding Ala in the β1 integrin compared to the analogous Asp residue in the β2 and β7 integrins. Interestingly, mutations in the adjacent to MIDAS (ADMIDAS) of integrins α4β7 and αLβ2 increased the binding and adhesion abilities compared to the wild-type, while the same mutations in the α2β1, α5β1, αVβ3, and αIIbβ3 integrins demonstrated decreased ligand binding and adhesion. We introduced a mutation in the αIIbβ3 to convert this MIDAS associated Ala²⁵² to Asp. By combination of this mutant with mutations of one or two ADMIDAS residues, we studied the effects of this residue on ligand binding and adhesion. Then, we performed molecular dynamics simulations on the wild-type and mutant αIIbβ3 integrin β I domains, and investigated the dynamics of metal ion binding sites in different integrin-RGD complexes. We found that the tendency of calculated binding free energies was in excellent agreement with the experimental results, suggesting that the variation in this MIDAS associated residue accounts for the differences in ligand binding and adhesion among different integrins, and it accounts for the conflicting results of ADMIDAS mutations within different integrins. This study sheds more light on the role of the MIDAS associated residue pertaining to ligand binding and adhesion and suggests that this residue may play a pivotal role in integrin-mediated cell rolling and firm adhesion.

Molecular dynamics simulation of tryptophan hydroxylase-1: binding modes and free energy analysis to phenylalanine derivative inhibitors

Serotonin is a neurotransmitter that modulates many central and peripheral functions. Tryptophan hydroxylase-1 (TPH1) is a key enzyme of serotonin synthesis. In the current study, the interaction mechanism of phenylalanine derivative TPH1 inhibitors was investigated using molecular dynamics (MD) simulations, free energy calculations, free energy decomposition analysis and computational alanine scanning. The predicted binding free energies of these complexes are consistent with the experimental data. The analysis of the individual energy terms indicates that although the van der Waals and electrostatics interaction contributions are important in distinguishing the binding affinities of these inhibitors, the electrostatic contribution plays a more crucial role in that. Moreover, it is observed that different configurations of the naphthalene substituent could form different binding patterns with protein, yet lead to similar inhibitory potency. The combination of different molecular modeling techniques is an efficient way to interpret the interaction mechanism of inhibitors and our work could provide valuable information for the TPH1 inhibitor design in the future.

Theoretical investigation on the diatomic ligand migration process and ligand binding properties in non-O2-binding H-NOX domain

The Nostoc sp (Ns) H-NOX (heme-nitric oxide or OXygen-binding) domain shares 35% sequence identity with soluble guanylate cyclase (sGC) and exhibits similar ligand binding property with the sGC. Previously, our molecular dynamic (MD) simulation work identified that there exists a Y-shaped tunnel system hosted in the Ns H-NOX interior, which servers for ligand migration. The tunnels were then confirmed by Winter et al. [PNAS 2011;108(43):E 881-889] recently using x-ray crystallography with xenon pressured conditions. In this work, to further investigate how the protein matrix of Ns H-NOX modulates the ligand migration process and how the distal residue composition affects the ligand binding prosperities, the free energy profiles for nitric oxide (NO), carbon monooxide (CO), and O2 migration are explored using the steered MDs simulation and the ligand binding energies are calculated using QM/MM schemes. The potential of mean force profiles suggest that the longer branch of the tunnel would be the most favorable route for NO migration and a second NO trapping site other than the distal heme pocket along this route in the Ns H-NOX was identified. On the contrary, CO and O2 would prefer to diffuse via the shorter branch of the tunnel. The QM/MM (quantum mechanics/molecular mechanics) calculations suggest that the hydrophobic distal pocket of Ns H-NOX would provide an approximately vacuum environment and the ligand discrimination would be determined by the intrinsic binding properties of the diatomic gas ligand to the heme group.