Molecular dynamics simulation of Axillaridine–A: A potent natural cholinesterase inhibitor

Molecular Dynamics (MD) simulations were carried out for human acetylcholinesterase (hAChE) and its complex with Axillaridine–A, in order to dynamically explore the active site of the protein and the behaviour of the ligand at the peripheral binding site. Simulation of the enzyme alone showed that the active site of AChE is located at the bottom of a deep and narrow cavity whose surface is lined with rings of aromatic residues while Tyr72 is almost perpendicular to the Trp286, which is responsible for stable π -π interactions. The complexation of AChE with Axillaridine-A, results in the reduction of gorge size due to interaction between the ligand and the active site residues. The gorge size was determined by the distance between the center of mass of Glu81 and Trp286. As far as the geometry of the active site is concerned, the presence of ligand in the active site alters its specific conformation, as revealed by stable hydrogen bondings established between amino acids. With the increasing interaction between ligand and the active amino acids, size of the active site of the complex decreases with respect to time. Axillaridine-A, forms stable π -π interactions with the aromatic ring of Tyr124 that results in inhibition of catalytic activity of the enzyme. This π -π interaction keeps the substrate stable at the edge of the catalytic gorge by inhibiting its catalytic activity. The MD results clearly provide an explanation for the binding pattern of bulky steroidal alkaloids at the active site of AChE.     

A computational study of the deacylation mechanism of human butyrylcholinesterase

To investigate the mechanism of the deacylation reaction in the active site of human butyrylcholinesterase (BuChE), we carried out quantum mechanical (QM) calculations on cluster models of the active site built from a crystallographic structure. The models consisted of the substrate butyrate moiety, the catalytic triad of residues (Ser198, Glu325, and His438), the "oxy-anion hole" (Gly116, Gly117, and Ala199), the side chain of Glu197, four water molecules, the side chain of Ser225, and the peptide linkage between Val321 and Asn322. Analyses of the equilibrium geometries, electronic properties, and energies of the QM models gave insights into the catalytic mechanism. In addition, the QM calculations provided the data required to build a molecular mechanics representation of the reactive BuChE region that was employed in molecular dynamics simulations followed by molecular-mechanics-Poisson-Boltzmann (MM-PB) calculations. Subsequently, we combined the QM energies with average MM-PB energies to estimate the free energy of the reactive structures in the enzyme. The rate-determining step corresponds to the formation of a tetrahedral intermediate with a free-energy barrier of approximately 14.0 kcal/mol. The modulation of the BuChE activity, exerted by either neutral molecules (glycerol, GOL) or a second butyrylcholine (CHO) molecule bound to the cation-pi site, does not involve any significant allosteric effect. Interestingly, the presence of GOL or CHO stabilizes a product complex formed between a butyric acid molecule and BuChE. These results are in consonance with the crystallographic structure of BuChE, in which the catalytic Ser198 interacts with a butyric fragment, while the cation-pi site is occupied by one GOL molecule.     

Kinetic analysis of butyrylcholinesterase-catalyzed hydrolysis of acetanilides

The aryl-acylamidase (AAA) activity of butyrylcholinesterase (BuChE) has been known for a long time. However, the kinetic mechanism of aryl-acylamide hydrolysis by BuChE has not been investigated. Therefore, the catalytic properties of human BuChE and its peripheral site mutant (D70G) toward neutral and charged aryl-acylamides were determined. Three neutral (o-nitroacetanilide, m-nitroacetanilide, o-nitrophenyltrifluoroacetamide) and one positively charged (3-(acetamido) N,N,N-trimethylanilinium, ATMA) acetanilides were studied. Hydrolysis of ATMA by wild-type and D70G enzymes showed a long transient phase preceding the steady state. The induction phase was characterized by a hysteretic "burst". This reflects the existence of two enzyme states in slow equilibrium with different catalytic properties. Steady-state parameters for hydrolysis of the three acetanilides were compared to catalytic parameters for hydrolysis of esters giving the same acetyl intermediate. Wild-type BuChE showed substrate activation while D70G displayed a Michaelian behavior with ATMA as with positively charged esters. Owing to the low affinity of BuChE for amide substrates, the hydrolysis kinetics of neutral amides was first order. Acylation was the rate-determining step for hydrolysis of aryl-acetylamide substrates. Slow acylation of the enzyme, relative to that by esters may, in part, be due suboptimal fit of the aryl-acylamides in the active center of BuChE. The hypothesis that AAA and esterase active sites of BuChE are non-identical was tested with mutant BuChE. It was found that mutations on the catalytic serine, S198C and S198D, led to complete loss of both activities. The silent variant (FS117) had neither esterase nor AAA activity. Mutation in the peripheral site (D70G) had the same effect on esterase and AAA activities. Echothiophate inhibited both activities identically. It was concluded that the active sites for esterase and AAA activities are identical, i.e. S198. This excludes any other residue present in the gorge for being the catalytic nucleophile pole.     

