Panita Decha

How amantadine and rimantadine inhibit proton transport in the M2 protein channel

To understand how antiviral drugs inhibit the replication of influenza A virus via the M2 ion channel, molecular dynamics simulations have been applied to the six possible protonation states of the M2 ion channel in free form and its complexes with two commercial drugs in a fully hydrated lipid bilayer. Among the six different states of free M2 tetramer, water density was present in the pore of the systems with mono-protonated, di-protonated at adjacent position, tri-protonated and tetra-protonated systems. In the presence of inhibitor, water density in the channel was considerably better reduced by rimantadine than amantadine, agreed well with the experimental IC(50) values. With the preferential position and orientation of the two drugs in all states, two mechanisms of action, where the drug binds to the opening pore and the histidine gate, were clearly explained, i.e., (i) inhibitor was detected to localize slightly closer to the histidine gate and can facilitate the orientation of His37 imidazole rings to lie in the close conformation and (ii) inhibitor acts as a blocker, binding at almost above the opening pore and interacts slightly with the three pore-lining residues, Leu26, Ala30 and Ser31. Here, the inhibitors were found to bind very weakly to the channel due to their allosteric hindrance while theirs side chains were strongly solvated.

Understanding of known drug-target interactions in the catalytic pocket of neuraminidase subtype N1.

To provide detailed information and insight into the drug‐target interaction, structure, solvation, and dynamic and thermodynamic properties, the three known‐neuraminidase inhibitors—oseltamivir (OTV), zanamivir (ZNV), and peramivir (PRV)—embedded in the catalytic site of neuraminidase (NA) subtype N1 were studied using molecular dynamics simulations. In terms of ligand conformation, there were major differences in the structures of the guanidinium and the bulky groups. The atoms of the guanidinium group of PRV were observed to form many more hydrogen bonds with the surrounded residues and were much less solvated by water molecules, in comparison with the other two inhibitors. Consequently, D151 lying on the 150‐loop (residues 147–152) of group‐1 neuraminidase (N1, N4, N5, and N8) was considerably shifted to form direct hydrogen bonds with the OH group of the PRV, which was located rather far from the 150‐loop. For the bulky group, direct hydrogen bonds were detected only between the hydrophilic side chain of ZNV and residues R224, E276, and E277 of N1 with rather weak binding, 20–70% occupation. This is not the case for OTV and PRV, in which flexibility and steric effects due to the hydrophobic side chain lead to the rearrangement of the surrounded residues, that is, the negatively charged side chain of E276 was shifted and rotated to form hydrogen bonds with the positively charged moiety of R224. Taking into account all the ligand‐enzyme interaction data, the gas phase MM interaction energy of −282.2 kcal/mol as well as the binding free energy (ΔGbinding) of −227.4 kcal/mol for the PRV‐N1 are significantly lower than those of the other inhibitors. The ordering of ΔGbinding of PRV < ZNV < OTV agrees well with the ordering of experimental IC50 value. Proteins 2008.

Evolution of Human Receptor Binding Affinity of H1N1 Hemagglutinins from 1918 to 2009 Pandemic Influenza A Virus

The recent outbreak of the novel 2009 H1N1 influenza in humans has focused global attention on this virus, which could potentially have introduced a more dangerous pandemic of influenza flu. In the initial step of the viral attachment, hemagglutinin (HA), a viral glycoprotein surface, is responsible for the binding to the human SIA alpha2,6-linked sialopentasaccharide host cell receptor (hHAR). Dynamical and structural properties, based on molecular dynamics simulations of the four different HAs of Spanish 1918 (H1-1918), swine 1930 (H1-1930), seasonal 2005 (H1-2005), and a novel 2009 (H1-2009) H1N1 bound to the hHAR were compared. In all four HA-hHAR complexes, major interactions with the receptor binding were gained from HA residue Y95 and the conserved HA residues of the 130-loop, 190-helix, and 220-loop. However, introduction of the charged HA residues K145 and E227 in the 2009 HA binding pocket was found to increase the HA-hHAR binding efficiency in comparison to the three previously recognized H1N1 strains. Changing of the noncharged HA G225 residue to a negatively charged D225 provides a larger number of hydrogen-bonding interactions. The increase in hydrophilicity of the receptor binding region is apparently an evolution of the current pandemic flu from the 1918 Spanish, 1930 swine, and 2005 seasonal strains. Detailed analysis could help the understanding of how different HAs effectively attach and bind with the hHAR.

Why amantadine loses its function in influenza m2 mutants: MD simulations.

Molecular dynamics simulations of the drug-resistant M2 mutants, A30T, S31N, and L26I, were carried out to investigate the inhibition of M2 activity using amantadine (AMT). The closed and open channel conformations were examined via non- and triply protonated H37. For the nonprotonated state, these mutants exhibited zero water density in the conducting region, and AMT was still bound to the channel pore. Thus, water transport is totally suppressed, similar to the wild-type channel. In contrast, the triply protonated states of the mutants exhibited a different water density and AMT position. A30T and L26I both have a greater water density compared to the wild-type M2, while for the A30T system, AMT is no longer inside the pore. Hydrogen bonding between AMT and H37 crucial for the bioactivity is entirely lost in the open conformation. The elimination of this important interaction of these mutations is responsible for the lost of AMTs function in influenza A M2. This is different for the S31N mutant in which AMT was observed to locate at the pore opening region and bond with V27 instead of S31.

Source of High Pathogenicity of an Avian Influenza Virus H5N1: Why H5 Is Better Cleaved by Furin

