Tertiary structure of thiopurine methyltransferase from Pseudomonas syringae, a bacterial orthologue of a polymorphic, drug-metabolizing enzyme

In humans, the enzyme thiopurine methyltransferase (TPMT) metabolizes 6-thiopurine (6-TP) medications, including 6-thioguanine, 6-mercaptopurine and azathioprine, commonly used for immune suppression and for the treatment of hematopoietic malignancies. S-Methylation by TPMT prevents the intracellular conversion of these drugs into active 6-thioguanine nucleotides (6-TGNs). Genetic polymorphisms in the TPMT protein sequence have been associated with decreased tissue enzymatic activities and an increased risk of life-threatening myelo-suppression from standard doses of 6-TP medications. Biochemical studies have demonstrated that TPMT deficiency is primarily associated with increased degradation of the polymorphic proteins through an ubiquitylation and proteasomal-dependent pathway. We have now determined the tertiary structure of the bacterial orthologue of TPMT from Pseudomonas syringae using NMR spectroscopy. Bacterial TPMT similarly catalyzes the S-adenosylmethionine (SAM)-dependent transmethylation of 6-TPs and shares 45% similarity (33% identity) with the human enzyme. Initial studies revealed an unstructured N terminus, which was removed for structural studies and subsequently determined to be required for enzymatic activity. Despite lacking sequence similarity to any protein of known three-dimensional structure, the tertiary structure of bacterial TPMT reveals a classical SAM-dependent methyltransferase topology, consisting of a seven-stranded beta-sheet flanked by alpha-helices on both sides. However, some deviations from the consensus topology, along with multiple insertions of structural elements, are evident. A review of the many experimentally determined tertiary structures of SAM-dependent methyltransferases demonstrates that such structural deviations from the consensus topology are common and often functionally important.     

Characterization of the Arabidopsis thaliana Arath;CDC25 dual-specificity tyrosine phosphatase

CDC25 enzymes are dual-specificity phosphatases involved in the regulation of the cell cycle. No CDC25 enzymes have been described in higher plant organisms. We report here the characterization of an Arabidopsis thaliana CDC25 enzyme, constituted by a sole catalytic domain and devoid of the N-terminal regulatory region found in the human CDC25. We describe the recombinant expression in Escherichia coli of the Arath;CDC25 and its purification for activity assay and structure determination by NMR. The recombinant enzyme has a tyrosine phosphatase activity towards an artificial substrate, a NMR characterization equally concludes to its correct folding. The secondary structure of the protein was predicted on the basis of the assigned chemical shift of (1)H, (15)N, and (13)C backbone atoms of the protein. The presence of a metal ion in the C-terminus of this new protein points to a zinc finger, and sequence homology indicates that this new structural element might be conserved in related plant homologs.     

Four human thiopurine s-methyltransferase alleles severely affect protein structure and dynamics

Thiopurine S-methyltransferase (TPMT) metabolizes cytotoxic thiopurine drugs used in the treatment of leukemia and inflammatory bowel disease. TPMT is a major pharmacogenomic target with 23 alleles identified to date. Several of these alleles cause rapid protein degradation and/or aggregation, making it extremely difficult to study the structural impact of the TPMT polymorphisms experimentally. We, therefore, have performed multiple molecular dynamics simulations of the four most common alleles [TPMT*2 (A80P), *3A (A154T/Y240C), *3B (A154T) and *3C (Y240C)] to investigate the molecular mechanism of TPMT inactivation at an atomic level. The A80P polymorphism in TPMT*2 disrupts helix α3 bordering the active site, which breaks several salt-bridge interactions and opens up a large cleft in the protein. The A154T polymorphism is located within the co-substrate binding site. The larger threonine alters the packing of substrate-binding residues (P68, L69, Y166), increasing the solvent exposure of the polymorphic site. This packing rearrangement may account for the complete lack of activity in the A154T mutant. The Y240C polymorphism is located in β-strand 9, distant from the active site. Side-chain contacts between residue 240 and helix α8 are lost in TPMT*3C. Residues 154 and 240 in TPMT*3A are connected through a hydrogen-bonding network. The dual polymorphisms result in a flattened, slightly distorted protein structure and an increase in the thiopurine-binding site solvent accessibility. The two variants that undergo the most rapid degradation in vivo, TPMT*2 and *3A, are also the most deformed in the simulations.     

