Structure-Based Design of a Dual-Targeted Covalent Inhibitor Against Papain-like and Main Proteases of SARS-CoV-2.

The two proteases, PLpro and Mpro, of SARS-CoV-2 are essential for replication of the virus. Using a structure-based co-pharmacophore screening approach, we developed a novel dual-targeted inhibitor that is equally potent in inhibiting PLpro and Mpro of SARS-CoV-2. The inhibitor contains a novel warhead, which can form a covalent bond with the catalytic cysteine residue of either enzyme. The maximum rate of the covalent inactivation is comparable to that of the most potent inhibitors reported for the viral proteases and covalent inhibitor drugs currently in clinical use. The covalent inhibition appears to be very specific for the viral proteases. The inhibitor has a potent antiviral activity against SARS-CoV-2 and is also well tolerated by mice and rats in toxicity studies. These results suggest that the inhibitor is a promising lead for development of drugs for treatment of COVID-19.

Combined QM(DFT)/MM molecular dynamics simulations of the deamination of cytosine by yeast cytosine deaminase (yCD).

Extensive combined quantum mechanical (B3LYP/6-31G*) and molecular mechanical (QM/MM) molecular dynamics simulations have been performed to elucidate the hydrolytic deamination mechanism of cytosine to uracil catalyzed by the yeast cytosine deaminase (yCD). Though cytosine has no direct binding to the zinc center, it reacts with the water molecule coordinated to zinc, and the adjacent conserved Glu64 serves as a general acid/base to shuttle protons from water to cytosine. The overall reaction consists of several proton-transfer processes and nucleophilic attacks. A tetrahedral intermediate adduct of cytosine and water binding to zinc is identified and similar to the crystal structure of yCD with the inhibitor 2-pyrimidinone. The rate-determining step with the barrier of 18.0 kcal/mol in the whole catalytic cycle occurs in the process of uracil departure where the proton transfer from water to Glu64 and nucleophilic attack of the resulting hydroxide anion to C2 of the uracil ring occurs synchronously. © 2016 Wiley Periodicals, Inc.

Host target modification as a strategy to counter pathogen hijacking of the jasmonate hormone receptor.

In the past decade, characterization of the host targets of pathogen virulence factors took a center stage in the study of pathogenesis and disease susceptibility in plants and humans. However, the impressive knowledge of host targets has not been broadly exploited to inhibit pathogen infection. Here, we show that host target modification could be a promising new approach to "protect" the disease-vulnerable components of plants. In particular, recent studies have identified the plant hormone jasmonate (JA) receptor as one of the common targets of virulence factors from highly evolved biotrophic/hemibiotrophic pathogens. Strains of the bacterial pathogen Pseudomonas syringae, for example, produce proteinaceous effectors, as well as a JA-mimicking toxin, coronatine (COR), to activate JA signaling as a mechanism to promote disease susceptibility. Guided by the crystal structure of the JA receptor and evolutionary clues, we succeeded in modifying the JA receptor to allow for sufficient endogenous JA signaling but greatly reduced sensitivity to COR. Transgenic Arabidopsis expressing this modified receptor not only are fertile and maintain a high level of insect defense, but also gain the ability to resist COR-producing pathogens Pseudomonas syringae pv. tomato and P. syringae pv. maculicola. Our results provide a proof-of-concept demonstration that host target modification can be a promising new approach to prevent the virulence action of highly evolved pathogens.

Molecular dynamics simulations of the Escherichia coli HPPK apo-enzyme reveal a network of conformational transitions.

6-Hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK) catalyzes the first reaction in the folate biosynthetic pathway. Comparison of its X-ray and nuclear magnetic resonance structures suggests that the enzyme undergoes significant conformational change upon binding to its substrates, especially in three catalytic loops. Experimental research has shown that even when confined by crystal contacts, loops 2 and 3 remain rather flexible when the enzyme is in its apo form, raising questions about the putative large-scale induced-fit conformational change of HPPK. To investigate the loop dynamics in a crystal-free environment, we performed conventional molecular dynamics simulations of the apo-enzyme at two different temperatures (300 and 350 K). Our simulations show that the crystallographic B-factors considerably underestimate the loop dynamics; multiple conformations of loops 2 and 3, including the open, semi-open, and closed conformations that an enzyme must adopt throughout its catalytic cycle, are all accessible to the apo-enzyme. These results revise our previous view of the functional mechanism of conformational change upon MgATP binding and offer valuable structural insights into the workings of HPPK. In this paper, conformational network analysis and principal component analysis related to the loops are discussed to support the presented conclusions.

