Arup Kumar Pathak

A comparative ab initio study of Br2•− and Br2 water clusters

The work presents ab initio results on structure and electronic properties of Br2*-.nH2O(n=1-10) and Br2.nH2O(n=1-8) hydrated clusters to study the effects of an excess electron on the microhydration of the halide dimer. A nonlocal density functional, namely, Beckes half-and-half hybrid exchange-correlation functional is found to perform well on the present systems with a split valence 6-31++G(d,p) basis function. Geometry optimizations for all the clusters are carried out with several initial guess structures and without imposing any symmetry restriction. Br2*-.nH2O clusters prefer to have symmetrical double hydrogen-bonding structures. Results on Br2.nH2O(n>or=2) cluster show that the O atom of one H2O is oriented towards one Br atom and the H atom of another H2O is directed to other Br atom making Br2 to exist as Br+-Br- entity in the cluster. The binding and solvation energies are calculated for the Br2*-.nH2O and Br2.nH2O clusters. Calculations of the vibrational frequencies show that the formation of Br2*- and Br2 water clusters induces significant shifts from the normal modes of isolated water. Excited-state calculations are carried out on Br2*-.nH2O clusters following configuration interaction with single electron excitation procedure and UV-VIS absorption profiles are simulated. There is an excellent agreement between the present theoretical UV-VIS spectra of Br2*-.10H2O cluster and the reported transient optical spectra for Br2*- in aqueous solution.

Microhydration shell structure in Cl2∙−∙nH2O clusters: A theoretical study

We present the results of a detailed study on structure and electronic properties of hydrated cluster Cl2*-.nH2O (n = 1-7) based on a nonlocal density functional, namely, Beckes [J. Chem. Phys. 98, 1372 (1993)] half and half hybrid exchange-correlation functional with a split valence 6-311++G(d,p) basis function. Geometry optimizations for all the clusters are carried out with various possible initial guess structures without any symmetry restriction. Several minimum energy structures (conformers) are predicted with a small difference in total energy. There is a competition between the binding of solvent H2O units with Cl2*- dimer radical anion directly through ion-molecule interaction and forming interwater hydrogen-bonding network in Cl2*-.nH2O (n > or = 2) hydrated cluster. Structure having interwater H-bonded network is more stable over the structure where H2O units are connected to the solute dimer radical anion Cl2*- rather independently either by single or double H bonding in a particular size (n) of hydrated cluster Cl2*-.nH2O. At the maximum four solvent H2O units reside in interwater H-bonding network present in these hydrated clusters. It is observed that up to six H2O units are independently linked to the anion having four double H bondings and two single H bondings suggesting the primary hydration number of Cl2*- to be 6. In all these clusters, the odd electron is found to be mostly localized over the two Cl atoms and these two atoms are bound by a three-electron hemibond. Calculated interaction (between solute and different water clusters) and vertical detachment energy profiles show saturation at n = 6 in the hydrated cluster Cl2*-.nH2O (n = 1-7). However, calculated solvation energy increases with the increase in number of solvent H2O molecules in the cluster. Interaction energy varies linearly with vertical detachment energy for the hydrated clusters Cl2*-.nH2O (n < or = 6). Calculation of the vibration frequencies show that the formation of Cl2*(-)-water clusters induces significant shifts from the normal stretching modes of isolated water. A clear difference in the pattern of IR spectra is observed in the O-H stretching region of water from hexa- to heptahydrated cluster.

Structure, energy, and IR spectra of I2- . nH2O clusters (n=1-8) : A theoretical study

The authors report theoretical results on structure, bonding, energy, and infrared spectra of iodine dimer radical anion hydrated clusters, I(2) (-).nH(2)O (n=1-8), based on a systematic study following density functional theory. Several initial guess structures are considered for each size cluster to locate minimum energy conformers with a Gaussian 6-311++G(d,p) split valence basis function (triple split valence 6-311 basis set is applied for iodine). It is observed that three different types of hydrogen bonded structures, namely, symmetrical double hydrogen bonding, single hydrogen bonding, and interwater hydrogen bonding structures, are possible in these hydrated clusters. But conformers having interwater hydrogen bonding arrangements are more stable compared to those of double or single hydrogen bonded structures. It is also noticed that up to four solvent H(2)O units can reside around the solute in interwater hydrogen bonding network. At the maximum six H(2)O units are independently linked to the dimer anion having four double hydrogen bonding and two single hydrogen bonding, suggesting the hydration number of I(2) (-) to be 6. However, conformers having H(2)O units independently linked to the iodine dimer anion are not the most stable structures. In all these hydrated clusters, the odd electron is found to be localized over two I atoms and the two atoms are bound by a three-electron hemi bond. The solvation, interaction, and vertical detachment energies are calculated for all I(2) (-).nH(2)O clusters. Energy of interaction and vertical detachment energy profiles show stepwise saturation, indicating geometrical shell closing in the hydrated clusters, but solvation energy profile fails to show such behavior. A linear correlation is observed between the calculated energy of interaction and vertical detachment energy. It is observed that formation of I(2) (-)-water cluster induces significant shifts from the normal O-H stretching modes of isolated H(2)O. However, bending mode of H(2)O remains insensitive to the successive addition of solvent H(2)O units. Weighted average energy profiles and IR spectra are reported for all the hydrated clusters based on the statistical population of individual conformers at room temperature.

