Active Site Water Molecules Research Articles

The β-subunit of tryptophan synthase is a latent tyrosine synthase.

Aromatic amino acids and their derivatives are diverse primary and secondary metabolites with critical roles in protein synthesis, cell structure and integrity, defense and signaling. All de novo aromatic amino acid production relies on a set of ancient and highly conserved chemistries. Here we introduce a new enzymatic transformation for L-tyrosine synthesis by demonstrating that the β-subunit of tryptophan synthase-which natively couples indole and L-serine to form L-tryptophan-can act as a latent 'tyrosine synthase'. A single substitution of a near-universally conserved catalytic residue unlocks activity toward simple phenol analogs and yields exclusive para carbon-carbon bond formation to furnish L-tyrosines. Structural and mechanistic studies show how a new active-site water molecule orients phenols for a nonnative mechanism of alkylation, with additional directed evolution resulting in a net >30,000-fold rate enhancement. This new biocatalyst can be used to efficiently prepare valuable L-tyrosine analogs at gram scales and provides the missing chemistry for a conceptually different pathway to L-tyrosine.

Structural and biochemical characterization of a nucleotide hydrolase from Streptococcus pneumonia

In this report, we structurally and biochemically characterized the unknown gene product SP1746 from Streptococcus pneumoniae serotype 4. Various crystal structures of SP1746 in the apo form and in complex with different nucleotides were determined. SP1746 is a globular protein, which belongs to the histidine-aspartate (HD) domain superfamily with two Fe3+ ions in the active site that are coordinated by key active site residues and water molecules. All nucleotides bind in a similar orientation in the active site with their phosphate groups anchored to the diiron cluster. Biochemically, SP1746 hydrolyzes different nucleotide substrates. SP1746 most effectively hydrolyzes diadenosine tetraphosphate (Ap4A) to two ADPs. Based on the aforementioned data, we annotated SP1746 as an Ap4A hydrolase, belonging to the YqeK family. Our invitro data indicate a potential role for SP1746 in regulating Ap4A homeostasis, which requires validation with invivo experiments in bacteria in the future.

PH-Dependent Structure and Dynamics of the Catalytic Domains of Human Carbonic Anhydrase II and IX.

Extensive computer simulation studies have been carried out to probe the pH-dependent structure and dynamics of the two most efficient isoenzymes II and IX of human carbonic anhydrase (HCA) that control the pH in the human body. The equilibrium structure and hydration of their catalytic domains are found to be largely unaffected by the variation of pH in the range studied, in close agreement with the known experimental results. In contrast, a significant effect of the change in pH is observed for the first time on the local electrostatic potential of the active site walls and the dynamics of active site water molecules. We also report for the first time the free energy and kinetics of coupled fluctuations of orientation and protonation states of the well-known His-mediated proton shuttle (His-64) in both isozymes at pH 7 and 8. The transitions between different tautomers of in or out conformations of His-64 side chain range between 109 and 106 s-1 depending on pH. Possible implications of these results on conformation-dependent pKa of His-64 side chain and its role in driving the catalysis toward hydration of CO2 or dehydration of HCO3- with varying pH are discussed.

Computational Insights on the Hydride and Proton Transfer Mechanisms of D-Arginine Dehydrogenase.

D-Arginine dehydrogenase from Pseudomonas aeruginosa (PaDADH) is an amine oxidase which catalyzes the conversion of D-arginine into iminoarginine. It contains a non-covalent FAD cofactor that is involved in the oxidation mechanism. Based on substrate, solvent, and multiple kinetic isotope effects studies, a stepwise hydride transfer mechanism is proposed. It was shown that D-arginine binds to the active site of enzyme as α-amino group protonated, and it is deprotonated before a hydride ion is transferred from its α-C to FAD. Based on a mutagenesis study, it was concluded that a water molecule is the most likely catalytic base responsible from the deprotonation of α-amino group. In this study, we formulated computational models based on ONIOM method to elucidate the oxidation mechanism of D-arginine into iminoarginine using the crystal structure of enzyme complexed with iminoarginine. The calculations showed that Arg222, Arg305, Tyr249, Glu87, His 48, and two active site water molecules play key roles in binding and catalysis. Model systems showed that the deprotonation step occurs prior to hydride transfer step, and active site water molecule(s) may have participated in the deprotonation process.

