Accurate Prediction of 1H NMR Chemical Shifts of Small Molecules Using Machine Learning.

NMR is widely considered the gold standard for organic compound structure determination. As such, NMR is routinely used in organic compound identification, drug metabolite characterization, natural product discovery, and the deconvolution of metabolite mixtures in biofluids (metabolomics and exposomics). In many cases, compound identification by NMR is achieved by matching measured NMR spectra to experimentally collected NMR spectral reference libraries. Unfortunately, the number of available experimental NMR reference spectra, especially for metabolomics, medical diagnostics, or drug-related studies, is quite small. This experimental gap could be filled by predicting NMR chemical shifts for known compounds using computational methods such as machine learning (ML). Here, we describe how a deep learning algorithm that is trained on a high-quality, "solvent-aware" experimental dataset can be used to predict 1H chemical shifts more accurately than any other known method. The new program, called PROSPRE (PROton Shift PREdictor) can accurately (mean absolute error of <0.10 ppm) predict 1H chemical shifts in water (at neutral pH), chloroform, dimethyl sulfoxide, and methanol from a user-submitted chemical structure. PROSPRE (pronounced "prosper") has also been used to predict 1H chemical shifts for >600,000 molecules in many popular metabolomic, drug, and natural product databases.

The plasma metabolome of long COVID patients two years after infection.

One of the major challenges currently faced by global health systems is the prolonged COVID-19 syndrome (also known as "long COVID") which has emerged as a consequence of the SARS-CoV-2 epidemic. It is estimated that at least 30% of patients who have had COVID-19 will develop long COVID. In this study, our goal was to assess the plasma metabolome in a total of 100 samples collected from healthy controls, COVID-19 patients, and long COVID patients recruited in Mexico between 2020 and 2022. A targeted metabolomics approach using a combination of LC-MS/MS and FIA MS/MS was performed to quantify 108 metabolites. IL-17 and leptin were measured in long COVID patients by immunoenzymatic assay. The comparison of paired COVID-19/long COVID-19 samples revealed 53 metabolites that were statistically different. Compared to controls, 27 metabolites remained dysregulated even after two years. Post-COVID-19 patients displayed a heterogeneous metabolic profile. Lactic acid, lactate/pyruvate ratio, ornithine/citrulline ratio, and arginine were identified as the most relevant metabolites for distinguishing patients with more complicated long COVID evolution. Additionally, IL-17 levels were significantly increased in these patients. Mitochondrial dysfunction, redox state imbalance, impaired energy metabolism, and chronic immune dysregulation are likely to be the main hallmarks of long COVID even two years after acute COVID-19 infection.

MagMet: A fully automated web server for targeted nuclear magnetic resonance metabolomics of plasma and serum.

Nuclear magnetic resonance (NMR) spectral analysis of biofluids can be a time-consuming process, requiring the expertise of a trained operator. With NMR becoming increasingly popular in the field of metabolomics, there is a growing need to change this paradigm and to automate the process. Here we introduce MagMet, an online web server, that automates the processing and quantification of 1D 1 H NMR spectra from biofluids-specifically, human serum/plasma metabolites, including those associated with inborn errors of metabolism (IEM). MagMet uses a highly efficient data processing procedure that performs automatic Fourier Transformation, phase correction, baseline optimization, chemical shift referencing, water signal removal, and peak picking/peak alignment. MagMet then uses the peak positions, linewidth information, and J-couplings from its own specially prepared standard metabolite reference spectral NMR library of 85 serum/plasma compounds to identify and quantify compounds from experimentally acquired NMR spectra of serum/plasma. MagMet employs linewidth adjustment for more consistent quantification of metabolites from higher field instruments and incorporates a highly efficient data processing procedure for more rapid and accurate detection and quantification of metabolites. This optimized algorithm allows the MagMet webserver to quickly detect and quantify 58 serum/plasma metabolites in 2.6 min per spectrum (when processing a dataset of 50-100 spectra). MagMet's performance was also assessed using spectra collected from defined mixtures (simulating other biofluids), with >100 previously measured plasma spectra, and from spiked serum/plasma samples simulating known IEMs. In all cases, MagMet performed with precision and accuracy matching the performance of human spectral profiling experts. MagMet is available at http://magmet.ca.