Structures of the FXYD regulatory proteins in lipid micelles and membranes

The FXYD membrane proteins constitute a family of conserved auxiliary subunits of the Na,K-ATPase, and have been the focus of recent attention due to their ability to finely regulate the activity of the enzyme complex in various physiological settings. In this review we describe the structures of the proteins, as well as their dynamics and their associations with the lipid bilayer membrane, which we have recently determined by NMR spectroscopy. Although the proteins are relatively small, their genes contain as many as six to nine small exons, and the coincidence of structured protein segments with their genetic elements suggests assembly from discrete structural modules through exon shuffling. The three-dimensional structures and backbone dynamics provide the foundation for understanding their intra-membrane association with the Na,K-ATPase α subunit, and the structure of FXYD1 suggests a mechanism whereby the phosphorylation of conserved Ser residues, by protein kinases A and C, could induce a conformational change in the cytoplasmic domain of the protein, to modulate its interaction with the α subunit.     

Solution structures, dynamics, and lipid-binding of the sterile alpha-motif domain of the deleted in liver cancer 2

The deleted in liver cancer 2 (DLC2) is a tumor suppressor gene, frequently found to be underexpressed in hepatocellular carcinoma. DLC2 is a multidomain protein containing a sterile alpha-motif (SAM) domain, a GTPase-activating protein (GAP) domain, and a lipid-binding StAR-related lipid-transfer (START) domain. The SAM domain of DLC2, DLC2-SAM, exhibits a low level of sequence homology (15-30%) with other SAM domains, and appears to be the prototype of a new subfamily of SAM domains found in DLC2-related proteins. In the present study, we have determined the three-dimensional solution structure of DLC2-SAM using NMR methods together with molecular dynamics simulated annealing. In addition, we performed a backbone dynamics study. The DLC2-SAM packed as a unique four alpha-helical bundle stabilized by interhelix hydrophobic interactions. The arrangement of the four helices is distinct from all other known SAM domains. In contrast to some members of the SAM domain family which form either dimers or oligomers, both biochemical analyses and rotational correlation time (tau(c)) measured by backbone 15N relaxation experiments indicated that DLC2-SAM exists as a monomer in solution. The interaction of DLC2-SAM domain with sodium dodecyl sulfate (SDS) micelles and 1,2-dimyristoyl-sn-glycerol-3-phosphatidylglycerol (DMPG) phospholipids was examined by CD and NMR spectroscopic techniques. The DLC2-SAM exhibits membrane binding properties accompanied by minor loss of the secondary structure of the protein. Deletion studies showed that the self-association of DLC2 in vivo does not require SAM domain, instead, a protein domain consisting of residues 120-672 mediates the self-association of DLC2.     

NMR structure of the p63 SAM domain and dynamical properties of G534V and T537P pathological mutants, identified in the AEC syndrome

The p63 protein is crucial for epidermal development, and its mutations cause the extrodactyly ectodermal dysplasia and cleft lip/palate syndrome. The three-dimensional solution structure of the p63 sterile α-motif (SAM) domain (residues 505–579), a region crucial to explaining the human genetic disease ankyloblepharon-ectodermal dysplasia-clefting syndrome (AEC), has been determined by nuclear magnetic resonance spectroscopy. The structure indicates that the domain is a monomer with the characteristic five-helix bundle topology observed in other SAM domains. It includes five tightly packed helices with an extended hydrophobic core to form a globular and compact structure. The dynamics of the backbone and the global correlation time of the molecule have also been investigated and compared with the dynamical properties obtained through molecular dynamics simulation. Attempts to purify the pathological G534V and T537P mutants, originally identified in AEC, were not successful because of the occurrence of unspecific proteolytic degradation of the mutated SAM domains. Analysis of the structural dynamic properties of the G534V and T537P mutants through molecular dynamics simulation and comparison with the wild type permits detection of differences in the degree of free-dom of individual residues and discussion of the possible causes for the pathology.     

