Hydrophilic Pocket Research Articles

Binding isotope effects as a tool for distinguishing hydrophobic and hydrophilic binding sites of HIV-1 RT.

The current treatment for HIV-1 infected patients consists of a cocktail of inhibitors, in an attempt to improve the potency of the drugs by adding the possible effects of each supplied compound. In this contribution, nine different inhibitors of HIV-1 RT, one of the three key proteins responsible for the virus replication, have been selected to develop and test a computational protocol that allows getting a deep insight into the inhibitors' binding mechanism. The interaction between the inhibitors and the protein have been quantified by computing binding free energies through FEP calculations, while a more detailed characterization of the kind of inhibitor-protein interactions is based on frequency analysis of the ligands in the initial and final state, i.e. in solution and binding the protein. QM/MM calculation of heavy atoms ((13)C, (15)N, and (18)O) binding isotope effects (BIE) have been used to identify the binding sites of the different inhibitors. Specific interactions between the isotopically labeled atoms of the inhibitors and polar residues and magnesium cations on the hydrophilic pocket of the protein are responsible for the frequencies shifting that can be detected when comparing the IR spectra of the compounds in solution and in the protein. On the contrary, it seems that changes in vdW interactions from solution to the final state when the ligand is interacting with residues of the hydrophobic cavity, does not influence frequency modes and then no BIE are observed. Our results suggest that a proper computational protocol can be a valuable tool which in turn can be used to increase the efficiency of anti AIDS drugs.

Relative Domain Folding and Stability of a Membrane Transport Protein

There is a limited understanding of the folding of multidomain membrane proteins. Lactose permease (LacY) of Escherichia coli is an archetypal member of the major facilitator superfamily of membrane transport proteins, which contain two domains of six transmembrane helices each. We exploit chemical denaturation to determine the unfolding free energy of LacY and employ Trp residues as site-specific thermodynamic probes. Single Trp LacY mutants are created with the individual Trps situated at mirror image positions on the two LacY domains. The changes in Trp fluorescence induced by urea denaturation are used to construct denaturation curves from which unfolding free energies can be determined. The majority of the single Trp tracers report the same stability and an unfolding free energy of approximately +2kcal mol−1. There is one exception; the fluorescence of W33 at the cytoplasmic end of helix I on the N domain is unaffected by urea. In contrast, the equivalent position on the first helix, VII, of the C-terminal domain exhibits wild-type stability, with the single Trp tracer at position 243 on helix VII reporting an unfolding free energy of +2kcal mol−1. This indicates that the region of the N domain of LacY at position 33 on helix I has enhanced stability to urea, when compared the corresponding location at the start of the C domain. We also find evidence for a potential network of stabilising interactions across the domain interface, which reduces accessibility to the hydrophilic substrate binding pocket between the two domains.

Nanoparticle and Gelation Stabilized Functional Composites of an Ionic Salt in a Hydrophobic Polymer Matrix

Polymer composites consisted of small hydrophilic pockets homogeneously dispersed in a hydrophobic polymer matrix are important in many applications where controlled release of the functional agent from the hydrophilic phase is needed. As an example, a release of biomolecules or drugs from therapeutic formulations or release of salt in anti-icing application can be mentioned. Here, we report a method for preparation of such a composite material consisted of small KCOOH salt pockets distributed in the styrene-butadiene-styrene (SBS) polymer matrix and demonstrate its effectiveness in anti-icing coatings. The mixtures of the aqueous KCOOH and SBS-cyclohexane solutions were firstly stabilized by adding silica nanoparticles to the emulsions and, even more, by gelation of the aqueous phase by agarose. The emulsions were observed in optical microscope to check its stability in time and characterized by rheological measurements. The dry composite materials were obtained via casting the emulsions onto the glass substrates and evaporations of the organic solvent. Composite polymer films were characterized by water contact angle (WCA) measurements. The release of KCOOH salt into water and the freezing delay experiments of water droplets on dry composite films demonstrated their anti-icing properties. It has been concluded that hydrophobic and thermoplastic SBS polymer allows incorporation of the hydrophilic pockets/phases through our technique that opens the possibility for controlled delivering of anti-icing agents from the composite.

