Unbinding of hACE2 and inhibitors from the receptor binding domain of SARS-CoV-2 spike protein.

The first direful biomolecular event leading to COVID-19 disease is the SARS-CoV-2 virus surface spike (S) protein-mediated interaction with the human transmembrane protein, angiotensin-converting enzyme 2 (hACE2). Prevention of this interaction presents an attractive alternative to thwart SARS-CoV-2 replications. The development of monoclonal antibodies (mAbs) in the convalescent plasma treatment, nanobody, and designer peptides, which recognizes epitopes that overlap with hACE2 binding sites in the receptor-binding domain (RBD) of S protein (S/RBD) and thereby blocking the infection has been the center stage of therapeutic research. Here we report atomistic and reliable in silico structure-energetic features of the S/RBD interactions with hACE2 and its two inhibitors (convalescent mAb, B38, and an alpaca nanobody, Ty1). The discovered potential of mean forces exhibits free energy basin and barriers along the interaction pathways, providing sufficient molecular insights to design a B38 mutant and a Ty1-based peptide with higher binding capacity. While the mutated B38 forms a 60-fold deeper free energy minimum, the designer peptide (Ty1-based) constitutes 38 amino acids and is found to form a 100-fold deeper free energy minimum in the first binding basin than their wild-type variants in complex with S/RBD. Our strategy may help to design more efficacious biologics towards therapeutic intervention against the current raging pandemic.Communicated by Ramaswamy H. Sarma.

Binding of human serum albumin with uranyl ion at various pH: an all atom molecular dynamics study.

Uranium is routinely handled in various stages of nuclear fuel cycle and its association with human serum albumin (HSA) has been reported in literature, however, their binding characteristics still remains obscure. The present study aims to understand interaction of uranium with HSA by employing all atom molecular dynamics simulation of the HSA-metal ion complex. His67, His247 and Asp249 residues constitute the major binding site of HSA, which capture the uranyl ion (UO22+). A total of six sets of initial coordinates are used for Zn2+-HSA and UO22+-HSA system at pH = 4, 7.4 and 9, respectively. Enhance sampling method, namely, well-tempered meta-dynamics (WT-MtD) is employed to study the binding and un-binding processes of UO22+ and Zn2+ ions. Potential of mean force (PMF) profiles are generated for all the six sets of complexes from the converged WT-MtD run. Various basins and barriers are observed along the (un)binding pathways. Hydrogen bond dynamics and short-range Coulomb interactions are evaluated from the equilibrium run at each basins and barriers for both the ions at all pH values. The binding of UO22+ ion with HSA is the result of the dynamical balance between UO22+-HSA and UO22+-water short range Coulomb interactions. Zn2+ ion interact more strongly than UO22+ at all pH through short range Coulomb interactions. PMF values further concludes that UO22+ cannot associate to the Zn2+ bound HSA protein but can be captured by free HSA at all pH values i.e. endosomal, alkaline and physiological pH.Communicated by Ramaswamy H. Sarma.

Heat-induced transitions of an empty minute virus of mice capsid in explicit water: all-atom MD simulation.

The capsid-like structure of the virus-based protein nanoparticles (NPs) can serve as bionanomaterials, with applications in biomedicines and nanotechnology. Release of packaged material from these nanocontainers is associated with subtle conformational changes of the NP structure, which in vitro, is readily accomplished by heating. Characterizing the structural changes as a function of temperature may provide fresh insights into nanomaterial/antiviral strategies. Here, we have calculated heat induced changes in the properties of an empty minute virus of mice particle using large-scale ≈ 3.0 × 106 all-atom molecular dynamics simulations. We focus on two heat induced structural changes of the NP, namely, dynamical transition (DT) and breathing transition (BT), both characterized by sudden and sharp change of measured parameters at temperatures, TDT and TBT, respectively. While DT is assessed by mean-square fluctuation of hydrogen atoms of the NP, BT is monitored through internal volume and permeation rate of water molecules through the NP. Both the transitions, resulting primarily from collective atomistic motion, are found to occur at temperatures widely separated from one another (TBT>TDT). The breathing motions, responsible for the translocation events of the packaged materials through the NP to kick off, are further probed by computing atomic resolution stresses from NVE simulations. Distribution of equilibrium atomistic stresses on the NP reveals a largely asymmetric nature and suggests structural breathing may actually represent large dynamic changes in the hotspot regions, far from the NP pores, which is in remarkable resemblance with recently conducted hydrogen-deuterium exchange coupled to mass spectrometry experiment. Communicated by Ramaswamy H. Sarma.

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.

