Modeling evolution of hydrogen bonding and stabilization of transition states in the process of cocaine hydrolysis catalyzed by human butyrylcholinesterase

Long-lasting blocking of interoceptive effects of cocaine by a highly efficient cocaine hydrolase in rats

Cocaine dependence is a serious world-wide public health problem without an FDA-approved pharmacotherapy. We recently designed and discovered a highly efficient long-acting cocaine hydrolase CocH5-Fc(M6). The present study examined the effectiveness and duration of CocH5-Fc(M6) in blocking interoceptive effects of cocaine by performing cocaine discrimination tests in rats, demonstrating that the duration of CocH5-Fc(M6) in blocking cocaine discrimination was dependent on cocaine dose and CocH5-Fc(M6) plasma concentration. Particularly, a dose of 3 mg/kg CocH5-Fc(M6) effectively attenuated discriminative stimulus effects of 10 mg/kg cocaine, cumulative doses of 10 and 32 mg/kg cocaine, and cumulative doses of 10, 32 and 56 mg/kg cocaine by ≥ 20% for 41, 19, and 10 days, and completely blocked the discriminative stimulus effects for 30, 13, and 5 days with corresponding threshold plasma CocH5-Fc(M6) concentrations of 15.9, 72.2, and 221 nM, respectively, under which blood cocaine concentration was negligible. Additionally, based on the data obtained, cocaine discrimination model is more sensitive than the locomotor activity to reveal cocaine effects and that CocH5-Fc(M6) itself has no long-term toxicity regarding behavioral activities such as lever pressing and food consumption in rats, further demonstrating that CocH5-Fc(M6) has the desired properties as a promising therapeutic candidate for prevenance of cocaine dependence.

Human Butyrylcholinesterase Mutants for (-)-Cocaine Hydrolysis: A Correlation Relationship between Catalytic Efficiency and Total Hydrogen Bonding Energy with an Oxyanion Hole

A combined computational and experimental study has been carried out to explore and test a quantitative correlation relationship between the relative catalytic efficiency (RCE) of human butyrylcholinesrase (BChE) mutant-catalyzed hydrolysis of substrate (-)-cocaine and the total hydrogen bonding energy (tHBE) of the carbonyl oxygen of the substrate with the oxyanion hole of the enzyme in the modeled transition-state structure (TS1), demonstrating a satisfactory linear correlation relationship between ln(RCE) and tHBE. The satisfactory correlation relationship has led us to computationally predict and experimentally confirm new human BChE mutants that have a further improved catalytic activity against (-)-cocaine, including the most active one (the A199S/F227S/S287G/A328W/Y332G mutant) with a 2790-fold improved catalytic efficiency (kcat/KM = 2.5 × 109 min-1 M-1) compared to the wild-type human BChE. Compared to the reference mutant (the A199S/S287G/A328W/Y332G mutant) tested in the reported clinical development of an enzyme therapy for cocaine dependence treatment, this new mutant (with a newly predicted additional F227S mutation) has an improved catalytic efficiency against (-)-cocaine by ∼2.6-fold. The good agreement between the computational and experimental ln(RCE) values suggests that the obtained correlation relationship is robust for computational prediction. A similar correlation relationship could also be explored in studying BChE or other serine hydrolases/esterases with an oxyanion hole stabilizing the carbonyl oxygen in the rate-determining reaction step of the enzymatic hydrolysis of other substrates.

A Comprehensive Review of Cholinesterase Modeling and Simulation

The present article reviews published efforts to study acetylcholinesterase and butyrylcholinesterase structure and function using computer-based modeling and simulation techniques. Structures and models of both enzymes from various organisms, including rays, mice, and humans, are discussed to highlight key structural similarities in the active site gorges of the two enzymes, such as flexibility, binding site location, and function, as well as differences, such as gorge volume and binding site residue composition. Catalytic studies are also described, with an emphasis on the mechanism of acetylcholine hydrolysis by each enzyme and novel mutants that increase catalytic efficiency. The inhibitory activities of myriad compounds have been computationally assessed, primarily through Monte Carlo-based docking calculations and molecular dynamics simulations. Pharmaceutical compounds examined herein include FDA-approved therapeutics and their derivatives, as well as several other prescription drug derivatives. Cholinesterase interactions with both narcotics and organophosphate compounds are discussed, with the latter focusing primarily on molecular recognition studies of potential therapeutic value and on improving our understanding of the reactivation of cholinesterases that are bound to toxins. This review also explores the inhibitory properties of several other organic and biological moieties, as well as advancements in virtual screening methodologies with respect to these enzymes.

Rationally engineered prolyl endopeptidases from Sphingomonas capsulata with improved hydrolytic activity towards pathogenic peptides of celiac diseases

Celiac disease affects approximately 1% of the population and is a major public health problem worldwide. It is trigged by gluten-derived peptides, which have unusually high proline-glutamine motif content and are highly resistant to proteolysis by digestive enzymes of the gastrointestinal tract. The only treatment for celiac disease is strict, lifelong adherence to a gluten-free diet, which is effective but costly and difficult to maintain. Therefore, novel non-dietary therapies for celiac disease are urgently needed. Gluten-degrading enzymes are promising non-dietary treatments, and some enzymes have been investigated in preclinical or clinical studies. A combination of prolyl endopeptidase from Sphingomonas capsulata (SC PEP) and a glutamine-specific endoprotease (EP-B2 from barley) known as latiglutenase showed insufficient benefits in phase II clinical trials, likely because of its low enzyme activity in the gastric environment. Therefore, improving enzyme activity is essential for the clinical application of SC PEP. Enzyme activity can be enhanced using computer-aided rational protein design tools. In this study, we combined molecular docking and molecular dynamics simulation to rationally design SC PEP mutants and experimentally evaluated their activities. We identified mutants with up to 90-103% increases in specific activity and up to 80-202% increases in the catalytic rate. We have investigated the mechanism underlying the enhanced activity of these mutants, and found that a conformational transition of the β-propeller domain and catalytic domain of SC PEP was important for enzyme activity, and this transition was affected by residues in the catalytic domain and at the domain interface; a shorter distance between the substrate Pro and the oxyanion holes was also crucial for improving SC PEP catalytic activity. Our results provide useful information for the rational design of highly active SC PEPs to accelerate the development of enzyme therapeutics candidates for Celiac disease.

Unexpected protonation state of Glu197 discovered from simulations of tacrine in butyrylcholinesterase

