Selectivity of a Singly Permeating Ion in Nonselective NaK Channel: Combined QM and MD Based Investigations

Achieving highly efficient and selective cesium extraction using 1,3-di-octyloxycalix[4]arene-crown-6 in n-octanol based solvent system: experimental and DFT investigation

Due to the long half-life of ¹³⁷Cs (t_1/2 ∼ 30 years), the selective extraction of cesium (Cs) from high level liquid waste is of paramount importance in the back end of the nuclear fuel cycle to avoid long term surveillance of high radiotoxic waste. As 1,3-di-octyloxycalix[4]arene-crown-6 (CC6) is suggested to be a promising candidate for selective Cs extraction, the improvement in the Cs extraction efficiency by CC6 has been investigated through the optimization of the effect of dielectric media on the extraction process. The effects of the feed acid (HNO₃, HCl, and HClO₄) and the composition of the diluents for the ligand in the organic phase on the extraction efficiency of Cs have been investigated systematically. In 100% n-octanol medium, Cs is found to form a 1 : 1 ion-pair complex with CC6 (0.03 M) providing a very high distribution ratio of D_Cs ∼ 22, suggesting n-octanol as the most suitable diluent for Cs extraction. No significant interference of other relevant cations such as Na, Mg and Sr was observed on the D_Cs value in the optimized solvent system. Density functional theory (DFT) based calculations have been carried out to elucidate the reason of ionic selectivity and enhanced Cs extraction efficiency of CC6 in the studied diluent systems. In addition to the ionic size-based selectivity of the crown-6 cavity, the polarity of the organic solvent system, the hydration energy of the ion, and the relative reorganization of CC6 upon complexation with Cs are understood to have roles in achieving the enhanced efficiency for the extraction of Cs by the CC6 extractant in nitrobenzene medium.

Separation scheme was developed for selective extraction of long-lived fission product ¹³⁷Cs using substituted calix crown 6 ether from aqueous acidic solution.

Revealing the cholinergic inhibition mechanism of Alzheimer's by galantamine: a metadynamics simulation study

Galantamine is one of the approved drugs based on the cholinergic hypothesis for the symptomatic treatment of mild to moderate Alzheimer's disease (AD). The etiology of AD is not fully known; however, the reported cholinergic hypothesis suggests the inadequate synthesis of the neurotransmitter acetylcholine (ACh) is responsible for this disease. The crystal structure of galantamine bound human acetylcholinesterase (hAChE) has been reported; however, the inhibition mechanism of hAChE by galantamine is not well understood. A Well-tempered metadynamics (WTMtD) simulation study has been performed with the crystal structure of galantamine bound hAChE. The reported mechanism for the degradation of ACh is suggested through a proton transfer process from a carboxylic group of Glu334 to the hydroxyl group of Ser203, which attacks ACh for the degradation to acetic acid and choline. Such proton transfer process is lowered in the presence of galantamine due to the separation of catalytic triad inside the gorge of AChE as observed with WTMtD. A docking study has been performed to examine the ACh's binding with the catalytic triad of galantamine bound hAChE. The docking results reveal that the approach of ACh to the catalytic triad is interrupted due to the galantamine's presence in the gorge of the enzyme.

MD simulation reveals differential binding of Cm(III) and Th(IV) with serum transferrin at acidic pH

The iron carrier human serum transferrin (sTf) is known to transport other metals, including some actinides (An). Radiotoxic An are routinely involved in the nuclear fuel cycle and the possibility of their accidental exposure cannot be ruled out. Understanding An interaction with sTf assumes a greater significance for the development of safe and efficacious chelators for their removal from the blood stream. Here we report several 100 ns equilibrium MD simulations of Cm(III)- and Th(IV)-loaded sTf at various protonation states of the protein to explore the possibility of the two An ions release and speciation. The results demonstrate variation in protonation state of dilysine pair (K206 and K296) and the tyrosine (Y188) residue is necessary for the opening of Cm(III)-bound protein and the release of the ion. For the tetravalent thorium, protonation of dilysine pair suffices to cause conformational changes of protein. However, in none of the protonation states, Th(IV) releases from sTf because of its strong electrostatic interaction with D63 in the first shell of the sTf binding cleft. Analysis of hydrogen bond, water bridge, and the evaluation of potential of mean forces of the An ions' release from sTf, substantiate the differential behavior of Cm(III) and Th(IV) at endosomal pH. The results provide insight in the regulation of Cm(III) and Th(IV) bioavailability that may prove useful for effective design of their decorporating agents and as well may help the future design of radiotherapy based on tetravalent ions.

Molecular dynamics simulations of plutonium binding and its decorporation from the binding-cleft of human serum transferrin

