Significance of structural changes in proteins: expected errors in refined protein structures

Insights on Rigidity and Flexibility at the Global and Local Levels of Protein Structures and Their Roles in Homologous Psychrophilic, Mesophilic, and Thermophilic Proteins: A Computational Study

The rigidity and flexibility of homologous psychrophilic (P), mesophilic (M), and thermophilic (T) proteins have been investigated at the global and local levels in terms of "packing factors" and "atomic fluctuations" obtained from B-factors. For comparison of atomic fluctuations, correction of errors by considering errors in B-factors from all sources in a consolidated manner and conversion of the fluctuations to the same temperature have been suggested and validated. The results indicate no differences in the global values like the average packing factor among the three classes of protein homologues, but at local levels there are differences. A comparison of homologous protein triplets show that the average atomic fluctuations at a given temperature mainly obey the order P > M > T. Packing factors and the atomic fluctuations are anti-correlated, suggesting that altering the rigidity of the active site might be a potential strategy to make tailor-made psychrophilic or thermophilic proteins from their mesophilic homologues. The computer codes developed and used in this work are available at the link https://github.com/Munna-Sarkar/proteins-rigidity-flexibility.git.

Structural perturbations of substrate binding and oxidation state changes in a lytic polysaccharide monooxygenase

LPMOs are enzymes which catalyse the oxidation of a C-H bond within polysaccharides, leading to their oxidative cleavage. To achieve this, LPMOs employ highly reactive oxidising intermediates, the generation of which is likely coupled to substrate binding to the enzyme. The nature of this coupling is unknown. Here we report a statistical comparison for four three-dimensional structures of an AA9 LPMO crystallised in the same space group but in different oxidation and substrate-binding states, to determine which significant structural perturbations occur at the enzyme upon either oxidation state change or the binding of substrate. In a novel step, we determine the global random error associated with the positional coordinates of atoms using the method of moments to ascertain the statistical estimators of Gaussian distributions of pairwise RMS differences between individual atoms in different structures. The results show that a change in the oxidation state of the copper leads to no significant structural changes, and that binding of the substrate leads to a single change in the conformation of a tryptophan residue. This tryptophan has previously been identified as part of a charge transfer pathway between the active site and the external surface of the protein, and the structural change identified herein may be part of the substrate-enzyme coupling mechanism.

Dynamics-Evolution Correspondence in Protein Structures

The genotype-phenotype mapping of proteins is a fundamental question in structural biology. In this Letter, with the analysis of a large dataset of proteins from hundreds of protein families, we quantitatively demonstrate the correlations between the noise-induced protein dynamics and mutation-induced variations of native structures, indicating the dynamics-evolution correspondence of proteins. Based on the investigations of the linear responses of native proteins, the origin of such a correspondence is elucidated. It is essential that the noise- and mutation-induced deformations of the proteins are restricted on a common low-dimensional subspace, as confirmed from the data. These results suggest an evolutionary mechanism of the proteins gaining both dynamical flexibility and evolutionary structural variability.

Adeno-Associated Virus (AAV-DJ)-Cryo-EM Structure at 1.56 Å Resolution

Adeno-associated virus is the leading viral vector for gene therapy. AAV-DJ is a recombinant variant developed for tropism to the liver. The AAV-DJ structure has been determined to 1.56 Å resolution through cryo-electron microscopy (cryo-EM). Only apoferritin is reported in preprints at 1.6 Å or higher resolution, and AAV-DJ nearly matches the highest resolutions ever attained through X-ray diffraction of virus crystals. However, cryo-EM has the advantage that most of the hydrogens are clear, improving the accuracy of atomic refinement, and removing ambiguity in hydrogen bond identification. Outside of secondary structures where hydrogen bonding was predictable a priori, the networks of hydrogen bonds coming from direct observation of hydrogens and acceptor atoms are quite different from those inferred even at 2.8 Å resolution. The implications for understanding viral assembly mean that cryo-EM will likely become the favored approach for high resolution structural virology.

Diversity in kinetics correlated with structure in nano body-stabilized LacY

