No abstract

The feasibility of using cucurbituril host molecule as a probable actinyl cation binders candidate is investigated through density functional theory based calculations. Various possible binding sites of the cucurbit[5]uril host molecule to uranyl are analyzed and based on the binding energy evaluations, μ(5)-binding is predicted to be favored. For this coordination, the structure, vibrational spectra, and binding energies are evaluated for the binding of three actinyls in hexa-valent and penta-valent oxidation states with functionalized cucurbiturils. Functionalizing cucurbituril with methyl and cyclohexyl groups increases the binding affinities of actinyls, whereas fluorination decreases the binding affinities as compared to the native host molecule. Surprisingly hydroxylation of the host molecule does not distinguish the oxidation state of the three actinyls.

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.

Aim: To performe a time-dependent topographical delineation of protein-drug interactions to gain molecular insight into the supremacy of Ortho-7 over HI-6 in reactivating tabun-conjugated mouse acetylcholinesterase (mAChE). Methods: We conducted all-atom steered molecular dynamics simulations of the two protein-drug complexes. Through a host of protein-drug interaction parameters (rupture force profiles, hydrogen bonds, water bridges, hydrophobic interactions), geometrical, and orientation ordering of the drugs, we monitored the enzyme's response during the release of the drugs from its active-site. Results: The results show the preferential binding of the drugs with the enzyme. The pyridinium ring of HI-6 shows excellent complementary binding with the peripheral anionic site, whereas one of two identical pyridinium rings of Ortho-7 has excellent binding compatibility in the enzyme active-site where it can orchestrate the reactivation process. We found that the active pyridinium ring of HI-6 undergoes a complete turn along the active site axis, directed away from the active-site region during the course of the simulation. Conclusion: Due to excellent cooperative binding of Ortho-7, as rendered by several cation-π interactions with the active-site gorge of the enzyme, Ortho-7 may be a more efficient reactivator than HI-6. Our work supports the growing body of evidence that the efficacy of the drugs is due to the differential bindings of the oximes with AChE and can aid to the rational design of oxime drugs.

Achieving an efficient separation of chemically similar Am(3+)/Eu(3+) pair in high level liquid waste treatment is crucial for managing the long-term nuclear waste disposal issues. The use of sophisticated supramolecules in a rigid framework could be the next step toward solving the long-standing problem. Here, we have investigated the possibility of separating Am(3+)/Eu(3+) pair with cucurbit-[5]-uril (CB[5]), a macrocycle from the cucurbit-[n]-uril family, using relativistic density functional theory (DFT) based calculations. We have explored the structures, binding, and energetics of metal-CB[5] complexation processes with and without the presence of counterions. Our study reveals an excellent selectivity of Eu(3+) over Am(3+) with CB[5] (ion exchange free energy, ΔΔGAm/Eu > 10 kcal mol(-1)). Both metals bind with the carbonyl portals via μ(5) coordination arrangement with the further involvement of three external water molecules. The presence of counterions, particularly nitrate, inside the hydrophobic cavity of CB[5], induces a cooperative cation-anion binding, resulting in enhancement of metal binding at the host. The overall binding process is found to be entropy driven resembling the recent experimental observations (Rawat et al. Dalton Trans. 2015, 44, 4246-4258). The optimized structural parameters for Eu(3+)-CB[5] complexes are found to be in excellent agreement with the available experimental information. To rationalize the computed selectivity trend, electronic structures are further scrutinized using energy decomposition analysis (EDA), quantum theory of atom in molecules (QTAIM), Mülliken population analysis (MPA), Nalewajski-Mrojek (NM) bond order, and molecular orbital analyses. Strong electrostatic ion-dipole interaction along with efficient charge transfer between CB[5] and Eu(3+) outweighs the better degree of covalency between CB[5] and Am(3+) leading to superior selectivity of Eu(3+) over Am(3+).

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.

Acid-related gastric diseases are associated with disorder of digestive tract acidification due to the acid secretion by gastric proton pump, H,K-ATPase. Omeprazole is one of the persuasive irreversible inhibitor of the proton pump H,K-ATPase. However, the reports on the mechanistic pathway of irreversible proton pump inhibitors (PPIs) on the acid activation and formation of disulfide complex are scarce in the literature. We have examined the acid activation PPIs, i.e., timoprazole, S-omeprazole and R-omeprazole using M062X/6-31++G(d,p) in aqueous phase with SMD solvation model. The proton pump inhibitor is a prodrug and activated in the acidic canaliculi of the gastric pump H,K-ATPase to sulfenic acid which can either form another acid activate intermediate sulfenamide or a disulfide complex with cysteine amino acid of H,K-ATPase. The quantum chemical calculations suggest that the transition state (TS5) for the disulfide complex formation is the rate-determining step of the multistep acid inhibition process by PPIs. The free energy barrier of TS5 is 5.5 kcal/mol higher for timoprazole compared to the S-omeprazole. The stability of the transition state for the formation of disulfide bond between S-omeprazole and cysteine amino acid of H,K-ATPase is governed by inter- and intramolecular hydrogen bonding. The disulfide complex for S-omeprazole is thermodynamically more stable by 4.5 kcal/mol in aqueous phase compared to disulfide complex of timoprazole, which corroborates the less efficacy of timoprazole as irreversible PPI for acid inhibition process. It has been speculated that sulfenic acid can either form sulfenamide or a stable disulfide complex with cysteine amino acid residue of H,K-ATPase. The M062X/6-31++G(d,p) level of theory calculated results reveal that the formation of tetra cyclic sulfenamide is unfavored by ∼17 kcal/mol for S-omeprazole and 11.5 kcal/mol for timoprazole compared to the disulfide complex formation in each case. The DFT calculations have further shed light on the acid activation process of R- and S-isomers of omeprazole. The calculated results suggest that the efficacy of these isomers lie on their metabolic pathway and excretion from human body.

We introduce a fractional single-file diffusion (FSFD) model describing the stochastic evolution of an assembly of nonpassing particles (single-file (SF) system), where in addition to the excluded-volume interaction between the particles, they individually follow sticky trajectories (fractional diffusion (FD)) on a fractal support. While the long-range correlations between the particles mediated through the former interaction is treated through the reflection principle, the sticky trajectories involve fractional time derivatives whose order is distributed over the interval from zero to one. A detailed study of the propagators of FSFD is made and compared with the normal SF diffusion through their time evolution in state space. Both the diffusion mechanisms, SF and FD, are known for anomalous subdiffusive behaviour; their combined entity, FSFD has been found to exhibit further anomality leading to ultraslow diffusion. Through an integral transformation of the Gaussian propagator, a generalised expression for the mean square displacement of a tagged particle undergoing FSFD is obtained.