The electrostatic driving force for nucleophilic catalysis in L-arginine deiminase: a combined experimental and theoretical study

L-arginine deiminase (ADI) catalyzes the hydrolysis of L-arginine to form L-citrulline and ammonia via two partial reactions. A working model of the ADI catalytic mechanism assumes nucleophilic catalysis by a stringently conserved active site Cys and general acid-general base catalysis by a stringently conserved active site His. Accordingly, in the first partial reaction, the Cys attacks the substrate guanidino C zeta atom to form a tetrahedral covalent adduct, which is protonated by the His at the departing ammonia group to facilitate the formation of the Cys- S-alkylthiouronium intermediate. In the second partial reaction, the His activates a water molecule for nucleophilic addition at the thiouronium C zeta atom to form the second tetrahedral intermediate, which eliminates the Cys in formation of the L-citrulline product. The absence of a basic residue near the Cys thiol suggested that the electrostatic environment of the Cys thiol, in the enzyme-substrate complex, stabilizes the Cys thiolate anion. The studies described in this paper explore the mechanism of stabilization of the Cys thiolate. First, the log(k(cat)/K(m)) and log k(cat) pH rate profiles were measured for several structurally divergent ADIs to establish the pH range for ADI catalysis. All ADIs were optimally active at pH 5, which suggested that the Cys pKa is strongly perturbed by the prevailing electrostatics of the ADI active site. The p K a of the Bacillus cereus ADI (BcADI) was determined by UV-pH titration to be 9.6. In contrast, the pKa determined by iodoacetamide Cys alkylation is 6.9. These results suggest that the negative electrostatic field from the two opposing Asp carboxylates perturbs the Cys pKa upward in the apoenzyme and that the binding of the iodoacetamide (a truncated analogue of the citrulline product) between the Cys thiol and the two Asp carboxylates shields the Cys thiol, thereby reducing its pKa. It is hypothesized that the bound positively charged guanidinium group of the L-arginine substrate further stabilizes the Cys thiolate. The so-called "substrate-assisted" Cys ionization, first reported by Fast and co-workers to operate in the related enzyme dimethylarginine dimethylaminohydrolase [Stone, E. M., Costello, A. L., Tierney, D. L., and Fast, W. (2006) Biochemistry 45, 5618-5630], was further explored computationally in ADI by using an ab initio quantum mechanics/molecular mechanics method. The energy profiles for formation of the tetrahedral intermediate in the first partial reaction were calculated for three different reaction scenarios. From these results, we conclude that catalytic turnover commences from the active configuration of the ADI(L-arginine) complex which consists of the Cys thiolate (nucleophile) and His imidazolium ion (general acid) and that the energy barriers for the nucleophilic addition of Cys thiolate to the thiouronium C zeta atom and His imidazolium ion-assisted elimination from the tetrahedral intermediate are small.     

Role of active site binding interactions in 4-chlorobenzoyl-coenzyme A dehalogenase catalysis.

4-Chlorobenzoyl-coenzyme A (4-CBA-CoA) dehalogenase catalyzes the hydrolytic dehalogenation of 4-CBA-CoA to 4-hydroxybenzoyl-CoA (4-HBA-CoA) via a multistep mechanism involving initial attack of Asp145 on C(4) of the substrate benzoyl ring to form a Meisenheimer intermediate (EMc), followed by expulsion of the chloride ion to form an arylated enzyme intermediate (EAr) and then ester hydrolysis in the EAr to form product. This study examines the role of binding interactions in dehalogenase catalysis. The enzyme and substrate groups positioned for favorable binding interaction were identified from the X-ray crystal structure of the enzyme-4-HBA-3'-dephospho-CoA complex. These groups were individually modified (via site-directed mutagenesis or chemical synthesis) for the purpose of disrupting the binding interaction. The changes in the Gibbs free energy of the enzyme-substrate complex (DeltaDeltaG(ES)) and enzyme-transition state complex (DeltaDeltaG) brought about by the modification were measured. Cases where DeltaDeltaG exceeds DeltaDeltaG(ES) are indicative of binding interactions used for catalysis. On the basis of this analysis, we show that the H-bond interactions between the Gly114 and Phe64 backbone amide NHs and the substrate benzoyl C=O group contribute an additional 3.1 kcal/mol of stabilization at the rate-limiting transition state. The binding interactions between the enzyme and the substrate CoA nucleotide moiety also intensify in the rate-limiting transition state, reducing the energy barrier to catalysis by an additional 3.3 kcal/mol. Together, these binding interactions contribute approximately 10(6) to the k(cat)/K(m).     

