Comparative study on the active sites of ficin and papain by the spin labeling method

Ficin was alkylated with a series of haloacetamide spin labels with various distances between the spin probes and reactive groups. From the relation of these distances to the tau c values of the labels incorporated into protein, it was estimated that the depth of the active site hole of ficin is ca. 8 A. The results are somewhat different from those reported previously for papain (S. Nakayama et al. (1981) Biochem. Biophys. Res. Commun. 98, 471-475). Examination of the pH dependence of the ESR spectra for ficin and papain alkylated with an iodoacetamide or a maleimide spin label suggested that these enzymes have an amino acid residue of pKa 4 (probably a histidine residue) around the active site cysteine and that the active site conformations change at around pH 5.     

Conformation of the active site of thiolsubtilisin: reaction with specific chloromethyl ketones and arylacryloylimidazoles.

The conformation of the active site of thiolsubtilisin, prepared from subtilisin by transformation of the active site Ser to Cys, was compared with that of subtilisin by kinetic and spectroscopic methods. Carbobenzyloxy-L-alanylglycyl-L-phenylalanine chloromethyl ketone inhibited thiolsubtilisin approximately 10(2) times faster than subtilisin; alkylation occurred at the sulfhydryl rather than the imidazolyl group of the active site. pH dependence of the inhibition is different from that of the reaction between a simple thiol with haloacetamide. Furthermore, several native chromophoric arylacryloyl-thiolsubtilisins and arylacryloyl-subtilisins showed similar red shifts when compared with their denatured forms. The rate of deacylation of arylacryloyl-thiolsubtilisins was faster than (or of the same order of magnitude as) the deacylation rate of the analogous arylacryloyl-subtilisins in 30% dioxane (v/v), pH 5--10. The deacylation rate--pH profiles of these arylacryloyl-thiolsubtilisins in 30% dioxane all give pK values of 7.7 which is identical with the pK in the deacylation of acyl-subtilisins. These facts strongly suggest that the active-site conformation remains intact on conversion from subtilisin to thiolsubtilisin. The low esterase and peptidase activities of thiolsubtilisin are most likely due to the relatively low basicity of -SH (compared with -OH).     

Mechanism of the binding of Z-L-tryptophan and Z-L-phenylalanine to thermolysin and stromelysin-1 in aqueous solutions

The chemical shift of the carboxylate carbon of Z-tryptophan is increased from 179.85 to 182.82 ppm and 182.87 ppm on binding to thermolysin and stromelysin-1 respectively. The chemical shift of Z-phenylalanine is also increased from 179.5 ppm to 182.9 ppm on binding to thermolysin. From pH studies we conclude that the pK_a of the inhibitor carboxylate group is lowered by at least 1.5 pK_a units when it binds to either enzyme. The signal at ~ 183 ppm is no longer observed when the active site zinc atom of thermolysin or stromelysin-1 is replaced by cobalt. We estimate that the distance of the carboxylate carbon of Z-[1-¹³C]-L-tryptophan is ≤ 3.71 Å from the active site cobalt atom of thermolysin. We conclude that the side chain of Z-[1-¹³C]-L-tryptophan is not bound in the S₂′ subsite of thermolysin. As the chemical shifts of the carboxylate carbons of the bound inhibitors are all ~ 183 ppm we conclude that they are all bound in a similar way most probably with the inhibitor carboxylate group directly coordinated to the active site zinc atom. Our spectrophotometric results confirm that the active site zinc atom is tetrahedrally coordinated when the inhibitors Z-tryptophan or Z-phenylalanine are bound to thermolysin.     

How the oxygen tolerance of a [NiFe]-hydrogenase depends on quaternary structure

