Dynamic phosphometabolomic profiling of human tissues and transgenic models by 18O-assisted ³¹P NMR and mass spectrometry.

Next-generation screening of disease-related metabolomic phenotypes requires monitoring of both metabolite levels and turnover rates. Stable isotope (18)O-assisted (31)P nuclear magnetic resonance (NMR) and mass spectrometry uniquely allows simultaneous measurement of phosphometabolite levels and turnover rates in tissue and blood samples. The (18)O labeling procedure is based on the incorporation of one (18)O into P(i) from [(18)O]H(2)O with each act of ATP hydrolysis and the distribution of (18)O-labeled phosphoryls among phosphate-carrying molecules. This enables simultaneous recording of ATP synthesis and utilization, phosphotransfer fluxes through adenylate kinase, creatine kinase, and glycolytic pathways, as well as mitochondrial substrate shuttle, urea and Krebs cycle activity, glycogen turnover, and intracellular energetic communication. Application of expanded (18)O-labeling procedures has revealed significant differences in the dynamics of G-6-P[(18)O] (glycolysis), G-3-P[(18)O] (substrate shuttle), and G-1-P[(18)O] (glycogenolysis) between human and rat atrial myocardium. In human atria, the turnover of G-3-P[(18)O], which defects are associated with the sudden death syndrome, was significantly higher indicating a greater importance of substrate shuttling to mitochondria. Phosphometabolomic profiling of transgenic hearts deficient in adenylate kinase (AK1-/-), which altered levels and mutations are associated to human diseases, revealed a stress-induced shift in metabolomic profile with increased CrP[(18)O] and decreased G-1-P[(18)O] metabolic dynamics. The metabolomic profile of creatine kinase M-CK/ScCKmit-/--deficient hearts is characterized by a higher G-6-[(18)O]P turnover rate, G-6-P levels, glycolytic capacity, γ/ß-phosphoryl of GTP[(18)O] turnover, as well as ß-[(18)O]ATP and ß-[(18)O]ADP turnover, indicating altered glycolytic, guanine nucleotide, and adenylate kinase metabolic flux. Thus, (18)O-assisted gas chromatography-mass spectrometry and (31)P NMR provide a suitable platform for dynamic phosphometabolomic profiling of the cellular energetic system enabling prediction and diagnosis of metabolic diseases states.

Electron spray ionization mass spectrometry and 2D 31P NMR for monitoring 18O/16O isotope exchange and turnover rates of metabolic oligophosphates.

A new method was here developed for the determination of (18)O-labeling ratios in metabolic oligophosphates, such as ATP, at different phosphoryl moieties (α-, ß-, and γ-ATP) using sensitive and rapid electrospray ionization mass spectrometry (ESI-MS). The ESI-MS-based method for monitoring of (18)O/(16)O exchange was validated with gas chromatography-mass spectrometry and 2D (31)P NMR correlation spectroscopy, the current standard methods in labeling studies. Significant correlation was found between isotopomer selective 2D (31)P NMR spectroscopy and isotopomer less selective ESI-MS method. Results demonstrate that ESI-MS provides a robust analytical platform for simultaneous determination of levels, (18)O-labeling kinetics and turnover rates of α-, ß-, and γ-phosphoryls in ATP molecule. Such method is advantageous for large scale dynamic phosphometabolomic profiling of metabolic networks and acquiring information on the status of probed cellular energetic system.

18O-assisted dynamic metabolomics for individualized diagnostics and treatment of human diseases.

Technological innovations and translation of basic discoveries to clinical practice drive advances in medicine. Today's innovative technologies enable comprehensive screening of the genome, transcriptome, proteome, and metabolome. The detailed knowledge, converged in the integrated "omics" (genomics, transcriptomics, proteomics, and metabolomics), holds an immense potential for understanding mechanism of diseases, facilitating their early diagnostics, selecting personalized therapeutic strategies, and assessing their effectiveness. Metabolomics is the newest "omics" approach aimed to analyze large metabolite pools. The next generation of metabolomic screening requires technologies for high throughput and robust monitoring of metabolite levels and their fluxes. In this regard, stable isotope 18O-based metabolite tagging technology expands quantitative measurements of metabolite levels and turnover rates to all metabolites that include water as a reactant, most notably phosphometabolites. The obtained profiles and turnover rates are sensitive indicators of energy and metabolic imbalances like the ones created by genetic deficiencies, myocardial ischemia, heart failure, neurodegenerative disorders, etc. Here we describe and discuss briefly the potential use of dynamic phosphometabolomic platform for disease diagnostics currently under development at Mayo Clinic.

