Dopamine transporter comparative molecular modeling and binding site prediction using the LeuT(Aa) leucine transporter as a template

Pharmacological and behavioral studies indicate that binding of cocaine and the amphetamines by the dopamine transporter (DAT) protein is principally responsible for initiating the euphoria and addiction associated with these drugs. The lack of an X-ray crystal structure for the DAT or any other member of the neurotransmitter:sodium symporter (NSS) family has hindered understanding of psychostimulant recognition at the atomic level; structural information has been obtained largely from mutagenesis and biophysical studies. The recent publication of a crystal structure for the bacterial leucine transporter LeuT(Aa), a distantly related NSS family homolog, provides for the first time a template for three-dimensional comparative modeling of NSS proteins. A novel computational modeling approach using the capabilities of the Molecular Operating Environment program MOE 2005.06 in conjunction with other comparative modeling servers generated the LeuT(Aa)-directed DAT model. Probable dopamine and amphetamine binding sites were identified within the DAT model using multiple docking approaches. Binding sites for the substrate ligands (dopamine and amphetamine) overlapped substantially with the analogous region of the LeuT(Aa) crystal structure for the substrate leucine. The docking predictions implicated DAT side chains known to be critical for high affinity ligand binding and suggest novel mutagenesis targets in elucidating discrete substrate and inhibitor binding sites. The DAT model may guide DAT ligand QSAR studies, and rational design of novel DAT-binding therapeutics.     

Insights into transport mechanism from LeuT engineered to transport tryptophan

LeuT is a bacterial homologue of the neurotransmitter:sodium symporter (NSS) family and, being the only NSS member to have been structurally characterized by X-ray crystallography, is a model protein for studying transporter structure and mechanism. Transport activity in LeuT was hypothesized to require structural transitions between open-to-out and occluded conformations dependent upon protein:ligand binding complementarity. Here, using crystallographic and functional analysis, we show that binding site modification produces changes in both structure and activity that are consistent with complementarity-dependent structural transitions to the occluded state. The mutation I359Q converts the activity of tryptophan from inhibitor to transportable substrate. This mutation changes the local environment of the binding site, inducing the bound tryptophan to adopt a different conformer than in the wild-type complex. Instead of trapping the transporter open, tryptophan binding now allows the formation of an occluded state. Thus, transport activity is correlated to the ability of the ligand to promote the structural transition to the occluded state, a step in the transport cycle that is dependent on protein:ligand complementarity in the central binding site.     

Threonine 67 is a key component in the coupling of the NSS amino acid transporter KAAT1

The crystallizations of the prokaryotic LeuT and of the eukaryotic DAT and SERT transporters represent important steps forward in the comprehension of the molecular physiology of Neurotransmitter: Sodium Symporters, although the molecular determinants of the coupling mechanism and of ion selectivity still remain to be fully elucidated. The insect NSS homologue KAAT1 exhibits unusual physiological features, such as the ability to use K⁺ as the driver ion, weak chloride dependence, and the ability of the driver ion to influence the substrate selectivity; these characteristics can help to define the molecular determinants of NSS function. Two non-conserved residues are present in the putative sodium binding sites of KAAT1: Ala 66, corresponding to Gly 20 in the Na2 site of LeuT, and Ser 68, corresponding to Ala 22 in the Na1 site. Thr 67 appears also to be significant since it is not conserved among NSS members, is present as threonine only in KAAT1 and in the paralogue CAATCH1 and, according to LeuT structure, is close to the amino acid binding site. Mutants of these residues were functionally characterized in Xenopus oocytes. The T67Y mutant exhibited uptake activity comparable to that of the wild type, but fully chloride-independent and with enhanced stereoselectivity. Interestingly, although dependent on the presence of sodium, the mutant showed reduced transport-associated currents, indicating uncoupling of the driver ion and amino acid fluxes. Thr 67 therefore appears to be a key component in the coupling mechanism, participating in a network that influences the cotransport of Na⁺ and the amino acid.     

