Effects of point mutations at the flexible loop alanine-145 of Escherichia coli dihydrofolate reductase on its stability and function

Exploring the potential and identifying Withania somnifera alkaloids as novel dihydrofolate reductase (DHFR) inhibitors by the AlteQ method

There is an urgent need to discover and develop novel drugs to combat Mycobacterium tuberculosis, the causative agent of tuberculosis (TB) in humans. Alkaloids have been shown to have wide-ranging therapeutic application and could be ideal candidates for drug development, and research is underway to develop new anti-tubercular drugs from natural sources. In this regard, the current research deals with finding novel lead compounds from the Withania somnifera (WS) plant. Broad health benefits of WS are due to the presence of diverse chemical constituents which include anaferine and anahygrine and which belong to the alkaloid family. In the present study, these two compounds have been theoretically studied to understand their electronic properties using the density functional theory (DFT) at the B3LYP/6-311 + G (d,p) level. HOMO and LUMO properties and molecular electrostatic potential (MEP) surface were calculated. Further, to understand the mechanism of action of these compounds and to identify their putative drug target, molecular docking and dynamics studies were employed against Mycobacterium tuberculosis dihydrofolate reductase (DHFR). It was determined that NADP+ affects stability of the complexes by reducing fluctuations of residues 14-23 and 117-126. It was also found that Ile5 and Gln28 play an important role in complexation. Electron density analysis (using the AlteQ method) of the intermolecular region, analyzing both the anaferin-NADP+ and anahygrin-NADP+ complexes showed that anaferin and anahygrin complexes are more stable in the presence of NADP+. It has been established that in most intermolecular contacts the contribution of the ligand to the electron density is greater than that of NADP+. The present study thus provides an excellent way to analyze the effect of anaferine and anahygrine in essential processes of M. tuberculosis.Communicated by Ramaswamy H. Sarma.

FtsH degrades dihydrofolate reductase by recognizing a partially folded species

AAA+ proteolytic machines play essential roles in maintaining and rebalancing the cellular proteome in response to stress, developmental cues, and environmental changes. Of the five AAA+ proteases in Escherichia coli, FtsH is unique in its attachment to the inner membrane and its function in degrading both membrane and cytosolic proteins. E. coli dihydrofolate reductase (DHFR) is a stable and biophysically well-characterized protein, which a previous study found resisted FtsH degradation despite the presence of an ssrA degron. By contrast, we find that FtsH degrades DHFR fused to a long peptide linker and ssrA tag. Surprisingly, we also find that FtsH degrades DHFR with shorter linkers and ssrA tag, and without any linker or tag. Thus, FtsH must be able to recognize a sequence element or elements within DHFR. We find that FtsH degradation of DHFR is noncanonical in the sense that it does not rely upon recognition of an unstructured polypeptide at or near the N-terminus or C-terminus of the substrate. Results using peptide-array experiments, mutant DHFR proteins, and fusion proteins suggest that FtsH recognizes an internal sequence in a species of DHFR that is partially unfolded. Overall, our findings provide insight into substrate recognition by FtsH and indicate that its degradation capacity is broader than previously reported.

Volume and Compressibility of Proteins

The partial specific (or molar) volume, expansibility, and compressibility of a protein are fundamental thermodynamic quantities for characterizing its structure in solution. We review the definitions, measurements, and implications of these volumetric quantities in relation to protein structural biology. The partial specific volumes under constant molality (isomolal) and chemical potential (isopotential) conditions of the cosolvent (multicomponent systems) are explained in terms of preferential solvent interactions relevant to the solubility and stability of proteins. The partial expansibility is briefly discussed in terms of the effects of temperature on protein-solvent interactions (hydration) and internal packing defects (cavities). We discuss the compressibility-structure-function relationships of proteins based on analyses of the correlations between the partial adiabatic compressibilities and the structures or functions of various globular proteins (including mutants), focusing on the roles of the internal cavities in structural fluctuations. The volume and compressibility changes associated with various conformational transitions are also discussed in terms of the changes in hydration and cavities in order to elucidate the nonnative structures and the transition mechanisms, especially those associated with pressure denaturation.

Vibrational softening of a protein on ligand binding