Structural Basis for the Assembly of Kinesin-5 into Bipolar Anti-Parallel Tetramers

Cell division is driven by organized mitotic spindles and requires diverse groups of molecular microtubule-based kinesin motor proteins. Kinesin-5 motors form a unique and conserved class of tetrameric bipolar motor proteins, consisting of dual motor domain sets oriented at opposite ends of a central rod-like structure. The bipolar organization is critical for kinesin-5 to crosslink and slide adjacent anti-parallel microtubules in opposite directions required for spindle elongation in anaphase. Here we present the atomic structure of Drosphila kinesin-5, KLF61F, bipolar assembly (BASS) domain, which explains how four molecules assemble to form a bipolar functional motor. Instead of the predicted dual traversing parallel coiled-coils, the BASS structure reveals a novel 26nm long antiparallel 4-helical bundle filament, where the basic assembly unit is an anti-parallel coiled-coil that originate from opposite ends of the bipolar structure, and folds onto a second unit in an anti-parallel manner. A striated pattern of hydrophobic and hydrophilic pockets precisely oriented monomers within the BASS tetramer to form this bipolar arrangement. Strikingly, the BASS N-termini undergo a helical exchange from anti-parallel bundle in the center to form homo-dimeric parallel helical coiled-coils at the poles. These parallel BASS N-termini are 100 degrees rotated compared to N-termini emerging from the opposite pole, providing a physical explanation for the known kinesin-5 preference for sliding anti-parallel versus parallel microtubules. We have generated a structural model for the kinesin-5 tetramer using the BASS structures and all electron microscopy that explains key features of the kinesin-5-driven sliding filament mechanism.

Molecular Simulation to Investigate the Cofactor Specificity for Pichia stipitis Xylose Reductase

Xylose is one of the most abundant carbohydrates in nature, and widely used to produce bioethanol via fermentation in industry. Xylulose can produce two key enzymes: xylose reductase and xylitol dehydrogenase. Owing to the disparate cofactor specificities of xylose reductase and xylitol dehydrogenase, intracellular redox imbalance is detected during the xylose fermentation, resulting in low ethanol yields. To overcome this barrier, a common strategy is applied to artificially modify the cofactor specificity of xylose reductase. In this study, we utilized molecular simulation approaches to construct a 3D (three-dimensional) structural model for the NADP-dependent Pichia stipitis xylose reductase (PsXR). Based on the 3D model, the favourable binding modes for both cofactors NAD and NADP were obtained using the flexible docking procedure and molecular dynamics simulation. Structural analysis of the favourable binding modes showed that the cofactor binding site of PsXR was composed of 3 major components: a hydrophilic pocket, a hydrophobic pocket as well as a linker channel between the aforementioned two pockets. The hydrophilic pocket could recognize the nicotinamide moiety of the cofactors by hydrogen bonding networks, while the hydrophobic pocket functioned to position the adenine moiety of the cofactors by hydrophobic and Π-Π stacking interactions. The linker channel contained some key residues for ligand-binding; their mutation could have impact to the specificity of PsXR. Finally, it was found that any of the two single mutations, K21A and K270N, might reverse the cofactor specificity of PsXR from major NADP- to NADdependent, which was further confirmed by the additional experiments. Our findings may provide useful insights into the cofactor specificity of PsXR, stimulating new strategies for better designing xylose reductase and improving ethanol production in industry.

Complex Structure of OspI and Ubc13: The Molecular Basis of Ubc13 Deamidation and Convergence of Bacterial and Host E2 Recognition

