Staphylococcal peroxidase inhibitor (SPIN): Residue-level investigation of the helical bundle domain

Myeloperoxidase is a critical component of the antibacterial arsenal of neutrophils, whereby it consumes H2O2 as an oxidant to convert halogen and pseudohalogen anions into cytotoxic hypohalous acids. Following phagocytosis by neutrophils, the human pathogen Staphylococcus aureus secretes a potent myeloperoxidase inhibitory protein, called SPIN, as part of its immune evasion repertoire. The matured S. aureus SPIN polypeptide consists of only 73 residues yet contains two functional domains: whereas the 60 residue C-terminal helical bundle domain is responsible for MPO binding, the 13 residue N-terminal domain is required to inhibit MPO. Previous studies have informed understanding of the SPIN N-terminal domain, but comparatively little is known about the helical domain insofar as the contribution of individual residues is concerned. To address this limitation, we carried out a residue-level structure/function investigation on the helical bundle domain of S. aureus SPIN. Using sequence conservation and existing structures of SPIN bound to human MPO as a guide, we selected residues L49, E50, H51, E52, Y55, and Y75 for interrogation by site-directed mutagenesis. We found that loss of L49 or E52 reduced SPIN activity by roughly an order of magnitude, but that loss of Y55 or H51 caused progressively greater loss of inhibitory potency. Direct binding studies by SPR showed that loss of inhibitory potency in these SPIN mutants resulted from a diminished initial interaction between the inhibitor and MPO. Together, our studies provide new insights into the structure/function relationships of SPIN and identify positions Y55 and H51 as critical determinants of SPIN function.

AI is a viable alternative to high throughput screening: a 318-target study

High throughput screening (HTS) is routinely used to identify bioactive small molecules. This requires physical compounds, which limits coverage of accessible chemical space. Computational approaches combined with vast on-demand chemical libraries can access far greater chemical space, provided that the predictive accuracy is sufficient to identify useful molecules. Through the largest and most diverse virtual HTS campaign reported to date, comprising 318 individual projects, we demonstrate that our AtomNet® convolutional neural network successfully finds novel hits across every major therapeutic area and protein class. We address historical limitations of computational screening by demonstrating success for target proteins without known binders, high-quality X-ray crystal structures, or manual cherry-picking of compounds. We show that the molecules selected by the AtomNet® model are novel drug-like scaffolds rather than minor modifications to known bioactive compounds. Our empirical results suggest that computational methods can substantially replace HTS as the first step of small-molecule drug discovery.

Inhibition of the C1s Protease and the Classical Complement Pathway by 6-(4-Phenylpiperazin-1-yl)Pyridine-3-Carboximidamide and Chemical Analogs

The classical pathway (CP) is a potent mechanism for initiating complement activity and is a driver of pathology in many complement-mediated diseases. The CP is initiated via activation of complement component C1, which consists of the pattern recognition molecule C1q bound to a tetrameric assembly of proteases C1r and C1s. Enzymatically active C1s provides the catalytic basis for cleavage of the downstream CP components, C4 and C2, and is therefore an attractive target for therapeutic intervention in CP-driven diseases. Although an anti-C1s mAb has been Food and Drug Administration approved, identifying small-molecule C1s inhibitors remains a priority. In this study, we describe 6-(4-phenylpiperazin-1-yl)pyridine-3-carboximidamide (A1) as a selective, competitive inhibitor of C1s. A1 was identified through a virtual screen for small molecules that interact with the C1s substrate recognition site. Subsequent functional studies revealed that A1 dose-dependently inhibits CP activation by heparin-induced immune complexes, CP-driven lysis of Ab-sensitized sheep erythrocytes, CP activation in a pathway-specific ELISA, and cleavage of C2 by C1s. Biochemical experiments demonstrated that A1 binds directly to C1s with a Kd of ∼9.8 μM and competitively inhibits its activity with an inhibition constant (Ki) of ∼5.8 μM. A 1.8-Å-resolution crystal structure revealed the physical basis for C1s inhibition by A1 and provided information on the structure-activity relationship of the A1 scaffold, which was supported by evaluating a panel of A1 analogs. Taken together, our work identifies A1 as a new class of small-molecule C1s inhibitor and lays the foundation for development of increasingly potent and selective A1 analogs for both research and therapeutic purposes.

Simultaneous inhibition of two neutrophil serine proteases by the S. aureus innate immune evasion protein EapH2

