Elizabeth E. Howell

Neutron diffraction studies of Escherichia coli dihydrofolate reductase complexed with methotrexate.

Hydrogen atoms play a central role in many biochemical processes yet are difficult to visualize by x-ray crystallography. Spallation neutron sources provide a new arena for protein crystallography with TOF measurements enhancing data collection efficiency and allowing hydrogen atoms to be located in smaller crystals of larger biological macromolecules. Here we report a 2.2-Å resolution neutron structure of Escherichia coli dihydrofolate reductase (DHFR) in complex with methotrexate (MTX). Neutron data were collected on a 0.3-mm3 D2O-soaked crystal at the Los Alamos Neutron Scattering Center. This study provides an example of using spallation neutrons to study protein dynamics, to identify protonation states directly from nuclear density maps, and to analyze solvent structure. Our structure reveals that the occluded loop conformation [monomer (mon.) A] of the DHFR·MTX complex undergoes greater H/D exchange compared with the closed-loop conformer (mon. B), partly because the Met-20 and β(F-G) loops readily exchange in mon. A. The eight-stranded β sheet of both DHFR molecules resists H/D exchange more than the helices and loops. However, the C-terminal strand, βH, in mon. A is almost fully exchanged. Several D2Os form hydrogen bonds with exchanged amides. At the active site, the N1 atom of MTX is protonated and thus charged when bound to DHFR. Several D2Os are observed at hydrophobic surfaces, including two pockets near the MTX-binding site. A previously unidentified D2O hydrogen bonds with the catalytic D27 in mon. B, stabilizing its negative charge.

Searching Sequence Space: Two Different Approaches to Dihydrofolate Reductase Catalysis

There are numerous examples of proteins that catalyze the same reaction while possessing different structures. This review focuses on two dihydrofolate reductases (DHFRs) that have disparate structures and discusses how the catalytic strategies of these two DHFRs are driven by their respective scaffolds. The two proteins are E. coli chromosomal DHFR (Ec DHFR) and a type II R‐plasmid‐encoded DHFR, typified by R67 DHFR. The former has been described as a very well evolved enzyme with an efficiency of 0.15, while the latter has been suggested to be a model for a “primitive” enzyme that has not yet been optimized by evolution. This comparison underlines what is important to catalysis in these two enzymes and concurrently highlights fundamental issues in enzyme catalysis.

STAAR: Statistical analysis of aromatic rings

The statistical analysis of aromatic rings program allows for an automated search for anion‐π interactions between phenylalanine residues and carboxylic acid moieties of neighboring aspartic acid or glutamic acid residues in protein data bank (PDB) structures. The program is written in C++ and is available both as a standalone code and through a web implementation that allows users to upload and analyze biomolecular structures in PDB format. The program outputs lists of Phe/Glu or Phe/Asp pairs involved in potential anion‐π interactions, together with geometrical (distance and angle between the Phes center of mass and Glu or Asps center of charge) and energetic (quantum mechanical Kitaura‐Morokuma interaction energy between the residues) descriptions of each anion‐π interaction. Application of the program on the latest content of the PDB shows that anion‐π interactions are present in thousands of protein structures and can possess strong energies, as low as −8.72 kcal/mol.

Toward resolving the catalytic mechanism of dihydrofolate reductase using neutron and ultrahigh-resolution X-ray crystallography.

