Molecular Basis of Substrate Polyspecificity of the Candida albicans Mdr1p Multidrug/H+ Antiporter

The molecular basis of polyspecificity of Mdr1p, a major drug/H⁺ antiporter of Candida albicans, is not elucidated. We have probed the nature of the drug-binding pocket by performing systematic mutagenesis of the 12 transmembrane segments. Replacement of the 252 amino acid residues with alanine or glycine yielded 2/3 neutral mutations while 1/3 led to the complete or selective loss of resistance to drugs or substrates transported by the pump. Using the GlpT-based 3D–model of Mdr1p, we roughly categorized these critical residues depending on their type and localization, 1°/ main structural impact (“S” group), 2°/ exposure to the lipid interface (“L” group), 3°/ buried but not facing the main central pocket, inferred as critical for the overall H⁺/drug antiport mechanism (“M” group) and finally 4°/ buried and facing the main central pocket (“B” group). Among “B” category, 13 residues were essential for the large majority of drugs/substrates, while 5 residues were much substrate-specific, suggesting a role in governing polyspecificity (P group). 3D superposition of the substrate-specific MFS Glut1 and XylE with the MDR substrate-polyspecific MdfA and Mdr1p revealed that the B group forms a common substrate interaction core while the P group is only found in the 2 MDR MFS transporters, distributed into 3 areas around the B core. This specific pattern has let us to propose that the structural basis for polyspecificity of MDR MFS transporters is the extended capacity brought by residues located at the periphery of a binding core to accomodate compounds differing in size and type.     

Enhanced enantioselectivity of a carboxyl esterase from Rhodobacter sphaeroides by directed evolution

The present work created an esterase variant from Rhodobacter sphaeroides (RspE) with enhanced selectivity in hydrolytic kinetic resolutions by directed evolution. A “model” substrate, methyl mandelate, was introduced in the high-throughput screening procedure. E values of a variant CH (Asn62Cys/Leu145His) for six different esters were 10–83, which were a relative improvement compared to 2–20 for the wild type. Our subsequent crystal structure interpretation and molecular dynamics simulations helped shed light on the source of enantioselectivity modified by directed evolution. Though mutations displayed no “direct” interaction with the substrate, they were hypothesized to strengthen the intramolecular interaction in the catalytic cavity of variant. Conformation analysis revealed that the enhanced enantioselectivity of variant CH for the seven substrates applied in this study was derived from the decrease in size of the substrate binding pocket.     

High-resolution structure of phosphoketolase from Bifidobacterium longum determined by cryo-EM single-particle analysis

In bifidobacteria, phosphoketolase (PKT) plays a key role in the central hexose fermentation pathway called “bifid shunt.” The three-dimensional structure of PKT from Bifidobacterium longum with co-enzyme thiamine diphosphate (ThDpp) was determined at 2.1 Å resolution by cryo-EM single-particle analysis using 196,147 particles to build up the structural model of a PKT octamer related by D₄ symmetry. Although the cryo-EM structure of PKT was almost identical to the X-ray crystal structure previously determined at 2.2 Å resolution, several interesting structural features were observed in the cryo-EM structure. Because this structure was solved at relatively high resolution, it was observed that several amino acid residues adopt multiple conformations. Among them, Q546–D547–H548–N549 (the QN-loop) demonstrate the largest structural change, which seems to be related to the enzymatic function of PKT. The QN-loop is at the entrance to the substrate binding pocket. The minor conformer of the QN-loop is similar to the conformation of the QN-loop in the crystal structure. The major conformer is located further from ThDpp than the minor conformer. Interestingly, the major conformer in the cryo-EM structure of PKT resembles the corresponding loop structure of substrate-bound Escherichia coli transketolase. That is, the minor and major conformers may correspond to “closed” and “open” states for substrate access, respectively. Moreover, because of the high-resolution analysis, many water molecules were observed in the cryo-EM structure of PKT. Structural features of the water molecules in the cryo-EM structure are discussed and compared with water molecules observed in the crystal structure.     