Conformational transition in signal transduction: metastable states and transition pathways in the activation of a signaling protein.

Signal transduction is of vital importance to the growth and adaptation of living organisms. The key to understand mechanisms of biological signal transduction is elucidation of the conformational dynamics of its signaling proteins, as the activation of a signaling protein is fundamentally a process of conformational transition from an inactive to an active state. A predominant form of signal transduction for bacterial sensing of environmental changes in the wild or inside their hosts is a variety of two-component systems, in which the conformational transition of a response regulator (RR) from an inactive to an active state initiates responses to the environmental changes. Here, RR activation has been investigated using RR468 as a model system by extensive unbiased all-atom molecular dynamics (MD) simulations in explicit solvent, starting from snapshots along a targeted MD trajectory that covers the conformational transition. Markov state modeling, transition path theory, and geometric analyses of the wealth of the MD data have provided a comprehensive description of the RR activation. It involves a network of metastable states, with one metastable state essentially the same as the inactive state and another very similar to the active state that are connected via a small set of intermediates. Five major pathways account for >75% of the fluxes of the conformational transition from the inactive to the active-like state. The thermodynamic stability of the states and the activation barriers between states are found, to identify rate-limiting steps. The conformal transition is initiated predominantly by movements of the ß3α3 loop, followed by movements of the ß4α4-loop and neighboring α4 helix region, and capped by additional movements of the ß3α3 loop. A number of transient hydrophobic and hydrogen bond interactions are revealed, and they may be important for the conformational transition.

Structural determinants of host specificity of complement Factor H recruitment by Streptococcus pneumoniae.

Many human pathogens have strict host specificity, which affects not only their epidemiology but also the development of animal models and vaccines. Complement Factor H (FH) is recruited to pneumococcal cell surface in a human-specific manner via the N-terminal domain of the pneumococcal protein virulence factor choline-binding protein A (CbpAN). FH recruitment enables Streptococcus pneumoniae to evade surveillance by human complement system and contributes to pneumococcal host specificity. The molecular determinants of host specificity of complement evasion are unknown. In the present study, we show that a single human FH (hFH) domain is sufficient for tight binding of CbpAN, present the crystal structure of the complex and identify the critical structural determinants for host-specific FH recruitment. The results offer new approaches to the development of better animal models for pneumococcal infection and redesign of the virulence factor for pneumococcal vaccine development and reveal how FH recruitment can serve as a mechanism for both pneumococcal complement evasion and adherence.

Crystallographic and molecular dynamics simulation analysis of Escherichia coli dihydroneopterin aldolase.

BACKGROUND: Dihydroneopterin aldolase (DHNA) catalyzes the conversion of 7,8-dihydroneopterin to 6-hydroxymethyl-7,8-dihydropterin and also the epimerization of DHNP to 7,8-dihydromonapterin. Previously, we determined the crystal structure of Staphylococcus aureus DHNA (SaDHNA) in complex with the substrate analogue neopterin (NP). We also showed that Escherichia coli DHNA (EcDHNA) and SaDHNA have significantly different binding and catalytic properties by biochemical analysis. On the basis of these structural and functional data, we proposed a catalytic mechanism involving two proton wires. RESULTS: To understand the structural basis for the biochemical differences and further investigate the catalytic mechanism of DHNA, we have determined the structure of EcDHNA complexed with NP at 1.07-Å resolution [PDB:2O90], built an atomic model of EcDHNA complexed with the substrate DHNP, and performed molecular dynamics (MD) simulation analysis of the substrate complex. EcDHNA has the same fold as SaDHNA and also forms an octamer that consists of two tetramers, but the packing of one tetramer with the other is significantly different between the two enzymes. Furthermore, the structures reveal significant differences in the vicinity of the active site, particularly in the loop that connects strands ß3 and ß4, mainly due to the substitution of nearby residues. The building of an atomic model of the complex of EcDHNA and the substrate DHNP and the MD simulation of the complex show that some of the hydrogen bonds between the substrate and the enzyme are persistent, whereas others are transient. The substrate binding model and MD simulation provide the molecular basis for the biochemical behaviors of the enzyme, including noncooperative substrate binding, indiscrimination of a pair of epimers as the substrates, proton wire switching during catalysis, and formation of epimerization product. CONCLUSIONS: The EcDHNA and SaDHNA structures, each in complex with NP, reveal the basis for the biochemical differences between EcDHNA and SaDHNA. The atomic substrate binding model and MD simulation offer insights into substrate binding and catalysis by DHNA. The EcDHNA structure also affords an opportunity to develop antimicrobials specific for Gram-negative bacteria, as DHNAs from Gram-negative bacteria are highly homologous and E. coli is a representative of this class of bacteria.

