Honggao Yan

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.

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.

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.

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. 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. 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.

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.

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.

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.

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 kcat and a dramatic increase in Km, indicating Glu64 is important for both binding and catalysis in the activation of 5FC. 19F 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. 1H and 15N 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 1H and 15N 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.

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.