Structural and dynamic requirements for optimal activity of the essential bacterial enzyme dihydrodipicolinate synthase

A metal ion-dependent conformational switch modulates activity of the Plasmodium M17 aminopeptidase

The metal-dependent M17 aminopeptidases are conserved throughout all kingdoms of life. This large enzyme family is characterized by a conserved binuclear metal center and a distinctive homohexameric arrangement. Recently, we showed that hexamer formation in Plasmodium M17 aminopeptidases was controlled by the metal ion environment, although the functional necessity for hexamer formation is still unclear. To further understand the mechanistic role of the hexameric assembly, here we undertook an investigation of the structure and dynamics of the M17 aminopeptidase from Plasmodium falciparum, PfA-M17. We describe a novel structure of PfA-M17, which shows that the active sites of each trimer are linked by a dynamic loop, and loop movement is coupled with a drastic rearrangement of the binuclear metal center and substrate-binding pocket, rendering the protein inactive. Molecular dynamics simulations and biochemical analyses of PfA-M17 variants demonstrated that this rearrangement is inherent to PfA-M17, and that the transition between the active and inactive states is metal dependent and part of a dynamic regulatory mechanism. Key to the mechanism is a remodeling of the binuclear metal center, which occurs in response to a signal from the neighboring active site and serves to moderate the rate of proteolysis under different environmental conditions. In conclusion, this work identifies a precise mechanism by which oligomerization contributes to PfA-M17 function. Furthermore, it describes a novel role for metal cofactors in the regulation of enzymes, with implications for the wide range of metalloenzymes that operate via a two-metal ion catalytic center, including DNA processing enzymes and metalloproteases.

Stochastic protein multimerization, activity, and fitness

Many proteins assemble into homomultimeric structures, with a number of subunits that can vary substantially among phylogenetic lineages. As protein-protein interactions require productive encounters among subunits, such variation might partially be explained by variation in cellular protein abundance. Protein abundance in turn depends on the intrinsic rates of production and decay of mRNA and protein molecules, as well as rates of cell growth and division. Using a stochastic framework for prediction of the multimeric state of a protein as a function of these processes and the free energy associated with interface-interface binding, we demonstrate agreement with a wide class of proteins using E. coli proteome data. As such, this platform, which links protein quaternary structure with biochemical rates governing gene expression, protein association and dissociation, and cell growth and division, can be extended to evolutionary models for the emergence and diversification of multimers. While it is tempting to think of multimerization as adaptive, the diversity of multimeric states raises the question of its functional role and impact on fitness. As a force driving selection, we consider the possible increase in enzymatic activity of proteins arising strictly as a consequence of interface-interface binding-namely, enhanced stability to degradation, substrate binding affinity, or catalytic rate of multimers with respect to monomers without invoking further conformational changes, as in allostery. For fixed cost of protein production, we find a benefit conferred by multimers that is dependent on context and can therefore become different in diverging lineages.

Substrate Locking Promotes Dimer-Dimer Docking of an Enzyme Antibiotic Target

Protein dynamics manifested through structural flexibility play a central role in the function of biological molecules. Here we explore the substrate-mediated change in protein flexibility of an antibiotic target enzyme, Clostridium botulinum dihydrodipicolinate synthase. We demonstrate that the substrate, pyruvate, stabilizes the more active dimer-of-dimers or tetrameric form. Surprisingly, there is little difference between the crystal structures of apo and substrate-bound enzyme, suggesting protein dynamics may be important. Neutron and small-angle X-ray scattering experiments were used to probe substrate-induced dynamics on the sub-second timescale, but no significant changes were observed. We therefore developed a simple technique, coined protein dynamics-mass spectrometry (ProD-MS), which enables measurement of time-dependent alkylation of cysteine residues. ProD-MS together with X-ray crystallography and analytical ultracentrifugation analyses indicates that pyruvate locks the conformation of the dimer that promotes docking to the more active tetrameric form, offering insight into ligand-mediated stabilization of multimeric enzymes.

Comparison of untagged and his-tagged dihydrodipicolinate synthase from the enteric pathogen Vibrio cholerae