The mechanism of a neurotransmitter:sodium symporter--inward release of Na+ and substrate is triggered by substrate in a second binding site

Eukaryotic neurotransmitter:sodium symporters (NSSs), targets for antidepressants and psychostimulants, terminate neurotransmission by sodium-driven reuptake. The crystal structure of LeuT(Aa), a prokaryotic NSS homolog, revealed an occluded state in which one leucine and two Na(+) ions are bound, but provided limited clues to the molecular mechanism of transport. Using steered molecular dynamics simulations, we explored the substrate translocation pathway of LeuT. We identified a second substrate binding site located in the extracellular vestibule comprised of residues shown recently to participate in binding tricyclic antidepressants. Binding and flux experiments showed that the two binding sites can be occupied simultaneously. The substrate in the secondary site allosterically triggers intracellular release of Na(+) and substrate from the primary site, thereby functioning as a "symport effector." Because tricyclic antidepressants bind differently to this secondary site, they do not promote substrate release from the primary site and thus act as symport uncouplers and inhibit transport.     

Elucidation of the glucose transport pathway in glucose transporter 4 via steered molecular dynamics simulations

Background

GLUT4 is a predominant insulin regulated glucose transporter expressed in major glucose disposal tissues such as adipocytes and muscles. Under the unstimulated state, GLUT4 resides within intracellular vesicles. Various stimuli such as insulin translocate this protein to the plasma membrane for glucose transport. In the absence of a crystal structure for GLUT4, very little is known about the mechanism of glucose transport by this protein. Earlier we proposed a homology model for GLUT4 and performed a conventional molecular dynamics study revealing the conformational rearrangements during glucose and ATP binding. However, this study could not explain the transport of glucose through the permeation tunnel.

Methodology/Principal Findings

To elucidate the molecular mechanism of glucose transport and its energetic, a steered molecular dynamics study (SMD) was used. Glucose was pulled from the extracellular end of GLUT4 to the cytoplasm along the pathway using constant velocity pulling method. We identified several key residues within the tunnel that interact directly with either the backbone ring or the hydroxyl groups of glucose. A rotation of glucose molecule was seen near the sugar binding site facilitating the sugar recognition process at the QLS binding site.

Conclusions/Significance

This study proposes a possible glucose transport pathway and aids the identification of several residues that make direct interactions with glucose during glucose transport. Mutational studies are required to further validate the observation made in this study.     

Ion-Controlled Conformational Dynamics in the Outward-Open Transition from an Occluded State of LeuT

Neurotransmitter:sodium symporter (NSS) proteins are secondary Na⁺-driven active transporters that terminate neurotransmission by substrate uptake. Despite the availability of high-resolution crystal structures of a bacterial homolog of NSSs—Leucine Transporter (LeuT)—and extensive computational and experimental structure-function studies, unanswered questions remain regarding the transport mechanisms. We used microsecond atomistic molecular-dynamics (MD) simulations and free-energy computations to reveal ion-controlled conformational dynamics of LeuT in relation to binding affinity and selectivity of the more extracellularly positioned Na⁺ binding site (Na1 site). In the course of MD simulations starting from the occluded state with bound Na⁺, but in the absence of substrate, we find a spontaneous transition of the extracellular vestibule of LeuT into an outward-open conformation. The outward opening is enhanced by the absence of Na1 and modulated by the protonation state of the Na1-associated Glu-290. Consistently, the Na⁺ affinity for the Na1 site is inversely correlated with the extent of outward-open character and is lower than in the occluded state with bound substrate; however, the Na1 site retains its selectivity for Na⁺ over K⁺ in such conformational transitions. To the best of our knowledge, our findings shed new light on the Na⁺-driven transport cycle and on the symmetry in structural rearrangements for outward- and inward-open transitions.     

Identification of the allosteric regulatory site of insulysin

