Exploration of the effects of sequence variations between HIV-1 and HIV-2 proteases on their three-dimensional structures

Abstract

HIV protease inhibitors (PIs) approved by the FDA (US Food and Drug Administration) are a major class of antiretroviral. HIV-2 protease (PR2) is naturally resistant to most of them as PIs were designed for HIV-1 protease (PR1). In this study, we explored the impact of amino-acid substitutions between PR1 and PR2 on the structure of protease (PR) by comparing the structural variability of 13 regions using 24 PR1 and PR2 structures complexed with diverse ligands. Our analyses confirmed structural rigidity of the catalytic region and highlighted the important role of three regions in the conservation of the catalytic region conformation. Surprisingly, we showed that the flap region, corresponding to a flexible region, exhibits similar conformations in PR1 and PR2. Furthermore, we identified regions exhibiting different conformations in PR1 and PR2, which could be explained by the intrinsic flexibility of these regions, by crystal packing, or by PR1 and PR2 substitutions. Some substitutions induce structural changes in the R2 and R4 regions that could have an impact on the properties of PI-binding site and could thus modify PI binding mode. Substitutions involved in structural changes in the elbow region could alter the flexibility of the PR2 flap regions relative to PR1, and thus play a role in the transition from the semi-open form to the closed form, and have an impact on ligand binding. These results improve the understanding of the impact of sequence variations between PR1 and PR2 on the natural resistance of HIV-2 to commercially available PIs.

Communicated by Ramaswamy H. Sarma     

Structural snapshots of the KMSKS loop rearrangement for amino acid activation by bacterial tyrosyl-tRNA synthetase

Tyrosyl-tRNA synthetase (TyrRS) has been studied extensively by mutational and structural analyses to elucidate its catalytic mechanism. TyrRS has the HIGH and KMSKS motifs that catalyze the amino acid activation with ATP. In the present study, the crystal structures of the Escherichia coli TyrRS catalytic domain, in complexes with l-tyrosine and a l-tyrosyladenylate analogue, Tyr-AMS, were solved at 2.0A and 2.7A resolution, respectively. In the Tyr-AMS-bound structure, the 2'-OH group and adenine ring of the Tyr-AMS are strictly recognized by hydrogen bonds. This manner of hydrogen-bond recognition is conserved among the class I synthetases. Moreover, a comparison between the two structures revealed that the KMSKS loop is rearranged in response to adenine moiety binding and hydrogen-bond formation, and the KMSKS loop adopts the more compact ("semi-open") form, rather than the flexible, open form. The HIGH motif initially recognizes the gamma-phosphate, and then the alpha and gamma-phosphates of ATP, with a slight rearrangement of the residues. The other residues around the substrate also accommodate the Tyr-AMS. This induced-fit form presents a novel "snapshot" of the amino acid activation step in the aminoacylation reaction by TyrRS. The present structures and the T.thermophilus TyrRS ATP-free and bound structures revealed that the extensive induced-fit conformational changes of the KMSKS loop and the local conformational changes within the substrate binding site form the basis for driving the amino acid activation step: the KMSKS loop adopts the open form, transiently shifts to the semi-open conformation according to the adenosyl moiety binding, and finally assumes the rigid ATP-bound, closed form. After the amino acid activation, the KMSKS loop adopts the semi-open form again to accept the CCA end of tRNA for the aminoacyl transfer reaction.     

The role of tyrosine 71 in modulating the flap conformations of BACE1

β-Site amyloid precursor protein cleaving enzyme 1 (BACE1) is a potential target for treating Alzheimer's disease. BACE1's binding site is partially covered by a flexible loop on its N-terminal domain, known as the "flap," which has been found in several conformations in crystal structures of BACE1 and other aspartyl proteases. The side chain of the invariant residue Tyr71 on the flap adopts several rotameric orientations, leading to our hypothesis that the orientation of this residue dictates the movement and conformations available to the flap. We investigated this hypothesis by performing 220 ns of molecular dynamics simulations of bound and unbound wild-type BACE1 as well as the unbound Y71A mutant. Our findings indicate that the flap exhibits various degrees of mobility and adopts different conformations depending on the Tyr71 orientation. Surprisingly, the "self-inhibited" form is stable in our simulations, making it a reasonable target for drug design. The alanine mutant, lacking a large side chain at position 71, displays significant differences in flap dynamics from wild type, freely sampling very open and closed conformations. Our simulations show that Tyr71, in addition to its previously determined functions in catalysis and substrate binding, has the important role of modulating flap conformations in BACE1.     