Calmodulin wraps around its binding domain in the plasma membrane Ca2+ pump anchored by a novel 18-1 motif.

Using solution NMR spectroscopy, we obtained the structure of Ca(2+)-calmodulin (holoCaM) in complex with peptide C28 from the binding domain of the plasma membrane Ca(2+)-ATPase (PMCA) pump isoform 4b. This provides the first atomic resolution insight into the binding mode of holoCaM to the full-length binding domain of PMCA. Structural comparison of the previously determined holoCaM.C20 complex with this holoCaM.C28 complex supports the idea that the initial binding step is represented by (holoCaM.C20) and the final bound complex by (holoCaM.C28). This affirms the existing multi-step kinetic model of PMCA4b activation by CaM. The complex exhibits a new binding motif in which holoCaM is wrapped around helical C28 peptide using two anchoring residues from the peptide at relative positions 18 and 1. The anchors correspond to Phe-1110 and Trp-1093, respectively, in full-length PMCA4b, and the peptide and CaM are oriented in an anti-parallel manner. This is a greater sequence distance between anchors than in any of the known holoCaM complexes with a helical peptide. Analysis of the geometry of holoCaM-peptide binding for the cases where the target peptide adopts an alpha(D)-helix with its anchors buried in the main hydrophobic pockets of the two CaM lobes establishes that only relative sequential positions of 10, 14, 17, and 18 are allowed for the second anchor.

(31)P NMR correlation maps of (18)O/ (16)O chemical shift isotopic effects for phosphometabolite labeling studies.

Intramolecular correlations among the (18)O-labels of metabolic oligophosphates, mapped by J-decoupled (31)P NMR 2D chemical shift correlation spectroscopy, impart stringent constraints to the (18)O-isotope distributions over the whole oligophosphate moiety. The multiple deduced correlations of isotopic labels enable determination of site-specific fractional isotope enrichments and unravel the isotopologue statistics. This approach ensures accurate determination of (18)O-labeling rates of phosphometabolites, critical in biochemical energy conversion and metabolic flux transmission. The biological usefulness of the J-decoupled (31)P NMR 2D chemical shift correlation maps was validated on adenosine tri-phosphate fractionally (18)O labeled in perfused mammalian hearts.

Structural basis of ubiquitin recognition by translesion synthesis DNA polymerase ι.

Cells have evolved mutagenic bypass mechanisms that prevent stalling of the replication machinery at DNA lesions. This process, translesion DNA synthesis (TLS), involves switching from high-fidelity DNA polymerases to specialized DNA polymerases that replicate through a variety of DNA lesions. In eukaryotes, polymerase switching during TLS is regulated by the DNA damage-triggered monoubiquitylation of PCNA. How the switch operates is unknown, but all TLS polymerases of the so-called Y-family possess PCNA and ubiquitin-binding domains that are important for their function. To gain insight into the structural mechanisms underlying the regulation of TLS by ubiquitylation, we have probed the interaction of ubiquitin with a conserved ubiquitin-binding motif (UBM2) of Y-family polymerase Polι. Using NMR spectroscopy, we have determined the structure of a complex of human Polι UBM2 and ubiquitin, revealing a novel ubiquitin recognition fold consisting of two α-helices separated by a central trans-proline residue conserved in all UBMs. We show that, different from the majority of ubiquitin complexes characterized to date, ubiquitin residue Ile44 only plays a modest role in the association of ubiquitin with Polι UBM2. Instead, binding of UBM2 is centered on the recognition of Leu8 in ubiquitin, which is essential for the interaction.

Solvent-induced differentiation of protein backbone hydrogen bonds in calmodulin.