The structure of lactose permease, stabilized in a periplasmic open conformation by two Gly to Trp replacements (LacYww) and complexed with a nanobody directed against this conformation, provides the highest resolution structure of the symporter. The nanobody binds in a different manner than two other nanobodies made against the same mutant, which also bind to the same general region on the periplasmic side. This region of the protein may represent an immune hotspot. The CDR3 loop of the nanobody is held by hydrogen bonds in a conformation that partially blocks access to the substrate-binding site. As a result, k_on and k_off for galactoside binding to either LacY or the double mutant complexed with the nanobody are lower than for the other two LacY/nanobody complexes though the K_d values are similar, reflecting the fact that the nanobodies rigidify structures along the pathway. While the wild-type LacY/nanobody complex clearly stabilizes a similar ‘extracellular open’ conformation in solution, judged by binding kinetics, the complex with wild-type LacY did not yet crystallize, suggesting the nanobody does not bind strongly enough to shift the equilibrium to stabilize a periplasmic side-open conformation suitable for crystallization. However, the similarity of the galactoside binding kinetics for the nanobody-bound complexes with wild type LacY and with LacY_WW indicates that they have similar structures, showing that the reported co-structures reliably show nanobody interactions with LacY.

Structure determination of a major facilitator peptide transporter: Inward facing PepTSt from Streptococcus thermophilus crystallized in space group P3121

Major facilitator superfamily (MFS) peptide transporters (typically referred to as PepT, POT or PTR transporters) mediate the uptake of di- and tripeptides, and so play an important dietary role in many organisms. In recent years, a better understanding of the molecular basis for this process has emerged, which is in large part due to a steep increase in structural information. Yet, the conformational transitions underlying the transport mechanism are still not fully understood, and additional data is therefore needed. Here we report in detail the detergent screening, crystallization, experimental MIRAS phasing, and refinement of the peptide transporter PepT_St from Streptococcus thermophilus. The space group is P3₁21, and the protein is crystallized in a monomeric inward facing form. The binding site is likely to be somewhat occluded, as the lobe encompassing transmembrane helices 10 and 11 is markedly bent towards the central pore of the protein, but the extent of this potential occlusion could not be determined due to disorder at the apex of the lobe. Based on structural comparisons with the seven previously determined P2₁2₁2₁ and C222₁ structures of inward facing PepT_St, the structural flexibility as well as the conformational changes mediating transition between the inward open and inward facing occluded states are discussed. In conclusion, this report improves our understanding of the structure and conformational cycle of PepT_St, and can furthermore serve as a case study, which may aid in supporting future structure determinations of additional MFS transporters or other integral membrane proteins.

B-factor Analysis and Conformational Rearrangement of Aldose Reductase

The NADPH-dependent reduction of glucose reaction that is catalyzed by Aldose Reductase (AR) follows a sequential ordered kinetic mechanism in which the co-factor NADPH binds to the enzyme prior to the aldehyde substrate. The kinetic/structural experiments have found a conformational change involving a hinge-like movement of a surface loop (residues 213-224) which is anticipated to take place upon the binding of the diphosphate moiety of NADPH. The reorientation of this loop, expected to permit the release of NADP⁺, represents the rate-limiting step of the catalytic mechanism. This study reveals: 1) The Translation/Libration/Screw (TLS) analysis of absolute B-factors of apo AR crystal structures indicates that the 212-224 loop might move as a rigid group. 2) Residues that make the flexible loop slide in the AR binary and ternary complexes. 3) The normalized B-factors separate this segment into three different clusters with fewer residues.

Methods of protein structure comparison

Despite its apparent simplicity, the problem of quantifying the differences between two structures of the same protein or complex is nontrivial and continues evolving. In this chapter, we described several methods routinely used to compare computational models to experimental answers in several modeling assessments. The two major classes of measures, positional distance-based and contact-based, are presented, compared, and analyzed. The most popular measure of the first class, the global RMSD, is shown to be the least representative of the degree of structural similarity because it is dominated by the largest error. Several distance-dependent algorithms designed to attenuate the drawbacks of RMSD are described. Measures of the second class, contact-based, are shown to be more robust and relevant. We also illustrate the importance of using combined measures, utility-based measures, and the role of the distributions derived from the pairs of experimental structures in interpreting the results.