Discrimination between Closely Related Cellular Metabolites by the SAM-I Riboswitch

The SAM-I riboswitch is a cis-acting element of genetic control found in bacterial mRNAs that specifically binds S-adenosylmethionine (SAM). We previously determined the 2.9-Å X-ray crystal structure of the effector-binding domain of this RNA element, revealing details of RNA-ligand recognition. To improve this structure, variations were made to the RNA sequence to alter lattice contacts, resulting in a 0.5-Å improvement in crystallographic resolution and allowing for a more accurate refinement of the crystallographic model. The basis for SAM specificity was addressed by a structural analysis of the RNA complexed to S-adenosylhomocysteine (SAH) and sinefungin and by measuring the affinity of SAM and SAH for a series of mutants using isothermal titration calorimetry. These data illustrate the importance of two universally conserved base pairs in the RNA that form electrostatic interactions with the positively charged sulfonium group of SAM, thereby providing a basis for discrimination between SAM and SAH.     

Local structural preferences of calpastatin, the intrinsically unstructured protein inhibitor of calpain

Calpain, the calcium-activated intracellular cysteine protease, is under the tight control of its intrinsically unstructured inhibitor, calpastatin. Understanding how potent inhibition by calpastatin can be reconciled with its unstructured nature provides deeper insight into calpain function and a more general understanding of how proteins devoid of a well-defined structure carry out their function. To this end, we performed a full NMR assignment of hCSD1 to characterize it in its solution state. Secondary chemical shift values and NMR relaxation data, R 1, R 2, and hetero-NOE, as well as spectral density function analysis have shown that conserved regions of calpastatin, subdomains A and C, which are responsible for calcium-dependent anchoring of the inhibitor to the enzyme, preferentially sample partially helical backbone conformations of a reduced flexibility. Moreover, the linker regions between subdomains are more flexible with no structural preference. The primary determinant of calpain inhibition, subdomain B, also has a non-fully random conformational preference, resembling a beta-turn structure also ascertained by prior studies of a 27-residue peptide encompassing the inhibitory region. This local structural preference is also confirmed by a deviation in chemical shift values between full-length calpastatin domain 1 and a truncated construct cut in the middle of subdomain B. At the C-terminal end of the molecule, a nascent helical region was found, which in contrast to the overall structural properties of the molecule may indicate a previously unknown functional region. Overall, these observations provide further evidence that supports previous suggestions that intrinsically unstructured proteins use preformed structural elements in efficient partner recognition.     

Structural mass spectrometry decodes domain interaction and dynamics of the full-length Human Histone Deacetylase 2

Human Histone Deacetylase 2 (HDAC2) belongs to a conserved enzyme superfamily that regulates deacetylation inside cells. HDAC2 is a drug target as it is known to be upregulated in cancers and neurodegenerative disorders. It consists of globular deacetylase and C-terminus intrinsically-disordered domains [1–3]. To date, there is no full-length structure of HDAC2 available due to the high intrinsic flexibility of its C-terminal domain. The intrinsically-disordered domain, however, is known to be important for the enzymatic function of HDAC2 [1, 4].Here we combine several structural Mass Spectrometry (MS) methodologies such as denaturing, native, ion mobility and chemical crosslinking, alongside biochemical assays and molecular modelling to study the structure and dynamics of the full-length HDAC2 for the first time. We show that MS can easily dissect heterogeneity inherent within the protein sample and at the same time probe the structural arrangement of the different conformers present.Activity assays combined with data from MS and molecular modelling suggest how the structural dynamics of the C-terminal domain, and its interactions with the catalytic domain, regulate the activity of this enzyme.     

Kinetic and spectroscopic studies of Tritrichomonas foetus low-molecular weight phosphotyrosyl phosphatase. Hydrogen bond networks and electrostatic effects