The origin of the high pathogenicity of an emerging avian influenza H5N1 due to the -RRRKK- insertion at the cleavage loop of the hemagglutinin H5, was studied using the molecular dynamics technique, in comparison with those of the noninserted H5 and H3 bound to the furin (FR) active site. The cleavage loop of the highly pathogenic H5 was found to bind strongly to the FR cavity, serving as a conformation suitable for the proteolytic reaction. With this configuration, the appropriate interatomic distances were found for all three reaction centers of the enzyme-substrate complex: the arrangement of the catalytic triad, attachment of the catalytic Ser(368) to the reactive S1-Arg, and formation of the oxyanion hole. Experimentally, the--RRRKK--insertion was also found to increase in cleavage of hemagglutinin by FR. The simulated data provide a clear answer to the question of why inserted H5 is better cleaved by FR than the other subtypes, explaining the high pathogenicity of avian influenza H5N1.

Combined QM/MM mechanistic study of the acylation process in furin complexed with the H5N1 avian influenza virus hemagglutinin's cleavage site.

Combined quantum mechanical/molecular mechanical (QM/MM) techniques have been applied to investigate the detailed reaction mechanism of the first step of the acylation process by furin in which the cleavage site of the highly pathogenic avian influenza virus subtype H5N1 (HPH5) acts as its substrate. The energy profile shows a simultaneous mechanism, known as a concerted reaction, of the two subprocesses: the proton transfer from Ser368 to His194 and the nucleophilic attack on the carbonyl carbon of the scissile peptide of the HPH5 cleavage site with a formation of tetrahedral intermediate (INT). The calculated energy barrier for this reaction is 16.2 kcal·mol−1 at QM/MM B3LYP/6‐31+G*//PM3‐CHARMM22 level of theory. Once the reaction proceeds, the ordering of the electrostatic stabilization by protein environment is of the enzyme‐substrate < transition state < INT complexes. Asp153 was found to play the most important role in the enzymatic reaction by providing the highest degree of intermediate complex stabilization. In addition, the negatively charged carbonyl oxygen of INT is well stabilized by the oxyanion hole constructed by Asn295s carboxamide and Ser368s backbone. Proteins 2009.

Theoretical studies on the molecular basis of HIV-1RT/NNRTIs interactions

Molecular dynamics simulations (MD) of the human immunodeficiency virus type 1 reverse transcriptase (HIV-1 RT) complexed with the four non-nucleoside reverse transcriptase inhibitors (NNRTIs): efavirenz (EFV), emivirine (EMV), etravirine (ETV) and nevirapine (NVP), were performed to examine the structures, binding free energies and the importance of water molecules in the binding site. The binding free energy, calculated using molecular mechanics Poisson-Boltzmann surface area (MM-PBSA), was found to decrease in the following order: EFV ∼ ETV > EMV > NVP. The decrease in stability of the HIV-1 RT/NNRTI complexes is in good agreement with the experimentally derived half maximal inhibitory concentration (IC50) values. The interaction energy of the protein-inhibitor complexes was found to be essentially associated with the cluster of seven hydrophobic residues, L100, V106, Y181, Y188, F227, W229 and P236, and two basic residues, K101 and K103. Moreover, these residues are considered to be the most frequently detected mutated amino acids during treatment by various NNRTIs and therefore, those most likely to have been selected in the population for resistance.

Evaluating how rimantadines control the proton gating of the influenza A M2-proton port via allosteric binding outside of the M2-channel: MD simulations

In order to understand how rimantadine (RMT) inhibits the proton conductance in the influenza A M2 channel via the recently proposed “allosteric mechanism”, molecular dynamics simulations were applied to the M2-tetrameric protein with four RMTs bound outside the channel at the three protonation states: the 0H-closed, 1H-intermediate and 3H-open situations. In the 0H-closed state, a narrow channel with the RMT-Asp44-Trp41 H-bond network was formed, therefore the water penetration through the channel was completely blocked. The Trp41-Asp44 interaction was absent in the 1H-intermediate state, whilst the binding of RMT to Asp44 remained, which resulted in a weakened helix-helix packing, therefore the channel was partially prevented. In the 3H-open state it was found that the electrostatic repulsion from the three charged His37 residues allowed the Trp41 gate to open, permitting water to penetrate through the channel. This agreed well with the potential of the means force which is in the following order: 0H > 1H > 3H.

Understanding of Drug‐Target Interactions: A case Study in Influenza Virus A Subtype H5N1

This study aims at gaining insight into molecular mechanisms of action of three drug targets of the life cycle of influenza virus A subtype H5N1, namely Hemagglutinin (H5), Neuraminidase (N1) and M2 ion channel (M2), using molecular mechanics and molecular dynamics techniques. In hemagglutinin, interest is focused on the high pathogenicity of the H5 due to the –RRRKK– insertion. MD simulations carried out for H5 in both high and low pathogenic forms (HPH5 and LPH5), aimed at understanding why HPH5 was experimentally observed to be 5‐fold better cleaved by furin relative to the non‐inserted sequence of LPH5. As the results, the cleavage loop of HPH5 was found to fit well and bind strongly into the catalytic site of human furin, serving as a conformation suitable for the proteolytic reaction. The second target, neuraminidase was studied by two different approaches. Firstly with MD simulations, rotation of the –NHAc and—OCHEt2 side chains of oseltamivir (OTV), leading directly to rearrangement of the catalyt...

Comment on "Cleavage mechanism of the H5N1 hemagglutinin by trypsin and furin" [Amino Acids 2008, January 31, Doi: 10.1007/s00726-007-0611-3].

Summary.Recently, Guo et al. have reported structural as well as the binding energy data of the particular interactions between the cleavage sites of hemagglutinin and serine proteases, trypsin and furin, using molecular docking approach. Due to a wrong assignment of protonation state on the histidine, one of the catalytic triad in the active site of both enzymes, their docking results are contradictory with the fundamental principle and previous theoretical studies of the known cleavage mechanism in serine proteases.