Structure and dynamics of the active site gorge of acetylcholinesterase: synergistic use of molecular dynamics simulation and X-ray crystallography.

The active site of acetylcholinesterase (AChE) from Torpedo californica is located 20 A from the enzyme surface at the bottom of a narrow gorge. To understand the role of this gorge in the function of AChE, we have studied simulations of its molecular dynamics. When simulations were conducted with pure water filling the gorge, residues in the vicinity of the active site deviated quickly and markedly from the crystal structure. Further study of the original crystallographic data suggests that a bis-quaternary decamethonium (DECA) ion, acquired during enzyme purification, residues in the gorge. There is additional electron density within the gorge that may represent small bound cations. When DECA and 2 cations are placed within the gorge, the simulation and the crystal structure are dramatically reconciled. The small cations, more so than DECA, appear to stabilize part of the gorge wall through electrostatic interactions. This part of the gorge wall is relatively thin and may regulate substrate, product, and water movement through the active site.     

Differential effect of pressure and temperature on the catalytic behaviour of wild-type human butyrylcholinesterase and its D70G mutant.

The combined action of temperature (10-35 degrees C) and pressure (0. 001-2 kbar) on the catalytic activity of wild-type human butyrylcholinesterase (BuChE) and its D70G mutant was investigated at pH 7.0 using butyrylthiocholine as the substrate. The residue D70, located at the mouth of the active site gorge, is an essential component of the peripheral substrate binding site of BuChE. Results showed a break in Arrhenius plots of wild-type BuChE (at Tt approximately 22 degrees C) whatever the pressure (dTt/dP = 1.6 +/- 1.5 degrees C.kbar-1), whereas no break was observed in Arrhenius plots of the D70G mutant. These results suggested a temperature-induced conformational change of the wild-type BuChE which did not occur for the D70G mutant. For the wild-type BuChE, at around a pressure of 1 kbar, an intermediate state, whose affinity for substrate was increased, appeared. This intermediate state was not seen for the mutant enzyme. The wild-type BuChE remained active up to a pressure of 2 kbar whatever the temperature, whereas the D70G mutant was found to be more sensitive to pressure inactivation (at pressures higher than 1.5 kbar the mutant enzyme lost its activity at temperatures lower than 25 degrees C). The results indicate that the residue D70 controls the conformational plasticity of the active site gorge of BuChE, and is involved in regulation of the catalytic activity as a function of temperature.     

Structure-activity relationships for inhibition of human cholinesterases by alkyl amide phenothiazine derivatives

Several lines of evidence indicate that inhibition of butyrylcholinesterase (BuChE) is important in the treatment of certain dementias. Further testing of this concept requires inhibitors that are both BuChE-selective and robust. N-alkyl derivatives (2, 3, 4) of phenothiazine (1) have previously been found to inhibit only BuChE in a mechanism involving pi-pi interaction between the phenothiazine tricyclic ring system and aromatic residues in the active site gorge. To explore features of phenothiazines that affect the selectivity and potency of BuChE inhibition, a series of N-carbonyl derivatives (5-25) was synthesized and examined for the ability to inhibit cholinesterases. Some of the synthesized derivatives also inhibited AChE through a different mechanism involving carbonyl interaction within the active site gorge. Binding of these derivatives takes place within the gorge, since this inhibition disappears when the molecular volume of the derivative exceeds the estimated active site gorge volume of this enzyme. In contrast, BuChE, with a much larger active site gorge, exhibited inhibition that increased directly with the molecular volumes of the derivatives. This study describes two distinct mechanisms for binding phenothiazine amide derivatives to BuChE and AChE. Molecular volume was found to be an important parameter for BuChE-specific inhibition.     