Virtually all of the eukaryotic low-molecular weight protein tyrosine phosphatases (LMW PTPases) studied to date contain a conserved, high-pK(a) histidine residue that is hydrogen bonded to a conserved active site asparagine residue of the phosphate binding loop. However, in the putative enzyme encoded by the genome of the trichomonad parasite Tritrichomonas foetus, this otherwise highly conserved histidine is replaced with a glutamine residue. We have cloned the gene, expressed the enzyme, demonstrated its catalytic activity, and examined the structural and functional roles of the glutamine residue using site-directed mutagenesis, kinetic measurements, and NMR spectroscopy. Titration studies of the two native histidine residues in the T. foetus enzyme as monitored by (1)H NMR revealed that H44 has a pK(a) of 6.4 and H143 has a pK(a) of 5.3. When a histidine residue was introduced in place of the native glutamine at position 67, a pK(a) of 8.2 was measured for this residue. Steady state kinetic methods were employed to study how mutation of the native glutamine to alanine, asparagine, and histidine affected the catalytic activity of the enzyme. Examination of k(cat)/K(m) showed that Q67H exhibits a substrate selectivity comparable to that of the wild-type (WT) enzyme, while Q67N and Q67A show reduced activity. The effect of pH on the reaction rate was examined. Importantly, the pH-rate profile of the WT TPTP enzyme revealed a much more clearly defined acidic limb than that which can be observed for other wild-type LMW PTPases. The pH-rate curve of the Q67H mutant shows a shift to a lower pH optimum relative to that seen for the wild-type enzyme. The Q67N and Q67A mutants showed curves that were shifted to higher pH optima. Although the active site of this enzyme is likely to be similar to that of other LMW PTPases, the hydrogen bonding and electrostatic changes afford new insight into factors affecting the pH dependence and catalysis by this family of enzymes.     

Characterisation of the conformational properties of urea-unfolded Im7: implications for the early stages of protein folding

The colicin immunity protein Im7 folds from its unfolded state in 6 M urea to its native four-helix structure through an on-pathway intermediate that lacks one of the helices of the native structure (helix III). In order to further characterize the folding mechanism of Im7, we have studied the conformational properties of the protein unfolded in 6 M urea in detail using heteronuclear NMR. Triple-resonance experiments with 13C/15N-labelled Im7 in 6 M urea provided almost complete resonance assignments for the backbone nuclei, and measurement of backbone 15N relaxation parameters allowed dynamic ordering of the unfolded polypeptide chain to be investigated. Reduced spectral density mapping and fitting backbone R2 relaxation rates to a polymer dynamics model identified four clusters of interacting residues, each predicted by the average area buried upon folding for each residue. Chemical shift analyses and measurement of NOEs detected with a long mixing-time 1H-1H-15N NOESY-HSQC spectrum confirmed the formation of four clusters. Each cluster of interacting side-chains in urea-unfolded Im7 occurs in a region of the protein that forms a helix in the protein, with the largest clusters being associated with the three long helices that are formed in the on-pathway folding intermediate, whilst the smallest cluster forms a helix only in the native state. NMR studies of a Phe15Ala Im7 variant and a protein in which residues 51-56 are replaced by three glycine residues (H3G3 Im7*), indicated that the clusters do not interact with each other, possibly because they are solvated by urea, as indicated by analysis of NOEs between the protein and the solvent. Based on these data, we suggest that dilution of the chaotrope to initiate refolding will result in collapse of the clusters, leading to the formation of persistent helical structure and the generation of the three-helix folding intermediate.     

Enzyme dynamics along the reaction coordinate: critical role of a conserved residue

