An Integrated Data-Driven Strategy for Safe-by-Design Nanoparticles: The FP7 MODERN Project.

The development and implementation of safe-by-design strategies is key for the safe development of future generations of nanotechnology enabled products. The safety testing of the huge variety of nanomaterials that can be synthetized is unfeasible due to time and cost constraints. Computational modeling facilitates the implementation of alternative testing strategies in a time and cost effective way. The development of predictive nanotoxicology models requires the use of high quality experimental data on the structure, physicochemical properties and bioactivity of nanomaterials. The FP7 Project MODERN has developed and evaluated the main components of a computational framework for the evaluation of the environmental and health impacts of nanoparticles. This chapter describes each of the elements of the framework including aspects related to data generation, management and integration; development of nanodescriptors; establishment of nanostructure-activity relationships; identification of nanoparticle categories; hazard ranking and risk assessment.

Prioritization of in silico models and molecular descriptors for the assessment of ready biodegradability.

Ready biodegradability is a key property for evaluating the long-term effects of chemicals on the environment and human health. As such, it is used as a screening test for the assessment of persistent, bioaccumulative and toxic substances. Regulators encourage the use of non-testing methods, such as in silico models, to save money and time. A dataset of 757 chemicals was collected to assess the performance of four freely available in silico models that predict ready biodegradability. They were applied to develop a new consensus method that prioritizes the use of each individual model according to its performance on chemical subsets driven by the presence or absence of different molecular descriptors. This consensus method was capable of almost eliminating unpredictable chemicals, while the performance of combined models was substantially improved with respect to that of the individual models.

Optimal descriptor as a translator of eclectic data into prediction of cytotoxicity for metal oxide nanoparticles under different conditions.

The Monte Carlo technique has been used to build up quantitative structure-activity relationships (QSARs) for prediction of dark cytotoxicity and photo-induced cytotoxicity of metal oxide nanoparticles to bacteria Escherichia coli (minus logarithm of lethal concentration for 50% bacteria pLC50, LC50 in mol/L). The representation of nanoparticles include (i) in the case of the dark cytotoxicity a simplified molecular input-line entry system (SMILES), and (ii) in the case of photo-induced cytotoxicity a SMILES plus symbol '^'. The predictability of the approach is checked up with six random distributions of available data into the visible training and calibration sets, and invisible validation set. The statistical characteristics of these models are correlation coefficient 0.90-0.94 (training set) and 0.73-0.98 (validation set).

In silico analysis of nanomaterials hazard and risk.

Because a variety of human-related activities, engineer-ed nanoparticles (ENMs) may be released to various environmental media and may cross environmental boundaries, and thus will be found in most media. Therefore, the potential environmental impacts of ENMs must be assessed from a multimedia perspective and with an integrated risk management approach that considers rapid developments and increasing use of new nanomaterials. Accordingly, this Account presents a rational process for the integration of in silico ENM toxicity and fate and transport analyses for environmental impact assessment. This approach requires knowledge of ENM toxicity and environmental exposure concentrations. Considering the large number of current different types of ENMs and that those numbers are likely to increase, there is an urgent need to accelerate the evaluation of their toxicity and the assessment of their potential distribution in the environment. Developments in high throughput screening (HTS) are now enabling the rapid generation of large data sets for ENM toxicity assessment. However, these analyses require the establishment of reliable toxicity metrics, especially when HTS includes data from multiple assays, cell lines, or organisms. Establishing toxicity metrics with HTS data requires advanced data processing techniques in order to clearly identify significant biological effects associated with exposure to ENMs. HTS data can form the basis for developing and validating in silico toxicity models (e.g., quantitative structure-activity relationships) and for generating data-driven hypotheses to aid in establishing and/or validating possible toxicity mechanisms. To correlate the toxicity of ENMs with their physicochemical properties, researchers will need to develop quantitative structure-activity relationships for nanomaterials (i.e., nano-SARs). However, as nano-SARs are applied in regulatory applications, researchers must consider their applicability and the acceptance level of false positive relative to false negative predictions and the reliability of toxicity data. To establish the environmental impact of ENMs identified as toxic, researchers will need to estimate the potential level of environmental exposure concentration of ENMs in the various media such as air, water, soil, and vegetation. When environmental monitoring data are not available, models of ENMs fate and transport (at various levels of complexity) serve as alternative approaches for estimating exposure concentrations. Risk management decisions regarding the manufacturing, use, and environmental regulations of ENMs would clearly benefit from both the assessment of potential ENMs exposure concentrations and suitable toxicity metrics. The decision process should consider the totality of available information: quantitative and qualitative data and the analysis of nanomaterials toxicity, and fate and transport behavior in the environment. Effective decision-making to address the potential impacts of nanomaterials will require considerations of the relevant environmental, ecological, technological, economic, and sociopolitical factors affecting the complete lifecycle of nanomaterials, while accounting for data and modeling uncertainties. Accordingly, researchers will need to establish standardized data management and analysis tools through nanoinformatics as a basis for the development of rational decision tools.