Structural insights into the substrate recognition of serine palmitoyltransferase from Sphingobacterium multivorum

Serine palmitoyltransferase (SPT) is a key enzyme of sphingolipid biosynthesis, which catalyzes the pyridoxal-5'-phosphate-dependent decarboxylative condensation reaction of L-serine (L-Ser) and palmitoyl-CoA (PalCoA) to form 3-ketodihydrosphingosine called long chain base (LCB). SPT is also able to metabolize L-alanine (L-Ala) and glycine (Gly), albeit with much lower efficiency. Human SPT is a membrane-bound large protein complex containing SPTLC1/SPTLC2 heterodimer as the core subunits, and it is known that mutations of the SPTLC1/SPTLC2 genes increase the formation of deoxy-type of LCBs derived from L-Ala and Gly to cause some neurodegenerative diseases. In order to study the substrate recognition of SPT, we examined the reactivity of Sphingobacterium multivorum SPT on various amino acids in the presence of PalCoA. The S. multivorum SPT could convert not only L-Ala and Gly but also L-homoserine, in addition to L-Ser, into the corresponding LCBs. Furthermore, we obtained high-quality crystals of the ligand-free form and the binary complexes with a series of amino acids, including a nonproductive amino acid, L-threonine, and determined the structures at 1.40-1.55 Å resolutions. The S. multivorum SPT accommodated various amino acid substrates through subtle rearrangements of the active-site amino acid residues and water molecules. It was also suggested that non-active-site residues mutated in the human SPT genes might indirectly influence the substrate specificity by affecting the hydrogen-bonding networks involving the bound substrate, water molecules, and amino acid residues in the active site of this enzyme. Collectively, our results highlight SPT structural features affecting substrate specificity for this stage of sphingolipid biosynthesis.

Mutagenic analysis of actin reveals the mechanism of His161 flipping that triggers ATP hydrolysis.

The dynamic assembly of actin is controlled by the hydrolysis of ATP, bound to the center of the molecule. Upon polymerization, actin undergoes a conformational change from the monomeric G-form to the fibrous F-form, which is associated with the flipping of the side chain of His161 toward ATP. His161 flipping from the gauche-minus to gauche-plus conformation leads to a rearrangement of the active site water molecules, including ATP attacking water (W1), into an orientation capable of hydrolysis. We previously showed that by using a human cardiac muscle α-actin expression system, mutations in the Pro-rich loop residues (A108G and P109A) and in a residue that was hydrogen-bonded to W1 (Q137A) affect the rate of polymerization and ATP hydrolysis. Here, we report the crystal structures of the three mutant actins bound to AMPPNP or ADP-Pi determined at a resolution of 1.35-1.55Å, which are stabilized in the F-form conformation with the aid of the fragmin F1 domain. In A108G, His161 remained non-flipped despite the global actin conformation adopting the F-form, demonstrating that the side chain of His161 is flipped to avoid a steric clash with the methyl group of A108. Because of the non-flipped His161, W1 was located away from ATP, similar to G-actin, which was accompanied by incomplete hydrolysis. In P109A, the absence of the bulky proline ring allowed His161 to be positioned near the Pro-rich loop, with a minor influence on ATPase activity. In Q137A, two water molecules replaced the side-chain oxygen and nitrogen of Gln137 almost exactly at their positions; consequently, the active site structure, including the W1 position, is essentially conserved. This seemingly contradictory observation to the reported low ATPase activity of the Q137A filament could be attributed to a high fluctuation of the active site water. Together, our results suggest that the elaborate structural design of the active site residues ensures the precise control of the ATPase activity of actin.