Given the emergence of multi drug resistant Vibrio cholerae strains, there is an urgent need to characterize new anti-cholera targets. One such target is the enzyme dihydrodipicolinate synthase (DHDPS; EC 4.3.3.7), which catalyzes the first committed step in the diaminopimelate pathway. This pathway is responsible for the production of two key metabolites in bacteria and plants, namely meso-2,6-diaminopimelate and L-lysine. Here, we report the cloning, expression and purification of untagged and His-tagged recombinant DHDPS from V. cholerae (Vc-DHDPS) and provide comparative structural and kinetic analyses. Structural studies employing circular dichroism spectroscopy and analytical ultracentrifugation demonstrate that the recombinant enzymes are folded and exist as dimers in solution. Kinetic analyses of untagged and His-tagged Vc-DHDPS show that the enzymes are functional with specific activities of 75.6 U/mg and 112 U/mg, K_M (pyruvate) of 0.14 mM and 0.15 mM, K_M (L-aspartate-4-semialdehyde) of 0.08 mM and 0.09 mM, and k_cat of 34 and 46 s^-1, respectively. These results demonstrate there are no significant changes in the structure and function of Vc-DHDPS upon the addition of an N-terminal His tag and, hence, the tagged recombinant product is suitable for future studies, including screening for new inhibitors as potential anti-cholera agents. Additionally, a polyclonal antibody raised against untagged Vc-DHDPS is validated for specifically detecting recombinant and native forms of the enzyme.

Dynamic Modelling Reveals 'Hotspots' on the Pathway to Enzyme-Substrate Complex Formation

Dihydrodipicolinate synthase (DHDPS) catalyzes the first committed step in the diaminopimelate pathway of bacteria, yielding amino acids required for cell wall and protein biosyntheses. The essentiality of the enzyme to bacteria, coupled with its absence in humans, validates DHDPS as an antibacterial drug target. Conventional drug design efforts have thus far been unsuccessful in identifying potent DHDPS inhibitors. Here, we make use of contemporary molecular dynamics simulation and Markov state models to explore the interactions between DHDPS from the human pathogen Staphylococcus aureus and its cognate substrate, pyruvate. Our simulations recover the crystallographic DHDPS-pyruvate complex without a priori knowledge of the final bound structure. The highly conserved residue Arg140 was found to have a pivotal role in coordinating the entry of pyruvate into the active site from bulk solvent, consistent with previous kinetic reports, indicating an indirect role for the residue in DHDPS catalysis. A metastable binding intermediate characterized by multiple points of intermolecular interaction between pyruvate and key DHDPS residue Arg140 was found to be a highly conserved feature of the binding trajectory when comparing alternative binding pathways. By means of umbrella sampling we show that these binding intermediates are thermodynamically metastable, consistent with both the available experimental data and the substrate binding model presented in this study. Our results provide insight into an important enzyme-substrate interaction in atomistic detail that offers the potential to be exploited for the discovery of more effective DHDPS inhibitors and, in a broader sense, dynamic protein-drug interactions.

Interactions between proteins and ligands underpin many important biological processes, such as binding of substrates to their cognate enzymes in the process of catalysis. These interactions are complex, often requiring several intermediate steps to fully transition into the bound state. Here, we have used computational simulation to study binding of pyruvate to Dihydrodipicolinate synthase (DHDPS), an enzyme in the bacterial diaminopimelate pathway. In bacteria, such as the human pathogen S. aureus, DHDPS functions to make building blocks necessary for protein and bacterial cell wall biosyntheses. As the enzyme is absent in humans, yet essential for bacterial growth, DHDPS is a valid target for broad-range antibiotics. However, known DHDPS inhibitors show poor potency. One avenue that has not yet been taken into consideration for inhibitor design is the dynamics of DHDPS’s interaction with its reaction substrates (e.g. pyruvate). Using molecular dynamics simulation, we find that pyruvate binding to DHDPS must pass through a transition intermediate ‘hotspot’ in which the substrate is held in place by a dense network of noncovalent bonds. Given that many of the protein residues involved in this interaction are also shared by DHDPS from many pathogenic bacteria, this binding intermediate ‘hotspot’ may help in development of better broad-range DHDPS inhibitors.