‘Oxygen-tolerant’ [NiFe]-hydrogenases can catalyze H₂ oxidation under aerobic conditions, avoiding oxygenation and destruction of the active site. In one mechanism accounting for this special property, membrane-bound [NiFe]-hydrogenases accommodate a pool of electrons that allows an O₂ molecule attacking the active site to be converted rapidly to harmless water. An important advantage may stem from having a dimeric or higher-order quaternary structure in which the electron-transfer relay chain of one partner is electronically coupled to that in the other. Hydrogenase-1 from E. coli has a dimeric structure in which the distal [4Fe-4S] clusters in each monomer are located approximately 12 Å apart, a distance conducive to fast electron tunneling. Such an arrangement can ensure that electrons from H₂ oxidation released at the active site of one partner are immediately transferred to its counterpart when an O₂ molecule attacks. This paper addresses the role of long-range, inter-domain electron transfer in the mechanism of O₂-tolerance by comparing the properties of monomeric and dimeric forms of Hydrogenase-1. The results reveal a further interesting advantage that quaternary structure affords to proteins.     

On the Permeation by Dioxygen of Urate Oxidase from Aspergillus flavus in Complex with Xanthine Anion: Dioxygen Pathways and a Portrait of the Enzyme Cavities from Molecular Dynamics Simulations in Water Solution

This work describes molecular dynamics (MD) simulations in aqueous media for the complex of the homotetrameric urate oxidase (UOX) from Aspergillus flavus with xanthine anion ( 5 ) in the presence of dioxygen (O₂). After 196.6 ns of trajectory from unrestrained MD, a O₂ molecule was observed leaving the bulk solvent to penetrate the enzyme between two subunits, A/C. From here, the same O₂ molecule was observed migrating, across subunit C, to the hydrophobic cavity that shares residue V227 with the active site. The latter was finally attained, after 378.3 ns of trajectory, with O₂ at a bonding distance from 5 . The reverse same O₂ pathway, from 5 to the bulk solvent, was observed as preferred pathway under random acceleration MD (RAMD), where an external, randomly oriented force was acting on O₂. Both MD and RAMD simulations revealed several cavities populated by O₂ during its migration from the bulk solvent to the active site or backwards. Paying attention to the last hydrophobic cavity that apparently serves as O₂ reservoir for the active site, it was noticed that its volume undergoes ample fluctuations during the MD simulation, as expected from the thermal motion of a flexible protein, independently from the particular subunit and no matter whether the cavity is filled or not by O₂.     

Studies of the mechanism of action and regulation of cAMP-dependent protein kinase

Magnetic resonance and kinetic studies of the catalytic subunit of a Type II cAMP-dependent protein kinase from bovine heart have established the active complex to be an enzyme-ATP-metal bridge. The metal ion is β,γ coordinated with Δ chirality at the β-phosphorous atom. The binding of a second metal ion at the active site which bridges the enzyme to the three phosphoryl groups of ATP, partially inhibits the reaction. Binding of the metal-ATP substrate to the enzyme occurs in a diffusion-controlled reaction followed by a 40 ° change in the glycosidic torsional angle. This conformational change results from strong interaction of the nucleotide base with the enzyme. NMR studies of four ATP-utilizing enzymes show a correlation between such conformational changes and high nucleotide base specificity. Heptapeptide substrates and substrate analogs bind to the active site of the catalytic subunit at a rate significantly lower than collision frequency indicating conformational selection by the enzyme or a subsequent slow conformational change. NMR studies of the conformation of the enzyme-bound peptide substrates have ruled out α-helical and β-pleated sheet structures. The results of kinetic studies of peptide substrates in which the amino acid sequence was systematically varied were used to rule out the obligatory requirement for all possible β-turn conformations within the heptapeptide although an enzymatic preference for a β_2–5 or β_3–6 turn could not be excluded. Hence if protein kinase has an absolute requirement for a specific secondary structure, then this structure must be a coil. In the enzyme-substrate complex the distance along the reaction coordinate between the γ-P of ATP and the serine oxygen of the peptide substrate (5.3 ± 0.7 Å) allows room for a metaphosphate intermediate. This finding together with kinetic observations as well as the location of the inhibitory metal suggest a dissociative mechanism for protein kinase, although a mechanism with some associative character remains possible. Regulation of protein kinase is accomplished by competition between the regulatory subunit and peptide or protein substrates at the active site of the catalytic subunit. Thus, the regulatory subunit is found by NMR to block the binding of the peptide substrate to the active site of protein kinase but allows the binding of the nucleotide substrate and divalent cations. The dissociation constant of the regulatory subunit from the active site (10^?10m) is increased ~10-fold by phosphorylation and ~10⁴-fold by the binding of cAMP, to a value (10^?5m) which exceeds the intracellular concentration of the R₂C₂ holoenzyme complex (10^?6m). The resulting dissociation of the holoenzyme releases the catalytic subunit, permitting the active site binding of peptide or protein substrates.     