Butyrylcholinesterase (BChE) has been actively involved in drug discoveries from many fields for decades. In the crystal structure of the BChE-tacrine complex, there is an unanticipated formyl-proline molecule resolved very close to tacrine, raising an essential question on how reliable it is to apply the binding pose in a crystal structure to analyze related experimental observations, in which no formyl-proline is actually involved. In this study, by performing a series of 100 ns molecular dynamics simulations, we demonstrate that it is safe to employ the structural information from this crystal structure to analyze related experimental observations. Surprisingly, Glu197 needs to be protonated to have the structures simulated appropriately. It should be noted that Glu197 has been commonly considered as deprotonated in diverse analyses due to its low pKa in aqueous solution, for which some interpretations are inconsistent or unclear. Our further investigation shows that the protonated Glu197 plays a very important role in preserving His438 within the catalytic triad through stabilizing a highly conserved water molecule. Interestingly, the catalytic triad and Glu197 have been long recognized for possibly deviating largely from the crystal structure, which might be catalytically deficient and is generally considered to result from the difference between the crystal and aqueous environment. Herein, our results suggest that the large deviations of the catalytic triad and Glu197 from the crystal structure are caused by the inappropriate protonation state of Glu197. This finding shall provide an important clue that has been long missing for a better understanding of BChE-related puzzles or even reconsideration of some BChE-catalyzed reaction mechanisms.

Plant-expressed cocaine hydrolase variants of butyrylcholinesterase exhibit altered allosteric effects of cholinesterase activity and increased inhibitor sensitivity

Butyrylcholinesterase (BChE) is an enzyme with broad substrate and ligand specificities and may function as a generalized bioscavenger by binding and/or hydrolyzing various xenobiotic agents and toxicants, many of which target the central and peripheral nervous systems. Variants of BChE were rationally designed to increase the enzyme’s ability to hydrolyze the psychoactive enantiomer of cocaine. These variants were cloned, and then expressed using the magnICON transient expression system in plants and their enzymatic properties were investigated. In particular, we explored the effects that these site-directed mutations have over the enzyme kinetics with various substrates of BChE. We further compared the affinity of various anticholinesterases including organophosphorous nerve agents and pesticides toward these BChE variants relative to the wild type enzyme. In addition to serving as a therapy for cocaine addiction-related diseases, enhanced bioscavenging against other harmful agents could add to the practicality and versatility of the plant-derived recombinant enzyme as a multivalent therapeutic.

Reaction Pathway for Cocaine Hydrolase-Catalyzed Hydrolysis of (+)-Cocaine

A recently designed and discovered cocaine hydrolase (CocH), engineered from human butyrylcholinesterase (BChE), has been proven promising as a novel enzyme therapy for treatment of cocaine overdose and addiction because it is highly efficient in catalyzing hydrolysis of naturally occurring (-)-cocaine. It has been known that the CocH also has a high catalytic efficiency against (+)-cocaine, a synthetic enantiomer of cocaine. Reaction pathway and the corresponding free energy profile for the CocH-catalyzed hydrolysis of (+)-cocaine have been determined, in the present study, by performing first-principles pseudobond quantum mechanical/molecular mechanical (QM/MM)-free energy (FE) calculations. Acordingt to the QM/MM-FE results, the catalytic hydrolysis process is initiated by the nucleophilic attack on carbonyl carbon of (-)-cocaine benzoyl ester via hydroxyl oxygen of S198 side chain, and the second reaction step (i.e. dissociation of benzoyl ester) is rate-determining. This finding for CocH-catalyzed hydrolysis of (+)-cocaine is remarkably different from that for the (+)-cocaine hydrolysis catalyzed by bacterial cocaine esterase in which the first reaction step of the deacylation is associated with the highest free energy barrier (~17.9 kcal/mol). The overall free energy barrier (~16.0 kcal/mol) calculated for the acylation stage of CocH-catalyzed hydrolysis of (+)-cocaine is in good agreement with the experimental free energy barrier of ~14.5 kcal/mol derivated from the experimental kinetic data.

Rational design, preparation, and characterization of a therapeutic enzyme mutant with improved stability and function for cocaine detoxification

Cocaine esterase (CocE) is known as the most efficient natural enzyme for cocaine hydrolysis. The major obstacle to the clinical application of wild-type CocE is the thermoinstability with a half-life of only ∼12 min at 37 °C. The previously designed T172R/G173Q mutant (denoted as enzyme E172-173) with an improved in vitro half-life of ∼6 h at 37 °C is currently in clinical trial Phase II for cocaine overdose treatment. Through molecular modeling and dynamics simulation, we designed and characterized a promising new mutant of E172-173 with extra L196C/I301C mutations (denoted as enzyme E196-301) to produce cross-subunit disulfide bonds that stabilize the dimer structure. The cross-subunit disulfide bonds were confirmed by X-ray diffraction. The designed L196C/I301C mutations have not only considerably extended the in vitro half-life at 37 °C to >100 days, but also significantly improved the catalytic efficiency against cocaine by ∼150%. In addition, the thermostable E196-301 can be PEGylated to significantly prolong the residence time in mice. The PEGylated E196-301 can fully protect mice from a lethal dose of cocaine (180 mg/kg, LD100) for at least 3 days, with an average protection time of ∼94h. This is the longest in vivo protection of mice from the lethal dose of cocaine demonstrated within all studies using an exogenous enzyme reported so far. Hence, E196-301 may be developed to become a more valuable therapeutic enzyme for cocaine abuse treatment, and it demonstrates that a general design strategy and protocol to simultaneously improve both the stability and function are feasible for rational protein drug design.

A highly efficient cocaine-detoxifying enzyme obtained by computational design

Compared to naturally occurring enzymes, computationally designed enzymes are usually much less efficient, with their catalytic activities being more than six orders of magnitude below the diffusion limit. Here we use a two-step computational design approach, combined with experimental work, to design a highly efficient cocaine hydrolising enzyme. We engineer E30-6 from human butyrylcholinesterase (BChE), which is specific for cocaine hydrolysis, and obtain a much higher catalytic efficiency for cocaine conversion than for conversion of the natural BChE substrate, acetylcholine (ACh). The catalytic efficiency of E30-6 for cocaine hydrolysis is comparable to that of the most efficient known naturally-occurring hydrolytic enzyme, acetylcholinesterase, the catalytic activity of which approaches the diffusion limit. We further show that E30-6 can protect mice from a subsequently administered lethal dose of cocaine, suggesting the enzyme may have therapeutic potential in the setting of cocaine detoxification or cocaine abuse.

Substrate selectivity of high-activity mutants of human butyrylcholinesterase

