Role of Asp102 in the catalytic relay system of serine proteases: a theoretical study

The role of Asp102 in the catalytic relay system of serine proteases is studied theoretically by calculating the free energy profiles of the single proton-transfer reaction by the Asn102 mutant trypsin and the concerted double proton-transfer reaction (so-called the charge-relay mechanism) of the wild-type trypsin. For each reaction, the reaction free energy profile of the rate-determining step (the tetrahedral intermediate formation step) is calculated by using ab initio QM/MM electronic structure calculations combined with molecular dynamics-free energy perturbation method. In the mutant reaction, the free energy monotonically increases along the reaction path. The rate-determining step of the mutant reaction is the formation of tetrahedral intermediate complex, not the base (His57) abstraction of the proton from Ser195. In contrast to the single proton-transfer reaction of the wild-type, MD simulations of the enzyme-substrate complex show that the catalytically favorable alignment of the relay system (the hydrogen bonding network between the mutant triad, His57, Asn102, and Ser195) is rarely observed even in the presence of a substrate at the active site. In the double proton-transfer reaction, the energy barrier is observed at the proton abstraction step, which corresponds to the rate-determining step of the single proton-transfer reaction of the wild-type. Although both reaction profiles show an increase of the activation barrier by several kcals/mol, these increases have different energetic origins: a large energetic loss of the electrostatic stabilization between His57 and Asn102 in the mutant reaction, while the lack of stabilization by the protein environment in the double proton-transfer reaction. Comparing the present results with the single proton transfer of the wild-type, Asp102 is proven to play two important roles in the catalytic process. One is to stabilize the protonated His57, or ionic intermediate, formed during the acylation, and the other is to fix the configuration around the active site, which is favorable to promote the catalytic process. These two factors are closely related to each other and are indispensable for the efficient catalysis. Also the present calculations suggest the importance of the remote site interaction between His57 and Val213-Ser214 at the catalytic transition state.     

Role of the catalytic triad and oxyanion hole in acetylcholinesterase catalysis: an ab initio QM/MM study

The initial step of the acylation reaction catalyzed by acetylcholinesterase (AChE) has been studied by a combined ab initio quantum mechanical/molecular mechanical (QM/MM) approach. The reaction proceeds through the nucleophilic addition of the Ser203 O to the carbonyl C of acetylcholine, and the reaction is facilitated by simultaneous proton transfer from Ser203 to His447. The calculated potential energy barrier at the MP2(6-31+G) QM/MM level is 10.5 kcal/mol, consistent with the experimental reaction rate. The third residue of the catalytic triad, Glu334, is found to be essential in stabilizing the transition state through electrostatic interactions. The oxyanion hole, formed by peptidic NH groups from Gly121, Gly122, and Ala204, is also found to play an important role in catalysis. Our calculations indicate that, in the AChE-ACh Michaelis complex, only two hydrogen bonds are formed between the carbonyl oxygen of ACh and the peptidic NH groups of Gly121 and Gly122. As the reaction proceeds, the distance between the carbonyl oxygen of ACh and NH group of Ala204 becomes smaller, and the third hydrogen bond is formed both in the transition state and in the tetrahedral intermediate.     

Protein engineering of nitrile hydratase activity of papain: molecular dynamics study of a mutant and wild-type enzyme

The mechanism of hydrolysis of the nitrile (N-acetyl-phenylalanyl-2-amino-propionitrile, I) catalyzed by Gln19Glu mutant of papain has been studied by nanosecond molecular dynamics (MD) simulations. MD simulations of the complex of mutant enzyme with I and of mutant enzyme covalently attached to both neutral (II) and protonated (III) thioimidate intermediates were performed. An MD simulation with the wild-type enzyme.I complex was undertaken as a reference. The ion pair between protonated His159 and thiolate of Cys25 is coplanar, and the hydrogen bonding interaction S(-)(25).HD1-ND1(159) is observed throughout MD simulation of the mutant enzyme.I complex. Such a sustained hydrogen bond is absent in nitrile-bound wild-type papain due to the flexibility of the imidazole ring of His159. The nature of the residue at position 19 plays a critical role in the hydrolysis of the covalent thioimidate intermediate. When position 19 represents Glu, the imidazolium ion of His159-ND1(+).Cys25-S(-) ion pair is distant, on average, from the nitrile nitrogen of substrate I. Near attack conformers (NACs) have been identified in which His159-ImH(+) is positioned to initiate a general acid-catalyzed addition of Cys-S(-) to nitrile. Though Glu19-CO(2)H is distant from nitrile nitrogen in the mutant.I structure, MD simulations of the mutant.II covalent adduct finds Glu19-CO(2)H hydrogen bonded to the thioimide nitrogen of II. This hydrogen bonded species is much less stable than the hydrogen bonded Glu19-CO(2)(-) with mutant-bound protonated thioimidate (III). This observation supports Glu19-CO(2)H general acid catalysis of the formation of mutant.III. This is the commitment step in the Gln19Glu mutant catalysis of nitrile hydrolysis.     

