Molecular modeling insights into the catalytic mechanism of the retaining galactosyltransferase LgtC

The bacterial enzyme lipopolysaccharyl alpha-galactosyltransferase C (EC 2.4.1.x, LgtC) is involved in the synthesis of lipooligosaccharides displayed on the cell surfaces of Neisseria meningitidis. LgtC catalyzes the transfer of a galactosyl residue from UDP-Gal to the terminal galactose residue of glycoconjugates with an overall retention of stereochemistry at the anomeric center. Several hypothetical catalytic mechanisms of the LgtC enzyme were examined herein using DFT quantum chemical methods up to the B3LYP/6-311++G**//B3LYP/6-31G* level. The computational model used to follow the reaction is based on the crystallographic structure of LgtC in complex with both the nucleotide-galactose donor and the oligosaccharide-acceptor analogues. The 136 atoms included in this model represent fragments of residues critical for the substrate binding and catalysis. From our calculations, the preferred pathway is predicted to be a one step mechanism with the nucleophilic attack of the acceptor oxygen onto the anomeric carbon and the proton transfer to a phosphate oxygen occurring simultaneously. This mechanism has an A(N)D(N)A(H)D(H) character, with the unique transition state structure in which the attacking galactose group is more closely bound to the anomeric carbon than to the UDP leaving group and where the hydrogen bond between the nucleophile and the leaving group oxygens facilitates the attack of the acceptor O4(') from the same side of the transferred galactose.     

Molecular modeling and its experimental verification for the catalytic mechanism of Candida antarctica lipase B

Quantum mechanical and molecular dynamics simulation analysis has been performed on the model system for CALB (Candida antarctica lipase B) with esters to study the reaction mechanism and conformational preference of catalytic hydrolysis and the esterification reaction. Using quantum mechanical analysis, the ping-pong bi-bi mechanism was applied and energies and 3-dimensional binding configurations of the whole reaction pathways were calculated. Further molecular dynamics simulation analysis was performed on the basis of the transition state obtained from quantum mechanical study to observe the effect of structures of,the substrates. Calculation results using substrates of different chain length and chiral configurations were compared for conformational preference. The calculated results showed very small influence on chain length, whereas chiral conformation showed big differences. Calculated results from molecular modeling studies have been compared qualitatively with the experimental data using racemic mixtures of (+/-)-cis-4-acetamido-cyclopent-2-ene-1-ethyl acetate as substrates.     

Exploration of the conformational landscape in pregnane X receptor reveals a new binding pocket

Ligand‐regulated pregnane X receptor (PXR), a member of the nuclear receptor superfamily, plays a central role in xenobiotic metabolism. Despite its critical role in drug metabolism, PXR activation can lead to adverse drug‐drug interactions and early stage metabolism of drugs. Activated PXR can induce cancer drug resistance and enhance the onset of malignancy. Since promiscuity in ligand binding makes it difficult to develop competitive inhibitors targeting PXR ligand binding pocket (LBP), it is essential to identify allosteric sites for effective PXR antagonism. Here, molecular dynamics (MD) simulation studies unravelled the existence of two different conformational states, namely “expanded” and “contracted”, in apo PXR ligand binding domain (LBD). Ligand binding events shifted this conformational equilibrium and locked the LBD in a single “ligand‐adaptable” conformational state. Ensemble‐based computational solvent mapping identified a transiently open potential small molecule binding pocket between α5 and α8 helices, named “α8 pocket”, whose opening‐closing mechanism directly correlated with the conformational shift in LBD. A virtual hit identified through structure‐based virtual screening against α8 pocket locks the pocket in its open conformation. MD simulations further revealed that the presence of small molecule at allosteric site disrupts the LBD dynamics and locks the LBD in a “tightly‐contracted” conformation. The molecular details provided here could guide new structural studies to understand PXR activation and antagonism.     

Role of protein flexibility in the catalytic cycle of p-hydroxybenzoate hydroxylase elucidated by the Pro293Ser mutant