Ubc13 is an important ubiquitin-conjugating (E2) enzyme in the NF-κB signaling pathway. The Shigella effector OspI targets Ubc13 and deamidates Gln100 of Ubc13 to a glutamic acid residue, leading to the inhibition of host inflammatory responses. Here we report the crystal structure of the OspI-Ubc13 complex at 2.3 Å resolution. The structure reveals that OspI uses two differently charged regions to extensively interact with the α1 helix, L1 loop and L2 loop of Ubc13. The Gln100 residue is bound within the hydrophilic catalytic pocket of OspI. A comparison between Ubc13-bound and wild-type free OspI structures revealed that Ubc13 binding induces notable structural reassembly of the catalytic pocket, suggesting that substrate binding might be involved in the catalysis of OspI. The OspI-binding sites in Ubc13 largely overlap with the binding residues for host ubiquitin E3 ligases and a deubiquitinating enzyme, which suggests that the bacterial effector and host proteins exploit the same surface on Ubc13 for specific recognition. Biochemical results indicate that both of the differently charged regions in OspI are important for the interaction with Ubc13, and the specificity determinants in Ubc13 for OspI recognition reside in the distinct residues in the α1 helix and L2 region. Our study reveals the molecular basis of Ubc13 deamidation by OspI, as well as a convergence of E2 recognition by bacterial and host proteins.

Beta‐arc 2 hydrogen bond network within hemolysin A facilitates hemolytic activation

In the current study, the hemolytic exoprotein from Proteus mirabilis, termed HpmA, was cloned in a truncated form (HpmA265) and assayed to characterize its biological activity and thermal stability. A template‐assisted hemolytic assay measured the ability of native HpmA265 and six mutant variants to facilitate the refolding of denatured full‐length HpmA into its active conformation, allowing it to lyse RBCs. Circular dichroism spectroscopy was employed to estimate the denaturation temperature (Tm) of native HpmA265 and the six mutant variants (each of which harbored one or more serine substitutions). The mutations were located in the region of the protein called beta‐arc 2, which forms the base of the beta‐helical architecture of the protein, a structure common to many proteins secreted via the two‐partner secretion pathway in Gram‐negative bacteria. Tm changes ranged from −4.0 to −12.8°C among the mutants, and hemolytic activation (by template‐assisted refolding of full‐length HpmA) was markedly decreased in the triple and quadruple mutants, but activity remained near the native parameters for the single and double mutants. This suggests that an integrated network of hydrogen bonds in beta‐arc 2, with cooperative interactions in the hydrophilic pocket and the exterior carboxyamide ladder, lend substantial structural stability to HpmA265. This research is supported in part by NSF‐RUI 1050435 (to TW).

Conformationally Controlled Oligocholate Membrane Transporters: Learning through Water Play

Controlled translocation of molecules and ions across lipid membranes is the basis of numerous biological functions. Because synthetic systems can help researchers understand the more complex biological ones, many chemists have developed synthetic mimics of biological transporters. Both systems need to deal with similar fundamental challenges. In addition to providing mechanistic insights into transport mechanisms, synthetic transporters are useful in a number of applications including separation, sensing, drug delivery, and catalysis. In this Account, we present several classes of membrane transporters constructed in our laboratory from a facially amphiphilic building block, cholic acid. Our "molecular baskets" can selectively shuttle glucose across lipid membranes without transporting smaller sodium ions. We have also built oligocholate foldamers that transiently fold into helices with internal hydrophilic binding pockets to transport polar guests. Lastly, we describe amphiphilic macrocycles, which form transmembrane nanopores in lipid bilayers through the strong associative interactions of encapsulated water molecules. In addition to presenting the different transport properties of these oligocholate transporters, we illustrate how fundamental studies of molecular behavior in solution facilitate the creation of new and useful membrane transporters, despite the large difference between the two environments. We highlight the strong conformational effect of transporters. Because the conformation of a molecule often alters its size and shape, and the distribution of functional groups, conformational control can be used rationally to tune the property of a transporter. Finally, we emphasize that, whenever water is the solvent, its unique properties--small size, strong solvation for ionic functionalities, and an extraordinary cohesive energy density (i.e., total intermolecular interactions per unit volume)--tend to become critical factors to be considered. Purposeful exploitation of these solvent properties may be essential to the success of the supramolecular process involved--this is also the reason for the "learning through water play" in the title of this Account.

The Three-Dimensional Structural Basis of Type II Hyperprolinemia