New insights into the mechanism of the Schiff base hydrolysis catalyzed by type I dehydroquinate dehydratase from S. enterica: a theoretical study

The reaction pathway of Schiff base hydrolysis catalyzed by type I dehydroquinate dehydratase (DHQD) from S. enterica has been studied by performing molecular dynamics (MD) simulations and density functional theory (DFT) calculations and the corresponding potential energy profile has also been identified. On the basis of the results, the catalytic hydrolysis process for the wild-type enzyme consists of three major reaction steps, including nucleophilic attack on the carbon atom involved in the carbon-nitrogen double bond of the Schiff base intermediate by a water molecule, deprotonation of the His143 residue, and dissociation between the product and the Lys170 residue of the enzyme. The remarkable difference between this and the previously proposed reaction mechanism is that the second step here, absent in the previously proposed reaction mechanism, plays an important role in facilitating the reaction through a key proton transfer by the His143 residue, resulting in a lower energy barrier. Comparison with our recently reported results on the Schiff base formation and dehydration processes clearly shows that the Schiff base hydrolysis is rate-determining in the overall reaction catalyzed by type I DHQD, consistent with the experimental prediction, and the calculated energy barrier of ～16.0 kcal mol(-1) is in good agreement with the experimentally derived activation free energy of ～14.3 kcal mol(-1). When the imidazole group of His143 residue is missing, the Schiff base hydrolysis is initiated by a hydroxide ion in the solution, rather than a water molecule, and both the reaction mechanism and the kinetics of Schiff base hydrolysis have been remarkably changed, clearly elucidating the catalytic role of the His143 residue in the reaction. The new mechanistic insights obtained here will be valuable for the rational design of high-activity inhibitors of type I DHQD as non-toxic antimicrobials, anti-fungals, and herbicides.     

Theoretical studies on the deacylation step of serine protease catalysis in the gas phase, in solution, and in elastase

The deacylation step of serine protease catalysis is studied using DFT and ab initio QM/MM calculations combined with MD/umbrella sampling calculations. Free energies of the entire reaction are calculated in the gas phase, in a continuum solvent, and in the enzyme elastase. The calculations show that a concerted mechanism in the gas phase is replaced by a stepwise mechanism when solvent effects or an acetate ion are added to the reference system, with the tetrahedral intermediate being a shallow minimum on the free energy surface. In the enzyme, the tetrahedral intermediate is a relatively stable species ( approximately 7 kcal/mol lower in energy than the transition state), mainly due to the electrostatic effects of the oxyanion hole and Asp102. It is formed in the first step of the reaction, as a result of a proton transfer from the nucleophilic water to His57 and of an attack of the remaining hydroxyl on the ester carbonyl. This is the rate-determining step of the reaction, which requires approximately 22 kcal/mol for activation, approximately 5 kcal/mol less than the reference reaction in water. In the second stage of the reaction, only small energy barriers are detected to facilitate the proton transfer from His57 to Ser195 and the breakdown of the tetrahedral intermediate. Those are attributed mainly to a movement of Ser195 and to a rotation of the His57 side chain. During the rotation, the imidazolium ion is stabilized by a strong H-bond with Asp102, and the C(epsilon)(1)-H...O H-bond with Ser214 is replaced by one with Thr213, suggesting that a "ring-flip mechanism" is not necessary as a driving force for the reaction. The movements of His57 and Ser195 are highly correlated with rearrangements of the binding site, suggesting that product release may be implicated in the deacylation process.     

Quantum-mechanical calculation of the reactions D + FH(ν = 0, 1, 2) → DF(ν′) + H and H + FD(ν = 0, 1, 2, 3) → HF(ν′) + D on a realistic potential energy surface

