Antibiotic Charge Profile Determines the Extent of L3 Dynamics in OmpF: An Expedited Passage for Molecules with a Positive Charge.

Efficient permeation into Gram-negative bacterial cells is a much-desired property in the design of antibacterial agents. The goal is to arrive at one or more chemical modifications of molecules that improve their uptake into the cell while maintaining a good binding affinity to the intracellular target. Previously, we proposed a mechanistic rationale for the fast permeation of bulky antibiotics that involves induced conformational dynamics in the constriction loop L3 of the OmpF channel. This flexibility is caused by the perturbation of a hydrogen bond network stabilizing the L3 loop due to the strong interactions of the positively charged moiety on the antibiotic with the residues of the L3 loop. In the present work, we examine how differences in the charge profile of antibiotic molecules can affect the permeation process, in particular, the L3 dynamics. To this end, we have performed all-atom molecular dynamics simulations to study the permeation process of molecules with differences in the net charge through the Escherichia coli OmpF channel. The results from these simulations suggest that a positively charged moiety on the antibiotic is responsible for strong interactions with the negatively charged residues of the L3 loop, promoting conformational dynamics in the L3 loop. In contrast, antibiotics without a positively charged moiety are unable to initiate such a dynamic response in the L3 loop. This distinct behavior of the L3 loop in the presence of molecules with different charge characteristics provides a plausible mechanism whereby large molecules with an appropriate charge distribution can leverage an L3 dynamic-dependent pathway to permeate efficiently. The results are relevant to the structure-based design of molecules with improved uptake properties achieved through systematic chemical modifications that effectively engage the L3 loop.

BtuB TonB-dependent transporters and BtuG surface lipoproteins form stable complexes for vitamin B12 uptake in gut Bacteroides.

Vitamin B12 (cobalamin) is required for most human gut microbes, many of which are dependent on scavenging to obtain this vitamin. Since bacterial densities in the gut are extremely high, competition for this keystone micronutrient is severe. Contrasting with Enterobacteria, members of the dominant genus Bacteroides often encode several BtuB vitamin B12 outer membrane transporters together with a conserved array of surface-exposed B12-binding lipoproteins. Here we show that the BtuB transporters from Bacteroides thetaiotaomicron form stable, pedal bin-like complexes with surface-exposed BtuG lipoprotein lids, which bind B12 with high affinities. Closing of the BtuG lid following B12 capture causes destabilisation of the bound B12 by a conserved BtuB extracellular loop, causing translocation of the vitamin to BtuB and subsequent transport. We propose that TonB-dependent, lipoprotein-assisted small molecule uptake is a general feature of Bacteroides spp. that is important for the success of this genus in colonising the human gut.

Fast prediction of antibiotic permeability through membrane channels using Brownian dynamics.

The efficient permeation across the Gram-negative bacterial membrane is an important step in the overall process of antibacterial action of a molecule and the one that has posed a significant hurdle on the way toward approved antibiotics. Predicting the permeability for a large library of molecules and assessing the effect of different molecular transformations on permeation rates of a given molecule is critical to the development of effective antibiotics. We present a computational approach for obtaining estimates of molecular permeability through a porin channel in a matter of hours using a Brownian dynamics approach. The fast sampling using a temperature acceleration scheme enables the approximate estimation of permeability using the inhomogeneous solubility diffusion model. Although the method is a significant approximation to similar all-atom approaches tested previously, we show that the present approach predicts permeabilities that correlate fairly well with the respective experimental permeation rates from liposome swelling experiments and accumulation rates from antibiotic accumulation assays, and is significantly, i.e., about 14 times, faster compared with a previously reported approach. The possible applications of the scheme in high-throughput screening for fast permeators are discussed.

Assembly of transmembrane pores from mirror-image peptides.