Type II hyperprolinemia is an autosomal recessive disorder caused by a deficiency in Δ1-pyrroline-5-carboxylate dehydrogenase (P5CDH; also known as ALDH4A1), the aldehyde dehydrogenase that catalyzes the oxidation of glutamate semialdehyde to glutamate. Here, we report the first structure of human P5CDH (HsP5CDH) and investigate the impact of the hyperprolinemia-associated mutation of Ser352 to Leu on the structure and catalytic properties of the enzyme. The 2. 5-Å-resolution crystal structure of HsP5CDH was determined using experimental phasing. Structures of the mutant enzymes S352A (2.4 Å) and S352L (2.85 Å) were determined to elucidate the structural consequences of altering Ser352. Structures of the 93% identical mouse P5CDH complexed with sulfate ion (1.3 Å resolution), glutamate (1.5 Å), and NAD+ (1.5 Å) were determined to obtain high-resolution views of the active site. Together, the structures show that Ser352 occupies a hydrophilic pocket and is connected via water-mediated hydrogen bonds to catalytic Cys348. Mutation of Ser352 to Leu is shown to abolish catalytic activity and eliminate NAD+ binding. Analysis of the S352A mutant shows that these functional defects are caused by the introduction of the nonpolar Leu352 side chain rather than the removal of the Ser352 hydroxyl. The S352L structure shows that the mutation induces a dramatic 8-Å rearrangement of the catalytic loop. Because of this conformational change, Ser349 is not positioned to interact with the aldehyde substrate, conserved Glu447 is no longer poised to bind NAD+, and Cys348 faces the wrong direction for nucleophilic attack. These structural alterations render the enzyme inactive.

Ligand Recognition of Lipocalin-Type Prostaglandin D Synthase

Lipocalin-type prostaglandin (PG) D synthase (L-PGDS) is a multi functional protein that functions as a PGD2 synthesizing enzyme, a scavenger of various lipophilic ligands, and an amyloid β chaperone in the brain. L-PGDS is a member of the lipocalin superfamily and has the ability to bind various lipophilic molecules, such as prostanoid, retinoid, bile pigment, and amyloid β peptide. However, the molecular mechanism for this wide variety of ligand binding is not fully understood. In this study, we determined the NMR structures of recombinant mouse L-PGDS and a L-PGDS/PGH2 analog. L-PGDS has a typical lipocalin fold, consisting of an eight-stranded β-barrel and a single α-helix. The interior of the barrel forms a hydrophobic cavity and the upper end of the barrel is open. The size of the barrel was found to be larger than those of other lipocalins and the cavity contained two pockets. The results of NMR titration, kinetic and molecular docking experiments revealed that PGH2 and retinoic acid occupied the hydrophilic pockets 1 and 2, respectively. A structural comparison of the L-PGDS/PGH2 analog complex with apo-L-PGDS showed that the H2-helix, CD-loop, and EF-loop located at the upper end of the β-barrel undergo a conformational change and cover the entrance of the cavity upon U-46619 binding. These results indicate that the two binding sites in the large cavity are responsible for the broad ligand specificity of L-PGDS in its induced-fit mechanism.

Design, synthesis and evaluation of polar head group containing 2-keto-oxazole inhibitors of FAAH

2-α-Keto oxazoles containing polar head groups in their C5-side chains were designed as fatty acid amide hydrolase (FAAH) inhibitors. Variation in the spacer length resulted in submicromolar α-keto-oxazole FAAH inhibitor (IC50=436nM) presenting electrostatic stabilizing interactions between its polar head group contained in the C5-side chain and the hydrophilic pocket of the enzyme.

RtcB, a Novel RNA Ligase, Can Catalyze tRNA Splicing and HAC1 mRNA Splicing in Vivo

