Acceleration of catalysis in dihydrofolate reductase by transient, site-specific photothermal excitation.

We have studied the role of protein dynamics in chemical catalysis in the enzyme dihydrofolate reductase (DHFR), using a pump-probe method that employs pulsed-laser photothermal heating of a gold nanoparticle (AuNP) to directly excite a local region of the protein structure and transient absorbance to probe the effect on enzyme activity. Enzyme activity is accelerated by pulsed-laser excitation when the AuNP is attached close to a network of coupled motions in DHFR (on the FG loop, containing residues 116-132, or on a nearby alpha helix). No rate acceleration is observed when the AuNP is attached away from the network (distal mutant and His-tagged mutant) with pulsed excitation, or for any attachment site with continuous wave excitation. We interpret these results within an energy landscape model in which transient, site-specific addition of energy to the enzyme speeds up the search for reactive conformations by activating motions that facilitate this search.

Multicomplex-based pharmacophore modeling in conjunction with multi-target docking and molecular dynamics simulations for the identification of PfDHFR inhibitors.

Plasmodium falciparum dihydrofolate reductase enzyme (PfDHFR) is counted as one of the attractive and validated antimalarial drug targets. However, the point mutations in the active site of wild-type PfDHFR have developed resistance against the well-known antifolates. Therefore, there is a dire need for the development of inhibitors that can inhibit both wild-type and mutant-type DHFR enzyme. In the present contribution, we have constructed the common feature pharmacophore models from the available PfDHFR. A representative hypothesis was prioritized and then employed for the screening of natural product library to search for the molecules with complementary features responsible for the inhibition. The screened candidates were processed via drug-likeness filters and molecular docking studies. The docking was carried out on the wild-type PfDHFR (3QGT); double-mutant PfDHFR (3UM5 and 1J3J) and quadruple-mutant PfDHFR (1J3K) enzymes. A total of eight common hits were obtained from the docking calculations that could be the potential inhibitors for both wild and mutant type DHFR enzymes. Eventually, the stability of these candidates with the selected proteins was evaluated via molecular dynamics simulations. Except for SPECS14, all the prioritized candidates were found to be stable throughout the simulation run. Overall, the strategy employed in the present work resulted in the retrieval of seven candidates that may show inhibitory activity against PfDHFR and could be further exploited as a scaffold to develop novel antimalarials. Communicated by Ramaswamy H. Sarma.

Thermal stabilization of dihydrofolate reductase using monte carlo unfolding simulations and its functional consequences.

Design of proteins with desired thermal properties is important for scientific and biotechnological applications. Here we developed a theoretical approach to predict the effect of mutations on protein stability from non-equilibrium unfolding simulations. We establish a relative measure based on apparent simulated melting temperatures that is independent of simulation length and, under certain assumptions, proportional to equilibrium stability, and we justify this theoretical development with extensive simulations and experimental data. Using our new method based on all-atom Monte-Carlo unfolding simulations, we carried out a saturating mutagenesis of Dihydrofolate Reductase (DHFR), a key target of antibiotics and chemotherapeutic drugs. The method predicted more than 500 stabilizing mutations, several of which were selected for detailed computational and experimental analysis. We find a highly significant correlation of r=0.65-0.68 between predicted and experimentally determined melting temperatures and unfolding denaturant concentrations for WT DHFR and 42 mutants. The correlation between energy of the native state and experimental denaturation temperature was much weaker, indicating the important role of entropy in protein stability. The most stabilizing point mutation was D27F, which is located in the active site of the protein, rendering it inactive. However for the rest of mutations outside of the active site we observed a weak yet statistically significant positive correlation between thermal stability and catalytic activity indicating the lack of a stability-activity tradeoff for DHFR. By combining stabilizing mutations predicted by our method, we created a highly stable catalytically active E. coli DHFR mutant with measured denaturation temperature 7.2°C higher than WT. Prediction results for DHFR and several other proteins indicate that computational approaches based on unfolding simulations are useful as a general technique to discover stabilizing mutations.

Protein dynamics and stability: the distribution of atomic fluctuations in thermophilic and mesophilic dihydrofolate reductase derived using elastic incoherent neutron scattering.

The temperature dependence of the dynamics of mesophilic and thermophilic dihydrofolate reductase is examined using elastic incoherent neutron scattering. It is demonstrated that the distribution of atomic displacement amplitudes can be derived from the elastic scattering data by assuming a (Weibull) functional form that resembles distributions seen in molecular dynamics simulations. The thermophilic enzyme has a significantly broader distribution than its mesophilic counterpart. Furthermore, although the rate of increase with temperature of the atomic mean-square displacements extracted from the dynamic structure factor is found to be comparable for both enzymes, the amplitudes are found to be slightly larger for the thermophilic enzyme. Therefore, these results imply that the thermophilic enzyme is the more flexible of the two.

X-ray crystallographic studies of Candida albicans dihydrofolate reductase. High resolution structures of the holoenzyme and an inhibited ternary complex.

The recent rise in systemic fungal infections has created a need for the development of new antifungal agents. As part of an effort to provide therapeutically effective inhibitors of fungal dihydrofolate reductase (DHFR), we have cloned, expressed, purified, crystallized, and determined the three-dimensional structure of Candida albicans DHFR. The 192-residue enzyme, which was expressed in Escherichia coli and purified by methotrexate affinity and cation exchange chromatography, was 27% identical to human DHFR. Crystals of C. albicans DHFR were grown as the holoenzyme complex and as a ternary complex containing a pyrroloquinazoline inhibitor. Both complexes crystallized with two molecules in the asymmetric unit in space group P21. The final structures had R-factors of 0.199 at 1.85-A resolution and 0.155 at 1.60-A resolution, respectively. The enzyme fold was similar to that of bacterial and vertebrate DHFR, and the binding of a nonselective diaminopyrroloquinazoline inhibitor and the interactions of NADPH with protein were typical of ligand binding to other DHFRs. However, the width of the active site cleft of C. albicans DHFR was significantly larger than that of the human enzyme, providing a basis for the design of potentially selective inhibitors.