Mapping the active site of factor Xa by selective inhibitors: An NMR and MD study

The structure of two selective inhibitors, Ac-Tyr-Ile-Arg-Ile-Pro-NH₂ and Ac-(4-Amino-Phe)-(Cyclohexyl-Gly)-Arg-NH₂, in the active site of the blood clotting enzyme factor Xa was determined by using transferred nuclear Overhauser effect nuclear magnetic resonance (NMR) spectroscopy. They represent a family of peptidic inhibitors obtained by the screening of a vast combinatorial library. Each structure was first calculated by using standard computational procedures (distance geometry, simulated annealing, energy minimization) and then further refined by systematic search of the conformation of the inhibitor docked in the active site and repeating the simulated annealing and energy minimization. The final structure was optimized by molecular dynamics simulations of the inhibitor-complex in water. The NMR restraints were kept throughout the refinement. The inhibitors assume a compact, very well defined conformation, embedded into the substrate binding site not in the same way as a substrate, blocking thus the catalysis. The model allows to explain the mode of action, affinity, and specificity of the peptides and to map the active site. Proteins 30:264–279, 1998. © 1998 Wiley-Liss, Inc.     

Influence of the protein structure surrounding the active site on the catalytic activity of [NiFeSe] hydrogenases

A combined experimental and theoretical study of the catalytic activity of a [NiFeSe] hydrogenase has been performed by H/D exchange mass spectrometry and molecular dynamics simulations. Hydrogenases are enzymes that catalyze the heterolytic cleavage or production of H₂. The [NiFeSe] hydrogenases belong to a subgroup of the [NiFe] enzymes in which a selenocysteine is a ligand of the nickel atom in the active site instead of cysteine. The aim of this research is to determine how much the specific catalytic properties of this hydrogenase are influenced by the replacement of a sulfur by selenium in the coordination of the bimetallic active site versus the changes in the protein structure surrounding the active site. The pH dependence of the D₂/H⁺ exchange activity and the high isotope effect observed in the Michaelis constant for the dihydrogen substrate and in the single exchange/double exchange ratio suggest that a “cage effect” due to the protein structure surrounding the active site is modulating the enzymatic catalysis. This “cage effect” is supported by molecular dynamics simulations of the diffusion of H₂ and D₂ from the outside to the inside of the protein, which show different accumulation of these substrates in a cavity next to the active site.     

Shifting the optimum pH of Bacillus circulans xylanase towards acidic side by introducing arginine

Electrostatic interactions are important in protein folding, binding, flexibility, stability and function. The pH at which the enzyme is maximally active is determined by the pK_as of the active site residues, which are modulated by several factors including the change in electrostatics in its vicinity. As the acidic xylanases are important in food and animal feed industries, electrostatic interactions are introduced in Bacillus circulans xylanase to shift their pH optima towards the acidic side. Arg substitutions are made to modulate the pKas of the active site residues. Neutral residues are substituted by Arg in such a way that the substituted residue can make direct interaction with the catalytic residues. However, the mutations with other titratable residues (Asp, Arg, Lys, His, Tyr, and Ser) present in between the catalytic sites and the substituted sites are avoided. Site directed mutagenesis was conducted to confirm the strategy. The results show the shift in pH optima of the mutants towards the acidic side by 0.5–1.5 unit. Molecular dynamics simulation of the mutant V37R reveals that the decrease in activity is due to the increase in distance between the substrate oxygen atoms and catalytic glutamates.     