Collinear (two-mathematical-dimensional (2MD)) coupled-channel quantum-mechanical calculations have been performed on the reactions D + FH(ν = 0, 1, 2) → DF(ν′) + H and H + FD(ν = 0, 1, 2, 3) → HF(ν′) + D on a potential energy surface with a 40 kcal/mole barrier to exchange. This barrier height is close to that predicted by ab initio calculations and suggested by experiments. The relative effectiveness of reagent vibrational and translational excitation to promote reaction is considered. A one-mathematical-dimensional (1MD) model for these reactions is constructed and is shown to work very well for the D + FH reaction at high temperatures, and less well for that reaction at lower temperatures as well as for the reverse H + FD reaction. Possible reasons for the breakdowns of this model are discussed.     

Molecular Dynamics Simulations and Steered Molecular Dynamics Simulations of Glabridin Bound to Wild Type and V30A Mutant Transthyretin: Ligand-linked Perturbation of Tertiary Conformation

Transthyretin(TTR), as a tetrameric protein, functions as a neuroprotector. The native TTR homotetramer dissociates into dimers and monomers. Dimers and monomers self-assemble into amyloid fibrils, and this process can lead to some diseases. Native TTR homotetramer is a widely accepted model for TTR amyloid formation. In this study, simulations using molecular dynamics(MD) and steered MD(SMD) were performed to explore the mechanisms for glabridin(Glab), a specific inhibitor for TTR binding, for V30A mutant and wild-type(WT) TTR. MD simulation results indicate that, compared with Glab binding to WT and V30A mutant, the WT TTR could lead to the collapse of β-strands from Ser52 to His56 at chain A. This phenomenon facilitated the easy dissociation of chains A and C. Calculations of the binding free energy between the two chains showed that the V30A-Glab TTR complex displayed a lower binding energy than other systems(WT TTR and WT-Glab TTR). Then, SMD simulation was performed to explore the unbinding pathway for Glab through the WT and V30A mutant TTR. The results show that Lys15(chain A) produced a hydrogen bond with Glab at the force peak via the WT TTR tunnel. Meanwhile, in the V30A TTR mutant, the hydrogen bond between Lys15(chain A) and Glab was broken at the force peak. This condition was beneficial for Glab to be taken off from the protein. Our theoretical results will be useful in designing a new specific inhibitor of TTR protein to control the TTR homotetramer dissociation.     

The Structural Basis of Mycobacterium tuberculosis RpoB Drug-Resistant Clinical Mutations on Rifampicin Drug Binding

Tuberculosis (TB), caused by the Mycobacterium tuberculosis infection, continues to be a leading cause of morbidity and mortality in developing countries. Resistance to the first-line anti-TB drugs, isoniazid (INH) and rifampicin (RIF), is a major drawback to effective TB treatment. Genetic mutations in the β-subunit of the DNA-directed RNA polymerase (rpoB) are reported to be a major reason of RIF resistance. However, the structural basis and mechanisms of these resistant mutations are insufficiently understood. In the present study, thirty drug-resistant mutants of rpoB were initially modeled and screened against RIF via a comparative molecular docking analysis with the wild-type (WT) model. These analyses prioritized six mutants (Asp441Val, Ser456Trp, Ser456Gln, Arg454Gln, His451Gly, and His451Pro) that showed adverse binding affinities, molecular interactions, and RIF binding hinderance properties, with respect to the WT. These mutant models were subsequently analyzed by molecular dynamics (MD) simulations. One-hundred nanosecond all-atom MD simulations, binding free energy calculations, and a dynamic residue network analysis (DRN) were employed to exhaustively assess the impact of mutations on RIF binding dynamics. Considering the global structural motions and protein–ligand binding affinities, the Asp441Val, Ser456Gln, and His454Pro mutations generally yielded detrimental effects on RIF binding. Locally, we found that the electrostatic contributions to binding, particularly by Arg454 and Glu487, might be adjusted to counteract resistance. The DRN analysis revealed that all mutations mostly distorted the communication values of the critical hubs and may, therefore, confer conformational changes in rpoB to perturb RIF binding. In principle, the approach combined fundamental molecular modeling tools for robust “global” and “local” level analyses of structural dynamics, making it well suited for investigating other similar drug resistance cases.     