Cocaine is one of the most addictive drugs, and there is still no FDA (Food and Drug Administration)-approved medication specific for cocaine abuse. A promising therapeutic strategy is to accelerate cocaine metabolism, producing biologically inactive metabolites via a route similar to the primary cocaine-metabolizing pathway, i.e. cocaine hydrolysis catalyzed by butyrylcholinesterase (BChE) in plasma. However, the native BChE has a low catalytic efficiency against the abused cocaine, i.e. (-)-cocaine. Our recently designed and discovered A199S/F227A/S287G/A328W/Y332G mutant and other mutants of human BChE have a considerably improved catalytic efficiency against (-)-cocaine. In the present study, we carried out both computational modeling and experimental kinetic analysis on the catalytic activities of these promising new BChE mutants against other known substrates, including neurotransmitter acetylcholine (ACh), acetylthiocholine (ATC), butyrylthiocholine (BTC), and (+)-cocaine, in comparison with the corresponding catalytic activity against (-)-cocaine. Both the computational modeling and kinetic analysis have consistently revealed that all the examined amino acid mutations only considerably improve the catalytic efficiency of human BChE against (-)-cocaine, without significantly improving the catalytic efficiency of the enzyme against any of the other substrates examined. In particular, all the examined BChE mutants have a slightly lower catalytic efficiency against neurotransmitter ACh compared to the wild-type BChE. This observation gives us confidence in developing an anti-cocaine enzyme therapy by using one of these BChE mutants, particularly the A199S/F227A/S287G/A328W/Y332G mutant.

Reaction Pathway and Free Energy Profile for Cocaine Hydrolase-Catalyzed Hydrolysis of (-)-Cocaine

Reaction pathway of (-)-cocaine hydrolysis catalyzed by our recently discovered most efficient cocaine hydrolase, which is the A199S/F227A/S287G/A328W/Y332G mutant of human butyrylcholinesterase (BChE), and the corresponding free energy profile have been studied by performing first-principles pseudobond quantum mechanical/molecular mechanical (QM/MM)-free energy (FE) calculations. Based on the QM/MM-FE results, the catalytic hydrolysis process consists of four major reaction steps, including the nucleophilic attack on carbonyl carbon of (-)-cocaine benzoyl ester by hydroxyl group of S198, dissociation of (-)-cocaine benzoyl ester, nucleophilic attack on carbonyl carbon of (-)-cocaine benzoyl ester by water, and finally the dissociation between (-)-cocaine benzoyl group and S198 of the enzyme. The second reaction step is rate-determining. The calculated free energy barrier associated with the transition state for the rate-determining step is ~15.0 kcal/mol, which is in excellent agreement with the experimentally-derived activation free energy of ~14.7 kcal/mol. The mechanistic insights obtained from the present study will be valuable for rational design of more active cocaine hydrolase against (-)-cocaine. In particular, future efforts aiming at further increasing the catalytic activity of the enzyme against (-)-cocaine should focus on stabilization of the transition state for the second reaction step in which the benzoyl ester of (-)-cocaine dissociates.

Reaction pathway and free energy profiles for butyrylcholinesterase-catalyzed hydrolysis of acetylthiocholine

The catalytic mechanism for butyrylcholineserase (BChE)-catalyzed hydrolysis of acetylthiocholine (ATCh) has been studied by performing pseudobond first-principles quantum mechanical/molecular mechanical-free energy (QM/MM-FE) calculations on both acylation and deacylation of BChE. Additional quantum mechanical (QM) calculations have been carried out, along with the QM/MM-FE calculations, to understand the known substrate activation effect on the enzymatic hydrolysis of ATCh. It has been shown that the acylation of BChE with ATCh consists of two reaction steps including the nucleophilic attack on the carbonyl carbon of ATCh and the dissociation of thiocholine ester. The deacylation stage includes nucleophilic attack of a water molecule on the carboxyl carbon of substrate and dissociation between the carboxyl carbon of substrate and hydroxyl oxygen of Ser198 side chain. QM/MM-FE calculation results reveal that the acylation of BChE is rate-determining. It has also been demonstrated that an additional substrate molecule binding to the peripheral anionic site (PAS) of BChE is responsible for the substrate activation effect. In the presence of this additional substrate molecule at PAS, the calculated free energy barrier for the acylation stage (rate-determining step) is decreased by ~1.7 kcal/mol. All of our computational predictions are consistent with available experimental kinetic data. The overall free energy barriers calculated for BChE-catalyzed hydrolysis of ATCh at regular hydrolysis phase and substrate activation phase are ~13.6 and ~11.9 kcal/mol, respectively, which are in reasonable agreement with the corresponding experimentally derived activation free energies of 14.0 kcal/mol (for regular hydrolysis phase) and 13.5 kcal/mol (for substrate activation phase).

Human butyrylcholinesterase-cocaine binding pathway and free energy profiles by molecular dynamics and potential of mean force simulations

In the present study, we have performed combined molecular dynamics and potential of mean force (PMF) simulations to determine the enzyme-substrate (ES) binding pathway and the corresponding free energy profiles for wild-type butyrylcholinesterase (BChE) binding with (-)/(+)-cocaine and for the A328W/Y332G mutant binding with (-)-cocaine. According to the PMF simulations, for each ES binding system, the substrate first binds with the enzyme at a peripheral anionic site around the entrance of the active-site gorge to form the first ES complex (ES1-like) during the binding process. Further evolution from the ES1-like complex to the nonprereactive ES complex is nearly barrierless, with a free energy barrier lower than 1.0 kcal/mol. So, the nonprereactive ES binding process should be very fast. The rate-determining step of the entire ES binding process is the subsequent evolution from the nonprereactive ES complex to the prereactive ES complex. Further accounting for the entire ES binding process, the PMF-based simulations qualitatively reproduced the relative order of the experimentally derived binding free energies (ΔG(bind)), although the simulations systematically overestimated the magnitude of the binding affinity and systematically underestimated the differences between the ΔG(bind) values. The obtained structural and energetic insights into the entire ES binding process provide a valuable base for future rational design of high-activity mutants of BChE as candidates for an enzyme therapy for cocaine overdose and abuse.

Recent progress in protein drug design and discovery with a focus on novel approaches to the development of anti-cocaine medications

Cocaine is highly addictive and no anti-cocaine medication is currently available. Accelerating cocaine metabolism, producing biologically inactive metabolites, is recognized as an ideal anti-cocaine medication strategy, especially for the treatment of acute cocaine toxicity. However, currently known wild-type enzymes have either too low a catalytic efficiency against the abused cocaine, in other words (-)-cocaine, or the in vivo half-life is too short. Novel computational strategies and design approaches have been developed recently to design and discover thermostable or high-activity mutants of enzymes based on detailed structures and catalytic/inactivation mechanisms. The structure- and mechanism-based computational design efforts have led to the discovery of high-activity mutants of butyrylcholinesterase and thermostable mutants of cocaine esterase as promising anti-cocaine therapeutics. The structure- and mechanism-based computational strategies and design approaches may be used to design high-activity and/or thermostable mutants of many other proteins that have clear therapeutic potentials and to design completely new therapeutic enzymes.