The distance between the anticodon loops of two tRNAs bound to the 70 S Escherichia coli ribosome.

Using singlet-singlet energy transfer, we have measured the distance between the anticodons of two transfer RNAs simultaneously bound to a messengerprogramed Escherichia coli 70 S ribosome. The fluorescent Y base adjacent to the anticodon of yeast tRNA_Y^Phe serves as a donor. A proflavine (Pf) chemically substituted for the Y base in tRNA_Pf^Phe serves as an acceptor. By exploiting the sequential binding properties of 70 S ribosomes for two deacylated tRNAs, we can fill the strong site with either tRNA_Y^Phe or tRNA_Pf^Phe and then the weak site with the other tRNA. In both cases donor quenching and sensitized emission of the acceptor are observed. Analysis of these results leads to an estimate for the Y-proflavine distance of 18 ± 2 Å. This distance is very short and suggests strongly that the two tRNAs are simultaneously in contact with adjacent codons of the message. Separate experiments show that binding of a tRNA to the weak site does not perturb the environment of the hypermodified base of a tRNA bound to the strong site. This supports the assignment of the strong site as the peptidyl site. It also indicates that binding of the second tRNA proceeds without a change in the anticodon structure of a pre-existing tRNA at the peptidyl site.     

Bioactivity of folding intermediates studied by the recovery of enzymatic activity during refolding

The peptide bond preceding proline residues realizes a cis/trans conformational switch with high switching resistance in native proteins and folding intermediates. Therefore, individual isomers have the potential to differ in bioactivity. However, information about isomer-specific bioactivities is difficult to obtain because of the risk of affecting isomeric distribution by bioactivity assay components.Here we present an approach that allows for the measurement of the recovery of enzymatic activities of wild-type RNase T₁ and RNase T₁ variants during refolding under conditions where the population of enzyme-substrate or enzyme-product complexes is negligible. Recovery of enzymatic activity was continuously monitored within the visible range of the spectrum by addition of a fluorescence-labeled nucleotide substrate to the refolding sample. We found that a nonnative trans conformation at Pro39 renders the RNase T₁ almost completely inactive. A folding intermediate having a nonnative trans conformation at Pro55 shows about 46% of the enzymatic activity referred to the native state. Pro55, in contrast to the active site located Pro39, is situated in a solvent-exposed loop region remote from active-site residues. In both cases, peptidyl prolyl cis/trans isomerases accelerate the regain of nucleolytic activity. Our findings show that even if there is a considerable distance between the site of isomerization and the active site, conformational control of the bioactivity of proteins is likely to occur, and that the surface location of prolyl bonds suffices for the control of buried active sites mediated by peptidyl prolyl cis/trans isomerases.     

Unveiling the Pathways of Dioxygen Through the C2 Component of the Environmentally Relevant Monooxygenase p‐Hydroxyphenylacetate Hydroxylase from Acinetobacter baumannii: A Molecular Dynamics Investigation

