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

Stabilization of the cyclodecadiene derivative isofuranodiene by silver (I) coordination. Mechanistic and biological aspects

The industrial extraction and further applications of isofuranodiene are limited because at room temperature it spontaneously converts to curzerene, a structurally less active isomer. This work definitively identified the structure of isofuranodiene in the solid state, showing the two methyl groups in syn position. In addition, two bioactive metal cations, namely, silver(I) and copper(II) ions, were used in the attempt to obtain the chemical stability of isofuranodiene: in the case of silver(I), a labile adduct was formed, while in the case of copper(II), a more stable 1:1 adduct was achieved. In the former, the presence of silver did not significantly affect the biological activities of isofuranodiene, while in the latter, the copper(II) coordination suppressed them. The biological activities of the isofuranodiene adducts were then evaluated as antiproliferative agents against human tumor cell lines (HCT116, MDA-MB 231, and T98G). In addition, for the first time, isofuranodiene was tested as an inhibitor of DHFR (DiHydroFolateReductase) from Escherichia coli. Anticancer activity was observed in the isofuranodiene with the AgCF₃SO₃ adduct, in the tested cell lines, with IC₅₀ values ranging from 4.89μM to 13.06μM, while inhibition assays highlighted a Ki of 6.22μM for isofuranodiene and of 0.17μM for the related silver adduct. Docking studies indicated a binding mode score of -6.83Kcal/mol for isofuranodiene, and an energy value of -11.82Kcal/mol for methotrexate (a classic DHFR inhibitor).

Aspects of Weak Interactions between Folate and Glycine Betaine

Folate, or vitamin B9, is an important compound in one-carbon metabolism. Previous studies have found weaker binding of dihydrofolate to dihydrofolate reductase in the presence of osmolytes. In other words, osmolytes are more difficult to remove from the dihydrofolate solvation shell than water; this shifts the equilibrium toward the free ligand and protein species. This study uses vapor-pressure osmometry to explore the interaction of folate with the model osmolyte, glycine betaine. This method yields a preferential interaction potential (μ23/RT value). This value is concentration-dependent as folate dimerizes. The μ23/RT value also tracks the deprotonation of folate's N3-O4 keto-enol group, yielding a pKa of 8.1. To determine which folate atoms interact most strongly with betaine, the interaction of heterocyclic aromatic compounds (as well as other small molecules) with betaine was monitored. Using an accessible surface area approach coupled with osmometry measurements, deconvolution of the μ23/RT values into α values for atom types was achieved. This allows prediction of μ23/RT values for larger molecules such as folate. Molecular dynamics simulations of folate show a variety of structures from extended to L-shaped. These conformers possess μ23/RT values from -0.18 to 0.09 m-1, where a negative value indicates a preference for solvation by betaine and a positive value indicates a preference for water. This range of values is consistent with values observed in osmometry and solubility experiments. As the average predicted folate μ23/RT value is near zero, this indicates folate interacts almost equally well with betaine and water. Specifically, the glutamate tail prefers to interact with water, while the aromatic rings prefer betaine. In general, the more protonated species in our small molecule survey interact better with betaine as they provide a source of hydrogens (betaine is not a hydrogen bond donor). Upon deprotonation of the small molecule, the preference swings toward water interaction because of its hydrogen bond donating capacities.

Ab initio structure determination from prion nanocrystals at atomic resolution by MicroED

Electrons, because of their strong interaction with matter, produce high-resolution diffraction patterns from tiny 3D crystals only a few hundred nanometers thick in a frozen-hydrated state. This discovery offers the prospect of facile structure determination of complex biological macromolecules, which cannot be coaxed to form crystals large enough for conventional crystallography or cannot easily be produced in sufficient quantities. Two potential obstacles stand in the way. The first is a phenomenon known as dynamical scattering, in which multiple scattering events scramble the recorded electron diffraction intensities so that they are no longer informative of the crystallized molecule. The second obstacle is the lack of a proven means of de novo phase determination, as is required if the molecule crystallized is insufficiently similar to one that has been previously determined. We show with four structures of the amyloid core of the Sup35 prion protein that, if the diffraction resolution is high enough, sufficiently accurate phases can be obtained by direct methods with the cryo-EM method microelectron diffraction (MicroED), just as in X-ray diffraction. The success of these four experiments dispels the concern that dynamical scattering is an obstacle to ab initio phasing by MicroED and suggests that structures of novel macromolecules can also be determined by direct methods.

Protein rethreading: A novel approach to protein design