Molecular dynamics study of naturally existing cavity couplings in proteins

Couplings between protein sub-structures are a common property of protein dynamics. Some of these couplings are especially interesting since they relate to function and its regulation. In this article we have studied the case of cavity couplings because cavities can host functional sites, allosteric sites, and are the locus of interactions with the cell milieu. We have divided this problem into two parts. In the first part, we have explored the presence of cavity couplings in the natural dynamics of 75 proteins, using 20 ns molecular dynamics simulations. For each of these proteins, we have obtained two trajectories around their native state. After applying a stringent filtering procedure, we found significant cavity correlations in 60% of the proteins. We analyze and discuss the structure origins of these correlations, including neighbourhood, cavity distance, etc. In the second part of our study, we have used longer simulations (≥100 ns) from the MoDEL project, to obtain a broader view of cavity couplings, particularly about their dependence on time. Using moving window computations we explored the fluctuations of cavity couplings along time, finding that these couplings could fluctuate substantially during the trajectory, reaching in several cases correlations above 0.25/0.5. In summary, we describe the structural origin and the variations with time of cavity couplings. We complete our work with a brief discussion of the biological implications of these results.

Use of a novel microtitration protocol to obtain diffraction-quality crystals of 4-hydroxy-2-oxoglutarate aldolase from Bos taurus

The enzyme 4-hydroxy-2-oxoglutarate aldolase (HOGA) catalyses the retro-aldol degradation of 4-hydroxy-2-oxoglutarate to pyruvate and glyoxylate as part of the hydroxyproline catabolic pathway in mammals. Mutations in the coding region of the human HOGA gene are associated with primary hyperoxaluria type 3, a disease characterized by excessive oxalate production and ultimately stone deposition. Native HOGA was purified from bovine kidney using an improved and streamlined purification protocol from which two crystal forms were obtained using two different approaches. Vapour diffusion using PEG 3350 as a precipitant produced monoclinic crystals that belonged to space group C2 and diffracted to 3.5 Å resolution. By comparison, orthorhombic crystals belonging to space group I222 or I2₁2₁2₁ and diffracting to beyond 2.25 Å resolution were obtained using a novel microtitration protocol with ammonium sulfate. The latter crystal form displayed superior diffraction quality and was suitable for structural determination by X-ray crystallography.

GroEL/ES chaperonin modulates the mechanism and accelerates the rate of TIM-barrel domain folding

The GroEL/ES chaperonin system functions as a protein folding cage. Many obligate substrates of GroEL share the (βα)8 TIM-barrel fold, but how the chaperonin promotes folding of these proteins is not known. Here, we analyzed the folding of DapA at peptide resolution using hydrogen/deuterium exchange and mass spectrometry. During spontaneous folding, all elements of the DapA TIM barrel acquire structure simultaneously in a process associated with a long search time. In contrast, GroEL/ES accelerates folding more than 30-fold by catalyzing segmental structure formation in the TIM barrel. Segmental structure formation is also observed during the fast spontaneous folding of a structural homolog of DapA from a bacterium that lacks GroEL/ES. Thus, chaperonin independence correlates with folding properties otherwise enforced by protein confinement in the GroEL/ES cage. We suggest that folding catalysis by GroEL/ES is required by a set of proteins to reach native state at a biologically relevant timescale, avoiding aggregation or degradation.

Crystal structure and in silico studies of dihydrodipicolinate synthase (DHDPS) from Aquifex aeolicus

Dihydrodipicolinate synthase (DHDPS, E.C.4.2.1.52) catalyzes the first committed step in the lysine biosynthetic pathway: the condensation of (S)-aspartate semialdehyde and pyruvate to form (4S)-4-hydroxy-2,3,4,5-tetrahydro-(2S)-dipicolinic acid. Since (S)-lysine biosynthesis does not occur in animals, DHDPS is an attractive target for rational antibiotic and herbicide design. Here, we report the crystal structure of DHDPS from a hyperthermophilic bacterium Aquifex aeolicus (AqDHDPS). L-Lysine is used as an important animal feed additive where the production is at the level of 1.5 million tons per year. The biotechnological manufacture of lysine has been going for more than 50 years which includes over synthesis and reverse engineering of DHDPS. AqDHDPS revealed a unique disulfide linkage which is not conserved in the homologues of AqDHDPS. In silico mutation of C139A and intermolecular ion-pair residues and the subsequent molecular dynamics simulation of the mutants showed that these residues are critical for the stability of AqDHDPS tetramer. MD simulations of AqDHDPS at three different temperatures (303, 363 and 393 K) revealed that the molecule is stable at 363 K. Thus, this structural and in silico study of AqDHDPS likely provides additional details towards the rational and structure-based design of hyper-L-lysine producing bacterial strains.