Structural enzymology and inhibition of the bi-functional folate pathway enzyme HPPK-DHPS from the biowarfare agent Francisella tularensis.

UNLABELLED: Two valid targets for antibiotic development, 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK) and dihydropteroate synthase (DHPS), catalyze consecutive reactions in folate biosynthesis. In Francisella tularensis (Ft), these two activities are contained in a single protein, FtHPPK-DHPS. Although Pemble et al. (PLoS One 5, e14165) determined the structure of FtHPPK-DHPS, they were unable to measure the kinetic parameters of the enzyme. In this study, we elucidated the binding and inhibitory activities of two HPPK inhibitors (HP-18 and HP-26) against FtHPPK-DHPS, determined the structure of FtHPPK-DHPS in complex with HP-26, and measured the kinetic parameters for the dual enzymatic activities of FtHPPK-DHPS. The biochemical analyses showed that HP-18 and HP-26 have significant isozyme selectivity, and that FtHPPK-DHPS is unique in that the catalytic efficiency of its DHPS activity is only 1/260,000 of that of Escherichia coli DHPS. Sequence and structural analyses suggest that HP-26 is an excellent lead for developing therapeutic agents for tularemia, and that the very low DHPS activity is due, at least in part, to the lack of a key residue that interacts with the substrate p-aminobenzoic acid (pABA). A BLAST search of the genomes of ten F. tularensis strains indicated that the bacterium contains a single FtHPPK-DHPS. The marginal DHPS activity and the single copy existence of FtHPPK-DHPS in F. tularensis make this bacterium more vulnerable to DHPS inhibitors. Current sulfa drugs are ineffective against tularemia; new inhibitors targeting the unique pABA-binding pocket may be effective and less subject to resistance because any mutations introducing resistance may make the marginal DHPS activity unable to support the growth of F. tularensis. DATABASE: The coordinates and structure factors have been deposited in the Protein Data Bank under accession code 4PZV.

Conformational transition of response regulator RR468 in a two-component system signal transduction process.

Signal transduction can be accomplished via a two-component system (TCS) consisting of a histidine kinase (HK) and a response regulator (RR). In this work, we simulate the response regulator RR468 from Thermotoga maritima, in which phosphorylation and dephosphorylation of a conserved aspartate residue acts as a switch via a large conformational change concentrated in three proximal loops. A detailed view of the conformational transition is obscured by the lack of stability of the intermediate states, which are difficult to detect using common structural biology techniques. Molecular dynamics (MD) trajectories of the inactive and active conformations were run, and show that the inactive (or active) trajectories do not exhibit sampling of the active (or inactive) conformations on this time scale. Targeted MD (TMD) was used to generate trajectories that span the inactive and active conformations and provide a view of how a localized event like phosphorylation can lead to conformational changes elsewhere in the protein, especially in the three proximal loops. The TMD trajectories are clustered to identify stages along the transition path. Residue interaction networks are identified that point to key residues having to rearrange in the process of transition. These are identified using both hydrogen bond analysis and residue interaction strength measurements. Potentials of mean force are generated for key residue rearrangements to ascertain their free energy barriers. We introduce methods that attempt to extrapolate from one conformation to the other and find that the most fluctuating proximal loop can transit part way from one to the other, suggesting that this conformational information is embedded in the sequence.

Structure of bacterial transcription factor SpoIIID and evidence for a novel mode of DNA binding.