The K-path entrance in cytochrome c oxidase is defined by mutation of E101 and controlled by an adjacent ligand binding domain

Three mutant forms of Rhodobacter sphaeroides cytochrome c oxidase (RsCcO) were created to test for multiple K-path entry sites (E101W), the existence of an “upper ligand site” (M350?W), and the nature and binding specificity of the “lower ligand site” (P315W/E101A) in the region of a crystallographically-defined deoxycholate at the K-path entrance. The effects of inhibitory and stimulatory detergents (dodecyl maltoside and Tween20) on these mutants are presented, as well as competition with other ligands, including the potentially physiologically relevant ligands cholesterol and retinoic acid. Ligands are shown to be able to compete with natural lipids to affect the activity of membrane-bound RsCcO. Results point to a single K-path entrance site at E101, with a single ligand binding pocket proximal to the entrance. The affinity of this pocket for amphipathic ligands is enhanced by removal of the E101 carboxyl and blocked by substituting a tryptophan in this area. A new crystal structure of the E101A mutant of RsCcO is presented that illustrates the structural basis of these results, showing that the loss of the E101 carboxyl creates a more hydrophobic groove consistent with altered ligand affinities.     

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

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

Free energy analysis of ω-transaminase reactions to dissect how the enzyme controls the substrate selectivity

ω-Transaminase (ω-TA) is the only naturally occurring enzyme allowing asymmetric amination of ketones for production of chiral amines. The active site of the enzyme was proposed to consist of two differently sized substrate binding pockets and the stringent steric constraint in the small pocket has presented a significant challenge to production of structurally diverse chiral amines. To provide a mechanistic understanding of how the (S)-specific ω-TA from Paracoccus denitrificans achieves the steric constraint in the small pocket, we developed a free energy analysis enabling quantification of individual contributions of binding and catalytic steps to changes in the total activation energy caused by structural differences in the substrate moiety that is to be accommodated by the small pocket. The analysis exploited kinetic and thermodynamic investigations using structurally similar substrates and the structural differences among substrates were regarded as probes to assess how much relative destabilizations of the reaction intermediates, i.e. the Michaelis complex and the transition state, were induced by the slight change of the substrate moiety inside the small pocket. We found that ≈80% of changes in the total activation energy resulted from changes in the enzyme-substrate binding energy, indicating that substrate selectivity in the small pocket is controlled predominantly by the binding step (K_M) rather than the catalytic step (k_cat). In addition, we examined the pH dependence of the kinetic parameters and the pH profiles of the K_M and k_cat values suggested that key active site residues involved in the binding and catalytic steps are decoupled. Taken together, these findings suggest that the active site residues forming the small pocket are mainly engaged in the binding step but not significantly involved in the catalytic step, which may provide insights into how to design a rational strategy for engineering of the small pocket to relieve the steric constraint toward bulky substituents.     

Disparity in productive binding mode of the slow-reacting enantiomer determines the novel catalytic behavior of Candida antarctica lipase B

In the context of specifying the origin of enzyme enantioselectivity, the present study explores the lipase enantioselectivity towards secondary alcohols of similar structure from the perspective of substrate binding. By carrying out molecular mechanics minimization as well as molecular dynamics simulation on tetrahedral reaction intermediates which are used as a model of transition state, we identify an unconventional productive binding mode (PBM)—M/H permutation type for Candida antarctica lipase B (CALB). The in silico results also indicate that different PBMs of the slow-reacting enantiomer do exist in one lipase even when there is little structural differences between substrates, e.g. compounds with Ph or CH₂CH₂Ph group display the M/H permutation type PBM while molecules with CH₂Ph show the M/L permutation type PBM. By relating the PBMs of substrates to the experimentally determined E-values obtained by Hoff et al. [16], we find that disparity in PBM of the slow-reacting enantiomer determines why E-values of substrates with CH₂Ph were lower than E-values of substrates with Ph or CH₂CH₂Ph group. The modeling results also suggest that the “pushed aside” effect of the F atom and Br atom accommodates the medium size substituent of the substrate better in the stereospecificity pocket of the enzyme.     