Tailored transmembrane alpha-helical pores with desired structural and functional versatility have promising applications in nanobiotechnology. Herein, we present a transmembrane pore DpPorA, based on the natural pore PorACj, built from D-amino acid α-helical peptides. Using single-channel current recordings, we show that DpPorA peptides self-assemble into uniform cation-selective pores in lipid membranes and exhibit properties distinct from their L-amino acid counterparts. DpPorA shows resistance to protease and acts as a functional nanopore sensor to detect cyclic sugars, polypeptides, and polymers. Fluorescence imaging reveals that DpPorA forms well-defined pores in giant unilamellar vesicles facilitating the transport of hydrophilic molecules. A second D-amino acid peptide based on the polysaccharide transporter Wza forms transient pores confirming sequence specificity in stable, functional pore formation. Finally, molecular dynamics simulations reveal the specific alpha-helical packing and surface charge conformation of the D-pores consistent with experimental observations. Our findings will aid the design of sophisticated pores for single-molecule sensing related technologies.

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.

The mechanism of conversion of substituted glycals to chiral acenes via Diels-Alder reaction: a computational study.

Synthesis of linearly fused aromatic systems using a glycal-based diene with an aryne is a long-standing topic of interest in glycal chemistry. We have examined the mechanistic pathways for the transformation of substituted glycals to chiral fused aromatic cores via Diels-Alder (DA) reaction using the SMDACN-M06-2X/6-31G(d) level of theory. The DA reactions of E (1a) and Z (1a') forms of C-2 alkenyl glycal and an aryl glycal (1b) as a diene were examined with a benzyne intermediate generated as a dienophile. The computational results reveal that 1a and 1b can only be transformed into the fused aromatic cores by the base-catalyzed reaction because a [1,5] sigmatropic hydrogen shift is not feasible. The activation free energy barrier for the base-catalyzed proton abstraction process is 4.2 kcal mol-1 and there is almost no barrier for stereoisomeric 1a DA-complexes. The activation free energy barrier values for stereoisomeric 1b DA-complexes for the base-catalyzed proton abstraction process are 10.8 and 12.4 kcal mol-1. The appropriate orientation of glycal-ring-oxygen and hydrogen at the 5th position of Z (1a') forms of C-2 alkenyl glycal facilitates the [1,5] sigmatropic hydrogen shift; however, the base-catalyzed reaction is energetically more favored than the former case. The rate-determining step for 1a and 1a' is the ring-opening step (18.2 and 19.5 kcal mol-1 for the S-stereoisomer), whereas the DA adduct formation step is the rate-determining step for 1b (16.1 kcal mol-1 for the S-stereoisomer). The structural analysis reveals the formation of the preferred S-stereoisomer over the R-stereoisomer with the respective dienes.

Common mechanism of thermostability in small α- and ß-proteins studied by molecular dynamics.

Protein thermostability is important to evolution, diseases, and industrial applications. Proteins use diverse molecular strategies to achieve stability at high temperature, yet reducing the entropy of unfolding seems required. We investigated five small α-proteins and five ß-proteins with known, distinct structures and thermostability (Tm ) using multi-seed molecular dynamics simulations at 300, 350, and 400 K. The proteins displayed diverse changes in hydrogen bonding, solvent exposure, and secondary structure with no simple relationship to Tm . Our dynamics were in good agreement with experimental B-factors at 300 K and insensitive to force-field choice. Despite the very distinct structures, the native-state (300 + 350 K) free-energy landscapes (FELs) were significantly broader for the two most thermostable proteins and smallest for the three least stable proteins in both the α- and ß-group and with both force fields studied independently (tailed t-test, 95% confidence level). Our results suggest that entropic ensembles stabilize proteins at high temperature due to reduced entropy of unfolding, viz., ΔG = ΔH - TΔS. Supporting this mechanism, the most thermostable proteins were also the least kinetically stable, consistent with broader FELs, typified by villin headpiece and confirmed by specific comparison to a mesophilic ortholog of Thermus thermophilus apo-pyrophosphate phosphohydrolase. We propose that molecular strategies of protein thermostabilization, although diverse, tend to converge toward highest possible entropy in the native state consistent with the functional requirements. We speculate that this tendency may explain why many proteins are not optimally structured and why molten-globule states resemble native proteins so much.

Flexible Fitting of Small Molecules into Electron Microscopy Maps Using Molecular Dynamics Simulations with Neural Network Potentials.