Extracellular adherence protein domain (EAP) proteins are high-affinity, selective inhibitors of neutrophil serine proteases (NSP), including cathepsin-G (CG) and neutrophil elastase (NE). Most Staphylococcus aureus isolates encode for two EAPs, EapH1 and EapH2, that contain a single functional domain and share 43% identity with one another. Although structure/function investigations from our group have shown that EapH1 uses a globally similar binding mode to inhibit CG and NE, NSP inhibition by EapH2 is incompletely understood due to a lack of NSP/EapH2 cocrystal structures. To address this limitation, we further studied NSP inhibition by EapH2 in comparison with EapH1. Like its effects on NE, we found that EapH2 is a reversible, time-dependent, and low nanomolar affinity inhibitor of CG. We characterized an EapH2 mutant which suggested that the CG binding mode of EapH2 is comparable to EapH1. To test this directly, we used NMR chemical shift perturbation to study EapH1 and EapH2 binding to CG and NE in solution. Although we found that overlapping regions of EapH1 and EapH2 were involved in CG binding, we found that altogether distinct regions of EapH1 and EapH2 experienced changes upon binding to NE. An important implication of this observation is that EapH2 might be capable of binding and inhibiting CG and NE simultaneously. We confirmed this unexpected feature by solving crystal structures of the CG/EapH2/NE complex and demonstrating their functional relevance through enzyme inhibition assays. Together, our work defines a new mechanism of simultaneous inhibition of two serine proteases by a single EAP protein.

Characterization of two distinct neutrophil serine protease-binding modes within a Staphylococcus aureus innate immune evasion protein family

Extracellular adherence protein domain (EAPs) proteins are a class of innate immune evasion proteins secreted by the human pathogen Staphylococcus aureus. EAPs are potent and selective inhibitors of cathepsin-G (CG) and neutrophil elastase (NE), which are the two most abundant neutrophil serine proteases (NSPs). Previous work from our group has shown that the prototypical EAP, EapH1, relies on plasticity within a single inhibitory site to block the activities of CG and NE. However, whether other EAPs follow similar structure-function relationships is unclear. To address this question, we studied the inhibitory properties of the first (Eap1) and second (Eap2) domains of the modular extracellular adherence protein of S. aureus and determined their structures when bound to CG and NE, respectively. We observed that both Eap1 and Eap2 displayed time-dependent inhibition of CG (on the order of 10^-9 M) and of NE (on the order of 10^-10 M). We also found that whereas the structures of Eap1 and Eap2 bound to CG showed an overall inhibitory mode like that seen previously for EapH1, the structures of Eap1 and Eap2 bound to NE revealed a new inhibitory mode involving a distal region of the EAP domain. Using site-directed mutagenesis of Eap1 and Eap2, along with enzyme assays, we confirmed the roles of interfacial residues in NSP inhibition. Taken together, our work demonstrates that EAPs can form structurally divergent complexes with two closely related serine proteases and further suggests that certain EAPs may be capable of inhibiting two NSPs simultaneously.

Local structural plasticity of the Staphylococcus aureus evasion protein EapH1 enables engagement with multiple neutrophil serine proteases

Members of the EAP family of Staphylococcus aureus immune evasion proteins potently inhibit the neutrophil serine proteases (NSPs) neutrophil elastase, cathepsin-G, and proteinase-3. Previously, we determined a 1.8 Å resolution crystal structure of the EAP family member EapH1 bound to neutrophil elastase. This structure revealed that EapH1 blocks access to the enzyme's active site by forming a noncovalent complex with this host protease. To determine how EapH1 inhibits other NSPs, we studied here the effects of EapH1 on cathepsin-G. We found that EapH1 inhibits cathepsin-G with a K_i of 9.8 ± 4.7 nm Although this K_i value is ∼466-fold weaker than the K_i for EapH1 inhibition of neutrophil elastase, the time dependence of inhibition was maintained. To define the physical basis for EapH1's inhibition of cathepsin-G, we crystallized EapH1 bound to this protease, solved the structure at 1.6 Å resolution, and refined the model to R_work and R_free values of 17.4% and 20.9%, respectively. This structure revealed a protease-binding mode for EapH1 with cathepsin-G that was globally similar to that seen in the previously determined EapH1-neutrophil elastase structure. The nature of the intermolecular interactions formed by EapH1 with cathepsin-G differed considerably from that with neutrophil elastase, however, with far greater contributions from the inhibitor backbone in the cathepsin-G-bound form. Together, these results reveal that EapH1's ability to form high-affinity interactions with multiple NSP targets is due to its remarkable level of local structural plasticity.

Staphylococcus aureus evasion proteins EapH1 and EapH2: Residue-level investigation of an alternative binding motif for human neutrophil elastase