In apo and holoCaM, almost half of the hydrogen bonds (H-bonds) at the protein backbone expected from the corresponding NMR or X-ray structures were not detected by h3JNC' couplings. The paucity of the h3JNC' couplings was considered in terms of dynamic features of these structures. We examined a set of seven proteins and found that protein-backbone H-bonds form two groups according to the h3JNC' couplings measured in solution. H-bonds that have h3JNC' couplings above the threshold of 0.2 Hz show distance/angle correlation among the H-bond geometrical parameters, and appear to be supported by the backbone dynamics in solution. The other H-bonds have no such correlation; they populate the water-exposed and flexible regions of proteins, including many of the CaM helices. The observed differentiation in a dynamical behavior of backbone H-bonds in apo and holoCaM appears to be related to protein functions.

Structural basis of BACH1 phosphopeptide recognition by BRCA1 tandem BRCT domains.

BRCT tandem domains, found in many proteins involved in DNA damage checkpoint and DNA repair pathways, were recently shown to be phosphopeptide binding motifs. Using solution nuclear magnetic resonance (NMR) spectroscopy and mutational analysis, we have characterized the interaction of BRCA1-BRCT domains with a phosphoserine-containing peptide derived from the DNA repair helicase BACH1. We show that a phenylalanine in the +3 position from the phosphoserine of BACH1 is bound to a conserved hydrophobic pocket formed between the two BRCT domains and that recognition of the phosphate group is mediated by lysine and serine side chains from the amino-terminal BRCT domain. Mutations that prevent phosphopeptide binding abolish BRCA1 function in DNA damage-induced checkpoint control. Our NMR data also reveal a dynamic interaction between BRCA1-BRCT and BACH1, where the bound phosphopeptide exists as an equilibrium of two conformations and where BRCA1-BRCT undergoes a transition to a more rigid conformation upon peptide binding.

Calcium-binding proteins afford calibration of dihedral-angle dependence of 3J(NC(gamma)) coupling constant in aspartate and asparagine residues.

Calibration of the 3J(NC(gamma)) couplings across the N-C(alpha)-C(beta)-C(gamma) fragment of aspartate and asparagine residues is afforded by two interactions that produce fixed conformations of the side chains in solution. One is the binding of these side chains to calcium ions; the other is the H-bond interaction of these side chains with a backbone amide.

H-bonding mediates polarization of peptide groups in folded proteins.

The carbon-nitrogen J-couplings in the hydrogen bonding chains of proteins show that H-bonding mediates peptide-group polarization, which results in the general reduction of peptide-group polarity of folded proteins in solution. The net effect is to make large regions of protein secondary structure, especially beta-sheets, intrinsically more hydrophobic, contributing thereby to overall stability of the tertiary structure.

Xanthones from Swertia punctata.

Isolation of 1-O-primeverosyl-3,8-dihydroxy-5-methoxyxanthone and 1-O-gentiobiosyl-3,7-dimethoxy-8-hydroxyxanthone, along with five known xanthones, isobellidifolin, methylbellidifolin, isoswertianin, methylswertianin and norswertianin-1-O-beta-D-glucoside, from the roots of Swertia punctata is reported. In the aerial parts four xanthones, bellidifolin, methylbellidifolin, swertianolin and mangiferin, and flavone-C-glucoside, isoorientin were identified. The chemotaxonomic and pharmacological significance of these results is discussed.

Molecular basis for the association of human E4B U box ubiquitin ligase with E2-conjugating enzymes UbcH5c and Ubc4.

Human E4B, also called UFD2a, is a U box-containing protein that functions as an E3 ubiquitin ligase and an E4 polyubiquitin chain elongation factor. E4B is thought to participate in the proteasomal degradation of misfolded or damaged proteins through association with chaperones. The U box domain is an anchor site for E2 ubiquitin-conjugating enzymes, but little is known of the binding mechanism. Using X-ray crystallography and NMR spectroscopy, we determined the structures of E4B U box free and bound to UbcH5c and Ubc4 E2s. Whereas previously characterized U box domains are homodimeric, we show that E4B U box is a monomer stabilized by a network of hydrogen bonds identified from scalar coupling measurements. These structural studies, complemented by calorimetry- and NMR-based binding assays, suggest an allosteric regulation of UbcH5c and Ubc4 by E4B U box and provide a molecular basis to understand how the ubiquitylation machinery involving E4B assembles.