The possibility of plutonium (Pu) intake by radiation workers can not be ruled out. Transportation of Pu(IV) to various organs/cells is mainly carried through iron-carrying protein, serum transferrin (sTf), by receptor-mediated endocytosis. Understanding the Pu-sTf interaction is a primary step toward future design of its decorporating agents. We report MD simulations of Pu(IV) binding with sTf and look out for its decorporation at extracellular pH using suitable ligands. MD simulations were carried out in polarizable water environment at different protonation states of the protein. Results unravel the binding motif of Pu(IV): (1) sTf binds the ion in closed conformation at extracellular serum pH with carbonate as synergistic anions, (2) change in protonation state of dilysine (K206 and K296)-trigger and that of the carbonate ion at acidic endosomal pH is found to cause conformational changes of protein, conducive for the heavy ion to be released, although; (3) strong electrostatic interaction between D63 in the binding-cleft and Pu(IV) is found not to ever set free the ion. In an endeavour to decorporate Pu(IV), fragmented molecular form of hydroxypyridinone (HOPO) and catechol (CAM)-based ligands are docked at the binding site (BS) of the protein and metadynamics simulations are conducted. Pu(IV) binding at BS is found to be so strong that it was not detached from BS with the docked HOPO. However, for the identical set of simulation parameters, CAM is found to facilitate dislodging the heavy ion from the protein's binding influence. Differential behaviour of the two chelators is further explored. Fragmented molecular form of hydroxy-pyridinone (HOPO) and catecholamide (CAM) ligands were docked at the binding-site (BS) of human serum transferrin (sTf) to explore their feasibility as plausible Pu(IV) decorporating agents by employing metadynamics method. CAM was found to dislodge Pu from the sTf BS, while HOPO could not.

Revealing the mechanistic pathway of cholinergic inhibition of Alzheimer's disease by donepezil: a metadynamics simulation study

Donepezil, an acetylcholinesterase inhibitor, is an approved drug for the symptomatic treatment of Alzheimer's disease (AD). The mechanistic pathway for the inhibition mechanism of acetylcholinesterase (AChE) by donepezil is not well explored. We report for the first time, the inhibition mechanism of AChE by the donepezil drug molecule for the hydrolysis of acetylcholine (ACh) with docking and well-tempered metadynamics (WTMtD) simulations with a human acetylcholinesterase (hAChE) crystal structure (). This study explored the orientation of the donepezil drug molecule inside the gorge of AChE. The 1D free energy surface obtained from WTMtD simulation studies reveals that the orientation of donepezil in the crystal donepezil (-87.25 kJ mol^-1) is energetically more favored than the other orientation of donepezil (-74.74 kJ mol^-1) for inhibition of AChE. The free energy landscape computation for the two sets of CVs further corroborates the 1D free energy surface. The WTMtD simulation performed with the crystal structure of donepezil bound hAChE gives the conformation of donepezil at Basin-I as similar to the conformation of donepezil observed in the crystal structure (). The WTMtD simulations further reveal that the bridged water molecules are more ordered near the catalytic triad of AChE to deter the nucleophilicity of Ser203 through intermolecular hydrogen bonding when donepezil approaches near to the active site gorge of AChE. The presence of donepezil near the active site of AChE can inhibit its approach for ACh hydrolysis; this is revealed through the docking study, where the drug molecule inside the active gorge of hAChE restricts the approach of ACh to Ser203 for the hydrolysis process.

Binding of Cm(III) and Th(IV) with Human Transferrin at Serum pH: Combined QM and MD Investigations

Human serum transferrin (sTf) can also function as a noniron metal transporter since only 30% of it is typically saturated with a ferric ion. While this function of sTf can be fruitfully utilized for targeted delivery of certain metal therapeutics, it also runs the risk of trafficking the lethal radionuclides into cells. A large number of actinide (An) ions are known to bind to the iron sites of sTf although molecular-level understanding of their binding is unclear. Understanding the radionuclide interaction with sTf is a primary step toward future design of their decorporating agents since irrespective of the means of contamination, the radionuclides are absorbed and transported by blood before depositing into target organs. Here, we report an extensive multiscale modeling approach of two An (curium(III) and thorium(IV)) ions' binding with sTf at serum physiological pH. We find that sTf binds both the heavy ions in a closed conformation with carbonate as synergistic anions and the An-loaded sTf maintains its closed conformation even after 100 ns of equilibrium molecular dynamics (MD) simulations. MD simulations are performed in a polarizable water environment, which also incorporates electronic continuum corrections for ions via charge rescaling. The molecular details of the An coordination and An exchange free energies with iron in the interdomain cleft of the protein are evaluated through a combination of quantum mechanical (QM) and MD studies. In line with reported experimental observations, well-tempered metadynamics results of the ions' binding energetics show that An-sTf complexes are less stable than Fe-sTf. Additionally, curium(III) is found to bind more weakly than thorium(IV). The latter result might suggest relative attenuation of thorium(IV) cytotoxicity when compared with curium(III).

H. R, Thamleena A. H, Joseph A, Ben A, V. U. K. Roles of different amino-acid residues towards binding and selective transport of K+ through KcsA K+-ion channel

Each amino acid in the selectivity filter plays a distinct role towards binding and transport of K⁺ ion through KcsA.

Designed inhibitors with hetero linkers for gastric proton pump H+,K+-ATPase: Steered molecular dynamics and metadynamics studies