Reaction mechanism for cocaine esterase-catalyzed hydrolyses of (+)- and (-)-cocaine: unexpected common rate-determining step

First-principles quantum mechanical/molecular mechanical free energy calculations have been performed to examine the catalytic mechanism for cocaine esterase (CocE)-catalyzed hydrolysis of (+)-cocaine in comparison with CocE-catalyzed hydrolysis of (-)-cocaine. It has been shown that the acylation of (+)-cocaine consists of nucleophilic attack of the hydroxyl group of Ser117 on the carbonyl carbon of (+)-cocaine benzoyl ester and the dissociation of (+)-cocaine benzoyl ester. The first reaction step of deacylation of (+)-cocaine, which is identical to that of (-)-cocaine, is rate-determining, indicating that CocE-catalyzed hydrolyses of (+)- and (-)-cocaine have a common rate-determining step. The computational results predict that the catalytic rate constant of CocE against (+)-cocaine should be the same as that of CocE against (-)-cocaine, in contrast with the remarkable difference between human butyrylcholinesterase-catalyzed hydrolyses of (+)- and (-)-cocaine. The prediction has been confirmed by experimental kinetic analysis on CocE-catalyzed hydrolysis of (+)-cocaine in comparison with CocE-catalyzed hydrolysis of (-)-cocaine. The determined common rate-determining step indicates that rational design of a high-activity mutant of CocE should be focused on the first reaction step of the deacylation. Furthermore, the obtained mechanistic insights into the detailed differences in the acylation between the (+)- and (-)-cocaine hydrolyses provide indirect clues for rational design of amino acid mutations that could more favorably stabilize the rate-determining transition state in the deacylation and, thus, improve the catalytic activity of CocE. This study provides a valuable mechanistic base for rational design of an improved esterase for therapeutic treatment of cocaine abuse.

Reaction pathway and free energy profile for prechemical reaction step of human butyrylcholinesterase-catalyzed hydrolysis of (-)-cocaine by combined targeted molecular dynamics and potential of mean force simulations

Combined targeted molecular dynamics (TMD) and potential of mean force (PMF) simulations have been carried out to uncover the detailed pathway and determine the corresponding free energy profile for the structural transformation from the nonprereactive butyrylcholinesterase (BChE)-(-)-cocaine binding to the prereactive BChE-(-)-cocaine binding associated with the (-)-cocaine rotation in the binding pocket of BChE. It has been shown that the structural transformation involves two transition states (TS1(rot) and TS2(rot)). TS1(rot) is mainly associated with the deformation of the nonprereactive complex, whereas TS2(rot) is mainly associated with the formation of the prereactive complex. It has also been demonstrated that the A328W/Y332G mutation significantly reduces the steric hindrance for (-)-cocaine rotation in the binding pocket of BChE and, thus, decreases the free energy barrier for the structural transformation from the nonprereactive binding to the prereactive binding. The calculated relative free energy barriers are all consistent with available experimental kinetic data. The new mechanistic insights obtained and the novel computational protocol tested in this study should be valuable for future computational design of high-activity mutants of BChE. The general computational strategy and approach based on the combined TMD and PMF simulations may be also valuable in computational studies of detailed pathways and free energy profiles for other similar mechanistic problems involving ligand rotation or another type of structural transformation in the binding pocket of a protein.

Design, preparation, and characterization of high-activity mutants of human butyrylcholinesterase specific for detoxification of cocaine

Cocaine is a widely abused drug without a U.S. Food and Drug Administration-approved medication. There is a recognized, promising anticocaine medication to accelerate cocaine metabolism, producing biologically inactive metabolites via a route similar to the primary cocaine-metabolizing pathway [i.e., cocaine hydrolysis catalyzed by butyrylcholinesterase (BChE) in plasma]. An ideal, therapeutically valuable mutant of human BChE should have not only a significantly improved catalytic activity against (-)-cocaine but also certain selectivity for (-)-cocaine over neurotransmitter acetylcholine (ACh), such that one would not expect systemic administration of the BChE mutant to interrupt cholinergic transmission. The present study accounting for the mutation-caused changes of the catalytic activities of BChE against both (-)-cocaine and ACh by means of molecular modeling and site-directed mutagenesis has led to identification of three BChE mutants that have not only a considerably improved catalytic efficiency against (-)-cocaine but also the desirable selectivity for (-)-cocaine over ACh. Two representative BChE mutants have been confirmed to be potent in actual protection of mice from acute toxicity (convulsion and lethality) of a lethal dose of cocaine (180 mg/kg). Pretreatment with the BChE mutant (i.e., 1 min before cocaine administration) dose-dependently protected mice against cocaine-induced convulsions and lethality. In particular, all mice pretreated with the mutant (e.g., 0.02 mg or more of A199S/F227A/S287G/A328W/E441D BChE) survived. The in vivo data reveal the primary factor (i.e., the relative catalytic efficiency), determining the efficacy in practical protection of mice from the acute cocaine toxicity and future direction for further improving the efficacy of the enzyme in the cocaine overdose treatment.

Reaction pathway and free energy profile for butyrylcholinesterase-catalyzed hydrolysis of acetylcholine

A catalytic mechanism for the butyrylcholinesterase (BChE)-catalyzed hydrolysis of acetylcholine (ACh) has been studied by performing pseudobond first-principles quantum mechanical/molecular mechanical-free energy calculations on both acylation and deacylation of BChE. It has been shown that the acylation with ACh includes two reaction steps, including nucleophilic attack on the carbonyl carbon of ACh and dissociation of choline ester. The deacylation stage includes nucleophilic attack of a water molecule on the carboxyl carbon of the substrate and dissociation between the carboxyl carbon of the substrate and the hydroxyl oxygen of the Ser198 side chain. Notably, despite the fact that acetylcholinesterase (AChE) and BChE are very similar enzymes, the acylation of BChE with ACh is rate-determining, which is remarkably different from the AChE-catalyzed hydrolysis of ACh, in which the deacylation is rate-determining. The computational prediction is consistent with available experimental kinetic data. The overall free energy barrier calculated for BChE-catalyzed hydrolysis of ACh is 13.8 kcal/mol, which is in good agreement with the experimentally derived activation free energy of 13.3 kcal/mol.

Design of high-activity mutants of human butyrylcholinesterase against (-)-cocaine: structural and energetic factors affecting the catalytic efficiency