Background

Insulin degrading enzyme (IDE) is responsible for the metabolism of insulin and plays a role in clearance of the Aβ peptide associated with Alzheimer''s disease. Unlike most proteolytic enzymes, IDE, which consists of four structurally related domains and exists primarily as a dimer, exhibits allosteric kinetics, being activated by both small substrate peptides and polyphosphates such as ATP.

Principal Findings

The crystal structure of a catalytically compromised mutant of IDE has electron density for peptide ligands bound at the active site in domain 1 and a distal site in domain 2. Mutating residues in the distal site eliminates allosteric kinetics and activation by a small peptide, as well as greatly reducing activation by ATP, demonstrating that this site plays a key role in allostery. Comparison of the peptide bound IDE structure (using a low activity E111F IDE mutant) with unliganded wild type IDE shows a change in the interface between two halves of the clamshell-like molecule, which may enhance enzyme activity by altering the equilibrium between closed and open conformations. In addition, changes in the dimer interface suggest a basis for communication between subunits.

Conclusions/Significance

Our findings indicate that a region remote from the active site mediates allosteric activation of insulysin by peptides. Activation may involve a small conformational change that weakens the interface between two halves of the enzyme.     

An intracellular interaction network regulates conformational transitions in the dopamine transporter

Neurotransmitter:sodium symporters (NSS)(1) mediate sodium-dependent reuptake of neurotransmitters from the synaptic cleft and are targets for many psychoactive drugs. The crystal structure of the prokaryotic NSS protein, LeuT, was recently solved at high resolution; however, the mechanistic details of regulation of the permeation pathway in this class of proteins remain unknown. Here we combine computational modeling and experimental probing in the dopamine transporter (DAT) to demonstrate the functional importance of a conserved intracellular interaction network. Our data suggest that a salt bridge between Arg-60 in the N terminus close to the cytoplasmic end of transmembrane segment (TM) 1 and Asp-436 at the cytoplasmic end of TM8 is stabilized by a cation-pi interaction between Arg-60 and Tyr-335 at the cytoplasmic end of TM6. Computational probing illustrates how the interactions may determine the flexibility of the permeation pathway, and mutagenesis within the network and results from assays of transport, as well as the state-dependent accessibility of a substituted cysteine in TM3, support the role of this network in regulating access between the substrate binding site and the intracellular milieu. The mechanism that emerges from these findings may be unique to the NSS family, where the local disruption of ionic interactions modulates the transition of the transporter between the outward- and inward-facing conformations.     

To be, or not to be two sites: that is the question about LeuT substrate binding

Transport proteins of the neurotransmitter sodium symporter (NSS) family regulate the extracellular concentration of several neurotransmitters in the central nervous system. The only member of this family for which atomic-resolution structural data are available is the prokaryotic homologue LeuT. This protein has been used as a model system to study the molecular mechanism of transport of the NSS family. In this Journal Club, we discuss two strikingly different LeuT transport mechanisms: one involving a single high-affinity substrate binding site and one recently proposed alternative involving two high-affinity substrate binding sites that are allosterically coupled.     

Coupled Global and Local Changes Direct Substrate Translocation by Neurotransmitter-Sodium Symporter Ortholog LeuT

Significant advances have been made in recent years in characterizing neurotransmitter:sodium symporter (NSS) family structure and function. Yet, many time-resolved events and intermediates that control the various stages of transport cycle remain to be elucidated. Whether NSSs harbor one or two sites for binding their substrates (neurotransmitters or amino acids), and what the role of the secondary site S2 is, if any, are still unresolved. Using molecular modeling and simulations for LeuT, a bacterial NSS, we present a comprehensive account of substrate-binding and -stabilization events, and subsequently triggered interactions leading to substrate (alanine) release. LeuT instantaneous conformation as it reconfigures from substrate-receiving (outward-facing) to -releasing (inward-facing) state appears to be a determinant of its affinity to bind substrate at site S2. In the outward-facing state, S1 robustly binds alanine and regulates subsequent redistribution of interactions to trigger extracellular gate closure; whereas S2 is only a transient binding site. The substrate-binding affinity at S2 increases in an intermediate close to inward-facing state. LeuT harbors the two substrate-binding sites, and small displacements of second substrate near S2 are observed to induce concerted small translocations in the substrate bound to primary site S1, although complete release requires collective structural rearrangements that fully expose the intracellular vestibule to the cytoplasm.     