Determination of ensemble-average pairwise root mean-square deviation from experimental B-factors

Root mean-square deviation (RMSD) after roto-translational least-squares fitting is a measure of global structural similarity of macromolecules used commonly. On the other hand, experimental x-ray B-factors are used frequently to study local structural heterogeneity and dynamics in macromolecules by providing direct information about root mean-square fluctuations (RMSF) that can also be calculated from molecular dynamics simulations. We provide a mathematical derivation showing that, given a set of conservative assumptions, a root mean-square ensemble-average of an all-against-all distribution of pairwise RMSD for a single molecular species, <RMSD(2)>(1/2), is directly related to average B-factors (<B>) and <RMSF(2)>(1/2). We show this relationship and explore its limits of validity on a heterogeneous ensemble of structures taken from molecular dynamics simulations of villin headpiece generated using distributed-computing techniques and the Folding@Home cluster. Our results provide a basis for quantifying global structural diversity of macromolecules in crystals directly from x-ray experiments, and we show this on a large set of structures taken from the Protein Data Bank. In particular, we show that the ensemble-average pairwise backbone RMSD for a microscopic ensemble underlying a typical protein x-ray structure is approximately 1.1 A, under the assumption that the principal contribution to experimental B-factors is conformational variability.

Decision-making in structure solution using Bayesian estimates of map quality: the PHENIX AutoSol wizard

Ten measures of experimental electron-density-map quality are examined and the skewness of electron density is found to be the best indicator of actual map quality. A Bayesian approach to estimating map quality is developed and used in the PHENIX AutoSol wizard to make decisions during automated structure solution.

Estimates of the quality of experimental maps are important in many stages of structure determination of macromolecules. Map quality is defined here as the correlation between a map and the corresponding map obtained using phases from the final refined model. Here, ten different measures of experimental map quality were examined using a set of 1359 maps calculated by re-analysis of 246 solved MAD, SAD and MIR data sets. A simple Bayesian approach to estimation of map quality from one or more measures is presented. It was found that a Bayesian estimator based on the skewness of the density values in an electron-density map is the most accurate of the ten individual Bayesian estimators of map quality examined, with a correlation between estimated and actual map quality of 0.90. A combination of the skewness of electron density with the local correlation of r.m.s. density gives a further improvement in estimating map quality, with an overall correlation coefficient of 0.92. The PHENIX AutoSol wizard carries out automated structure solution based on any combination of SAD, MAD, SIR or MIR data sets. The wizard is based on tools from the PHENIX package and uses the Bayesian estimates of map quality described here to choose the highest quality solutions after experimental phasing.

The role of protein dynamics in thymidylate synthase catalysis: variants of conserved 2'-deoxyuridine 5'-monophosphate (dUMP)-binding Tyr-261

The enzyme thymidylate synthase (TS) catalyzes the reductive methylation of 2'-deoxyuridine 5'-monophosphate (dUMP) to 2'-deoxythymidine 5'-monophosphate. Using kinetic and X-ray crystallography experiments, we have examined the role of the highly conserved Tyr-261 in the catalytic mechanism of TS. While Tyr-261 is distant from the site of methyl transfer, mutants at this position show a marked decrease in enzymatic activity. Given that Tyr-261 forms a hydrogen bond with the dUMP 3'-O, we hypothesized that this interaction would be important for substrate binding, orientation, and specificity. Our results, surprisingly, show that Tyr-261 contributes little to these features of the mechanism of TS. However, the residue is part of the structural core of closed ternary complexes of TS, and conservation of the size and shape of the Tyr side chain is essential for maintaining wild-type values of kcat/Km. Moderate increases in Km values for both the substrate and cofactor upon mutation of Tyr-261 arise mainly from destabilization of the active conformation of a loop containing a dUMP-binding arginine. Besides binding dUMP, this loop has a key role in stabilizing the closed conformation of the enzyme and in shielding the active site from the bulk solvent during catalysis. Changes to atomic vibrations in crystals of a ternary complex of Escherichia coli Tyr261Trp are associated with a greater than 2000-fold drop in kcat/Km. These results underline the important contribution of dynamics to catalysis in TS.