The Staphylococcus aureusExtracellular Adherence Protein (Eap) and its homologs, EapH1 and EapH2, are a family of secreted proteins that potently inhibit the neutrophil serine proteases Neutrophil Elastase (hNE), Cathepsin G, and Proteinase 3. Similarly to EapH1, inhibition of hNE by EapH2 is characterized by a rapid association rate (2.9 × 10⁵ M^-1s^-1) coupled with a very slow dissociation rate (5.9 × 10^-4 s^-1), yielding an apparent inhibition constant of 2.11 nM. As with EapH1, inhibition of hNE by EapH2 is also time-dependent in character. A phenylalanine in EapH2 replaces the leucine in EapH1 that sits over the hNE catalytic serine and creates a potential steric clash. Indeed, the EapH1 L59F mutant is severely decreased in its ability to inhibit hNE (~9500-fold). When compared to the EapH1:hNE co-crystal structure, a model of the EapH2:hNE complex predicts an alternative binding motif comprised of EapH2 residues 120-127. These putative interfacing residues were individually mutated and kinetically interrogated. The EapH2 N127A mutant resulted in the largest decrease in hNE inhibition (~200-fold) and loss of the time-dependent characteristic. Surprisingly, the time-dependent characteristic was still abolished in the EapH2 T125A mutant, even though it was less perturbed in hNE inhibition (~25-fold). T125 forms an intra-molecular hydrogen bond to the carbonyl oxygen of N127 in the EapH2 crystal structure. Given these observations, we conclude (i) that EapH2 has an altogether distinct hNE binding motif than EapH1, (ii) that N127 is the main functional determinant in EapH2, and (iii) that T125 serves an ancillary role aiding in the optimal orientation of N127.

Investigation of Human Neutrophil Elastase Inhibition by Staphylococcus aureus EapH1: The Key Role Played by Arginine 89

Staphylococcus aureus secretes a family of potent, noncovalent inhibitory proteins that selectively target the neutrophil serine proteases neutrophil elastase, cathepsin-G, and proteinase-3. A majority of our understanding of these so-called EAP domain proteins has come from structure/function studies on EapH1 and its effects on human neutrophil elastase (hNE). Inspection of the EapH1/hNE cocrystal structure suggested that EapH1 residues R89, E94, and K95 are positioned near the EapH1/hNE interface and might contribute to the potent inhibition of hNE by EapH1. In this study, we used site directed mutagenesis, kinetic evaluation, and surface plasmon resonance to probe the individual contributions of R89, E94, and K95 to EapH1 function. We found that the wild-type EapH1/hNE complex is characterized by a fast association rate (2.0 × 10⁶ M^-1 s^-1) and a very slow dissociation rate (4.3 × 10^-5 s^-1), yielding an apparent inhibition constant of 21 pM. The slow dissociation rate of EapH1 from hNE resulted in a time-dependent inhibition pattern. Although conservative mutants E94Q and K95M, as well as the E94Q/K95M double mutant, had on- and off-rates comparable to wild-type EapH1, mutation of R89 to methionine resulted in a 15,000-fold decrease in inhibition (321 nM) and loss of the time-dependent inhibition characteristic. The double mutants R89M/E94Q and R89M/K95M, as well as the triple mutant R89M/E94Q/K95M were similarly perturbed. Mutation of R89 to lysine restored a portion of the inhibition of hNE (27 nM). Given these observations, we conclude that R89 is a primary contributor to EapH1 function vis-à-vis time-dependent inhibition of hNE.

Functional Analysis of the Bacteriophage T4 Rad50 Homolog (gp46) Coiled-coil Domain

Rad50 and Mre11 form a complex involved in the detection and processing of DNA double strand breaks. Rad50 contains an anti-parallel coiled-coil with two absolutely conserved cysteine residues at its apex. These cysteine residues serve as a dimerization domain and bind a Zn(2+) cation in a tetrathiolate coordination complex known as the zinc-hook. Mutation of the zinc-hook in bacteriophage T4 is lethal, indicating the ability to bind Zn(2+) is critical for the functioning of the MR complex. In vitro, we found that complex formation between Rad50 and a peptide corresponding to the C-terminal domain of Mre11 enhances the ATPase activity of Rad50, supporting the hypothesis that the coiled-coil is a major conduit for communication between Mre11 and Rad50. We constructed mutations to perturb this domain in the bacteriophage T4 Rad50 homolog. Deletion of the Rad50 coiled-coil and zinc-hook eliminates Mre11 binding and ATPase activation but does not affect its basal activity. Mutation of the zinc-hook or disruption of the coiled-coil does not affect Mre11 or DNA binding, but their activation of Rad50 ATPase activity is abolished. Although these mutants excise a single nucleotide at a normal rate, they lack processivity and have reduced repetitive exonuclease rates. Restricting the mobility of the coiled-coil eliminates ATPase activation and repetitive exonuclease activity, but the ability to support single nucleotide excision is retained. These results suggest that the coiled-coiled domain adopts at least two conformations throughout the ATPase/nuclease cycle, with one conformation supporting enhanced ATPase activity and processivity and the other supporting nucleotide excision.