The present study was aimed to explore the correlation between the protein structure and catalytic efficiency of butyrylcholinesterase (BChE) mutants against (-)-cocaine by modeling the rate-determining transition state (TS1), i.e., the transition state for the first step of chemical reaction process, of (-)-cocaine hydrolysis catalyzed by various mutants of human BChE in comparison with the wild type. Molecular modeling of the TS1 structures revealed that mutations on certain nonactive site residues can indirectly affect the catalytic efficiency of the enzyme against (-)-cocaine through enhancing or weakening the overall hydrogen bonding between the carbonyl oxygen of (-)-cocaine benzoyl ester and the oxyanion hole of the enzyme. Computational insights and predictions were supported by the catalytic activity data obtained from wet experimental tests on the mutants of human BChE, including five new mutants reported for the first time. The BChE mutants with at least ∼1000-fold improved catalytic efficiency against (-)-cocaine compared to the wild-type BChE are all associated with the TS1 structures having stronger overall hydrogen bonding between the carbonyl oxygen of (-)-cocaine benzoyl ester and the oxyanion hole of the enzyme. The combined computational and experimental data demonstrate a reasonable correlation relationship between the hydrogen-bonding distances in the TS1 structure and the catalytic efficiency of the enzyme against (-)-cocaine.

Characterization of a high-activity mutant of human butyrylcholinesterase against (-)-cocaine

Cocaine addiction and overdose are a well-known public health problem. There is no approved medication available for cocaine abuse treatment. Our recently designed and discovered high-activity mutant (A199S/S287G/A328W/Y332G) of human butyrylcholinesterase (BChE) has been recognized to be worth exploring for clinical application in humans as a potential anti-cocaine medication. The catalytic rate constant (k(cat)) and Michaelis-Menten constant (K(M)) for (-)-cocaine hydrolysis catalyzed by A199S/S287G/A328W/Y332G BChE (without fusion with any other peptide) have been determined to be 3,060 min(-1) and 3.1 microM, respectively, in the present study. The determined kinetic parameters reveal that the un-fused A199S/S287G/A328W/Y332G mutant has a approximately 1,080-fold improved catalytic efficiency (k(cat)/K(M)) against (-)-cocaine compared to the wild-type BChE. The approximately 1,080-fold improvement in the catalytic efficiency of the un-fused A199S/S287G/A328W/Y332G mutant is very close to the previously reported the approximately 1,000-fold improvement in the catalytic efficiency of the A199S/S287G/A328W/Y332G mutant fused with human serum albumin. These results suggest that the albumin fusion did not significantly change the catalytic efficiency of the BChE mutant while extending the plasma half-life. In addition, we have also examined the catalytic activities of the A199S/S287G/A328W/Y332G mutant against two other substrates, acetylthiocholine (ATC) and butyrylthiocholine (BTC). It has been shown that the A199S/S287G/A328W/Y332G mutations actually decreased the catalytic efficiencies of BChE against ATC and BTC, while considerably improving the catalytic efficiency of BChE against (-)-cocaine.

Free energy perturbation simulation on transition states and high-activity mutants of human butyrylcholinesterase for (-)-cocaine hydrolysis

A unified computational approach based on free energy perturbation (FEP) simulations of transition states has been employed to calculate the mutation-caused shifts of the free energy change from the free enzyme to the rate-determining transition state for (-)-cocaine hydrolysis catalyzed by the currently most promising series of mutants of human butyrylcholinesterase (BChE) that contain the A199S/A328W/Y332G mutations. The FEP simulations were followed by Michaelis-Menten kinetics analysis determining the individual k(cat) and K(M) values missing for the A199S/F227A/A328W/Y332G mutant in this series. The calculated mutation-caused shifts of the free energy change from the free enzyme to the rate-determining transition state are in good agreement with the experimental kinetic data, demonstrating that the unified computational approach based on the FEP simulations of the transition states may be valuable for future computational design of new BChE mutants with a further improved catalytic efficiency against (-)-cocaine.

Mechanism for cocaine blocking the transport of dopamine: insights from molecular modeling and dynamics simulations

Molecular modeling and dynamics simulations have been performed to study how cocaine inhibits dopamine transporter (DAT) for the transport of dopamine. The computationally determined DAT-ligand binding mode is totally different from the previously proposed overlap binding mode in which cocaine- and dopamine-binding sites are the same (Beuming, T.; et al. Nat. Neurosci. 2008, 11, 780-789). The new cocaine-binding site does not overlap with, but is close to, the dopamine-binding site. Analysis of all results reveals that when cocaine binds to DAT, the initial binding site is likely the one modeled in this study because this binding site can naturally accommodate cocaine. Then cocaine may move to the dopamine-binding site after DAT makes some necessary conformational change and expands the binding site cavity. It has been demonstrated that cocaine may inhibit the transport of dopamine through both blocking the initial DAT-dopamine binding and reducing the kinetic turnover of the transporter following the DAT-dopamine binding. The relative contributions to the phenomenological inhibition of the transport of dopamine from blocking the initial binding and reducing the kinetic turnover can be different in different types of assays. The obtained general structural and mechanistic insights are consistent with available experimental data and could be valuable for guiding future studies toward understanding cocaine's inhibiting of other transporters.

Reaction pathway and free-energy barrier for reactivation of dimethylphosphoryl-inhibited human acetylcholinesterase

The dephosphorylation/reactivation mechanism and the corresponding free-energy profile of the dimethylphosphoryl-inhibited conjugate of human acetylcholinesterase (AChE) has been studied by performing first-principles quantum mechanical/molecular mechanical free-energy (QM/MM-FE) calculations. On the basis of the QM/MM-FE results, for the favorable reaction pathway, the entire dephosphorylation/reactivation process consists of three reaction steps, including the nucleophilic water attack on the P atom, the spatial reorganization of the dimethylphosphoryl group, and the dissociation between the dimethylphosphoryl group and Ser203 of AChE. The overall free-energy barrier for the entire dephosphorylation/reactivation reaction is found to be the free-energy change from the initial reactant to the transition state associated with the spatial reorganization step, and the calculated overall free-energy barrier (20.1 to 23.5 kcal/mol) is reasonably close to the experimentally derived activation free energy of 22.3 kcal/mol. In addition, key amino acid residues and their specific roles in the reaction process have been identified.

Reaction pathway and free-energy barrier for reactivation of dimethylphosphoryl-inhibited human acetylcholinesterase