Quantum chemical and well-tempered metadynamics study to design adenine analogs for orthogonal Preq1 riboswitch.

Communicated by Ramaswamy H. Sarma.

Temperature Induced Dynamical Transition of Biomolecules in Polarizable and Nonpolarizable TIP3P Water.

Temperature induced dynamical transition (DT), associated with a sharp rise in molecular flexibility, is well-known to be exhibited between 270 and 280 K in glycerol to 200-230 K in hydrated biomolecules and is controlled by diffusivity (viscosity) of the solvation layer. In the molecular dynamics (MD) community, especially for water as a solvent, this has been an intense area of research despite decades of investigations. However, in general, water in these studies is described by empirical nonpolarizable force fields in which electronic polarizability is treated implicitly with effective charges and related parameters. This might have led to the present trait of discovery that DTs of biomolecules, irrespective of the potential functions for water models used, occur within a narrow band of temperature variation (30-40 K). Whereas a water molecule in a biomolecular surface and one in bulk are polarized differently, therefore explicit treatment of water polarizability would be a powerful approach toward the treatment of hydration water, believed to cause the DT manifestation. Using MD simulations, we investigated the effects of polarizable water on the DT of biomolecules and the dynamic properties of hydration water. We chose two types of solutes: globular protein (lysozyme) and more open and flexible RNAs (a hairpin and a riboswitch) with different natures of hydrophilic sites than proteins in general. We found that the characteristic temperature of DT ( TDT) for the solutes in polarizable water is always higher than that in its nonpolarizable counterpart. In particular, for RNAs, the variations are found to be ∼45 K between the two water models, whereas for the more compact lysozyme, it is only ∼4 K. The results are discussed in light of the enormous increase in relaxation times of a liquid upon cooling in the paradigm of dynamic switchover in hydration water with liquid-liquid phase transition, derived from the existence of the second critical point. Our result supports the idea that structures of biomolecules and their interactions with the hydration water determines TDT and provides evidence for the decisive role of polarizable water on the onset of DT, which has been hitherto ignored.

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).

Stereoselective Metabolism of Omeprazole by Cytochrome P450 2C19 and 3A4: Mechanistic Insights from DFT Study.