Theoretical perspectives on the reaction mechanism of serine proteases: the reaction free energy profiles of the acylation process

The reaction mechanism of serine proteases (trypsin), which catalyze peptide hydrolysis, is studied theoretically by ab initio QM/MM electronic structure calculations combined with Molecular Dynamics-Free Energy Perturbation calculations. We have calculated the entire reaction free energy profiles of the first reaction step of this enzyme (acylation process). The present calculations show that the rate-determining step of the acylation is the formation of the tetrahedral intermediate, and the breakdown of this intermediate has a small energy barrier. The calculated activation free energy for the acylation is approximately 17.8 kcal/mol at QM/MM MP2/(aug)-cc-pVDZ//HF/6-31(+)G/AMBER level, and this reaction is an exothermic process. MD simulations of the enzyme-substrate (ES) complex and the free enzyme in aqueous phase show that the substrate binding induces slight conformational changes around the active site, which favor the alignment of the reactive fragments (His57, Asp102, and Ser195) together in a reactive orientation. It is also shown that the proton transfer from Ser195 to His57 and the nucleophilic attack of Ser195 to the carbonyl carbon of the scissile bond of the substrate occur in a concerted manner. In this reaction, protein environment plays a crucial role to lowering the activation free energy by stabilizing the tetrahedral intermediate compared to the ES complex. The polarization energy calculations show that the enzyme active site is in a very polar environment because of the polar main chain contributions of protein. Also, the ground-state destabilization effect (steric strain) is not a major catalytic factor. The most important catalytic factor of stabilizing the tetrahedral intermediate is the electrostatic interaction between the active site and particular regions of protein: the main chain NH groups in Gly193 and Ser195 (so-called oxyanion hole region) stabilize negative charge generated on the carbonyl oxygen of the scissile bond, and the main chain carbonyl groups in Ile212 approximately Ser214 stabilize a positive charge generated on the imidazole ring of His57.     

DFT investigation on the mechanism of the deacetylation reaction catalyzed by LpxC

LpxC is a key enzyme in the biochemical synthesis of Lipid A, an important outer cell-membrane component found in a number of pathogenic bacteria. Using DFT, we have investigated the binding of the substrate within its active site as well as the deacetylation mechanism it catalyzes. The substrate is found to preferentially coordinate to the active site Zn2+ via its carbonyl oxygen between a Zn2+-bound H2O and an adjacent threonine (Thr191). Furthermore, upon substrate binding a nearby Glu78 residue is found to readily deprotonate the remaining Zn2+-bound H2O. Unlike several related metallopeptidases, the mechanism of LpxC is found to proceed via four steps; (i) initial hydroxylation of the substrates' carbonyl carbon to give a gem-diolate intermediate, (ii) protonation of the amide nitrogen by the histidine His265-H+, (iii) a barrier-less change in the active site-intermediate hydrogen-bond network and finally, (iv) C-N bond cleavage. Notably, the rate-determining step of the mechanism of LpxC is found to be the initial hydroxylation while the final C-N bond cleavage occurs with an overall barrier of 23.6 kJ mol-1. Furthermore, LpxC uses a general acid/base pair mechanism as indicated by the fact that both His265-H+ and Glu78 are accordingly involved.     

Influence of metal cofactors and water on the catalytic mechanism of creatininase-creatinine in aqueous solution from molecular dynamics simulation and quantum study

The reaction mechanism of creatinine-creatininase binding to form creatine as a final product has been investigated by using a combined ab initio quantum mechanical/molecular mechanical approach and classical molecular dynamics (MD) simulations. In MD simulations, an X-ray crystal structure of the creatininase/creatinine was modified for creatininase/creatinine complexes and the MD simulations were run for free creatininase and creatinine in water. MD results reveal that two X-ray water molecules can be retained in the active site as catalytic water. The binding free energy from Molecular Mechanics Poisson-Boltzmann Surface Area calculation predicted the strong binding of creatinine with Zn²⁺, Asp45 and Glu183. Two step mechanisms via Mn²⁺/Zn²⁺ (as in X-ray structure) and Zn²⁺/Zn²⁺ were proposed for water adding step and ring opening step with two catalytic waters. The pathway using synchronous transit methods with local density approximations with PWC functional for the fragment in the active region were obtained. Preferable pathway Zn²⁺/Zn²⁺ was observed due to lower activation energy in water adding step. The calculated energy in the second step for both systems were comparable with the barrier of 26.03 and 24.44 kcal/mol for Mn²⁺/Zn²⁺ and Zn²⁺/Zn²⁺, respectively.     