Conformational flexibility of the enzyme architecture is essential for biological function. These structural transitions often encompass significant portions of the enzyme molecule. Here, we present a detailed study of functionally relevant RNase A dynamics in the wild type and a D121A mutant form by NMR spin-relaxation techniques. In the wild-type enzyme, the dynamic properties are largely conserved in the apo, enzyme-substrate, and enzyme-product complexes. In comparison, mutation of aspartic acid 121 to alanine disrupts the timing of active-site dynamics, the product-release step, and global conformational changes, indicating that D121 plays a significant role in coordinating the dynamic events in RNase A. In addition, this mutation results in 90% loss of catalytic activity despite the absence of direct participation of D121 in the chemical reaction or in interactions with the substrate. These data suggest that one role of this conserved residue is to facilitate important millisecond protein dynamics.     

PPM: a side-chain and backbone chemical shift predictor for the assessment of protein conformational ensembles

The combination of the wide availability of protein backbone and side-chain NMR chemical shifts with advances in understanding of their relationship to protein structure makes these parameters useful for the assessment of structural-dynamic protein models. A new chemical shift predictor (PPM) is introduced, which is solely based on physical?Cchemical contributions to the chemical shifts for both the protein backbone and methyl-bearing amino-acid side chains. To explicitly account for the effects of protein dynamics on chemical shifts, PPM was directly refined against 100?ns long molecular dynamics (MD) simulations of 35 proteins with known experimental NMR chemical shifts. It is found that the prediction of methyl-proton chemical shifts by PPM from MD ensembles is improved over other methods, while backbone C??, C??, C??, N, and H^N chemical shifts are predicted at an accuracy comparable to the latest generation of chemical shift prediction programs. PPM is particularly suitable for the rapid evaluation of large protein conformational ensembles on their consistency with experimental NMR data and the possible improvement of protein force fields from chemical shifts.     

Protein dynamics detected in a membrane-embedded potassium channel using two-dimensional solid-state NMR spectroscopy

We report longitudinal ¹⁵N relaxation rates derived from two-dimensional (¹⁵N, ¹³C) chemical shift correlation experiments obtained under magic angle spinning for the potassium channel KcsA-Kv1.3 reconstituted in multilamellar vesicles. Thus, we demonstrate that solid-state NMR can be used to probe residue-specific backbone dynamics in a membrane-embedded protein. Enhanced backbone mobility was detected for two glycine residues within the selectivity filter that are highly conserved in potassium channels and that are of core relevance to the filter structure and ion selectivity.     

Solution structure of the Vts1 SAM domain in the presence of RNA

The yeast Vts1 SAM (sterile alpha motif) domain is a member of a new class of SAM domains that specifically bind RNA. To elucidate the structural basis for RNA binding, the solution structure of the Vts1 SAM domain, in the presence of a specific target RNA, has been solved by multidimensional heteronuclear NMR spectroscopy. The Vts1 SAM domain retains the "core" five-helix-bundle architecture of traditional SAM domains, but has additional short helices at N and C termini, comprising a small substructure that caps the core helices. The RNA-binding surface of Vts1, determined by chemical shift perturbation, maps near the ends of three of the core helices, in agreement with mutational data and the electrostatic properties of the molecule. These results provide a structural basis for the versatility of the SAM domain in protein and RNA-recognition.     

Solution structure and backbone dynamics of human liver fatty acid binding protein: fatty acid binding revisited

Liver fatty acid binding protein (L-FABP), a cytosolic protein most abundant in liver, is associated with intracellular transport of fatty acids, nuclear signaling, and regulation of intracellular lipolysis. Among the members of the intracellular lipid binding protein family, L-FABP is of particular interest as it can i), bind two fatty acid molecules simultaneously and ii), accommodate a variety of bulkier physiological ligands such as bilirubin and fatty acyl CoA. To better understand the promiscuous binding and transport properties of L-FABP, we investigated structure and dynamics of human L-FABP with and without bound ligands by means of heteronuclear NMR. The overall conformation of human L-FABP shows the typical β-clam motif. Binding of two oleic acid (OA) molecules does not alter the protein conformation substantially, but perturbs the chemical shift of certain backbone and side-chain protons that are involved in OA binding according to the structure of the human L-FABP/OA complex. Comparison of the human apo and holo L-FABP structures revealed no evidence for an "open-cap" conformation or a "swivel-back" mechanism of the K90 side chain upon ligand binding, as proposed for rat L-FABP. Instead, we postulate that the lipid binding process in L-FABP is associated with backbone dynamics.     