Local flexibility in molecular function paradigm

It is generally accepted that the functional activity of biological macromolecules requires tightly packed three-dimensional structures. Recent theoretical and experimental evidence indicates, however, the importance of molecular flexibility for the proper functioning of some proteins. We examined high resolution structures of proteins in various functional categories with respect to the secondary structure assessment. The latter was considered as a characteristic of the inherent flexibility of a polypeptide chain. We found that the proteins in functionally competent conformational states might be comprised of 20-70% flexible residues. For instance, proteins involved in gene regulation, e.g. transcription factors, are on average largely disordered molecules with over 60% of amino acids residing in "coiled" configurations. In contrast, oxygen transporters constitute a class of relatively rigid molecules with only 30% of residues being locally flexible. Phylogenic comparison of a large number of protein families with respect to the propagation of secondary structure illuminates the growing role of the local flexibility in organisms of greater complexity. Furthermore the local flexibility in protein molecules appears to be dependent on the molecular confinement and is essentially larger in extracellular proteins.

An information theoretic approach to macromolecular modeling: I. Sequence alignments

We are interested in applying the principles of information theory to structural biology calculations. In this article, we explore the information content of an important computational procedure: sequence alignment. Using a reference state developed from exhaustive sequences, we measure alignment statistics and evaluate gap penalties based on first-principle considerations and gap distributions. We show that there are different gap penalties for different alphabet sizes and that the gap penalties can depend on the length of the sequences being aligned. In a companion article, we examine the information content of molecular force fields.

The channel architecture of aquaporin 0 at a 2.2-A resolution

We determined the x-ray structure of bovine aquaporin 0 (AQP0) to a resolution of 2.2 A. The structure of this eukaryotic, integral membrane protein suggests that the selectivity of AQP0 for water transport is based on the identity and location of signature amino acid residues that are hallmarks of the water-selective arm of the AQP family of proteins. Furthermore, the channel lumen is narrowed only by two, quasi-2-fold related tyrosine side chains that might account for reduced water conductance relative to other AQPs. The channel is functionally open to the passage of water because there are eight discreet water molecules within the channel. Comparison of this structure with the recent electron-diffraction structure of the junctional form of sheep AQP0 at pH 6.0 that was interpreted as closed shows no global change in the structure of AQP0 and only small changes in side-chain positions. We observed no structural change to the channel or the molecule as a whole at pH 10, which could be interpreted as the postulated pH-gating mechanism of AQP0-mediated water transport at pH >6.5. Contrary to the electron-diffraction structure, the comparison shows no evidence of channel gating induced by association of the extracellular domains of AQP0 at pH 6.0. Our structure aids the analysis of the interaction of the extracellular domains and the possibility of a cell-cell adhesion role for AQP0. In addition, our structure illustrates the basis for formation of certain types of cataracts that are the result of mutations.

The process of structure-based drug design

The field of structure-based drug design is a rapidly growing area in which many successes have occurred in recent years. The explosion of genomic, proteomic, and structural information has provided hundreds of new targets and opportunities for future drug lead discovery. This review summarizes the process of structure-based drug design and includes, primarily, the choice of a target, the evaluation of a structure of that target, the pivotal questions to consider in choosing a method for drug lead discovery, and evaluation of the drug leads. Key principles in the field of structure-based drug design will be illustrated through a case study that explores drug design for AmpC beta-lactamase.

Local complexity of amino acid interactions in a protein core

Atomic resolution structures of proteins indicate that the core is typically well packed, suggesting a densely connected network of interactions between amino acid residues. The combinatorial complexity of energetic interactions in such a network could be enormous, a problem that limits our ability to relate structure and function. Here, we report a case study of the complexity of amino acid interactions in a localized region within the core of the GFP, a particularly stable and tightly packed molecule. Mutations at three sites within the chromophore-binding pocket display an overlapping pattern of conformational change and are thermodynamically coupled, seemingly consistent with the dense network model. However, crystallographic and energetic analyses of coupling between mutations paint a different picture; pairs of mutations couple through independent "hotspots" in the region of structural overlap. The data indicate that, even in highly stable proteins, the core contains sufficient plasticity in packing to uncouple high-order energetic interactions of residues, a property that is likely general in proteins.