Computational mutation analysis of hydrogen abstraction and radical rearrangement steps in the catalysis of coenzyme B12-dependent diol dehydratase

A mutation analysis of the catalytic functions of active-site residues of coenzyme B(12)-dependent diol dehydratase in the conversion of 1,2-propanediol to 1,1-propanediol has been carried out by using QM/MM computations. Mutants His143Ala, Glu170Gln, Glu170Ala, and Glu170Ala/Glu221Ala were considered to estimate the impact of the mutations of His143 and Glu170. In the His143Ala mutant the activation energy for OH migration increased to 16.4 from 11.5 kcal mol(-1) in the wild-type enzyme. The highest activation energy, 19.6 kcal mol(-1), was measured for hydrogen back-abstraction in this reaction. The transition state for OH migration is not sufficiently stabilized by the hydrogen-bonding interaction formed between the spectator OH group and Gln170 in the Glu170Gln mutant, which demonstrates that a strong proton acceptor is required to promote OH migration. In the Glu170Ala mutant, a new strong hydrogen bond is formed between the spectator OH group and Glu221. A computed activation energy of 13.6 kcal mol(-1) for OH migration in the Glu170Ala mutant is only 2.1 kcal mol(-1) higher than the corresponding barrier in the wild-type enzyme. Despite the low activation barrier, the Glu170Ala mutant is inactive because the subsequent hydrogen back-abstraction is energetically demanding in this mutant. OH migration is not feasible in the Glu170Ala/Glu221Ala mutant because the activation barrier for OH migration is greatly increased by the loss of COO(-) groups near the spectator OH group. This result indicates that the effect of partial deprotonation of the spectator OH group is the most important factor in reducing the activation barrier for OH migration in the conversion of 1,2-propanediol to 1,1-propanediol catalyzed by diol dehydratase.     

Wild type and mutant forms of recombinant horseradish peroxidase C expressed inEscherichia coli