HIV-1 protease substrate binding and product release pathways explored with coarse-grained molecular dynamics

We analyze the encounter of a peptide substrate with the native HIV-1 protease, the mechanism of substrate incorporation in the binding cleft, and the dissociation of products after substrate hydrolysis. To account for the substrate, we extend a coarse-grained model force field, which we previously developed to study the flap opening dynamics of HIV-1 protease on a microsecond timescale. Molecular and Langevin dynamics simulations show that the flaps need to open for the peptide to bind and that the protease interaction with the substrate influences the flap opening frequency and interval. On the other hand, release of the products does not require flap opening because they can slide out from the binding cleft to the sides of the enzyme. Our data show that in the protease-substrate complex the highest fluctuations correspond to the 17- and 39-turns and the substrate motion is anticorrelated with the 39-turn. Moreover, the active site residues and the flap tips move in phase with the peptide. We suggest some mechanistic principles for how the flexibility of the protein may be involved in ligand binding and release.     

Dynamic flaps in HIV-1 protease adopt unique ordering at different stages in the catalytic cycle

The flexibility of HIV-1 protease flaps is known to be essential for the enzymatic activity. Here we attempt to capture a multitude of conformations of the free and substrate-bound HIV-1 protease that differ drastically in their flap arrangements. The substrate binding process suggests the opening of active site gate in conjunction with a reversal of flap tip ordering, from the native semiopen state. The reversed-flap, open-gated enzyme readily transforms to a closed conformation after proper placement of the substrate into the binding cleft. After substrate processing, the closed state protease which possessed opposite flap ordering relative to the semiopen state, encounters another flap reversal via a second open conformation that facilitates the evolution of native semiopen state of correct flap ordering. The complicated transitional pathway, comprising of many high and low energy states, is explored by combining standard and activated molecular dynamics (MD) simulation techniques. The study not only complements the existing findings from X-ray, NMR, EPR, and MD studies but also provides a wealth of detailed information that could help the structure-based drug design process.     

A molecular dynamics study comparing a wild-type with a multiple drug resistant HIV protease: differences in flap and aspartate 25 cavity dimensions

HIV proteases can develop resistance to therapeutic drugs by mutating specific residues, but still maintain activity with their natural substrates. To gain insight into why mutations confer such resistance, long ( approximately 70 ns) Molecular Dynamics simulations in explicit solvent were performed on a multiple drug resistant (MDR) mutant (with Asn25 in the crystal structure mutated in silico back to the catalytically active Asp25) and a wild type (WT) protease. HIV proteases are homodimers, with characteristic flap tips whose conformations and dynamics are known to be important influences of ligand binding to the aspartates that form the catalytic center. The WT protease undergoes a transition between 25 and 35 ns that is absent in the MDR protease. The origin of this distinction is investigated using principal component analysis, and is related to differences in motion mainly in the flap region of each monomer. Trajectory analysis suggests that the WT transition arises from a concerted motion of the flap tip distances to their catalytic aspartate residues, and the distance between the two flap tips. These distances form a triangle that in the WT expands the active site from an initial (semi-open) form to an open form, in a correlated manner. In contrast, the MDR protease remains in a more closed configuration, with uncorrelated fluctuations in the distances defining the triangle. This contrasting behavior suggests that the MDR mutant achieves its resistance to drugs by making its active site less accessible to inhibitors. The migration of water to the active site aspartates is monitored. Water molecules move in and out of the active site and individual waters hydrogen bond to both aspartate carboxylate oxygens, with residence times in the ns time regime.     

Probing conformational plasticity of the activation domain of trypsin: the role of glycine hinges

Trypsin-like serine proteases play essential roles in diverse physiological processes such as hemostasis, apoptosis, signal transduction, reproduction, immune response, matrix remodeling, development, and differentiation. All of these proteases share an intriguing activation mechanism that involves the transition of an unfolded domain (activation domain) of the zymogen to a folded one in the active enzyme. During this conformational change, activation domain segments move around highly conserved glycine hinges. In the present study, hinge glycines were replaced by alanine residues via site directed mutagenesis. The effects of these mutations on the interconversion of the zymogen-like and active conformations as well as on catalytic activity were studied. Mutant trypsins showed zymogen-like structures to varying extents characterized by increased flexibility of some activation domain segments, a more accessible N-terminus and a deformed substrate binding site. Our results suggest that the trypsinogen to trypsin transition is hindered by the mutations, which results in a shift of the equilibrium between the inactive zymogen-like and active enzyme conformations toward the inactive state. Our data also showed, however, that the inactive conformations of the various mutants differ from each other. Binding of substrate analogues shifted the conformational equilibrium toward the active enzyme since inhibited forms of the trypsin mutants showed similar structural features as the wild-type enzyme. The catalytic activity of the mutants correlated with the proper conformation of the active site, which could be supported by varying conformations of the N-terminus and the autolysis loop. Transient kinetic measurements confirmed the existence of an inactive to active conformational transition occurring prior to substrate binding.     