RtcB enzymes are novel RNA ligases that join 2',3'-cyclic phosphate and 5'-OH ends. The phylogenetic distribution of RtcB points to its candidacy as a tRNA splicing/repair enzyme. Here we show that Escherichia coli RtcB is competent and sufficient for tRNA splicing in vivo by virtue of its ability to complement growth of yeast cells that lack the endogenous "healing/sealing-type" tRNA ligase Trl1. RtcB also protects yeast trl1Δ cells against a fungal ribotoxin that incises the anticodon loop of cellular tRNAs. Moreover, RtcB can replace Trl1 as the catalyst of HAC1 mRNA splicing during the unfolded protein response. Thus, RtcB is a bona fide RNA repair enzyme with broad physiological actions. Biochemical analysis of RtcB highlights the uniqueness of its active site and catalytic mechanism. Our findings draw attention to tRNA ligase as a promising drug target.

Peptide aldehyde inhibitors challenge the substrate specificity of the SARS-coronavirus main protease

SARS coronavirus main protease (SARS-CoV M pro) is essential for the replication of the virus and regarded as a major antiviral drug target. The enzyme is a cysteine protease, with a catalytic dyad (Cys-145/His-41) in the active site. Aldehyde inhibitors can bind reversibly to the active-site sulfhydryl of SARS-CoV M pro. Previous studies using peptidic substrates and inhibitors showed that the substrate specificity of SARS-CoV M pro requires glutamine in the P1 position and a large hydrophobic residue in the P2 position. We determined four crystal structures of SARS-CoV M pro in complex with pentapeptide aldehydes (Ac-ESTLQ-H, Ac-NSFSQ-H, Ac-DSFDQ-H, and Ac-NSTSQ-H). Kinetic data showed that all of these aldehydes exhibit inhibitory activity towards SARS-CoV M pro, with K i values in the μM range. Surprisingly, the X-ray structures revealed that the hydrophobic S2 pocket of the enzyme can accommodate serine and even aspartic-acid side-chains in the P2 positions of the inhibitors. Consequently, we reassessed the substrate specificity of the enzyme by testing the cleavage of 20 different tetradecapeptide substrates with varying amino-acid residues in the P2 position. The cleavage efficiency for the substrate with serine in the P2 position was 160-times lower than that for the original substrate (P2 = Leu); furthermore, the substrate with aspartic acid in the P2 position was not cleaved at all. We also determined a crystal structure of SARS-CoV M pro in complex with aldehyde Cm-FF-H, which has its P1-phenylalanine residue bound to the relatively hydrophilic S1 pocket of the enzyme and yet exhibits a high inhibitory activity against SARS-CoV M pro, with K i = 2.24 ± 0.58 μM. These results show that the stringent substrate specificity of the SARS-CoV M pro with respect to the P1 and P2 positions can be overruled by the highly electrophilic character of the aldehyde warhead, thereby constituting a deviation from the dogma that peptidic inhibitors need to correspond to the observed cleavage specificity of the target protease.

Ligand-based 3-D pharmacophore generation and molecular docking of mTOR kinase inhibitors

The 3-D structure of the human mTOR kinase domain was modeled based on the crystal structure of PI3Kγ using comparative modeling methods, and the ATP-binding site of mTOR was characterized. The mTOR kinase 3-D model structure is similar to the structure of the PI3Kγ kinase domain, and exhibits great similarity to PI3Kγ at the active site of the kinase. Pharmacophore generation, the docking of mTOR inhibitors, and molecular dynamics (MD) simulations of mTOR-inhibitor docked complexes were carried out in this work. The best pharmacophore model generated from 27 ATP-competitive mTOR inhibitors comprised two hydrogen-bond acceptors, one aromatic ring, and one hydrophobic feature. These 27 inhibitors were docked into the ATP-binding site comprising the DFG motif, and the interactions in each protein-inhibitor complex were characterized. Mapping the pharmacophore model onto the docked inhibitors explained the specificity of the features of the pharmacophore and how they were arranged inside the active site of mTOR kinase. MD studies revealed important structural features, such as the large hydrophobic pocket "HP" and hydrophilic pocket "A1," and the solvent-exposed hydrophilic pocket "A2" at the active site of mTOR. Our results provide structural models of mTOR-inhibitor complexes and clues that should aid in the design of better mTOR kinase inhibitors.