Association rule mining of cellular responses induced by metal and metal oxide nanoparticles.

Relationships among fourteen different biological responses (including ten signaling pathway activities and four cytotoxicity effects) of murine macrophage (RAW264.7) and bronchial epithelial (BEAS-2B) cells exposed to six metal and metal oxide nanoparticles (NPs) were analyzed using both statistical and data mining approaches. Both the pathway activities and cytotoxicity effects were assessed using high-throughput screening (HTS) over an exposure period of up to 24 h and concentration range of 0.39-200 mg L(-1). HTS data were processed by outlier removal, normalization, and hit-identification (for significantly regulated cellular responses) to arrive at reliable multiparametric bioactivity profiles for the NPs. Association rule mining was then applied to the bioactivity profiles followed by a pruning process to remove redundant rules. The non-redundant association rules indicated that "significant regulation" of one or more cellular responses implies regulation of other (associated) cellular response types. Pairwise correlation analysis (via Pearson's χ(2) test) and self-organizing map clustering of the different cellular response types indicated consistency with the identified non-redundant association rules. Furthermore, in order to explore the potential use of association rules as a tool for data-driven hypothesis generation, specific pathway activity experiments were carried out for ZnO NPs. The experimental results confirmed the association rule identified for the p53 pathway and mitochondrial superoxide levels (via MitoSox reagent) and further revealed that blocking of the transcriptional activity of p53 lowered the MitoSox signal. The present approach of using association rule mining for data-driven hypothesis generation has important implications for streamlining multi-parameter HTS assays, improving the understanding of NP toxicity mechanisms, and selection of endpoints for the development of nanomaterial structure-activity relationships.

A US perspective on closing the carbon cycle to defossilize difficult-to-electrify segments of our economy.

Electrification to reduce or eliminate greenhouse gas emissions is essential to mitigate climate change. However, a substantial portion of our manufacturing and transportation infrastructure will be difficult to electrify and/or will continue to use carbon as a key component, including areas in aviation, heavy-duty and marine transportation, and the chemical industry. In this Roadmap, we explore how multidisciplinary approaches will enable us to close the carbon cycle and create a circular economy by defossilizing these difficult-to-electrify areas and those that will continue to need carbon. We discuss two approaches for this: developing carbon alternatives and improving our ability to reuse carbon, enabled by separations. Furthermore, we posit that co-design and use-driven fundamental science are essential to reach aggressive greenhouse gas reduction targets.

Nano-SAR development for bioactivity of nanoparticles with considerations of decision boundaries.

The development of classification nano-structure-activity Relationships (nano-SARs) of nanoparticle (NP) bioactivity is presented with the aim of demonstrating the integration of multiparametric toxicity/bioactivity assays to arrive at statistically meaningful class definitions (i.e., bioactivity/inactivity endpoints), as well as the implications of nano-SAR applicability domains and decision boundaries. Nano-SARs are constructed based on a dataset of 44 iron oxide core nanoparticles (NPs), used in molecular imaging and nano-sensing, containing bioactivity profiles for four cell types and four different assays. Class definitions are developed on the basis of 'hit' (i.e., significant bioactivity) identification analysis and self-organizing map based consensus clustering; these class definitions enable construction of nano-SARs of a high classification accuracy (>78%) with different NP descriptor combinations that include primary size, spin-lattice and spin-spin relaxivities, and zeta potentials. Analysis of the nano-SAR performance for different class definitions suggests that H4 (i.e., class with at least four hits) is a reasonable endpoint (from a 'regulatory' viewpoint) for keeping the level of false negatives (i.e., incorrect labeling of bioactive NPs as inactive) low. The establishment of a quantitative nano-SAR applicability domain is demonstrated, making use of a probability density with the H4 class definition and naive Bayesian classifier (NBC) model (with spin-lattice relaxivity and zeta potential as descriptors). Decision boundaries are determined for the above H4/NBC nano-SAR for different acceptance levels of false negative to false positive predictions, illustrating a practical approach that may assist in regulatory decision making with a consideration of reducing the likelihood of identifying bioactive NPs as being inactive.

A comparative survey of chemistry-driven in silico methods to identify hazardous substances under REACH.

This paper presents an inventory of in silico screening tools to identify substance properties of concern under the European chemicals' legislation REACH. The objective is to support the selection and implementation of appropriate tools as building blocks within integrated testing strategies (ITS). The relevant concerns addressed are persistence, bioaccumulation potential, acute and long-term aquatic toxicity, PBT/vPvB properties ((very) persistent, (very) bioaccumulative, toxic), CMR (carcinogenicity, mutagenicity, reproductive toxicity), endocrine disruption and skin sensitisation. The inventory offers a comparative evaluation of methods with respect to the underlying algorithms (how does the method work?) and the applicability domains (when does the method work?) as well as their limitations (when does the method not work?). The inventory explicitly addresses the reliability of predictions of different in silico models for diverse chemicals by applicability domain considerations. The confidence in predictions can be greatly improved by consensus modelling that allows for taking conflicting results into account. The inventory is complemented by a brief discussion of socio-economic tools for assessing the potential efficiency gains of using in silico methods compared to traditional in vivo testing of chemical hazards.

Molecular structural dataset of lignin macromolecule elucidating experimental structural compositions.

Lignin is one of the most abundant biopolymers in nature and has great potential to be transformed into high-value chemicals. However, the limited availability of molecular structure data hinders its potential industrial applications. Herein, we present the Lignin Structural (LGS) Dataset that includes the molecular structure of milled wood lignin focusing on two major monomeric units (coniferyl and syringyl), and the six most common interunit linkages (phenylpropane ß-aryl ether, resinol, phenylcoumaran, biphenyl, dibenzodioxocin, and diaryl ether). The dataset constitutes a unique resource that covers a part of lignin's chemical space characterized by polymer chains with lengths in the range of 3 to 25 monomer units. Structural data were generated using a sequence-controlled polymer generation approach that was calibrated to match experimental lignin properties. The LGS dataset includes 60 K newly generated lignin structures that match with high accuracy (~90%) the experimentally determined structural compositions available in the literature. The LGS dataset is a valuable resource to advance lignin chemistry research, including computational simulation approaches and predictive modelling.

Preparing for the next pandemic via transfer learning from existing diseases with hierarchical multi-modal BERT: a study on COVID-19 outcome prediction.

Developing prediction models for emerging infectious diseases from relatively small numbers of cases is a critical need for improving pandemic preparedness. Using COVID-19 as an exemplar, we propose a transfer learning methodology for developing predictive models from multi-modal electronic healthcare records by leveraging information from more prevalent diseases with shared clinical characteristics. Our novel hierarchical, multi-modal model ([Formula: see text]) integrates baseline risk factors from the natural language processing of clinical notes at admission, time-series measurements of biomarkers obtained from laboratory tests, and discrete diagnostic, procedure and drug codes. We demonstrate the alignment of [Formula: see text]'s predictions with well-established clinical knowledge about COVID-19 through univariate and multivariate risk factor driven sub-cohort analysis. [Formula: see text]'s superior performance over state-of-the-art methods shows that leveraging patient data across modalities and transferring prior knowledge from similar disorders is critical for accurate prediction of patient outcomes, and this approach may serve as an important tool in the early response to future pandemics.

Classification NanoSAR development for cytotoxicity of metal oxide nanoparticles.

A classification-based cytotoxicity nanostructure-activity relationship (nanoSAR) is presented based on a set of nine metal oxide nanoparticles to which transformed bronchial epithelial cells (BEAS-2B) were exposed over a range of concentrations (0.375-200 mg L(-1) ) and exposure times up to 24 h. The nanoSAR is developed using cytotoxicity data from a high-throughput screening assay that was processed to identify and label toxic (in terms of the propidium iodide uptake of BEAS-2B cells) versus nontoxic events relative to an unexposed control cell population. Starting with a set of fourteen intuitive but fundamental physicochemical nanoSAR input parameters, a number of models were identified which had a classification accuracy above 95%. The best-performing model had a 100% classification accuracy in both internal and external validations. This model is based on three descriptors: atomization energy of the metal oxide, period of the nanoparticle metal, and nanoparticle primary size, in addition to nanoparticle volume fraction (in solution). Notwithstanding the success of the present modeling approach with a relatively small nanoparticle library, it is important to recognize that a significantly larger data set would be needed in order to expand the applicability domain and increase the confidence and reliability of data-driven nanoSARs.

Analysis of nanoparticle agglomeration in aqueous suspensions via constant-number Monte Carlo simulation.

A constant-number direct simulation Monte Carlo (DSMC) model was developed for the analysis of nanoparticle (NP) agglomeration in aqueous suspensions. The modeling approach, based on the "particles in a box" simulation method, considered both particle agglomeration and gravitational settling. Particle-particle agglomeration probability was determined based on the classical Derjaguin-Landau-Verwey-Overbeek (DLVO) theory and considerations of the collision frequency as impacted by Brownian motion. Model predictions were in reasonable agreement with respect to the particle size distribution and average agglomerate size when compared with dynamic light scattering (DLS) measurements for aqueous TiO(2), CeO(2), and C(60) nanoparticle suspensions over a wide range of pH (3-10) and ionic strength (0.01-156 mM). Simulations also demonstrated, in quantitative agreement with DLS measurements, that nanoparticle agglomerate size increased both with ionic strength and as the solution pH approached the isoelectric point (IEP). The present work suggests that the DSMC modeling approach, along with future use of an extended DLVO theory, has the potential for becoming a practical environmental analysis tool for predicting the agglomeration behavior of aqueous nanoparticle suspensions.

Self-organizing map analysis of toxicity-related cell signaling pathways for metal and metal oxide nanoparticles.

The response of a murine macrophage cell line exposed to a library of seven metal and metal oxide nanoparticles was evaluated via High Throughput Screening (HTS) assay employing luciferase-reporters for ten independent toxicity-related signaling pathways. Similarities of toxicity response among the nanoparticles were identified via Self-Organizing Map (SOM) analysis. This analysis, applied to the HTS data, quantified the significance of the signaling pathway responses (SPRs) of the cell population exposed to nanomaterials relative to a population of untreated cells, using the Strictly Standardized Mean Difference (SSMD). Given the high dimensionality of the data and relatively small data set, the validity of the SOM clusters was established via a consensus clustering technique. Analysis of the SPR signatures revealed two cluster groups corresponding to (i) sublethal pro-inflammatory responses to Al2O3, Au, Ag, SiO2 nanoparticles possibly related to ROS generation, and (ii) lethal genotoxic responses due to exposure to ZnO and Pt nanoparticles at a concentration range of 25-100 µg/mL at 12 h exposure. In addition to identifying and visualizing clusters and quantifying similarity measures, the SOM approach can aid in developing predictive quantitative-structure relations; however, this would require significantly larger data sets generated from combinatorial libraries of engineered nanoparticles.

Visualizing biomolecular electrostatics in virtual reality with UnityMol-APBS.

Virtual reality is a powerful tool with the ability to immerse a user within a completely external environment. This immersion is particularly useful when visualizing and analyzing interactions between small organic molecules, molecular inorganic complexes, and biomolecular systems such as redox proteins and enzymes. A common tool used in the biomedical community to analyze such interactions is the Adaptive Poisson-Boltzmann Solver (APBS) software, which was developed to solve the equations of continuum electrostatics for large biomolecular assemblages. Numerous applications exist for using APBS in the biomedical community including analysis of protein ligand interactions and APBS has enjoyed widespread adoption throughout the biomedical community. Currently, typical use of the full APBS toolset is completed via the command line followed by visualization using a variety of two-dimensional external molecular visualization software. This process has inherent limitations: visualization of three-dimensional objects using a two-dimensional interface masks important information within the depth component. Herein, we have developed a single application, UnityMol-APBS, that provides a dual experience where users can utilize the full range of the APBS toolset, without the use of a command line interface, by use of a simple graphical user interface (GUI) for either a standard desktop or immersive virtual reality experience.

Molecular dynamics simulations of zinc oxide solubility: From bulk down to nanoparticles.

The solubility of metal oxides is one of the key descriptors for the evaluation of their potential toxic effects, both in the bulk form and in nanoparticulated aggregates. Current work presents a new methodology for the in silico assessment of the solubility of metal oxides, which is demonstrated using a well-studied system, ZnO. The calculation of the solubility is based on statistical thermodynamics tools combined with Density Functional Tight Binding theory for the evaluation of the free energy exchange during the dissolution process. Models of small ZnO clusters are used for describing the final dissolved material, since the complete ionic dissolution of ZnO is hindered by the formation of O2- anions in solution, which are highly unstable. Results show very good agreement between the computed solubility values and experimental data for ZnO bulk, up to 0.5 mg L-1 and equivalents of 50 µg L-1 for the free Zn2+ cation in solution. However, the reference model for solid nanoparticles formed by free space nanoparticles can only give a limited quantitative solubility evaluation for ZnO nanoparticles.

Perspectives from the NanoSafety Modelling Cluster on the validation criteria for (Q)SAR models used in nanotechnology.

Nanotechnology and the production of nanomaterials have been expanding rapidly in recent years. Since many types of engineered nanoparticles are suspected to be toxic to living organisms and to have a negative impact on the environment, the process of designing new nanoparticles and their applications must be accompanied by a thorough risk analysis. (Quantitative) Structure-Activity Relationship ([Q]SAR) modelling creates promising options among the available methods for the risk assessment. These in silico models can be used to predict a variety of properties, including the toxicity of newly designed nanoparticles. However, (Q)SAR models must be appropriately validated to ensure the clarity, consistency and reliability of predictions. This paper is a joint initiative from recently completed European research projects focused on developing (Q)SAR methodology for nanomaterials. The aim was to interpret and expand the guidance for the well-known "OECD Principles for the Validation, for Regulatory Purposes, of (Q)SAR Models", with reference to nano-(Q)SAR, and present our opinions on the criteria to be fulfilled for models developed for nanoparticles.

In Silico Design of Optimal Dissolution Kinetics of Fe-Doped ZnO Nanoparticles Results in Cancer-Specific Toxicity in a Preclinical Rodent Model.

Cancer cells have unique but widely varying characteristics that have proven them difficult to be treated by classical therapeutics and calls for novel and selective treatment options. Nanomaterials (NMs) have been shown to display biological effects as a function of their chemical composition, and the extent and exact nature of these effects can vary between different biological environments. Here, ZnO NMs are doped with increasing levels of Fe, which allows to finely tune their dissolution rate resulting in significant differences in their biological behavior on cancer or normal cells. Based on in silico analysis, 2% Fe-doped ZnO NMs are found to be optimal to cause selective cancer cell death, which is confirmed in both cultured cells and syngeneic tumor models, where they also reduce metastasis formation. These results show that upon tuning NM chemical composition, NMs can be designed as a targeted selective anticancer therapy.

CompNanoTox2015: novel perspectives from a European conference on computational nanotoxicology on predictive nanotoxicology.

A first European Conference on Computational Nanotoxicology, CompNanoTox, was held in November 2015 in Benahavís, Spain with the objectives to disseminate and integrate results from the European modeling and database projects (NanoPUZZLES, ModENPTox, PreNanoTox, MembraneNanoPart, MODERN, eNanoMapper and EU COST TD1204 MODENA) as well as to create synergies within the European NanoSafety Cluster. This conference was supported by the COST Action TD1204 MODENA on developing computational methods for toxicological risk assessment of engineered nanoparticles and provided a unique opportunity for cross fertilization among complementary disciplines. The efforts to develop and validate computational models crucially depend on high quality experimental data and relevant assays which will be the basis to identify relevant descriptors. The ambitious overarching goal of this conference was to promote predictive nanotoxicology, which can only be achieved by a close collaboration between the computational scientists (e.g. database experts, modeling experts for structure, (eco) toxicological effects, performance and interaction of nanomaterials) and experimentalists from different areas (in particular toxicologists, biologists, chemists and material scientists, among others). The main outcome and new perspectives of this conference are summarized here.

Use of Quasi-SMILES and Monte Carlo Optimization to Develop Quantitative Feature Property/Activity Relationships (QFPR/QFAR) for Nanomaterials.

The CORAL software (http://www.insilico.eu/coral) has been used to develop quantitative feature-property/activity relationships (QFPRs/QFARs) for the prediction of endpoints related to different categories of nanomaterials. In contrast to previous models built up by using CORAL from a representation of the molecular structure by using simplified molecular input-line entry system (SMILES), the current QFPR/QFARs are based on an integrated representation of acting conditions (i.e., a combination of physicochemical and/or biochemical factors) of nanomaterials via the so-called quasi-SMILES notation. In contrast to traditional quantitative structure - property / activity relationships (QSPRs/QSARs), the new models are able to provide new insight on the conditions of acting of substances (e.g., chemicals and nanomaterials) independently of their molecular structure. The development and validation of the QFPR/QFAR models was carried out following the OECD principles. The statistical quality of models developed from quasi-SMILES is acceptable, with values for the determination coefficient in the range of 0.70 to 0.85 for various endpoints of environmental and human health relevance. Perspectives of the QFPR/QFAR and their interaction and overlapping with traditional QSPR/QSAR are also discussed.

Predicting Cell Association of Surface-Modified Nanoparticles Using Protein Corona Structure - Activity Relationships (PCSAR).

Nanoparticles are likely to interact in real-case application scenarios with mixtures of proteins and biomolecules that will absorb onto their surface forming the so-called protein corona. Information related to the composition of the protein corona and net cell association was collected from literature for a library of surface-modified gold and silver nanoparticles. For each protein in the corona, sequence information was extracted and used to calculate physicochemical properties and statistical descriptors. Data cleaning and preprocessing techniques including statistical analysis and feature selection methods were applied to remove highly correlated, redundant and non-significant features. A weighting technique was applied to construct specific signatures that represent the corona composition for each nanoparticle. Using this basic set of protein descriptors, a new Protein Corona Structure-Activity Relationship (PCSAR) that relates net cell association with the physicochemical descriptors of the proteins that form the corona was developed and validated. The features that resulted from the feature selection were in line with already published literature, and the computational model constructed on these features had a good accuracy (R(2)LOO=0.76 and R(2)LMO(25%)=0.72) and stability, with the advantage that the fingerprints based on physicochemical descriptors were independent of the specific proteins that form the corona.