Open and closed states of Mrlip1 DAG lipase revealed by molecular dynamics simulation

ABSTRACT

The lid and flap domains control the catalytic activity of lipase through the opening and closing motion. However, this gating mechanism of diacylglycerol (DAG) lipase is poorly understood due to the lack of 3D structures in open conformations. In this study, the opening and closing states of Mrlip1 DAG lipase are revealed by the homology modelling and molecular dynamic simulations. It was found that the active residues (Ser₁₇₁, His₂₈₁ and Asp₂₂₈) in the catalytic pocket of Mrlip1 DAG lipase are covered by the lid domain in the closed conformation, and exposed to the solvent in the open conformation. The role of residues Phe₂₇₈ and Gln₂₈₂ in the flap domain, as well as that of Thr₁₀₁ and Thr₁₀₇ in the lid domains are also identified in gating mechanism. The site-directed mutagenesis have been carried out to illustrate the putative alterations of enzyme specificity. Our results suggest that the substrate specificity is achieved by these two key residues Phe₂₇₈ and Gln₂₈₂, and the irreversible conversion from DAG to TAG (Triacylglycerol) lipase are enabled by the two-point mutations.     

Molecular docking based virtual screening of natural compounds as potential BACE1 inhibitors: 3D QSAR pharmacophore mapping and molecular dynamics analysis

Beta-site APP cleaving enzyme1 (BACE1) catalyzes the rate determining step in the generation of Aβ peptide and is widely considered as a potential therapeutic drug target for Alzheimer’s disease (AD). Active site of BACE1 contains catalytic aspartic (Asp) dyad and flap. Asp dyad cleaves the substrate amyloid precursor protein with the help of flap. Currently, there are no marketed drugs available against BACE1 and existing inhibitors are mostly pseudopeptide or synthetic derivatives. There is a need to search for a potent inhibitor with natural scaffold interacting with flap and Asp dyad. This study screens the natural database InterBioScreen, followed by three-dimensional (3D) QSAR pharmacophore modeling, mapping, in silico ADME/T predictions to find the potential BACE1 inhibitors. Further, molecular dynamics of selected inhibitors were performed to observe the dynamic structure of protein after ligand binding. All conformations and the residues of binding region were stable but the flap adopted a closed conformation after binding with the ligand. Bond oligosaccharide interacted with the flap as well as catalytic dyad via hydrogen bond throughout the simulation. This led to stabilize the flap in closed conformation and restricted the entry of substrate. Carbohydrates have been earlier used in the treatment of AD because of their low toxicity, high efficiency, good biocompatibility, and easy permeability through the blood–brain barrier. Our finding will be helpful in identify the potential leads to design novel BACE1 inhibitors for AD therapy.     

Molecular Dynamics Study of the Opening Mechanism for DNA Polymerase I

During DNA replication, DNA polymerases follow an induced fit mechanism in order to rapidly distinguish between correct and incorrect dNTP substrates. The dynamics of this process are crucial to the overall effectiveness of catalysis. Although X-ray crystal structures of DNA polymerase I with substrate dNTPs have revealed key structural states along the catalytic pathway, solution fluorescence studies indicate that those key states are populated in the absence of substrate. Herein, we report the first atomistic simulations showing the conformational changes between the closed, open, and ajar conformations of DNA polymerase I in the binary (enzyme∶DNA) state to better understand its dynamics. We have applied long time-scale, unbiased molecular dynamics to investigate the opening process of the fingers domain in the absence of substrate for B. stearothermophilis DNA polymerase in silico. These simulations are biologically and/or physiologically relevant as they shed light on the transitions between states in this important enzyme. All closed and ajar simulations successfully transitioned into the fully open conformation, which is known to be the dominant binary enzyme-DNA conformation from solution and crystallographic studies. Furthermore, we have detailed the key stages in the opening process starting from the open and ajar crystal structures, including the observation of a previously unknown key intermediate structure. Four backbone dihedrals were identified as important during the opening process, and their movements provide insight into the recognition of dNTP substrate molecules by the polymerase binary state. In addition to revealing the opening mechanism, this study also demonstrates our ability to study biological events of DNA polymerase using current computational methods without biasing the dynamics.     