Helix-coil transition theory including long-range electrostatic interactions: application to globular proteins

An extension of the Zimm–Bragg two-state theory for the helix–coil transition in polypeptides, which takes into account the effect of peptide charge–dipole interactions on helix stability, is presented. This new theory incorporates these interactions in an expression that is parameterized on recently obtained experimental data on polypeptides for which electrostatic effects are known to influence helix content. Unlike previous two-state or multistate models, which are parameterized on protein x-ray data, the present theoretical treatment in independent of such protein data. The theoretical model is applied to a series of peptides derived from the C-peptide of ribonuclease A, which have been the object of recent spectroscopic studies. The new theoretical approach can account for most of the structural information derived from studies of these C-peptides, and for overall average helix probabilities that are close in magnitude to those observed for these polypeptides in solution. An application of this new formulation for the prediction of the locations of α-helices in globular proteins from their amino acid sequence is also presented.     

Spermine-nucleic acid interactions: a theoretical study

The interaction of spermine with nucleic acids is simulated theoretically using refined semi-empirical energy formulae and an advanced minimization procedure. Various nucleic acids are considered: model homopolymeric DNA's, a dodecamer (CGCGAATTCGCG) of type B-DNA, as well as a transfer RNA, tRNA^Phe. The dominant role of electrostatic potential in determining the preferential binding sites of spermine is demonstrated in each of these cases and the role of counterions, nucleic acid structure, and base-pair sequence is analyzed.     

Comparison of different schemes to treat long-range electrostatic interactions in molecular dynamics simulations of a protein crystal

Eight molecular dynamics simulations of a ubiquitin crystal unit cell were performed to investigate the effect of different schemes to treat the long-range electrostatic interactions as well as the need to include counter ions. A crystal system was chosen as the test system, because the higher charge density compared with a protein in solution makes it more sensitive to the way of treating the electrostatic interactions. Three different schemes of treating the long-range interactions were compared: straight cutoff, reaction-field approximation, and a lattice-sum method (P3M). For each of these schemes, two simulations were performed, one with and one without the counter ions. Two additional simulations with a reaction-field force and different initial placements of the counter ions were performed to examine the effect of the initial positions of the ions. The inclusion of long-range electrostatic interactions using either a reaction-field or a lattice-sum method proved to be necessary for the simulation of crystals. These two schemes did not differ much in their ability to reproduce the crystallographic structure. The inclusion of counter ions, on the other hand, seems not necessary for obtaining a stable simulation. The initial positions of the ions have a visible but small effect on the simulation.     

sp C-H and sp C-H agostic ruthenium complexes: a combined experimental and theoretical study

Halide abstraction from the 18 electron Ru(II) complex RuCl(CO)₂[2,6-(CH₂P^tBu₂)₂C₆H₃] (2) with AgPF₆ results in the exclusive formation of the cationic complex {Ru(CO)₂[2,6-(CH₂P^tBu₂)₂C₆H₃]}⁺PF₆⁻ (3). The molecular structures of 2 and 3 were determined by complete single-crystal diffraction studies. X-ray crystallographic analysis of 3 reveals that the “open” coordination site is occupied by an agostic interaction between the metal center and an sp³ C-H bond of a tert-butyl substituent. DFT gas phase calculations (B97-1/SDD) show the necessity of two sterically demanding tert-butyl substituents on one P donor atom for the agostic interaction to occur. The reaction of 3 with H₂ results in the quantitative conversion to {Ru(H)(CO)₂[2,6-(CH₂P^tBu₂)₂C₆H₄]}⁺PF₆⁻ (4) where the aromatic C_ipso-H bond is η²-coordinated to the metal center. Treatment of the agostic complex 4 with Et₃N results in the formation of the neutral complex Ru(H)(CO)₂[2,6-(CH₂P^tBu₂)₂C₆H₃] (5). The mechanistic details of 3 + H₂ → 4 were investigated by DFT calculations at the B97-1/SDB-cc-pVDZ//B97-1/SDD level of theory.     

Cumulative site-directed charge-change replacements in bacteriophage T4 lysozyme suggest that long-range electrostatic interactions contribute little to protein stability