Despite significant advances in resolution, the potential for cryo-electron microscopy (EM) to be used in determining the structures of protein-drug complexes remains unrealized. Determination of accurate structures and coordination of bound ligands necessitates simultaneous fitting of the models into the density envelopes, exhaustive sampling of the ligand geometries, and, most importantly, concomitant rearrangements in the side chains to optimize the binding energy changes. In this article, we present a flexible-fitting pipeline where molecular dynamics flexible fitting (MDFF) is used to refine structures of protein-ligand complexes from 3 to 5 Å electron density data. Enhanced sampling is employed to explore the binding pocket rearrangements. To provide a model that can accurately describe the conformational dynamics of the chemically diverse set of small-molecule drugs inside MDFF, we use QM/MM and neural-network potential (NNP)/MM models of protein-ligand complexes, where the ligand is represented using the QM or NNP model, and the protein is represented using established molecular mechanical force fields (e.g., CHARMM). This pipeline offers structures commensurate to or better than recently submitted high-resolution cryo-EM or X-ray models, even when given medium to low-resolution data as input. The use of the NNPs makes the algorithm more robust to the choice of search models, offering a radius of convergence of 6.5 Å for ligand structure determination. The quality of the predicted structures was also judged by density functional theory calculations of ligand strain energy. This strain potential energy is found to systematically decrease with better fitting to density and improved ligand coordination, indicating correct binding interactions. A computationally inexpensive protocol for computing strain energy is reported as part of the model analysis protocol that monitors both the ligand fit as well as model quality.

Structural insights into the interactions of flavin mononucleotide (FMN) and riboflavin with FMN riboswitch: a molecular dynamics simulation study.

Antibiotics resistance is becoming a serious problem associated with fatalities and suffering patients. New antibiotics that can target the broader spectrum of cellular processes are warranted. One of the recent approaches in this regard is to target the special type of RNA riboswitches in bacteria. In this report, we have explored the mechanistic pathways of ligand-dependent conformational changes of flavin mononucleotide (FMN) riboswitch using molecular dynamics (MD) simulation studies. Cognate ligands FMN and riboflavin (RBF) have shown very different behavior with FMN riboswitch in terms of their role in the gene regulation process. These two ligands have similar scaffold, except the terminal phosphate group in FMN ligand. The MD simulations reveal that the binding of FMN ligand with the riboswitch does not lead to global folding of structure, rather lead to local changes in riboswitch structure. The binding free energy calculated with molecular mechanics Poisson-Boltzmann surface area method suggests the stronger binding of FMN than RBF to the riboswitch and electrostatic energy contributes chiefly to stabilize the complex. Further, the hydrogen bonding analysis identified the key binding site residues G11, G32, G62 of the riboswitch with FMN and RBF. The critical role of the phosphate group in the FMN ligand for binding with the active site of a riboswitch is also borne out in this study. These results unravel the importance of functional groups in natural ligands on designing newer ligands for FMN riboswitch as new antibiotics in the future.Communicated by Ramaswamy H. Sarma.

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

Communicated by Ramaswamy H. Sarma.

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.

Indium-Catalyzed Denitrogenative Transannulation of Pyridotriazoles: Synthesis of Pyrido[1,2- a]indoles.

Pyrido[1,2- a]indoles are known for medicinally and pharmaceutically important compounds; however, the efficient synthetic routes are scarce in the literature. We report herein a convenient and efficient route to synthesize these molecules through indium-catalyzed transannulation of pyridotriazoles with arenes. A library of compounds have been synthesized employing the method developed with various substituted pyrido[1,2- a]indole derivatives in moderate to good yields. The density functional theory study using SMDDCB-M06/6-31++G(d,p)/LANL2DZ//B3LYP/6-31G(d)/LANL2DZ method suggests that the reactions proceed via indium-carbenoid complex.

DFT studies on quantum mechanical tunneling in tautomerization of three-membered rings.