Acid suppressant SCH28080 and its derivatives reversibly reduce acid secretion activity of the H⁺,K⁺-ATPase in a K⁺ competitive manner. The results on homologation of the SCH28080 by varying the linker chain length suggested the improvement in efficacy. However, the pharmacokinetic studies reveal that the hydrophobic nature of the CH₂ linker units may not help it to function as a better acid suppressant. We have exploited the role of linker unit to enhance the efficacy of such reversible acid suppressant drug molecules using hetero linker, i.e., disulfide and peroxy linkers. The logarithm of partition coefficient defined for a drug molecule relates to the partition coefficient, which allows the optimum solubility characteristics to reach the active site. The logarithm of partition coefficient calculated for the designed inhibitors suggests that inhibitors would possibly reach the active site in sufficient concentration like in the case of SCH28080. The steered molecular dynamics studies have revealed that the Inhibitor-1 with disulfide linker unit is more stable at the active site due to greater noncovalent interactions compared to the SCH28080. Centre of mass distance analysis suggests that the Cysteine-813 amino acid residue selectively plays an important role in the inhibition of H⁺,K⁺-ATPase for Inhibitor-1. Furthermore, the quantum chemical calculations with M11L/6-31+G(d,p) level of theory have been performed to account the noncovalent interactions responsible for the stabilization of inhibitor molecules in the active site gorge of the gastric proton pump at different time scale. The hydrogen bonding and hydrophobic interaction studies corroborate the center of mass distance analysis as well. Well-tempered metadynamics free energy surface and center of mass separation analysis for the Inhibitor-1 is in good agreement with the steered molecular dynamics results. The torsional angle of the linker units seems to be crucial for better efficacy of drug molecules. The torsional angle of linker units of SCH28080 (COCH₂C) and of Inhibitor 1 (CSSC) prefers to lie within ∼60°-90° for a longer time during the simulations, whereas, the peroxy linker (COOC) of Inhibitor 2 prefers to adopt ∼120-160°. Therefore, it appears that the smaller torsion angle of linker units can achieve better interactions with the active site residues of H⁺,K⁺-ATPase to inhibit the acid secretion activity. The reversible drug molecules with disulfide linker unit would be a promising candidate as proton pump antagonist to H⁺,K⁺-ATPase.

Divalent ions are potential permeating blockers of the non-selective NaK ion channel: combined QM and MD based investigations

The interplay between the hydration energy barrier and optimum site-selectivity controls the permeation and blocking features of divalent metal ions.

Asn47 and Phe114 modulate the inner sphere reorganization energies of type zero copper proteins

The geometric structures and electron transfer properties of type 1 Cu proteins are reasonably understood at the molecular level (E. I. Solomon and R. G. Hadt, Coord. Chem. Rev., 2011, 255, 774-789, J. J. Warren, K. M. Lancaster, J. H. Richards and H. B. Gray, J. Inorg. Biochem., 2012, 115, 119-126). Much understanding of type 1 copper electron transfer reactivity has come from site directed mutagenesis studies. For example, artificial "type zero" Cu-centres constructed in cupredoxin-azurin have showcased the capacity of outer-sphere hydrogen bonding networks to enhance Cu II/I electron transfer reactivity. In this paper, we have elaborated on earlier kinetics and electronic structural studies of type zero Cu by calculating the inner sphere reorganization energies of type 1, type 2, and type zero Cu proteins using density functional theory (DFT). Although the choice of density functionals for copper systems is not straightforward, we have benchmarked the density functionals against the recently reported ESI-PES data for two synthetic copper models (S. Niu, D.-L. Huang, P. D. Dau, H.-T. Liu, L.-S. Wang and T. J. Ichiye, Chem. Theory Comput., 2014, 10, 1283). For the Cu proteins, our calculations predict that changes in the coordination number upon metal reduction lead to large inner sphere reorganization energies for type 2 Cu sites, whereas retention in the coordination number is observed for type zero Cu sites. These variations in the coordination number are modulated by the outer-sphere coordinating residues Asn47 and Phe114, which are involved in hydrogen bonding with the Asp112 side chain.

Computational studies of transport in ion channels using metadynamics

Molecular dynamics simulations have played a fundamental role in numerous fields of science by providing insights into the structure and dynamics of complex systems at the atomistic level. However, exhaustive sampling by standard molecular dynamics is in most cases computationally prohibitive, and the time scales accessible remain significantly shorter than many biological processes of interest. In particular, in the study of ion channels, realistic models to describe permeation and gating require accounting for large numbers of particles and accurate interaction potentials, which severely limits the length of the simulations. To overcome such limitations, several advanced methods have been proposed among which is metadynamics. In this algorithm, an external bias potential to accelerate sampling along selected collective variables is introduced. This bias potential discourages visiting regions of the configurational space already explored. In addition, the bias potential provides an estimate of the free energy as a function of the collective variables chosen once the simulation has converged. In this review, recent contributions of metadynamics to the field of ion channels are discussed, including how metadynamics has been used to search for transition states, predict permeation pathways, treat conformational flexibility that underlies the coupling between gating and permeation, or compute free energy of permeation profiles. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.