Coupled Global and Local Changes Direct Substrate Translocation by?Neurotransmitter-Sodium Symporter Ortholog LeuT

Significant advances have been made in recent years in characterizing neurotransmitter:sodium symporter (NSS) family structure and function. Yet, many time-resolved events and intermediates that control the various stages of transport cycle remain to be elucidated. Whether NSSs harbor one or two sites for binding their substrates (neurotransmitters or amino acids), and what the role of the secondary site S2 is, if any, are still unresolved. Using molecular modeling and simulations for LeuT, a bacterial NSS, we present a comprehensive account of substrate-binding and -stabilization events, and subsequently triggered interactions leading to substrate (alanine) release. LeuT instantaneous conformation as it reconfigures from substrate-receiving (outward-facing) to -releasing (inward-facing) state appears to be a determinant of its affinity to bind substrate at site S2. In the outward-facing state, S1 robustly binds alanine and regulates subsequent redistribution of interactions to trigger extracellular gate closure; whereas S2 is only a transient binding site. The substrate-binding affinity at S2 increases in an intermediate close to inward-facing state. LeuT harbors the two substrate-binding sites, and small displacements of second substrate near S2 are observed to induce concerted small translocations in the substrate bound to primary site S1, although complete release requires collective structural rearrangements that fully expose the intracellular vestibule to the cytoplasm.     

Molecular dynamics of leucine and dopamine transporter proteins in a model cell membrane lipid bilayer

The dopamine transporter (DAT) operates via facilitated diffusion, harnessing an inward Na⁺ gradient to drive dopamine from the extracellular synaptic cleft to the neuron interior. The DAT is relevant to central nervous system disorders such as Parkinson disease and attention‐deficit hyperactivity disorder and is the primary site of action for the abused psychostimulants cocaine and amphetamines. Crystallization of a DAT homolog, the bacterial leucine transporter LeuT, provided the first reliable 3‐D DAT template. Here, the LeuT crystal structure and the DAT molecular model have been combined with their respective substrates, leucine and dopamine, in lipid bilayer molecular dynamics simulations toward tracking substrate movement along the protein's substrate/ion permeation pathway. Specifically, movement of residue pairs that comprise the “external gate” was followed as a function of substrate presence. The transmembrane (TM) 1 arginine‐TM 10 aspartate strut formed less readily in DAT compared with LeuT, with or without substrate present. For LeuT but not DAT, the addition of substrate enhanced the chances of forming the TM 1‐10 bridge. Also, movement of the fourth extracellular loop EL‐4 in the presence of substrate was more pronounced for DAT, the EL‐4 unwinding to a degree. The overall similarity between the LeuT and DAT molecular dynamics simulations indicated that LeuT was a legitimate model to guide DAT structure‐function predictions. There were, nevertheless, differences significant enough to allow for DAT‐unique insights, which may include how cocaine, methylphenidate (Ritalin, NIDA Drug Supply, Rockville, MD), and other DAT blockers are not recognized as substrates even though they can access the primary substrate binding pocket. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

Dopamine transporter binding is unaffected by L-DOPA administration in normal and MPTP-treated monkeys

Background