A Multiple Proton Transfer Mechanism for the Charging Step of the Aminoacylation Reaction at the Active Site of Aspartyl tRNA Synthetase.

Aspartyl-tRNA synthetase catalyzes the attachment of aspartic acid to its cognate tRNA by the aminoacylation reaction during the initiation of the protein biosynthesis process. In the second step of the aminoacylation reaction, known as the charging step, the aspartate moiety is transferred from aspartyl-adenylate to the 3'-OH of A76 of tRNA through a proton transfer process. We have investigated different pathways for the charging step through three separate QM/MM simulations combined with the enhanced sampling method of well-sliced metadynamics and found out the most feasible pathway for the reaction at the active site of the enzyme. In the charging reaction, both the phosphate group and the ammonium group after deprotonation can potentially act as a base for proton transfer in the substrate-assisted mechanism. We have considered three possible mechanisms involving different pathways of proton transfer, and only one of them is determined to be enzymatically feasible. The free energy landscape along reaction coordinates where the phosphate group acts as the general base showed that, in the absence of water, the barrier height is 52.6 kcal/mol. The free energy barrier is reduced to 39.7 kcal/mol when the active site water molecules are also treated quantum mechanically, thus allowing a water mediated proton transfer. The charging reaction involving the ammonium group of the aspartyl adenylate is found to follow a path where first a proton from the ammonium group moves to a water in the vicinity forming a hydronium ion (H3O+) and NH2 group. The hydronium ion subsequently passes the proton to the Asp233 residue, thus minimizing the chance of back proton transfer from hydronium to the NH2 group. The neutral NH2 group subsequently takes the proton from the O3' of A76 with a free energy barrier of 10.7 kcal/mol. In the next step, the deprotonated O3' makes a nucleophilic attack to the carbonyl carbon forming a tetrahedral transition state with a free energy barrier of 24.8 kcal/mol. Thus, the present work shows that the charging step proceeds through a multiple proton transfer mechanism where the amino group formed after deprotonation acts as the base to capture a proton from O3' of A76 rather than the phosphate group. The current study also shows the important role played by Asp233 in the proton transfer process.

Insight into the Role of Active Site Protonation States and Water Molecules in the Catalytic Inhibition of DPP4 by Vildagliptin.

Vildagliptin (VIL) is an antidiabetic drug that inhibits dipeptidyl peptidase-4 (DPP4) through a covalent mechanism. The molecular bases for this inhibitory process have been addressed experimentally and computationally. Nevertheless, relevant issues remain unknown such as the roles of active site protonation states and conserved water molecules nearby the catalytic center. In this work, molecular dynamics simulations were applied to examine the structures of 12 noncovalent VIL-DPP4 complexes encompassing all possible protonation states of three noncatalytic residues (His126, Asp663, Asp709) that were inconclusively predicted by different computational tools. A catalytically competent complex structure was only achieved in the system with His126 in its ε-form and nonconventional neutral states for Asp663/Asp709. This complex suggested the involvement of one water molecule in the catalytic process of His740/Ser630 activation, which was confirmed by QM/MM simulations. Our findings support the suitability of a novel water-mediated mechanism in which His740/Ser630 activation occurs concertedly with the nucleophilic attack on VIL and the imidate protonation by Tyr547. Then, the restoration of His740/ Tyr547 protonation states occurs via a two-water hydrogen bonding network in a low-barrier process, thus describing the final step of the catalytic cycle for the first time. Additionally, two hydrolytic mechanisms were proposed based on the hydrogen bonding networks formed by water molecules and the catalytic residues along the inhibitory mechanism. These findings are valuable to unveil the molecular features of the covalent inhibition of DPP4 by VIL and support the future development of novel derivatives with improved structural or mechanistic profiles.

Probing the general base for DNA polymerization in telomerase: a molecular dynamics investigation.