Butyrylcholinesterase-catalyzed hydrolysis of N-methylindoxyl acetate: analysis of volume changes upon reaction and hysteretic behavior

Hydrolysis of the neutral substrate N-methylindoxyl acetate (NMIA) by wild-type human butyrylcholinesterase (BuChE) and peripheral site mutants (D70G, Y332A, D70G/Y332A) was found to follow the Michaelis-Menten kinetics. K(m) was 0.14 mM for wild-type, and 0.07-0.16 mM for D70G, Y332A and D70G/Y332A, indicating that the peripheral site is not involved in NMIA binding. The values of k(cat) were of the same order for all enzymes: 12,000-18,000 min(-1). Volume changes upon substrate binding (-DeltaV(K(m))) and the activation volumes (DeltaV++(k(cat)) associated with hydrolysis of NMIA were calculated from the pressure dependence of the catalytic constants. Values of -DeltaV(K(m)) indicate that NMIA binds to an aromatic residue, presumed to be W82, the active site binding locus. Binding is accompanied by a release of water molecules from the gorge. Residue 70 controls the number of water molecules that are released upon substrate binding. The values of DeltaV++(k(cat)), which are positive for wild-type and faintly positive for D70G, clearly indicate that the catalytic steps are accompanied by re-entry of water into the gorge. Results support the premise that residue D70 is involved in the conformational stabilization of the active site gorge and in control of its hydration. A slow transient, preceding the steady state, was seen on a time scale of several minutes. The induction time rapidly increased with NMIA concentration to reach a limit at substrate saturation. Much shorter induction times (<1 min) were seen for hydrolysis of benzoylcholine (BzCh) by wild-type BuChE and for hydrolysis of butyrylthiocholine (BuSCh) by the active site mutants E197Q and E197Q/G117H. This slow transient was interpreted in terms of hysteresis without kinetic cooperativity. The hysteretic behavior of BuChE results from a slow conformational equilibrium between two enzyme states E and E'. NMIA binds only to the primed form E'. Kosmotropic salts and hydrostatic pressure were found to shift the equilibrium toward E'. The E-->E' transition is accompanied by a negative activation volume (DeltaV++(0)= -45+/-10 ml/mol), and the E' form is more compact than E. Hydration water in the gorge of E' appears to be more structured than in the unprimed form.     

Analysis of the activation of acetylcholinesterase by carbon nanoparticles using a monolithic immobilized enzyme microreactor: role of the water molecules in the active site gorge

A biochromatographic system was used to study the direct effect of carbon nanoparticles (CNPs) on the acetylcholinesterase (AChE) activity. The AChE enzyme was covalently immobilized on a monolithic CIM-disk via its NH₂ residues. Our results showed an increase in the AChE activity in presence of CNPs. The catalytic constant (k_cat) was increased while the Michaelis constant (K_m) was slightly decreased. This indicated an increase in the enzyme efficiency with increase of the substrate affinity to the active site. The thermodynamic data of the activation mechanism of the enzyme, i.e. ΔH* and ΔS*, showed no change in the substrate interaction mechanism with the anionic binding site. The increase of the enthalpy (ΔH*) and the entropy (ΔS*) with decrease in the free energy of activation (Ea) was related to structural conformation change in the active site gorge. This affected the stability of water molecules in the active site gorge and facilitated water displacement by substrate for entering to the active site of the enzyme.     

Novel isosorbide di-ester compounds as inhibitors of acetylcholinesterase

We report herein that a variety of isosorbide di-esters, previously reported to be novel substrates for butyrylcholinesterase (BuChE, EC 3.1.1.8), are in fact inhibitors of the homologous enzyme acetylcholinesterase (AChE), with IC(50) values in the micromolar range. In vitro studies show that they are mixed inhibitors of the enzyme, and thus the ternary enzyme-inhibitor-substrate complex can form in acetylcholinesterase. This is rationalised by molecular modelling which shows that the compounds bind in the mid-gorge area. In this position, simultaneous substrate binding might be possible, but the hydrolysis of this substrate is prevented. The di-esters dock within the butyrylcholinesterase gorge in a very different manner, with the ester sidechain at the 5-position occupying the acyl pocket at residues Leu286 and Val288, and the 2-ester binding to Trp82. The carbonyl group of the 2-ester is susceptible to nucleophilic attack by Ser198 of the catalytic triad. The larger residues of the acyl pocket in acetylcholinesterase prevent binding in this manner. The results complement each other and explain the differing behaviours of the esters in the cholinesterase enzymes. These findings may prove very significant for future work.     

Nanosecond dynamics of acetylcholinesterase near the active center gorge

To delineate the role of peptide backbone flexibility and rapid molecular motion in acetylcholinesterase catalysis and inhibitor association, we investigated the decay of fluorescence anisotropy at three sites of fluorescein conjugation to cysteine-substitution mutants of the enzyme. One cysteine was placed in a loop at the peripheral site near the rim of the active center gorge (H287C); a second was in a helical region outside of the active center gorge (T249C); a third was at the tip of a small, flexible omega loop well separated from the gorge (A262C). Mutation and fluorophore conjugation did not appreciably alter catalytic or inhibitor binding parameters of the enzyme. The results show that each site examined was associated with a high degree of segmental motion; however, the A262C and H287C sites were significantly more flexible than the T249C site. Association of the active center inhibitor, tacrine, and the peripheral site peptide inhibitor, fasciculin, had no effect on the anisotropy decay of fluorophores at positions 249 and 262. Fasciculin, but not tacrine, on the other hand, dramatically altered the decay profile of the fluorophore at the 287 position, in a manner consistent with fasciculin reducing the segmental motion of the peptide chain in this local region. The results suggest that the motions of residues near the active center gorge and across from the Cys(69)-Cys(96) omega loop are uncoupled and that ligand binding at the active center or the peripheral site does not influence acetylcholinesterase conformational dynamics globally, but induces primarily domain localized decreases in flexibility proximal to the bound ligand.     

How an Inhibitor Bound to Subunit Interface Alters Triosephosphate Isomerase Dynamics

The tunnel region at triosephosphate isomerase (TIM)’s dimer interface, distant from its catalytic site, is a target site for certain benzothiazole derivatives that inhibit TIM’s catalytic activity in Trypanosoma cruzi, the parasite that causes Chagas disease. We performed multiple 100-ns molecular-dynamics (MD) simulations and elastic network modeling (ENM) on both apo and complex structures to shed light on the still unclear inhibitory mechanism of one such inhibitor, named bt10. Within the time frame of our MD simulations, we observed stabilization of aromatic clusters at the dimer interface and enhancement of intersubunit hydrogen bonds in the presence of bt10, which point to an allosteric effect rather than destabilization of the dimeric structure. The collective dynamics dictated by the topology of TIM is known to facilitate the closure of its catalytic loop over the active site that is critical for substrate entrance and product release. We incorporated the ligand’s effect on vibrational dynamics by applying mixed coarse-grained ENM to each one of 54,000 MD snapshots. Using this computationally efficient technique, we observed altered collective modes and positive shifts in eigenvalues due to the constraining effect of bt10 binding. Accordingly, we observed allosteric changes in the catalytic loop’s dynamics, flexibility, and correlations, as well as the solvent exposure of catalytic residues. A newly (to our knowledge) introduced technique that performs residue-based ENM scanning of TIM revealed the tunnel region as a key binding site that can alter global dynamics of the enzyme.     

Molecular dynamics of mouse acetylcholinesterase complexed with huperzine A

Two molecular dynamics simulations were performed for a modeled complex of mouse acetylcholinesterase liganded with huperzine A (HupA). Analysis of these simulations shows that HupA shifts in the active site toward Tyr 337 and Phe 338, and that several residues in the active site area reach out to make hydrogen bonds with the inhibitor. Rapid fluctuations of the gorge width are observed, ranging from widths that allow substrate access to the active site, to pinched structures that do not allow access of molecules as small as water. Additional openings or channels to the active site are found. One opening is formed in the side wall of the active site gorge by residues Val 73, Asp 74, Thr 83, Glu 84, and Asn 87. Another opening is formed at the base of the gorge by residues Trp 86, Val 132, Glu 202, Gly 448, and Ile 451. Both of these openings have been observed separately in the Torpedo californica form of the enzyme. These channels could allow transport of waters and ions to and from the bulk solution. © 1999 John Wiley & Sons, Inc. Biopoly 50: 347–359, 1999     

Crystallographic studies on Ascaris suum NAD-malic enzyme bound to reduced cofactor and identification of an effector site

The crystal structure of the mitochondrial NAD-malic enzyme from Ascaris suum, in a quaternary complex with NADH, tartronate, and magnesium has been determined to 2.0-A resolution. The structure closely resembles the previously determined structure of the same enzyme in binary complex with NAD. However, a significant difference is observed within the coenzyme-binding pocket of the active site with the nicotinamide ring of NADH molecule rotating by 198 degrees over the C-1-N-1 bond into the active site without causing significant movement of the other catalytic residues. The implications of this conformational change in the nicotinamide ring to the catalytic mechanism are discussed. The structure also reveals a binding pocket for the divalent metal ion in the active site and a binding site for tartronate located in a highly positively charged environment within the subunit interface that is distinct from the active site. The tartronate binding site, presumably an allosteric site for the activator fumarate, shows striking similarities and differences with the activator site of the human NAD-malic enzyme that has been reported recently. Thus, the structure provides additional insights into the catalytic as well as the allosteric mechanisms of the enzyme.     

Rate-determining step of butyrylcholinesterase-catalyzed hydrolysis of benzoylcholine and benzoylthiocholine. Volumetric study of wild-type and D70G mutant behavior.

The rate-limiting step for hydrolysis of the positively charged oxoester benzoylcholine (BzCh) by human butyrylcholinesterase (BuChE) is deacylation (k(3)), whereas it is acylation (k(2)) for hydrolysis of the homologous thioester benzoylthiocholine (BzSCh). Steady-state hydrolysis of BzCh and BzSCh by wild-type BuChE and its peripheral anionic site mutant D70G was investigated at different hydrostatic pressures, which allowed determination of volume changes associated with substrate binding, and the activation volumes for the chemical steps. A differential nonlinear pressure-dependence of the catalytic parameters for hydrolysis of both substrates by both enzymes was shown. Nonlinearity of the plots may be explained in terms of compressibility changes or rate-limiting changes. To distinguish between these two possibilities, enzyme phosphorylation by diisopropylfluorophosphate (DFP) in the presence of substrate (BzSCh) under pressure was studied. There was no pressure dependence of volume changes for DFP binding or for phosphorylation of either wild-type or D70G. Analysis of the pressure dependence for steady-state hydrolysis of substrates, and for phosphorylation by DFP provided evidence that no enzyme compressibility changes occurred during the catalyzed reactions. Thus, the nonlinear pressure dependence of substrate hydrolysis reflects changes in the rate-limiting step with pressure. Change in rate-determining step occurred at a pressure of 100 MPa for hydrolysis of BzCh by wild-type and at 75 MPa for D70G. For hydrolysis of BzSCh the change occurred at higher pressures because k(2) < k(3) at atmospheric pressure for this substrate. Elementary volume change contributions upon initial binding, productive binding, acylation and deacylation were calculated from the pressure differentiation of kinetic constants. This analysis shed light on the molecular events taking place along the hydrolysis pathways of BzCh and BzSCh by wild-type BuChE and the D70G mutant. In addition, volume change differences between wild-type and D70G provided new evidence that residue D70 in the peripheral site controls hydration of the active site gorge and the dynamics of the water molecule network during catalysis. Finally, a steady-state kinetic study of the oxyanion hole mutant (G117H) showed that substitution of the ethereal sulfur for oxygen in the substrate alters the final adjustment of substrate in the active site and stabilization of the acylation transition state.     

Crystal structure analysis,covalent docking,and molecular dynamics calculations reveal a conformational switch in PhaZ7 PHB depolymerase

An open and a closed conformation of a surface loop in PhaZ7 extracellular poly(3‐hydroxybutyrate) depolymerase were identified in two high‐resolution crystal structures of a PhaZ7 Y105E mutant. Molecular dynamics (MD) simulations revealed high root mean square fluctuations (RMSF) of the 281–295 loop, in particular at residue Asp289 (RMSF 7.62 Å). Covalent docking between a 3‐hydroxybutyric acid trimer and the catalytic residue Ser136 showed that the binding energy of the substrate is significantly more favorable in the open loop conformation compared to that in the closed loop conformation. MD simulations with the substrate covalently bound depicted 1 Å RMSF higher values for the residues 281–295 in comparison to the apo (substrate‐free) form. In addition, the presence of the substrate in the active site enhanced the ability of the loop to adopt a closed form. Taken together, the analysis suggests that the flexible loop 281–295 of PhaZ7 depolymerase can act as a lid domain to control substrate access to the active site of the enzyme. Proteins 2017; 85:1351–1361. © 2017 Wiley Periodicals, Inc.     

Dynamics of the native and the ligand-bound structures of eosinophil cationic protein: network of hydrogen bonds at the catalytic site

Eosinophil Cationic Protein (ECP) is sequentially and structurally similar to ribonuclease A (RNase A). It belongs to the RNase A family of proteins and the RNA catalysis is essential to its biological function. In the present study, we have generated the dinucleotide-bound structures of ECP by docking the dinucleotides to a number of molecular dynamics (MD) generated ECP structures. The stability of the docked enzyme-ligand complexes was ascertained by extensive MD simulations. The modes of ligand binding are explored by essential dynamics studies. The role of water molecules in the stability of the complex and in the catalysis was investigated. The active site residues form a complex network of connections with the ligand and with a water molecule. The catalytic mechanism of the RNA cleavage is examined on the basis of the active site geometry obtained by the simulations.     

Mouse acetylcholinesterase unliganded and in complex with huperzine A: A comparison of molecular dynamics simulations

A 1 ns molecular dynamics simulation of unliganded mouse acetylcholinesterase (AChE) is compared to a previous simulation of mouse AChE complexed with huperzine A (HupA). Several common features are observed. In both simulations, the active site gorge fluctuates in size during the 1 ns trajectory and is completely pinched off several times. Many of the residues in the gorge that formed hydrogen bonds with HupA in the simulation of the complex now form hydrogen bonds with other protein residues and water molecules in the gorge. The opening of a “backdoor” entrance to the active site that was found in the simulation of the complex is also observed in the unliganded simulation. Differences between the two simulations include overall lower structural rms deviations for residues in the gorge in the unliganded simulation, a smaller diameter of the gorge in the absence of HupA, and the disappearance of a side channel that was frequently present in the liganded simulation. The differences between the two simulations can be attributed, in part, to the interaction of AChE with HupA. © 1999 John Wiley & Sons, Inc. Biopoly 50: 35–43, 1999     

Cofactor assisted gating mechanism in the active site of NADH oxidase from Thermus thermophilus

NADH oxidase (NOX) from Thermus thermophilus is a member of a structurally homologous flavoprotein family of nitroreductases and flavin reductases. The importance of local conformational dynamics in the active site of NOX has been recently demonstrated. The enzyme activity was increased by 250% in the presence of 1 M urea with no apparent perturbation of the native structure of the protein. The present in silico results correlate with the in vitro data and suggest the possible explanation about the effect of urea on NOX activity at the molecular level. Both, X-ray structure and molecular dynamics (MD) simulations, show open conformation of the active site represented by approximately 0.9 nm distance between the indole ring of Trp47 and the isoalloxazine ring of FMN412. In this conformation, the substrate molecule can bind in the active site without sterical restraints. MD simulations also indicate more stable conformation of the active site called "closed" conformation. In this conformation, Trp47 and the isoalloxazine ring of FMN412 are so close to each other (approximately 0.5 nm) that the substrate molecule is unable to bind between them without perturbing this conformation. The open/close transition of the active site between Trp47 and the flavin ring is accompanied by release of the "tightly" bound water molecule from the active site--cofactor assisted gating mechanism. The presence of urea in aqueous solutions of NOX prohibits closing of the active site and even unlocks the closed active site because of the concomitant binding of a urea molecule in the active site cavity. The binding of urea in the active site is stabilized by formation of one/two persistent hydrogen bonds involving the carbonyl group of the urea molecule. Our report represents the first MD study of an enzyme from the novel flavoprotein family of nitroreductases and flavin reductases. The common occurrence of aromatic residues covering the active sites in homologous enzymes suggests the possibility of a general gating mechanism and the importance of local dynamics within this flavoprotein family.     

Specific potassium ion interactions facilitate homocysteine binding to betaine‐homocysteine S‐methyltransferase

Betaine‐homocysteine S‐methyltransferase (BHMT) is a zinc‐dependent methyltransferase that uses betaine as the methyl donor for the remethylation of homocysteine to form methionine. This reaction supports S‐adenosylmethionine biosynthesis, which is required for hundreds of methylation reactions in humans. Herein we report that BHMT is activated by potassium ions with an apparent K_M for K⁺ of about 100 µM. The presence of potassium ions lowers the apparent K_M of the enzyme for homocysteine, but it does not affect the apparent K_M for betaine or the apparent k_cat for either substrate. We employed molecular dynamics (MD) simulations to theoretically predict and protein crystallography to experimentally localize the binding site(s) for potassium ion(s). Simulations predicted that K⁺ ion would interact with residues Asp26 and/or Glu159. Our crystal structure of BHMT bound to homocysteine confirms these sites of interaction and reveals further contacts between K⁺ ion and BHMT residues Gly27, Gln72, Gln247, and Gly298. The potassium binding residues in BHMT partially overlap with the previously identified DGG (Asp26‐Gly27‐Gly28) fingerprint in the Pfam 02574 group of methyltransferases. Subsequent biochemical characterization of several site‐specific BHMT mutants confirmed the results obtained by the MD simulations and crystallographic data. Together, the data herein indicate that the role of potassium ions in BHMT is structural and that potassium ion facilitates the specific binding of homocysteine to the active site of the enzyme. Proteins 2014; 82:2552–2564. © 2014 Wiley Periodicals, Inc.     

Insights into the structure and dynamics of the dinuclear zinc beta-lactamase site from Bacteroides fragilis

Herein, we report quantum chemical calculations and molecular dynamics (MD) simulations of the dinuclear form of the Bacteroides fragilis zinc beta-lactamase. We studied four different configurations which differ in the protonation state of the Asp103 residue and in the presence or absence of a Zn1-OH-Zn2 bridge. The flexibility of the Zn1-OH-Zn2 bridge was studied by means of quantum mechanical (QM) calculations on cluster models while the relative stabilities of the different configurations were estimated from QM linear scaling calculations on the enzyme. Contacts between important residues (Cys104, Asp69, Lys185, etc.), the solvation of the zinc ions, and the conformation of the active site beta-hairpin loop were characterized by the MD analyses. The influence of the buried sodium ion close to the Zn2 position was investigated by carrying out a secondary simulation where the sodium ion was replaced with an internal water molecule. The comparative structural analyses among the different MD trajectories augmented with energetic calculations have demonstrated that the B. fragilis protein efficiently binds the internal Na(+) ion observed crystallographically. Moreover, we found that when Asp103 is unprotonated, a rigid Zn1-OH-Zn2 bridge results, while for neutral Asp103, a fluctuating Zn1-Zn2 distance was possible via the breaking and formation of the Zn1-OH-Zn2 bridge. The mechanistic implications of these observations are discussed in detail.