Dynamic Mechanism of a Fluorinated Oxime Reactivator Unbinding from AChE Gorge in Polarizable Water.

A well-tempered metadynamics simulation is performed to study the unbinding process of a fluorinated oxime (FHI-6) drug from the active-site gorge of acetylcholinesterase enzyme in a polarizable water medium. Cation-π interactions and water bridge and hydrogen bridge formations between the protein and the drug molecule are found to strongly influence the unbinding process, forming basins and barriers along the gorge pathway. Distinct unbinding pathways are found when FHI-6 was compared with its recently reported nonfluorinated analogue, HI-6. For example, because of permanent positive charges on both the pyridinium rings of HI-6, it exhibits the minimum in the potential of mean force of the unbinding process in the gorge mouth (where the peripheral anion site, PAS, of the enzyme is located), which is largely caused by cation-π interactions. However, the same interaction, both in the catalytic active-site (CAS) and PAS regions, is found to be greatly enhanced in its lipophilic fluorinated analogue, FHI-6, causing a deep potential energy minimum in the bound state. This may render FHI-6 to be held more firmly in the CAS region of the gorge, as is also evidenced from the microkinetics of unbinding transitions, measured through a combination of metadynamics and hyperdynamics simulations.

Water isotope effect on the thermostability of a polio viral RNA hairpin: A metadynamics study.

Oral polio vaccine is considered to be the most thermolabile of all the common childhood vaccines. Despite heavy water (D2O) having been known for a long time to stabilise attenuated viral RNA against thermodegradation, the molecular underpinnings of its mechanism of action are still lacking. Whereas, understanding the basis of D2O action is an important step that might reform the way other thermolabile drugs are stored and could possibly minimize the cold chain problem. Here using a combination of parallel tempering and well-tempered metadynamics simulation in light water (H2O) and in D2O, we have fully described the free energy surface associated with the folding/unfolding of a RNA hairpin containing a non-canonical basepair motif, which is conserved within the 3'-untranslated region of poliovirus-like enteroviruses. Simulations reveal that in heavy water (D2O) there is a considerable increase of the stability of the folded basin as monitored through an intramolecular hydrogen bond (HB), size, shape, and flexibility of RNA structures. This translates into a higher melting temperature in D2O by 41 K when compared with light water (H2O). We have explored the hydration dynamics of the RNA, hydration shell around the RNA surface, and spatial dependence of RNA-solvent collective HB dynamics in the two water systems. Simulation in heavy water clearly showed that D2O strengthens the HB network in the solvent, lengthens inter-residue water-bridge lifetime, and weakens dynamical coupling of the hairpin to its solvation environment, which enhances the rigidity of solvent exposed sites of the native configurations. The results might suggest that like other added osmoprotectants, D2O can act as a thermostabilizer when used as a solvent.

Protein-Drug Interactions with Effective Polarization in Polarizable Water: Oxime Unbinding from AChE Gorge.

Despite the fact that polarizability of water is different in the bulk and in protein, simulations of protein-ligand complexes are mostly carried out in nonpolarizable water media. We present oxime (HI-6) unbinding from the active site gorge of AChE, known to be strongly influenced by intermolecular cation-π, hydrogen bridge (HB) and water bridge (WB) interactions and by molecular simulations with effective polarization in polarizable mean-field model of TIP3P water. Enabled by the recent availability of a method of obtaining microkinetics of rare events, we set out to investigate the rate constants of unbinding transitions from one basin to the other through a combination of metadynamics and hyperdynamics simulations. The results underpin the importance of electronic polarization effects on the pathways, potential of mean force, rate constants, and HB and WB dynamics of unbinding transitions of a drug molecule ligated to protein interior. The method is also applicable to unravel the binding mechanisms.

Effect of a buried ion pair in the hydrophobic core of a protein: An insight from constant pH molecular dynamics study.