Catalytic mechanism of bacteriophage T4 Rad50 ATP hydrolysis

Spontaneous double-strand breaks (DSBs) are one of the most deleterious forms of DNA damage, and their improper repair can lead to cellular dysfunction. The Mre11 and Rad50 proteins, a nuclease and an ATPase, respectively, form a well-conserved complex that is involved in the initial processing of DSBs. Here we examine the kinetic and catalytic mechanism of ATP hydrolysis by T4 Rad50 (gp46) in the presence and absence of Mre11 (gp47) and DNA. Single-turnover and pre-steady state kinetics on the wild-type protein indicate that the rate-limiting step for Rad50, the MR complex, and the MR-DNA complex is either chemistry or a conformational change prior to catalysis. Pre-steady state product release kinetics, coupled with viscosity steady state kinetics, also supports that the binding of DNA to the MR complex does not alter the rate-limiting step. The lack of a positive deuterium solvent isotope effect for the wild type and several active site mutants, combined with pH-rate profiles, implies that chemistry is rate-limiting and the ATPase mechanism proceeds via an asymmetric, dissociative-like transition state. Mutation of the Walker A/B and H-loop residues also affects the allosteric communication between Rad50 active sites, suggesting possible routes for cooperativity between the ATP active sites.

Molecular docking and NMR binding studies to identify novel inhibitors of human phosphomevalonate kinase

Phosphomevalonate kinase (PMK) phosphorylates mevalonate-5-phosphate (M5P) in the mevalonate pathway, which is the sole source of isoprenoids and steroids in humans. We have identified new PMK inhibitors with virtual screening, using autodock. Promising hits were verified and their affinity measured using NMR-based (1)H-(15)N heteronuclear single quantum coherence (HSQC) chemical shift perturbation and fluorescence titrations. Chemical shift changes were monitored, plotted, and fitted to obtain dissociation constants (K(d)). Tight binding compounds with K(d)'s ranging from 6-60 μM were identified. These compounds tended to have significant polarity and negative charge, similar to the natural substrates (M5P and ATP). HSQC cross peak changes suggest that binding induces a global conformational change, such as domain closure. Compounds identified in this study serve as chemical genetic probes of human PMK, to explore pharmacology of the mevalonate pathway, as well as starting points for further drug development.

Disruption of the bacteriophage T4 Mre11 dimer interface reveals a two-state mechanism for exonuclease activity

The Mre11-Rad50 (MR) complex is a central player in DNA repair and is implicated in the processing of DNA ends caused by double strand breaks. Recent crystal structures of the MR complex suggest that several conformational rearrangements occur during its ATP hydrolysis cycle. A comparison of the Mre11 dimer interface from these structures suggests that the interface is dynamic in nature and may adopt several different arrangements. To probe the functional significance of the Mre11 dimer interface, we have generated and characterized a dimer disruption Mre11 mutant (L101D-Mre11). Although L101D-Mre11 binds to Rad50 and dsDNA with affinity comparable with the wild-type enzyme, it does not activate the ATP hydrolysis activity of Rad50, suggesting that the allosteric communication between Mre11 and Rad50 has been interrupted. Additionally, the dsDNA exonuclease activity of the L101D-MR complex has been reduced by 10-fold under conditions where processive exonuclease activity is required. However, we unexpectedly found that under steady state conditions, the nuclease activity of the L101D-MR complex is significantly greater than that of the wild-type complex. Based on steady state and single-turnover nuclease assays, we have assigned the rate-determining step of the steady state nuclease reaction to be the productive assembly of the complex at the dsDNA end. Together, our data suggest that the Mre11 dimer interface adopts at least two different states during the exonuclease reaction.

Functional evaluation of bacteriophage T4 Rad50 signature motif residues