Radiotracer imaging of the presynaptic nigrostriatal dopaminergic system is used to assess disease progression in Parkinson''s disease (PD) and may provide a useful adjunct to clinical assessment during therapeutic trials of potential neuroprotective agents. Several clinical trials comparing dopamine agonists to L-DOPA or early vs. late L-DOPA have revealed differences between clinical assessment and imaging of the presynaptic dopaminergic system, hence questioning the comparability of these measures as neuroprotection outcome variables. Thus, results of these studies may have been affected by factors other than the primary biological process investigated.

Methodology/Principal Findings

We tested the possibility that L-DOPA might interfere with DAT binding. Post-mortem DAT binding was conducted in normal and MPTP-treated macaque monkeys that were administered L-DOPA, acutely or chronically. In parallel, DAT SPECT was conducted in MPTP-treated animals that were administered chronic L-DOPA. [99mTc]TRODAT-1 SPECT binding was similarly reduced in all MPTP monkeys regardless of L-DOPA treatment. L-DOPA had no significant effect on post-mortem DAT binding either in saline or in MPTP-lesioned animals.

Conclusions/Significance

These data indicate that L-DOPA does not induce modifications of DAT expression detectable by SPECT of by DAT binding autoradiography, suggesting that differences between clinical assessment and radiotracer imaging in clinical trials may not be specifically related to L-DOPA treatment.     

Molecular mechanism of ion-ion and ion-substrate coupling in the Na+-dependent leucine transporter LeuT

Ion-coupled transport of neurotransmitter molecules by neurotransmitter:sodium symporters (NSS) play an important role in the regulation of neuronal signaling. One of the major events in the transport cycle is ion-substrate coupling and formation of the high-affinity occluded state with bound ions and substrate. Molecular mechanisms of ion-substrate coupling and the corresponding ion-substrate stoichiometry in NSS transporters has yet to be understood. The recent determination of a high-resolution structure for a bacterial homolog of Na⁺/Cl⁻-dependent neurotransmitter transporters, LeuT, offers a unique opportunity to analyze the functional roles of the multi-ion binding sites within the binding pocket. The binding pocket of LeuT contains two metal binding sites. The first ion in site NA1 is directly coupled to the bound substrate (Leu) with the second ion in the neighboring site (NA2) only ∼7 Å away. Extensive, fully atomistic, molecular dynamics, and free energy simulations of LeuT in an explicit lipid bilayer are performed to evaluate substrate-binding affinity as a function of the ion load (single versus double occupancy) and occupancy by specific monovalent cations. It was shown that double ion occupancy of the binding pocket is required to ensure substrate coupling to Na⁺ and not to Li⁺ or K⁺ cations. Furthermore, it was found that presence of the ion in site NA2 is required for structural stability of the binding pocket as well as amplified selectivity for Na⁺ in the case of double ion occupancy.     

The structure of the NPC1L1 N-terminal domain in a closed conformation

Background

NPC1L1 is the molecular target of the cholesterol lowering drug Ezetimibe and mediates the intestinal absorption of cholesterol. Inhibition or deletion of NPC1L1 reduces intestinal cholesterol absorption, resulting in reduction of plasma cholesterol levels.

Principal Findings

Here we present the 2.8 Å crystal structure of the N-terminal domain (NTD) of NPC1L1 in the absence of cholesterol. The structure, combined with biochemical data, reveals the mechanism of cholesterol selectivity of NPC1L1. Comparison to the cholesterol free and bound structures of NPC1(NTD) reveals that NPC1L1(NTD) is in a closed conformation and the sterol binding pocket is occluded from solvent.

Conclusion

The structure of NPC1L1(NTD) reveals a degree of flexibility surrounding the entrance to the sterol binding pocket, suggesting a gating mechanism that relies on multiple movements around the entrance to the sterol binding pocket.     

Structural and biochemical studies of human 4-hydroxy-2-oxoglutarate aldolase: implications for hydroxyproline metabolism in primary hyperoxaluria

Background