Mutational Mapping and Modeling of the Binding Site for (S)-Citalopram in the Human Serotonin Transporter

The serotonin transporter (SERT) regulates extracellular levels of the neurotransmitter serotonin (5-hydroxytryptamine) in the brain by facilitating uptake of released 5-hydroxytryptamine into neuronal cells. SERT is the target for widely used antidepressant drugs, including imipramine, fluoxetine, and (S)-citalopram, which are competitive inhibitors of the transport function. Knowledge of the molecular details of the antidepressant binding sites in SERT has been limited due to lack of structural data on SERT. Here, we present a characterization of the (S)-citalopram binding pocket in human SERT (hSERT) using mutational and computational approaches. Comparative modeling and ligand docking reveal that (S)-citalopram fits into the hSERT substrate binding pocket, where (S)-citalopram can adopt a number of different binding orientations. We find, however, that only one of these binding modes is functionally relevant from studying the effects of 64 point mutations around the putative substrate binding site. The mutational mapping also identify novel hSERT residues that are crucial for (S)-citalopram binding. The model defines the molecular determinants for (S)-citalopram binding to hSERT and demonstrates that the antidepressant binding site overlaps with the substrate binding site.     

Structural basis for the substrate specificity of PepA from Streptococcus pneumoniae, a dodecameric tetrahedral protease

Regulated cytosolic proteolysis is one of the key cellular processes ensuring proper functioning of a cell. M42 family proteases show a broad spectrum of substrate specificities, but the structural basis for such diversity of the substrate specificities is lagging behind biochemical data. Here we report the crystal structure of PepA from Streptococcus pneumoniae, a glutamyl aminopeptidase belonging to M42 family (SpPepA). We found that Arg-257 in the substrate binding pocket is strategically positioned so that Arg-257 can make electrostatic interactions with the acidic residue of a substrate at its N-terminus. Structural comparison of the substrate binding pocket of the M42 family proteases, along with the structure-based multiple sequence alignment, argues that the appropriate electrostatic interactions contribute to the selective substrate specificity of SpPepA.     

Bovine factor Xa and bovine trypsin: A comparison with respect to inhibition by some amidines and guanidines

The enzymatic activity of activated bovine blood clotting factor X toward the synthetic substrate N α-benzoyl-l-arginine ethyl ester and the inhibitory effects of a series of low molecular weight synthetic aromatic amidine and guanidine compounds on that activity were studied using the steady-state kinetic method. The kinetic parameters, K_m and κ_cat, and the apparent dissociation constant K_i^′ for each inhibitor, were determined for activated factor X hydrolysis of Bz-Arg-OEt at 37 °C, pH 7.8 in 0.1 n NaCl and 0.001 m CaCl₂. The same constants were determined for bovine β-trypsin under identical conditions. Comparison of kinetic constants determined for both enzymes shows that activated factor X binds the substrate Bz-Arg-OEt less efficiently than β-trypsin by several orders of magnitude. However, binding of the inhibitors benzamidine, p-aminobenzamidine, pentamidine, M&B 4596, phenylguanidine, and p-guanidinobenzoic acid is similar for both enzymes. The results indicate that these two closely related serine proteases differ little in the structural arrangement and accessibility of the anionic “pocket” at which these inhibitors bind. The large differences observed with respect to substrate binding activity probably reflect substantial structural differences between the two enzymes at secondary sites adjacent to the primary anionic site.     

Functional characterization of a novel tropinone reductase-like gene in Dendrobium nobile Lindl