Bacteriophage T4 lysozyme is a basic molecule with an isoelectric point above 9.0, and an excess of nine positive charges at neutral pH. It might be expected that it would be energetically costly to bring these out-of-balance charges from the extended, unfolded, form of the protein into the compact folded state. To determine the contribution of such long-range electrostatic interactions to the stability of the protein, five positively charged surface residues, Lys16, Arg119, Lys135, Lys147 and Arg154, were individually replaced with glutamic acid. Eight selected double, triple and quadruple mutants were also constructed so as to sequentially reduce the out-of-balance formal charge on the molecule from +9 to +1 units. Each of the five single variant proteins was crystallized and high-resolution X-ray analysis confirmed that each mutant structure was, in general, very similar to the wild-type. In the case of R154E, however, the Arg154 to Glu replacement caused a rearrangement in which Asp127 replaced Glu128 as the capping residue of a nearby alpha-helix. The thermal stabilities of all 13 variant proteins were found to be fairly similar, ranging from 0.5 kcal/mol more stable than wild-type to 1.7 kcal/mol less stable than wild-type. In the case of the five single charge-change variants, for which the structures were determined, the changes in stability can be rationalized in terms of changes in local interactions at the site of the replacement. There is no evidence that the reduction in the out-of-balance charge on the molecule increases the stability of the folded relative to the unfolded form, either at pH 2.8 or at pH 5.3. This indicates that long-range electrostatic interactions between the substituted amino acid residues and other charged groups on the surface of the molecule are weak or non-existent. Furthermore, the relative stabilities of the multiple charge replacement mutant proteins were found to be almost exactly equal to the sums of the relative stabilities of the constituent single mutant proteins. This also clearly indicates that the electrostatic interactions between the replaced charges are negligibly small. The activities of the charge-change mutant lysozymes, as measured by the rate of hydrolysis of cell wall suspensions, are essentially equal to that of the wild-type lysozyme, but on a lysoplate assay the mutant enzymes appear to have higher activity.(ABSTRACT TRUNCATED AT 400 WORDS)     

Contributions to catalysis and potential interactions of the three catalytic domains in a contiguous trimeric creatine kinase

Three separate creatine kinase (CK) isoform families exist in animals. Two of these (cytoplasmic and mitochondrial) are obligate oligomers. A third, flagellar, is monomeric but contains the residues for three complete CK domains. It is not known whether the active sites in each of the contiguous flagellar domains are catalytically competent, and, if so, whether they are capable of acting independently. Here we have utilized site-directed mutagenesis to selectively disable individual active sites and all possible combinations thereof. Kinetic studies showed that these mutations had minimal impact on substrate binding and synergism. Interestingly, the active sites were not catalytically equivalent, and were in fact interdependent, a phenomenon that has previously been reported only in the oligomeric CK isoforms.     

Novel derivatives of ansa-titanocenes procured from 6-phenylfulvene: a combined experimental and theoretical study

The previously prepared trans-[(1,2-diphenyl-1,2-dicyclopentadienyl)ethanediyl] titanium(IV) dichloride, [1,2-(Ph)₂C₂H₂{η⁵-C₅H₄}₂]Ti(Cl)₂, was synthesised using an alternative procedure, from which its crystal structure was determined. Using this compound, a variety of other ansa-titanocene derivatives were synthesised by replacement of the chloride ligands with selected substituents. Thus RTi(X)(Y) systems where R=1,2-(Ph)₂C₂H₂η⁵-C₅H₄₂; X=Y=CH₃; X=CH₃, Y=Cl; X=Y=NCS; X=Y=NCO; X=Y=OPh and (X/Y)=O have been synthesised and characterised. DFT calculations were performed on the complexes trans-[(1,2-diphenyl-1,2-dicyclopentadienyl)-ethanediyl] titanium(IV) dichloride, bis-(6,6-diphenylfulvene)titanium and bis-(6,6-diphenylfulvene)iron. This demonstrated the role that the metal centre plays in their formation, generating either an ansa-metallocene, a ‘tucked in’ fulvene complex or a metallocene coordinating fulvene anions.     

Wobble base-pairing in codon-anticodon interactions: a theoretical modelling study