SpoIIID is evolutionarily conserved in endospore-forming bacteria, and it activates or represses many genes during sporulation of Bacillus subtilis. An SpoIIID monomer binds DNA with high affinity and moderate sequence specificity. In addition to a predicted helix-turn-helix motif, SpoIIID has a C-terminal basic region that contributes to DNA binding. The nuclear magnetic resonance (NMR) solution structure of SpoIIID in complex with DNA revealed that SpoIIID does indeed have a helix-turn-helix domain and that it has a novel C-terminal helical extension. Residues in both of these regions interact with DNA, based on the NMR data and on the effects on DNA binding in vitro of SpoIIID with single-alanine substitutions. These data, as well as sequence conservation in SpoIIID binding sites, were used for information-driven docking to model the SpoIIID-DNA complex. The modeling resulted in a single cluster of models in which the recognition helix of the helix-turn-helix domain interacts with the major groove of DNA, as expected. Interestingly, the C-terminal extension, which includes two helices connected by a kink, interacts with the adjacent minor groove of DNA in the models. This predicted novel mode of binding is proposed to explain how a monomer of SpoIIID achieves high-affinity DNA binding. Since SpoIIID is conserved only in endospore-forming bacteria, which include important pathogenic Bacilli and Clostridia, whose ability to sporulate contributes to their environmental persistence, the interaction of the C-terminal extension of SpoIIID with DNA is a potential target for development of sporulation inhibitors.

Feedback inhibition of deoxy-D-xylulose-5-phosphate synthase regulates the methylerythritol 4-phosphate pathway.

The 2-C-methyl-D-erythritol 4-phosphate (MEP) pathway leads to the biosynthesis of isopentenyl diphosphate (IDP) and dimethylallyl diphosphate (DMADP), the precursors for isoprene and higher isoprenoids. Isoprene has significant effects on atmospheric chemistry, whereas other isoprenoids have diverse roles ranging from various biological processes to applications in commercial uses. Understanding the metabolic regulation of the MEP pathway is important considering the numerous applications of this pathway. The 1-deoxy-D-xylulose-5-phosphate synthase (DXS) enzyme was cloned from Populus trichocarpa, and the recombinant protein (PtDXS) was purified from Escherichia coli. The steady-state kinetic parameters were measured by a coupled enzyme assay. An LC-MS/MS-based assay involving the direct quantification of the end product of the enzymatic reaction, 1-deoxy-D-xylulose 5-phosphate (DXP), was developed. The effect of different metabolites of the MEP pathway on PtDXS activity was tested. PtDXS was inhibited by IDP and DMADP. Both of these metabolites compete with thiamine pyrophosphate for binding with the enzyme. An atomic structural model of PtDXS in complex with thiamine pyrophosphate and Mg(2+) was built by homology modeling and refined by molecular dynamics simulations. The refined structure was used to model the binding of IDP and DMADP and indicated that IDP and DMADP might bind with the enzyme in a manner very similar to the binding of thiamine pyrophosphate. The feedback inhibition of PtDXS by IDP and DMADP constitutes an important mechanism of metabolic regulation of the MEP pathway and indicates that thiamine pyrophosphate-dependent enzymes may often be affected by IDP and DMADP.

Bisubstrate analog inhibitors of 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase: new lead exhibits a distinct binding mode.

6-Hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK), a key enzyme in the folate biosynthesis pathway catalyzing the pyrophosphoryl transfer from ATP to 6-hydroxymethyl-7,8-dihydropterin, is an attractive target for developing novel antimicrobial agents. Previously, we studied the mechanism of HPPK action, synthesized bisubstrate analog inhibitors by linking 6-hydroxymethylpterin to adenosine through phosphate groups, and developed a new generation of bisubstrate inhibitors by replacing the phosphate bridge with a piperidine-containing linkage. To further improve linker properties, we have synthesized a new compound, characterized its protein binding/inhibiting properties, and determined its structure in complex with HPPK. Surprisingly, this inhibitor exhibits a new binding mode in that the adenine base is flipped when compared to previously reported structures. Furthermore, the side chain of amino acid residue E77 is involved in protein-inhibitor interaction, forming hydrogen bonds with both 2' and 3' hydroxyl groups of the ribose moiety. Residue E77 is conserved among HPPK sequences, but interacts only indirectly with the bound MgATP via water molecules. Never observed before, the E77-ribose interaction is compatible only with the new inhibitor-binding mode. Therefore, this compound represents a new direction for further development.

Role of glutamate 64 in the activation of the prodrug 5-fluorocytosine by yeast cytosine deaminase.