Protein engineering is an important tool for the design of proteins with novel and desirable features. Templates from the protein databank (PDB) are often used as initial models that can be modified to introduce new properties. We examine whether it is possible to reconnect a protein in a manner that generates a new topology yet preserves its structural integrity. Here, we describe the rethreading of dihydrofolate reductase (DHFR) from E. coli (wtDHFR). The rethreading process involved the removal of three native loops, and the introduction of three new loops with alternate connections. The structure of the rethreaded DHFR (rDHFR-1) was determined to 1.6 Å, demonstrating the success of the rethreading process. Both wtDHFR and rDHFR-1 exhibited similar affinities towards methotrexate. However, rDHFR-1 showed no reducing activity towards dihydrofolate, and exhibited about ~6-fold lower affinity towards NADPH than wtDHFR. This work demonstrates that protein rethreading can be a powerful tool for the design of a large array of proteins with novel structures and topologies, and that by careful rearrangement of a protein sequence, the sequence to structure relationship can be expanded substantially.

Long-Range Electrostatics-Induced Two-Proton Transfer Captured by Neutron Crystallography in an Enzyme Catalytic Site

Neutron crystallography was used to directly locate two protons before and after a pH-induced two-proton transfer between catalytic aspartic acid residues and the hydroxy group of the bound clinical drug darunavir, located in the catalytic site of enzyme HIV-1 protease. The two-proton transfer is triggered by electrostatic effects arising from protonation state changes of surface residues far from the active site. The mechanism and pH effect are supported by quantum mechanics/molecular mechanics (QM/MM) calculations. The low-pH proton configuration in the catalytic site is deemed critical for the catalytic action of this enzyme and may apply more generally to other aspartic proteases. Neutrons therefore represent a superb probe to obtain structural details for proton transfer reactions in biological systems at a truly atomic level.

Structural Insights into Mycobacterium tuberculosis Rv2671 Protein as a Dihydrofolate Reductase Functional Analogue Contributing to para-Aminosalicylic Acid Resistance

Mycobacterium tuberculosis (Mtb) Rv2671 is annotated as a 5-amino-6-ribitylamino-2,4(1H,3H)-pyrimidinedione 5'-phosphate (AROPP) reductase (RibD) in the riboflavin biosynthetic pathway. Recently, a strain of Mtb with a mutation in the 5' untranslated region of Rv2671, which resulted in its overexpression, was found to be resistant to dihydrofolate reductase (DHFR) inhibitors including the anti-Mtb drug para-aminosalicylic acid (PAS). In this study, a biochemical analysis of Rv2671 showed that it was able to catalyze the reduction of dihydrofolate (DHF) to tetrahydrofolate (THF), which explained why the overexpression of Rv2671 was sufficient to confer PAS resistance. We solved the structure of Rv2671 in complex with the NADP(+) and tetrahydrofolate (THF), which revealed the structural basis for the DHFR activity. The structures of Rv2671 complexed with two DHFR inhibitors, trimethoprim and trimetrexate, provided additional details of the substrate binding pocket and elucidated the differences between their inhibitory activities. Finally, Rv2671 was unable to catalyze the reduction of AROPP, which indicated that Rv2671 and its closely related orthologues are not involved in riboflavin biosynthesis.

Vacuum-Ultraviolet Circular Dichroism Spectra of Escherichia coli Dihydrofolate Reductase and Its Mutants: Contributions of Phenylalanine and Tyrosine Side Chains and Exciton Coupling of Two Tryptophan Side Chains

Vacuum-ultraviolet (VUV) circular dichroism (CD) spectroscopy has recently been used for secondary structure analysis of proteins; however, the contribution of aromatic side chains to protein VUV CD spectra is unresolved. In this report, VUV CD spectra of 10 Escherichia coli dihydrofolate reductase (DHFR) mutants, in which each phenylalanine or tyrosine residue was mutated to leucine, were measured down to 175 nm at 25 °C and pH 8.0 to elucidate the contributions of these aromatic side chains to the high-energy transitions of peptide bonds. The VUV CD spectra of these mutants were different from the spectrum of the wild-type protein, indicating that the contribution of the phenylalanine and tyrosine side chains of DHFR extends to the VUV region. Furthermore, the VUV CD spectrum and the folate- or NADP(+)-induced spectral change of F103L mutant DHFR indicated a modification and regeneration of exciton coupling between the Trp47 and Trp74 side chains, respectively, suggesting that exciton coupling may also contribute to the CD spectrum of DHFR in the VUV region. These results should be useful for theoretically characterizing the contribution of aromatic side chains to protein CD spectra, leading to the improvement of protein secondary-structure analysis by VUV CD spectroscopy.

Structural Characterization of Dihydrofolate Reductase Complexes by Top-Down Ultraviolet Photodissociation Mass Spectrometry