Significance There is immense difficulty in mapping out the complete details of an enzyme’s mechanism, especially those that catalyze an acid-base reaction, owing to the simple fact that hydrogen atom positions are rarely known with any confidence. Ultrahigh-resolution X-ray and, better still, neutron crystallography can provide this crucial layer of information. We paired these techniques to reveal the catalytic mechanism of dihydrofolate reductase (DHFR), an enzyme necessary for nucleotide biosynthesis and a classical drug target. In a complex that closely resembles the catalytically active state, DHFR stabilizes a particular substrate conformer and likely elevates the pKa of the substrate atom that is protonated. This protonation occurs directly via water, with its access to the substrate regulated by structural fluctuation of the enzyme. Dihydrofolate reductase (DHFR) catalyzes the NADPH-dependent reduction of dihydrofolate (DHF) to tetrahydrofolate (THF). An important step in the mechanism involves proton donation to the N5 atom of DHF. The inability to determine the protonation states of active site residues and substrate has led to a lack of consensus regarding the catalytic mechanism involved. To resolve this ambiguity, we conducted neutron and ultrahigh-resolution X-ray crystallographic studies of the pseudo-Michaelis ternary complex of Escherichia coli DHFR with folate and NADP+. The neutron data were collected to 2.0-Å resolution using a 3.6-mm3 crystal with the quasi-Laue technique. The structure reveals that the N3 atom of folate is protonated, whereas Asp27 is negatively charged. Previous mechanisms have proposed a keto-to-enol tautomerization of the substrate to facilitate protonation of the N5 atom. The structure supports the existence of the keto tautomer owing to protonation of the N3 atom, suggesting that tautomerization is unnecessary for catalysis. In the 1.05-Å resolution X-ray structure of the ternary complex, conformational disorder of the Met20 side chain is coupled to electron density for a partially occupied water within hydrogen-bonding distance of the N5 atom of folate; this suggests direct protonation of substrate by solvent. We propose a catalytic mechanism for DHFR that involves stabilization of the keto tautomer of the substrate, elevation of the pKa value of the N5 atom of DHF by Asp27, and protonation of N5 by water that gains access to the active site through fluctuation of the Met20 side chain even though the Met20 loop is closed.

Mechanistic Studies of R67 Dihydrofolate Reductase EFFECTS OF pH AND AN H62C MUTATION

R67 dihydrofolate reductase (DHFR) is encoded by an R-plasmid, and expression of this enzyme in bacteria confers resistance to the antibacterial drug, trimethoprim. This DHFR variant is not homologous in either sequence or structure with chromosomal DHFRs. The crystal structure of tetrameric R67 DHFR indicates a single active site pore that traverses the length of the molecule (Narayana, N., Matthews, D. A., Howell, E. E., and Xuong, N.-H. (1995) Nat. Struct. Biol. 2, 1018-1025). A pH profile of enzyme activity in R67 DHFR displays an acidic pKa that is protein concentration-dependent. This pKa describes dissociation of active tetramer into two relatively inactive dimers upon protonation of His-62 and the symmetry-related His-162, His-262, and His-362 residues at the dimer-dimer interfaces. Construction of an H62C mutation results in stabilization of the active tetramer via disulfide bond formation at the dimer-dimer interfaces. The oxidized, tetrameric form of H62C R67 DHFR is quite active at pH 7, and a pH profile displays increasing activity at low pH. These results indicate protonated dihydrofolate (pKa = 2.59) is the productive substrate and that R67 DHFR does not possess a proton donor.

Isothermal Titration Calorimetry for Measuring Macromolecule-Ligand Affinity

Isothermal titration calorimetry (ITC) is a useful tool for understanding the complete thermodynamic picture of a binding reaction. In biological sciences, macromolecular interactions are essential in understanding the machinery of the cell. Experimental conditions, such as buffer and temperature, can be tailored to the particular binding system being studied. However, careful planning is needed since certain ligand and macromolecule concentration ranges are necessary to obtain useful data. Concentrations of the macromolecule and ligand need to be accurately determined for reliable results. Care also needs to be taken when preparing the samples as impurities can significantly affect the experiment. When ITC experiments, along with controls, are performed properly, useful binding information, such as the stoichiometry, affinity and enthalpy, are obtained. By running additional experiments under different buffer or temperature conditions, more detailed information can be obtained about the system. A protocol for the basic setup of an ITC experiment is given.

A Balancing Act between Net Uptake of Water during Dihydrofolate Binding and Net Release of Water upon NADPH Binding in R67 Dihydrofolate Reductase

R67 dihydrofolate reductase (DHFR) catalyzes the reduction of dihydrofolate (DHF) to tetrahydrofolate using NADPH as a cofactor. This enzyme is a homotetramer possessing 222 symmetry, and a single active site pore traverses the length of the protein. A promiscuous binding surface can accommodate either DHF or NADPH, thus two nonproductive complexes can form (2NADPH or 2DHF) as well as a productive complex (NADPH·DHF). The role of water in binding was monitored using a number of different osmolytes. From isothermal titration calorimetry (ITC) studies, binding of NADPH is accompanied by the net release of 38 water molecules. In contrast, from both steady state kinetics and ITC studies, binding of DHF is accompanied by the net uptake of water. Although different osmolytes have similar effects on NADPH binding, variable results are observed when DHF binding is probed. Sensitivity to water activity can also be probed by an in vivo selection using the antibacterial drug, trimethoprim, where the water content of the media is decreased by increasing concentrations of sorbitol. The ability of wild type and mutant clones of R67 DHFR to allow host Escherichia coli to grow in the presence of trimethoprim plus added sorbitol parallels the catalytic efficiency of the DHFR clones, indicating water content strongly correlates with the in vivo function of R67 DHFR.