Proline 293 of p-hydroxybenzoate hydroxylase from Pseudomonas aeruginosa is in a highly conserved region of the flavoprotein aromatic hydroxylases. It is thought to impart rigidity to the backbone, as it partially cradles the FAD in these hydroxylases. Thus, this residue has been substituted with serine by site-directed mutagenesis to investigate the importance of flexibility of the peptide segment in catalysis. Differential scanning calorimetry demonstrated that the mutation has decreased the stability of the folded mutant protein compared to the wild-type PHBH. The increased flexibility in the protein backbone enhanced the accessibility of the flavin hydroperoxide intermediate to the solvent, causing an increase in the elimination of H(2)O(2) from this labile intermediate and, consequently, a decrease in the efficiency of substrate hydroxylation. Additionally, the increased accessibility of this mutant form of the enzyme makes it more susceptible than the wild-type enzyme to being trapped in the hydroxyflavin intermediate form in the presence of high levels of p-hydroxybenzoate. The mutation also lowers the pK(a) of the phenolic oxygen of bound p-hydroxybenzoate, and eliminates the pH dependence of the rate constant for flavin reduction by NADPH. These experimental observations lead to a model that explains how the wild-type protein can sense the charge of the 4-substituent of the aromatic ligand and link this charge to a flavin conformational change that is required for reaction with NADPH: (i) The peptide oxygen of Pro 293 is repelled by the negative charge of the phenolic oxygen of p-hydroxybenzoate. (ii) This repulsion is transmitted through the peptide backbone, causing the movement of Asn 300. (iii) The change in the position of Asn 300 triggers the movement of the flavin from the largely buried "in" conformation to the exposed, reactive "out" conformation.     

Molecular cloning of a ferret angiotensin II AT(1) receptor reveals the importance of position 163 for Losartan binding

A complementary DNA for the angiotensin II (AngII) type 1 (AT(1)) receptor from Mustela putorius furo (ferret) was isolated from a ferret atria cDNA library. The cDNA encodes a protein (fAT(1)) of 359 amino acids having high homologies (93-99%) to other mammalian AT(1) receptor counterparts. When fAT(1) was expressed in COS-7 cells and photoaffinity labeled with the photoactive analogue (125)I-?Sar(1), Bpa(8)AngII, a protein of 100 kDa was detected by autoradiography. The formation of this complex was specific since it was abolished in the presence of the AT(1) non-peptidic antagonist L-158,809. Functional analysis indicated that the fAT(1) receptor efficiently coupled to phospholipase C as demonstrated by an increase in inositol phosphate production following stimulation with AngII. Binding studies revealed that the fAT(1) receptor had a high affinity for the peptide antagonist ?Sar(1), Ile(8)AngII (K(d) of 5. 8+/-1.4 nM) but a low affinity for the AT(1) selective non-peptidic antagonist DuP 753 (K(d) of 91+/-15.6 nM). Interestingly, when we substituted Thr(163) with an Ala residue, which occupies this position in many mammalian AT(1) receptors, we restored the high affinity of this receptor for Dup 753 (11.7+/-5.13 nM). These results suggest that position 163 of the AT(1) receptor does not contribute to the overall binding of peptidic ligands but that certain non-peptidic antagonists such as Dup 753 are clearly dependent on this position for efficient binding.     

Molecular modeling of Mycobacterium tuberculosis dUTpase: docking and catalytic mechanism studies

Mycobacterium tuberculosis is a leading cause of infectious disease in the world today. This outlook is aggravated by a growing number of M. tuberculosis infections in individuals who are immunocompromised as a result of HIV infections. Thus, new and more potent anti-TB agents are necessary. Therefore, dUTpase was selected as a target enzyme to combat M. tuberculosis. In this work, molecular modeling methods involving docking and QM/MM calculations were carried out to investigate the binding orientation and predict binding affinities of some potential dUTpase inhibitors. Our results suggest that the best potential inhibitor investigated, among the compounds studied in this work, is the compound dUPNPP. Regarding the reaction mechanism, we concluded that the decisive stage for the reaction is the stage 1. Furthermore, it was also observed that the compounds with a -1 electrostatic charge presented lower activation energy in relation to the compounds with a -2 charge.     

Systems biology modeling reveals a possible mechanism of the tumor cell death upon oncogene inactivation in EGFR addicted cancers