Molecular Dynamics Simulation Studies of dTTP Binding and Catalysis Mediated by YhdE Dimerization

YhdE is a Maf-like (multicopy associated filamentation) protein that primarily acts as dTTPase to hydrolyze dTTP into dTMP and two phosphate molecules in cell metabolism pathway. Two crystal structures of YhdE have been previously determined, representing the open and closed active site conformations, respectively. Based on the structures, we have carried out molecular dynamics simulations and free energy calculations to investigate dTTP binding to and hydrolysis by YhdE. Our results suggest that YhdE closed state is structurally more compact than its open state at room temperature. YhdE open state is a favorable conformation for dTTP binding and closed state is a structurally favorable conformation for catalytic reaction. This observation is supported by the structure of YhdE homolog in complex with a nucleotide analog. Free energy calculations reveal that YhdE dimerization occurs preferentially in dTTP binding and is favorable for successive cooperative reaction. The key residues R11, R12 and K80, are found to contribute to the substrate stabilization. Further, YhdE dimerization and binding of dTTP induce the cooperative effect through a direct allosteric communication network in YhdE from the dTTP binding sites in the catalytic center to the intermolecular β-strand in YhdE dimer.     

Charge neutralization in the active site of the catalytic trimer of aspartate transcarbamoylase promotes diverse structural changes

A classical model for allosteric regulation of enzyme activity posits an equilibrium between inactive and active conformations. An alternative view is that allosteric activation is achieved by increasing the potential for conformational changes that are essential for catalysis. In the present study, substitution of a basic residue in the active site of the catalytic (C) trimer of aspartate transcarbamoylase with a non‐polar residue results in large interdomain hinge changes in the three chains of the trimer. One conformation is more open than the chains in both the wild‐type C trimer and the catalytic chains in the holoenzyme, the second is closed similar to the bisubstrate‐analog bound conformation and the third hinge angle is intermediate to the other two. The active‐site 240s loop conformation is very different between the most open and closed chains, and is disordered in the third chain, as in the holoenzyme. We hypothesize that binding of anionic substrates may promote similar structural changes. Further, the ability of the three catalytic chains in the trimer to access the open and closed active‐site conformations simultaneously suggests a cyclic catalytic mechanism, in which at least one of the chains is in an open conformation suitable for substrate binding whereas another chain is closed for catalytic turnover. Based on the many conformations observed for the chains in the isolated catalytic trimer to date, we propose that allosteric activation of the holoenzyme occurs by release of quaternary constraint into an ensemble of active‐site conformations.     

Exploring the conformational states and rearrangements of Yarrowia lipolytica Lipase

We report the 1.7 Å resolution crystal structure of the Lip2 lipase from Yarrowia lipolytica in its closed conformation. The Lip2 structure is highly homologous to known structures of the fungal lipase family (Thermomyces lanuginosa, Rhizopus niveus, and Rhizomucor miehei lipases). However, it also presents some unique features that are described and discussed here in detail. Structural differences, in particular in the conformation adopted by the so-called lid subdomain, suggest that the opening mechanism of Lip2 may differ from that of other fungal lipases. Because the catalytic activity of lipases is strongly dependent on structural rearrangement of this mobile subdomain, we focused on elucidating the molecular mechanism of lid motion. Using the x-ray structure of Lip2, we carried out extensive molecular-dynamics simulations in explicit solvent environments (water and water/octane interface) to characterize the major structural rearrangements that the lid undergoes under the influence of solvent or upon substrate binding. Overall, our results suggest a two-step opening mechanism that gives rise first to a semi-open conformation upon adsorption of the protein at the water/organic solvent interface, followed by a further opening of the lid upon substrate binding.     

A kinetic alignment of orthologous inosine-5'-monophosphate dehydrogenases