Neutron scattering experiments have demonstrated that binding of the cancer drug methotrexate softens the low-frequency vibrations of its target protein, dihydrofolate reductase (DHFR). Here, this softening is fully reproduced using atomic detail normal-mode analysis. Decomposition of the vibrational density of states demonstrates that the largest contributions arise from structural elements of DHFR critical to stability and function. Mode-projection analysis reveals an increase of the breathing-like character of the affected vibrational modes consistent with the experimentally observed increased adiabatic compressibility of the protein on complexation.

Coupling effects of distal loops on structural stability and enzymatic activity of Escherichia coli dihydrofolate reductase revealed by deletion mutants

Residues distal from the active site in dihydrofolate reductase (DHFR) have regulatory roles in catalytic reaction and also folding stability. The couplings of the distal residues to the ones in the active site have been analyzed using site-directed mutants. To expand our understanding of the structural and functional influences of distal residue mutation, we explored the structural stability and enzymatic activity of deletion mutants. Deletion has greater structural and dynamical impacts on the corresponding part than site-directed mutation does. Thus, deletion amplifies the effects caused by distal mutations, which should make the mutual couplings among the distant residues more apparent. We focused on residues 52, 67, 121, and 145 in the four distinct loops of DHFR. All the single-residue deletion mutants showed marked reduction in stability, except for Delta52 in an alphaC-betaC loop. Double deletion mutants showed that the loop alphaC-betaC has nonadditive couplings with the betaF-betaG and betaG-betaH loops regarding stability. Single deletion to the loops alphaC-betaC or betaC-betaD resulted in considerable activity reduction, demonstrating that the loops couple to the residues near the active site. The four loops were shown to be functionally interdependent from the double deletion experiments.

Effects of pressure on enzyme function of Escherichia coli dihydrofolate reductase

To elucidate the effects of pressure on the function of Escherichia coli dihydrofolate reductase (DHFR), the enzyme activity and the dissociation constants of substrates and cofactors were measured at pressures up to 250 MPa at 25 degrees C and pH 7.0. The enzyme activity decreased with increasing pressure, accompanying the activation volume of 7.8 ml mol(-1). The values of the Michaelis constant (K(m)) for dihydrofolate and NADPH were slightly higher at 200 MPa than at atmospheric pressure. The hydride-transfer step was insensitive to pressure, as monitored by the effects of the deuterium isotope of NADPH on the reaction velocity. The dissociation constants of substrates and cofactors increased with pressure, producing volume reductions from 6.5 ml mol(-1) (tetrahydrofolate) to 33.5 ml mol(-1) (NADPH). However, the changes in Gibbs free energy with dissociation of many ligands showed different pressure dependences below and above 50 MPa, suggesting conformational changes of the enzyme at high pressure. The enzyme function at high pressure is discussed based on the volume levels of the intermediates and the candidates for the rate-limiting process.

Ensemble-based signatures of energy propagation in proteins: a new view of an old phenomenon

The ability of a protein to transmit the energetic effects of binding from one site to another constitutes the underlying basis for allosterism and signal transduction. Despite clear experimental evidence indicating the ability of proteins to transmit the effects of binding, the means by which this propagation is facilitated is not well understood. Using our previously developed ensemble-based description of the equilibrium, we investigated the physical basis of energy propagation and identified several fundamental and general aspects of energetic coupling between residues in a protein. First, partitioning of a conformational ensemble into four distinct sub-ensembles allows for explanation of the range of experimentally observed coupling behaviors (i.e., positive, neutral, and negative coupling between various regions of the protein structure). Second, the relative thermodynamic properties of these four sub-ensembles define the energetic coupling between residues as either positive, neutral, or negative. Third, analysis of the structural and thermodynamic features of the states within each sub-ensemble reveals significant variability. This third result suggests that a quantitative description of energy propagation in proteins requires an understanding of the structural and energetic features of more than just one or a few low-energy states, but also of many high-energy states. Such findings illuminate the difficulty in interpreting energy propagation in proteins in terms of a structural pathway that physically links coupled sites.

Evolutional Design of a Hyperactive Cysteine- and Methionine-free Mutant of Escherichia coli Dihydrofolate Reductase