Structural and Functional Analysis of Two New Positive Allosteric Modulators of GluA2 Desensitization and Deactivation

At the dimer interface of the extracellular ligand-binding domain of α-amino-3-hydroxy-5-methylisoxazole-4-propionic acid (AMPA) receptors a hydrophilic pocket is formed that is known to interact with two classes of positive allosteric modulators, represented by cyclothiazide and the ampakine 2H,3H,6aH-pyrrolidino(2,1-3',2')1,3-oxazino(6',5'-5,4)benzo(e)1,4-dioxan-10-one (CX614). Here, we present structural and functional data on two new positive allosteric modulators of AMPA receptors, phenyl-1,4-bis-alkylsulfonamide (CMPDA) and phenyl-1,4-bis-carboxythiophene (CMPDB). Crystallographic data show that these compounds bind within the modulator-binding pocket and that substituents of each compound overlap with distinct moieties of cyclothiazide and CX614. The goals of the present study were to determine 1) the degree of modulation by CMPDA and CMPDB of AMPA receptor deactivation and desensitization; 2) whether these compounds are splice isoform-selective; and 3) whether predictions of mechanism of action could be inferred by comparing molecular interactions between the ligand-binding domain and each compound with those of cyclothiazide and CX614. CMPDB was found to be more isoform-selective than would be predicted from initial binding assays. It is noteworthy that these new compounds are both more potent and more effective and may be more clinically relevant than the AMPA receptor modulators described previously.

Bilirubin oxidase from Myrothecium verrucaria: X-ray determination of the complete crystal structure and a rational surface modification for enhanced electrocatalytic O2 reduction

The blue multi-copper oxidase bilirubin oxidase (BOx) from the ascomycete plant pathogen Myrothecium verrucaria (Mv) efficiently catalyses the oxidation of bilirubin to biliverdin, with the concomitant reduction of O(2) to water, a reaction of considerable interest for low-temperature bio-fuel cell applications. We have solved the complete X-ray determined structure of Mv BOx at 2.4 Å resolution, using molecular replacement with the Spore Coat Protein A (CotA) enzyme from Bacillus subtilis (PDB code 1GSK) as a template. The structure reveals an unusual environment around the blue type 1 copper (T1 Cu) that includes two non-coordinating hydrophilic amino acids, asparagine and threonine. The presence of a long, narrow and hydrophilic pocket near the T1 Cu suggests that structure of the substrate-binding site is dynamically determined in vivo. We show that the interaction between the binding pocket of Mv BOx and its highly conjugated natural organic substrate, bilirubin, can be used to stabilise the enzyme on a pyrolytic graphite electrode, more than doubling its electrocatalytic activity relative to the current obtained by simple adsorption of the protein to the carbon surface.

On the mechanisms of bananin activity against severe acute respiratory syndrome coronavirus

In a previous study, severe acute respiratory syndrome coronavirus (SARS‐CoV) was cultured in the presence of bananin, an effective adamantane‐related molecule with antiviral activity. In the present study, we show that all bananin‐resistant variants exhibit mutations in helicase and membrane protein, although no evidence of bananin interference on their mutual interaction has been found. A structural analysis on protein sequence mutations found in SARS‐CoV bananin‐resistant variants was performed. The S259/L mutation of SARS‐CoV helicase is always found in all the identified bananin‐resistant variants, suggesting a primary role of this mutation site for bananin activity. From a structural analysis of SARS‐CoV predicted helicase structure, S259 is found in a hydrophilic surface pocket, far from the enzyme active sites and outside the helicase dimer interface. The S/L substitution causes a pocket volume reduction that weakens the interaction between bananin and SARS‐CoV mutated helicase, suggesting a possible mechanism for bananin antiviral activity.

Structural insight into human variegate porphyria disease