The Crick wobble hypothesis attributes the phenomenon of codon degeneracy to a certain impreciseness of pairing between the third base of the codon and the first base of the anticodon. This theoretical study investigates the pairing properties of some wobble bases, including both, observed and unobserved pairs. Some wobble base-pairs are predicted to follow the Watson-Crick pairs in configuration and pairing facility, while others deviate from this norm. The observed U:V pair is unique in that a pairing configuration may be suggested for it wherein the hydrogen-bonding involves the exocyclic 5-carboxymethoxy group of V. By comparing the theoretical data on the configurations of these pairs with the evidence for their existence/non-existence in nature, some guidelines emerge for differentiating between observed and unobserved base pairs on the basis of the pairing configuration.     

Conformation of d(CpG) modified by the carcinogen 4-aminobiphenyl: a combined experimental and theoretical analysis

The change in DNA conformation produced by the attachment of a reactive substance is likely to be a vital factor in determining the biological consequences of the reaction. We have prepared a deoxydinucleoside monophosphate containing the major adduct derived from the carcinogenic amine 4-aminobiphenyl and analyzed its conformation by theoretical and experimental methods. Reaction of d(CpG) with N-acetoxy-N-(trifluoroacetyl)-4-aminobiphenyl afforded the product modified at C-8 of guanine with 4-aminobiphenyl. After purification by reverse-phase high-performance liquid chromatography, milligram amounts of product were obtained. It was analyzed by circular dichroism, proton magnetic resonance, and minimized potential-energy calculations. A flexible molecule with a mixture of conformers is indicated. Both carcinogen-base-stacked states and base-base-stacked states, with guanine both syn anti, contribute to the population mixture on the dimer level. The global minimum-energy conformation has syn-guanine and carcinogen-base stacking. Forms of this type are calculated to represent roughly 58% of the conformer population. Because of the twisted nature of the biphenyl moiety, carcinogen-base stacking inherently involves less overlap than that in the planar and rigid three-ringed aminofluorene analogue. This difference might relate to the diminished effectiveness of the aminobiphenyl vs. the aminofluorene adduct as a frameshift mutagen in Salmonella typhimurium 1538.     

Nutrient utilization by bovine articular chondrocytes: a combined experimental and theoretical approach

A combined experimental-numerical approach was adopted to characterize glucose and oxygen uptake and lactate production by bovine articular chondrocytes in a model system. For a wide range of cell concentrations, cells in agarose were supplemented with either low or high glucose medium. During an initial culture phase of 48 h, oxygen was monitored noninvasively using a biosensor system. Glucose and lactate were determined by medium sampling. In order to quantify glucose and oxygen uptake, a finite element approach was adopted to describe diffusion and uptake in the experimental model. Numerical predictions of lactate, based on simple relations for cell metabolism, were found to agree well for low glucose, but not for high glucose medium. Oxygen did not play a role in either case. Given the close association between chondrocyte energy metabolism and matrix synthesis, a quantifiable prediction of utilization can present a valuable contribution in the optimization of tissue engineering conditions.     

A combined theoretical and experimental study of the interaction of metrizamide with proteins

A theory of isopycnic sedimentation of interacting systems provides a theoretical explanation for the isopycnic behavior of catalase and other proteins in metrizamide gradients of pH 7.0–7.5, particularly the bimodal banding patterns reported by some investigators. As for whether or not the model is realistic, a spectroscopic investigation provides an unambiguous demonstration of a reversible interaction between metrizamide and the catalase dimers formed by dissociation of catalase tetramer at pH 3. The metrizamide-catalase dimer complex(es) is of high bouyant density. In addition, metrizamide at pH 7.4 interacts with heat-denatured ovalbumin. These observations suggest that the denser band in the experimental bimodal isopycnic patterns might be a complex(es) of metrizamide with a comformationally altered state of the protein.     

Adsorption of DNA to mica mediated by divalent counterions: a theoretical and experimental study