In this work, models of the homotetrameric C2 component of the monooxygenase p‐hydroxyphenylacetate hydroxylase from Acinetobacter baumannii, in complex with dioxygen (O₂) and, or not, the substrate p‐hydroxyphenylacetate (HPA) were built. Both models proved to be amenable to random‐acceleration molecular dynamics (RAMD) simulations, whereby a tiny randomly oriented external force, acting on O₂ at the active site in front of flavin mononucleotide (FMNH^?), accelerated displacement of O₂ toward the bulk solvent. This allowed us to carry out a sufficiently large number of RAMD simulations to be of statistical significance. The two systems behaved very similarly under RAMD, except for O₂ leaving the active site more easily in the absence of HPA, but then finding similar obstacles in getting to the gate as when the active site was sheltered by HPA. This challenges previous conclusions that HPA can only reach the active center after that the C4aOOH derivative of FMNH^? is formed, requiring uptake of O₂ at the active site before HPA. According to these RAMD simulations, O₂ could well get to FMNH^? also in the presence of the substrate at the active site.     

Rotational resonance NMR study of the active site structure in bacteriorhodopsin: conformation of the Schiff base linkage.

Rotational resonance, a new solid-state NMR technique for determining internuclear distances, is used to measure a distance in the active site of bacteriorhodopsin (bR) that changes in different states of the protein. The experiments are targeted to the active site of bR through 13C labeling of both the retinal chromophore and the Lys side chains of the protein. The time course of the rotor-driven magnetization exchange between a pair of 13C nuclei is then observed to determine the dipolar coupling and therefore the internuclear distance. Using this approach, we have measured the distance from [14-13C]retinal to [epsilon-13C]Lys216 in dark-adapted bR in order to examine the structure of the retinal-protein linkage and its role in coupling the isomerizations of retinal to unidirectional proton transfer. This distance depends on the configuration of the intervening C=N bond. The 3.0 +/- 0.2 A distance observed in bR555 demonstrates that the C=N bond is syn, and the 4.1 +/- 0.3 A distance observed in bR568 demonstrates that the C=N bond is anti. These direct distance determinations independently confirm the configurations previously deduced from solid-state NMR chemical shift and resonance Raman vibrational spectra. The spectral selectivity of rotational resonance allows these two distances to be measured independently in a sample containing both bR555 and bR568; the presence of both states and of 25% lipid in the sample demonstrates the use of rotational resonance to measure an active site distance in a membrane protein with an effective molecular mass of about 85 kDa.(ABSTRACT TRUNCATED AT 250 WORDS)     

Crystal structure of vespid phospholipase A1 reveals insights into the mechanism for cause of membrane dysfunction

Vespid phospholipase A₁ (vPLA₁) from the black-bellied hornet (Vespa basalis) catalyzes the hydrolysis of emulsified phospholipids and shows potent hemolytic activity that is responsible for its lethal effect. To investigate the mechanism of vPLA₁ towards its function such as hemolysis and emulsification, we isolated vPLA₁ from V. basalis venom and determined its crystal structure at 2.5 Å resolution. vPLA₁ belongs to the α/β hydrolase fold family. It contains a tightly packed β-sheet surrounded by ten α-helices and a Gly-X-Ser-X-Gly motif, characteristic of a serine hydrolyase active site. A bound phospholipid was modeled into the active site adjacent to the catalytic Ser-His-Asp triad indicating that Gln95 is located at hydrogen-bonding distance from the substrate's phosphate group. Moreover, a hydrophobic surface comprised by the side chains of Phe53, Phe62, Met91, Tyr99, Leu197, Ala167 and Pro169 may serve as the acyl chain-binding site. vPLA₁ shows global similarity to the N-terminal domain of human pancreatic lipase (HPL), but with some local differences. The lid domain and the β9 loop responsible for substrate selectivity in vPLA₁ are shorter than in HPL. Thus, solvent-exposed hydrophilic residues can easily accommodate the polar head groups of phospholipids, thereby accounting for the high activity level of vPLA₁. Our result provides a potential explanation for the ability of vPLA₁ to hydrolyze phospholipids of cell membrane.     

Rabbit liver fructose 1,6-bisphosphatase: labeling of the active and allosteric sites with pyridoxal 5-phosphate and sequence of a nonapeptide from the active site