We developed a strategy for finding out the adapted variants of enzymes, and we applied it to an enzyme, dihydrofolate reductase (DHFR), in terms of its catalytic activity so that we successfully obtained several hyperactive cysteine- and methionine-free variants of DHFR in which all five methionyl and two cysteinyl residues were replaced by other amino acid residues. Among them, a variant (M1A/M16N/M20L/M42Y/C85A/M92F/C152S), named as ANLYF, has an approximately seven times higher k(cat) value than wild type DHFR. Enzyme kinetics and crystal structures of the variant were investigated for elucidating the mechanism of the hyperactivity. Steady-state and transient binding kinetics of the variant indicated that the kinetic scheme of the catalytic cycle of ANLYF was essentially the same as that of wild type, showing that the hyperactivity was brought about by an increase of the dissociation rate constants of tetrahydrofolate from the enzyme-NADPH-tetrahydrofolate ternary complex. The crystal structure of the variant, solved and refined to an R factor of 0.205 at 1.9-angstroms resolution, indicated that an increased structural flexibility of the variant and an increased size of the N-(p-aminobenzoyl)-L-glutamate binding cleft induced the increase of the dissociation constant. This was consistent with a large compressibility (volume fluctuation) of the variant. A comparison of folding kinetics between wild type and the variant showed that the folding of these two enzymes was similar to each other, suggesting that the activity enhancement of the enzyme can be attained without drastic changes of the folding mechanism.

Effects of Mutation at Methionine-42 of Escherichia coli Dihydrofolate Reductase on Stability and Function: Implication of Hydrophobic Interactions

Methionine-42, distal to the active site of Escherichia coli dihydrofolate reductase, was substituted by site-directed mutagenesis with 14 amino acids (Ala, Cys, Glu, Gln, Gly, His, Ile, Leu, Pro, Ser, Thr, Trp, Tyr, and Val) to elucidate its role in the stability and function of this enzyme. Far-ultraviolet circular dichroism spectra of these mutants showed a distinctive negative peak at around 230 nm beside 220 nm, depending on the hydrophobicity of the amino acids introduced. The fluorescence intensity also increased in an order similar to that of the amino acids. These spectroscopic data suggest that the mutations do not affect the secondary structure, but strongly perturb the exciton coupling between Trp47 and Trp74. The free energy of urea unfolding, deltaG(o)u, increased with increases in the side-chain hydrophobicity in the range 2.96-6.40 kcal x mol(-1), which includes the value for the wild-type enzyme (6.08 kcal x mol(-1)). The steady-state kinetic parameters, Km and kcat, also increased with increases in the side-chain hydrophobicity, with the M42W mutant showing the largest increases in Km (35-fold) and kcat (4.3-fold) compared with the wild-type enzyme. These results demonstrate that site 42 distal to the active site plays an important role in the stability and function of this enzyme, and that the main effect of the mutations is to modify of hydrophobic interactions with the residues surrounding this position.

Compressibility gives new insight into protein dynamics and enzyme function

The adiabatic compressibility of enzyme is largely influenced by binding of coenzyme and substrate, due to the changes in atomic packing. Amino acid substitution also induces large changes in compressibility parallel to enzyme activity. These results demonstrate that a small alteration of local structure by ligand binding and mutation is dramatically magnified in the flexibility of protein molecule to affect the function. Compressibility gives new insight into protein dynamics and enzyme function from the aspect of atomic packing or cavity which cannot be obtained by other techniques.

Molecular dynamics simulation of Escherichia coli dihydrofolate reductase and its protein fragments: relative stabilities in experiment and simulations

We have carried out molecular dynamics simulations of the native dihydrofolate reductase from Escherichia coli and several of its folded protein fragments at standard temperature. The simulations have shown fragments 1--36, 37--88, and 89--159 to be unstable, with a C(alpha)RMSD (C(alpha) root mean squared deviation) >5 A after 3.0 nsec of simulation. The unfolding of fragment 1--36 was immediate, whereas fragments 37--88 and 89--159 gradually unfolded because of the presence of the beta-sheet core structure. In the absence of residues 1--36, the two distinct domains comprising fragment 39--159 associated with each other, resulting in a stable conformation. This conformation retained most of its native structural elements. We have further simulated fragments derived from computational protein cutting. These were also found to be unstable, with the exception of fragment 104--159. In the absence of alpha(4), the loose loop region of residues 120--127 exhibited a beta-strand-like behavior, associating itself with the beta-sheet core of the protein fragment. The current study suggests that the folding of dihydrofolate reductase involves cooperative folding of distinct domains which otherwise would have been unstable as independent folded units in solution. Finally, the critical role of residues 1--36 in allowing the two distinct domains of fragment 104--159 to fold into the final native conformation is discussed.