Telomerase is an RNA-dependent DNA polymerase that plays a role in the maintenance of the 3' end of the eukaryotic chromosome, known as a telomere, by catalyzing the DNA polymerization reaction in cancer and embryonic stem cells. The detailed molecular details of the DNA polymerization by telomerase, especially the general base for deprotonating the terminal 3'-hydroxyl, which triggers the chemical reaction, remain elusive. We conducted a computational investigation using hybrid quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) simulations to probe the detailed mechanism of the reaction. Our simulations started with the telomerase:RNA:DNA:dNTP ternary complex, and by using enhanced sampling QM/MM MD simulations, we probed the general base involved directly in the polymerization. We report the participation of an aspartate (Asp344) coordinated to Mg and an active site water molecule, jointly acting as a base during nucleic acid addition. The Asp344 residue remains transiently protonated during the course of the reaction, and later it deprotonates by transferring its proton to the water at the end of the reaction.

Structural and Energetic Origin of Different Product Specificities and Activities for SETD3 and Its Mutants on the Methylation of the β-Actin H73K Peptide: Insights from a QM/MM Study.

The methylation of the lysine residue can affect some fundamental biological processes, and specific biological effects of the methylations are often related to product specificity of methyltransferases. The question remains concerning how active-site structural features and dynamics control the activity as well as the number (1, 2, or 3) of methyl groups on methyl lysine products. SET domain containing protein 3 (SETD3) has been identified recently as the β-actin histidine73-N3 methyltransferase, and also, it has a weak methylation activity on the H73K β-actin peptide for which the target H73 residue is mutated into K73. Interestingly, the K73 methylation activity of SETD3 increases significantly as a result of the N255 → A or N255 → F/W273 → A mutation, and the N255A product specificity also differs from that of wild-type. Here, we performed QM/MM molecular dynamics and potential of mean force (PMF) simulations for SETD3 and its mutants (N255A and N255F/W273A) to study how SETD3 and its mutants could have different product specificities and activities for the K73 methylation. The PMF simulations show that the barrier for the first methylation of K73 is higher compared to the barrier of the H73 methylation in SETD3. Moreover, the second methylation of K73 has been found to have a barrier from the free energy simulation that is higher by 2.2 kcal/mol compared to the barrier of the first methyl transfer to K73, agreeing with the suggestion that SETD3 is a monomethylase. For the first, second, and third methylations of K73 in the N255A mutant, the barriers obtained from the PMF simulations for transferring the second and third methyl groups are found to be lower relative to the barrier for the first methyl transfer. Thus, N255A can be considered as a trimethyl lysine methyltransferase. In addition, for the first K73 methylation, the activities from the PMF simulations follow the order of N255F/W273A > N255A > WT, in agreement with experiments. The examination of the structural and dynamic results at the active sites provides better understanding of different product specificities and activities for the K73 methylations in SETD3 and its mutants. It is demonstrated that the existence of well-balanced interactions at the active site leading to the near attack conformation is of crucial importance for the efficient methyl transfers. Moreover, the presence of potential interactions (e.g., the C-H···O and cation-π interactions) that are strengthening at the transition state can also be important. Furthermore, the activity as well as product specificity of the K73 methylation also seems to be controlled by certain active-site water molecules which may be released to provide extra space for the addition of more methyl groups on K73.

Monofunctional Heme-Catalases.

The review focuses on four issues that are critical for the understanding of monofunctional catalases. How hydrogen peroxide (H2O2) reaches the active site and outcompetes water molecules to be able to function at a very high rate is one of the issues examined. Part of the answer is a gate valve system that is instrumental to drive out solvent molecules from the final section of the main channel. A second issue relates to how the enzyme deals with an unproductive reactive compound I (Cpd I) intermediate. Peroxidatic two and one electron donors and the transfer of electrons to the active site from NADPH and other compounds are reviewed. The new ascribed catalase reactions are revised, indicating possible measurement pitfalls. A third issue concerns the heme b to heme d oxidation, why this reaction occurs only in some large-size subunit catalases (LSCs), and the possible role of singlet oxygen in this and other modifications. The formation of a covalent bond between the proximal tyrosine with the vicinal residue is analyzed. The last issue refers to the origin and function of the additional C-terminal domain (TD) of LSCs. The TD has a molecular chaperone activity that is traced to a gene fusion between a Hsp31-type chaperone and a small-size subunit catalase (SSC).