Despite many evidences supporting the concept of "oncogene addiction" and many hypotheses rationalizing it, there is still a lack of detailed understanding to the precise molecular mechanism underlying oncogene addiction. In this account, we developed a mathematic model of epidermal growth factor receptor (EGFR) associated signaling network, which involves EGFR-driving proliferation/pro-survival signaling pathways Ras/extracellular-signal-regulated kinase (ERK) and phosphoinositol-3 kinase (PI3K)/AKT, and pro-apoptotic signaling pathway apoptosis signal-regulating kinase 1 (ASK1)/p38. In the setting of sustained EGFR activation, the simulation results show a persistent high level of proliferation/pro-survival effectors phospho-ERK and phospho-AKT, and a basal level of pro-apoptotic effector phospho-p38. The potential of p38 activation (apoptotic potential) due to the elevated level of reactive oxygen species (ROS) is largely suppressed by the negative crosstalk between PI3K/AKT and ASK1/p38 pathways. Upon acute EGFR inactivation, the survival signals decay rapidly, followed by a fast increase of the apoptotic signal due to the release of apoptotic potential. Overall, our systems biology modeling together with experimental validations reveals that inhibition of survival signals and concomitant release of apoptotic potential jointly contribute to the tumor cell death following the inhibition of addicted oncogene in EGFR addicted cancers.     

Model assembly study of the ligand binding by p-hydroxybenzoate hydroxylase: Correlation between the calculated binding energies and the experimental dissociation constants

The energies of binding of seven ligands by p-hydroxybenzoate hydroxylase (PHBH) were calculated theoretically. Direct enzyme–ligand interaction energies were calculated using the ab initio quantum mechanical model assembly of the active site at the 3-21G level. Solvation energies of the ligands needed in the evaluation of the binding energies were calculated with the semiempirical AM1–SM2 method and the long-range electrostatic interaction energies between the ligands and the protein matrix classically using the static charge distributions of the ligands and the protein. Energies for proton-transfer between the ligands OH or SH substituent at position 4 and the active-site tyrosine within the ab initio model assemblies were calculated and compared to the corresponding pK_as in aqueous solution. Excluding 3,4-dihydroxybenzoate, the natural product of PHBH, a linear relationship between the calculated binding energies and the experimental binding free energies was found with a correlation coefficient of 0.90. Contributions of the direct enzyme–ligand interaction energies, solvation energies and the long-range electrostatic interaction energies to the calculated binding energies were analyzed. The proton-transfer energies of the ligands with substituents ortho to the ionized OH were found to be perturbed less in the model calculations than the energies of their meta isomers as deduced from the corresponding pK_as. © 1995 Wiley-Liss, Inc.     

Modelling flavin and substrate substituent effects on the activation barrier and rate of oxygen transfer by p-hydroxybenzoate hydroxylase

The simulation of enzymatic reactions, using computer models, is becoming a powerful tool in the most fundamental challenge in biochemistry: to relate the catalytic activity of enzymes to their structure. In the present study, various computed parameters were correlated with the natural logarithm of experimental rate constants for the hydroxylation of various substrate derivatives catalysed by wild-type para-hydroxybenzoate hydroxylase (PHBH) as well as for the hydroxylation of the native substrate (p-hydroxybenzoate) by PHBH reconstituted with a series of 8-substituted flavins. The following relative parameters have been calculated and tested: (a) energy barriers from combined quantum mechanical/molecular mechanical (QM/MM) (AM1/CHARMM) reaction pathway calculations, (b) gas-phase reaction enthalpies (AM1) and (c) differences between the HOMO and LUMO energies of the isolated substrate and cofactor molecules (AM1 and B3LYP/6-31+G(d)). The gas-phase approaches yielded good correlations, as long as similarly charged species are involved. The QM/MM approach resulted in a good correlation, even including differently charged species. This indicates that the QM/MM model accounts quite well for the solvation effects of the active site surroundings, which vary for differently charged species. The correlations obtained demonstrate quantitative structure activity relationships for an enzyme-catalysed reaction including, for the first time, substitutions on both substrate and cofactor.     

Molecular dissection of the inhibitor binding pocket of mitotic kinesin Eg5 reveals mutants that confer resistance to antimitotic agents