4-hydroxy-2-oxoglutarate (HOG) aldolase is a unique enzyme in the hydroxyproline degradation pathway catalyzing the cleavage of HOG to pyruvate and glyoxylate. Mutations in this enzyme are believed to be associated with the excessive production of oxalate in primary hyperoxaluria type 3 (PH3), although no experimental data is available to support this hypothesis. Moreover, the identity, oligomeric state, enzymatic activity, and crystal structure of human HOGA have not been experimentally determined.

Methodology/Principal Findings

In this study human HOGA (hHOGA) was identified by mass spectrometry of the mitochondrial enzyme purified from bovine kidney. hHOGA performs a retro-aldol cleavage reaction reminiscent of the trimeric 2-keto-3-deoxy-6-phosphogluconate aldolases. Sequence comparisons, however, show that HOGA is related to the tetrameric, bacterial dihydrodipicolinate synthases, but the reaction direction is reversed. The 1.97 Å resolution crystal structure of hHOGA bound to pyruvate was determined and enabled the modeling of the HOG-Schiff base intermediate and the identification of active site residues. Kinetic analyses of site-directed mutants support the importance of Lys196 as the nucleophile, Tyr168 and Ser77 as components of a proton relay, and Asn78 and Ser198 as unique residues that facilitate substrate binding.

Conclusions/Significance

The biochemical and structural data presented support that hHOGA utilizes a type I aldolase reaction mechanism, but employs novel residue interactions for substrate binding. A mapping of the PH3 mutations identifies potential rearrangements in either the active site or the tetrameric assembly that would likely cause a loss in activity. Altogether, these data establish a foundation to assess mutant forms of hHOGA and how their activity could be pharmacologically restored.     

Molecular Dynamics Simulation of the Allosteric Regulation of eIF4A Protein from the Open to Closed State,Induced by ATP and RNA Substrates

Background

Eukaryotic initiation factor 4A (eIF4A) plays a key role in the process of protein translation initiation by facilitating the melting of the 5′ proximal secondary structure of eukaryotic mRNA for ribosomal subunit attachment. It was experimentally postulated that the closed conformation of the eIF4A protein bound by the ATP and RNA substrates is coupled to RNA duplex unwinding to promote protein translation initiation, rather than an open conformation in the absence of ATP and RNA substrates. However, the allosteric process of eIF4A from the open to closed state induced by the ATP and RNA substrates are not yet fully understood.

Methodology

In the present work, we constructed a series of diplex and ternary models of the eIF4A protein bound by the ATP and RNA substrates to carry out molecular dynamics simulations, free energy calculations and conformation analysis and explore the allosteric properties of eIF4A.

Results

The results showed that the eIF4A protein completes the conformational transition from the open to closed state via two allosteric processes of ATP binding followed by RNA and vice versa. Based on cooperative allosteric network analysis, the ATP binding to the eIF4A protein mainly caused the relative rotation of two domains, while the RNA binding caused the proximity of two domains via the migration of RNA bases in the presence of ATP. The cooperative binding of ATP and RNA for the eIF4A protein plays a key role in the allosteric transition.     

Substrate binding and formation of an occluded state in the leucine transporter

Translocation through the extracellular vestibule and binding of leucine in the leucine transporter (LeuT) have been studied with molecular dynamics simulations. More than 0.1 μs of all-atom molecular dynamics simulations have been performed on different combinations of LeuT, bound substrate, and bound structural Na⁺ ions to describe molecular events involved in substrate binding and in the formation of the occluded state and to investigate the dynamics of this state. Three structural features are found to be directly involved in the initial steps of leucine transport: a Na⁺ ion directly coordinated to leucine (Na-1), two aromatic residues closing the binding site toward the extracellular vestibule (Tyr-108 and Phe-253), and a salt bridge in the extracellular vestibule (Arg-30 and Asp-404). These features account for observed differences between simulations of LeuT with and without bound substrate and for a possible pathway for leucine binding and thereby formation of the occluded LeuT binding site.     

