The Superfund Research Program Analytics Portal: linking environmental chemical exposure to biological phenotypes.

The OSU/PNNL Superfund Research Program (SRP) represents a longstanding collaboration to quantify Polycyclic Aromatic Hydrocarbons (PAHs) at various superfund sites in the Pacific Northwest and assess their potential impact on human health. To link the chemical measurements to biological activity, we describe the use of the zebrafish as a high-throughput developmental toxicity model that provides quantitative measurements of the exposure to chemicals. Toward this end, we have linked over 150 PAHs found at Superfund sites to the effect of these same chemicals in zebrafish, creating a rich dataset that links environmental exposure to biological response. To quantify this response, we have implemented a dose-response modelling pipeline to calculate benchmark dose parameters which enable potency comparison across over 500 chemicals and 12 of the phenotypes measured in zebrafish. We provide a rich dataset for download and analysis as well as a web portal that provides public access to this dataset via an interactive web site designed to support exploration and re-use of these data by the scientific community at http://srp.pnnl.gov .

De novo sequencing and native mass spectrometry revealed hetero-association of dirigent protein homologs and potential interacting proteins in Forsythia × intermedia.

The discovery of dirigent proteins (DPs) and their functions in plant phenol biochemistry was made over two decades ago with Forsythia × intermedia. Stereo-selective, DP-guided, monolignol-derived radical coupling in vitro was then reported to afford the optically active lignan, (+)-pinoresinol from coniferyl alcohol, provided one-electron oxidase/oxidant capacity was present. It later became evident that DPs have several distinct sub-families, presumably with different functions. Some known DPs require other essential enzymes/proteins (e.g. oxidases) for their functions. However, the lack of a fully sequenced genome for Forsythia × intermedia made it difficult to profile other components co-purified with the (+)-pinoresinol forming DP. Herein, we used an integrated bottom-up, top-down, and native mass spectrometry (MS) approach to de novo sequence the extracted proteins via adaptation of our initial report of DP solubilization and purification. Using publicly available transcriptome and genomic data from closely related species, we identified 14 proteins that were putatively associated with either DP function or the cell wall. Although their co-occurrence after extraction and chromatographic separation is suggestive for potential protein-protein interactions, none were found to form stable protein complexes with DPs in native MS under the specific experimental conditions we have explored. Interestingly, two new DP homologs were found and they formed hetero-trimers. Molecular dynamics simulations suggested that similar hetero-trimers were possible between Arabidopsis DP homologs with comparable sequence similarities. Nevertheless, our integrated mass spectrometry method development helped prepare for future investigations directed to the discovery of novel proteins and protein-protein interactions. These advantages can be highly beneficial for plant and microbial research where fully sequenced genomes may not be readily available.

Quantum Mechanical Methods Predict Accurate Thermodynamics of Biochemical Reactions.