Absolute Binding Free Energy Calculations for Buried Water Molecules.

Water often plays a key role in mediating protein-ligand interactions. Understanding contributions from active-site water molecules to binding thermodynamics of a ligand is important in predicting binding free energies for ligand optimization. In this work, we tested a non-equilibrium switching method for absolute binding free energy calculations on water molecules in binding sites of 13 systems. We discuss the lessons we learned about identified issues that affected our calculations and ways to address them. This work fits with our larger focus on how to do accurate ligand binding free energy calculations when water rearrangements are very slow, such as rearrangements due to ligand modification (as in relative free energy calculations) or ligand binding (as in absolute free energy calculations). The method studied in this work can potentially be used to account for limited water sampling via providing endpoint corrections to free energy calculations using our calculated binding free energy of water.

Elucidation of the tyrosinase/O2/monophenol ternary intermediate that dictates the monooxygenation mechanism in melanin biosynthesis

Melanins are highly conjugated biopolymer pigments that provide photoprotection in a wide array of organisms, from bacteria to humans. The rate-limiting step in melanin biosynthesis, which is the ortho-hydroxylation of the amino acid L-tyrosine to L-DOPA, is catalyzed by the ubiquitous enzyme tyrosinase (Ty). Ty contains a coupled binuclear copper active site that binds O2 to form a μ:η2:η2-peroxide dicopper(II) intermediate (oxy-Ty), capable of performing the regioselective monooxygenation of para-substituted monophenols to catechols. The mechanism of this critical monooxygenation reaction remains poorly understood despite extensive efforts. In this study, we have employed a combination of spectroscopic, kinetic, and computational methods to trap and characterize the elusive catalytic ternary intermediate (Ty/O2/monophenol) under single-turnover conditions and obtain molecular-level mechanistic insights into its monooxygenation reactivity. Our experimental results, coupled with quantum-mechanics/molecular-mechanics calculations, reveal that the monophenol substrate docks in the active-site pocket of oxy-Ty fully protonated, without coordination to a copper or cleavage of the μ:η2:η2-peroxide O-O bond. Formation of this ternary intermediate involves the displacement of active-site water molecules by the substrate and replacement of their H bonds to the μ:η2:η2-peroxide by a single H bond from the substrate hydroxyl group. This H-bonding interaction in the ternary intermediate enables the unprecedented monooxygenation mechanism, where the μ-η2:η2-peroxide O-O bond is cleaved to accept the phenolic proton, followed by substrate phenolate coordination to a copper site concomitant with its aromatic ortho-hydroxylation by the nonprotonated μ-oxo. This study provides insights into O2 activation and reactivity by coupled binuclear copper active sites with fundamental implications in biocatalysis.

Different Enzyme Conformations Induce Different Mechanistic Traits in HIV-1 Protease.

The influence of the dynamical flexibility of enzymes on reaction mechanisms is a cornerstone in biological sciences. In this study, we aim to 1) study the convergence of the activation free energy by using the first step of the reaction catalysed by HIV-1 protease as a case study, and 2) provide further evidence for a mechanistic divergence in this enzyme, as two different reaction pathways were seen to contribute to this step. We used quantum mechanics/molecular mechanics molecular dynamics simulations, on four different initial conformations that led to different barriers in a previous study. Despite the sampling, the four activation free energies still spanned a range of 5.0 kcal ⋅ mol-1 . Furthermore, the new simulations did confirm the occurrence of an unusual mechanistic divergence, with two different mechanistic pathways displaying equivalent barriers. An active-site water molecule is proposed to influence the mechanistic pathway.

The mechanisms of catalysis and ligand binding for the SARS-CoV-2 NSP3 macrodomain from neutron and x-ray diffraction at room temperature.