The efficacy of S-omeprazole as a proton pump inhibitor compared with that of its enantiomer R-omeprazole is studied using density functional theoretical calculations. The pharmacokinetic studies suggest that the efficacy of S-omeprazole presumably depends on metabolic pathway and excretion from the human body. The density functional theory calculations at SMDwater-B3LYP-D3/6-311+G(d,p)/LANL2DZ//B3LYP/6-31G(d)/LANL2DZ with triradicaloid model active species, [Por•+FeIV(SH)O], of CYP2C19 enzyme with high-spin quartet and low-spin doublet states demonstrate C-H bond activation mechanism through a two-state rebound process for the hydroxylation of R-omeprazole and S-omeprazole. The calculated activation free energy barriers for the hydrogen abstraction are 15.7 and 17.5 kcal/mol for R-omeprazole and S-omeprazole, respectively. The hydroxylation of R-omeprazole and S-omeprazole is thermodynamically favored; however, the hydroxylated intermediate of S-omeprazole further disintegrates to metabolite 5- O-desmethylomeprazole with a higher kinetic barrier. We have examined the sulfoxidation of S-omeprazole to omeprazole sulfone metabolite by CYP3A4, and the observed activation free energy barrier is 9.9 kcal/mol. The computational results reveal that CYP2C19 exclusively metabolizes R-omeprazole to hydroxyomeprazole, which is hydrophilic and can easily excrete, whereas CYP3A4 metabolizes S-omeprazole to lipophilic sulfone; hence, the excretion of this metabolite would be relatively slower from the body. The spin density analysis and molecular orbital analysis performed using biorthogonalization calculations indicate that R-omeprazole favors high-spin pathway for metabolism process whereas S-omeprazole prefers the low-spin pathway.

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.

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 CH2 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 (COCH2C) 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 bacterial NaK ion channel is distinctly different from other known ion channels due to its inherent non-selective feature. One of the unexplored and rather interesting features is its ability to permeate divalent metal ions (such as Ca2+ and Ba2+) and not monovalent alkali metal ions. Several intriguing questions about the energetics and structural aspects still remain unanswered. For instance, what causes Ca2+ to permeate as well as block the selectivity filter (SF) of the NaK ion channel and act as a "permeating blocker"? How and at what energetic cost does another chemical congener, Sr2+, as well as Ba2+, a potent blocker of the K+ ion channel, permeate through the SF of the NaK ion channel? Finally, how do their translocation energetics differ from those of monovalent ions such as K+? Here, in an attempt to address these outstanding issues, we elucidate the structure, binding and selectivity of divalent ions (Ca2+, Sr2+ and Ba2+) as they permeate through the SF of the NaK ion channel using all-atom molecular dynamics simulations and density functional theory based calculations. We unveil mechanistic insight into this translocation event using well-tempered metadynamics simulations in a polarizable environment using the mean-field model of water and incorporating electronic continuum corrections for ions via charge rescaling. The results show that, akin to K+ coordination, Sr2+ and Ba2+ bind at the SF in a very similar fashion and remain octa-coordinated at all sites. Interestingly, differing from its local hydration structure, Ca2+ interacts with eight carbonyls to remain at the middle of the S3 site. Furthermore, the binding of divalent metals at SF binding sites is more favorable than the binding of K+. However, their permeation through the extracellular entrance faces a considerably higher energetic barrier compared to that for K+, which eventually manifests their inherent blocking feature.

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.

Unbinding of fluorinated oxime drug from the AChE gorge in polarizable water: a well-tempered metadynamics study.

Despite the fact that fluorination makes a drug more lipophilic, the molecular level understanding of protein-fluorinated drug interactions is very poor. Due to their enhanced ability to penetrate the blood brain barrier, they are suitable for reactivation of organophosphorus inactivated acetylcholinesterase (AChE) in the central nervous system. We systematically studied the unbinding of fluorinated obidoxime (FOBI) and non-fluorinated obidoxime (OBI) from the active site gorge of the serine hydrolase AChE in mean field polarizable water by employing all atom molecular dynamics simulations. It is observed that the unbinding process is strongly influenced by cation-π, hydrogen bond (HB) and water bridge interactions. The FOBI drug interacts more strongly with the protein residues than OBI and this is also verified from quantum mechanical calculations. Distinct unbinding pathways for FOBI and OBI are observed as evident from the 1D and 2D potential of mean force of the unbinding profiles. The present study suggests that the FOBI drug is held more firmly in the gorge of AChE in comparison to OBI and may lead to higher reactivation efficiency of the inactivated enzyme.

Revealing the Mechanistic Pathway of Acid Activation of Proton Pump Inhibitors To Inhibit the Gastric Proton Pump: A DFT Study.

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.

Efficient Separation of Europium Over Americium Using Cucurbit-[5]-uril Supramolecule: A Relativistic DFT Based Investigation.

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+).

Ortho-7 bound to the active-site gorge of free and OP-conjugated acetylcholinesterase: cation-π interactions.