The adsorption of DNA molecules onto a flat mica surface is a necessary step to perform atomic force microscopy studies of DNA conformation and observe DNA-protein interactions in physiological environment. However, the phenomenon that pulls DNA molecules onto the surface is still not understood. This is a crucial issue because the DNA/surface interactions could affect the DNA biological functions. In this paper we develop a model that can explain the mechanism of the DNA adsorption onto mica. This model suggests that DNA attraction is due to the sharing of the DNA and mica counterions. The correlations between divalent counterions on both the negatively charged DNA and the mica surface can generate a net attraction force whereas the correlations between monovalent counterions are ineffective in the DNA attraction. DNA binding is then dependent on the fractional surface densities of the divalent and monovalent cations, which can compete for the mica surface and DNA neutralizations. In addition, the attraction can be enhanced when the mica has been pretreated by transition metal cations (Ni(2+), Zn(2+)). Mica pretreatment simultaneously enhances the DNA attraction and reduces the repulsive contribution due to the electrical double-layer force. We also perform end-to-end distance measurement of DNA chains to study the binding strength. The DNA binding strength appears to be constant for a fixed fractional surface density of the divalent cations at low ionic strength (I < 0.1 M) as predicted by the model. However, at higher ionic strength, the binding is weakened by the screening effect of the ions. Then, some equations were derived to describe the binding of a polyelectrolyte onto a charged surface. The electrostatic attraction due to the sharing of counterions is particularly effective if the polyelectrolyte and the surface have nearly the same surface charge density. This characteristic of the attraction force can explain the success of mica for performing single DNA molecule observation by AFM. In addition, we explain how a reversible binding of the DNA molecules can be obtained with a pretreated mica surface.     

Towards numerical forecasting of long-range air transport of birch pollen: theoretical considerations and a feasibility study

This paper considers the feasibility of numerical simulation of large-scale atmospheric transport of allergenic pollen. It is shown that at least small grains, such as birch pollen, can stay in the air for a few days, which leads to a characteristic scale for their transport of ∼10³ km. The analytical consideration confirmed the applicability of existing dispersion models to the pollen transport task and provided some reference parameterizations of the key processes, including dry and wet deposition. The results were applied to the Finnish Emergency Dispersion Modelling System (SILAM), which was then used to analyze pollen transport to Finland during spring time in 2002–2004. Solutions of the inverse problems (source apportionment) showed that the main source areas, from which the birch flowering can affect Finnish territory, are the Baltic States, Russia, Germany, Poland, and Sweden—depending on the particular meteorological situation. Actual forecasting of pollen dispersion required a birch forest map of Europe and a unified European model for birch flowering, both of which were nonexistent before this study. A map was compiled from the national forest inventories of Western Europe and satellite images of broadleaf forests. The flowering model was based on the mean climatological dates for the onset of birch forests rather than conditions of any specific year. Utilization of probability forecasting somewhat alleviated the problem, but the development of a European-wide flowering model remains the main obstacle for real-time forecasting of large-scale pollen distribution.     

Nonlinear mechanical property of tracheal cartilage: a theoretical and experimental study

BACKGROUND: Despite being the stiffest airway of the bronchial tree, the trachea undergoes significant deformation due to intrathoracic pressure during breathing. The mechanical properties of the trachea affect the flow in the airway and may contribute to the biological function of the lung. METHOD: A Fung-type strain energy density function was used to investigate the nonlinear mechanical behavior of tracheal cartilage. A bending test on pig tracheal cartilage was performed and a mathematical model for analyzing the deformation of tracheal cartilage was developed. The constants included in the strain energy density function were determined by fitting the experimental data. RESULT: The experimental data show that tracheal cartilage is a nonlinear material displaying higher strength in compression than in tension. When the compression forces varied from -0.02 to -0.03N and from -0.03 to -0.04N, the deformation ratios were 11.03+/-2.18% and 7.27+/-1.59%, respectively. Both were much smaller than the deformation ratios (20.01+/-4.49%) under tension forces of 0.02 to 0.01N. The Fung-type strain energy density function can capture this nonlinear behavior very well, whilst the linear stress-strain relation cannot. It underestimates the stability of trachea by exaggerating the displacement in compression. This study may improve our understanding of the nonlinear behavior of tracheal cartilage and it may be useful for the future study on tracheal collapse behavior under physiological and pathological conditions.     

The magnitude of electrostatic interactions in inhibitor binding and during catalysis by ribonuclease A.