The dephosphorylation/reactivation mechanism and the corresponding free-energy profile of the dimethylphosphoryl-inhibited conjugate of human acetylcholinesterase (AChE) has been studied by performing first-principles quantum mechanical/molecular mechanical free-energy (QM/MM-FE) calculations. On the basis of the QM/MM-FE results, for the favorable reaction pathway, the entire dephosphorylation/reactivation process consists of three reaction steps, including the nucleophilic water attack on the P atom, the spatial reorganization of the dimethylphosphoryl group, and the dissociation between the dimethylphosphoryl group and Ser203 of AChE. The overall free-energy barrier for the entire dephosphorylation/reactivation reaction is found to be the free-energy change from the initial reactant to the transition state associated with the spatial reorganization step, and the calculated overall free-energy barrier (20.1 to 23.5 kcal/mol) is reasonably close to the experimentally derived activation free energy of 22.3 kcal/mol. In addition, key amino acid residues and their specific roles in the reaction process have been identified.

Fundamental reaction mechanism and free energy profile for (-)-cocaine hydrolysis catalyzed by cocaine esterase

The fundamental reaction mechanism of cocaine esterase (CocE)-catalyzed hydrolysis of (-)-cocaine and the corresponding free energy profile have been studied by performing pseudobond first-principles quantum mechanical/molecular mechanical free energy (QM/MM-FE) calculations. On the basis of the QM/MM-FE results, the entire hydrolysis reaction consists of four reaction steps, including the nucleophilic attack on the carbonyl carbon of (-)-cocaine benzoyl ester by the hydroxyl group of Ser117, dissociation of (-)-cocaine benzoyl ester, nucleophilic attack on the carbonyl carbon of (-)-cocaine benzoyl ester by water, and finally dissociation between the (-)-cocaine benzoyl group and Ser117 of CocE. The third reaction step involving the nucleophilic attack of a water molecule was found to be rate-determining, which is remarkably different from (-)-cocaine hydrolysis catalyzed by wild-type butyrylcholinesterase (BChE; where the formation of the prereactive BChE-(-)-cocaine complex is rate-determining) or its mutants containing Tyr332Gly or Tyr332Ala mutation (where the first chemical reaction step is rate-determining). Besides, the role of Asp259 in the catalytic triad of CocE does not follow the general concept of the "charge-relay system" for all serine esterases. The free energy barrier calculated for the rate-determining step of CocE-catalyzed hydrolysis of (-)-cocaine is 17.9 kcal/mol, which is in good agreement with the experimentally derived activation free energy of 16.2 kcal/mol. In the present study, where many sodium ions are present, the effects of counterions are found to be significant in determining the free energy barrier. The finding of the significant effects of counterions on the free energy barrier may also be valuable in guiding future mechanistic studies on other charged enzymes.

Model of human butyrylcholinesterase tetramer by homology modeling and dynamics simulation

A mutant of human butyrylcholinesterase (BChE) with high activity against cocaine would be highly promising as a drug for therapeutic treatment of cocaine abuse and overdose. It is desirable to design a recombinant BChE mutant with a long half-life in human circulation. Studies showed that BChE subunits can be assembled by a peptide containing the proline-rich attachment domain (PRAD) to form a stable tetramer. The models of BChE tetramer complexed with PRAD with various sequences have been constructed, in the present study, on the basis of homology modeling and molecular dynamics simulation of explicit water-solvated systems. The 3D models enable us to understand how the BChE subunits are arranged in the tetramer and how the tetramerization domain of BChE is associated with PRAD to form a stable tetramer of human BChE. It has been shown that the six conserved hydrophobic residues located on the C-terminal of BChE are responsible for the key electrostatic and hydrophobic interactions between the tetramerization domain of BChE and PRAD. The simulated tetramer structures suggest that mutation of three residues, i.e., Phe547, Met554, and Phe561, to other hydrophobic residues may be beneficial for increasing the binding between the tetramerization domain of BChE and PRAD. Thus, the detailed structural insights obtained from this study may be valuable for rational design of a recombinant BChE tetramer with a longer residence time in circulation.

Understanding the mechanism of drug resistance due to a codon deletion in protoporphyrinogen oxidase through computational modeling

Protoporphyrinogen oxidase (PPO; EC 1.3.3.4) is the last common enzyme for the enzymatic transformation of protoporphyrinogen-IX to protoporphyrin-IX, which is the key common intermediate leading to heme and chlorophyll. Hence, PPO has been identified as one of the most importance action targets for the treatment of some important diseases including cancer and variegated porphyria (VP). In the agricultural field, PPO inhibitors have been used as herbicides for many years. Recently, a unique drug resistance was found to be associated with a nonactive site residue (Gly210) deletion rather than substitution in A. tuberculatus PPO. In the present study, extensive computational simulations, including homology modeling, molecular dynamics (MD) simulations, and molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) calculations, have been carried out to uncover the detailed molecular mechanism of drug resistance associated with Gly210 deletion. Although Gly210 in the wild-type A. tuberculatus PPO has no direct interaction with the inhibitors, all the computational models and energetic results indicated that Gly210 deletion has great effects on the hydrogen-bonding network and the conformational change of the binding pocket. An interchain hydrogen bond between Gly210 with Ser424, playing an important role in stabilizing the local conformation of the wild-type enzyme, disappeared after Gly210 deletion. As a result, the mutant-type PPO has a lower affinity than the wild-type enzyme, which accounts for the molecular mechanism of drug resistance. The structural and mechanistic insights obtained from the present study provide a new starting point for future rational design of novel PPO inhibitors to overcome drug resistance associated with Gly210 deletion.

Free-energy perturbation simulation on transition states and redesign of butyrylcholinesterase

It is recognized that an ideal anti-cocaine treatment is to accelerate cocaine metabolism by producing biologically inactive metabolites via a route similar to the primary cocaine-metabolizing pathway, i.e., butyrylcholinesterase (BChE)-catalyzed hydrolysis of cocaine. BChE mutants with a higher catalytic activity against (-)-cocaine are highly desired for use as an exogenous enzyme in humans. To develop a rational design for high-activity mutants, we carried out free-energy perturbation (FEP) simulations on various mutations of the transition-state structures in addition to the corresponding free-enzyme structures by using an extended FEP procedure. The FEP simulations on the mutations of both the free-enzyme and transition-state structures allowed us to calculate the mutation-caused shift of the free-energy change from the free enzyme (BChE) to the transition state, and thus to theoretically predict the mutation-caused shift of the catalytic efficiency (k(cat)/K(M)). The computational predictions are supported by the kinetic data obtained from the wet experiments, demonstrating that the FEP-based computational design approach is promising for rational design of high-activity mutants of an enzyme. One of the BChE mutants designed and discovered in this study has an approximately 1800-fold improved catalytic efficiency against (-)-cocaine compared to wild-type BChE. The high-activity mutant may be therapeutically valuable.

QM/MM methods for biomolecular systems