The nonstructural protein 3 (NSP3) macrodomain of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) (Mac1) removes adenosine diphosphate (ADP) ribosylation posttranslational modifications, playing a key role in the immune evasion capabilities of the virus responsible for the coronavirus disease 2019 pandemic. Here, we determined neutron and x-ray crystal structures of the SARS-CoV-2 NSP3 macrodomain using multiple crystal forms, temperatures, and pHs, across the apo and ADP-ribose–bound states. We characterize extensive solvation in the Mac1 active site and visualize how water networks reorganize upon binding of ADP-ribose and non-native ligands, inspiring strategies for displacing waters to increase the potency of Mac1 inhibitors. Determining the precise orientations of active site water molecules and the protonation states of key catalytic site residues by neutron crystallography suggests a catalytic mechanism for coronavirus macrodomains distinct from the substrate-assisted mechanism proposed for human MacroD2. These data provoke a reevaluation of macrodomain catalytic mechanisms and will guide the optimization of Mac1 inhibitors.

Structure-affinity relationships of reversible proline analog inhibitors targeting proline dehydrogenase.

Proline dehydrogenase (PRODH) catalyzes the first step of proline catabolism, the FAD-dependent oxidation of L-proline to Δ1-pyrroline-5-carboxylate. PRODH plays a central role in the metabolic rewiring of cancer cells, which has motivated the discovery of inhibitors. Here, we studied the inhibition of PRODH by 18 proline-like compounds to understand the structural and chemical features responsible for the affinity of the best-known inhibitor, S-(-)-tetrahydro-2-furoic acid (1). The compounds were screened, and then six were selected for more thorough kinetic analysis: cyclobutane-1,1-dicarboxylic acid (2), cyclobutanecarboxylic acid (3), cyclopropanecarboxylic acid (4), cyclopentanecarboxylic acid (16), 2-oxobutyric acid (17), and (2S)-oxetane-2-carboxylic acid (18). These compounds are competitive inhibitors with inhibition constants in the range of 1.4-6 mM, compared to 0.3 mM for 1. Crystal structures of PRODH complexed with 2, 3, 4, and 18 were determined. All four inhibitors bind in the proline substrate site, but the orientations of their rings differ from that of 1. The binding of 3 and 18 is accompanied by compression of the active site to enable nonpolar contacts with Leu513. Compound 2 is unique in that the additional carboxylate displaces a structurally conserved water molecule from the active site. Compound 18 also destabilizes the conserved water, but by an unexpected non-steric mechanism. The results are interpreted using a chemical double mutant thermodynamic cycle. This analysis revealed unanticipated synergism between ring size and hydrogen bonding to the conserved water. These structure-affinity relationships provide new information relevant to the development of new inhibitor design strategies targeting PRODH.

Conformational exploration of RbgA using molecular dynamics: Possible implications in ribosome maturation and GTPase activity in different nucleotide bound states

Ribosome biogenesis GTPase A (RbgA) is involved in the late steps of the 50S ribosomal subunit maturation by binding into the 45S pre-ribosomal subunit. The association of RbgA to the 45S intermediate subunit depends on its bound nucleotide (GTP/GDP), probably because of the conformational shifts that occur between the GTP and GDP bound states. However, the available crystal structures of Staphylococcus aureus RbgA (SaRbgA) do not show any significant variations between different nucleotide bound states. Therefore, conformational exploration of SaRbgA in different nucleotide bound states was carried out using all-atom molecular dynamics (MD) simulations. Exploration of conformational distribution using cluster analysis and principal component analysis (PCA) revealed that GDP and pppGpp bound systems exhibit a larger distribution. This is majorly due to the fluctuations of the C-terminal tail (C-tail) as a result of the unwinding of α-helical secondary conformations into loop conformations which are observed from RMSF and DSSP analyses. Further investigation of the network of interactions revealed that the GTP and GMPPNP bound systems hold the C-tail in an α-helical form through stronger interactions between the active-site and C-tail. We also find that the presence of Mg2+ positions Sw-I loop away from the bound nucleotide and stabilizes the active-site water molecules. This seems to assist SaRbgA GTPase activity. In addition, mutations at the C-terminal and Sw-II conserved residues exhibit a larger conformational distribution majorly due to the C-tail fluctuations suggesting that the C-tail of SaRbgA probably interacts with the rRNA or rprotein in the process of ribosome biogenesis.