Constant pH molecular dynamics (CpHMD) is a commonly used sampling method, which incorporates the coupling of conformational flexibility and protonation state of a protein during the simulation by using pH as an external parameter. The effects on the structure and stability of a hyperstable variant of staphylococcal nuclease (Δ+PHS) protein of an artificial charge pair buried in its hydrophobic core are investigated by applying both CpHMD and accelerated molecular dynamics coupled with constant pH (CpHaMD) methods. Generalized Born electrostatics is used to model the solvent water. Two sets of starting coordinates of V23E/L36K variant of Δ+PHS, namely, Maestro generated coordinates from Δ+PHS and crystal structure coordinates of the same are considered for detail investigations. On the basis of root mean square displacement (RMSD) and root mean square fluctuations (RMSF) calculations, it is observed that this variant is stable over a wide range of pH. The calculated pKa values for aspartate and glutamate residues based on both CpHMD and CpHaMD simulations are consistent with the reported experimental values (within ± 0.5 to ± 1.5 pH unit), which clearly indicates that the local chemical environment of the carboxylic acids in V23E/L36K variant are comparable to the parent form. The strong salt bridge interaction between the mutated pair, E23/K36 and additional hydrogen bonds formed in the V23E/L36K variant, may help to compensate for the unfavorable self-energy experienced by the burial of these residues in the hydrophobic core. However, from RMSD, RMSF, and pKa analysis, no significant change in the global conformation of V23E/L36K variant with respect to the parent form, Δ+PHS is noticed.

Unbinding free energy of acetylcholinesterase bound oxime drugs along the gorge pathway from metadynamics-umbrella sampling investigation.

Because of the pivotal role that the nerve enzyme, acetylcholinesterase plays in terminating nerve impulses at cholinergic synapses. Its active site, located deep inside a 20 Å gorge, is a vulnerable target of the lethal organophosphorus compounds. Potent reactivators of the intoxicated enzyme are nucleophiles, such as bispyridinium oxime that binds to the peripheral anionic site and the active site of the enzyme through suitable cation-π interactions. Atomic scale molecular dynamics and free energy calculations in explicit water are used to study unbinding pathways of two oxime drugs (Ortho-7 and Obidoxime) from the gorge of the enzyme. The role of enzyme-drug cation-π interactions are explored with the metadynamics simulation. The metadynamics discovered potential of mean force (PMF) of the unbinding events is refined by the umbrella sampling (US) corrections. The bidimensional free energy landscape of the metadynamics runs are further subjected to finite temperature string analysis to obtain the transition tube connecting the minima and bottlenecks of the unbinding pathway. The PMF is also obtained from US simulations using the biasing potential constructed from the transition tube and are found to be consistent with the metadynamics-US corrected results. Although experimental structural data clearly shows analogous coordination of the two drugs inside the gorge in the bound state, the PMF of the drug trafficking along the gorge pathway point, within an equilibrium free energy context, to a multistep process that differs from one another. Routes, milestones and subtlety toward the unbinding pathway of the two oximes at finite temperature are identified.

Structure and photoelectron spectral properties of I(-)·nCO2 clusters: an ab initio study.

Structures and photoelectron spectral properties of I(-)·nCO(2) (n=1-7) clusters are presented at the level of second-order Møller-Plesset perturbation theory with relativistic corrections. Triple split-valence 6-311++G(d,p) basis set functions are employed herein. It is observed that the CO(2) molecules approach the I(-) anion from one side in all the clusters and that I(-)·nCO(2) clusters prefer the surface structure. The calculated vertical detachment energy of these clusters is in excellent agreement with the reported experimentally measured values (within 4%). Efforts are also made to extract vertical detachment energy of large size of clusters, including the bulk. The extracted vertical detachment energy values for larger clusters (n=8-13) by employing the microscopic theory-based expression are also close (within 4%) to that of the experimentally measured values.

Theoretical study on the spectroscopic properties of CO3(*-).nH2O clusters: extrapolation to bulk.

Vertical detachment energies (VDE) and UV/Vis absorption spectra of hydrated carbonate radical anion clusters, CO(3)(*-).nH(2)O (n=1-8), are determined by means of ab initio electronic structure theory. The VDE values of the hydrated clusters are calculated with second-order Moller-Plesset perturbation (MP2) and coupled cluster theory using the 6-311++G(d,p) set of basis functions. The bulk VDE value of an aqueous carbonate radical anion solution is predicted to be 10.6 eV from the calculated weighted average VDE values of the CO(3)(*-).nH(2)O clusters. UV/Vis absorption spectra of the hydrated clusters are calculated by means of time-dependent density functional theory using the Becke three-parameter nonlocal exchange and the Lee-Yang-Parr nonlocal correlation functional (B3LYP). The simulated UV/Vis spectrum of the CO(3)(*-).8H(2)O cluster is in excellent agreement with the reported experimental spectrum for CO(3)(*-) (aq), obtained based on pulse radiolysis experiments.