Assessment of the role of computations in structural biology

Computations are now an integrated part of structural biology and are used in data gathering, data processing, and data storage as well as in a full spectrum of theoretical pursuits. In this review, we focus on areas of great promise and call attention to important issues of internal consistency and error analysis.

Improved amino acid flexibility parameters

Protein molecules exhibit varying degrees of flexibility throughout their three-dimensional structures, with some segments showing little mobility while others may be so disordered as to be unresolvable by techniques such as X-ray crystallography. Atomic displacement parameters, or B-factors, from X-ray crystallographic studies give an experimentally determined indication of the degree of mobility in a protein structure. To provide better estimators of amino acid flexibility, we have examined B-factors from a large set of high-resolution crystal structures. Because of the differences among structures, it is necessary to normalize the B-factors. However, many proteins have segments of unusually high mobility, which must be accounted for before normalization can be performed. Accordingly, a median-based method from quality control studies was used to identify outliers. After removal of outliers from, and normalization of, each protein chain, the B-factors were collected for each amino acid in the set. It was found that the distribution of normalized B-factors followed a Gumbel, or extreme value distribution, and the location parameter, or mode, of this distribution was used as an estimator of flexibility for the amino acid. These new parameters have a higher correlation with experimentally determined B-factors than parameters from earlier methods.

The only active mutant of thymidylate synthase D169, a residue far from the site of methyl transfer, demonstrates the exquisite nature of enzyme specificity

Cysteine is the only variant of D169, a cofactor-binding residue in thymidylate synthase, that shows in vivo activity. The 2.4 A crystal structure of Escherichia coli thymidylate synthase D169C in a complex with dUMP and the antifolate CB3717 shows it to be an asymmetric dimer, with only one active site covalently bonded to dUMP. At the active site with covalently bound substrate, C169 S gamma adopts the roles of both carboxyl oxygens of D169, making a 3.6 A S...H[bond]N hydrogen bond to 3-NH of CB3717 and a 3.4 A water-mediated hydrogen bond to H212. Analogous hydrogen bonds formed during the enzyme reaction are important for cofactor binding and are postulated to contribute to catalysis. The C169 side chain is likely to be ionized, making it a better hydrogen bond acceptor than a neutral sulfhydryl group. At the second active site, C169 S gamma makes a shorter (3 A) hydrogen bond to the 3-NH of CB3717, CB3717 is approximately 1.5 A out of its binding site and there is no covalent bond between dUMP and the catalytic cysteine. Changes to partitioning among productive and non-productive conformations of reaction intermediates may contribute as much, if not more, to the diminished activity of this mutant than reduced stabilization of transition states.

Crystal structure of the HIV-1 integrase catalytic core and C-terminal domains: a model for viral DNA binding

Insolubility of full-length HIV-1 integrase (IN) limited previous structure analyses to individual domains. By introducing five point mutations, we engineered a more soluble IN that allowed us to generate multidomain HIV-1 IN crystals. The first multidomain HIV-1 IN structure is reported. It incorporates the catalytic core and C-terminal domains (residues 52-288). The structure resolved to 2.8 A is a Y-shaped dimer. Within the dimer, the catalytic core domains form the only dimer interface, and the C-terminal domains are located 55 A apart. A 26-aa alpha-helix, alpha6, links the C-terminal domain to the catalytic core. A kink in one of the two alpha6 helices occurs near a known proteolytic site, suggesting that it may act as a flexible elbow to reorient the domains during the integration process. Two proteins that bind DNA in a sequence-independent manner are structurally homologous to the HIV-1 IN C-terminal domain, suggesting a similar protein-DNA interaction in which the IN C-terminal domain may serve to bind, bend, and orient viral DNA during integration. A strip of positively charged amino acids contributed by both monomers emerges from each active site of the dimer, suggesting a minimally dimeric platform for binding each viral DNA end. The crystal structure of the isolated catalytic core domain (residues 52-210), independently determined at 1.6-A resolution, is identical to the core domain within the two-domain 52-288 structure.

