Partitioning roles of side chains in affinity, orientation, and catalysis with structures for mutant complexes: asparagine-229 in thymidylate synthase

Evolution of the genetic code by incorporation of amino acids that improved or changed protein function

Fifty years have passed since the genetic code was deciphered, but how the genetic code came into being has not been satisfactorily addressed. It is now widely accepted that the earliest genetic code did not encode all 20 amino acids found in the universal genetic code as some amino acids have complex biosynthetic pathways and likely were not available from the environment. Therefore, the genetic code evolved as pathways for synthesis of new amino acids became available. One hypothesis proposes that early in the evolution of the genetic code four amino acids-valine, alanine, aspartic acid, and glycine-were coded by GNC codons (N = any base) with the remaining codons being nonsense codons. The other sixteen amino acids were subsequently added to the genetic code by changing nonsense codons into sense codons for these amino acids. Improvement in protein function is presumed to be the driving force behind the evolution of the code, but how improved function was achieved by adding amino acids has not been examined. Based on an analysis of amino acid function in proteins, an evolutionary mechanism for expansion of the genetic code is described in which individual coded amino acids were replaced by new amino acids that used nonsense codons differing by one base change from the sense codons previously used. The improved or altered protein function afforded by the changes in amino acid function provided the selective advantage underlying the expansion of the genetic code. Analysis of amino acid properties and functions explains why amino acids are found in their respective positions in the genetic code.

Mapping of ligand-binding cavities in proteins

The complex interactions between proteins and small organic molecules (ligands) are intensively studied because they play key roles in biological processes and drug activities. Here, we present a novel approach to characterize and map the ligand-binding cavities of proteins without direct geometric comparison of structures, based on Principal Component Analysis of cavity properties (related mainly to size, polarity, and charge). This approach can provide valuable information on the similarities and dissimilarities, of binding cavities due to mutations, between-species differences and flexibility upon ligand-binding. The presented results show that information on ligand-binding cavity variations can complement information on protein similarity obtained from sequence comparisons. The predictive aspect of the method is exemplified by successful predictions of serine proteases that were not included in the model construction. The presented strategy to compare ligand-binding cavities of related and unrelated proteins has many potential applications within protein and medicinal chemistry, for example in the characterization and mapping of "orphan structures", selection of protein structures for docking studies in structure-based design, and identification of proteins for selectivity screens in drug design programs.

Mechanism of influence of phosphorylation on serine 124 on a decrease of catalytic activity of human thymidylate synthase

Regulation by phosphorylation is a well-established mechanism for controlling biological activity of proteins. Recently, phosphorylation of serine 124 in human thymidylate synthase (hTS) has been shown to lower the catalytic activity of the enzyme. To clarify a possible mechanism of the observed influence, molecular dynamics (MD), essential dynamics (ED) and MM-GBSA studies were undertaken. Structures derived from the MD trajectories reveal incorrect binding alignment between the pyrimidine ring of the substrate, dUMP, and the pterine ring of the cofactor analogue, THF, in the active site of the phosphorylated enzyme. The ED analysis indicates changes in the behavior of collective motions in the phosphorylated enzyme, suggesting that the formation of the closed ternary complex is hindered. Computed free energies, in agreement with structural analysis, predict that the binding of dUMP and THF to hTS is favored in the native compared to phosphorylated state of the enzyme. The paper describes at the structural level how phosphorylation at the distant site influences the ligand binding. We propose that the 'phosphorylation effect' is transmitted from the outside loop of Ser 124 into the active site via a subtle mechanism initiated by the long-range electrostatic repulsion between the phosphate groups of dUMP and Ser124. The mechanism can be described in terms of the interplay between the two groups of amino acids: the link (residues 125-134) and the patch (residues 189-192), resulting in the change of orientation of the pyrimidine ring of dUMP, which, in turn, prevents the correct alignment between the latter ring and the pterin ring of THF.

A molecular modeling study of the interaction of 2'-fluoro-substituted analogues of dUMP/FdUMP with thymidylate synthase

Molecular dynamics simulations and free energy calculations are presented, exploring previously described experimentally studied interactions of a series of 2'-fluoro-substituted dUMP/FdUMP analogues with thymidylate synthase (TS). The results show the inhibitory behaviors of 2'-F-ara-UMP, 2',2''-diF-dUMP and 2',5-diF-ara-UMP to be dependent upon the binding positions and orientations adopted by the molecules of these compounds in the active site of TS. The binding mode of 2',5-diF-ara-UMP suggests a novel role of the active site residue Trp 80, stabilizing through hydrophobic stacking the binding position of the pyrimidine ring in 2',5-diF-ara-UMP.

The effect of 5-substitution in the pyrimidine ring of dUMP on the interaction with thymidylate synthase: molecular modeling and QSAR

Thymidylate synthase (TS) is a target enzyme for a number of anticancer agents including the 5-fluorouracil metabolite, FdUMP. The present paper reports on molecular modeling studies of the effect of substitution at C(5) position in the pyrimidine ring of the TS substrate, dUMP, on the binding affinity for the enzyme. The results of molecular dynamics simulations show that the binding of C(5) analogues of dUMP to TS in the binary complexes does not undergo changes, unless a substituent with a large steric effect, such as the propyl group, is involved. On the other hand, apparent differences in the binding of the TS cofactor, resulting from varying substitution at dUMP C(5), are observed in the modeled structures of the ternary complexes of TS. These binding characteristics are supplemented with a classical QSAR model quantifying the relation between the affinity for TS and the substituent electronic and steric effects of C(5) analogues of dUMP. Based on the findings from the present work, the perspectives for finding promising new C(5) analogues of dUMP as potential agents targeted against TS are discussed.

Two active site asparagines are essential for the reaction mechanism of the class III anaerobic ribonucleotide reductase from bacteriophage T4