Yeast cytosine deaminase (yCD) catalyzes the hydrolytic deamination of cytosine to uracil as well as the deamination of the prodrug 5-fluorocytosine (5FC) to the anticancer drug 5-fluorouracil. In this study, the role of Glu64 in the activation of the prodrug 5FC was investigated by site-directed mutagenesis, biochemical, nuclear magnetic resonance (NMR), and computational studies. Steady-state kinetics studies showed that the mutation of Glu64 causes a dramatic decrease in k(cat) and a dramatic increase in K(m), indicating Glu64 is important for both binding and catalysis in the activation of 5FC. (19)F NMR experiments showed that binding of the inhibitor 5-fluoro-1H-pyrimidin-2-one (5FPy) to the wild-type yCD causes an upfield shift, indicating that the bound inhibitor is in the hydrated form, mimicking the transition state or the tetrahedral intermediate in the activation of 5FC. However, binding of 5FPy to the E64A mutant enzyme causes a downfield shift, indicating that the bound 5FPy remains in an unhydrated form in the complex with the mutant enzyme. (1)H and (15)N NMR analysis revealed trans-hydrogen bond D/H isotope effects on the hydrogen of the amide of Glu64, indicating that the carboxylate of Glu64 forms two hydrogen bonds with the hydrated 5FPy. ONIOM calculations showed that the wild-type yCD complex with the hydrated form of the inhibitor 1H-pyrimidin-2-one is more stable than the initial binding complex, and in contrast, with the E64A mutant enzyme, the hydrated inhibitor is no longer favored and the conversion has a higher activation energy, as well. The hydrated inhibitor is stabilized in the wild-type yCD by two hydrogen bonds between it and the carboxylate of Glu64 as revealed by (1)H and (15)N NMR analysis. To explore the functional role of Glu64 in catalysis, we investigated the deamination of cytosine catalyzed by the E64A mutant by ONIOM calculations. The results showed that without the assistance of Glu64, both proton transfers before and after the formation of the tetrahedral reaction intermediate become partially rate-limiting steps. The results of the experimental and computational studies together indicate that Glu64 plays a critical role in both the binding and the chemical transformation in the conversion of the prodrug 5FC to the anticancer drug 5-fluorouracil.

Bisubstrate analogue inhibitors of 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase: New design with improved properties.

6-Hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK), a key enzyme in the folate biosynthetic pathway, catalyzes the pyrophosphoryl transfer from ATP to 6-hydroxymethyl-7,8-dihydropterin. The enzyme is essential for microorganisms, is absent from humans, and is not the target for any existing antibiotics. Therefore, HPPK is an attractive target for developing novel antimicrobial agents. Previously, we characterized the reaction trajectory of HPPK-catalyzed pyrophosphoryl transfer and synthesized a series of bisubstrate analog inhibitors of the enzyme by linking 6-hydroxymethylpterin to adenosine through 2, 3, or 4 phosphate groups. Here, we report a new generation of bisubstrate analog inhibitors. To improve protein binding and linker properties of such inhibitors, we have replaced the pterin moiety with 7,7-dimethyl-7,8-dihydropterin and the phosphate bridge with a piperidine linked thioether. We have synthesized the new inhibitors, measured their K(d) and IC(50) values, determined their crystal structures in complex with HPPK, and established their structure-activity relationship. 6-Carboxylic acid ethyl ester-7,7-dimethyl-7,8-dihydropterin, a novel intermediate that we developed recently for easy derivatization at position 6 of 7,7-dimethyl-7,8-dihydropterin, offers a much high yield for the synthesis of bisubstrate analogs than that of previously established procedure.

Probing single-molecule enzyme active-site conformational state intermittent coherence.

The relationship between protein conformational dynamics and enzymatic reactions has been a fundamental focus in modern enzymology. Using single-molecule fluorescence resonance energy transfer (FRET) with a combined statistical data analysis approach, we have identified the intermittently appearing coherence of the enzymatic conformational state from the recorded single-molecule intensity-time trajectories of enzyme 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK) in catalytic reaction. The coherent conformational state dynamics suggests that the enzymatic catalysis involves a multistep conformational motion along the coordinates of substrate-enzyme complex formation and product releasing, presenting as an extreme dynamic behavior intrinsically related to the time bunching effect that we have reported previously. The coherence frequency, identified by statistical results of the correlation function analysis from single-molecule FRET trajectories, increases with the increasing substrate concentrations. The intermittent coherence in conformational state changes at the enzymatic reaction active site is likely to be common and exist in other conformation regulated enzymatic reactions. Our results of HPPK interaction with substrate support a multiple-conformational state model, being consistent with a complementary conformation selection and induced-fit enzymatic loop-gated conformational change mechanism in substrate-enzyme active complex formation.