The repair of DNA double-strand breaks (DSBs) is essential to maintaining the integrity of the genome, and organisms have evolved a conserved mechanism to facilitate their repair. In eukaryotes, archaea, and some bacteriophage, a complex made up of Mre11 and Rad50 (MR complex), which are a nuclease and ATPase, respectively, is involved in the initial processing of DSBs. Rad50 is a member of the ATP Binding Cassette (ABC) protein superfamily, the members of which contain an important Signature motif that acts in trans to complete the dimeric ATP binding site. To explore the functional relevance of this motif, four of its five residues were mutated in bacteriophage T4 Rad50, and their respective ATPase and nuclease activities were evaluated. The mutations reveal the functional roles of the Signature motif in ATP binding, hydrolysis, and cooperativity. In several mutants, the degree of DNA activation of ATP hydrolysis activity is reduced, indicating that the Signature motif is involved in allosteric signal transmission between the DNA and ATP binding sites of the MR complex. ATP hydrolysis is not required for nuclease activity when the probe is near the beginning of the DNA substrate; however, when an internal probe is used, decreases in ATPase activity have substantial effects on nuclease activity, suggesting that ATP hydrolysis is involved in translocation of the complex. Unexpectedly, the ATP hydrolysis and nuclease activities are not directly correlated with each other, and each mutation appears to differentially affect the exonuclease activity of Mre11.

Crystal structures of Staphylococcus epidermidis mevalonate diphosphate decarboxylase bound to inhibitory analogs reveal new insight into substrate binding and catalysis

The polyisoprenoid compound undecaprenyl phosphate is required for biosynthesis of cell wall peptidoglycans in gram-positive bacteria, including pathogenic Enterococcus, Streptococcus, and Staphylococcus spp. In these organisms, the mevalonate pathway is used to produce the precursor isoprenoid, isopentenyl 5-diphosphate. Mevalonate diphosphate decarboxylase (MDD) catalyzes formation of isopentenyl 5-diphosphate in an ATP-dependent irreversible reaction and is therefore an attractive target for inhibitor development that could lead to new antimicrobial agents. To facilitate exploration of this possibility, we report the crystal structure of Staphylococcus epidermidis MDD (1.85 Å resolution) and, to the best of our knowledge, the first structures of liganded MDD. These structures include MDD bound to the mevalonate 5-diphosphate analogs diphosphoglycolyl proline (2.05 Å resolution) and 6-fluoromevalonate diphosphate (FMVAPP; 2.2 Å resolution). Comparison of these structures provides a physical basis for the significant differences in K(i) values observed for these inhibitors. Inspection of enzyme/inhibitor structures identified the side chain of invariant Ser(192) as making potential contributions to catalysis. Significantly, Ser → Ala substitution of this side chain decreases k(cat) by ∼10(3)-fold, even though binding interactions between FMVAPP and this mutant are similar to those observed with wild type MDD, as judged by the 2.1 Å cocrystal structure of S192A with FMVAPP. Comparison of microbial MDD structures with those of mammalian counterparts reveals potential targets at the active site periphery that may be exploited to selectively target the microbial enzymes. These studies provide a structural basis for previous observations regarding the MDD mechanism and inform future work toward rational inhibitor design.

Biochemical characterization of bacteriophage T4 Mre11-Rad50 complex

The Mre11-Rad50 complex (MR) from bacteriophage T4 (gp46/47) is involved in the processing of DNA double-strand breaks. Here, we describe the activities of the T4 MR complex and its modulation by proteins involved in homologous recombination. T4 Mre11 is a Rad50- and Mn(2+)-dependent dsDNA exonuclease and ssDNA endonuclease. ATP hydrolysis is required for the removal of multiple nucleotides via dsDNA exonuclease activity but not for the removal of the first nucleotide or for ssDNA endonuclease activity, indicating ATP hydrolysis is only required for repetitive nucleotide removal. By itself, Rad50 is a relatively inefficient ATPase, but the presence of Mre11 and dsDNA increases ATP hydrolysis by 20-fold. The ATP hydrolysis reaction exhibits positive cooperativity with Hill coefficients ranging from 1.4 for Rad50 alone to 2.4 for the Rad50-Mre11-DNA complex. Kinetic assays suggest that approximately four nucleotides are removed per ATP hydrolyzed. Directionality assays indicate that the prevailing activity is a 3' to 5' dsDNA exonuclease, which is incompatible with the proposed role of MR in the production of 3' ssDNA ends. Interestingly, we found that in the presence of a recombination mediator protein (UvsY) and ssDNA-binding protein (gp32), Mre11 is capable of using Mg(2+) as a cofactor for its nuclease activity. Additionally, the Mg(2+)-dependent nuclease activity, activated by UvsY and gp32, results in the formation of endonuclease reaction products. These results suggest that gp32 and UvsY may alter divalent cation preference and facilitate the formation of a 3' ssDNA overhang, which is a necessary intermediate for recombination-mediated double-strand break repair.