Global Minimum‐Energy Structure and Spectroscopic Properties of I2.−⋅n H2O Clusters: A Monte Carlo Simulated Annealing Study

The vibrational (IR and Raman) and photoelectron spectral properties of hydrated iodine-dimer radical-anion clusters, I(2)(*-) x n H(2)O (n=1-10), are presented. Several initial guess structures are considered for each size of cluster to locate the global minimum-energy structure by applying a Monte Carlo simulated annealing procedure including spin-orbit interaction. In the Raman spectrum, hydration reduces the intensity of the I-I stretching band but enhances the intensity of the O-H stretching band of water. Raman spectra of more highly hydrated clusters appear to be simpler than the corresponding IR spectra. Vibrational bands due to simultaneous stretching vibrations of O-H bonds in a cyclic water network are observed for I(2)(*-) x n H(2)O clusters with n > or = 3. The vertical detachment energy (VDE) profile shows stepwise saturation that indicates closing of the geometrical shell in the hydrated clusters on addition of every four water molecules. The calculated VDE of finite-size small hydrated clusters is extrapolated to evaluate the bulk VDE value of I(2)(*-) in aqueous solution as 7.6 eV at the CCSD(T) level of theory. Structure and spectroscopic properties of these hydrated clusters are compared with those of hydrated clusters of Cl(2)(*-) and Br(2)(*-).

Theoretical studies on photoelectron and IR spectral properties of Br2∙−(H2O)n clusters

We report vertical detachment energy (VDE) and IR spectra of Br2.-.(H2O)n clusters (n=1-8) based on first principles electronic structure calculations. Cluster structures and IR spectra are calculated at Beckes half-and-half hybrid exchange-correlation functional (BHHLYP) with a triple split valence basis function, 6-311++G(d,p). VDE for the hydrated clusters is calculated based on second order Moller-Plesset perturbation (MP2) theory with the same set of basis function. On full geometry optimization, it is observed that conformers having interwater hydrogen bonding among solvent water molecules are more stable than the structures having double or single hydrogen bonded structures between the anionic solute, Br2.-, and solvent water molecules. Moreover, a conformer having cyclic interwater hydrogen bonded network is predicted to be more stable for each size hydrated cluster. It is also noticed that up to four solvent H2O units can reside around the solute in a cyclic interwater hydrogen bonded network. The excess electron in these hydrated clusters is localized over the solute atoms. Weighted average VDE is calculated for each size (n) cluster based on statistical population of the conformers at 150 K. A linear relationship is obtained for VDE versus (n+3)(-1/3) and bulk VDE of Br2.- aqueous solution is calculated as 10.01 eV at MP2 level of theory. BHHLYP density functional is seen to make a systematic overestimation in VDE values by approximately 0.5 eV compared to MP2 data in all the hydrated clusters. It is observed that hydration increases VDE of bromine dimer anion system by approximately 6.4 eV. Calculated IR spectra show that the formation of Br2.--water clusters induces large shifts from the normal O-H stretching bands of isolated water keeping bending modes rather insensitive. Hydrated clusters, Br2.-.(H2O)n, show characteristic sharp features of O-H stretching bands of water in the small size clusters.

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. Proteins 2014; 82:1799–1818.

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.

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.

A theoretical study on structures, energetics, and spectra of Br–.nCO2 clusters: Towards bridging the gap between micro-domain and macro-domain

Structures, energetics, and spectra of Br(-).nCO(2) (n = 1-8) clusters are studied based on ab initio electronic structure theory. The geometry of each size of clusters is evaluated by employing second-order Moller-Plesset (MP2) perturbation theory. It is observed that the solvent CO(2) molecules approach the bromide moiety from one side in an asymmetric fashion except for the Br(-).8CO(2) cluster. Simple electrostatic model for charge-quadrupole interactions is valid for the Br(-).nCO(2) clusters. Reduced variational space based energy decomposition method shows that the electrostatic interaction is the major component and polarization and charge transfer energies are the other significant components of the total interaction energy. Both adiabatic and vertical electron detachment energies and solvation energies are calculated at MP2 level of theory. We have observed an excellent agreement between theory and experiment for the vertical detachment and solvation energies. Calculated quantities based on the analytical expression which connects the finite domain to macroscopic one are found to be very good in agreement with the available experimental results. The present study reveals a 2.6 eV increase in the detachment energy of bromide anion due to the solvation effect of CO(2), which is relatively small compared to that of the corresponding 4.7 eV increase in detachment energy in water.