Active site water molecules revealed in the 2.1 A resolution structure of a site-directed mutant of isocitrate dehydrogenase

Isocitrate dehydrogenase catalyses the two step, acid base, oxidative decarboxylation of isocitrate to alpha-ketoglutarate. Lysine 230 was suggested to act as proton donor based on geometry and spatial proximity to isocitrate. To clarify further the role of lysine 230, we co-crystallized the lysine-to-methionine mutant (K230M) with isocitrate and with alpha-ketoglutarate. Crystals were flash-frozen and the two structures were determined and refined to 2. 1 A. Several new features were identified relative to the wild-type structure. Seven side-chains previously unplaced in the wild-type structure were identified and included in the model, and the amino acid terminus was extended by an alanine residue. Many additional water molecules were identified. Examination of the K230M active sites (K230M isocitrate and K230M-ketoglutarate) revealed that tyrosine 160 protrudes further into the active site in the presence of either isocitrate or alpha-ketoglutarate in K230 M than it does in the wild-type structure. Also, methionine 230 was not as fully extended, and asparagine 232 rotates approximately 30 degrees toward the ligand permitting polar interactions. Outside the active site cleft a tetragonal volume of density was identified as a sulfate molecule. Its location and interactions suggest it may influence the equilibrium between the tetragonal and the orthorhombic forms of isocitrate dehydrogenase. Differences observed in the active site water structure between the wild-type and K230M structures were due to a single point mutation. A water molecule was located in the position equivalent to that occupied by the wild-type epsilon-amine of lysine 230; a water molecule in that location in K230M suggests it may influence catalysis in the mutant. Comparison of K230M complexed with isocitrate and alpha-ketoglutarate illuminates the influence a ligand has on active site water structure.

Native protein fluctuations: the conformational-motion temperature and the inverse correlation of protein flexibility with protein stability

We study the fluctuations of native proteins by exact enumeration using the HP lattice model. The model fluctuations increase with temperature. We observe a low-temperature point, below which large fluctuations are frozen out. This prediction is consistent with the observation by Tilton et al. [R. F. Tilton, Jr., J. C. Dewan, and G. A. Petsko, Biochemistry 31, 2469 (1992)], that the thermal motions of ribonuclease A increase sharply above about 200 K. We also explore protein "flexibility" as defined by Debye-Waller-like factors and solvent accessibilities of core residues to hydrogen exchange. We find that proteins having greater stability tend to have fewer large fluctuations, and hence lower flexibilities. If flexibility is necessary for enzyme catalysis, this could explain why proteins from thermophilic organisms, which are exceptionally stable, may be catalytically inactive at normal temperatures.

Contributions of orientation and hydrogen bonding to catalysis in Asn229 mutants of thymidylate synthase

We have determined structures of binary and ternary complexes of five Asn229 variants of thymidylate synthase (TS) and related their structures to the kinetic constants measured previously. Asn229 forms two hydrogen bonds to the pyrimidine ring of the substrate 2'-deoxyuridine-5'-monophosphate (dUMP). These hydrogen bonds constrain the orientation of dUMP in binary complexes with dUMP, and in ternary complexes with dUMP and the TS cofactor, 5,10-methylene-5,6,7,8-tetrahydrofolate. In N229 mutants, where these hydrogen bonds cannot be made, dUMP binds in a misoriented or more disordered fashion. Most N229 mutants exhibit no activity for the dehalogenation of 5-bromo-dUMP, which requires correct orientation of dUMP against Cys198. Since bound dUMP forms the binding surface against which the pterin ring of cofactor binds, misorientation of dUMP results in higher Km values for cofactor. At the same time, binding of the cofactor aids in ordering and positioning dUMP for catalysis. Hydrophobic mutants, such as N229I, favor an arrangement of solvent molecules and side-chains around the ligands similar to that in a proposed transition state for ternary complex formation in wild-type TS, and kcat values are similar to the wild-type value. Smaller, more hydrophilic mutants favor arrangements of the solvent and side-chains surrounding the ligands that do not resemble the proposed transition state. These changes correspond to decreases in kcat of up to 2000-fold, with only modest increases in Km or Kd. These results are consistent with the proposal that the hydrogen-bonding network between water, dUMP and side-chains in the active-site cavity contributes to catalysis in TS. Asn229 has the unique ability to maintain this critical network, without sterically interfering with dUMP binding.