Binding sites in Escherichia coli dihydrofolate reductase communicate by modulating the conformational ensemble

To explore how distal mutations affect binding sites and how binding sites in proteins communicate, an ensemble-based model of the native state was used to define the energetic connectivities between the different structural elements of Escherichia coli dihydrofolate reductase. Analysis of this model protein has allowed us to identify two important aspects of intramolecular communication. First, within a protein, pair-wise couplings exist that define the magnitude and extent to which mutational effects propagate from the point of origin. These pair-wise couplings can be identified from a quantity we define as the residue-specific connectivity. Second, in addition to the pair-wise energetic coupling between residues, there exists functional connectivity, which identifies energetic coupling between entire functional elements (i.e., binding sites) and the rest of the protein. Analysis of the energetic couplings provides access to the thermodynamic domain structure in dihydrofolate reductase as well as the susceptibility of the different regions of the protein to both small-scale (e.g., point mutations) and large-scale perturbations (e. g., binding ligand). The results point toward a view of allosterism and signal transduction wherein perturbations do not necessarily propagate through structure via a series of conformational distortions that extend from one active site to another. Instead, the observed behavior is a manifestation of the distribution of states in the ensemble and how the distribution is affected by the perturbation.

Effect of ligand binding on the flexibility of dihydrofolate reductase as revealed by compressibility

The partial specific volume, v, and adiabatic compressibility, beta(s), of Escherichia coli dihydrofolate reductase were measured at 30 degrees C in the presence of various ligands (folate, dihydrofolate, tetrahydrofolate, NADPH, NADP, methotrexate, and KCl). Binding of these ligands (binary and ternary complexes) brought about large changes of v (0.734-0.754 cm(3) g(-1)) and beta(s) (6. 6x10(-6)-9.8x10(-6) bar(-1)), keeping a linear relationship between the two parameters. The values of v and beta(s) increased with an increase in internal cavity, V(cav), and a decrease in accessible surface area, ASA, which were calculated from the X-ray crystal structures of the complexes. A large variation of V(cav) relative to ASA by ligand binding suggested that the cavity is a dominant factor and the effect of hydration might be small for the ligand-induced changes of v and beta(s). The beta(s) values of the binary and ternary complexes suggested a characteristic conformational flexibility of the kinetic intermediates in the enzyme reaction coordinate. Comparison of beta(s) with the cavity distribution in the crystal structures revealed that the flexibility of the intermediates was mainly determined by the total cavity volume with minor contributions of the number, position, and size of cavities. These results demonstrate that the compressibility is a useful measure of the conformational flexibility of the intermediates in the enzyme reaction and that the combined study of compressibility and X-ray crystallography gives new insight into the protein dynamics through the behavior of the cavities.

Circular permutation analysis as a method for distinction of functional elements in the M20 loop of Escherichia coli dihydrofolate reductase

A functional element of an enzyme can be defined as the smallest unit of the local peptide backbone of which the connectivity is crucial for the catalytic activity. In order to elucidate the distribution of functional elements in an active site flexible loop (the M20 loop) of Escherichia coli dihydrofolate reductase, systematic cleavage of main chain connectivity was performed using circular permutation. Our analysis is based on the assumption that a permutation within a functional element would significantly reduce enzyme function, whereas ones outside or at the boundaries of the elements would affect the function only slightly. Thus, a functional element would be assigned as the minimum peptide chain between the identified boundaries. Comparison of the activities of the circularly permuted variants revealed that the peptide chain around the M20 loop could be divided into four regions (regions 1-4). Region 1 was found to play an important role in overall tertiary fold because most variants permuted at region 1 did not accumulate in E. coli cells stably. A distinction between region 2 and region 3 was in agreement with the extent of movements calculated from the coordinates of alpha carbons, supporting the idea that the movement of peptide backbone is a key feature of enzyme function. The boundary between region 3 and region 4 coincided with that between the M20 loop and the following alpha helix. From equilibrium binding studies, region 2 was found to be involved in the binding of nicotinamide substrates, whereas region 4 appeared to be very important for the binding of pterin substrates.