NMR dynamics investigation of ligand-induced changes of main and side-chain arginine N-H's in human phosphomevalonate kinase

Phosphomevalonate kinase (PMK) catalyzes phosphoryl transfer from adenosine triphosphate (ATP) to mevalonate 5-phosphate (M5P) on the pathway for synthesizing cholesterol and other isoprenoids. To permit this reaction, its substrates must be brought proximal, which would result in a significant and repulsive buildup of negative charge. To facilitate this difficult task, PMK contains 17 arginines and eight lysines. However, the way in which this charge neutralization and binding is achieved, from a structural and dynamics perspective, is not known. More broadly, the role of arginine side-chain dynamics in binding of charged substrates has not been experimentally defined for any protein to date. Herein we report a characterization of changes to the dynamical state of the arginine side chains in PMK due to binding of its highly charged substrates, ATP and M5P. These studies were facilitated by the use of arginine-selective labeling to eliminate spectral overlap. Model-free analysis indicated that while substrate binding has little effect on the arginine backbone dynamics, binding of either substrate leads to significant rigidification of the arginine side chains throughout the protein, even those that are >8 A from the binding site. Such a global rigidification of arginine side chains is unprecedented and suggests that there are long-range electrostatic interactions of sufficient strength to restrict the motion of arginine side chains on the picosecond-to-nanosecond time scale. It will be interesting to see whether such effects are general for arginine residues in proteins that bind highly charged substrates, once additional studies of arginine side-chain dynamics are reported.

Substrate induced structural and dynamics changes in human phosphomevalonate kinase and implications for mechanism

Phosphomevalonate kinase (PMK) catalyzes an essential step in the mevalonate pathway, which is the only pathway for synthesis of isoprenoids and steroids in humans. PMK catalyzes transfer of the gamma-phosphate of ATP to mevalonate 5-phosphate (M5P) to form mevalonate 5-diphosphate. Bringing these phosphate groups in proximity to react is especially challenging, given the high negative charge density on the four phosphate groups in the active site. As such, conformational and dynamics changes needed to form the Michaelis complex are of mechanistic interest. Herein, we report the characterization of substrate induced changes (Mg-ADP, M5P, and the ternary complex) in PMK using NMR-based dynamics and chemical shift perturbation measurements. Mg-ADP and M5P K(d)'s were 6-60 microM in all complexes, consistent with there being little binding synergy. Binding of M5P causes the PMK structure to compress (tau(c) = 13.5 nsec), whereas subsequent binding of Mg-ADP opens the structure up (tau(c) = 15.6 nsec). The overall complex seems to stay very rigid on the psec-nsec timescale with an average NMR order parameter of S(2) approximately 0.88. Data are consistent with addition of M5P causing movement around a hinge region to permit domain closure, which would bring the M5P domain close to ATP to permit catalysis. Dynamics data identify potential hinge residues as H55 and R93, based on their low order parameters and their location in extended regions that connect the M5P and ATP domains in the PMK homology model. Likewise, D163 may be a hinge residue for the lid region that is homologous to the adenylate kinase lid, covering the "Walker-A" catalytic loop. Binding of ATP or ADP appears to cause similar conformational changes; however, these observations do not indicate an obvious role for gamma-phosphate binding interactions. Indeed, the role of gamma-phosphate interactions may be more subtle than suggested by ATP/ADP comparisons, because the conservative O to NH substitution in the beta-gamma bridge of ATP causes a dramatic decrease in affinity and induces few chemical shift perturbations. In terms of positioning of catalytic residues, binding of M5P induces a rigidification of Gly21 (adjacent to the catalytically important Lys22), although exchange broadening in the ternary complex suggests some motion on a slower timescale does still occur. Finally, the first nine residues of the N-terminus are highly disordered, suggesting that they may be part of a cleavable signal or regulatory peptide sequence.

Human mevalonate diphosphate decarboxylase: characterization, investigation of the mevalonate diphosphate binding site, and crystal structure