Directly observed hydrogen bonds at calcium-binding-sites of calmodulin in solution relate to affinity of the calcium-binding.

Molecular tuning to calcium-binding in the EF-hand motif of holo-calmodulin was studied in solution by NMR (h3)J(NC') H-bond couplings. In the N-terminus lobe of holo-calmodulin, the glutamate crucial for Ca(2+) coordination has network of H-bonds weaker than inferred from the X-ray crystal structure. This glutamate at position 12 appears shifted away from the Ca(2+) preferred coordination, which can explain the lower affinity of the calcium-binding to the N-terminus with respect to C-terminus EF hands.

Structural dependencies of protein backbone 2JNC' couplings.

Protein folding can introduce strain in peptide covalent geometry, including deviations from planarity that are difficult to detect, especially for a protein in solution. We have found dependencies in protein backbone (2)J(NC') couplings on the planarity and the relative orientation of the sequential peptide planes. These dependences were observed in experimental (2)J(NC') couplings from seven proteins, and also were supported by DFT calculations for a model tripeptide. Findings indicate that elevated (2)J(NC') couplings may serve as reporters of structural strain in the protein backbone imposed by protein folds. Such information, supplemented with the H-bond strengths derived from (h3)J(NC') couplings, provides useful insight into the overall energy profile of the protein backbone in solution.

Three-dimensional structure of the water-insoluble protein crambin in dodecylphosphocholine micelles and its minimal solvent-exposed surface.

We chose crambin, a hydrophobic and water-insoluble protein originally isolated from the seeds of the plant Crambe abyssinica, as a model for NMR investigations of membrane-associated proteins. We produced isotopically labeled crambin(P22,L25) (variant of crambin containing Pro22 and Leu25) as a cleavable fusion with staphylococcal nuclease and refolded the protein by an approach that has proved successful for the production of proteins with multiple disulfide bonds. We used NMR spectroscopy to determine the three-dimensional structure of the protein in two membrane-mimetic environments: in a mixed aqueous-organic solvent (75%/25%, acetone/water) and in DPC micelles. With the sample in the mixed solvent, it was possible to determine (>NH...OC<) hydrogen bonds directly by the detection of (h3)J(NC)' couplings. H-bonds determined in this manner were utilized in the refinement of the NMR-derived protein structures. With the protein in DPC (dodecylphosphocholine) micelles, we used manganous ion as an aqueous paramagnetic probe to determine the surface of crambin that is shielded by the detergent. With the exception of the aqueous solvent exposed loop containing residues 20 and 21, the protein surface was protected by DPC. This suggests that the protein may be similarly embedded in physiological membranes. The strategy described here for the expression and structure determination of crambin should be applicable to structural and functional studies of membrane active toxins and small membrane proteins.

Electron delocalization mediates the metal-dependent capacity for CH/pi interactions of acetylacetonato chelates.

CH/pi interactions between the coordinated acetylacetonato ligand and phenyl rings were analyzed in the crystal structures from the Cambridge Structural Database and by quantum chemical calculations. The acetylacetonato ligand may engage in two types of interactions: it can be hydrogen atom donor or acceptor. The analysis of crystal structures and calculations show that interactions with the acetylacetonato ligand acting as hydrogen atom donor depend on the metal in an acetylacetonato chelate ring; the chelate rings with soft metals make stronger interactions. The same trend was not observed in the interactions where the acetylacetonato chelate ring acts as the hydrogen atom acceptor.

Prediction of volatile anesthetic binding sites in proteins.

Computational methods designed to predict and visualize ligand protein binding interactions were used to characterize volatile anesthetic (VA) binding sites and unoccupied pockets within the known structures of VAs bound to serum albumin, luciferase, and apoferritin. We found that both the number of protein atoms and methyl hydrogen, which are within approximately 8 A of a potential ligand binding site, are significantly greater in protein pockets where VAs bind. This computational approach was applied to structures of calmodulin (CaM), which have not been determined in complex with a VA. It predicted that VAs bind to [Ca(2+)](4)-CaM, but not to apo-CaM, which we confirmed with isothermal titration calorimetry. The VA binding sites predicted for the structures of [Ca(2+)](4)-CaM are located in hydrophobic pockets that form when the Ca(2+) binding sites in CaM are saturated. The binding of VAs to these hydrophobic pockets is supported by evidence that halothane predominantly makes contact with aliphatic resonances in [Ca(2+)](4)-CaM (nuclear Overhauser effect) and increases the Ca(2+) affinity of CaM (fluorescence spectroscopy). Our computational analysis and experiments indicate that binding of VA to proteins is consistent with the hydrophobic effect and the Meyer-Overton rule.

Cation-pi interaction in a folded polypeptide.

Cationic and aromatic side chains from protein residues interact to stabilize tertiary structure. The stabilization energy originates in part from electrostatic attraction between the cation, and regions of high electron density in pi-orbitals of the aromatic group, leading to the name cation-pi interaction. The lysine and tyrosine containing peptide, N-acetyl-Pro-Pro-Lys-Tyr-Asp-Lys-NH(2), has near uv CD characteristic of tyrosine in a structured environment. Nuclear Overhauser effect (NOE), coupling constant, and ring current chemical shift constraints obtained with (1)H NMR confirm that the peptide (t6p) folds. Simulated annealing consistent with all NMR constraints produces a 40-structure ensemble for t6p with potential energies within one standard deviation of the lowest value observed. Calculated binding energies indicate that cation-pi and cation-phenolic OH interactions exists between the Lys3 and Tyr4 side chains in most of the structures. The t6p peptide in solution is a model for these interactions in a protein. A perturbing electric field from the cationic ground state charge intermingles the excited states of the aromatic group. This intermingling effect may provide a cation-pi signature effect in the tyrosine spectroscopy. The absorption and CD for the lowest energy electronic transitions of the tyrosine phenol were computed for the ensemble. Red-shifted peak energy and hypochromicity in the absorbance band, and decreasing rotational strength, correlates with increasing binding energy of the complex indicating the cation-pi spectroscopic signature. The ensemble average spectroscopic signature effects in t6p are small and in agreement with observation.

Structural dependencies of h3JNC' scalar coupling in protein H-bond chains.

The H-bond ((h3)J(NC')) and peptide bond ((1)J(NC')) scalar couplings establish connectivity of the electronic structure in the H-bond chains of proteins. The correlated changes of (h3)J(NC') and (1)J(NC') couplings extend over several peptide groups in the chains. Consequently, the electronic structure of the H-bond chains can affect (h3)J(NC') in a manner that is independent of the local H-bond geometry. By taking this into account, and by using a more complete set of H-bond geometry parameters, we have predicted (h3)J(NC') couplings in the H-bond chains with deviations commensurate to the standard deviations of the experimentally determined values. We have created a comprehensive database of (h3)J(NC') and (1)J(NC') couplings by measuring the coupling constants in ubiquitin (alphabeta-fold) intestinal fatty acid binding protein (beta-barrel) and carp parvalbumin (alpha-helical).

Saturation transfer difference nuclear magnetic resonance spectroscopy as a method for screening proteins for anesthetic binding.

The effects of anesthetics on cellular function may result from direct interactions between anesthetic molecules and proteins. These interactions have a low affinity and are difficult to characterize. To identify proteins that bind anesthetics, we used nuclear magnetic resonance saturation transfer difference (STD) spectroscopy. The method is based on the nuclear Overhauser effect between bound anesthetic protons and all protein protons. To establish STD as a method for testing anesthetic binding to proteins, we conducted measurements on a series of protein/anesthetic solutions studied before by other methods. STD was able to identify that volatile anesthetics bind to bovine serum albumin, oleic acid reduces halothane binding to bovine serum albumin, and halothane binds to apomyoglobin but not lysozyme. Using STD, we found that halothane binding to calmodulin is Ca2+ -dependent, which demonstrates anesthetic specificity for a protein conformation. Thus, STD is a powerful tool for investigating anesthetic-protein interactions because of its abilities to detect weak binding, to screen a single protein for binding of multiple anesthetics simultaneously, and to detect a change in anesthetic binding caused by conformational changes or competition with other ligands.