Tuning of the H-transfer coordinate in primitive versus well-evolved enzymes.

The nature of an H-transfer reaction catalyzed by a primitive enzyme is examined and compared to the same reaction catalyzed by a mature (highly evolved) enzyme. The findings are evaluated using two different theoretical models. The tunneling correction model[1–3] suggests that the reaction catalyzed by the mature enzyme involves extensive tunneling, while that of the primitive enzyme involves no tunneling contribution. Marcus-like models,[2–5] on the other hand, suggest that the reaction catalyzed by the primitive enzyme has a poorly reorganized reaction coordinate, while the mature enzyme has tuned the reaction coordinate to near perfect reorganization. The latter interpretation does not indicate the degree of tunneling, but it does address the level of system preparation that brings the reaction coordinate to the tunneling conformation. Importantly, the findings indicate that, in contrast to the primitive enzyme, the mature one has evolved to catalyze a reaction with asignificant tunneling contribution or with a perfectly reorganized reaction coordinate for H-tunneling (using tunneling-correction or the Marcus-like models, respectively).

Coated tube enzyme immunoassay: factors affecting sensitivity and effects of reversible protein binding to polystyrene.

Coated tube enzyme immunoassay using alkaline phosphatase conjugated to rabbit (anti-human IgG) antiserum was studied to determine conditions of maximum sensitivity. The competitive binding assay utilized showed a large increase in sensitivity with immobilized antigen levels below the levels giving rise to the maximum in the coating-antigen dilution series. The effects of reversible antigen binding to the solid phase were investigated by comparison of untreated polystyrene tubes, polystyrene tubes treated with glutaraldehyde and glass tubes activated with an aminosilane. The use of glutaraldehyde treated tubes reduced, and the use of activated glass tubes prevented the time dependent release of immobilized antigen seen with the untreated polystyrene tubes. By comparison of these solid phases, it is shown that reversible antigen immobilized in a competitive binding assay gives rise to poorer conjugate binding (three-fold), and poorer sensitivity (six-fold). A noncompetitive response was found to occur at high free antibody levels and low competing antigen concentrations. This binding behavior is moderated by the minimization of the reversible antigen immobilization.

Determining anion-quadrupole interactions among protein, DNA, and ligand molecules

Background An extensive search through the Protein Databank (about 4500 nonredundant structures) was previously completed within our lab to analyze the energetic and geometric characteristics of an understudied molecular interaction known as an anion-quadrupole (AQ) interaction. Such an interaction occurs when the positively charged edge of an aromatic ring, resulting from a quadruple moment (i.e., a dual dipole moment), renders the aromatic molecule noncovalently bound to a nearby anionic molecule. The study considered a very limited scenario of molecules that can participate in AQ interactions, consisting of the phenyl group of a phenylalanine (phe) amino acid as the aromatic participant and the carboxylate group of an aspartate (asp) or glutamate (glu) amino acid as the anionic participant. The results revealed anion-quadrupole pairs to be prevalent within most of the protein structures. It was also observed that the interaction energy for AQ pairs was heavily dependent on the angle between the anion and plane of the aromatic ring, favoring a more planar interaction. In light of these critical observations being made from such a limited scenario, only phe-glu and phe-asp pairs and in a reduced sample set of the PDB, we are now continuing this work of identifying AQ interactions using a greatly expanded strategy. We are following these four aims: 1. Optimizing the AQ-search program to run in a semi-parallel fashion and on a large cluster of processors in order to handle larger analyses, 2. Adding to our search additional anionic participants which will include non-protein structures such as DNA and small ligands, 3. Studying a subset of the AQ pairs with molecular dynamics simulations in buried and solvent exposed environments to observe non-static behavioral traits as well as the reproducibility of AQ interactions by force field parameters. 4. Building an online database for public access to our data and search program.