Combined quantum-mechanics/molecular-mechanics (QM/MM) approaches have become the method of choice for modeling reactions in biomolecular systems. Quantum-mechanical (QM) methods are required for describing chemical reactions and other electronic processes, such as charge transfer or electronic excitation. However, QM methods are restricted to systems of up to a few hundred atoms. However, the size and conformational complexity of biopolymers calls for methods capable of treating up to several 100,000 atoms and allowing for simulations over time scales of tens of nanoseconds. This is achieved by highly efficient, force-field-based molecular mechanics (MM) methods. Thus to model large biomolecules the logical approach is to combine the two techniques and to use a QM method for the chemically active region (e.g., substrates and co-factors in an enzymatic reaction) and an MM treatment for the surroundings (e.g., protein and solvent). The resulting schemes are commonly referred to as combined or hybrid QM/MM methods. They enable the modeling of reactive biomolecular systems at a reasonable computational effort while providing the necessary accuracy.

Most efficient cocaine hydrolase designed by virtual screening of transition states

Cocaine is recognized as the most reinforcing of all drugs of abuse. There is no anticocaine medication available. The disastrous medical and social consequences of cocaine addiction have made the development of an anticocaine medication a high priority. It has been recognized that an ideal anticocaine medication is one that accelerates cocaine metabolism producing biologically inactive metabolites via a route similar to the primary cocaine-metabolizing pathway, i.e., cocaine hydrolysis catalyzed by plasma enzyme butyrylcholinesterase (BChE). However, wild-type BChE has a low catalytic efficiency against the abused cocaine. Design of a high-activity enzyme mutant is extremely challenging, particularly when the chemical reaction process is rate-determining for the enzymatic reaction. Here we report the design and discovery of a high-activity mutant of human BChE by using a novel, systematic computational design approach based on transition-state simulations and activation energy calculations. The novel computational design approach has led to discovery of the most efficient cocaine hydrolase, i.e., a human BChE mutant with an approximately 2000-fold improved catalytic efficiency, promising for therapeutic treatment of cocaine overdose and addiction as an exogenous enzyme in human. The encouraging discovery resulted from the computational design not only provides a promising anticocaine medication but also demonstrates that the novel, generally applicable computational design approach is promising for rational enzyme redesign and drug discovery.

Modeling the catalysis of anti-cocaine catalytic antibody: competing reaction pathways and free energy barriers

The competing reaction pathways and the corresponding free energy barriers for cocaine hydrolysis catalyzed by an anti-cocaine catalytic antibody, mAb15A10, were studied by using a novel computational strategy based on the binding free energy calculations on the antibody binding with cocaine and transition states. The calculated binding free energies were used to evaluate the free energy barrier shift from the cocaine hydrolysis in water to the antibody-catalyzed cocaine hydrolysis for each reaction pathway. The free energy barriers for the antibody-catalyzed cocaine hydrolysis were predicted to be the corresponding free energy barriers for the cocaine hydrolysis in water plus the calculated free energy barrier shifts. The calculated free energy barrier shift of -6.87 kcal/mol from the dominant reaction pathway of the cocaine benzoyl ester hydrolysis in water to the dominant reaction pathway of the antibody-catalyzed cocaine hydrolysis is in good agreement with the experimentally derived free energy barrier shift of -5.93 kcal/mol. The calculated mutation-caused shifts of the free energy barrier are also reasonably close to the available experimental activity data. The good agreement suggests that the protocol for calculating the free energy barrier shift from the cocaine hydrolysis in water to the antibody-catalyzed cocaine hydrolysis may be used in future rational design of possible high-activity mutants of the antibody as anti-cocaine therapeutics. The general strategy of the free energy barrier shift calculation may also be valuable in studying a variety of chemical reactions catalyzed by other antibodies or proteins through noncovalent bonding interactions with the substrates.

Structure-and-mechanism-based design and discovery of therapeutics for cocaine overdose and addiction

(-)-Cocaine is a widely abused drug and there is currently no available anti-cocaine therapeutic. Promising agents, such as anti-cocaine catalytic antibodies and high-activity mutants of human butyrylcholinesterase (BChE), for therapeutic treatment of cocaine overdose have been developed through structure-and-mechanism-based design and discovery. In particular, a unique computational design strategy based on the modeling and simulation of the rate-determining transition state has been developed and used to design and discover desirable high-activity mutants of BChE. One of the discovered high-activity mutants of BChE has a approximately 456-fold improved catalytic efficiency against (-)-cocaine. The encouraging outcome of the structure-and-mechanism-based design and discovery effort demonstrates that the unique computational design approach based on transition state modeling and simulation is promising for rational enzyme redesign and drug discovery. The general approach of the structure-and-mechanism-based design and discovery may be used to design high-activity mutants of any enzyme or catalytic antibody.

Rational design of an enzyme mutant for anti-cocaine therapeutics

(-)-Cocaine is a widely abused drug and there is no available anti-cocaine therapeutic. The disastrous medical and social consequences of cocaine addiction have made the development of an effective pharmacological treatment a high priority. An ideal anti-cocaine medication would be to accelerate (-)-cocaine metabolism producing biologically inactive metabolites. The main metabolic pathway of cocaine in body is the hydrolysis at its benzoyl ester group. Reviewed in this article is the state-of-the-art computational design of high-activity mutants of human butyrylcholinesterase (BChE) against (-)-cocaine. The computational design of BChE mutants have been based on not only the structure of the enzyme, but also the detailed catalytic mechanisms for BChE-catalyzed hydrolysis of (-)-cocaine and (+)-cocaine. Computational studies of the detailed catalytic mechanisms and the structure-and-mechanism-based computational design have been carried out through the combined use of a variety of state-of-the-art techniques of molecular modeling. By using the computational insights into the catalytic mechanisms, a recently developed unique computational design strategy based on the simulation of the rate-determining transition state has been employed to design high-activity mutants of human BChE for hydrolysis of (-)-cocaine, leading to the exciting discovery of BChE mutants with a considerably improved catalytic efficiency against (-)-cocaine. One of the discovered BChE mutants (i.e., A199S/S287G/A328W/Y332G) has a approximately 456-fold improved catalytic efficiency against (-)-cocaine. The encouraging outcome of the computational design and discovery effort demonstrates that the unique computational design approach based on the transition-state simulation is promising for rational enzyme redesign and drug discovery.

Free energy perturbation (FEP) simulation on the transition states of cocaine hydrolysis catalyzed by human butyrylcholinesterase and its mutants