We have examined keto-enol and amino-imino tautomerization in a set of three-membered ring systems (1-5) in the absence and presence of water molecules. Aromaticity governs the keto-enol and amino-imino tautomerization processes in (1-5), which lead to the formation of enol and imine derivatives. The possibility of quantum mechanical tunneling (QMT) has not been reported in the tautomerization processes of three-membered ring systems. M062X/6-311+G(d,p) level of theory QMT calculations reveal that tunneling is not possible in the water unassisted processes because of very high free energy activation barriers. The activation free energy barriers for the amino-imino tautomerization of 5, aziridine-2,3-diimine, and one water assisted, 5-W, are 58.1 kcal mol-1 and 14.8 kcal mol-1, respectively and the lowest among the 3-membered rings examined. The classical over the barrier rate constant (kCVT) obtained by QMT calculation for 5-W→5-W-P is 10.6 s-1 at 273 K. Inclusion of small curvature tunneling (SCT) enhances the classical over the barrier rate constant by 15.1 times at 273 K, i.e., kCVT+SCT is 160 s-1 and reveals nonclassical behaviour for the tautomerization of 5-W. A higher kinetic isotope effect in the tautomerization process of 4-W and 5-W also indicates a pronounced contribution of tunneling toward the tautomerization process. The two-water assisted tautomerization of 3 has the highest activation free energy barrier in the series indicating a nonclassical contribution to 3-2W→3-2W-P. These results suggest that the tautomerization processes of 1-5 are experimentally feasible by tunneling and aromaticity.

Influence of gauche effect on uncharged oxime reactivators for the reactivation of tabun-inhibited AChE: quantum chemical and steered molecular dynamics studies.

The neutral oxime reactivator RS194B with a seven-membered ring has shown better efficacy towards the tabun-inhibited AChE than that of RS69N with a six-membered ring and RS41A with a five-membered ring. The difference in the efficacy of these reactivators has remained unexplored. We report here the origin of the difference of efficacy of these reactivators based on the conformational analysis, quantum chemical calculations and steered molecular dynamics (SMD) simulations. The conformational analysis using B3LYP/6-31G(d) level of theory revealed that RS41A and RS194B are more stable in gauche conformation due to the gauche effect (-N-C-C-N- bonds) whereas RS69N prefers anti-conformation. The SMD simulations show that RS194B retains in more stable gauche conformation inside the active gorge of AChE during different time intervals that experiences more hydrogen bonding, hydrophobic interactions with the catalytic anionic site (CAS) residues and weaker interactions with the peripheral anionic site (PAS) residues compared to RS41A and RS69N. In an effort to design an even superior reactivator, RS194B-S has been chosen with a subtle change in the geometry of RS194B by replacing the carbonyl oxygen with the sulfur atom. The newly designed reactivator RS194B-S can also be a promising candidate to reactivate tabun-inhibited AChE.

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.

Instant Detection of Hydrogen Cyanide Gas and Cyanide Salts in Solid Matrices and Water by using CuII and NiII Complexes of Intramolecularly Hydrogen Bonded Zwitterions.

A series of intramolecularly hydrogen-bonded zwitterionic compartmental ligands HL1-HL4, containing a pendent diamine arm that is monoprotonated and an aldehyde functionality at two different ortho-positions of a 4-halophenoxide, is reported herein. Single-crystal X-ray diffraction (SXRD) provides persuasive evidence for the identification of this class of proton-transferred zwitterions at room temperature. The solid-state photoluminescent nature of these zwitterions remains intact in aqueous and organic solutions. Grinding of HL1 and HL2 with Cu2+ /Ni2+ salts develop turn-on probes 1-4. Compounds 1 and 4 are dinuclear CuII and NiII species, respectively. Compound 2 is a tetranuclear CuII complex. Interestingly, compound 3 is a mononuclear NiII species in which both nitrogen atoms in the pendant diamine arm are protonated and, therefore, not coordinated to the NiII center. All these probes (1-4) display an instant response to the poison gas hydrogen cyanide (HCN) and cyanide salts present in both solid matrices and aqueous (100 % water) solution. Selective and rapid sensing of HCN gas and cyanide salts in solid/soil/water phases, without any interference, by the mechanosynthesized complexes 1-4 can be perceived easily by the naked eye under a hand-held UV lamp.