Auto-FACE: An NMR Based Binding Site Mapping Program for Fast Chemical Exchange Protein-Ligand Systems

Background

Nuclear Magnetic Resonance (NMR) spectroscopy offers a variety of experiments to study protein-ligand interactions at atomic resolution. Among these experiments, N Heteronuclear Single Quantum Correlation (HSQC) experiment is simple, less time consuming and highly informative in mapping the binding site of the ligand. The interpretation of N HSQC becomes ambiguous when the chemical shift perturbations are caused by non-specific interactions like allosteric changes and local structural rearrangement. Under such cases, detailed chemical exchange analysis based on chemical shift perturbation will assist in locating the binding site accurately.

Methodology/Principal Findings

We have automated the mapping of binding sites for fast chemical exchange systems using information obtained from N HSQC spectra of protein serially titrated with ligand of increasing concentrations. The automated program Auto-FACE (Auto-FAst Chemical Exchange analyzer) determines the parameters, e.g. rate of change of perturbation, binding equilibrium constant and magnitude of chemical shift perturbation to map the binding site residues. Interestingly, the rate of change of perturbation at lower ligand concentration is highly sensitive in differentiating the binding site residues from the non-binding site residues. To validate this program, the interaction between the protein and the ligand BH3I-1 was studied. Residues in the hydrophobic BH3 binding groove of were easily identified to be crucial for interaction with BH3I-1 from other residues that also exhibited perturbation. The geometrically averaged equilibrium constant () calculated for the residues present at the identified binding site is consistent with the values obtained by other techniques like isothermal calorimetry and fluorescence polarization assays (). Adjacent to the primary site, an additional binding site was identified which had an affinity of 3.8 times weaker than the former one. Further NMR based model fitting for individual residues suggest single site model for residues present at these binding sites and two site model for residues present between these sites. This implies that chemical shift perturbation can represent the local binding event much more accurately than the global binding event.

Conclusion/Significance

Detail NMR chemical shift perturbation analysis enabled binding site residues to be distinguished from non-binding site residues for accurate mapping of interaction site in complex fast exchange system between small molecule and protein. The methodology is automated and implemented in a program called “Auto-FACE”, which also allowed quantitative information of each interaction site and elucidation of binding mechanism.     

Catalytic Water Co-Existing with a Product Peptide in the Active Site of HIV-1 Protease Revealed by X-Ray Structure Analysis

Background

It is known that HIV-1 protease is an important target for design of antiviral compounds in the treatment of Acquired Immuno Deficiency Syndrome (AIDS). In this context, understanding the catalytic mechanism of the enzyme is of crucial importance as transition state structure directs inhibitor design. Most mechanistic proposals invoke nucleophilic attack on the scissile peptide bond by a water molecule. But such a water molecule coexisting with any ligand in the active site has not been found so far in the crystal structures.

Principal Findings

We report here the first observation of the coexistence in the active site, of a water molecule WAT1, along with the carboxyl terminal product (Q product) peptide. The product peptide has been generated in situ through cleavage of the full-length substrate. The N-terminal product (P product) has diffused out and is replaced by a set of water molecules while the Q product is still held in the active site through hydrogen bonds. The position of WAT1, which hydrogen bonds to both the catalytic aspartates, is different from when there is no substrate bound in the active site. We propose WAT1 to be the position from where catalytic water attacks the scissile peptide bond. Comparison of structures of HIV-1 protease complexed with the same oligopeptide substrate, but at pH 2.0 and at pH 7.0 shows interesting changes in the conformation and hydrogen bonding interactions from the catalytic aspartates.

Conclusions/Significance

The structure is suggestive of the repositioning, during substrate binding, of the catalytic water for activation and subsequent nucleophilic attack. The structure could be a snap shot of the enzyme active site primed for the next round of catalysis. This structure further suggests that to achieve the goal of designing inhibitors mimicking the transition-state, the hydrogen-bonding pattern between WAT1 and the enzyme should be replicated.