The stepwise reduction of dihydrofolate to tetrahydrofolate entails significant conformational changes of dihydrofolate reductase (DHFR). Binary and ternary complexes of DHFR containing cofactor NADPH, inhibitor methotrexate (MTX), or both NADPH and MTX were characterized by 193 nm ultraviolet photodissociation (UVPD) mass spectrometry. UVPD yielded over 80% sequence coverage of DHFR and resulted in production of fragment ions that revealed the interactions between DHFR and each ligand. UVPD of the binary DHFR·NADPH and DHFR·MTX complexes led to an unprecedented number of fragment ions containing either an N- or C-terminal protein fragment still bound to the ligand via retention of noncovalent interactions. In addition, holo-fragments retaining both ligands were observed upon UVPD of the ternary DHFR·NADPH·MTX complex. The combination of extensive holo and apo fragment ions allowed the locations of the NADPH and MTX ligands to be mapped, with NADPH associated with the adenosine binding domain of DHFR and MTX interacting with the loop domain. These findings are consistent with previous crystallographic evidence. Comparison of the backbone cleavage propensities for apo DHFR and its holo counterparts revealed significant variations in UVPD fragmentation in the regions expected to experience conformational changes upon binding NADPH, MTX, or both ligands. In particular, the subdomain rotation and loop movements, which are believed to occur upon formation of the transition state of the ternary complex, are reflected in the UVPD mass spectra. The UVPD spectra indicate enhanced backbone cleavages in regions that become more flexible or show suppressed backbone cleavages for those regions either shielded by the ligand or involved in new intramolecular interactions. This study corroborates the versatility of 193 nm UVPD mass spectrometry as a sensitive technique to track enzymatic cycles that involve conformational rearrangements.

Sub-atomic resolution X-ray crystallography and neutron crystallography: promise, challenges and potential

The International Year of Crystallography saw the number of macromolecular structures deposited in the Protein Data Bank cross the 100000 mark, with more than 90000 of these provided by X-ray crystallography. The number of X-ray structures determined to sub-atomic resolution (i.e. ≤1 Å) has passed 600 and this is likely to continue to grow rapidly with diffraction-limited synchrotron radiation sources such as MAX-IV (Sweden) and Sirius (Brazil) under construction. A dozen X-ray structures have been deposited to ultra-high resolution (i.e. ≤0.7 Å), for which precise electron density can be exploited to obtain charge density and provide information on the bonding character of catalytic or electron transfer sites. Although the development of neutron macromolecular crystallography over the years has been far less pronounced, and its application much less widespread, the availability of new and improved instrumentation, combined with dedicated deuteration facilities, are beginning to transform the field. Of the 83 macromolecular structures deposited with neutron diffraction data, more than half (49/83, 59%) were released since 2010. Sub-mm(3) crystals are now regularly being used for data collection, structures have been determined to atomic resolution for a few small proteins, and much larger unit-cell systems (cell edges >100 Å) are being successfully studied. While some details relating to H-atom positions are tractable with X-ray crystallography at sub-atomic resolution, the mobility of certain H atoms precludes them from being located. In addition, highly polarized H atoms and protons (H(+)) remain invisible with X-rays. Moreover, the majority of X-ray structures are determined from cryo-cooled crystals at 100 K, and, although radiation damage can be strongly controlled, especially since the advent of shutterless fast detectors, and by using limited doses and crystal translation at micro-focus beams, radiation damage can still take place. Neutron crystallography therefore remains the only approach where diffraction data can be collected at room temperature without radiation damage issues and the only approach to locate mobile or highly polarized H atoms and protons. Here a review of the current status of sub-atomic X-ray and neutron macromolecular crystallography is given and future prospects for combined approaches are outlined. New results from two metalloproteins, copper nitrite reductase and cytochrome c', are also included, which illustrate the type of information that can be obtained from sub-atomic-resolution (∼0.8 Å) X-ray structures, while also highlighting the need for complementary neutron studies that can provide details of H atoms not provided by X-ray crystallography.

Escherichia coli dihydrofolate reductase catalyzed proton and hydride transfers: temporal order and the roles of Asp27 and Tyr100

The reaction catalyzed by Escherichia coli dihydrofolate reductase (ecDHFR) has become a model for understanding enzyme catalysis, and yet several details of its mechanism are still unresolved. Specifically, the mechanism of the chemical step, the hydride transfer reaction, is not fully resolved. We found, unexpectedly, the presence of two reactive ternary complexes [enzyme:NADPH:7,8-dihydrofolate (E:NADPH:DHF)] separated by one ionization event. Furthermore, multiple kinetic isotope effect (KIE) studies revealed a stepwise mechanism in which protonation of the DHF precedes the hydride transfer from the nicotinamide cofactor (NADPH) for both reactive ternary complexes of the WT enzyme. This mechanism was supported by the pH- and temperature-independent intrinsic KIEs for the C-H→C hydride transfer between NADPH and the preprotonated DHF. Moreover, we showed that active site residues D27 and Y100 play a synergistic role in facilitating both the proton transfer and subsequent hydride transfer steps. Although D27 appears to have a greater effect on the overall rate of conversion of DHF to tetrahydrofolate, Y100 plays an important electrostatic role in modulating the pKa of the N5 of DHF to enable the preprotonation of DHF by an active site water molecule.