NMR structure of a heterodimeric SAM:SAM complex: characterization and manipulation of EphA2 binding reveal new cellular functions of SHIP2

The sterile alpha motif (SAM) for protein-protein interactions is encountered in over 200 proteins, but the structural basis for its interactions is just becoming clear. Here we solved the structure of the?EphA2-SHIP2 SAM:SAM heterodimeric complex by use of NMR restraints from chemical shift perturbations, NOE and RDC experiments. Specific contacts between the protein surfaces differ significantly from a previous model and other SAM:SAM complexes. Molecular dynamics and docking simulations indicate fluctuations in the complex toward alternate, higher energy conformations. The interface suggests that EphA family members bind to SHIP2 SAM, whereas EphB members may not; correspondingly, we demonstrate binding of EphA1, but not of EphB2, to SHIP2. A variant of EphB2 SAM was designed that binds SHIP2. Functional characterization of a mutant EphA2 compromised in SHIP2 binding reveals two previously unrecognized functions of SHIP2 in suppressing ligand-induced activation of EphA2 and in promoting receptor coordinated chemotactic cell migration.     

SAM1 domain of SASH1 harbors distinctive structural heterogeneity

The sterile alpha motif (SAM) domains are among the most versatile protein domains in biology, and the variety of the oligomerization states contribute to their diverse roles in many diseases. A better understanding of the structure and dynamics of various SAM domains will provide a scientific basis for drug development targeting them. Here, we used SEC-MALS, HPLC, NMR, and other biophysical techniques to characterize the structural features and dynamics of the SAM1 domain in SASH1. SASH1 is a scaffold protein belonging to the same family as SASH3. Unlike the dimerization seen in SASH3′s SAM domain, our SEC-MALS and SE-HPLC showed that SAM1 exists primarily as a less compact monomer with a minor oligomer. NMR assignment, relaxation, and exchange experiments revealed the presence of both a disordered monomer and a more structured oligomer with multiple timescale exchange regimes in solution. Mutagenesis and SE-HPLC showed that D663A/T664K substitutions in SAM1 increased its oligomerization. In sum, this study is the first to characterize a disordered structure for a SAM domain, provides additional evidence and framework for the diversity of SAM domains, and identifies a region in SAM1 as a potential starting point to further characterize the structural mechanism of oligomerization of the domain.     

Identification of a residue critical for maintaining the functional conformation of BPTI

The effects of amino acid replacements on the backbone dynamics of bovine pancreatic trypsin inhibitor (BPTI) were examined using 15N NMR relaxation experiments. Previous studies have shown that backbone amide groups within the trypsin-binding region of the wild-type protein undergo conformational exchange processes on the micros time scale, and that replacement of Tyr35 with Gly greatly increases the number of backbone atoms involved in such motions. In order to determine whether these mutational effects are specific to the replacement of this residue with Gly, six additional replacements were examined in the present study. In two of these, Tyr35 was replaced with either Ala or Leu, and the other four were single replacements of Tyr23, Phe33, Asn43 or Asn44, all of which are highly buried in the native structure and conserved in homologous proteins. The Y35A and Y35L mutants displayed dynamic properties very similar to those of the Y35G mutant, with the backbone segments including residues 10-19 and 32-44 undergoing motions revealed by enhanced 15N transverse relaxation rates. On the other hand, the Y23L, N43G and N44A substitutions caused almost no detectable changes in backbone dynamics, on either the ns-ps or ms-micros time scales, even though each of these replacements significantly destabilizes the native conformation. Replacement of Phe33 with Leu caused intermediate effects, with several residues that have previously been implicated in motions in the wild-type protein displaying enhanced transverse relaxation rates. These results demonstrate that destabilizing amino acid replacements can be accommodated in a native protein with dramatically different effects on conformational dynamics and that Tyr35 plays a particularly important role in defining the conformation of the trypsin-binding site of BPTI.