The basicity of an active-site water molecule discriminates between tyrosinase and catechol oxidase activity

Tyrosinase (Ty) and catechol oxidase (CO) are members of type-3 copper enzymes. While Ty catalyzes both phenolase and catecholase reactions, CO catalyzes only the latter reaction. In the present study, Ty was found to catalyze the catecholase reaction, but hardly the phenolase reaction in the presence of the metallochaperon called “caddie protein (Cad)”. The ability of the substrates to dissociate the motif shielding the active-site pocket seems to contribute critically to the substrate specificity of Ty. In addition, a mutation at the N191 residue, which forms a hydrogen bond with a water molecule near the active center, decreased the inherent ratio of phenolase versus catecholase activity. Unlike the wild-type complex, reaction intermediates were not observed when the catalytic reaction toward the Y98 residue of Cad was progressed in the crystalline state. The increased basicity of the water molecule may be necessary to inhibit the proton transfer from the conjugate acid to a hydroxide ion bridging the two copper ions. The deprotonation of the substrate hydroxyl by the bridging hydroxide seems to be significant for the efficient catalytic cycle of the phenolase reaction.

Assessment of CYP2C9 Structural Models for Site of Metabolism Prediction.

Structure-based prediction of a compound's potential sites of metabolism (SOMs) mediated by cytochromes P450 (CYPs) is highly advantageous in the early stage of drug discovery. However, the accuracy of the SOMs prediction can be influenced by several factors. CYP2C9 is one of the major drug-metabolizing enzymes in humans and is responsible for the metabolism of ∼13 % of clinically used drugs. In this study, we systematically evaluated the effects of protein crystal structure models, scoring functions, heme forms, conserved active-site water molecules, and protein flexibility on SOMs prediction of CYP2C9 substrates. Our results demonstrated that, on average, ChemScore and GlideScore outperformed four other scoring functions: Vina, GoldScore, ChemPLP, and ASP. The performance of the crystal structure models with pentacoordinated heme was generally superior to that of the hexacoordinated iron-oxo heme (referred to as Compound I) models. Inclusion of the conserved active-site water molecule improved the prediction accuracy of GlideScore, but reduced the accuracy of ChemScore. In addition, the effect of the conserved water on SOMs prediction was found to be dependent on the receptor model and the substrate. We further found that one of snapshots from molecular dynamics simulations on the apo form can improve the prediction accuracy when compared to the crystal structural model.

Structural, enzymatic and pharmacological profiles of AplTX-II - A basic sPLA2 (D49) isolated from the Agkistrodon piscivorus leucostoma snake venom

A basic sPLA2 (D49) from the venom of snake Agkistrodon piscivorus leucostoma (AplTX-II) was isolated, purified and characterized. We determined the enzymatic and pharmacological profiles of this toxin. AplTX-II was isolated with a high level of purity through reverse phase chromatography and molecular exclusion. The enzyme showed pI 9.48 and molecular weight of 14,003 Da. The enzymatic activity of the AplTX-II depended on Ca2+ pH and temperature. The comparison of the primary structure with other sPLA2s revealed that AplTX-II presented all the structural reasons expected for a basic sPLA2s. Additionally, we have resolved its structure with the docked synthetic substrate NOBA (4-nitro-3-octanoyloxy benzoic acid) by homology modeling, and performed MD simulations with explicit solvent. Structural similarities were found between the enzyme's modeled structure and other snake sPLA2 X-Ray structures, available in the PDB database. NOBA and active-site water molecules spontaneously adopted stable positions and established interactions in full agreement with the reaction mechanism, proposed for the physiological substrate, suggesting that NOBA hydrolysis is an excellent model to study phospholipid hydrolysis.