Dendrobium nobile, a herbal medicine plant, contains many important alkaloids and other secondary metabolites with pharmacological and clinical effects. However, the biosynthetic pathway of these secondary metabolites is largely unknown. In present study, a cDNA sequence (DnTR2) that encodes a peptide with high similarity to known tropinone reductase (TR) was cloned from D. nobile Lindl. Sequence comparison and phylogenetic analysis showed that DnTR2 was evolutionarily distant from those well-characterized subgroups of TRs. qRT-PCR revealed that DnTR2 was expressed constitutively in all three vegetative organs (leaves, stems and roots) and was regulated by methyl jasmonate (MeJA), salicylic acid (SA) and nitrogen oxide (NO). Catalytic activity analysis using recombinant protein found that DnTR2 was not able to reduce tropinone, but reduced the two structural analogs of tropinone, 3-quinuclidinone hydrochloride and 4-methylcyclohexanone. Structural modeling and comparison suggested that the substrate specificity of TRs may not be determined by their phylogenetic relationships but by the amino acids that compose the substrate binding pocket. To verify this hypothesis, a site-directed mutagenesis was performed and it successfully restored the DnTR2 with tropinone reduction activity. Our results also showed that the substrate specificity of TRs was determined by a few residues that compose the substrate binding pocket which may have an important role for directed selecting of TRs with designated substrate specificities.     

Lipid-dependent sequential allosteric activation of heat-sensing TRPV1 channels by anchor-stereoselective “hot” vanilloid compounds and analogs

Both a silent resident phosphatidylinositol lipid and a “hot” vanilloid agonist capsaicin or resiniferatoxin have been shown to share the same inter-subunit binding pocket between a voltage sensor like domain and a pore domain in TRPV1. However, how the vanilloid competes off the resident lipid for allosteric TRPV1 activation is unknown. Here, the in silico research suggested that anchor-stereoselective sequential cooperativity between an initial recessive transient silent weak ligand binding site and a subsequent dominant steady-state strong ligand binding site in the vanilloid pocket may facilitate the lipid release for allosteric activation of TRPV1 by vanilloids or analogs upon non-covalent interactions. Thus, the resident lipid may play a critical role in allosteric activation of TRPV1 by vanilloid compounds and analogs.     

A crystal structure of the human protein kinase Mps1 reveals an ordered conformation of the activation loop

Monopolar spindle 1 (Mps1) is a dual-specificity protein kinase, orchestrating faithful chromosome segregation during mitosis. All reported structures of the Mps1 kinase adopt the hallmarks of an inactive conformation, which includes a mostly disordered activation loop. Here, we present a 2.4 Å resolution crystal structure of an “extended” version of the Mps1 kinase domain, which shows an ordered activation loop. However, the other structural characteristics of an active kinase are not present. Our structure shows that the Mps1 activation loop can fit to the ATP binding pocket and interferes with ATP, but less so with inhibitors binding, partly explain the potency of various Mps1 inhibitors.     

Substrate Binding Mechanism of a Type I Extradiol Dioxygenase

A meta-cleavage pathway for the aerobic degradation of aromatic hydrocarbons is catalyzed by extradiol dioxygenases via a two-step mechanism: catechol substrate binding and dioxygen incorporation. The binding of substrate triggers the release of water, thereby opening a coordination site for molecular oxygen. The crystal structures of AkbC, a type I extradiol dioxygenase, and the enzyme substrate (3-methylcatechol) complex revealed the substrate binding process of extradiol dioxygenase. AkbC is composed of an N-domain and an active C-domain, which contains iron coordinated by a 2-His-1-carboxylate facial triad motif. The C-domain includes a β-hairpin structure and a C-terminal tail. In substrate-bound AkbC, 3-methylcatechol interacts with the iron via a single hydroxyl group, which represents an intermediate stage in the substrate binding process. Structure-based mutagenesis revealed that the C-terminal tail and β-hairpin form part of the substrate binding pocket that is responsible for substrate specificity by blocking substrate entry. Once a substrate enters the active site, these structural elements also play a role in the correct positioning of the substrate. Based on the results presented here, a putative substrate binding mechanism is proposed.     