The mitotic kinesin Eg5 plays an essential role in establishing the bipolar spindle. Recently, several antimitotic inhibitors have been shown to share a common binding region on Eg5. Considering the importance of Eg5 as a potential drug target for cancer chemotherapy it is essential to understand the molecular mechanism, by which these agents block Eg5 activity, and to determine the "key residues" crucial for inhibition. Eleven residues in the inhibitor binding pocket were mutated and the effects were monitored by kinetic analysis and mass spectrometry. Mutants R119A, D130A, P131A, I136A, V210A, Y211A and L214A abolish the inhibitory effect of monastrol. Results for W127A and R221A are less striking, but inhibitor constants are still considerably modified compared to wild-type Eg5. Only one residue, Leu214, was found to be essential for inhibition by STLC. W127A, D130A, V210A lead to increased K(i)(app) values, but binding of STLC is still tight. R119A, P131A, Y211A and R221A convert STLC into a classical rather than a tight-binding inhibitor with increased inhibitor constants. These results demonstrate that monastrol and STLC interact with different amino acids within the same binding region, suggesting that this site is highly flexible to accommodate different types of inhibitors. The drug specificity is due to multiple interactions not only with loop L5, but also with residues located in helices alpha2 and alpha3. These results suggest that tumour cells might develop resistance to Eg5 inhibitors, by expressing Eg5 point mutants that retain the enzyme activity, but prevent inhibition, a feature that is observed for certain tubulin inhibitors.     

Crystal structure and substrate binding modeling of the uroporphyrinogen-III decarboxylase from Nicotiana tabacum. Implications for the catalytic mechanism

The enzymatic catalysis of many biological processes of life is supported by the presence of cofactors and prosthetic groups originating from the common tetrapyrrole precursor uroporphyrinogen-III. Uroporphyrinogen-III decarboxylase catalyzes its conversion into coproporphyrinogen-III, leading in plants to chlorophyll and heme biosynthesis. Here we report the first crystal structure of a plant (Nicotiana tabacum) uroporphyrinogen-III decarboxylase, together with the molecular modeling of substrate binding in tobacco and human enzymes. Its structural comparison with the homologous human protein reveals a similar catalytic cleft with six invariant polar residues, Arg(32), Arg(36), Asp(82), Ser(214) (Thr in Escherichia coli), Tyr(159), and His(329) (tobacco numbering). The functional relationships obtained from the structural and modeling analyses of both enzymes allowed the proposal for a refined catalytic mechanism. Asp(82) and Tyr(159) seem to be the catalytic functional groups, whereas the other residues may serve in substrate recognition and binding, with Arg(32) steering its insertion. The crystallographic dimer appears to represent the protein dimer under physiological conditions. The dimeric arrangement offers a plausible mechanism at least for the first two (out of four) decarboxylation steps.     

Site-directed mutagenesis reveals transition-state stabilization as a general catalytic mechanism for aminoacyl-tRNA synthetases

Some aminoacyl-tRNA synthetases of almost negligible homology do have a small region of similarity around four-residue sequence His-Ile(or Leu or Met)-Gly-His(or Asn), the HIGH sequence. The first histidine in this sequence in the tyrosyl-tRNA synthetase, His-45, has been shown to form part of a binding site for the gamma-phosphate of ATP in the transition state for the reaction as does Thr-40. Residue His-56 in the valyl-tRNA synthetase begins a HIGH sequence, and there is a threonine at position 52, one position closer to the histidine than in the tyrosyl-tRNA synthetase. The mutants Thr----Ala-52 and His----Asn-56 have been made and their complete free energy profiles for the formation of valyl adenylate determined. Difference energy diagrams have been constructed by comparison with the reaction of wild-type enzyme. The difference energy profiles are very similar to those for the mutants Thr----Ala-40 and His----Asn-45 of the tyrosyl-tRNA synthetase. Thr-52 and His-56 of the valyl-tRNA synthetase contribute little binding energy to valine, ATP, and Val-AMP. Instead, the wild-type enzyme binds the transition state and pyrophosphate some 6 kcal/mol more tightly than do the mutants. Preferential transition-state stabilization is thus an important component of catalysis by the valyl-tRNA synthetase. Further, by analogy to the tyrosyl-tRNA synthetase, the valyl-tRNA synthetase has a binding site for the gamma-phosphate of ATP in the transition state, and this is likely to be a general feature of aminoacyl-tRNA synthetases that have a HIGH region.     