From knock-out phenotype to three-dimensional structure of a promising antibiotic target from Streptococcus pneumoniae

Given the rise in drug-resistant Streptococcus pneumoniae, there is an urgent need to discover new antimicrobials targeting this pathogen and an equally urgent need to characterize new drug targets. A promising antibiotic target is dihydrodipicolinate synthase (DHDPS), which catalyzes the rate-limiting step in lysine biosynthesis. In this study, we firstly show by gene knock out studies that S. pneumoniae (sp) lacking the DHDPS gene is unable to grow unless supplemented with lysine-rich media. We subsequently set out to characterize the structure, function and stability of the enzyme drug target. Our studies show that sp-DHDPS is folded and active with a k(cat) = 22 s(-1), K(M)(PYR) = 2.55 ± 0.05 mM and K(M)(ASA) = 0.044 ± 0.003 mM. Thermal denaturation experiments demonstrate sp-DHDPS exhibits an apparent melting temperature (T(M)(app)) of 72 °C, which is significantly greater than Escherichia coli DHDPS (Ec-DHDPS) (T(M)(app) = 59 °C). Sedimentation studies show that sp-DHDPS exists in a dimer-tetramer equilibrium with a K(D)(4→2) = 1.7 nM, which is considerably tighter than its E. coli ortholog (K(D)(4→2) = 76 nM). To further characterize the structure of the enzyme and probe its enhanced stability, we solved the high resolution (1.9 Å) crystal structure of sp-DHDPS (PDB ID 3VFL). The enzyme is tetrameric in the crystal state, consistent with biophysical measurements in solution. Although the sp-DHDPS and Ec-DHDPS active sites are almost identical, the tetramerization interface of the s. pneumoniae enzyme is significantly different in composition and has greater buried surface area (800 Å(2)) compared to its E. coli counterpart (500 Å(2)). This larger interface area is consistent with our solution studies demonstrating that sp-DHDPS is considerably more thermally and thermodynamically stable than Ec-DHDPS. Our study describe for the first time the knock-out phenotype, solution properties, stability and crystal structure of DHDPS from S. pneumoniae, a promising antimicrobial target.

Crystal, solution and in silico structural studies of dihydrodipicolinate synthase from the common grapevine

Dihydrodipicolinate synthase (DHDPS) catalyzes the rate limiting step in lysine biosynthesis in bacteria and plants. The structure of DHDPS has been determined from several bacterial species and shown in most cases to form a homotetramer or dimer of dimers. However, only one plant DHDPS structure has been determined to date from the wild tobacco species, Nicotiana sylvestris (Blickling et al. (1997) J. Mol. Biol. 274, 608-621). Whilst N. sylvestris DHDPS also forms a homotetramer, the plant enzyme adopts a 'back-to-back' dimer of dimers compared to the 'head-to-head' architecture observed for bacterial DHDPS tetramers. This raises the question of whether the alternative quaternary architecture observed for N. sylvestris DHDPS is common to all plant DHDPS enzymes. Here, we describe the structure of DHDPS from the grapevine plant, Vitis vinifera, and show using analytical ultracentrifugation, small-angle X-ray scattering and X-ray crystallography that V. vinifera DHDPS forms a 'back-to-back' homotetramer, consistent with N. sylvestris DHDPS. This study is the first to demonstrate using both crystal and solution state measurements that DHDPS from the grapevine plant adopts an alternative tetrameric architecture to the bacterial form, which is important for optimizing protein dynamics as suggested by molecular dynamics simulations reported in this study.