New insights into the structural and functional involvement of the gate loop in AcrB export activity

AcrB is a major multidrug exporter in Escherichia coli and other Gram-negative bacteria. Its gate loop, located between the proximal and the distal pockets, have been reported to play important role in the export of many antibiotics. This loop location, rigidity and interactions with substrates have led recent reports to suggest that AcrB export mechanism operates in a sequential manner. First the substrate binds the proximal pocket in the access monomer, then it moves to bind the distal pocket in the binding monomer and subsequently it is extruded in the extrusion monomer. Recently, we have demonstrated that the gate loop is not required for the binding of Erythromycin but the integrity of this loop is important for an efficient export of this substrate. However, here we show that the antibiotic susceptibilities of the same AcrB gate loop mutants for Doxorubicin were unaffected, suggesting that this loop is not required for its export, and we demonstrate that this substrate may use principally the tunnel-1, located between transmembranes 8 and 9, more often than previously reported. To further explain our findings, here we address the gate loop mutations effects on AcrB solution energetics (fold, stability, molecular dynamics) and on the in vivo efflux of Erythromycin and Doxorubicin. Finally, we discuss the efflux and the discrepancy between the structural and the functional experiments for Erythromycin in these gate loop mutants.     

A model of pancreatic lipase and the orientation of enzymes at interfaces

A model for the interfacial orientation and the mode of action of lipase is proposed. Lipase is oriented so that its active site is near the oil-water boundary. This orientation is achieved by oil-enzyme bonding at the “hydrophobic head” of the enzyme, a region free of electric charges and relatively resistant to unfolding. The measured K_M is a complex constant including the dissociation constant of this oil-enzyme “complex”. The interfacial orientation of lipase is further aided by hydrophilic negative charges on the “back” of the enzyme and by a hydrophilic carbohydrate “tail”.It is suggested that similar hydrophobic heads and hydrophilic tails and asymmetric charge distributions establish the orientation of many enzymes which act at interfaces. Many phospholipases, for instance, appear to be charge-oriented, and the carbohydrate residues of ribonucleases and many other glycoproteins may be hydrophilic tails.Lipase is probably a serine enzyme with a catalytic center similar to that of chymotrypsin, but more hindered, perhaps owing to the presence of a leucine residue, and there is no binding of substrate lipid chains in the “active complex”. The substrate molecule is fixated on the enzyme in a two-dimensional orientation, because its leaving alkoxy group must be received by the serine hydroxyl hydrogen which is directed towards the imidazol ring of the reactive histidine through a hydrogen bond. The high turnover rate of lipolysis, 5 × 10⁵/min, exceptional even for an enzyme, results from the extremely high substrate concentration near the active site, and from an almost complete extrusion of water because of the hydrophobicity of both the active site and the substrate. In addition, both substrate and enzyme, because of their polarity, are already so favorably positioned at the interface that the formation of the “active complex” requires only a proper two-dimensional alignment, perhaps with partial extraction of the substrate molecule from the lipid phase.     

New insight into the architecture of oxy‐anion pocket in unliganded conformation of GAT domains: A MD‐simulation study

Human Guanine Monophosphate Synthetase (hGMPS) converts XMP to GMP, and acts as a bifunctional enzyme with N‐terminal “glutaminase” (GAT) and C‐terminal “synthetase” domain. The enzyme is identified as a potential target for anti‐cancer and immunosuppressive therapies. GAT domain of enzyme plays central role in metabolism, and contains conserved catalytic residues Cys104, His190, and Glu192. MD simulation studies on GAT domain suggest that position of oxyanion in unliganded conformation is occupied by one conserved water molecule (W1), which also stabilizes that pocket. This position is occupied by a negatively charged atom of the substrate or ligand in ligand bound crystal structures. In fact, MD simulation study of Ser75 to Val indicates that W1 conserved water molecule is stabilized by Ser75, while Thr152, and His190 also act as anchor residues to maintain appropriate architecture of oxyanion pocket through water mediated H‐bond interactions. Possibly, four conserved water molecules stabilize oxyanion hole in unliganded state, but they vacate these positions when the enzyme (hGMPS)‐substrate complex is formed. Thus this study not only reveals functionally important role of conserved water molecules in GAT domain, but also highlights essential role of other non‐catalytic residues such as Ser75 and Thr152 in this enzymatic domain. The results from this computational study could be of interest to experimental community and provide a testable hypothesis for experimental validation. Conserved sites of water molecules near and at oxyanion hole highlight structural importance of water molecules and suggest a rethink of the conventional definition of chemical geometry of inhibitor binding site. Proteins 2016; 84:360–373. © 2016 Wiley Periodicals, Inc.     