Distribution analysis of the variation of B-factors of X-ray crystal structures; temperature and structural variations in lysozyme

The B-factor (isotropic temperature factor) data for X-ray structures of hen egg-white lysozyme from the study of Young et al. (Young, Dewan, Nave, and Tilton J. Appl. Cryst. 1993, 26, 309-319) potentially contain information about the relative contributions of static and dynamic variation to these factors. The six structures of the protein were obtained at two widely different temperatures (100 and 298 K), with two crystal forms (monoclinic and tetragonal) and other experimental differences. In addition, the monoclinic lysozyme crystals with two molecules per asymmetric unit allow direct examination of variation between structures determined under identical conditions at both temperatures. The B-factors from these structures all have complex distribution functions as might be expected considering all of the influences that these values must reflect. The empirical cumulative distribution functions (eCDF's) of these data show that they are representative of complex, multicomponent distributions. Distribution analysis using the DANFIP procedure (Wampler, Anal. Biochemistry 1990, 186, 209-218) of the data sets reveals that they can be modeled as four to six Gaussian subpopulations, that these subpopulations do not correlate with specific atom types, specific amino acid residues or fixed locations in the structure. While they do seem to correlate with localized groupings of atoms, these grouping vary from structure to structure even within the same crystal under the same conditions. Temperature seems to have a global effect in this case, but it is clear that other factors including experimental error influence the distribution of B-factors within a given structure. This analysis also helps explain the oft observed lack of atomic level correlation between experimental B-factors and calculated mean square displacements from molecular dynamics simulations.

Binding of the anticancer drug ZD1694 to E. coli thymidylate synthase: assessing specificity and affinity

BACKGROUND

Thymidylate synthase (TS) catalyzes the reductive methylation of deoxyuridine monophosphate (dUMP) by 5, 10-methylenetetrahydrofolate (CH2H4folate) to form deoxythymidine monophosphate (dTMP) and dihydrofolate (H2folate). The essential role of TS in the cell life cycle makes it an attractive target for the development of substrate and cofactor-based inhibitors that may find efficacy as anticancer and antiproliferative drugs. Antifolates that compete specifically with the binding of CH2H4 folate include the cofactor analog CB3717 (10-propargyl-5,8-dideazafolate). However, the development of potent cofactor analog inhibitors of TS, such as CB3717, as drugs has been slowed by their toxicity, which often becomes apparent as hepatic and renal toxicity mediated by the specific chemistry of the compound. Attempts to abolish toxicity in human patients while preserving potency against the target enzyme, have led to the development of ZD1694, which has already shown significant activity against colorectal tumours.

RESULTS

The three dimensional crystallographic structure of ZD1694 in complex with dUMP and Escherichia coli TS has been determined to a resolution of 2.2 . This was used to evaluate the specific structural determinants of ZD1694 potency and to correlate structure/activity relationships between it and the closely related ligand, CB3717. ZD1694 binds to TS in the same manner as CB3717 and H2 folate, but a methyl group on its quinazoline ring, its thiophene ring and the methyl group at N10 are compensated for by plastic accommodation of the enzyme active site coupled with specific rearrangement in the solvent structure. A specific hydrogen bond between the protein and the inhibitor CB3717 is absent in the case of ZD1694 whose monoglutamate tail is reoriented and more well ordered.

CONCLUSIONS

The binding mode of ZD1694 to thymidylate synthase has been determined at atomic resolution. ZD1694 forms a ternary complex with dUMP and participates in the multi-step TS reaction through the covalent bond formation between dUMP and Cys146 thereby competing with CH2H4 folate at the active site. Analysis of this inhibitor ternary complex structure and comparison with that of CB3717 reveals that the enzyme accommodates the differences between the two inhibitors with small shifts in the positions of key active site residues and by repositioning an active site water molecule, thereby preserving a general binding mode of these inhibitors.