Practical Aspects of NMR-Based Metabolomics.

While NMR-based metabolomics is only about 20 years old, NMR has been a key part of metabolic and metabolism studies for >40 years. Historically, metabolic researchers used NMR because of its high level of reproducibility, superb instrument stability, facile sample preparation protocols, inherently quantitative character, non-destructive nature, and amenability to automation. In this chapter, we provide a short history of NMR-based metabolomics. We then provide a detailed description of some of the practical aspects of performing NMR-based metabolomics studies including sample preparation, pulse sequence selection, and spectral acquisition and processing. The two different approaches to metabolomics data analysis, targeted vs. untargeted, are briefly outlined. We also describe several software packages to help users process NMR spectra obtained via these two different approaches. We then give several examples of useful or interesting applications of NMR-based metabolomics, ranging from applications to drug toxicology, to identifying inborn errors of metabolism to analyzing the contents of biofluids from dairy cattle. Throughout this chapter, we will highlight the strengths and limitations of NMR-based metabolomics. Additionally, we will conclude with descriptions of recent advances in NMR hardware, methodology, and software and speculate about where NMR-based metabolomics is going in the next 5-10 years.

HMDB 5.0: the Human Metabolome Database for 2022.

The Human Metabolome Database or HMDB (https://hmdb.ca) has been providing comprehensive reference information about human metabolites and their associated biological, physiological and chemical properties since 2007. Over the past 15 years, the HMDB has grown and evolved significantly to meet the needs of the metabolomics community and respond to continuing changes in internet and computing technology. This year's update, HMDB 5.0, brings a number of important improvements and upgrades to the database. These should make the HMDB more useful and more appealing to a larger cross-section of users. In particular, these improvements include: (i) a significant increase in the number of metabolite entries (from 114 100 to 217 920 compounds); (ii) enhancements to the quality and depth of metabolite descriptions; (iii) the addition of new structure, spectral and pathway visualization tools; (iv) the inclusion of many new and much more accurately predicted spectral data sets, including predicted NMR spectra, more accurately predicted MS spectra, predicted retention indices and predicted collision cross section data and (v) enhancements to the HMDB's search functions to facilitate better compound identification. Many other minor improvements and updates to the content, the interface, and general performance of the HMDB website have also been made. Overall, we believe these upgrades and updates should greatly enhance the HMDB's ease of use and its potential applications not only in human metabolomics but also in exposomics, lipidomics, nutritional science, biochemistry and clinical chemistry.

NP-MRD: the Natural Products Magnetic Resonance Database.

The Natural Products Magnetic Resonance Database (NP-MRD) is a comprehensive, freely available electronic resource for the deposition, distribution, searching and retrieval of nuclear magnetic resonance (NMR) data on natural products, metabolites and other biologically derived chemicals. NMR spectroscopy has long been viewed as the 'gold standard' for the structure determination of novel natural products and novel metabolites. NMR is also widely used in natural product dereplication and the characterization of biofluid mixtures (metabolomics). All of these NMR applications require large collections of high quality, well-annotated, referential NMR spectra of pure compounds. Unfortunately, referential NMR spectral collections for natural products are quite limited. It is because of the critical need for dedicated, open access natural product NMR resources that the NP-MRD was funded by the National Institute of Health (NIH). Since its launch in 2020, the NP-MRD has grown quickly to become the world's largest repository for NMR data on natural products and other biological substances. It currently contains both structural and NMR data for nearly 41,000 natural product compounds from >7400 different living species. All structural, spectroscopic and descriptive data in the NP-MRD is interactively viewable, searchable and fully downloadable in multiple formats. Extensive hyperlinks to other databases of relevance are also provided. The NP-MRD also supports community deposition of NMR assignments and NMR spectra (1D and 2D) of natural products and related meta-data. The deposition system performs extensive data enrichment, automated data format conversion and spectral/assignment evaluation. Details of these database features, how they are implemented and plans for future upgrades are also provided. The NP-MRD is available at https://np-mrd.org.

MarkerDB: an online database of molecular biomarkers.

MarkerDB is a freely available electronic database that attempts to consolidate information on all known clinical and a selected set of pre-clinical molecular biomarkers into a single resource. The database includes four major types of molecular biomarkers (chemical, protein, DNA [genetic] and karyotypic) and four biomarker categories (diagnostic, predictive, prognostic and exposure). MarkerDB provides information such as: biomarker names and synonyms, associated conditions or pathologies, detailed disease descriptions, detailed biomarker descriptions, biomarker specificity, sensitivity and ROC curves, standard reference values (for protein and chemical markers), variants (for SNP or genetic markers), sequence information (for genetic and protein markers), molecular structures (for protein and chemical markers), tissue or biofluid sources (for protein and chemical markers), chromosomal location and structure (for genetic and karyotype markers), clinical approval status and relevant literature references. Users can browse the data by conditions, condition categories, biomarker types, biomarker categories or search by sequence similarity through the advanced search function. Currently, the database contains 142 protein biomarkers, 1089 chemical biomarkers, 154 karyotype biomarkers and 26 374 genetic markers. These are categorized into 25 560 diagnostic biomarkers, 102 prognostic biomarkers, 265 exposure biomarkers and 6746 predictive biomarkers or biomarker panels. Collectively, these markers can be used to detect, monitor or predict 670 specific human conditions which are grouped into 27 broad condition categories. MarkerDB is available at https://markerdb.ca.

Serum Metabolite Biomarkers for Predicting Residual Feed Intake (RFI) of Young Angus Bulls.

Residual feed intake (RFI) is a feed efficiency measure commonly used in the livestock industry to identify animals that efficiently/inefficiently convert feed into meat or body mass. Selection for low-residual feed intake (LRFI), or feed efficient animals, is gaining popularity among beef producers due to the fact that LRFI cattle eat less and produce less methane per unit weight gain. RFI is a difficult and time-consuming measure to perform, and therefore a simple blood test that could distinguish high-RFI (HRFI) from LRFI animals (early on) would potentially benefit beef farmers in terms of optimizing production or selecting which animals to cull or breed. Using three different metabolomics platforms (nuclear magnetic resonance (NMR) spectrometry, liquid chromatography-tandem mass spectrometry (LC-MS/MS), and inductively coupled plasma mass spectrometry (ICP-MS)) we successfully identified serum biomarkers for RFI that could potentially be translated to an RFI blood test. One set of predictive RFI biomarkers included formate and leucine (best for NMR), and another set included C4 (butyrylcarnitine) and LysoPC(28:0) (best for LC-MS/MS). These serum biomarkers have high sensitivity and specificity (AUROC > 0.85), for distinguishing HRFI from LRFI animals. These results suggest that serum metabolites could be used to inexpensively predict and categorize bovine RFI values. Further validation using a larger, more diverse cohort of cattle is required to confirm these findings.

Candidate serum metabolite biomarkers of residual feed intake and carcass merit in sheep.

Mutton and lamb sales continue to grow globally at a rate of 5% per year. However, sheep farming struggles with low profit margins due to high feed costs and modest carcass yields. Selecting those sheep expected to convert feed efficiently and have high carcass merit, as early as possible in their life cycle, could significantly improve the profitability of sheep farming. Unfortunately, direct measurement of feed conversion efficiency (via residual feed intake [RFI]) and carcass merit is a labor-intensive and expensive procedure. Thus, indirect, marker-assisted evaluation of these traits has been explored as a means of reducing the cost of its direct measurement. One promising and potentially inexpensive route to discover biomarkers of RFI and/or carcass merit is metabolomics. Using quantitative metabolomics, we profiled the blood serum metabolome (i.e., the sum of all measurable metabolites) associated with sheep RFI and carcass merit and identified candidate biomarkers of these traits. The study included 165 crossbred ram-lambs that underwent direct measurement of feed consumption to determine their RFI classification (i.e., low vs. high) using the GrowSafe System over a period 40 d. Carcass merit was evaluated after slaughter using standardized methods. Prior to being sent to slaughter, one blood sample was drawn from each animal, and serum prepared and frozen at -80 °C to limit metabolite degradation. A subset of the serum samples was selected based on divergent RFI and carcass quality for further metabolomic analyses. The analyses were conducted using three analytical methods (nuclear magnetic resonance spectroscopy, liquid chromatography mass spectrometry, and inductively coupled mass spectrometry), which permitted the identification and quantification of 161 unique metabolites. Biomarker analyses identified three significant (P < 0.05) candidate biomarkers of sheep RFI (AUC = 0.80), seven candidate biomarkers of carcass yield grade (AUC = 0.77), and one candidate biomarker of carcass muscle-to-bone ratio (AUC = 0.74). The identified biomarkers appear to have roles in regulating energy metabolism and protein synthesis. These results suggest that serum metabolites could be used to categorize and predict sheep for their RFI and carcass merit. Further validation using a larger (3×) and more diverse cohort of sheep is required to confirm these findings.

Automated Tools for the Analysis of 1D-NMR and 2D-NMR Spectra.

Nuclear magnetic resonance (NMR) spectroscopy is becoming increasingly automated. Most modern NMR spectrometers are now equipped with auto-tune/auto-match probes along with automated locking and shimming systems. Likewise, more and more instruments, especially for NMR-based metabolomics applications, are equipped with automated sample changers. All this instrumental automation allows NMR data to be collected at a rate of >100 samples/day. However, a continuing bottleneck in NMR-based metabolomics has been the time required to manually analyze and annotate the collected NMR spectra. In many cases, manual spectral annotation and analysis can take one or more hours per spectrum. Fortunately, over the past few years, several software tools have been developed that largely automate the spectral deconvolution or spectral annotation process. Using these tools requires that the samples must be prepared and the NMR spectra must be acquired in a very specific manner. In this chapter, we will describe the step-by-step preparation of biofluid samples along with the required protocols for acquiring optimal spectra for automated NMR metabolomics analysis. We will also discuss the use of three common tools (Chenomx NMR Suite, Bayesil, and COLMARm) for (semi-) automated profiling, and annotation of 1D- and 2D-NMR spectra of biofluids.

Rapid and reliable protein structure determination via chemical shift threading.

Protein structure determination using nuclear magnetic resonance (NMR) spectroscopy can be both time-consuming and labor intensive. Here we demonstrate how chemical shift threading can permit rapid, robust, and accurate protein structure determination using only chemical shift data. Threading is a relatively old bioinformatics technique that uses a combination of sequence information and predicted (or experimentally acquired) low-resolution structural data to generate high-resolution 3D protein structures. The key motivations behind using NMR chemical shifts for protein threading lie in the fact that they are easy to measure, they are available prior to 3D structure determination, and they contain vital structural information. The method we have developed uses not only sequence and chemical shift similarity but also chemical shift-derived secondary structure, shift-derived super-secondary structure, and shift-derived accessible surface area to generate a high quality protein structure regardless of the sequence similarity (or lack thereof) to a known structure already in the PDB. The method (called E-Thrifty) was found to be very fast (often < 10 min/structure) and to significantly outperform other shift-based or threading-based structure determination methods (in terms of top template model accuracy)-with an average TM-score performance of 0.68 (vs. 0.50-0.62 for other methods). Coupled with recent developments in chemical shift refinement, these results suggest that protein structure determination, using only NMR chemical shifts, is becoming increasingly practical and reliable. E-Thrifty is available as a web server at http://ethrifty.ca .

HMDB 4.0: the human metabolome database for 2018.

The Human Metabolome Database or HMDB (www.hmdb.ca) is a web-enabled metabolomic database containing comprehensive information about human metabolites along with their biological roles, physiological concentrations, disease associations, chemical reactions, metabolic pathways, and reference spectra. First described in 2007, the HMDB is now considered the standard metabolomic resource for human metabolic studies. Over the past decade the HMDB has continued to grow and evolve in response to emerging needs for metabolomics researchers and continuing changes in web standards. This year's update, HMDB 4.0, represents the most significant upgrade to the database in its history. For instance, the number of fully annotated metabolites has increased by nearly threefold, the number of experimental spectra has grown by almost fourfold and the number of illustrated metabolic pathways has grown by a factor of almost 60. Significant improvements have also been made to the HMDB's chemical taxonomy, chemical ontology, spectral viewing, and spectral/text searching tools. A great deal of brand new data has also been added to HMDB 4.0. This includes large quantities of predicted MS/MS and GC-MS reference spectral data as well as predicted (physiologically feasible) metabolite structures to facilitate novel metabolite identification. Additional information on metabolite-SNP interactions and the influence of drugs on metabolite levels (pharmacometabolomics) has also been added. Many other important improvements in the content, the interface, and the performance of the HMDB website have been made and these should greatly enhance its ease of use and its potential applications in nutrition, biochemistry, clinical chemistry, clinical genetics, medicine, and metabolomics science.

Unraveling the meaning of chemical shifts in protein NMR.

Chemical shifts are among the most informative parameters in protein NMR. They provide wealth of information about protein secondary and tertiary structure, protein flexibility, and protein-ligand binding. In this report, we review the progress in interpreting and utilizing protein chemical shifts that has occurred over the past 25years, with a particular focus on the large body of work arising from our group and other Canadian NMR laboratories. More specifically, this review focuses on describing, assessing, and providing some historical context for various chemical shift-based methods to: (1) determine protein secondary and super-secondary structure; (2) derive protein torsion angles; (3) assess protein flexibility; (4) predict residue accessible surface area; (5) refine 3D protein structures; (6) determine 3D protein structures and (7) characterize intrinsically disordered proteins. This review also briefly covers some of the methods that we previously developed to predict chemical shifts from 3D protein structures and/or protein sequence data. It is hoped that this review will help to increase awareness of the considerable utility of NMR chemical shifts in structural biology and facilitate more widespread adoption of chemical-shift based methods by the NMR spectroscopists, structural biologists, protein biophysicists, and biochemists worldwide. This article is part of a Special Issue entitled: Biophysics in Canada, edited by Lewis Kay, John Baenziger, Albert Berghuis and Peter Tieleman.

A robust algorithm for optimizing protein structures with NMR chemical shifts.

Over the past decade, a number of methods have been developed to determine the approximate structure of proteins using minimal NMR experimental information such as chemical shifts alone, sparse NOEs alone or a combination of comparative modeling data and chemical shifts. However, there have been relatively few methods that allow these approximate models to be substantively refined or improved using the available NMR chemical shift data. Here, we present a novel method, called Chemical Shift driven Genetic Algorithm for biased Molecular Dynamics (CS-GAMDy), for the robust optimization of protein structures using experimental NMR chemical shifts. The method incorporates knowledge-based scoring functions and structural information derived from NMR chemical shifts via a unique combination of multi-objective MD biasing, a genetic algorithm, and the widely used XPLOR molecular modelling language. Using this approach, we demonstrate that CS-GAMDy is able to refine and/or fold models that are as much as 10 Å (RMSD) away from the correct structure using only NMR chemical shift data. CS-GAMDy is also able to refine of a wide range of approximate or mildly erroneous protein structures to more closely match the known/correct structure and the known/correct chemical shifts. We believe CS-GAMDy will allow protein models generated by sparse restraint or chemical-shift-only methods to achieve sufficiently high quality to be considered fully refined and "PDB worthy". The CS-GAMDy algorithm is explained in detail and its performance is compared over a range of refinement scenarios with several commonly used protein structure refinement protocols. The program has been designed to be easily installed and easily used and is available at http://www.gamdy.ca.

A simple method to measure protein side-chain mobility using NMR chemical shifts.

Protein side-chain motions are involved in many important biological processes including enzymatic catalysis, allosteric regulation, and the mediation of protein-protein, protein-DNA, protein-RNA, and protein-cofactor interactions. NMR spectroscopy has long been used to provide insights into the motions of side-chain groups. Currently, the method of choice for studying side-chain dynamics by NMR is the measurement of methyl (2)H autorelaxation. Methyl (2)H autorelaxation exhibits simple relaxation mechanisms and can be straightforwardly converted to meaningful dynamic parameters. However, methyl groups can only be found in 6 of 19 side-chain bearing amino acids. Consequently, only a sparse assessment of protein side-chain dynamics is possible. Therefore, there is a significant interest in developing novel methods of studying side-chain motions that can be applied to all types of side-chains. Here, we show how side-chain chemical shifts can be used to determine the magnitude of fast side-chain motions in proteins. The chemical shift method is applicable to all side-chain bearing residues and does not require any additional measurements beyond standard NMR experiments for backbone and side-chain assignments.

Resolution-by-proxy: a simple measure for assessing and comparing the overall quality of NMR protein structures.

In protein X-ray crystallography, resolution is often used as a good indicator of structural quality. Diffraction resolution of protein crystals correlates well with the number of X-ray observables that are used in structure generation and, therefore, with protein coordinate errors. In protein NMR, there is no parameter identical to X-ray resolution. Instead, resolution is often used as a synonym of NMR model quality. Resolution of NMR structures is often deduced from ensemble precision, torsion angle normality and number of distance restraints per residue. The lack of common techniques to assess the resolution of X-ray and NMR structures complicates the comparison of structures solved by these two methods. This problem is sometimes approached by calculating "equivalent resolution" from structure quality metrics. However, existing protocols do not offer a comprehensive assessment of protein structure as they calculate equivalent resolution from a relatively small number (<5) of protein parameters. Here, we report a development of a protocol that calculates equivalent resolution from 25 measurable protein features. This new method offers better performance (correlation coefficient of 0.92, mean absolute error of 0.28 Å) than existing predictors of equivalent resolution. Because the method uses coordinate data as a proxy for X-ray diffraction data, we call this measure "Resolution-by-Proxy" or ResProx. We demonstrate that ResProx can be used to identify under-restrained, poorly refined or inaccurate NMR structures, and can discover structural defects that the other equivalent resolution methods cannot detect. The ResProx web server is available at http://www.resprox.ca.

Exploring the essential collective dynamics of interacting proteins: application to prion protein dimers.

Essential collective dynamics is a promising and robust approach for analysing the slow motions of macromolecules from short molecular dynamics trajectories. In this study, an extension of the method to treat a collection of interacting protein molecules is presented. The extension is applied to investigate the effects of dimerization on the collective dynamics of ovine prion protein molecules in two different arrangements. Examination of the structural plasticity shows that aggregation has a restricting effect on the local mobility of the prion protein molecules in the interfacial regions. Domain motions of the two dimeric ovine prion protein conformations are distinctly different and can be related to interatomic correlations at the interface. Correlated motions are among the slow collective modes extensively analysed by considering both main-chain and side-chain atoms. Correlation maps reveal the existence of a vast network of dynamically correlated side groups, which extends beyond individual subunits via interfacial interconnections. The network is formed by a core of hydrophobic side chains surrounded by hydrophilic groups at the periphery. The relevance of these findings are discussed in the context of mutations associated with prion diseases. The binding free energy of the dimeric conformations is evaluated to probe their thermodynamic stability. The descriptions afforded by the analysis of the essential collective dynamics of the prion dimers are consistent with their binding free energies. The agreement validates the extension of the methodology and provides a means of interpreting the collective dynamics in terms of the thermodynamic stability of ovine prion proteins.

Use of proteinase K nonspecific digestion for selective and comprehensive identification of interpeptide cross-links: application to prion proteins.

Chemical cross-linking combined with mass spectrometry is a rapidly developing technique for structural proteomics. Cross-linked proteins are usually digested with trypsin to generate cross-linked peptides, which are then analyzed by mass spectrometry. The most informative cross-links, the interpeptide cross-links, are often large in size, because they consist of two peptides that are connected by a cross-linker. In addition, trypsin targets the same residues as amino-reactive cross-linkers, and cleavage will not occur at these cross-linker-modified residues. This produces high molecular weight cross-linked peptides, which complicates their mass spectrometric analysis and identification. In this paper, we examine a nonspecific protease, proteinase K, as an alternative to trypsin for cross-linking studies. Initial tests on a model peptide that was digested by proteinase K resulted in a "family" of related cross-linked peptides, all of which contained the same cross-linking sites, thus providing additional verification of the cross-linking results, as was previously noted for other post-translational modification studies. The procedure was next applied to the native (PrP(C)) and oligomeric form of prion protein (PrPß). Using proteinase K, the affinity-purifiable CID-cleavable and isotopically coded cross-linker cyanurbiotindipropionylsuccinimide and MALDI-MS cross-links were found for all of the possible cross-linking sites. After digestion with proteinase K, we obtained a mass distribution of the cross-linked peptides that is very suitable for MALDI-MS analysis. Using this new method, we were able to detect over 60 interpeptide cross-links in the native PrP(C) and PrPß prion protein. The set of cross-links for the native form was used as distance constraints in developing a model of the native prion protein structure, which includes the 90-124-amino acid N-terminal portion of the protein. Several cross-links were unique to each form of the prion protein, including a Lys(185)-Lys(220) cross-link, which is unique to the PrPß and thus may be indicative of the conformational change involved in the formation of prion protein oligomers.

Solution NMR of a 463-residue phosphohexomutase: domain 4 mobility, substates, and phosphoryl transfer defect.

Phosphomannomutase/phosphoglucomutase contributes to the infectivity of Pseudomonas aeruginosa, retains and reorients its intermediate by 180°, and rotates domain 4 to close the deep catalytic cleft. Nuclear magnetic resonance (NMR) spectra of the backbone of wild-type and S108C-inactivated enzymes were assigned to at least 90%. (13)C secondary chemical shifts report excellent agreement of solution and crystallographic structure over the 14 α-helices, C-capping motifs, and 20 of the 22 ß-strands. Major and minor NMR peaks implicate substates affecting 28% of assigned residues. These can be attributed to the phosphorylation state and possibly to conformational interconversions. The S108C substitution of the phosphoryl donor and acceptor slowed transformation of the glucose 1-phosphate substrate by impairing k(cat). Addition of the glucose 1,6-bisphosphate intermediate accelerated this reaction by 2-3 orders of magnitude, somewhat bypassing the defect and apparently relieving substrate inhibition. The S108C mutation perturbs the NMR spectra and electron density map around the catalytic cleft while preserving the secondary structure in solution. Diminished peak heights and faster (15)N relaxation suggest line broadening and millisecond fluctuations within four loops that can contact phosphosugars. (15)N NMR relaxation and peak heights suggest that domain 4 reorients slightly faster in solution than domains 1-3, and with a different principal axis of diffusion. This adds to the crystallographic evidence of domain 4 rotations in the enzyme, which were previously suggested to couple to reorientation of the intermediate, substrate binding, and product release.

Comparative analysis of essential collective dynamics and NMR-derived flexibility profiles in evolutionarily diverse prion proteins.

Collective motions on ns-µs time scales are known to have a major impact on protein folding, stability, binding and enzymatic efficiency. It is also believed that these motions may have an important role in the early stages of prion protein misfolding and prion disease. In an effort to accurately characterize these motions and their potential influence on the misfolding and prion disease transmissibility we have conducted a combined analysis of molecular dynamic simulations and NMR-derived flexibility measurements over a diverse range of prion proteins. Using a recently developed numerical formalism, we have analyzed the essential collective dynamics (ECD) for prion proteins from 8 different species including human, cow, elk, cat, hamster, chicken, turtle and frog. We also compared the numerical results with flexibility profiles generated by the random coil index (RCI) from NMR chemical shifts. Prion protein backbone flexibility derived from experimental NMR data and from theoretical computations show strong agreement with each other, demonstrating that it is possible to predict the observed RCI profiles employing the numerical ECD formalism. Interestingly, flexibility differences in the loop between second beta strand (S2) and the second alpha helix (HB) appear to distinguish prion proteins from species that are susceptible to prion disease and those that are resistant. Our results show that the different levels of flexibility in the S2-HB loop in various species are predictable via the ECD method, indicating that ECD may be used to identify disease resistant variants of prion proteins, as well as the influence of prion proteins mutations on disease susceptibility or misfolding propensity.