Structural Basis for Substrate Recognition and Specificity in Aklavinone-11-Hydroxylase from Rhodomycin Biosynthesis

In the biosynthesis of several anthracyclines, aromatic polyketides produced by many Streptomyces species, the aglycone core is modified by a specific flavin adenine dinucleotide (FAD)- and NAD(P)H-dependent aklavinone-11-hydroxylase. Here, we report the crystal structure of a ternary complex of this enzyme from Streptomyces purpurascens, RdmE, with FAD and the substrate aklavinone. The enzyme is built up of three domains, a FAD-binding domain, a domain involved in substrate binding, and a C-terminal thioredoxin-like domain of unknown function. RdmE exhibits structural similarity to aromatic hydroxylases from the p-hydroxybenzoate hydroxylase family, but unlike most other related enzymes, RdmE is a monomer. The substrate is bound in a hydrophobic pocket in the interior of the enzyme, and access to this pocket is provided through a different route than for the isoalloxazine ring of FAD—the backside of the ligand binding cleft. The architecture of the substrate binding pocket and the observed enzyme-aklavinone interactions provide a structural explanation for the specificity of the enzyme for non-glycosylated substrates with C9-R stereochemistry. The isoalloxazine ring of the flavin cofactor is bound in the “out” conformation but can be modeled in the “in” conformation without invoking large conformational changes of the enzyme. This model places the flavin ring in a position suitable for catalysis, almost perpendicular to the tetracyclic ring system of the substrate and with a distance of the C4a carbon atom of the isoalloxazine ring to the C-11 carbon atom of the substrate of 4.8 Å. The structure suggested that a Tyr224-Arg373 pair might be involved in proton abstraction at the C-6 hydroxyl group, thereby increasing the nucleophilicity of the aromatic ring system and facilitating electrophilic attack by the perhydroxy-flavin intermediate. Replacement of Tyr224 by phenylalanine results in inactive enzyme, whereas mutants at position Arg373 retain catalytic activity close to wild-type level. These data establish an essential role of residue Tyr224 in catalysis, possibly in aligning the substrate in a position suitable for catalysis.     

Controlling electron transfer in Acyl-CoA oxidases and dehydrogenases: a structural view

Plants produce a unique peroxisomal short chain-specific acyl-CoA oxidase (ACX4) for beta-oxidation of lipids. The short chain-specific oxidase has little resemblance to other peroxisomal acyl-CoA oxidases but has an approximately 30% sequence identity to mitochondrial acyl-CoA dehydrogenases. Two biochemical features have been linked to structural properties by comparing the structures of short chain-specific Arabidopsis thaliana ACX4 with and without a substrate analogue bound in the active site to known acyl-CoA oxidases and dehydrogenase structures: (i) a solvent-accessible acyl binding pocket is not required for oxygen reactivity, and (ii) the oligomeric state plays a role in substrate pocket architecture but is not linked to oxygen reactivity. The structures indicate that the acyl-CoA oxidases may encapsulate the electrons for transfer to molecular oxygen by blocking the dehydrogenase substrate interaction site with structural extensions. A small binding pocket observed adjoining the flavin adenine dinucleotide N5 and C4a atoms could increase the number of productive encounters between flavin adenine dinucleotide and O2.