Two horseradish peroxidase C (HRPC) mutants with substitutions in the active center, i.e., Phe41→ His and Phel43→ Glu, were compared to the wild-type recombinant enzyme expressed in Escherichia coli in terms of the enzymatic activity and stability under irradiation. Both mutations caused a significant decrease in activity, but it was still possible to follow the effect of mutations on the key steps of the reaction mechanism. Phe41 can be considered a nonpolar barrier separating histidine residues in the active center and providing a firm noncovalent binding with the highly hydrophobic porphyrin ring. The replacement of Phe41 with the ionizable His residue destabilizes the enzyme. The Phel43→ Glu replacement creates a negative charge at the entrance of the heme-binding pocket, and protects the latter from both donor substrates and free radicals. The radiolytic inactivation of the wild-type and mutant forms of recombinant HRP suggested different binding sites for iodide, 2,2′-bis(3-ethylbenzothiasoline-6-sulfonate (ABTS), guaiacol, and o-phenylene diamine. The study of kinetics and inactivation is in agreement with the direct binding of iodide to the heme porphyrin ring. The results also suggest that the ABTS binding site is less accessible than that for o-phenylene diamine.     

On the formation of radical dications of protonated amino acids in a "microsolution" of water or acetonitrile and their reactivity towards the solvent

In high-energy collisions (50 keV) between O2 and protonated amino acids AH+, radical dications AH2+* are formed for A = Phe, His, Met, Tyr, and Trp. When solvated by water or acetonitrile (S), AH2+*(S)1,2 are formed for A = Arg, His, Met, Tyr, and Trp. The stability of the hydrogen-deficient AH2+* in the "microsolution" depends on the energetics of the electron transfer reaction AH2+* +S --> AH++S+*, the hydrogen abstraction reaction AH2+*+S --> AH2(2+)+[S-H]*, and the proton transfer reaction AH2+* + S --> A+*+SH+. Using B3LYP/ 6-311+G(2d,p)//B3LYP/6-31+G(d) model chemistry, we describe these three reactions in detail for A=Tyr and find that the first two reactions are unfavorable whereas the third one is favorable. However, energy is required for the formation of Tyr+* and SH+ from TyrH2+*(S) to overcome the Coulomb barrier, which renders the complex observable with a life-time larger than 5 micros. The ionization energy, IE, of TyrH+ is calculated to be 11.1 eV in agreement with an experimental measurement of 10.1+/-2.1 eV ([IE(CH3CN)+IE(Tyr)]/ 2); hydration further lowers the IE by 0.3 eV [IE(TyrH+(H2O) = 10.8 eV, calculated]. We estimate the ionization energies of TrpH+, HisH+, and MetH+ to be 10.1+/-2.1 eV, 12.4+/-0.2 eV, and 12.4+/-0.2 eV, and that of PheH+ to be larger than 12.6 eV.     

Mechanism of reaction of hydrogen peroxide with horseradish peroxidase: identification of intermediates in the catalytic cycle.

The mechanism of the reaction of horseradish peroxidase isoenzyme C (HRPC) with hydrogen peroxide to form the reactive enzyme intermediate compound I has been studied using electronic absorbance, rapid-scan stopped-flow, and electron paramagnetic resonance (EPR) spectroscopies at both acid and basic pH. The roles of the active site residues His42 and Arg38 in controlling heterolytic cleavage of the H(2)O(2) oxygen-oxygen bond have been probed with site-directed mutant enzymes His42 --> Leu (H42L), Arg38 --> Leu (R38L), and Arg38 --> Gly (R38G). The biphasic reaction kinetics of H42L with H(2)O(2) suggested the presence of an intermediate species and, at acid pH, a reversible second step, probably due to a neutral enzyme-H(2)O(2) complex and the ferric-peroxoanion-containing compound 0. EPR also indicated the formation of a protein radical situated more than approximately 10 A from the heme iron. The stoichiometry of the reaction of the H42L/H(2)O(2) reaction product and 2,2'-azinobis(3-ethylbenzothiazolinesulfonic acid) (ABTS) was concentration dependent and fell from a value of 2 to 1 above 0.7 mM ABTS. These data can be explained if H(2)O(2) undergoes homolytic cleavage in H42L. The apparent rate of compound I formation by H42L, while low, was pH independent in contrast to wild-type HRPC where the rate falls at acid pH, indicating the involvement of an ionizable group with pK(a) approximately 4. In R38L and R38G, the apparent pK(a) was shifted to approximately 8 but there is no evidence that homolytic cleavage of H(2)O(2) occurs. These data suggest that His42 acts initially as a proton acceptor (base catalyst) and then as a donor (acid catalyst) at neutral pH and predict the observed slower rate and lower efficiency of heterolytic cleavage observed at acid pH. Arg38 is influential in lowering the pK(a) of His42 and additionally in aligning H(2)O(2) in the active site, but it does not play a direct role in proton transfer.     

An experimental and theoretical study of NSCl decomposition in the presence of trace amounts of water

An experimental and density functional theory study of the decomposition reaction of NSCl with water molecule(s) into NH 3 and SO 2 is described. Experiments involving the decomposition of NSCl in the presence of trace amounts of water showed that the formation of both HNSO and SO 2 is faster when the cell and the cell surface are fresh and more likely to be coated in adsorbed water. The density functional theory computational results suggest that the decomposition involves two main steps: reactions 1 NSCl + nH 2O --> HNSO + ( n - 1) H 2O + HCl and 2 HNSO + nH 2O --> NH 3 + SO 2 + ( n - 1)H 2O. The barrier for reaction 1 in which NSCl decomposes to form HNSO becomes substantially lower when the explicit hydrogen bonding of water molecules is considered. The water assisted decomposition reaction mechanism can help account for the experimentally observed conversion of NSCl to produce HNSO in the presence of water. The calculated barrier of reaction 2 (the decomposition of HNSO into NH 3 and SO 2) can also explain the experimental observation of SO 2 bands from the decomposition of NSCl in the presence of water.     

Density functional calculations of ATP systems. 2. ATP hydrolysis at the active site of actin

The hydrolysis of adenosine 5'-triphosphate (ATP) at the active site of actin has been studied using density functional calculations. The active site is modeled by the triphosphate tail of ATP, an Mg cation, surrounding water molecules, and the nearby protein residues. Four reaction paths have been followed by constraining coordinates that allow phosphate stretching, nucleophilic attack of the catalytic water, and OH(-) formation via water deprotonation. The lowest-energy barrier (21.0 kcal/mol) is obtained for a dissociative reaction where the terminal phosphate breaks on approaching the catalytic water, followed by proton release via a proton wire mechanism. A higher barrier (39.6 kcal/mol) results for an associative reaction path where OH(-) is formed first, with a pentacoordinated phosphorus atom (P-O distances 2.1 A). Stretching the terminal bridging P-O bond results in bond rupture at 2.8 A with an energy barrier of 28.8 kcal/mol. The residues Gln137 and His161 are not important in the reactions, but insight into their roles in vivo has been obtained. The favored coordination of the end products H(2)PO(4)(-) and ADP(3-) includes a hydrogen bond and an O-Mg-O bridge between the phosphates as well as a hydrogen bond between H(2)PO(4)(-) and the Ser14 side chain. The total energy is 2.1 kcal/mol lower than in the initial reactants. Classical simulations of ATP- and ADP.P(i)-actin show few hydrolysis-induced differences in the protein structure, indicating that phosphate migration is necessary for a change in conformation.     

Molecular dynamic and free energy studies of primary resistance mutations in HIV-1 protease-ritonavir complexes

To understand the basis of drug resistance of the HIV-1 protease, molecular dynamic (MD) and free energy calculations of the wild-type and three primary resistance mutants, V82F, I84V, and V82F/I84V, of HIV-1 protease complexed with ritonavir were carried out. Analysis of the MD trajectories revealed overall structures of the protein and the hydrogen bonding of the catalytic residues to ritonavir were similar in all four complexes. Substantial differences were also found near the catalytic binding domain, of which the double mutant complex has the greatest impact on conformational changes of the protein and the inhibitor. The tip of the HIV-1 protease flap of the double mutant has the greater degree of opening with respect to that of the others. Additionally, the phenyl ring of Phe82 moves away from the binding pocket S1', and the conformational change of ritonavir subsite P1' consequently affects the cavity size of the protein and the conformational energy of the inhibitor. Calculations of binding free energy using the solvent continuum model were able to reproduce the same trend of the experimental inhibition constant. The results show that the resistance mutants require hydrophobic residues to maintain the interactions in the binding pocket. Changes of the cavity volume correlate well with free energy penalties due to the mutation and are responsible for the loss of drug susceptibility.     

Molecular dynamics simulation of cocaine binding with human butyrylcholinesterase and its mutants

Molecular dynamics (MD) simulations were carried out to study cocaine binding with wild-type human butyrylcholinesterase (BChE) and its mutants based on a recently reported X-ray crystal structure of human BChE. For each BChE-cocaine system, we simulated both the nonprereactive and prereactive complexes in water. Despite the significant difference found at the acyl binding pocket, the simulated structures confirm the fundamental structural and mechanistic insights obtained from earlier computational studies of wild-type BChE with cocaine based on a homology model, e.g. the rate-determining step for BChE-catalyzed hydrolysis of biologically active (-)-cocaine is the (-)-cocaine rotation in the active site from the nonprereactive BChE-(-)-cocaine complex to the prereactive complex. It has been demonstrated that the MD simulations on both the nonprereactive and prereactive BChE-cocaine complexes can clearly reveal whether specific mutations produce the desired BChE-(-)-cocaine binding structures in which the (-)-cocaine rotation is less hindered while the required prereactive BChE-(-)-cocaine binding is maintained. Based on the MD simulations, both A328W/Y332A and A328W/Y332G BChE's are expected to have catalytic activity for (-)-cocaine hydrolysis higher than that of wild-type BChE and the activity of A328W/Y332G BChE should be slightly higher than that of A328W/Y332A BChE due to the less-hindered (-)-cocaine rotation in the mutant BChE's. However, the less-hindered (-)-cocaine rotation is only a necessary condition for a higher activity mutant BChE. The (-)-cocaine rotation is also less hindered in A328W/Y332A/Y419S BChE, but (-)-cocaine binds with A328W/Y332A/Y419S BChE in a way that is not suitable for the catalysis. Thus, A328W/Y332A/Y419S BChE is expected to lose the catalytic activity. The computational predictions were confirmed by our experimental kinetic data, demonstrating that the MD simulation-based computational protocol used in this study is reliable in prediction of the catalytic activity of BChE mutants for (-)-cocaine hydrolysis.