IMP dehydrogenase (IMPDH) catalyzes two very different chemical transformations, a dehydrogenase reaction and a hydrolysis reaction. The enzyme toggles between the open conformation required for the dehydrogenase reaction and the closed conformation of the hydrolase reaction by moving a mobile flap into the NAD site. Despite these multiple functional constraints, the residues of the flap and NAD site are highly diverged, and the equilibrium between open and closed conformations ( K c ) varies widely. In order to understand how differences in the dynamic properties of the flap influence the catalytic cycle, we have delineated the kinetic mechanism of IMPDH from the pathogenic protozoan parasite Cryptosporidium parvum ( CpIMPDH), which was obtained from a bacterial source through horizontal gene transfer, and its host counterpart, human IMPDH type 2 (hIMPDH2). Interestingly, the intrinsic binding energy of NAD (+) differentially distributes across the dinucleotide binding sites of these two enzymes as well as in the previously characterized IMPDH from Tritrichomonas foetus ( TfIMPDH). Both the dehydrogenase and hydrolase reactions display significant differences in the host and parasite enzymes, in keeping with the phylogenetic and structural divergence of their active sites. Despite large differences in K c , the catalytic power of both the dehydrogenase and hydrolase conformations are similar in CpIMPDH and TfIMPDH. This observation suggests that the closure of the flap simply sets the stage for catalysis rather than plays a more active role in the chemical transformation. This work provides the essential mechanistic framework for drug discovery.     

Structural, kinetic, and thermodynamic studies of specificity designed HIV-1 protease

HIV-1 protease recognizes and cleaves more than 12 different substrates leading to viral maturation. While these substrates share no conserved motif, they are specifically selected for and cleaved by protease during viral life cycle. Drug resistant mutations evolve within the protease that compromise inhibitor binding but allow the continued recognition of all these substrates. While the substrate envelope defines a general shape for substrate recognition, successfully predicting the determinants of substrate binding specificity would provide additional insights into the mechanism of altered molecular recognition in resistant proteases. We designed a variant of HIV protease with altered specificity using positive computational design methods and validated the design using X-ray crystallography and enzyme biochemistry. The engineered variant, Pr3 (A28S/D30F/G48R), was designed to preferentially bind to one out of three of HIV protease's natural substrates; RT-RH over p2-NC and CA-p2. In kinetic assays, RT-RH binding specificity for Pr3 increased threefold compared to the wild-type (WT), which was further confirmed by isothermal titration calorimetry. Crystal structures of WT protease and the designed variant in complex with RT-RH, CA-p2, and p2-NC were determined. Structural analysis of the designed complexes revealed that one of the engineered substitutions (G48R) potentially stabilized heterogeneous flap conformations, thereby facilitating alternate modes of substrate binding. Our results demonstrate that while substrate specificity could be engineered in HIV protease, the structural pliability of protease restricted the propagation of interactions as predicted. These results offer new insights into the plasticity and structural determinants of substrate binding specificity of the HIV-1 protease.     

Crystal structures of mouse class II alcohol dehydrogenase reveal determinants of substrate specificity and catalytic efficiency

The structure of mouse class II alcohol dehydrogenase (ADH2) has been determined in a binary complex with the coenzyme NADH and in a ternary complex with both NADH and the inhibitor N-cyclohexylformamide to 2.2 A and 2.1 A resolution, respectively. The ADH2 dimer is asymmetric in the crystal with different orientations of the catalytic domains relative to the coenzyme-binding domains in the two subunits, resulting in a slightly different closure of the active-site cleft. Both conformations are about half way between the open apo structure and the closed holo structure of horse ADH1, thus resembling that of ADH3. The semi-open conformation and structural differences around the active-site cleft contribute to a substantially different substrate-binding pocket architecture as compared to other classes of alcohol dehydrogenase, and provide the structural basis for recognition and selectivity of alcohols and quinones. The active-site cleft is more voluminous than that of ADH1 but not as open and funnel-shaped as that of ADH3. The loop with residues 296-301 from the coenzyme-binding domain is short, thus opening up the pocket towards the coenzyme. On the opposite side, the loop with residues 114-121 stretches out over the inter-domain cleft. A cavity is formed below this loop and adds an appendix to the substrate-binding pocket. Asp301 is positioned at the entrance of the pocket and may control the binding of omega-hydroxy fatty acids, which act as inhibitors rather than substrates. Mouse ADH2 is known as an inefficient ADH with a slow hydrogen-transfer step. By replacing Pro47 with His, the alcohol dehydrogenase activity is restored. Here, the structure of this P47H mutant was determined in complex with NADH to 2.5 A resolution. His47 is suitably positioned to act as a catalytic base in the deprotonation of the substrate. Moreover, in the more closed subunit, the coenzyme is allowed a position closer to the catalytic zinc. This is consistent with hydrogen transfer from an alcoholate intermediate where the Pro/His replacement focuses on the function of the enzyme.     