Transformation of Substituted Glycals to Chiral Fused Aromatic Cores via Annulative π-Extension Reactions with Arynes.

The Diels-Alder addition of arynes to appropriately substituted vinyl/aryl glycals followed by π-extension via pyran ring opening smoothly furnished meta-disubstituted fused aromatic cores containing a stereodefined orthogonally protected chiral side chain. The method is broad in terms of aryl homologation, affording benzene, naphthalene, and phenanthrene derivatives. Base-induced deprotonation followed by cleavage of the allylic C-O bond appear to be the crucial steps leading to the development of aromaticity, which is the driving force behind the annulative π-extension process. The present protocol can be used for the synthesis of meta-disubstituted naphthalene aldehydes and substrates for aldolases.

DFT Study To Explore the Importance of Ring Size and Effect of Solvents on the Keto-Enol Tautomerization Process of α- and ß-Cyclodiones.

We have explored the effect of ring size on keto-enol tautomerization of α- and ß-cyclodiones using the M062X-SMDaq/6-31+G(d,p)//M062X/6-31+G(d,p) level of theory. The calculated results show that the activation free energy barrier for the keto-enol tautomerization process of α-cyclopropanedione (1) is 54.9 kcal/mol, which is lower compared to that of the other cyclic diketo systems studied here. The four-membered α- and ß-cyclobutanedione (2 and 6) do not favor keto-enol tautomerization unlike other studied cyclic systems because of the ring strain developed in the transition-state geometries and their corresponding products. Water-assisted keto-enol tautomerization with one molecule reveals that the free energy activation barriers reduce almost half compared to those for the uncatalyzed systems. The two-water-assisted process is favorable as the activation free energy barriers lowered by ∼10 kcal/mol compared to those of the one-water-assisted process. The ion-pair formation seems to govern the lowering of activation barriers of α- and ß-cyclodiones with two water molecules during the keto-enol tautomerization process, which however also overcomes the favorable aromatization in the three-membered ring system. The free energy activation barriers calculated with the M062X-SMDaq/6-31+G(d,p) level predicted that the keto-enol tautomerization process for the α-cyclodiones follows the following trend: 2 > 3 > 4 > 5 > 1. Water-assisted tautomerization of α-cyclodiones also predicted 1-W and 1-2W as the most favored processes; however, 5-W and 5-2W were found to be disfavored in this case. The ß-cyclodione systems also showed similar trends as obtained with α-diketone systems. The influence of bulk solvent on the keto-enol tautomerization process favors the formation of the enol form in a more polar solvent medium even under mixed solvent conditions in acetonitrile and hexane at M062X-SMDacetonitrile/6-31+G(d,p) and M062X-SMDhexane/6-31+G(d,p) levels of theory.

Probing the Structural and Electronic Effects on the Origin of π-Facial Stereoselectivity in 1-Methylphosphole 1-Oxide Cycloadditions and Cyclodimerization.

We have examined the π-facial stereoselectivity in the Diels-Alder reactions of phosphole oxides computationally. The experimentally observed syn-cycloadditions have been rationalized with the Cieplak model and distortion-interaction model. The natural bond orbital analysis suggests that the hyperconjugative interactions are energetically preferred between the antiperiplanar methyl group present in the -P=O unit and the developing incipient (-C-C-) bond in syn-adducts in accordance with the Cieplak model. The distortion-interaction analysis carried out for syn and anti transition states of Diels-Alder reactions of 1-substituted phosphole 1-oxide with different dienophiles reveals that the syn selectivity is favored by distortions and interaction energies compared with the anti selectivity. The formation of a syn adduct is also stabilized by the πCC-σ*PO orbital interaction, and the repulsive n-π interaction destabilizes the anti adduct that leads to the 7.0 kcal/mol thermodynamic preference for the former adduct. Furthermore, the distortion-interaction model rationalizes the formation of stereospecific products in these Diels-Alder reactions, which however is not explicable with the much-debated Cieplak model.

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.