Role of protein conformational dynamics in the catalysis by 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase.

Enzymatic catalysis has conflicting structural requirements of the enzyme. In order for the enzyme to form a Michaelis complex, the enzyme must be in an open conformation so that the substrate can get into its active center. On the other hand, in order to maximize the stabilization of the transition state of the enzymatic reaction, the enzyme must be in a closed conformation to maximize its interactions with the transition state. The conflicting structural requirements can be resolved by a flexible active center that can sample both open and closed conformational states. For a bisubstrate enzyme, the Michaelis complex consists of two substrates in addition to the enzyme. The enzyme must remain flexible upon the binding of the first substrate so that the second substrate can get into the active center. The active center is fully assembled and stabilized only when both substrates bind to the enzyme. However, the side-chain positions of the catalytic residues in the Michaelis complex are still not optimally aligned for the stabilization of the transition state, which lasts only approximately 10(-13) s. The instantaneous and optimal alignment of catalytic groups for the transition state stabilization requires a dynamic enzyme, not an enzyme which undergoes a large scale of movements but an enzyme which permits at least a small scale of adjustment of catalytic group positions. This review will summarize the structure, catalytic mechanism, and dynamic properties of 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase and examine the role of protein conformational dynamics in the catalysis of a bisubstrate enzymatic reaction.

Productive versus unproductive nucleotide binding in yeast guanylate kinase mutants: comparison of R41M with K14M by proton two dimensional transferred NOESY.

The R41M and K14M mutant enzymes of yeast guanylate kinase (GKy) were studied to investigate the effects of these site-directed mutations on bound-substrate conformations. Published X-ray crystal structures of yeast guanylate kinase indicate that K14 is part of the "P" loop involved in ATP and ADP binding, while R41 is suggested as a hydrogen bonding partner for the phosphoryl moiety of GMP. Both of these residues might be involved in transition state stabilization. Adenosine conformations of ATP and ADP and guanosine conformations of GMP bound to R41M and K14M mutant yeast guanylate kinase in the complexes GKy.MgATP, GKy.MgADP, and GKy.MgADP.[u-(13)C]GMP were determined by two-dimensional transferred nuclear Overhauser effect (TRNOESY) measurements combined with molecular dynamics simulations, and these conformations were compared with previously published conformations for the wild type. In the fully constrained, two substrate complexes, GKy.MgADP.[u-(13)C]GMP, the guanyl glycosidic torsion angle, chi, is 51 +/- 5 degrees for R41M and 47 +/- 5 degrees for K14M. Both are similar to the published 50 +/- 5 degrees published for wild type. For R41M with adenyl nucleotides, the glycosidic torsion angle, chi, was 55 +/- 5 degrees with MgATP, and 47 +/- 5 degrees with MgADP, which compares well to the 54 +/- 5 degrees published for wild type. However, for K14M with adenyl nucleotides, the glycosidic torsion angle was 30 +/- 5 degrees with MgATP and 28 +/- 5 degrees with MgADP. The results indicate that bound adenyl-nucleotides have significantly different conformations in the wild-type and K14M mutant enzymes, suggesting that K14 plays an important role in orienting the triphosphate of MgATP for catalysis.

Dynamics in Thermotoga neapolitana adenylate kinase: 15N relaxation and hydrogen-deuterium exchange studies of a hyperthermophilic enzyme highly active at 30 degrees C.

Backbone conformational dynamics of Thermotoga neapolitana adenylate kinase in the free form (TNAK) and inhibitor-bound form (TNAK*Ap5A) were investigated at 30 degrees C using (15)N NMR relaxation measurements and NMR monitored hydrogen-deuterium exchange. With kinetic parameters identical to those of Escherichia coli AK (ECAK) at 30 degrees C, TNAK is a unique hyperthermophilic enzyme. These catalytic properties make TNAK an interesting and novel model to study the interplay between protein rigidity, stability, and activity. Comparison of fast time scale dynamics (picosecond to nanosecond) in the open and closed states of TNAK and ECAK at 30 degrees C reveals a uniformly higher rigidity across all domains of TNAK. Within this framework of a rigid TNAK structure, several residues located in the AMP-binding domain and in the core-lid hinge regions display high picosecond to nanosecond time scale flexibility. Together with the recent comparison of ECAK dynamics with those of hyperthermophilic Aquifex aeolicus AK (AAAK), our results provide strong evidence for the role of picosecond to nanosecond time scale fluctuations in both stability and activity. In the slow time scales, TNAK's increased rigidity is not uniform but localized in the AMP-binding and lid domains. The core domain amides of ECAK and TNAK in the open and closed states show comparable protection against exchange. Significantly, the hinges framing the lid domain show similar exchange data in ECAK and TNAK open and closed forms. Our NMR relaxation and hydrogen-deuterium exchange studies therefore suggest that TNAK maintains high activity at 30 degrees C by localizing flexibility to the hinge regions that are key to facilitating conformational changes.