Class III ribonucleotide reductase is an anaerobic enzyme that uses a glycyl radical to catalyze the reduction of ribonucleotides to deoxyribonucleotides and formate as ultimate reductant. The reaction mechanism of class III ribonucleotide reductases requires two cysteines within the active site, Cys-79 and Cys-290 in bacteriophage T4 NrdD numbering. Cys-290 is believed to form a transient thiyl radical that initiates the reaction with substrate and Cys-79 to take part as a transient thiyl radical in later steps of the reductive reaction. The recently solved three-dimensional structure of class III ribonucleotide reductase (RNR) from bacteriophage T4 shows that two highly conserved asparagines, Asn-78 and Asn-311, are positioned close to the essential Cys-79. We have investigated the function of Asn-78 and Asn-311 by site-directed mutagenesis and measured enzyme activity and glycyl radical formation in five single (N78(A/C/D) and N311(A/C)) and one double (N78A/N311A) mutant proteins. Our results suggest that both asparagines are important for the catalytic mechanism of class III RNR and that one asparagine can partially compensate for the lack of the other functional group in the single Asn --> Ala mutant proteins. A plausible role for these two asparagines could be in positioning formate in the active site to orient it toward the proposed thiyl radical of Cys-79. This would also control the highly reactive carbon dioxide radical anion form of formate within the active site before it is released as carbon dioxide. A detailed reaction scheme including the function of the two asparagines and two formate molecules is proposed for class III RNRs.

Crystal structures of 5-fluoro-dUrd and its 2 and/or 4-thio analogues: models of substituted dUMP pyrimidine ring interacting with thymidylate synthase

In order to understand the influence on thymidylate synthase interactions with dUMP analogues of the pyrimidine ring 2- and/or 4-thio, and 5-fluoro substitutions, X-ray diffractions by crystals of 5-fluoro-dUrd and its 2- and 4-thio, and 2,4-dithio analogues were measured, the four structures solved and refined. The following conclusions were suggested by results of comparative analyses of structural parameters (bond lengths, valence angles), followed by theoretical considerations based on calculated resonance structure distributions and aromaticity indices of the uracil, thiouracil, fluorouracil and fluorothiouracil rings. The effect of 4-thio substitution of FdUMP, altering specificity of inactivation of thymidylate synthases from various sources, is probably due to weaker proton acceptor power of the 4-thio substituent and increasing acidity (enhanced proton-donor power) of the N(3)-H moiety, resulting in an impaired fitness into the network of hydrogen bonds in the enzyme active center cleft. 2,4-Dithio substitution results in (i) impaired pyrimidine ring recognition by the enzyme active center, due to the 4-thio substituent (ii) increased pyrimidine ring aromaticity in dUMP, leading to resistance of C(6) to nucleophilic attack by the enzyme active center cysteine and (iii) altered planarity of the pyrimidine ring and deflections, with respect to the ring plane, of substituents at C(2), C(4) and C(5). 5-Fluoro substitution apparently activates the pyrimidine ring towards the interaction with thymidylate synthase by producing local strain, which results in an increased reactivity as predicted by the Walsh-Bent rule.

Thymidylate synthase inhibition: a structure-based rationale for drug design

Thymidylate synthase (TS) is a very interesting target in antiproliferative diseases. Its inhibition causes thimineless death of the cells and compounds inhibiting TS are widely used in anticancer therapy. The classical antifolate TS inhibitors are structural analogs of the folate cofactor; they often share the same metabolic pathways and this causes the development of resistance inside the cells. A detailed analysis of the available x-ray crystal structures of the complexes of the enzyme with different substrates and inhibitors support the finding of a structural basis of their biological activity. TS inhibitors nonstructural analog of folate, non-analog antifolate inhibitors (NAAI), are welcome as a new interesting research topic. Among the most recent and interesting ones, compounds from Agouron related to the indole structure, are independent on the folate metabolism, highly active and specific for human TS. Other compounds, phthalein derivatives, can inhibit TS enzymes from various sources and show an interesting biological activity profile: they inhibit better bacterial and fungal TS than human TS. The x-ray crystal structures of some of these inhibitors with TS show that they bind in a different binding site from that of the classical folate TS inhibitors. This indicates a potential allosteric binding site useful for future drug discovery studies.

Thermodynamic characterization of the binding of dCMP to the Asn229Asp mutant of thymidylate synthase

Isothermal titration microcalorimetry and equilibrium dialysis have been used to characterize the binding of 2'-deoxycytidine 5'-monophosphate (dCMP) to the Asn229Asp mutant of Lactobacillus casei recombinant thymidylate synthase at pH 7.4 over a temperature range of 15 degrees C to 35 degrees C. Equilibrium dialysis analysis shows that dCMP binds to two sites in the dimer of both wild-type and mutant thymidylate synthase. A concomitant net uptake of protons with binding of dCMP to both enzymes, was detected carrying out calorimetric experiments in various buffer systems with different heats of ionization. The change in protonation for binding of dCMP to wild-type enzyme is lower than that obtained for binding of this nucleotide to TS N229D, which suggests that the pK value of Asp-229 is increased upon dCMP binding to the mutant enzyme. At 25 degrees C, although the binding of dCMP to wild-type and N229D TS is favoured by both enthalpy and entropy changes, the enthalpy change is more negative for the mutant protein. Thus, the substitution of Asn 229 for Asp results in a higher affinity of TS for dCMP due to a more favourable enthalpic contribution. The Gibbs energy change of binding of dCMP to the mutant enzyme is weakly temperature-dependent, because of the enthalpy-entropy compensation arising from a negative heat capacity change of binding equal to -0.83 +/- 0.02 kJ K(-1) per mol of dCMP bound.