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.     

Characteristics of CO3(2-)-water hydrogen bonds in aqueous solution: insights from HF/MM and B3LYP/MM MD simulations

Two combined quantum mechanics/molecular mechanics (QM/MM) molecular dynamics simulations, namely HF/MM and B3LYP/MM, have been performed to investigate the local hydration structure and dynamics of carbonate (CO(3)(2-)) in dilute aqueous solution. With respect to the QM/MM scheme, the QM region, which contains the CO(3)(2-) and its surrounding water molecules, was treated at HF and B3LYP levels of accuracy, respectively, using the DZV+ basis set, while the rest of the system is described by classical MM potentials. For both the HF/MM and B3LYP/MM simulations, it is observed that the hydrogen bonds between CO(3)(2-) oxygens and their nearest-neighbor waters are relatively strong, i.e., compared to water-water hydrogen bonds in the bulk, and that the first shell of each CO(3)(2-) oxygen atom somewhat overlaps with the others, which allows migration of water molecules among the coordinating sites to exist. In addition, it is observed that first-shell waters are either "loosely" or "tightly" bound to the respective CO(3)(2-) oxygen atoms, leading to large fluctuations in the number of first-shell waters, ranging from 1 to 6 (HF/MM) and 2 to 7 (B3LYP/MM), with the prevalent value of 3. Upon comparing the HF and B3LYP methods in describing this hydrated ion, the latter is found to overestimate the hydrogen-bond strength in the CO(3)(2-)-water complexes, resulting in a slightly more compact hydration structure at each of the CO(3)(2-) oxygens.     

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.     

Molecular insights into the binding affinity and specificity of the hemagglutinin cleavage loop from four highly pathogenic H5N1 isolates towards the proprotein convertase furin

Abstract  

The furin (FR) complex with each of four different sequences of hemagglutinin from the highly pathogenic H5N1 strains (HPH5), which were identified during the 2004–2010 influenza outbreaks in Thailand, were evaluated by molecular dynamics simulations, so as to compare the specificity and recognition of the enzyme–substrate binding. Relative to the conventional HPH5 inserted (H5Sq1, RERRRKKR), the S5-R or S6-R arginine residue is replaced by the smaller lysine in the H5Sq2 (RERKRKKR) and H5Sq3 (REKRRKKR) strains, respectively, whereas the S3-K lysine residue is deleted in H5Sq4 (RERRR_KR). The molecular dynamics results of the intermolecular interactions, in terms of hydrogen bonds and per-residue decomposition energy, between the substrate and furin revealed that the deletion of the positively charged amino acid at the S3 position in H5Sq4 leads to a notably weaker binding and specificity with the furin active site compared with that of FR–H5Sq1. A slight change in the substrate binding was found in the FR–H5Sq2 and FR–H5Sq3 complexes as a result of the replacement of the arginine with the shorter side-chained lysine (same positive charge). Altogether, the predicted binding free energy of the enzyme–substrate complexes was found to be in the following order: FR–H5Sq1 < FR–H5Sq2 ~ FR–H5Sq3 ≪ FR–H5Sq4.