A novel computational protocol based on free energy perturbation (FEP) simulations on both the free enzyme and transition state structures has been developed and tested to predict the mutation-caused shift of the free energy change from the free enzyme to the rate-determining transition state for human butyrylcholinesterase (BChE)-catalyzed hydrolysis of (-)-cocaine. The calculated shift, denoted by DeltaDeltaG(1 --> 2), of such kind of free energy change determines the catalytic efficiency (kcat/KM) change caused by the simulated mutation transforming enzyme 1 to enzyme 2. By using the FEP-based computational protocol, the DeltaDeltaG(1 --> 2) values for the mutations A328W/Y332A --> A328W/Y332G and A328W/Y332G --> A328W/Y332G/A199S were calculated to be -0.22 and -1.94 kcal/mol, respectively. The calculated DeltaDeltaG(1 --> 2) values predict that the change from the A328W/Y332A mutant to the A328W/Y332G mutant should slightly improve the catalytic efficiency and that the change from the A328W/Y332G mutant to the A328W/Y332G/A199S mutant should significantly improve the catalytic efficiency of the enzyme for the (-)-cocaine hydrolysis. The predicted catalytic efficiency increases are supported by the experimental data showing that kcat/KM = 8.5 x 10(6), 1.4 x 10(7), and 7.2 x 10(7) min(-1) M(-1) for the A328W/Y332A, A328W/Y332G, and A328W/Y332G/A199S mutants, respectively. The qualitative agreement between the computational and experimental data suggests that the FEP simulations may provide a promising protocol for rational design of high-activity mutants of an enzyme. The general computational strategy of the FEP simulation on a transition state can be used to study the effects of a mutation on the activation free energy for any enzymatic reaction.

First-principles calculation of pKa for cocaine, nicotine, neurotransmitters, and anilines in aqueous solution

The absolute pKa values of 24 representative amine compounds, including cocaine, nicotine, 10 neurotransmitters, and 12 anilines, in aqueous solution were calculated by performing first-principles electronic structure calculations that account for the solvent effects using four different solvation models, i.e., the surface and volume polarization for electrostatic interaction (SVPE) model, the standard polarizable continuum model (PCM), the integral equation formalism for the polarizable continuum model (IEFPCM), and the conductor-like screening solvation model (COSMO). Within the examined computational methods, the calculations using the SVPE model lead to the absolute pKa values with the smallest root-mean-square-deviation (rmsd) value (1.18). When the SVPE model was replaced by the PCM, IEFPCM, and COSMO, the rmsd value of the calculated absolute pKa values became 3.21, 2.72, and 3.08, respectively. All types of calculated pKa values linearly correlate with the experimental pKa values very well. With the empirical corrections using the linear correlation relationships, the theoretical pKa values are much closer to the corresponding experimental data and the rmsd values become 0.51-0.83. The smallest rmsd value (0.51) is also associated with the SVPE model. All of the results suggest that the first-principles electronic structure calculations using the SVPE model are a reliable approach to the pKa prediction for the amine compounds.

How dopamine transporter interacts with dopamine: insights from molecular modeling and simulation

By performing homology modeling, molecular docking, and molecular dynamics simulations, we have developed three-dimensional (3D) structural models of both dopamine transporter and dopamine transporter-dopamine complex in the environment of lipid bilayer and solvent water. According to the simulated structure of dopamine transporter-dopamine complex, dopamine was orientated in a hydrophobic pocket at the midpoint of the membrane. The modeled 3D structures provide some detailed structural and mechanistic insights concerning how dopamine transporter (DAT) interacts with dopamine at atomic level, extending our mechanistic understanding of the dopamine reuptake with the help of Na(+) ions. The general features of the modeled 3D structures are consistent with available experimental data. Based on the modeled structures, our calculated binding free energy (DeltaG(bind) = -6.4 kcal/mol) for dopamine binding with DAT is also reasonably close to the experimentally derived DeltaG(bind) value of -7.4 kcal/mol. Finally, a possible dopamine-entry pathway, which involves formation and breaking of the salt bridge between side chains of Arg(85) and Asp(476), is proposed based on the results obtained from the modeling and molecular dynamics simulation. The new structural and mechanistic insights obtained from this computational study are expected to stimulate future, further biochemical and pharmacological studies on the detailed structures and mechanisms of DAT and other homologous transporters.

How can (-)-epigallocatechin gallate from green tea prevent HIV-1 infection? Mechanistic insights from computational modeling and the implication for rational design of anti-HIV-1 entry inhibitors

Possible inhibitors preventing human immunodeficiency virus type 1 (HIV-1) entry into the cells are recognized as hopeful next-generation anti-HIV-1 drugs. It is highly desirable to develop a potent inhibitor blocking binding of glycoprotein CD4 of the cell with glycoprotein gp120 of HIV-1, because the gp120-CD4 binding is the initial step of HIV-1 entry into the cells. It has been recently reported that (-)-epigallocatechin gallate (EGCG) from green tea is an inhibitor blocking gp120-CD4 binding. But the inhibitory mechanism remains unknown. For understanding the inhibitory mechanism, extensive molecular docking, molecular dynamics simulations, and binding free-energy calculations have been performed in this study to predict the most favorable structures of CD4-EGCG, gp120-CD4, and gp120-CD4-EGCG binding complexes in water. The results reveal that EGCG binds with CD4 in such a way that the calculated binding affinity of gp120 with the CD4-EGCG complex is negligible. So, the favorable binding of EGCG with CD4 can effectively block gp120-CD4 binding. The calculated CD4-EGCG binding affinity (DeltaG(bind) = -5.5 kcal/mol, K(d) = 94 microM) is in excellent agreement with available experimental data suggesting IC(50) approximately 100 microM for EGCG-blocking CD4-gp120 binding. These results and insights provide a rational basis for future design of novel, more potent inhibitors to block gp120-CD4 binding.

Bioactive Permethrin/β-Cyclodextrin Inclusion Complex

Permethrin is popularly used in a variety of therapeutic areas. However, the poor water solubility of permethrin seriously limits its wider clinical applications. The present study demonstrates that solubility of permethrin in aqueous solution can considerably increase in the presence of beta-cyclodextrin (beta-CD). Extensive experimental data along with computational modeling reveal the formation of stable permethrin/beta-CD inclusion complexes, including permethrin(beta-CD) and permethrin(beta-CD)2, through hydrophobic binding. Both permethrin(beta-CD) and permethrin(beta-CD)2 complexes coexisted in aqueous solution, and the ratio of the concentration of permethrin(beta-CD) complex to that of permethrin(beta-CD)2 complex was dependent on the concentration of beta-CD. The complexation of permethrin with beta-CD significantly improved the bioavailability of permethrin and, therefore, increased the bioactivity. The significant increase of the bioactivity of permethrin in the presence of beta-CD provides an effective approach to improve the practical use of permethrin in public health and agriculture.