Scoring ligand similarity in structure-based virtual screening.

Scoring to identify high-affinity compounds remains a challenge in virtual screening. On one hand, protein-ligand scoring focuses on weighting favorable and unfavorable interactions between the two molecules. Ligand-based scoring, on the other hand, focuses on how well the shape and chemistry of each ligand candidate overlay on a three-dimensional reference ligand. Our hypothesis is that a hybrid approach, using ligand-based scoring to rank dockings selected by protein-ligand scoring, can ensure that high-ranking molecules mimic the shape and chemistry of a known ligand while also complementing the binding site. Results from applying this approach to screen nearly 70 000 National Cancer Institute (NCI) compounds for thrombin inhibitors tend to support the hypothesis. EON ligand-based ranking of docked molecules yielded the majority (4/5) of newly discovered, low to mid-micromolar inhibitors from a panel of 27 assayed compounds, whereas ranking docked compounds by protein-ligand scoring alone resulted in one new inhibitor. Since the results depend on the choice of scoring function, an analysis of properties was performed on the top-scoring docked compounds according to five different protein-ligand scoring functions, plus EON scoring using three different reference compounds. The results indicate that the choice of scoring function, even among scoring functions measuring the same types of interactions, can have an unexpectedly large effect on which compounds are chosen from screening. Furthermore, there was almost no overlap between the top-scoring compounds from protein-ligand versus ligand-based scoring, indicating the two approaches provide complementary information. Matchprint analysis, a new addition to the SLIDE (Screening Ligands by Induced-fit Docking, Efficiently) screening toolset, facilitated comparison of docked molecules' interactions with those of known inhibitors. The majority of interactions conserved among top-scoring compounds for a given scoring function, and from the different scoring functions, proved to be conserved interactions in known inhibitors. This was particularly true in the S1 pocket, which was occupied by all the docked compounds.

Dynamics of the conformational transitions in the assembling of the Michaelis complex of a bisubstrate enzyme: a (15)N relaxation study of Escherichia coli 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase.

6-Hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK) catalyzes the transfer of pyrophosphate from ATP to 6-hydroxymethyl-7,8-dihydropterin (HP), which follows an ordered bi-bi kinetic mechanism with ATP binding to the enzyme first. HPPK undergoes dramatic conformational changes during its catalytic cycle as revealed by X-ray crystallography, and the conformational changes are essential for the enzymatic catalysis as shown by site-directed mutagenesis and biochemical and crystallographic analysis of the mutants. However, the dynamic properties of the enzyme have not been measured experimentally. Here, we report a (15)N NMR relaxation study of the dynamic properties of Escherichia coli HPPK from the apo form to the binary substrate complex with MgATP (represented by MgAMPCPP, an ATP analogue) to the Michaelis complex (ternary substrate complex) with MgATP (represented by MgAMPCPP) and HP (represented by 7,7-dimethyl-6-hydroxypterin, an HP analogue). The results show that the binding of the nucleotide to HPPK does not cause major changes in the dynamic properties of the enzyme. Whereas enzymes are often more rigid when bound to the ligand or the substrate, the internal mobility of HPPK is not reduced and is even moderately increased in the binary complex, particularly in the catalytic loops. The internal mobility of the catalytic loops is significantly quenched upon the formation of the ternary complex, but some mobility remains. The enhanced motions in the catalytic loops of the binary substrate complex may be required for the assembling of the ternary complex. On the other hand, some degrees of mobility in the catalytic loops of the ternary complex may be required for the optimal stabilization of the transition state, which may need the instantaneous adjustment and alignment of the side-chain positions of catalytic residues. Such dynamic behaviors may be characteristic of bisubstrate enzymes.