Expression in Escherichia coli of his-tagged human mevalonate diphosphate decarboxylase (hMDD) has expedited enzyme isolation, characterization, functional investigation of the mevalonate diphosphate binding site, and crystal structure determination (2.4A resolution). hMDD exhibits V(max)=6.1+/-0.5 U/mg; K(m) for ATP is 0.69+/-0.07 mM and K(m) for (R,S) mevalonate diphosphate is 28.9+/-3.3 microM. Conserved polar residues predicted to be in the hMDD active site were mutated to test functional importance. R161Q exhibits a approximately 1000-fold diminution in specific activity, while binding the fluorescent substrate analog, TNP-ATP, comparably to wild-type enzyme. Diphosphoglycolyl proline (K(i)=2.3+/-0.3 uM) and 6-fluoromevalonate 5-diphosphate (K(i)=62+/-5 nM) are competitive inhibitors with respect to mevalonate diphosphate. N17A exhibits a V(max)=0.25+/-0.0 2U/mg and a 15-fold inflation in K(m) for mevalonate diphosphate. N17A's K(i) values for diphosphoglycolyl proline and fluoromevalonate diphosphate are inflated (>70-fold and 40-fold, respectively) in comparison with wild-type enzyme. hMDD structure indicates the proximity (2.8A) between R161 and N17, which are located in an interior pocket of the active site cleft. The data suggest the functional importance of R161 and N17 in the binding and orientation of mevalonate diphosphate.

Biochemical and structural basis for feedback inhibition of mevalonate kinase and isoprenoid metabolism

Mevalonate kinase (MK), which catalyzes a key reaction in polyisoprenoid and sterol metabolism in many organisms, is subject to feedback regulation by farnesyl diphosphate and related compounds. The structures of human mevalonate kinase and a binary complex of the rat enzyme incubated with farnesyl thiodiphosphate (FSPP) are reported. Significant FSPP hydrolysis occurs under crystallization conditions; this results in detection of farnesyl thiophosphate (FSP) in the structure of the binary complex. Farnesyl thiodiphosphate competes with substrate ATP to produce feedback inhibition of mevalonate kinase. The binding sites for these metabolites overlap, with the phosphate of FSP nearly superimposed on ATP's beta-phosphate and FSP's polyisoprenoid chain overlapping ATP's adenosine moiety. Several hydrophobic amino acid side chains are positioned near the polyisoprenoid chain of FSP and their functional significance has been evaluated in mutagenesis experiments with human MK, which exhibits the highest reported sensitivity to feedback inhibition. Results suggest that single and double mutations at T104 and I196 produce a significant inflation of the K(i) for FSPP (approximately 40-fold for T104A/I196A). Such an effect persists when K(i) values are normalized for effects on the K(m) for ATP, suggesting that it may be possible to engineer MK proteins with altered sensitivity to feedback inhibition. Comparison of animal MK protein alignments and structures with those of a MK protein from Streptococcus pneumoniae indicates that sequence differences between N- and C-terminal domains correlate with differences in interdomain angles. Bacterial MK proteins exhibit more solvent exposure of feedback inhibitor binding sites and, consequently, weaker binding of these inhibitors.

Functional evaluation of conserved basic residues in human phosphomevalonate kinase

Phosphomevalonate kinase (PMK) catalyzes the cation-dependent reaction of mevalonate 5-phosphate with ATP to form mevalonate 5-diphosphate and ADP, a key step in the mevalonate pathway for isoprenoid/sterol biosynthesis. Animal PMK proteins belong to the nucleoside monophosphate (NMP) kinase family. For many NMP kinases, multiple basic residues contribute to the neutralization of the negatively charged pentacoordinate phosphate reaction intermediate. Loss of basicity can result in catalytically impaired enzymes. On the basis of this precedent, conserved basic residues of human PMK have been mutated, and purified forms of the mutated proteins have been kinetically and biophysically characterized. K48M and R73M mutants exhibit diminished Vmax values in both reaction directions (>1000-fold) with only slight Km perturbations (<10-fold). In both forward and reverse reactions, R110M exhibits a large (>10,000-fold) specific activity diminution. R111M exhibits substantially inflated Km values for mevalonate 5-phosphate and mevalonate 5-diphosphate (60- and 30-fold, respectively) as well as decreases [50-fold (forward) and 85-fold (reverse)] in Vmax. R84M also exhibits inflated Km values (50- and 33-fold for mevalonate 5-phosphate and mevalonate 5-diphosphate, respectively). The Ki values for R111M and R84M product inhibition by mevalonate 5-diphosphate are inflated by 45- and 63-fold; effects are comparable to the 30- and 38-fold inflations in Km for mevalonate 5-diphosphate. R141M exhibits little perturbation in Vmax [14-fold (forward) and 10-fold (reverse)] but has inflated Km values for ATP and ADP (48- and 136-fold, respectively). The Kd of ATP for R141M, determined by changes in tryptophan fluorescence, is inflated 27-fold compared to wt PMK. These data suggest that R110 is important to PMK catalysis, which is also influenced by K48 and R73. R111 and R84 contribute to binding of mevalonate 5-phosphate and R141 to binding of ATP.

Phosphomevalonate Kinase: Functional Investigation of the Recombinant Human Enzyme