It is demonstrated that a model of nucleotide binding to ribonuclease A similar to that proposed by Hammes and coworkers (G. G. Hammes (1968), Adv. Protein Chem. 23, 1) is at least, approximately applicable for both cyclic nucleotide substrates and mononucleotide inhibitors at pH values less than or equal to 6.5 and as a function of ionic strength. Calorimetric data on various inhibitors show that the binding reaction can be thermodynamically dissected into a contribution arising from van der Waal's interaction of the nucleoside moiety, characterized by a large negative enthalpy change, and a contribution arising from electrostatic interactions between the negatively charged phosphate group of the inhibitor and the positively charged protein fabric, characterized by a large positive unitary entropy change. Assuming a catalytic mechanism involving the formation of a dianionic pentacoordinated phosphate transition state intermediate, the magnitude of the effect of electrostatic interactions on the overall rate enhancement by the enzyme is estimated to be 2 times 10(2) to 10(6). It is suggested that this effect, along with substrate approximation effects, is sufficient to "explain" the catalytic behavior of the enzyme.     

Mechano-electric interactions in heterogeneous myocardium: development of fundamental experimental and theoretical models

The heart is structurally and functionally a highly non-homogenous organ, yet its main function as a pump can only be achieved by the co-ordinated contraction of millions of ventricular cells. This apparent contradiction gives rise to the hypothesis that ‘well-organised’ inhomogeneity may be a pre-requisite for normal cardiac function. Here, we present a set of novel experimental and theoretical tools for the study of this concept. Heterogeneity, in its most condensed form, can be simulated using two individually controlled, mechanically interacting elements (duplex). We have developed and characterised three different types of duplexes: (i) biological duplex, consisting of two individually perfused biological samples (like thin papillary muscles or a trabeculae), (ii) virtual duplex, made-up of two interacting mathematical models of cardiac muscle, and (iii) hybrid duplex, containing a biological sample that interacts in real-time with a virtual muscle. In all three duplex types, in-series or in-parallel mechanical interaction of elements can be studied during externally isotonic, externally isometric, and auxotonic modes of contraction and relaxation. Duplex models, therefore, mimic (patho-)physiological mechano-electric interactions in heterogeneous myocardium at the multicellular level, and in an environment that allows one to control mechanical, electrical and pharmacological parameters. Results obtained using the duplex method show that: (i) contractile elements in heterogeneous myocardium are not ‘independent’ generators of tension/shortening, as their ino- and lusitropic characteristics change dynamically during mechanical interaction—potentially matching microscopic contractility to macroscopic demand, (ii) mechanical heterogeneity contributes differently to action potential duration (APD) changes, depending on whether mechanical coupling of elements is in-parallel or in-series, which may play a role in mechanical tuning of distant tissue regions, (iii) electro-mechanical activity of mechanically interacting contractile elements is affected by their activation sequence, which may optimise myocardial performance by smoothing intrinsic differences in APD. In conclusion, we present a novel set of tools for the experimental and theoretical investigation of cardiac mechano-electric interactions in healthy and/or diseased heterogeneous myocardium, which allows for the testing of previously inaccessible concepts.     

Cyanolysis and azidolysis of epoxides by haloalcohol dehalogenase: theoretical study of the reaction mechanism and origins of regioselectivity

Haloalcohol dehalogenase HheC catalyzes the reversible dehalogenation of vicinal haloalcohols to form epoxides and free halides. In addition, HheC is able to catalyze the irreversible and highly regioselective ring-opening of epoxides with nonhalide nucleophiles, such as CN (-) and N 3 (-). For azidolysis of aromatic epoxides, the regioselectivity observed with HheC is opposite to the regioselectivity of the nonenzymatic epoxide-opening. This, together with a relatively broad substrate specificity, makes HheC a promising tool for biocatalytic applications. We have designed large quantum chemical models of the HheC active site and used density functional theory to study the reaction mechanism of the HheC-catalyzed ring-opening of ( R)-styrene oxide with the nucleophiles CN (-) and N 3 (-). Both the cyanolysis and the azidolysis reactions are shown to take place in a single concerted step. The results support the suggested role of the putative Ser132-Tyr145-Arg149 catalytic triad, where Tyr145 acts as a general acid, donating a proton to the substrate, and Arg149 interacts with Tyr145 and facilitates proton abstraction, while Ser132 positions the substrate and reduces the barrier for epoxide opening through interaction with the emerging oxyanion of the substrate. We have also studied the regioselectivity of ( R)-styrene oxide opening for both the cyanolysis and the azidolysis reactions. The employed active site model was shown to be able to reproduce the experimentally observed beta-regioselectivity of HheC. In silico mutations of various groups in the HheC active site model were performed to elucidate the important factors governing the regioselectivity.