Mechanism of substrate recognition by drug-resistant human immunodeficiency virus type 1 protease variants revealed by a novel structural intermediate

Human immunodeficiency virus type 1 (HIV-1) protease processes and cleaves the Gag and Gag-Pol polyproteins, allowing viral maturation, and therefore is an important target for antiviral therapy. Ligand binding occurs when the flaps open, allowing access to the active site. This flexibility in flap geometry makes trapping and crystallizing structural intermediates in substrate binding challenging. In this study, we report two crystal structures of two HIV-1 protease variants bound with their corresponding nucleocapsid-p1 variant. One of the flaps in each of these structures exhibits an unusual "intermediate" conformation. Analysis of the flap-intermediate and flap-closed crystal structures reveals that the intermonomer flap movements may be asynchronous and that the flap which wraps over the P3 to P1 (P3-P1) residues of the substrate might close first. This is consistent with our hypothesis that the P3-P1 region is crucial for substrate recognition. The intermediate conformation is conserved in both the wild-type and drug-resistant variants. The structural differences between the variants are evident only when the flaps are closed. Thus, a plausible structural model for the adaptability of HIV-1 protease to recognize substrates in the presence of drug-resistant mutations has been proposed.     

Investigation of flap flexibility of β-secretase using molecular dynamic simulations

Flap motif and its dynamics were extensively reported in aspartate proteases, e.g. HIV proteases and plasmepsins. Herein, we report the first account of flap dynamics amongst different conformations of β-secretase using molecular dynamics simulation. Various parameters were proposed and a selected few were picked which could appropriately describe the flap motion. Three systems were studied, namely Free (BACE_Free) and two ligand-bound conformations, which belonged to space groups P6₁22 (BACE_Bound1) and C222₁ (BACE_Bound2), respectively and four parameters (distance between the flaps tip residue, Thr72 and Ser325, d₁; dihedral angle, ? (Thr72-Asp32-Asp228-Ser325); TriCα angles, θ₁ (Thr72-Asp32-Ser325), and θ₂ (Thr72-Asp228-Ser325)) were proposed to understand the change in dynamics of flap domain and the extent of flap opening and closing. Analysis of, θ₂, d₁, θ₁ and ? confirmed that the BACE_Free adopted semi-open, open and closed conformations with slight twisting during flap opening. However, BACE_Bound1 (P6122) showed an adaptation to open conformation due to lack of hydrogen bond interaction between the ligand and flap tip residue. A slight flap twisting, ? (lateral twisting) was observed for BACE_Bound1 during flap opening which correlates with the opening of BACE_Free. Contradictory to the BACE_Bound1, the BACE_Bound2 locked the flap in a closed conformation throughout the simulation due to formation of a stable hydrogen bond interaction between the flap tip residue and ligand. Analyses of all three systems highlight that d₁, θ₂ and ? can be precisely used to describe the extent of flap opening and closing concurrently with snapshots along the molecular dynamics trajectory across several conformations of β-secretase.     

Crystal structures and biochemical studies of human lysophosphatidic acid phosphatase type 6

Lysophosphatidic acid (LPA) is an important bioactive phospholipid involved in cell signaling through Gprotein- coupled receptors pathways. It is also involved in balancing the lipid composition inside the cell, and modulates the function of lipid rafts as an intermediate in phospholipid metabolism. Because of its involvement in these important processes, LPA degradation needs to be regulated as precisely as its production. Lysophosphatidic acid phosphatase type 6 (ACP6) is an LPA-specific acid phosphatase that hydrolyzes LPA to monoacylglycerol (MAG) and phosphate. Here, we report three crystal structures of human ACP6 in complex with malonate, L- (+)-tartrate and tris, respectively. Our analyses revealed that ACP6 possesses a highly conserved Rossmann-foldlike body domain as well as a less conserved cap domain. The vast hydrophobic substrate-binding pocket, which is located between those two domains, is suitable for accommodating LPA, and its shape is different from that of other histidine acid phosphatases, a fact that is consistent with the observed difference in substrate preferences. Our analysis of the binding of three molecules in the active site reveals the involvement of six conserved and crucial residues in binding of the LPA phosphate group and its catalysis. The structure also indicates a water-supplying channel for substrate hydrolysis. Our structural data are consistent with the fact that the enzyme is active as a monomer. In combination with additional mutagenesis and enzyme activity studies, our structural data provide important insights into substrate recognition and the mechanism for catalytic activity of ACP6.