Analysis of the active site of the flavoprotein p-hydroxybenzoate hydroxylase and some ideas with respect to its reaction mechanism

The flavoprotein p-hydroxybenzoate hydroxylase has been studied extensively by biochemical techniques by others and in our laboratory by X-ray crystallography. As a result of the latter investigations, well-refined crystal structures are known of the enzyme complexed (i) with its substrate p-hydroxybenzoate and (ii) with its reaction product 3,4-dihydroxybenzoate and (iii) the enzyme with reduced FAD. Knowledge of these structures and the availability of the three-dimensional structure of a model compound for the reactive flavin 4a-hydroperoxide intermediate has allowed a detailed analysis of the reaction with oxygen. In the model of this reaction intermediate, fitted to the active site of p-hydroxybenzoate hydroxylase, all possible positions of the distal oxygen were surveyed by rotating this oxygen about the single bond between the C4a and the proximal oxygen. It was found that the distal oxygen is free to sweep an arc of about 180 degrees in the active site. The flavin 4a-peroxide anion, which is formed after reaction of molecular oxygen with reduced FAD, might accept a proton from an active-site water molecule or from the hydroxyl group of the substrate. The position of the oxygen to be transferred with respect to the substrate appears to be almost ideal for nucleophilic attack of the substrate onto this oxygen. The oxygen is situated above the 3-position of the substrate where the substitution takes place, at an angle of about 60 degrees with the aromatic plane, allowing strong interactions with the pi electrons of the substrate. Polarization of the peroxide oxygen-oxygen bond by the enzyme may enhance the reactivity of flavin 4a-peroxide.     

Molecular mechanism for pterin-mediated inactivation of tyrosine hydroxylase: formation of insoluble aggregates of tyrosine hydroxylase

Tyrosine hydroxylase (TH), an iron-containing enzyme, catalyzes the first and rate-limiting step of catecholamine biosynthesis, and requires tetrahydrobiopterin (BH4) as a cofactor. We found that preincubation of recombinant human TH with BH4 results in the irreversible inactivation of the enzyme at a concentration far less than the Km value toward BH4 in spite of its cofactor role, whereas oxidized biopterin, which has no cofactor activity, does not affect the enzyme activity. We show that TH is inactivated by BH4 in competition with the binding of dopamine. The sequential addition of BH4 to TH results in a gradual decrease in the intensity of the fluorescence and CD spectra without changing their overall profiles. Sedimentation velocity analysis demonstrated an association of TH molecules with each other in the presence of BH4, and studies using gel-permeation chromatography, turbidity measurements, and transmission electron microscopy demonstrated the formation of amorphous aggregates with large molecular weights following the association of the TH proteins. These results suggest that BH4 not only acts as a cofactor, but also accelerates the aggregation of TH. We propose a novel mechanism for regulating the amount of TH protein, and discuss its physiological significance.     

Absorption spectra of radical forms of 2,4-dihydroxybenzoic acid, a substrate for p-hydroxybenzoate hydroxylase.

Combined optical and conductimetric measurements in aqueous solution indicate that at high pH (greater than or equal to 10).OH radicals react with the phenoxide form of 2,4-dihydroxybenzoic acid to form transiently phenoxyl radicals and a small amount of hydroxyeyclohexadienyl (HCHD) radicals by 150 ns. The respective yields of 88 and 12% of the total.OH radical yield were deduced from conductance and optical changes as well as from studies using a low potential reductant. The HCHD radical possesses a pKa of 8.0 +/- 0.1 and the constructed spectrum of the deprotonated forms of HCHD has a lambda max at 420 nm with a minimum extinction coefficient of approximately 7250 M-1 cm-1. The red shift in lambda max and increase in extinction coefficient compared to the revised spectral properties of the protonated form of the HCHD radical (lambda max at 390 nm with extinction coefficient of approximately 4500 M-1 cm-1), together with the pKa of the HCHD radical, provide an explanation for the pH-dependent spectral changes of the so-called highly absorbing intermediate II species, observed in the functioning of the enzyme p-hydroxybenzoate hydroxylase. These results add further to the evidence in support of the proposal that intermediate II is composed of species which absorb similarly to the flavin 4(a)-hydroxide and a form of the substrate/product such as the HCHD radical (Anderson, R. F., Patel, K. B., and Stratford, M. R. L. (1987) J. Biol. Chem. 262, 17475-17479).     