The active and allosteric sites of fructose 1,6-bisphosphatase (Fru-P₂ase, EC 3.1.3.11) were labeled by reaction with pyridoxal phosphate and sodium borohydride in the presence of the allosteric inhibitor AMP or the substrate, Fru-P₂, respectively. Modification of the active site results in loss of activity. Modification of the allosteric site decreases the sensitivity of the enzyme to inhibition by AMP and alters its ability to bind to blue dextran-Sepharose. The allosteric and active sites have been located on different cyanogen bromide peptides; the sequence of a nonapeptide from the active site is (H)GlyLysLeuArgLeuLeu TyrGluCys(OH). The lysyl residue is modified by pyridoxal phosphate.     

Structural Evidence for a Copper-Bound Carbonate Intermediate in the Peroxidase and Dismutase Activities of Superoxide Dismutase

Copper-zinc superoxide dismutase (SOD) is of fundamental importance to our understanding of oxidative damage. Its primary function is catalysing the dismutation of superoxide to O₂ and H₂O₂. SOD also reacts with H₂O₂, leading to the formation of a strong copper-bound oxidant species that can either inactivate the enzyme or oxidise other substrates. In the presence of bicarbonate (or CO₂) and H₂O₂, this peroxidase activity is enhanced and produces the carbonate radical. This freely diffusible reactive oxygen species is proposed as the agent for oxidation of large substrates that are too bulky to enter the active site. Here, we provide direct structural evidence, from a 2.15 Å resolution crystal structure, of (bi)carbonate captured at the active site of reduced SOD, consistent with the view that a bound carbonate intermediate could be formed, producing a diffusible carbonate radical upon reoxidation of copper. The bound carbonate blocks direct access of substrates to Cu(I), suggesting that an adjunct to the accepted mechanism of SOD catalysed dismutation of superoxide operates, with Cu(I) oxidation by superoxide being driven via a proton-coupled electron transfer mechanism involving the bound carbonate rather than the solvent. Carbonate is captured in a different site when SOD is oxidised, being located in the active site channel adjacent to the catalytically important Arg143. This is the probable route of diffusion from the active site following reoxidation of the copper. In this position, the carbonate is poised for re-entry into the active site and binding to the reduced copper.     

Effects of protonation state of Asp181 and position of active site water molecules on the conformation of PTP1B

In protein tyrosine phosphatase 1B (PTP1B), the flexible WPD loop adopts a closed conformation (WPD_closed) in the active state of PTP1B, bringing the catalytic Asp181 close to the active site pocket, while WPD loop is in an open conformation (WPD_open) in the inactive state. Previous studies showed that Asp181 may be protonated at physiological pH, and ordered water molecules exist in the active site. In the current study, molecular dynamics simulations are employed at different Asp181 protonation states and initial positions of active site water molecules, and compared with the existing crystallographic data of PTP1B. In WPD_closed conformation, the active site is found to maintain its conformation only in the protonated state of Asp181 in both free and liganded states, while Asp181 is likely to be deprotonated in WPD_open conformation. When the active site water molecule network that is a part of the free WPD_closed crystal structure is disrupted, intermediate WPD loop conformations, similar to that in the PTPRR crystal structure, are sampled in the MD simulations. In liganded PTP1B, one active site water molecule is found to be important for facilitating the orientation of Cys215 and the phosphate ion, thus may play a role in the reaction. In conclusion, conformational stability of WPD loop, and possibly catalytic activity of PTP1B, is significantly affected by the protonation state of Asp181 and position of active site water molecules, showing that these aspects should be taken into consideration both in MD simulations and inhibitor design. © Proteins 2013. © 2012 Wiley Periodicals, Inc.     

嗜热丁烷氧化古菌Candidatus Syntrophoarchaeum烷基转移酶MtaA催化丁基辅酶M的分子机制研究