Conditions of forming protein complexes with GroEL can influence the mechanism of chaperonin-assisted refolding.

The interaction of GroEL with urea-unfolded dihydrofolate reductase (DHFR) has been studied in the presence of DHFR substrates by investigating the ability of GroES to release enzyme under conditions where a stable GroES-GroEL-DHFR ternary complex can be formed. In these circumstances, GroES could only partially discharge the DHFR if ADP was present in the solution and approximately half of the DHFR remained bound on the chaperonin. This bound DHFR could be rescued by addition of ATP and KCl into the refolding mixture. The stable ternary complex did not show any significant protection of bound DHFR against proteolysis by Proteinase K. These results are in contrast to those observed with the GroEL-DHFR complex formed by thermal inactivation of DHFR at 45 degrees C in which GroES addition leads to partial protection of bound DHFR. Thus, the method of presentation influences the properties of the bound intermediates. It is suggested that the ability of GroES to bind on the same side of the GroEL double toroid as the target protein and displace it into the central cavity depends on the way the protein-substrate is presented to the GroEL molecule. Therefore, the compact folding intermediate formed by thermal unfolding can be protected against proteolysis after GroES binds to form a ternary complex. In addition, structural changes within GroEL induced by the experimental conditions may contribute to differences in the properties of the complexes. The more open urea-unfolded DHFR binds on the surface of chaperonin and can be displaced into solution by the tighter binding GroES molecule. It is suggested that the state of the unfolded protein when it is presented to GroEL determines the detailed mechanism of its assisted refolding. It follows that individual proteins, having characteristic folding intermediates, can have different detailed mechanisms of chaperonin-assisted folding.

Comparison of protein surfaces using a genetic algorithm.

A genetic algorithm (GA) is described which is used to compare the solvent-accessible surfaces of two proteins or fragments of proteins, represented by a dot surface calculated using the Connolly algorithm. The GA is used to move one surface relative to the other to locate the most similar surface region between the two. The matching process is enhanced by the use of the surface normals and shape terms provided by the Connolly program and also by a simple hydrogen-bonding descriptor and an additional shape descriptor. The algorithm has been tested in applications ranging from the comparison of small surface patches to the comparison of whole protein surfaces, and it has performed correctly in all cases. Examples of the matches are given and a quantitative analysis of the quality of the matches is performed. A number of possible future enhancements to the program are described which would allow the GA to be used for more complex surface comparisons.

Development of nonpolar surfaces in the folding of Escherichia coli dihydrofolate reductase detected by 1-anilinonaphthalene-8-sulfonate binding.

The development of nonpolar surfaces during the folding of Escherichia coli dihydrofolate reductase (DHFR) was studied by monitoring the time-dependent fluorescence of 1-anilinonaphthalene-8-sulfonate (ANS) included in the refolding solution. Stopped-flow refolding experiments demonstrated a rapid increase in fluorescence intensity within the dead time of mixing (5 ms), indicating that the earliest detectable folding intermediate contains hydrophobic surfaces which are capable of binding ANS. A further increase in fluorescence intensity over the next 300 ms coincides with the formation of a set of four intermediates which are known to contain a specific tertiary contact [Kuwajima, K., Garvey, E. P., Finn, B. E., Matthews, C. R., & Sugai, S. (1991) Biochemistry 30, 7693-7703]. Experiments performed in the presence of polar fluorescence quenching agents indicate that the binding sites for ANS in the burst phase species are more exposed to solvent than those in the subsequent set of intermediates. When considered along with the above study of the formation of secondary structure by stopped-flow circular dichroism, these results imply that DHFR initially forms a molten globule intermediate. Subdomains containing specific tertiary structure and more solvent-excluded ANS binding sites then form before ultimately being converted to native or native-like conformations during the rate-limiting steps in the folding of DHFR. The occurrence of similar kinetic phases observed by ANS binding during the folding of a number of other proteins suggests that this may be a common scheme for protein folding reactions.

A polypeptide bound by the chaperonin groEL is localized within a central cavity.

Chaperonins are oligomeric protein complexes that play an essential role in the cell, mediating ATP-dependent polypeptide chain folding in a variety of cellular compartments. They appear to bind early folding intermediates, preventing their aggregation; in the presence of MgATP and a cochaperonin, bound polypeptides are released in a stepwise manner, associated with folding to the native state. Chaperonin complexes appear in the electron microscope as cylindrical structures, usually composed of two stacked rings, each containing, by negative staining, an electron dense central "hole" approximately 6.0 nm in diameter. We sought to identify the site on the Escherichia coli chaperonin groEL, where the "molten globule"-like intermediate of dihydrofolate reductase (DHFR) becomes bound, by examining in the scanning transmission electron microscope complexes formed between groEL and DHFR molecules bearing covalently crosslinked 1.4-nm gold clusters. In top views of the groEL complexes, gold densities were observed in the central region; in side views, the densities were seen at the end portions of the cylinders, corresponding to positions within the individual rings. In some cases, two gold densities were observed in the same groEL complex. We conclude that folding intermediates are bound inside central cavities within individual chaperonin rings. In this potentially sequestered location, folding intermediates with a compact conformation can be bound at multiple sites by surrounding monomeric members of the ring; localization of folding within the cavity could also facilitate rebinding of structures that initially fail to incorporate properly into the folding protein.