Human protoporphyrinogen IX oxidase (hPPO), a mitochondrial inner membrane protein, converts protoporphyrinogen IX to protoporphyrin IX in the heme biosynthetic pathway. Mutations in the hPPO gene cause the inherited human disease variegate porphyria (VP). In this study, we report the crystal structure of hPPO in complex with the coenzyme flavin adenine dinucleotide (FAD) and the inhibitor acifluorfen at a resolution of 1.9 Å. The structural and biochemical analyses revealed the molecular details of FAD and acifluorfen binding to hPPO as well as the interactions of the substrate with hPPO. Structural analysis and gel chromatography indicated that hPPO is a monomer rather than a homodimer in vitro. The founder-effect mutation R59W in VP patients is most likely caused by a severe electrostatic hindrance in the hydrophilic binding pocket involving the bulky, hydrophobic indolyl ring of the tryptophan. Forty-seven VP-causing mutations were purified by chromatography and kinetically characterized in vitro. The effect of each mutation was demonstrated in the high-resolution crystal structure.

Binding and Orientation of Tricyclic Antidepressants within the Central Substrate Site of the Human Serotonin Transporter

Tricyclic antidepressants (TCAs) have been used for decades, but their orientation within and molecular interactions with their primary target is yet unsettled. The recent finding of a TCA binding site in the extracellular vestibule of the bacterial leucine transporter 11 A above the central site has prompted debate about whether this vestibular site in the bacterial transporter is applicable to binding of antidepressants to their relevant physiological target, the human serotonin transporter (hSERT). We present an experimentally validated structural model of imipramine and analogous TCAs in the central substrate binding site of hSERT. Two possible binding modes were observed from induced fit docking calculations. We experimentally validated a single binding mode by combining mutagenesis of hSERT with uptake inhibition studies of different TCA analogs according to the paired mutation ligand analog complementation paradigm. Using this experimental method, we identify a salt bridge between the tertiary aliphatic amine and Asp(98). Furthermore, the 7-position of the imipramine ring is found vicinal to Phe(335), and the pocket lined by Ala(173) and Thr(439) is utilized by 3-substituents. These protein-ligand contact points unambiguously orient the TCA within the central binding site and reveal differences between substrate binding and inhibitor binding, giving important clues to the inhibition mechanism. Consonant with the well established competitive inhibition of uptake by TCAs, the resulting binding site for TCAs in hSERT is fully overlapping with the serotonin binding site in hSERT and dissimilar to the low affinity noncompetitive TCA site reported in the leucine transporter (LeuT).

Ultrafast Dynamics of Diatomic Ligand Binding to Nitrophorin 4

Nitrophorin 4 (NP4) is a heme protein that stores and delivers nitric oxide (NO) through pH-sensitive conformational change. This protein uses the ferric state of a highly ruffled heme to bind NO tightly at low pH and release it at high pH. In this work, the rebinding kinetics of NO and CO to NP4 are investigated as a function of iron oxidation state and the acidity of the environment. The geminate recombination process of NO to ferrous NP4 at both pH 5 and pH 7 is dominated by a single approximately 7 ps kinetic phase that we attribute to the rebinding of NO directly from the distal pocket. The lack of pH dependence explains in part why NP4 cannot use the ferrous state to fulfill its function. The kinetic response of ferric NP4NO shows two distinct phases. The relative geminate amplitude of the slower phase increases dramatically as the pH is raised from 5 to 8. We assign the fast phase of NO rebinding to a conformation of the ferric protein with a closed hydrophobic pocket. The slow phase is assigned to the protein in an open conformation with a more hydrophilic heme pocket environment. Analysis of the ultrafast kinetics finds the equilibrium off-rate of NO to be proportional to the open state population as well as the pH-dependent amplitude of escape from the open pocket. When both factors are considered, the off-rate increases by more than an order of magnitude as the pH changes from 5 to 8. The recombination of CO to ferrous NP4 is observed to have a large nonexponential geminate amplitude with rebinding time scales of approximately 10(-11)-10(-9) s at pH 5 and approximately 10(-10)-10(-8) s at pH 7. The nonexponential CO rebinding kinetics at both pH 5 and pH 7 are accounted for using a simple model that has proven effective for understanding CO binding in a variety of other heme systems (Ye, X.; et al. Proc. Natl. Acad. Sci. U.S.A. 2007, 104, 14682).