【背景】嗜热古菌Candidatus Syntrophoarchaeum可以与硫酸盐还原细菌共生,通过逆转产甲烷途径进行正丁烷的氧化,但在该过程中负责催化丁基辅酶M氧化的酶尚未确定。【目的】利用分子动力学模拟证明Ca.Syntrophoarchaeum中mta A基因编码的蛋白可以特异性催化丁基辅酶M中丁基的转移,并非转移甲基。【方法】使用Methanosarcina mazei辅酶M甲基转移酶Mta A的晶体结构(PDB ID:4ay8)作为模板,对Mta A_1 (Gen Bank登录号OFV65993.1)和Mta A_2 (Gen Bank登录号OFV65678.1)进行同源建模。使用分子对接得到两者分别结合CH_3-Co M和C_4H_9-Co M时的结构,并用AMBER18进行分子动力学模拟。【结果】当Mta A_1和Mta A_2分别结合C_4H_9-Co M时,表现出与4ay8晶体结构类似的TIM-Barrel折叠三维结构,但在活性中心形状、Zn~(2+)与底物距离以及活性位点附近氨基酸配位方式等方面存在差异,这可能是导致Ca.Syntrophoarchaeum中mta A基因编码的蛋白催化丁基辅酶M氧化的原因。其中Mta A_2与4ay8结构更相似,活性中心氨基酸配位更完整,暗示其更可能具备催化活性。然而当Mta A_1和Mta A_2分别结合CH_3-Co M时,整体结构不合实际,活性中心Zn~(2+)与底物距离过远,表明底物几乎不可能与酶结合。【结论】Ca.Syntrophoarchaeum中的Mta A_1和Mta A_2很可能是特异性的丁基转移酶,而非催化甲基的转移,其中Mta A_2具备活性的可能性更高。     

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.     

Simulation of structural and functional properties of mevalonate diphosphate decarboxylase (MVD)

Mevalonate 5-diphosphate decarboxylase (MVD) is an important enzyme in the mevalonate pathway catalyzing the ATP-dependent decarboxylation of mevalonate 5-diphosphate (MDP) to yield isopentynyl diphosphate (IPP) which is an ubiquitous precursor for isoprenoids and sterols. Although there are studies to show the involvement of certain amino acid residues in MVD activity, the structure and the function of the active site is yet to be investigated. Therefore the objectives of this study were to elucidate the active site of Saccharomyces cerevisiae MVD (scMVD) using a molecular docking and simulation-based approach. The Cartesian coordinates of scMVD retrieved from the PDB database were used in the docking procedure. 3D atomic coordinates of MDP, ATP and an inhibitor trifluoromevalonate (TFMDP) were generated using Gaussian 98. ATP, MDP and TFMDP were docked into the potential active site identified by sequence analyses using Hex 4.2. The complexes obtained from docking procedure were subjected to 1.5 ns simulation by GROMACS 3.2. Investigation of complexes revealed that Ala15, Lys18, Ser121 &; Ser155; Lys22, Ser153 &; Ser155 and Tyr19, Ser121, Ser153, Gly154 &; Thr209 of MVD are within hydrogen bond forming distances of MDP, ATP and TFMDP, respectively indicating their possible involvement in active site formation through H-bond formation. The presence of a water molecule between the carboxyl group of Asp302, a previously characterized active site residue and C3 region of MDP at a distance of 3 Å suggests that deprotonation of the hydroxyl of the C3 takes place via a water molecule. Conjunction with reported crucial catalytic activity of Ser121 of MVD and our finding of the presence of this residue in hydrogen bond forming distance to MDP suggests that this hydrogen bond helps in proper orienting of MDP for phosphorylation /decarboxylation. We further suggest that the reported greater RMS deviation of Pro₇₉- Leu mutated MVD with respect to native MVD of temperature sensitive mutant phenotype of S. cerevisiae is due to partial unfolding of MVD as a result of mutation. Finally, this study provides a tantalizing glimpse about hitherto unknown structural and functional properties of the active site of MVD.