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 (IC(50)) 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.

The activated state of a sodium channel voltage sensor in a membrane environment.

Direct structural insights on the fundamental mechanisms of permeation, selectivity, and gating remain unavailable for the Na(+) and Ca(2+) channel families. Here, we report the spectroscopic structural characterization of the isolated Voltage-Sensor Domain (VSD) of the prokaryotic Na(+) channel NaChBac in a lipid bilayer. Site-directed spin-labeling and EPR spectroscopy were carried out for 118 mutants covering all of the VSD. EPR environmental data were used to unambiguously assign the secondary structure elements, define membrane insertion limits, and evaluate the activated conformation of the isolated-VSD in the membrane using restrain-driven molecular dynamics simulations. The overall three-dimensional fold of the NaChBac-VSD closely mirrors those seen in KvAP, Kv1.2, Kv1.2-2.1 chimera, and MlotiK1. However, in comparison to the membrane-embedded KvAP-VSD, the structural dynamics of the NaChBac-VSD reveals a much tighter helix packing, with subtle differences in the local environment of the gating charges and their interaction with the rest of the protein. Using cell complementation assays we show that the NaChBac-VSD can provide a conduit to the transport of ions in the resting or "down" conformation, a feature consistent with our EPR water accessibility measurements in the activated or "up" conformation. These results suggest that the overall architecture of VSD's is remarkably conserved among K(+) and Na(+) channels and that pathways for gating-pore currents may be intrinsic to most voltage-sensors. Cell complementation assays also provide information about the putative location of the gating charges in the "down/resting" state and hence a glimpse of the extent of conformational changes during activation.

How does each substituent functional group of oseltamivir lose its activity against virulent H5N1 influenza mutants?

To reveal the source of oseltamivir-resistance in influenza (A/H5N1) mutants, the drug-target interactions at each functional group were investigated using MD/LIE simulations. Oseltamivir in the H274Y mutation primarily loses the electrostatic and the vdW interaction energies at the -NH(3)(+) and -OCHEt(2) moieties corresponding to the weakened hydrogen-bonds and changed distances to N1 residues. Differentially, the N294S mutation showed small changes of binding energies and intermolecular interactions. Interestingly, the presence of different conformations of E276 positioned between the -OCHEt(2) group and the mutated residue is likely to play an important role in oseltamivir-resistant identification. In the H274Y mutant, it moves towards the -OCHEt(2) group leading to a reduction in hydrophobicity and pocket size, whilst in the N294S mutant it acts as the hydrogen network center bridging with R224 and the mutated residue S294. The molecular details have answered a question of how the H274Y and N294S mutations confer the high- and medium-level of oseltamivir-resistance to H5N1.

Susceptibility of antiviral drugs against 2009 influenza A (H1N1) virus.

The recent outbreak of the novel strain of influenza A (H1N1) virus has raised a global concern of the future risk of a pandemic. To understand at the molecular level how this new H1N1 virus can be inhibited by the current anti-influenza drugs and which of these drugs it is likely to already be resistant to, homology modeling and MD simulations have been applied on the H1N1 neuraminidase complexed with oseltamivir, and the M2-channel with adamantanes bound. The H1N1 virus was predicted to be susceptible to oseltamivir, with all important interactions with the binding residues being well conserved. In contrast, adamantanes are not predicted to be able to inhibit the M2 function and have completely lost their binding with the M2 residues. This is mainly due to the fact that the M2 transmembrane of the new H1N1 strain contains the S31N mutation which is known to confer resistance to adamantanes.

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 AMT's 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 oseltamivir resistance in avian influenza H5N1 virus with the H274Y mutation.

Molecular dynamics simulations were carried out for the mutant oseltamivir-NA complex, to provide detailed information on the oseltamivir-resistance resulting from the H274Y mutation in neuraminidase (NA) of avian influenza H5N1 viruses. In contrast with a previous proposal, the H274Y mutation does not prevent E276 and R224 from forming the hydrophobic pocket for the oseltamivir bulky group. Instead, reduction of the hydrophobicity and size of pocket in the area around an ethyl moiety at this bulky group were found to be the source of the oseltamivir-resistance. These changes were primarily due to the dramatic rotation of the hydrophilic -COO(-) group of E276 toward the ethyl moiety. In addition, hydrogen-bonding interactions with N1 residues at the -NH(3) (+) and -NHAc groups of oseltamivir were replaced by a water molecule. The calculated binding affinity of oseltamivir to NA was significantly reduced from -14.6 kcal mol(-1) in the wild-type to -9.9 kcal mol(-1) in the mutant-type.

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 Asn295's carboxamide and Ser368's backbone.

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.

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.

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 (DeltaG(binding)) of -227.4 kcal/mol for the PRV-N1 are significantly lower than those of the other inhibitors. The ordering of DeltaG(binding) of PRV < ZNV < OTV agrees well with the ordering of experimental IC(50) value.

Ab initio QM/MM dynamics of H3O+ in water.

A molecular dynamics (MD) simulation based on a combined ab initio quantum mechanics/molecular mechanics (QM/MM) method has been performed to investigate the solvation structure and dynamics of H3O+ in water. The QM region is a sphere around the central H3O+ ion, and contains about 6-8 water molecules. It is treated at the Hartree-Fock (HF) level, while the rest of the system is described by means of classical pair potentials. The Eigen complex (H9O4+) is found to be the most prevalent species in the aqueous solution, partly due to the selection scheme of the center of the QM region. The QM/MM results show that the Eigen complex frequently converts back and forth into the Zundel (H5O2+) structure. Besides the three nearest-neighbor water molecules directly hydrogen-bonded to H3O+, other neighbor waters, such as a fourth water molecule which interacts preferentially with the oxygen atom of the hydronium ion, are found occasionally near the ion. Analyses of the water exchange processes and the mean residence times of water molecules in the ion's hydration shell indicate that such next-nearest neighbor water molecules participate in the rearrangement of the hydrogen bond network during fluctuative formation of the Zundel ion and, thus, contribute to the Grotthuss transport of the proton.

Structure and dynamics of hydrated NH(4) (+): an ab initio QM/MM molecular dynamics simulation.

A combined ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulation has been performed to investigate solvation structure and dynamics of NH(4) (+) in water. The most interesting region, the sphere includes an ammonium ion and its first hydration shell, was treated at the Hartree-Fock level using DZV basis set, while the rest of the system was described by classical pair potentials. On the basis of detailed QM/MM simulation results, the solvation structure of NH(4) (+) is rather flexible, in which many water molecules are cooperatively involved in the solvation shell of the ion. Of particular interest, the QM/MM results show fast translation and rotation of NH(4) (+) in water. This phenomenon has resulted from multiple coordination, which drives the NH(4) (+) to translate and rotate quite freely within its surrounding water molecules. In addition, a "structure-breaking" behavior of the NH(4) (+) is well reflected by the detailed analysis on the water exchange process and the mean residence times of water molecules surrounding the ion.