Phosphomevalonate kinase (PMK) catalyzes a key step in isoprenoid/sterol biosynthesis, converting mevalonate 5-phosphate and ATP to mevalonate 5-diphosphate and ADP. To expedite functional and structural study of this enzyme, an expression plasmid encoding His-tagged human PMK has been constructed and recombinant enzyme isolated in an active, stable form. PMK catalyzes a reversible reaction; kinetic constants of human PMK have been determined for both forward (formation of mevalonate 5-diphosphate) and reverse (formation of mevalonate 5-phosphate) reactions. Animal and invertebrate PMKs are not orthologous to plant, fungal, or bacterial PMKs, limiting the information available from sequence alignment analysis. A homology model for the structure of human PMK has been generated. The model conforms to a nucleoside monophosphate kinase family fold. This result, together with sequence comparisons of animal and invertebrate PMKs, suggests an N-terminal basic residue rich sequence as a possible "Walker A" ATP binding motif. The functions of four basic (K17, R18, K19, K22) residues and one acidic (D23) residue in the conserved sequence have been tested by mutagenesis and characterization of isolated mutant proteins. Substrate K(m) values for K17M, R18Q, K19M, and D23N have been measured for forward and reverse reactions; in comparison with wild-type PMK values, only modest (<12-fold) changes are observed. In contrast, R18Q exhibits a V(max) decrease of 100/300-fold (forward/reverse reaction). K22M activity is too low for measurement at nonsaturating substrate concentration; specific activity is decreased by >10000-fold in both forward/reverse reactions, suggesting an active site location and an important role in phosphoryl transfer.

Functional dissection of a conserved motif within the pilus retraction protein PilT

PilT is a hexameric ATPase required for type IV pilus retraction in gram-negative bacteria. Retraction of type IV pili mediates intimate attachment to and signaling in host cells, surface motility, biofilm formation, natural transformation, and phage sensitivity. We investigated the in vivo and in vitro roles of each amino acid of the distinct, highly conserved C-terminal AIRNLIRE motif in PilT. Substitution of amino acids A288, I289, L292, and I293 as well as a double substitution of R290 and R294 abolished Pseudomonas aeruginosa PilT function in vivo, as measured by a loss of surface motility and phage sensitivity. When introduced into purified Aquifex aeolicus PilT, substitutions in the AIRNLIRE motif did not disrupt ATPase activity or oligomerization. In contrast, a K136Q substitution in the broadly conserved nucleotide binding motif prevented PilT function in vivo as well as in vitro. We propose that the AIRNLIRE motif forms an amphipathic alpha helix which transmits signals between a surface-exposed protein interaction site and the ATPase core of PilT, and we recognize a potential functional homology in other type II secretion ATPases.

The pilus-retraction protein PilT: ultrastructure of the biological assembly

PilT is a biological motor required for the retraction of bacterial type IV pili. Nesseria gonorrhoeae PilT has been purified and its ultrastructure has been examined by freeze-etch electron microscopy, revealing a 115 A outer diameter, 15-35 A inner diameter ring. Aquifex aeolicus PilT crystals were obtained in a primitive hexagonal space group (unit-cell parameters a = b = 107.3, c = 68.5 A) and diffract to a minimum Bragg spacing of 2.8 A when PilT is co-crystallized with adenine nucleotides. Initial phases to 3.5 A resolution have been determined by multiwavelength anomalous dispersion and density modification. Resulting electron-density maps show a hexameric A. aeolicus PilT ring 105 A wide by 55 A high, with an inner cavity that varies in shape and width from 20 to 40 A over the height of the complex. Both PilT ultrastructures are very similar to type II and type IV secretion ATPases in overall shape, size and assembly.

Aquifex aeolicus PilT, homologue of a surface motility protein, is a thermostable oligomeric NTPase

Bacterial surface motility works by retraction of surface-attached type IV pili. This retraction requires the PilT protein, a member of a large family of putative NTPases from type II and IV secretion systems. In this study, the PilT homologue from the thermophilic eubacterium Aquifex aeolicus was cloned, overexpressed, and purified. A. aeolicus PilT was shown to be a thermostable ATPase with a specific activity of 15.7 nmol of ATP hydrolyzed/min/mg of protein. This activity was abolished when a conserved lysine in the nucleotide-binding motif was altered. The substrate specificity was low; UTP, CTP, ATP, GTP, dATP, and dGTP served as substrates, UTP having the highest activity of these in vitro. Based on sedimentation equilibrium and size exclusion chromatography, PilT was identified as a approximately equal 5- to 6-subunit oligomer. Potential implications of the NTPase activity of PilT in pilus retraction are discussed.