Thermodynamics plays a crucial role in regulating the metabolic processes in all living organisms. Accurate determination of biochemical and biophysical properties is important to understand, analyze, and synthetically design such metabolic processes for engineered systems. In this work, we extensively performed first-principles quantum mechanical calculations to assess its accuracy in estimating free energy of biochemical reactions and developed automated quantum-chemistry (QC) pipeline (https://appdev.kbase.us/narrative/45710) for the prediction of thermodynamics parameters of biochemical reactions. We benchmark the QC methods based on density functional theory (DFT) against different basis sets, solvation models, pH, and exchange-correlation functionals using the known thermodynamic properties from the NIST database. Our results show that QC calculations when combined with simple calibration yield a mean absolute error in the range of 1.60-2.27 kcal/mol for different exchange-correlation functionals, which is comparable to the error in the experimental measurements. This accuracy over a diverse set of metabolic reactions is unprecedented and near the benchmark chemical accuracy of 1 kcal/mol that is usually desired from DFT calculations.

Exploring the Impacts of Conformer Selection Methods on Ion Mobility Collision Cross Section Predictions.

The prediction of structure dependent molecular properties, such as collision cross sections as measured using ion mobility spectrometry, are crucially dependent on the selection of the correct population of molecular conformers. Here, we report an in-depth evaluation of multiple conformation selection techniques, including simple averaging, Boltzmann weighting, lowest energy selection, low energy threshold reductions, and similarity reduction. Generating 50 000 conformers each for 18 molecules, we used the In Silico Chemical Library Engine (ISiCLE) to calculate the collision cross sections for the entire data set. First, we employed Monte Carlo simulations to understand the variability between conformer structures as generated using simulated annealing. Then we employed Monte Carlo simulations to the aforementioned conformer selection techniques applied on the simulated molecular property: the ion mobility collision cross section. Based on our analyses, we found Boltzmann weighting to be a good trade-off between precision and theoretical accuracy. Combining multiple techniques revealed that energy thresholds and root-mean-squared deviation-based similarity reductions can save considerable computational expense while maintaining property prediction accuracy. Molecular dynamic conformer generation tools like AMBER can continue to generate new lowest energy conformers even after tens of thousands of generations, decreasing precision between runs. This reduced precision can be ameliorated and theoretical accuracy increased by running density functional theory geometry optimization on carefully selected conformers.

The multi-dimensional embryonic zebrafish platform predicts flame retardant bioactivity.

Flame retardant chemicals (FRCs) commonly added to many consumer products present a human exposure burden associated with adverse health effects. Under pressure from consumers, FRC manufacturers have adopted some purportedly safer replacements for first-generation brominated diphenyl ethers (BDEs). In contrast, second and third-generation organophosphates and other alternative chemistries have limited bioactivity data available to estimate their hazard potential. In order to evaluate the toxicity of existing and potential replacement FRCs, we need efficient screening methods. We built a 61-FRC library in which we systemically assessed developmental toxicity and potential neurotoxicity effects in the embryonic zebrafish model. Data were compared to publicly available data generated in a battery of cell-based in vitro assays from ToxCast, Tox21, and other alternative models. Of the 61 FRCs, 19 of 45 that were tested in the ToxCast assays were bioactive in our zebrafish model. The zebrafish assays detected bioactivity for 10 of the 12 previously classified developmental neurotoxic FRCs. Developmental zebrafish were sufficiently sensitive at detecting FRC structure-bioactivity impacts that we were able to build a classification model using 13 physicochemical properties and 3 embryonic zebrafish assays that achieved a balanced accuracy of 91.7%. This work illustrates the power of a multi-dimensional in vivo platform to expand our ability to predict the hazard potential of new compounds based on structural relatedness, ultimately leading to reliable toxicity predictions based on chemical structure.

Pterocarpan synthase (PTS) structures suggest a common quinone methide-stabilizing function in dirigent proteins and proteins with dirigent-like domains.

The biochemical activities of dirigent proteins (DPs) give rise to distinct complex classes of plant phenolics. DPs apparently began to emerge during the aquatic-to-land transition, with phylogenetic analyses revealing the presence of numerous DP subfamilies in the plant kingdom. The vast majority (>95%) of DPs in these large multigene families still await discovery of their biochemical functions. Here, we elucidated the 3D structures of two pterocarpan-forming proteins with dirigent-like domains. Both proteins stereospecifically convert distinct diastereomeric chiral isoflavonoid precursors to the chiral pterocarpans, (-)- and (+)-medicarpin, respectively. Their 3D structures enabled comparisons with stereoselective lignan- and aromatic terpenoid-forming DP orthologs. Each protein provides entry into diverse plant natural products classes, and our experiments suggest a common biochemical mechanism in binding and stabilizing distinct plant phenol-derived mono- and bis-quinone methide intermediates during different C-C and C-O bond-forming processes. These observations provide key insights into both their appearance and functional diversification of DPs during land plant evolution/adaptation. The proposed biochemical mechanisms based on our findings provide important clues to how additional physiological roles for DPs and proteins harboring dirigent-like domains can now be rationally and systematically identified.

Block Design with Common Reference Samples Enables Robust Large-Scale Label-Free Quantitative Proteome Profiling.

Label-free quantitative proteomics has become an increasingly popular tool for profiling global protein abundances. However, one major limitation is the potential performance drift of the LC-MS platform over time, which, in turn, limits its utility for analyzing large-scale sample sets. To address this, we introduce an experimental and data analysis scheme based on a block design with common references within each block for enabling large-scale label-free quantification. In this scheme, a large number of samples (e.g., >100 samples) are analyzed in smaller and more manageable blocks, minimizing instrument drift and variability within individual blocks. Each designated block also contains common reference samples (e.g., controls) for normalization across all blocks. We demonstrated the robustness of this approach by profiling the proteome response of human macrophage THP-1 cells to 11 engineered nanomaterials at two different doses. A total of 116 samples were analyzed in six blocks, yielding an average coverage of 4500 proteins per sample. Following a common reference-based correction, 2537 proteins were quantified with high reproducibility without any imputation of missing values from 116 data sets. The data revealed the consistent quantification of proteins across all six blocks, as illustrated by the highly consistent abundances of house-keeping proteins in all samples and the high levels of correlation among samples from different blocks. The data also demonstrated that label-free quantification is robust and accurate enough to quantify even very subtle abundance changes as well as large fold-changes. Our streamlined workflow is easy to implement and can be readily adapted to other large cohort studies for reproducible label-free proteome quantification.

A proteome-wide assessment of the oxidative stress paradigm for metal and metal-oxide nanomaterials in human macrophages.

Responsible implementation of engineered nanomaterials (ENMs) into commercial applications is an important societal issue, driving demand for new approaches for rapid and comprehensive evaluation of their bioactivity and safety. An essential part of any research focused on identifying potential hazards of ENMs is the appropriate selection of biological endpoints to evaluate. Herein, we use a tiered strategy employing both targeted biological assays and untargeted quantitative proteomics to elucidate the biological responses of human THP-1 derived macrophages across a library of metal/metal oxide ENMs, raised as priority ENMs for investigation by NIEHS's Nanomaterial Health Implications Research (NHIR) program. Our results show that quantitative cellular proteome profiles readily distinguish ENM types based on their cytotoxic potential according to induction of biological processes and pathways involved in the cellular antioxidant response, TCA cycle, oxidative stress, endoplasmic reticulum stress, and immune responses as major processes impacted. Interestingly, bioinformatics analysis of differentially expressed proteins also revealed new biological processes that were influenced by all ENMs independent of their cytotoxic potential. These included biological processes that were previously implicated as mechanisms cells employ as adaptive responses to low levels of oxidative stress, including cell adhesion, protein translation and protein targeting. Unsupervised clustering revealed the most striking proteome changes that differentiated ENM classes highlight a small subset of proteins involved in the oxidative stress response (HMOX1), protein chaperone functions (HS71B, DNJB1), and autophagy (SQSTM), providing a potential new panel of markers of ENM-induced cellular stress. To our knowledge, the results represent the most comprehensive profiling of the biological responses to a library of ENMs conducted using quantitative mass spectrometry-based proteomics. The results provide a basis to identify the patterns of a diverse set of cellular pathways and biological processes impacted by ENM exposure in an important immune cell type, laying the foundation for multivariate, pathway-level structure activity assessments of ENMs in the future.

Time-dependent behavioral data from zebrafish reveals novel signatures of chemical toxicity using point of departure analysis.

High-content imaging of larval zebrafish behavior can be used as a screening approach to rapidly evaluate the relative potential for chemicals to cause toxicity. However, most statistical methods applied to these data transform movement values to incidence-based "hits" and calculate lowest effect levels (LELs), which loses individual fish resolution of behavior and defies hazard ranking due to reliance on applied dose levels. We developed a parallelizable workflow to calculate benchmark dose (BMD) values from dynamic, high-content zebrafish behavior data that scales for high-throughput chemical screening. To capture the zebrafish movement response from light to dark stimulus, we summarized time-dependent data using both area under the curve and the immediate change at the transition point into two novel metrics that characterized abnormal behavior as a function of chemical concentration. The BMD workflow was applied to calculate BMD10 values of 1,060 ToxCast chemicals for 24 zebrafish endpoints, including behavior, mortality and morphology. The BMD10 values provided better precision and separation than LELs for clustering chemicals since they were derived from models that best-fit their concentration-response curves. Analysis of BMD10 values revealed behavioral signatures as the most sensitive endpoints. High concordance in chemical activity was observed between ToxCast in vitro data and zebrafish in vivo behavioral data, however ToxPi analysis indicated that rankings based on in vitro data were not a reliable predictor of in vivo rankings for lower potency chemicals. This analysis method will enable the use of high-content zebrafish behavioral screening data for BMD analysis in toxicological hazard assessment.

Evaluation of In Silico Multifeature Libraries for Providing Evidence for the Presence of Small Molecules in Synthetic Blinded Samples.

The current gold standard for unambiguous molecular identification in metabolomics analysis is comparing two or more orthogonal properties from the analysis of authentic reference materials (standards) to experimental data acquired in the same laboratory with the same analytical methods. This represents a significant limitation for comprehensive chemical identification of small molecules in complex samples. The process is time consuming and costly, and the majority of molecules are not yet represented by standards. Thus, there is a need to assemble evidence for the presence of small molecules in complex samples through the use of libraries containing calculated chemical properties. To address this need, we developed a Multi-Attribute Matching Engine (MAME) and a library derived in part from our in silico chemical library engine (ISiCLE). Here, we describe an initial evaluation of these methods in a blinded analysis of synthetic chemical mixtures as part of the U.S. Environmental Protection Agency's (EPA) Non-Targeted Analysis Collaborative Trial (ENTACT, Phase 1). For molecules in all mixtures, the initial blinded false negative rate (FNR), false discovery rate (FDR), and accuracy were 57%, 77%, and 91%, respectively. For high evidence scores, the FDR was 35%. After unblinding of the sample compositions, we optimized the scoring parameters to better exploit the available evidence and increased the accuracy for molecules suspected as present. The final FNR, FDR, and accuracy were 67%, 53%, and 96%, respectively. For high evidence scores, the FDR was 10%. This study demonstrates that multiattribute matching methods in conjunction with in silico libraries may one day enable reduced reliance on experimentally derived libraries for building evidence for the presence of molecules in complex samples.

An integrated experimental-computational approach for predicting virulence in New Zealand white rabbits and humans following inhalation exposure to Bacillus anthracis spores.

Inhalation of Bacillus anthracis spores can lead to an anthrax infection that can be fatal. Previously published mathematical models have extrapolated kinetic rates associated with bacterial growth in New Zealand White (NZW) rabbits to humans, but to date, actual measurements of the underlying processes associated with anthrax virulence between species have not been conducted. To address this knowledge gap, we have quantified species-specific rate constants associated with germination, proliferation, and immune cell inactivation of B. anthracis Sterne using an in vitro test platform that includes primary lung epithelial and immune cells. The generated data was then used to develop a physiologically based biokinetic model (PBBK) which quantitatively compares bacterial growth and mean time to death under lethal conditions in rabbits and humans. Simulations based upon our in vitro data and previously published in vivo data from rabbits indicate that disease progression is likely to be faster in humans than in NZW rabbits under comparable total deposited dose conditions. With the computational framework established, PBBK parameters can now be refined using experimental data for lethal B. anthracis strains (e.g. Ames) under identical conditions in future studies. The PBBK model can also be linked to existing aerosol dosimetry models that account for species-specific differences in aerosol deposition patterns to further improve the human health risk assessment of inhalation anthrax.

ISiCLE: A Quantum Chemistry Pipeline for Establishing in Silico Collision Cross Section Libraries.

High-throughput, comprehensive, and confident identifications of metabolites and other chemicals in biological and environmental samples will revolutionize our understanding of the role these chemically diverse molecules play in biological systems. Despite recent technological advances, metabolomics studies still result in the detection of a disproportionate number of features that cannot be confidently assigned to a chemical structure. This inadequacy is driven by the single most significant limitation in metabolomics, the reliance on reference libraries constructed by analysis of authentic reference materials with limited commercial availability. To this end, we have developed the in silico chemical library engine (ISiCLE), a high-performance computing-friendly cheminformatics workflow for generating libraries of chemical properties. In the instantiation described here, we predict probable three-dimensional molecular conformers (i.e., conformational isomers) using chemical identifiers as input, from which collision cross sections (CCS) are derived. The approach employs first-principles simulation, distinguished by the use of molecular dynamics, quantum chemistry, and ion mobility calculations, to generate structures and chemical property libraries, all without training data. Importantly, optimization of ISiCLE included a refactoring of the popular MOBCAL code for trajectory-based mobility calculations, improving its computational efficiency by over 2 orders of magnitude. Calculated CCS values were validated against 1983 experimentally measured CCS values and compared to previously reported CCS calculation approaches. Average calculated CCS error for the validation set is 3.2% using standard parameters, outperforming other density functional theory (DFT)-based methods and machine learning methods (e.g., MetCCS). An online database is introduced for sharing both calculated and experimental CCS values ( metabolomics.pnnl.gov ), initially including a CCS library with over 1 million entries. Finally, three successful applications of molecule characterization using calculated CCS are described, including providing evidence for the presence of an environmental degradation product, the separation of molecular isomers, and an initial characterization of complex blinded mixtures of exposure chemicals. This work represents a method to address the limitations of small molecule identification and offers an alternative to generating chemical identification libraries experimentally by analyzing authentic reference materials. All code is available at github.com/pnnl .

Isolation of Tryptanthrin and Reassessment of Evidence for Its Isobaric Isostere Wrightiadione in Plants of the Wrightia Genus.

A series of Wrightia hanleyi extracts was screened for activity against Mycobacterium tuberculosis H37Rv. One active fraction contained a compound that initially appeared to be either the isoflavonoid wrightiadione or the alkaloid tryptanthrin, both of which have been previously reported in other Wrightia species. Characterization by NMR and MS, as well as evaluation of the literature describing these compounds, led to the conclusion that wrightiadione (1) was misidentified in the first report of its isolation from W. tomentosa in 1992 and again in 2015 when reported in W. pubescens and W. religiosa. Instead, the molecule described in these reports and in the present work is almost certainly the isobaric (same nominal mass) and isosteric (same number of atoms, valency, and shape) tryptanthrin (2), a well-known quinazolinone alkaloid found in a variety of plants including Wrightia species. Tryptanthrin (2) is also accessible synthetically via several routes and has been thoroughly characterized. Wrightiadione (1) has been synthesized and characterized and may have useful biological activity; however, this compound can no longer be said to be known to exist in Nature. To our knowledge, this misidentification of wrightiadione (1) has heretofore been unrecognized.

All that is silver is not toxic: silver ion and particle kinetics reveals the role of silver ion aging and dosimetry on the toxicity of silver nanoparticles.

BACKGROUND: When suspended in cell culture medium, nano-objects composed of soluble metals such as silver can dissolve resulting in ion formation, altered particle properties (e.g. mass, morphology, etc.), and modulated cellular dose. Cultured cells are exposed not just to nanoparticles but to a complex, dynamic mixture of altered nanoparticles, unbound ions, and ion-ligand complexes. Here, three different cell types (RAW 264.7 macrophages and bone marrow derived macrophages from wild-type C57BL/6 J mice and Scavenger Receptor A deficient (SR-A(-/-)) mice) were exposed to 20 and 110 nm silver nanoparticles, and RAW 264.7 cells were exposed to freshly mixed silver ions, aged silver ions (ions incubated in cell culture medium), and ions formed from nanoparticle dissolution. The In Vitro Sedimentation, Diffusion, Dissolution, and Dosimetry Model (ISD3) was used to predict dose metrics for each exposure scenario. RESULTS: Silver nanoparticles, freshly mixed ions, and ions from nanoparticle dissolution were toxic, while aged ions were not toxic. Macrophages from SR-A(-/-) mice did not take up 20 nm silver nanoparticles as well as wild-types but demonstrated no differences in silver levels after exposure to 110 nm nanoparticles. Dose response modeling with ISD3 predicted dose metrics suggest that amount of ions in cells and area under the curve (AUC) of ion amount in cells are the most predictive of cell viability after nanoparticle and combined nanoparticle/dissolution-formed-ions exposures, respectively. CONCLUSIONS: Results of this study suggest that the unbound silver cation is the ultimate toxicant, and ions formed extracellularly drive toxicity after exposure to nanoparticles. Applying computational modeling (ISD3) to better understand dose metrics for soluble nanoparticles allows for better interpretation of in vitro hazard assessments.

Efficient discrimination of natural stereoisomers of chicoric acid, an HIV-1 integrase inhibitor.

Plants from the Asteraceae family are known to contain a wide spectrum of phytochemicals with various nutraceutical properties. One important phytochemical, chicoric acid (CA), is reported to exist in plants, such as Sonchus oleraceus and Bidens pilosa, as stereoisomers. These CA molecules occur either as the naturally abundant RR-chicoric acid (RR-CA), or the less abundant RS-chicoric acid (RS-CA), also known as meso-chicoric acid. To date, little is known about the biological activity of RS-CA, but there is evidence of its anti-human immunodeficiency virus (HIV) properties. In this study, a reliable analytical method was developed to distinguish between the two stereoisomers detected in S. oleraceus and B. pilosa. For structure identification and characterization of CA molecules, liquid chromatography-mass spectrometry (LC-MS) was used in combination with ultraviolet radiation (UV)-induced geometrical isomerization, molecular dynamics (MD) simulations, and density functional theory (DFT) models. Optimized structures from DFT calculations were used for docking studies against the HIV-1 integrase enzyme. Different retention times on the reverse phase chromatograms revealed that the plants produce two different CA stereoisomers: S. oleraceus produced the RR-CA isomer, while B. pilosa produced the RS-CA isomer. DFT results demonstrated the RR-CA molecule was more stable than RS-CA due to the stabilizing force of intra-molecular hydrogen bonding. Differences in the HIV-1 integrase enzyme binding modes were observed, with the RR-CA being a more potent inhibitor than the RS-CA molecule. The results highlight the significance of plant metabolite structural complexity from both chemical and biological perspectives. Furthermore, the study demonstrates that induced-formation of geometrical isomers, in combination with the predictive ability of DFT models and the resolving power of the LC-MS, can be exploited to distinguish structurally closely related compounds, such as stereoisomers.

Explicit ions/implicit water generalized Born model for nucleic acids.

The ion atmosphere around highly charged nucleic acid molecules plays a significant role in their dynamics, structure, and interactions. Here we utilized the implicit solvent framework to develop a model for the explicit treatment of ions interacting with nucleic acid molecules. The proposed explicit ions/implicit water model is based on a significantly modified generalized Born (GB) model and utilizes a non-standard approach to define the solute/solvent dielectric boundary. Specifically, the model includes modifications to the GB interaction terms for the case of multiple interacting solutes-disconnected dielectric boundary around the solute-ion or ion-ion pairs. A fully analytical description of all energy components for charge-charge interactions is provided. The effectiveness of the approach is demonstrated by calculating the potential of mean force for Na+-Cl- ion pair and by carrying out a set of Monte Carlo (MC) simulations of mono- and trivalent ions interacting with DNA and RNA duplexes. The monovalent (Na+) and trivalent (CoHex3+) counterion distributions predicted by the model are in close quantitative agreement with all-atom explicit water molecular dynamics simulations used as reference. Expressed in the units of energy, the maximum deviations of local ion concentrations from the reference are within k B T. The proposed explicit ions/implicit water GB model is able to resolve subtle features and differences of CoHex distributions around DNA and RNA duplexes. These features include preferential CoHex binding inside the major groove of the RNA duplex, in contrast to CoHex biding at the "external" surface of the sugar-phosphate backbone of the DNA duplex; these differences in the counterion binding patters were earlier shown to be responsible for the observed drastic differences in condensation propensities between short DNA and RNA duplexes. MC simulations of CoHex ions interacting with the homopolymeric poly(dA·dT) DNA duplex with modified (de-methylated) and native thymine bases are used to explore the physics behind CoHex-thymine interactions. The simulations suggest that the ion desolvation penalty due to proximity to the low dielectric volume of the methyl group can contribute significantly to CoHex-thymine interactions. Compared to the steric repulsion between the ion and the methyl group, the desolvation penalty interaction has a longer range and may be important to consider in the context of methylation effects on DNA condensation.

An automated framework for NMR chemical shift calculations of small organic molecules.

When using nuclear magnetic resonance (NMR) to assist in chemical identification in complex samples, researchers commonly rely on databases for chemical shift spectra. However, authentic standards are typically depended upon to build libraries experimentally. Considering complex biological samples, such as blood and soil, the entirety of NMR spectra required for all possible compounds would be infeasible to ascertain due to limitations of available standards and experimental processing time. As an alternative, we introduce the in silico Chemical Library Engine (ISiCLE) NMR chemical shift module to accurately and automatically calculate NMR chemical shifts of small organic molecules through use of quantum chemical calculations. ISiCLE performs density functional theory (DFT)-based calculations for predicting chemical properties-specifically NMR chemical shifts in this manuscript-via the open source, high-performance computational chemistry software, NWChem. ISiCLE calculates the NMR chemical shifts of sets of molecules using any available combination of DFT method, solvent, and NMR-active nuclei, using both user-selected reference compounds and/or linear regression methods. Calculated NMR chemical shifts are provided to the user for each molecule, along with comparisons with respect to a number of metrics commonly used in the literature. Here, we demonstrate ISiCLE using a set of 312 molecules, ranging in size up to 90 carbon atoms. For each, calculation of NMR chemical shifts have been performed with 8 different levels of DFT theory, and with solvation effects using the implicit solvent Conductor-like Screening Model. The DFT method dependence of the calculated chemical shifts have been systematically investigated through benchmarking and subsequently compared to experimental data available in the literature. Furthermore, ISiCLE has been applied to a set of 80 methylcyclohexane conformers, combined via Boltzmann weighting and compared to experimental values. We demonstrate that our protocol shows promise in the automation of chemical shift calculations and, ultimately, the expansion of chemical shift libraries.

Integration among databases and data sets to support productive nanotechnology: Challenges and recommendations.

Many groups within the broad field of nanoinformatics are already developing data repositories and analytical tools driven by their individual organizational goals. Integrating these data resources across disciplines and with non-nanotechnology resources can support multiple objectives by enabling the reuse of the same information. Integration can also serve as the impetus for novel scientific discoveries by providing the framework to support deeper data analyses. This article discusses current data integration practices in nanoinformatics and in comparable mature fields, and nanotechnology-specific challenges impacting data integration. Based on results from a nanoinformatics-community-wide survey, recommendations for achieving integration of existing operational nanotechnology resources are presented. Nanotechnology-specific data integration challenges, if effectively resolved, can foster the application and validation of nanotechnology within and across disciplines. This paper is one of a series of articles by the Nanomaterial Data Curation Initiative that address data issues such as data curation workflows, data completeness and quality, curator responsibilities, and metadata.

The Role of Specific Chemokines in the Amelioration of Colitis by Appendicitis and Appendectomy.

The appendix contains abundant lymphoid tissue and is constantly exposed to gut flora. When completed at a young age, appendicitis followed by appendectomy (AA) prevents or significantly ameliorates Inflammatory Bowel Diseases (IBDs) in later life. Inflammatory bowel disease comprises Crohn's disease and ulcerative colitis. Our murine AA model is the only existing experimental model of AA. In our unique model, AA performed in the most proximal colon limits colitis pathology in the most distal colon by curbing T-helper 17 cell activity, diminishing autophagy, modulating interferon activity-associated molecules, and suppressing endothelin vaso-activity-mediated immunopathology. In the research presented in this paper, we have examined the role of chemokines in colitis pathology with our murine AA model. Chemokines are a family of small cytokines with four conserved cysteine residues. Chemokines induce chemotaxis in adjacent cells with corresponding receptors. All 40 known chemokine genes and 24 chemokine receptor genes were examined for gene expression levels in distal colons three days post-AA and 28 days post-AA. At 28 days post-AA, the chemokine gene CCL5 was significantly upregulated. Furthermore, Gene Set Enrichment Analysis (GSEA) showed upregulation of seven CCL5-associated gene-sets 28 days post-AA in contrast to just one gene-set downregulated at the same time-point. The chemokine gene CXCL11 was significantly upregulated three days post-AA and 28 days post-AA. Evaluation using GSEA showed upregulation of six CXCL11-associated gene sets but no downregulation of any gene set. At 28 days post-AA, CCL17 gene expression was significantly downregulated. There was no expression of any chemokine receptor gene three days post-AA, but CCR10 was the only chemokine receptor gene that displayed differential gene expression (upregulation) 28 days post-AA. No CCR10-associated gene set was upregulated in GSEA in contrast to one downregulated gene set. Our analysis resulted in identifying three new therapeutic targets towards ameliorating colitis: CCL5, CXCL11, and CCL17. While CCL5 and CXCL11 are good therapeutic chemokine candidates to be exogenously administered, CCL17 is a good candidate chemokine to competitively inhibit or limit colitis pathology.

Atomic Radius and Charge Parameter Uncertainty in Biomolecular Solvation Energy Calculations.

Atomic radii and charges are two major parameters used in implicit solvent electrostatics and energy calculations. The optimization problem for charges and radii is underdetermined, leading to uncertainty in the values of these parameters and in the results of solvation energy calculations using these parameters. This paper presents a new method for quantifying this uncertainty in implicit solvation calculations of small molecules using surrogate models based on generalized polynomial chaos (gPC) expansions. There are relatively few atom types used to specify radii parameters in implicit solvation calculations; therefore, surrogate models for these low-dimensional spaces could be constructed using least-squares fitting. However, there are many more types of atomic charges; therefore, construction of surrogate models for the charge parameter space requires compressed sensing combined with an iterative rotation method to enhance problem sparsity. We demonstrate the application of the method by presenting results for the uncertainties in small molecule solvation energies based on these approaches. The method presented in this paper is a promising approach for efficiently quantifying uncertainty in a wide range of force field parametrization problems, including those beyond continuum solvation calculations. The intent of this study is to provide a way for developers of implicit solvent model parameter sets to understand the sensitivity of their target properties (solvation energy) on underlying choices for solute radius and charge parameters.