Despite the immense importance of cation-π interactions prevailing in bispyridinium drug acetylcholinesterase (AChE) complexes, a precise description of cation-π interactions at molecular level has remained elusive. Here, we consider a bispyridinium drug, namely, ortho-7 in three different structures of AChE, with and without complexation with organophosphorus (OP) compounds for detailed investigation using all atom molecular dynamics simulation. By quantum mechanical calculations, Y72, W86, Y124, W286, Y337, and Y341 aromatic residues of the enzyme are investigated for possible cation-π interactions with ortho-7. The cation-π interactions in each of the protein-drug complexes are studied using distance, angle, a suitable functional form of them, and electrostatic criteria. The variation of cation-π functional is remarkably consistent with that of the Columbic variation. It is clearly observed that cation-π interactions for some of the residues in the catalytic active site (CAS) and peripheral anionic site (PAS) of the enzyme are either enhanced or reduced based on the nature of OP conjugation (i.e., nerve gas, tabun or pesticide, fenamiphos) when compared with the OP-free enzyme. The strength of cation-π interaction is strongly dependent on the type OP conjugation. The effect of conjugation at CAS is also seen to influence the cation-π interaction at the PAS region. The variation of cation-π interactions on the type of conjugating OP compounds might be suggestive of a reason as to why wide spectrum drug against any OP poisoning is yet to arrive in the market.

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.

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

Ion channels, such as potassium channels are known to discriminate ions to achieve remarkable selective transportation of K(+) over Na(+) through the membrane. The recently reported NaK ion channel, on the contrary, seems to be an exception, as it is observed to permeate most of the group IA alkali metal cations and hence is suggested to be nonselective in nature. However, does that correspond to a complete annihilation of selectivity inside the selectivity filter (SF) of the channel? What is the origin of such nonselectivity/selectivity, if any? The present computational study is an extensive multiscale modeling approach to find the probable answers to these intriguing questions. Here, we have used density functional theory (DFT) based calculations using a realistic truncated model of SF from the crystal structures of the NaK ion channel to evaluate the binding of various alkali metal ions (Na(+), K(+) and Cs(+)), free from "contamination" due to the absence any other "rivalry" cations, in its different binding sites. Among all of the possible binding sites, a vestibule is noticed to be nonselective and seen to act as a probable binding site only in the presence of multiple ions. Binding sites S3 and S4 are found to be selective for K(+) and Na(+), respectively. As an important observation, we find that calculations on oversimplified models using an isolated ion binding site may lead to an erroneous selectivity trend as it neglects the synergetics of consecutive binding sites on the final outcome. Energy decomposition analysis revealed ion-dipole electrostatics as the major contributing interaction in metal-bound binding sites. Our investigations find that although NaK is permeable to monovalent alkali metal ions, strongly "site specific" selectivity does exist at the three well-defined noncontiguous binding sites of the SF. Different important physicomechanical parameters (such as ligating environment, synergistic influence of binding sites, and topological constraints) are found to be the determining factor to induce the "site specific" selectivity of ions during translocation. Wherever possible, our computed results are compared with the available experimental findings. We finally conduct a detailed umbrella sampling-corrected metadynamics simulation in order to obtain an ion permeation free energy landscape within the SF that corroborates well with the "site specific" selectivity trend.

Water-Mediated Differential Binding of Strontium and Cesium Cations in Fulvic Acid.

The migration of potentially harmful radionuclides, such as cesium ((137)Cs) and strontium ((90)Sr), in soil is governed by the chemical and biological reactivity of soil components. Soil organic matter (SOM) that can be modeled through fulvic acid (FA) is known to alter the mobility of radionuclide cations, Cs(+) and Sr(2+). Shedding light on the possible interaction mechanisms at the atomic level of these two ions with FA is thus vital to explain their transport behavior and for the design of new ligands for the efficient extraction of radionuclides. Here we have performed molecular dynamics, metadynamics simulations, and density-functional-theory-based calculations to understand the binding mechanism of Sr(2+) and Cs(+) cations with FA. Our studies predict that interaction of Cs(+) to FA is very weak as compared with Sr(2+). While the water-FA interaction is largely responsible for the weak binding of Cs(+) to FA, leading to the outer sphere complexation of the ion with FA, the interaction between Sr(2+) and FA is stronger and thus can surpass the existing secondary nonbonding interaction between coordinated waters and FA, leading to inner sphere complexation of the ion with FA. We also find that entropy plays a dominant role for Cs(+) binding to FA, whereas Sr(2+) binding is an enthalpy-driven process. Our predicted results are found to be in excellent agreement with the available experimental data on complexation of Cs(+) and Sr(2+) with SOM.