Analysis of the Staphylococcus aureus DgkB structure reveals a common catalytic mechanism for the soluble diacylglycerol kinases

Soluble diacylglycerol (DAG) kinases function as regulators of diacylglycerol metabolism in cell signaling and intermediary metabolism. We report the structure of a DAG kinase, DgkB from Staphylococcus aureus, both as the free enzyme and in complex with ADP. The molecule is a tight homodimer, and each monomer comprises two domains with the catalytic center located within the interdomain cleft. Two distinctive features of DkgB are a structural Mg2+ site and an associated Asp*water*Mg2+ network that extends toward the active site locale. Site-directed mutagenesis revealed that these features play important roles in the catalytic mechanism. The key active site residues and the components of the Asp*water*Mg2+ network are conserved in the catalytic cores of the mammalian signaling DAG kinases, indicating that these enzymes use the same mechanism and have similar structures as DgkB.     

Farnesol-induced apoptosis in Aspergillus nidulans reveals a possible mechanism for antagonistic interactions between fungi

The dimorphic fungus Candida albicans secretes farnesol, which acts as a quorum-sensing molecule and prevents the yeast to mycelium conversion. In this study we examined the effect of farnesol in the filamentous fungus Aspergillus nidulans. We show that externally added farnesol has no effect on hyphal morphogenesis; instead, it triggers morphological features characteristic of apoptosis. Additional experiments suggest that mitochondria and reactive oxygen species (ROS) participate in farnesol-induced apoptosis. Moreover, the effects of farnesol appear to be mediated by the FadA heterotrimeric G protein complex. Because A. nidulans does not secrete detectable amounts of farnesol, we propose that it responds to farnesol produced by other fungi. In agreement with this notion, growth and development were impaired in a farnesol-dependent manner when A. nidulans was co-cultivated with C. albicans. Taken together, our data suggest that farnesol, in addition to its quorum-sensing function that regulates morphogenesis, is also employed by C. albicans to reduce competition from other microbes.     

Kinetic analysis of autotaxin reveals substrate-specific catalytic pathways and a mechanism for lysophosphatidic acid distribution

Autotaxin (ATX) is a secreted lysophospholipase D that hydrolyzes lysophosphatidylcholine (LPC) into lysophosphatidic acid (LPA), initiating signaling cascades leading to cancer metastasis, wound healing, and angiogenesis. Knowledge of the pathway and kinetics of LPA synthesis by ATX is critical for developing quantitative physiological models of LPA signaling. We measured the individual rate constants and pathway of the LPA synthase cycle of ATX using the fluorescent lipid substrates FS-3 and 12-(N-methyl-N-(7-nitrobenz-2-oxa-1,3-diazol-4-yl))-LPC. FS-3 binds rapidly (k(1) ≥500 μm(-1) s(-1)) and is hydrolyzed slowly (k(2) = 0.024 s(-1)). Release of the first hydrolysis product is random and rapid (≥1 s(-1)), whereas release of the second is slow and rate-limiting (0.005-0.007 s(-1)). Substrate binding and hydrolysis are slow and rate-limiting with LPC. Product release is sequential with choline preceding LPA. The catalytic pathway and kinetics depend strongly on the substrate, suggesting that ATX kinetics could vary for the various in vivo substrates. Slow catalysis with LPC reveals the potential for LPA signaling to spread to cells distal to the site of LPC substrate binding by ATX. An ATX mutant in which catalytic threonine at position 210 is replaced with alanine binds substrate weakly, favoring a role for Thr-210 in binding as well as catalysis. FTY720P, the bioactive form of a drug currently used to treat multiple sclerosis, inhibits ATX in an uncompetitive manner and slows the hydrolysis reaction, suggesting that ATX inhibition plays a significant role in lymphocyte immobilization in FTY720P-based therapeutics.