Next Article in Journal
Surface Plasmon Enhancement of Eu3+ Emission Intensity in LaPO4/Ag Nanoparticles
Next Article in Special Issue
Controlling the Porosity and Biocidal Properties of the Chitosan-Hyaluronate Matrix Hydrogel Nanocomposites by the Addition of 2D Ti3C2Tx MXene
Previous Article in Journal
Autogenous Dentin Graft in Bone Defects after Lower Third Molar Extraction: A Split-Mouth Clinical Trial
Previous Article in Special Issue
Engineering of 2D Ti3C2 MXene Surface Charge and its Influence on Biological Properties
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Materials
Volume 13
Issue 14
10.3390/ma13143083
materials-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Surface-Related Features Responsible for Cytotoxic Behavior of MXenes Layered Materials Predicted with Machine Learning Approach

by
Maciej E. Marchwiany
1,†,
Magdalena Birowska
2,*,†,
Mariusz Popielski
2,
Jacek A. Majewski
2 and
Agnieszka M. Jastrzębska
3
1
Interdisciplinary Centre for Mathematical and Computational Modelling (ICM), University of Warsaw, Pawińskiego 5a, 02-106 Warsaw, Poland
2
Faculty of Physics, University of Warsaw, Pasteura 5, 00-092 Warsaw, Poland
3
Faculty of Materials Science and Engineering, Warsaw University of Technology, Wołoska 141, 02-507 Warsaw, Poland
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Materials 2020, 13(14), 3083; https://doi.org/10.3390/ma13143083
Submission received: 19 June 2020 / Revised: 6 July 2020 / Accepted: 8 July 2020 / Published: 10 July 2020
(This article belongs to the Special Issue New Findings of MXenes: Preparation, Properties and Applications in Biotechnology and Catalysis)
Browse Figures
Versions Notes

Abstract

:
To speed up the implementation of the two-dimensional materials in the development of potential biomedical applications, the toxicological aspects toward human health need to be addressed. Due to time-consuming and expensive analysis, only part of the continuously expanding family of 2D materials can be tested in vitro. The machine learning methods can be used—by extracting new insights from available biological data sets, and provide further guidance for experimental studies. This study identifies the most relevant highly surface-specific features that might be responsible for cytotoxic behavior of 2D materials, especially MXenes. In particular, two factors, namely, the presence of transition metal oxides and lithium atoms on the surface, are identified as cytotoxicity-generating features. The developed machine learning model succeeds in predicting toxicity for other 2D MXenes, previously not tested in vitro, and hence, is able to complement the existing knowledge coming from in vitro studies. Thus, we claim that it might be one of the solutions for reducing the number of toxicological studies needed, and allows for minimizing failures in future biological applications.

Graphical Abstract

1. Introduction

Experimental development of two-dimensional (2D) materials is booming and the industrial applications are not only envisioned but also are becoming a reality. On the other hand, safety verification by systematic and in-depth studies pose a great challenge. The first stage of research into the safety of individual materials are the in vitro studies on mammalian cells. If material is nontoxic, they are moving forward into more advanced studies. In this regard, toxicity appears when the viability of tested cells decreases from 100% to below 70%. This means reduction over 30% which is significantly below the natural 20% variance of survival rate assumed for the living organisms.
One of the most perspective 2D materials that can be utilized in technological applications are the MXenes. During last years their development has risen exponentially. MXenes are defined as early transition metal carbides, nitrides, and carbonitrides. They have received much attention as their unique 2D crystal structure can be easily tuned to produce dramatic improvement in material properties [1]. Therefore, it is not surprising that they have been successfully applied in many fields of materials science and technology [2]. The term ‘MXene’ reflects the unique 2D structure of the material in which the formula Mn+1XnTz perfectly matches the arrangement of its layered features. In this formula—M is early transition metal, X is carbon and/or nitrogen, n = 1, 2, 3, or 4, whereas Tz corresponds to functional groups terminating the surface (e.g., -OH, = -O, -F) [3]. The family of MXenes has expanded rapidly since the discovery of their first representative—the Ti3C2Tz phase in 2011 by Naguib et al. [4]. Note, that the first package of MXenes included only several phases with 19+ successfully synthesized in subsequent years [5]. Since that time, ten years have passed and now researchers are able to predict new MXene phases theoretically [6], and new phases have been successfully synthesized (see Reference [7]).
In the case of MXenes cytotoxicity, the first studies concerned in vitro testing of multilayered Ti3C2Tz MXene and showed a potential threat related to the generation of reactive oxygen species (ROS) [8]. Further studies showed differences in toxicological effects in view of MXenes stoichiometry (i.e., Ti3C2Tz or Ti2CTz) [9]. Moreover, the importance of flake thickness was highlighted by us not only in the case of material stability but also potential toxicity [10]. In addition, cytotoxicity was also examined in respect to the number of layers constituting MXene compounds [11].
In this regard, we have come to the moment that innovative solutions are needed to extract the most promising representatives of MXenes with the highest potential for application and the lowest cytotoxicological threats. It becomes obvious that it is impossible to carry out screening investigations for all MXenes phases in reasonable time and manageable costs. The most time- and cost-consuming analyses are undoubtedly the biological studies, which are also inevitable to push through MXenes applications in industry. What is more, many certification procedures involve verification of the safety of market products containing the claimed nanomaterials. Accordingly, there is a strong demand for theoretical solutions that could overcome the problem of so many complicated analyses.
One such solutions might be the machine learning (ML) procedure. Machine learning has so far proved its applicability for cytotoxicity studies of large number of various chemicals (see Reference [12]) as well as recent predictions of synthesis of various MXenes compounds [13]. In addition, it can be also used to effectively analyse complex surface science data [14]. Thus, we assume that it might be one of the solutions for reducing the number of toxicological studies needed, and allows for minimizing failures in future biological applications. Machine learning studies concerning toxicity of drug, molecules have been carried out previously, by using deep learning and XGBoost [15] approaches, and using Atomic Fingerprints [16,17,18,19,20,21]. However, to the best of our knowledge, there is a lack of toxicological research based on ML methods concerning layered materials, in particular, 2D structures. Thus, we present the first approach of predicting the cytotoxicity of experimentally synthesized MXene compounds.
The aim of this study is to provide the prediction of the potential cytotoxic behavior of MXenes materials based on ML model with some elemental information provided from experiments. We first determine the biological and physico-chemical features that describe each 2D material in relation to the tested cytotoxicity. On this basis, we determine the possible descriptors for predicting the cytotoxic behavior of MXenes compounds. In particular, the surface characteristics, morphology, and structure have been taken into account as inputs for our theoretical model. Next, the Random Forest [22] approach has been applied to identify the most important features, that might have an impact on MXenes cytotoxicity. Then, we apply Principle Component Analysis (PCA) [23] as feature engineering to improve our model. We use the key features to train machine learning models. The models are checked by a 10-fold cross-validation scheme, with their performance measured by accuracy score, data set contains experimentally measured samples provided from our own research studies, as well as from screening experiments. Then, we use this model to predict the potential cytotoxic behavior of 19 experimentally examined MXene compounds not tested in vitro. It is crucial to note, that this study relies on small size of dataset, thus, we include variety of external methods to validate our predictions.
Note, that the theoretical methodology developed here can be further applied to other types of 2D materials, which exhibit diverse surface characteristics. The detailed information about the experimental data and elemental features used in this study are presented in Appendices Appendix A and Appendix B, respectively.

2. Materials and Methods

In order to determine a good quality ML model, it is crucial to determine potential descriptors that can characterize the cytotoxic behavior and identify ML algorithms well suited for a given dataset.

2.1. The Choice of Descriptors

Descriptors, which are representative of the compounds properties play a crucial role in ML, thus, we decided to describe in details the rationale behind the particular choice of the descriptors used in our studies.
At first, we choose only the MXenes compounds tested in vitro which are well experimentally characterised and appear cytotoxic or not. In this regard, the particular MXene exhibits cytotoxicity if it causes the reduction of viability of tested cells in vitro below 70% within the concentration range up to 250 m g L . In addition, the ROS level above 120% at a concentration of 250 m g L is also an indicative of the presence of cytotoxicity. The rest of MXenes appears as non-cytotoxic. Furthermore, we gather and sort the information from material synthesis methodology as well as results from characterisation. This includes the chemical composition of MXene (M, X), surface modification with external compounds, lateral size, thickness, etching agent, delaminating agent, elements and oxides (MxOy) present on the surface which are studied by the X-ray photoelectron spectroscopy (XPS). Based on these considerations, we aim to select particular experimental descriptors listed in Table A2.
Based on the previous experiments [9] that linked the cytotoxic behavior of MXenes with the stoichiometry of the structures and thickness of the flakes, the data concerning the geometry information (lattice and atomic structures) about the known 2D MXenes compounds have been taken into account. The theoretical descriptors have been build on the basis of the theoretical dataset, namely atomic types and positions using Atom-Centered Symmetry Functions (ACSF) [24]. The theoretical descriptors are collected in the Table A3.
Note, that the choice of the descriptors which represent the cytotoxic behavior is also limited to the experimental data available in literature. Therefore, the prediction might change whenever more experimental data are provided.

2.2. Datasets

The toxicological in vitro data for 2D MXenes is taken from recently published high-throughput screening experiments, therefore, it is reliable and convenient for comparison. We have decided to test independently three datatsets in order to answer following questions:
  • What kind of datasets and descriptors: theoretical, experimental or combined ones are sufficient to build a good quality ML model?
  • Could the cytotoxic behavior be predicted based on purely theoretical descriptors (type of atoms, stoichiometry, etc.), and hence, does it mean that no experimental data need to be provided to predict the cytotoxic behavior?
  • Does inclusion of the geometrical descriptors in the ML model improve qualitatively predictions of the cytotoxicity?
Three datatsets are listed below. For each of the datasets the ML models have been built separately and discussed in details in the corresponding subsections of Results.
  • Dataset I (experimental set)—the experimental data have been selected based on detailed description of the experiment, as well as detailed information about the structure, surface modification, in-depth characterization of morphology, and the cytotoxic effect of the MXenes compound on the cells tested in vitro. Those information have been collected from the literature listed in Table A1 and presented in Appendix A. It consists of 71 records and elemental features (descriptors) listed in Table A2.
  • Dataset II (theoretical set)—data taken from the two-dimensional database [25] concerning the geometry information about the known 2D MXenes compounds (61 records). The elemental features are collected in Table A3.
  • Dataset III (combined set)—dataset consists of both Dataset I and Dataset II. The number of the records and elemental features are combined from those two datasets correspondingly.
It is worth mentioning, that the first two datasets overlapped partially, namely, the geometry of each of the compounds in the first set is known. In addition, we do not know anything about the class of function describing cytotoxicity, so it is not possible to point out the class of algorithms that should be used. Thus, only the size and the type of the variables in the dataset determine our choice of ML algorithms. Detailed analysis of the applicability of machine learning algorithms can be found elsewhere [26]. Below, we present briefly the ML algorithms used in the present study:
  • Logistic regression [27] with regularization L1 and L2 (regLOG-L1, regLOG-L2). This approach allows avoiding over-learning a model even for a large number of variables. The algorithm removes unimportant features for the model.
  • Random Forest (RF) [22] is commonly used for a small dataset, and must be used with care regarding over-learning. It allows for selecting the most important features.
  • Support Vector Machine (SVM) [28] uses only part of the dataset, thus, it can be easily applied to a small size of dataset. The key point of prediction in the SVM algorithm is the choice of kernel. In this study, we have tested the commonly used kernels such as—linear, rbf, and sigmoid, denoted by us regSVM-lin, regSVM-rbf, regSVN-sig, respectively, throughout this paper.
  • Extreme Random Tree (ERT) [29] is an extension of a Random Forest algorithm, and is known to be computationally faster than RF. Both ERT and RF are known to work well for any dataset.
Parametric models such as linear regression are used to help us understand a phenomenon by determining the functional dependences. In the case of non-parametric models such as Random Forest, the crucial issue is to identify the importance of features, and thus, it allows us to understand the studied phenomenon. Note, that other commonly used ML methods such as Kernel Ridge Regression (KKR) [30] or Neural Networks (NN) [31], are well suited for large datasets, thus, are not applicable in our case.
In addition, our datasets face a commonly known issue, namely, the class imbalance problem. Significantly, this problem is widely reported for the toxicity of many other materials, where the size of the positive data (toxic samples) is considerably smaller than the negative data (non-toxic samples) (see Reference [32]). To solve this problem, we make use of various data-balancing algorithms such as—Weight classifier (Weight), and generating synthetic samples (SMOTE). The other commonly used algorithms such as oversample minority class or undersample majority class are not applicable in the case of MXenes materials, due to the small number of toxic records for which proper statistics cannot be built.
We have used the Python programming language (version 3.6.4) with the scikit-learn [33] and XRT [34] libraries for data analysis and machine learning. The Pandas [35] library was adopted in order to read and process the data, whereas the NumPy package (version 1.16.1) [36] was used to construct the features.

3. Results

In this study, we make use of state-of-the-art machine learning (ML) methods to identify the cytotoxicity of experimentally synthesized as well as deeply characterized 2D MXenes.
This section is divided into three subsections related to the prediction for the three datasets: experimental set, theoretical set and combined one, as it was mentioned in previous section. Then, the models are simplified, by selecting the most important features based on Random Forest algorithm or by the construction of new features from the given ones by the use of the Principle Component Analysis (PCA). The models are tested by 10-fold cross-validation, with the performance measured by class balanced accuracy score of correct predictions [37]. The accuracy score metric is defined in the range of [0, 1].

3.1. Dataset I—Experimental Set

Our theoretical investigations (see Table 1) reveal that the accuracy score for balanced data shows a good level of precision, greater than 0.72 (except for regSVM-sig) for all of the algorithms employed in this paper. Moreover, note that data balancing techniques improve the results approximately by a few percent (see the results for balanced versus unbalanced data collected in Table 1). The largest values are obtained in the case of the SVM with rbf kernel (regSVN-rbf) for Weight and SMOTE techniques, and are equal to 0.92 and 0.93 respectively. The high accuracy obtained in the case of regRF and regSVM-rbf manifest the non-linearity of the studied problem and the need for using non-parametric models. Unfortunately, this has a negative effect on understanding of the results and the underlying phenomenon. It is worth noting, that most of the variables used in this study are categorical variables described by one-hot-encoded or labeling methods. By use of the SMOTE data balancing algorithm, the results do not include pure features but contain mixed ones, which results in losing the physical interpretation of the outcomes. Thus, we have decided to use the Weight data balancing algorithms further in this study.
In order to understand the cytotoxicity issue and its dependency on selected descriptors, feature selection and feature engineering techniques are applied further in this study. The crucial premise of feature selection is that the data contains some variables that are either redundant or irrelevant, and can thus be removed without much loss of information. Feature selection and feature engineering techniques are methodological tools that allow simplification of the model, and hence, they can facilitate the interpretation of the studied phenomenon. In order to select the most important features, the feature importance score is obtained from a random forest analysis. Feature importance shows how much weight the model assigns to the given descriptor during predictions, and thus, gives insight into which variables are crucial for predicting the cytotoxicity of MXenes materials. For comparison purposes, we also use feature engineering from PCA.
The results reveal that three features are already sufficient for a good level of prediction accuracy, for all employed algorithms except regSVM-sig (see Figure 1).
Moreover, all the models based on PCA show lower accuracy than RF, which means that there is a low correlation between the features. The PCA approach is based on correlation between the features, while feature importance from RF selects the most important, original variables. PCA does not improve the results, thus, the original variables selected by RF have been chosen.
In addition, the feature importance score shows that there are two crucial features, four equally important features, and the rest seem to be unimportant from the RF analysis (see Figure 2). The most important features are the presence of the MxOy and the Li atoms on the MXenes surfaces.
Note that the results presented here show that experimental data and descriptors are sufficient to build a good quality ML model for cytotoxicity predictions in MXenes materials.

3.2. Dataset II—Experimental Set

Here, the ML model is built based on theoretical data. We have tested the dataset with structural information of the compounds included, namely position and type of atoms. There are many methods available for building descriptors such as Atom-Centered Symmetry Functions (ACSF) [24], Coulomb Matrix [38], or Ewald Sum Matrix [39], which convert the atomic positions into variables that can be used in machine learning. We have used the Weighted Atom-Centered Symmetry Functions (wACSF) [40] as descriptors, in order to substantially decrease the number of variables. The parameters of this model have been adopted from Reference [24]. Note that only 10% records of dataset have been labeled as toxic or non-toxic. In order to effectively elucidate the information contained in this dataset, clustering technique has been applied. Therefore, we have used unsupervised ML techniques. All these methods do not allow us to correctly distinguish toxic and non-toxic records in two separate clusters. In other words, the results reveal that the clustering technique cannot be viewed as a mechanism for toxicity prediction of MXenes.
All the results presented in this subsection reveal that taking into account geometrical descriptors from theoretical dataset of MXenes materials are not sufficient to build a ML model for cytotoxicity prediction of MXenes compounds. However, this theoretical set can be combined with experimental set, to determine whether such enlarged database improves model predictions, which is a subject of study in the next subsection.

3.3. Dataset III—Combined Experimental and Theoretical Sets

Our study reveals that including the information related to the geometry of the compounds (see Table 2) does not improve the results, and gives a similar level of accuracy score of correct predictions as obtained for dataset I (see Table 1), for all of the methods employed here.
Then, we have built models that include from one to five features by use of the Random Forest algorithm and PCA method, similarly to the approach for dataset I. The analysis demonstrates that two features are sufficient to describe the toxicity of MXenes compounds with high accuracy of predictions (see Figure 3). From the feature importance ranking, we find out that the topmost descriptors are the presence of MxOy, Li on the surface of MXenes, and surface modification with external compounds (see Figure 4). Note, that the order of the six top important features is the same as in the case for dataset I.

3.4. Discussion of Datasets

Here, we present the discussion concerning the results of employed datatsets in response to the following questions (previously stated in Section 2.2):
  • What kind of datasets and descriptors: theoretical, experimental or combined ones are sufficient to build a good quality ML model?
    Our results indicate that experimental data are sufficient to build an effective model with high accuracy score of correct cytotoxicity predictions based on 10-fold cross-validation scheme.
  • Could the cytotoxic behavior be predicted based on purely theoretical descriptors (type of atoms, stoichiometry etc.) and hence, does it mean that no experimental data need to be provided to predict the cytotoxic behavior?
    Our results reveal that taking into account only geometrical descriptors from theoretical dataset of MXenes materials are not sufficient to build a ML model for cytotoxicity predictions. Thus, the experimental information such as material synthesis methodology and characterisations play a crucial role in building ML model concerning cytotoxic behavior of MXenes compounds.
  • Does inclusion of the geometrical descriptors in the ML model improve qualitatively predictions of the cytotoxicity?
    The inclusion of geometrical features of MXenes do not qualitatively change the ML results based on experimental dataset.

4. Model Predictions

After successful verification of the ML model based on experimental dataset, we are able to predict toxicity for 2D MXenes not tested in vitro. To do so, we have searched the available literature covering the MXenes compounds for which all the elemental features listed in Table A2 have been provided, but with no in-vitro studies carried out. Despite the fact that there are around a hundred phases synthesized so far, we have only found 19 MXenes compounds, for which comprehensive data on the material are available (see Table 3).
The ML models predict two of 2D MXenes can exhibit cytotoxic properties with a high probability of prediction equal to 0.9, while the rest of them are predicted to be non-toxic (see Table 3). It is worth mentioning that for the non-toxic ones, no presence of MxOy on the surface has been reported. The presence of MxOy is the key toxicity-generating feature obtained from our studies.
Our results demonstrate that our ML model is able to complement the existing knowledge coming from in vitro studies. However, note that this prediction has to be viewed with some caution, knowing that traditional k-fold cross-validation is highly optimistic when evaluating machine learning models, due to the fact that materials datasets are rarely uniformly distributed.

5. Discussion

Our results indicate that knowledge of the surface and its modification might be crucial issue concerning the toxicity of the studied layered 2D materials, whereas geometrical descriptors may have little impact on the outcomes. These results are in the line with recent experimental findings concerning the presence of TixOy on the surface [8,58,59], and the biological knowledge of cytotoxicity mechanisms [60], as well as physical and chemical intuition. It should be stressed that this conclusion is much more definitive than we expected at the beginning of our studies. The reason is that the chemical diversity and inhomogeneity of MXenes are already widely known and pose a major challenge in such complex analysis. The second corresponding aspect is the divided surface characteristics. The primary strictly depends on the starting materials (MAX phases). The latter is undoubtedly far more problematic, if it comes into interactions with highly sensitive systems such as biological ones. Basically, the chemical composition of the surface of MXenes, almost certainly, is closely related to the type of ‘M’ element and a resulting chemical composition of the MxOy passivation layer that occurs as a result of M reaction with oxygen and/or water [61]. In fact, every surface of the MXene exposed to the air can naturally react with oxygen because the freshly exposed metallic surface is energetically unsaturated and possesses high reactivity. This can happen immediately after the delamination process (but certainly must also depend on the MXene stability). What’s more, the freshly exposed surface of the MXene also acquires bonding with products of chemical reactions that occur during acidic etching of the ‘A’ element from the MAX phase. As can be seen, the aforementioned surface-related features influence at first the material itself, but finally, may result in different biological effects such as the appearance or lack of cytotoxicity.
From this regard, MXenes are the most interesting, because their surface is highly unstable and susceptible to oxidation in the ambient conditions [52]. It has been shown that the surface effects highly affect the living cells and it is closely linked to the oxidation process and decomposition to toxic oxides. In such a case, mechanisms of cytotoxic action refer to cell cycle, DNA synthesis, and cellular membrane integrity [62]. From this perspective, surface chemistry of MXenes should be managed to avoid further toxic effects using redirecting into more safe surface compositions [59]. This is indeed a right direction when the surface modification cannot be used. In other cases, polyanionic salts can be effectively used for MXene flakes edges capping which results in highly decreased decomposition and stability [50]. In other approach, natural antioxidants are good in diminishing surface oxidation [49]. In addition, approaches that involve surface modification with bio-organic moieties, specifically designed for cytotoxicity mitigation, are also highly recommended [62,63].

6. Conclusions

Here, we present the first theoretical study concerning the toxicological aspects of 2D MXene materials by employing various machine learning models. Our work demonstrate that the most important features potentially responsible for the toxicological properties are related to the presence of transition metal oxides MxOy and Lithium atoms on the surface, as well as surface modification with external compounds. Our detailed analysis reveals, that the crucial issue is what happens on the surface, while the structural information of the systems might have minimal impact on cytotoxicological aspects of MXenes materials.
Our ML model successfully complement existing experimental studies, for which no cytotoxicological measurements have been carried out. In particular, we have predicted the cytotoxicity of 19 MXenes compounds, for which two of them are predicted to be cytotoxic with 0.9 probability. The rest of the compounds are predicted to be non-toxic and can be potentially applied in many technological areas [2,64,65,66,67,68].
Moreover, our results show that the cytotoxic prediction of MXenes can be examined on the materials that are well experimentally characterised in terms of surface chemistry and the presence of oxides on the MXenes surface. Thus, we claim that it might be one of the solutions for reducing the number of toxicological studies needed, and allows for minimizing failures in future biological applications.
In addition, the theoretical research methodology based on ML models developed here can be further applied to other types of 2D materials exhibiting complex structure and diverse surface characteristics, such as for example, novel 2D transition metal borides, so called MBenes [69] as well as van der Waals heterostructures [70,71]. We expect that the predictions presented here will facilitate the experimental efforts by providing the information that might accelerate time consuming and expensive cytotoxical experimental studies, by reducing the large number of compounds, and hence, speeding up potential future applications.

Author Contributions

Conceptualization, A.M.J.; methodology, M.E.M.; software, M.E.M.; validation, M.B. and A.M.J.; formal analysis, M.E.M. and M.B.; investigation, M.E.M. and M.P.; resources, M.P. and A.M.J.; data curation, M.P. and M.E.M.; writing—original draft preparation, M.B.; writing–review and editing, M.B., A.M.J. and J.A.M.; visualization, M.B. and M.P.; supervision, M.B. and J.A.M.; project administration, A.M.J.; funding acquisition, A.M.J. All authors have read and agreed to the published version of the manuscript.

Funding

The study was accomplished thanks to the funds allotted by the National Science Centre on the basis of decision no. DEC-2017/26/E/ST8/01073, within the framework of the research project ‘SONATA BIS 7’ no. UMO-2017/26/E/ST8/01073. M.B. is funded by the National Science Centre grant no. UMO-2016/23/D/ST3/03446. J.A.M acknowledges the support of NCN through the grant OPUS-16 (UMO-2018/31/B/ST3/03758). Access to computing facilities of PL-Grid Polish Infrastructure for Supporting Computational Science in the European Research Space and of the Interdisciplinary Center of Modeling (ICM), University of Warsaw is gratefully acknowledged. We made use of computing facilities of TU Dresden ZIH within the project TransPheMat.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Abbreviations

The following abbreviations are used in this manuscript:
SPSoybean phospholipid
PVPpolyvinylpyrrolidone
HAhyaluronic acid
PEIpolyethylene imine
PEGpolyethylene glycol
PVApolyvinyl alcohol
NMPN-methyl-2-pyrrolidone
TMAOHtetramethylammonium hydroxide
TBAOHtetrabuthylammonium hydroxide
TPAOHtetrapropylammonium hydroxide
DMFdimethylformamide
DMSOdimethyl sulfoxide
CTACcetanecyltrimethylammonium chloride
c(RGDyC)cyclic arginine-glycine-aspartic pentapeptide
PLLpoly-L-lysine
APTES(3-aminopropyl) triethoxysilane
CEAcarcinoembryonic antigen
L-ACCL-ascorbic acid
PASpolyanionic salts

Appendix A. Detailed Information about the Experimental Data Used in the Machine Learning Models

Here we present the list of the high-throughput experimental data that we used to build the database (see Table A1).
Table
Table A1. The types of delaminated 2D MXenes compounds used in this study. The information about the particular compound are presented in details in the following references.
Table A1. The types of delaminated 2D MXenes compounds used in this study. The information about the particular compound are presented in details in the following references.
MXenes Compound:
Ti3C2References [8,41,42,43,44,45,46,47,48,49,50,51,52,72,73,74,75,76,77,78,79,80,81]
Ta4C3References [82,83,84]
Nb2CReferences [54,55,85,86]
Ti2CReference [67]
Mo2CReferences [87]
Mo1.33CReference [56]
Nb4C3References [62]
V2CReferences [50,53]
Ti2NReference [10]
Ti4N3Reference [57]

Appendix B. Detailed Information about the Elemental Features Used in the Machine Learning Models

Here we present the detailed information about the elemental features taken into account in the machine learning predictions, listed in Table A2 and Table A3 respectively.
Table
Table A2. Detailed description of the elemental features used in ML scheme and applied for the dataset I and dataset III.
Table A2. Detailed description of the elemental features used in ML scheme and applied for the dataset I and dataset III.
Detailed Description of the Feature
Surface modificationPVP, SP, MnOx+SP, HA, FexOy+SP, PEI, PEG, DNA+Pt+Pd, Ag, L-ACC, PAS
with external compoundsCTAC+PEG+SiO2+c(RGDyC)+SiO2, Au+Fe3O4, Au+PEG, PLL, APTES+CEA, PVA, Au
Lateral sizefrom few to hundredths of nm, from few to hundredths of μ m
Thicknessfrom few to tens of nm
Etching agentHF, LiF+HCl, LiF+HCl+AlCl3, NaF+HCl, KF+HCl, KF+LiF+NaF,
Delaminating agentUltrasounds, DMF+high pressure+high Temp., high pressure+high Temp., TBAOH, TBAOH+ultrasounds, TPAOH, TMAOH, DMSO+ultrasounds, no additional treatment
Carbon (C) on a surface1-Yes, 0-No
Oxygen (O) on a surface1-Yes, 0-No
Fluor (F) on a surface1-Yes, 0-No
Aluminium (Al) on a surface1-Yes, 0-No
Titanium (Ti) on a surface1-Yes, 0-No
Nitrogen (N) on a surface1-Yes, 0-No
Chloride (Cl) on a surface1-Yes, 0-No
Silicon (Si) on a surface1-Yes, 0-No
Lithium (Li) on a surface1-Yes, 0-No
Other on a surface1-Yes, 0-No
presence of MxOy1-Yes, 0-No
Toxic1-Yes, 0-No
Table
Table A3. Detailed description of the elemental features used in ML scheme for the dataset II and constituting the part of the dataset III.
Table A3. Detailed description of the elemental features used in ML scheme for the dataset II and constituting the part of the dataset III.
G 2 ∖ w G 2 see Reference [24]; cutoff RC = 6, shift Rs = ( 0 ⋯ N ) , number of the functions N = 10, η = ( 1 , N R s N )
G 4 ∖ w G 4 see Reference [24]; cutoff RC = 6, shift Rs = ( 0 ⋯ N ) , number of the functions N = 10, λ = − 1 , 1 , ζ = 1 , 2 , 4 , 16

References

  1. Zhang, J.F.; Cao, H.Y.; Wang, H.B. Research Progress of Novel Two-dimensional Material MXene. J. Inorg. Mater. 2017, 32, 561. [Google Scholar] [CrossRef] [Green Version]
  2. Hong Ng, V.M.; Huang, H.; Zhou, K.; Lee, P.S.; Que, W.; Xu, J.Z.; Kong, L.B. Recent progress in layered transition metal carbides and/or nitrides (MXenes) and their composites: Synthesis and applications. J. Mater. Chem. A 2017, 5, 3039–3068. [Google Scholar] [CrossRef]
  3. Verger, L.; Xu, C.; Natu, V.; Cheng, H.M.; Ren, W.; Barsoum, M.W. Overview of the synthesis of MXenes and other ultrathin 2D transition metal carbides and nitrides. Curr. Opin. Solid State Mater. Sci. 2019, 23, 149–163. [Google Scholar] [CrossRef]
  4. Naguib, M.; Mashtalir, O.; Carle, J.; Presser, V.; Lu, J.; Hultman, L.; Gogotsi, Y.; Barsoum, M.W. Two-Dimensional Transition Metal Carbides. ACS Nano 2012, 6, 1322–1331. [Google Scholar] [CrossRef]
  5. Naguib, M.; Halim, J.; Lu, J.; Cook, K.M.; Hultman, L.; Gogotsi, Y.; Barsoum, M.W. New Two-Dimensional Niobium and Vanadium Carbides as Promising Materials for Li-Ion Batteries. J. Am. Chem. Soc. 2013, 135, 15966–15969. [Google Scholar] [CrossRef]
  6. Khazaei, M.; Arai, M.; Sasaki, T.; Chung, C.Y.; Venkataramanan, N.S.; Estili, M.; Sakka, Y.; Kawazoe, Y. Novel Electronic and Magnetic Properties of Two-Dimensional Transition Metal Carbides and Nitrides. Adv. Funct. Mater. 2013, 23, 2185–2192. [Google Scholar] [CrossRef]
  7. Anasori, B.; Lukatskaya, M.R.; Gogotsi, Y. 2D metal carbides and nitrides (MXenes) for energy storage. Nat. Rev. Mater. 2017, 2. [Google Scholar] [CrossRef]
  8. Jastrzębska, A.; Szuplewska, A.; Wojciechowski, T.; Chudy, M.; Ziemkowska, W.; Chlubny, L.; Rozmysłowska, A.; Olszyna, A. In vitro studies on cytotoxicity of delaminated Ti3C2 MXene. J. Hazard. Mater. 2017, 339, 1–8. [Google Scholar] [CrossRef]
  9. Jastrzębska, A.M.; Karwowska, E.; Wojciechowski, T.; Ziemkowska, W.; Rozmysłowska, A.; Chlubny, L.; Olszyna, A.R. The Atomic Structure of Ti2C and Ti3C2 MXenes is Responsible for Their Antibacterial Activity Toward E. coli Bacteria. J. Mater. Eng. Perform. 2018, 28, 1272–1277. [Google Scholar] [CrossRef]
  10. Soundiraraju, B.; George, B.K. Two-Dimensional Titanium Nitride (Ti2N) MXene: Synthesis, Characterization, and Potential Application as Surface-Enhanced Raman Scattering Substrate. ACS Nano 2017, 11, 8892–8900. [Google Scholar] [CrossRef]
  11. Scheibe, B.; Wychowaniec, J.K.; Scheibe, M.; Peplińska, B.; Jarek, M.; Nowaczyk, G.; Przysiecka, Ł. Cytotoxicity Assessment of Ti-Al-C Based MAX Phases and Ti3C2Tx MXenes on Human Fibroblasts and Cervical Cancer Cells. ACS Biomater. Sci. Eng. 2019, 5, 6557–6569. [Google Scholar] [CrossRef]
  12. Yin, Z.; Ai, H.; Zhang, L.; Ren, G.; Wang, Y.; Zhao, Q.; Liu, H. Predicting the cytotoxicity of chemicals using ensemble learning methods and molecular fingerprints. J. Appl. Toxicol. 2019, 39, 1366–1377. [Google Scholar] [CrossRef] [PubMed]
  13. Frey, N.C.; Wang, J.; Vega Bellido, G.I.; Anasori, B.; Gogotsi, Y.; Shenoy, V.B. Prediction of Synthesis of 2D Metal Carbides and Nitrides (MXenes) and Their Precursors with Positive and Unlabeled Machine Learning. ACS Nano 2019, 13, 3031–3041. [Google Scholar] [CrossRef] [PubMed]
  14. Madiona, R.M.; Winkler, D.A.; Muir, B.W.; Pigram, P.J. Optimal machine learning models for robust materials classification using ToF-SIMS data. Appl. Surf. Sci. 2019, 487, 773–783. [Google Scholar] [CrossRef]
  15. Breiman, L. Arcing The Edge; Technical Report 486; University of California: Berkeley, CA, USA, 1997. [Google Scholar]
  16. Wu, Y.; Wang, G. Machine Learning Based Toxicity Prediction: From Chemical Structural Description to Transcriptome Analysis. Int. J. Mol. Sci. 2018, 19, 2358. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Grenet, I.; Yin, Y.; Comet, J.P.; Gelenbe, E. Machine Learning to Predict Toxicity of Compounds. In Proceedings of the ICANN 2018: Artificial Neural Networks and Machine Learning, Rhodes, Greece, 4–7 October 2018; Springer: Berlin, Germany, 2018; pp. 335–345. [Google Scholar]
  18. Mayr, A.; Klambauer, G.; Unterthiner, T.; Hochreiter, S. DeepTox: Toxicity Prediction using Deep Learning. Front. Environ. Sci. 2016, 3, 80. [Google Scholar] [CrossRef] [Green Version]
  19. Lei, T.; Sun, H.; Kang, Y.; Zhu, F.; Liu, H.; Zhou, W.; Wang, Z.; Li, D.; Li, Y.; Hou, T. ADMET Evaluation in Drug Discovery. 18. Reliable Prediction of Chemical-Induced Urinary Tract Toxicity by Boosting Machine Learning Approaches. Mol. Pharm. 2017, 14, 3935–3953. [Google Scholar] [CrossRef]
  20. Grenet, I.; Yin, Y.; Comet, J.P. G-Networks to Predict the Outcome of Sensing of Toxicity. Sensors 2018, 18, 3483. [Google Scholar] [CrossRef] [Green Version]
  21. Grenet, I.; Merlo, K.; Comet, J.P.; Tertiaux, R.; Rouquié, D.; Dayan, F. Stacked Generalization with Applicability Domain Outperforms Simple QSAR on in Vitro Toxicological Data. J. Chem. Inf. Model. 2019, 59, 1486–1496. [Google Scholar] [CrossRef]
  22. Ho, T.K. Random Decision Forest. In Proceedings of the 3rd International Conference on Document Analysis and Recognition, Montreal, QC, Canada, 14–16 August 1995. [Google Scholar]
  23. Jolliffe, I. Principal Component Analysis, Second Edition; Springer: Berlin, Germany, 1986. [Google Scholar]
  24. Behler, J. Atom-centered symmetry functions for constructing high-dimensional neural network potentials. J. Chem. Phys. 2011, 134, 074106. [Google Scholar] [CrossRef]
  25. Jain, A.; Ong, S.P.; Hautier, G.; Chen, W.; Richards, W.D.; Dacek, S.; Cholia, S.; Gunter, D.; Skinner, D.; Ceder, G.; et al. The Materials Project: A materials genome approach to accelerating materials innovation. APL Mater. 2013, 1, 011002. [Google Scholar] [CrossRef] [Green Version]
  26. Forman, G.; Cohen, I. Learning from Little: Comparison of Classifiers Given Little Training. In Knowledge Discovery in Databases: PKDD 2004; Boulicaut, J.F., Esposito, F., Giannotti, F., Pedreschi, D., Eds.; Springer: Berlin/Heidelberg, Germany, 2004; pp. 161–172. [Google Scholar]
  27. Bishop, C.M. Pattern Recognition and Machine Learning; Springer: Berlin, Germany, 2006. [Google Scholar]
  28. Boser, B.E.; Guyon, I.M.; Vapnik, V.N. A Training Algorithm for Optimal Margin Classifiers. In Proceedings of the Fifth Annual Workshop on Computational Learning Theory, COLT ’92, Pittsburgh, PA, USA, 27–29 July 1992; ACM: New York, NY, USA, 1992; pp. 144–152. [Google Scholar] [CrossRef]
  29. Geurts, P.; Ernst, D.; Wehenkel, L. Extremely randomized trees. Mach. Learn. 2006, 63, 3–42. [Google Scholar] [CrossRef] [Green Version]
  30. Zhang, Y.; Duchi, J.; Wainwright, M. Divide and Conquer Kernel Ridge Regression. In Machine Learning Research, Proceedings of the 26th Annual Conference on Learning Theory, Princeton, NJ, USA, 12–14 June 2013; Shalev-Shwartz, S., Steinwart, I., Eds.; PMLR: Princeton, NJ, USA, 2013; Volume 30, pp. 592–617. [Google Scholar]
  31. McCulloch, W.S.; Pitts, W. A logical calculus of the ideas immanent in nervous activity. Bull. Math. Biophys. 1943, 5, 115–133. [Google Scholar] [CrossRef]
  32. Sun, Y.; Wong, A.K.C.; Kamel, M.S. Classification of imbalanced data: A review. Int. J. Pattern Recognit. Artif. Intell. 2009, 23, 687–719. [Google Scholar] [CrossRef]
  33. Pedregosa, F.; Varoquaux, G.; Gramfort, A.; Michel, V.; Thirion, B.; Grisel, O.; Blondel, M.; Prettenhofer, P.; Weiss, R.; Dubourg, V.; et al. Scikit-learn: Machine Learning in Python. J. Mach. Learn. Res. 2011, 12, 2825–2830. [Google Scholar]
  34. Klementiev, K.; Chernikov, R. Powerful scriptable ray tracing package xrt. In Advances in Computational Methods for X-Ray Optics III; Del Rio, M.S., Chubar, O., Eds.; International Society for Optics and Photonics, SPIE: Bellingham, WA, USA, 2014; Volume 9209, pp. 60–75. [Google Scholar] [CrossRef]
  35. Available online: https://pandas.pydata.org/ (accessed on 9 July 2020).
  36. Available online: https://numpy.org/ (accessed on 9 July 2020).
  37. Mosley, L. QA Balanced Approach to the Multi-Class Imbalance Problem. Master’s Thesis, Iowa State University, Ames, IA, USA, 2013. [Google Scholar]
  38. Rupp, M.; Tkatchenko, A.; Müller, K.R.; Von Lilienfeld, O.A. Fast and Accurate Modeling of Molecular Atomization Energies with Machine Learning. Phys. Rev. Lett. 2012, 108, 058301. [Google Scholar] [CrossRef]
  39. Faber, F.; Lindmaa, A.; Von Lilienfeld, O.A.; Armiento, R. Crystal structure representations for machine learning models of formation energies. Int. J. Quantum Chem. 2015, 115, 1094–1101. [Google Scholar] [CrossRef] [Green Version]
  40. Gastegger, M.; Schwiedrzik, L.; Bittermann, M.; Berzsenyi, F.; Marquetand, P. wACSF—Weighted atom-centered symmetry functions as descriptors in machine learning potentials. J. Chem. Phys. 2018, 148, 241709. [Google Scholar] [CrossRef] [Green Version]
  41. Ren, C.E.; Hatzell, K.B.; Alhabeb, M.; Ling, Z.; Mahmoud, K.A.; Gogotsi, Y. Charge- and Size-Selective Ion Sieving Through Ti3C2Tx MXene Membranes. J. Phys. Chem. Lett. 2015, 6, 4026–4031. [Google Scholar] [CrossRef]
  42. Rakhi, R.B.; Nayak, P.; Xia, C.; Alshareef, H.N. Novel amperometric glucose biosensor based on MXene nanocomposite. Sci. Rep. 2016, 6, 36422. [Google Scholar] [CrossRef] [Green Version]
  43. Shahzad, A.; Rasool, K.; Miran, W.; Nawaz, M.; Jang, J.; Mahmoud, K.A.; Lee, D.S. Two-Dimensional Ti3C2Tx MXene Nanosheets for Efficient Copper Removal from Water. ACS Sustain. Chem. Eng. 2017, 5, 11481–11488. [Google Scholar] [CrossRef]
  44. Wu, L.; Lu, X.; Dhanjai; Wu, Z.S.; Dong, Y.; Wang, X.; Zheng, S.; Chen, J. 2D transition metal carbide MXene as a robust biosensing platform for enzyme immobilization and ultrasensitive detection of phenol. Biosens. Bioelectron. 2018, 107, 69–75. [Google Scholar] [CrossRef] [PubMed]
  45. Arabi Shamsabadi, A.; Sharifian Gh., M.; Anasori, B.; Soroush, M. Antimicrobial Mode-of-Action of Colloidal Ti3C2Tx MXene Nanosheets. ACS Sustain. Chem. Eng. 2018, 6, 16586–16596. [Google Scholar] [CrossRef]
  46. Kumar, S.; Lei, Y.; Alshareef, N.H.; Quevedo-Lopez, M.; Salama, K.N. Biofunctionalized two-dimensional Ti3C2 MXenes for ultrasensitive detection of cancer biomarker. Biosens. Bioelectron. 2018, 121, 243–249. [Google Scholar] [CrossRef] [Green Version]
  47. Zheng, J.; Wang, B.; Ding, A.; Weng, B.; Chen, J. Synthesis of MXene/DNA/Pd/Pt nanocomposite for sensitive detection of dopamine. J. Electroanal. Chem. 2018, 816, 189–194. [Google Scholar] [CrossRef]
  48. Pandey, R.P.; Rasool, K.; Madhavan, V.E.; Aïssa, B.; Gogotsi, Y.; Mahmoud, K.A. Ultrahigh-flux and fouling-resistant membranes based on layered silver/MXene (Ti3C2Tx) nanosheets. J. Mater. Chem. A 2018, 6, 3522–3533. [Google Scholar] [CrossRef]
  49. Zhao, X.; Vashisth, A.; Prehn, E.; Sun, W.; Shah, S.A.; Habib, T.; Chen, Y.; Tan, Z.; Lutkenhaus, J.L.; Radovic, M.; et al. Antioxidants Unlock Shelf-Stable Ti3C2Tx (MXene) Nanosheet Dispersions. Matter 2019, 1, 513–526. [Google Scholar] [CrossRef] [Green Version]
  50. Natu, V.; Hart, J.L.; Sokol, M.; Chiang, H.; Taheri, M.L.; Barsoum, M.W. Edge Capping of 2D-MXene Sheets with Polyanionic Salts To Mitigate Oxidation in Aqueous Colloidal Suspensions. Angew. Chem. Int. Ed. 2019, 58, 12655–12660. [Google Scholar] [CrossRef] [PubMed]
  51. Zhu, J.; Tang, Y.; Yang, C.; Wang, F.; Cao, M. Composites of TiO2 Nanoparticles Deposited on Ti3C2 MXene Nanosheets with Enhanced Electrochemical Performance. J. Electrochem. Soc. 2016, 163, A785–A791. [Google Scholar] [CrossRef]
  52. Zhang, C.J.; Pinilla, S.; McEvoy, N.; Cullen, C.P.; Anasori, B.; Long, E.; Park, S.H.; Seral-Ascaso, A.; Shmeliov, A.; Krishnan, D.; et al. Oxidation Stability of Colloidal Two-Dimensional Titanium Carbides (MXenes). Chem. Mater. 2017, 29, 4848–4856. [Google Scholar] [CrossRef]
  53. Liu, F.; Zhou, J.; Wang, S.; Wang, B.; Shen, C.; Wang, L.; Hu, Q.; Huang, Q.; Zhou, A. Preparation of High-Purity V2C MXene and Electrochemical Properties as Li-Ion Batteries. J. Electrochem. Soc. 2017, 164, 709–713. [Google Scholar] [CrossRef]
  54. Mashtalir, O.; Lukatskaya, M.R.; Zhao, M.Q.; Barsoum, M.W.; Gogotsi, Y. Amine-Assisted Delamination of Nb2C MXene for Li-Ion Energy Storage Devices. Adv. Mater. 2015, 27, 3501–3506. [Google Scholar] [CrossRef] [PubMed]
  55. Peng, C.; Wei, P.; Chen, X.; Zhang, Y.; Zhu, F.; Cao, Y.; Wang, H.; Yu, H.; Peng, F. A hydrothermal etching route to synthesis of 2D MXene (Ti3C2, Nb2C): Enhanced exfoliation and improved adsorption performance. Ceram. Int. 2018, 44, 18886–18893. [Google Scholar] [CrossRef]
  56. Tao, Q.; Dahlqvist, M.; Lu, J.; Kota, S.; Meshkian, R.; Halim, J.; Palisaitis, J.; Hultman, L.; Barsoum, M.W.; Persson, P.O.; et al. Two-dimensional Mo1.33C MXene with divacancy ordering prepared from parent 3D laminate with in-plane chemical ordering. Nat. Commun. 2017, 8. [Google Scholar] [CrossRef]
  57. Urbankowski, P.; Anasori, B.; Makaryan, T.; Er, D.; Kota, S.; Walsh, P.L.; Zhao, M.; Shenoy, V.B.; Barsoum, M.W.; Gogotsi, Y. Synthesis of two-dimensional titanium nitride Ti4N3 (MXene). Nanoscale 2016, 8, 11385–11391. [Google Scholar] [CrossRef]
  58. Szuplewska, A.; Rozmysłowska-Wojciechowska, A.; Poźniak, S.; Wojciechowski, T.; Birowska, M.; Popielski, M.; Chudy, M.; Ziemkowska, W.; Chlubny, L.; Moszczyńska, D.; et al. Multilayered stable 2D nano-sheets of Ti2NTx MXene: Synthesis, characterization, and anticancer activity. J. Nanobiotechnol. 2019, 17. [Google Scholar] [CrossRef]
  59. Jastrzębska, A.M.; Szuplewska, A.; Rozmysłowska-Wojciechowska, A.; Chudy, M.; Olszyna, A.; Birowska, M.; Popielski, M.; Majewski, J.A.; Scheibe, B.; Natu, V.; et al. On tuning the cytotoxicity of Ti3C2 (MXene) flakes to cancerous and benign cells by post-delamination surface modifications. 2D Mater. 2020, 7, 025018. [Google Scholar] [CrossRef]
  60. Franqui, L.; De Luna, L.; Loret, T.; Martinez, D.; Bussy, C. Assessing the Adverse Effects of Two-Dimensional Materials Using Cell Culture-Based Models.; Springer: Berlin, Germany, 2019. [Google Scholar] [CrossRef]
  61. Habib, T.; Zhao, X.; Shah, S.A.; Chen, Y.; Sun, W.; An, H.; Lutkenhaus, J.L.; Radovic, M.; Green, M.J. Oxidation stability of Ti3C2Tx MXene nanosheets in solvents and composite films. NPJ 2D Mater. Appl. 2019, 3. [Google Scholar] [CrossRef]
  62. Jastrzebska, A.; Szuplewska, A.; Rozmysłowska-Wojciechowska, A.; Mitrzak, J.; Wojciechowski, T.; Chudy, M.; Moszczyńska, D.; Wójcik, A.; Prenger, K.; Naguib, M. Juggling Surface Charges of 2D Niobium Carbide MXenes for a Reactive Oxygen Species Scavenging and Effective Targeting of the Malignant Melanoma Cell Cycle into Programmed Cell Death. ACS Sustain. Chem. Eng. 2020, 8, 7942–7951. [Google Scholar] [CrossRef]
  63. Rozmysłowska-Wojciechowska, A.; Szuplewska, A.; Wojciechowski, T.; Poźniak, S.; Mitrzak, J.; Chudy, M.; Ziemkowska, W.; Chlubny, L.; Olszyna, A.; Jastrzębska, A. A simple, low-cost and green method for controlling the cytotoxicity of MXenes. Mater. Sci. Eng. C 2020, 111, 110790. [Google Scholar] [CrossRef]
  64. Li, X.; Wang, C.; Cao, Y.; Wang, G. Functional MXene Materials: Progress of Their Applications. Chem. Asian J. 2018, 13, 2742–2757. [Google Scholar] [CrossRef] [PubMed]
  65. Liu, F.; Zhou, A.; Chen, J.; Jia, J.; Zhou, W.; Wang, L.; Hu, Q. Preparation of Ti3C2 and Ti2C MXenes by fluoride salts etching and methane adsorptive properties. Appl. Surf. Sci. 2017, 416, 781–789. [Google Scholar] [CrossRef]
  66. Zhu, J.; Ha, E.; Zhao, G.; Zhou, Y.; Huang, D.; Yue, G.; Hu, L.; Sun, N.; Wang, Y.; Lee, L.Y.S.; et al. Functional MXene Materials: Progress of Their Applications. Coord. Chem. Rev. 2017, 352, 306–327. [Google Scholar] [CrossRef]
  67. Szuplewska, A.; Kulpińska, D.; Dybko, A.; Jastrzębska, A.M.; Wojciechowski, T.; Rozmysłowska, A.; Chudy, M.; Grabowska-Jadach, I.; Ziemkowska, W.; Brzózka, Z.; et al. 2D Ti2C (MXene) as a novel highly efficient and selective agent for photothermal therapy. Mater. Sci. Eng. C 2019, 98, 874–886. [Google Scholar] [CrossRef] [PubMed]
  68. Rozmysłowska-Wojciechowska, A.; Wojciechowski, T.; Ziemkowska, W.; Chlubny, L.; Olszyna, A.; Jastrzębska, A. Surface interactions between 2D Ti3C2/Ti2C MXenes and lysozyme. Appl. Surf. Sci. 2019, 473, 409–418. [Google Scholar] [CrossRef]
  69. Khazaei, M.; Wang, J.; Estili, M.; Ranjbar, A.; Suehara, S.; Arai, M.; Esfarjani, K.; Yunoki, S. Novel MAB phases and insights into their exfoliation into 2D MBenes. Nanoscale 2019, 11, 11305–11314. [Google Scholar] [CrossRef]
  70. Geim, A.K.; Grigorieva, I.V. Van der Waals heterostructures. Nature 2013, 499, 419–425. [Google Scholar] [CrossRef]
  71. Birowska, M.; Urban, J.; Baranowski, M.; Maude, D.K.; Plochocka, P.; Szwacki, N.G. The impact of hexagonal boron nitride encapsulation on the structural and vibrational properties of few layer black phosphorus. Nanotechnology 2019, 30, 195201. [Google Scholar] [CrossRef] [Green Version]
  72. Xue, Q.; Zhang, H.; Zhu, M.; Pei, Z.; Li, H.; Wang, Z.; Huang, Y.; Huang, Y.; Deng, Q.; Zhou, J.; et al. Photoluminescent Ti3C2 MXene Quantum Dots for Multicolor Cellular Imaging. Adv. Mater. 2017, 29, 1604847. [Google Scholar] [CrossRef]
  73. Yu, X.; Cai, X.; Cui, H.; Lee, S.W.; Yu, X.F.; Liu, B. Fluorine-free preparation of titanium carbide MXene quantum dots with high near-infrared photothermal performances for cancer therapy. Nanoscale 2017, 9, 17859–17864. [Google Scholar] [CrossRef]
  74. Zhou, L.; Wu, F.; Yu, J.; Deng, Q.; Zhang, F.; Wang, G. Titanium carbide (Ti3C2Tx) MXene: A novel precursor to amphiphilic carbide-derived graphene quantum dots for fluorescent ink, light-emitting composite and bioimaging. Carbon 2017, 118, 50–57. [Google Scholar] [CrossRef]
  75. Lin, H.; Wang, X.; Yu, L.; Chen, Y.; Shi, J. Two-Dimensional Ultrathin MXene Ceramic Nanosheets for Photothermal Conversion. Nano Lett. 2017, 17, 384–391. [Google Scholar] [CrossRef] [PubMed]
  76. Dai, C.; Lin, H.; Xu, G.; Liu, Z.; Wu, R.; Chen, Y. Biocompatible 2D Titanium Carbide (MXenes) Composite Nanosheets for pH-Responsive MRI-Guided Tumor Hyperthermia. Chem. Mater. 2017, 29, 8637–8652. [Google Scholar] [CrossRef]
  77. Liu, G.; Zou, J.; Tang, Q.; Yang, X.; Zhang, Y.; Zhang, Q.; Huang, W.; Chen, P.; Shao, J.; Dong, X. Surface Modified Ti3C2 MXene Nanosheets for Tumor Targeting Photothermal/Photodynamic/Chemo Synergistic Therapy. ACS Appl. Mater. Interfaces 2017, 9, 40077–40086. [Google Scholar] [CrossRef]
  78. Chen, X.; Sun, X.; Xu, W.; Pan, G.; Zhou, D.; Zhu, J.; Wang, H.; Bai, X.; Dong, B.; Song, H. Ratiometric photoluminescence sensing based on Ti3C2 MXene quantum dots as an intracellular pH sensor. Nanoscale 2018, 10, 1111–1118. [Google Scholar] [CrossRef] [PubMed]
  79. Han, X.; Huang, J.; Lin, H.; Wang, Z.; Li, P.; Chen, Y. 2D Ultrathin MXene-Based Drug-Delivery Nanoplatform for Synergistic Photothermal Ablation and Chemotherapy of Cancer. Adv. Healthc. Mater. 2018, 7, 1701394. [Google Scholar] [CrossRef]
  80. Hussein, E.A.; Zagho, M.M.; Rizeq, B.R.; Younes, N.N.; Pintus, G.; Mahmoud, K.A.; Nasrallah, G.K.; Elzatahry, A.A. Plasmonic MXene-based nanocomposites exhibiting photothermal therapeutic effects with lower acute toxicity than pure MXene. Int. J. Nanomed. 2019, 14, 4529–4539. [Google Scholar] [CrossRef] [Green Version]
  81. Tang, W.; Dong, Z.; Zhang, R.; Yi, X.; Yang, K.; Jin, M.; Yuan, C.; Xiao, Z.; Liu, Z.; Cheng, L. Multifunctional Two-Dimensional Core–Shell MXene@Gold Nanocomposites for Enhanced Photo–Radio Combined Therapy in the Second Biological Window. ACS Nano 2019, 13, 284–294. [Google Scholar] [CrossRef]
  82. Dai, C.; Chen, Y.; Jing, X.; Xiang, L.; Yang, D.; Lin, H.; Liu, Z.; Han, X.; Wu, R. Two-Dimensional Tantalum Carbide (MXenes) Composite Nanosheets for Multiple Imaging-Guided Photothermal Tumor Ablation. ACS Nano 2017, 11, 12696–12712. [Google Scholar] [CrossRef]
  83. Lin, H.; Wang, Y.; Gao, S.; Chen, Y.; Shi, J. Theranostic 2D Tantalum Carbide (MXene). Adv. Mater. 2018, 30, 1703284. [Google Scholar] [CrossRef]
  84. Liu, Z.; Lin, H.; Zhao, M.; Dai, C.; Zhang, S.; Peng, W.; Chen, Y. 2D Superparamagnetic Tantalum Carbide Composite MXenes for Efficient Breast-Cancer Theranostics. Theranostics 2018, 8, 1648–1664. [Google Scholar] [CrossRef] [PubMed]
  85. Lin, H.; Gao, S.; Dai, C.; Chen, Y.; Shi, J. A Two-Dimensional Biodegradable Niobium Carbide (MXene) for Photothermal Tumor Eradication in NIR-I and NIR-II Biowindows. J. Am. Chem. Soc. 2017, 139, 16235–16247. [Google Scholar] [CrossRef] [PubMed]
  86. Han, X.; Jing, X.; Yang, D.; Lin, H.; Wang, Z.; Ran, H.; Li, P.; Chen, Y. Therapeutic mesopore construction on 2D Nb2C MXenes for targeted and enhanced chemo-photothermal cancer therapy in NIR-II biowindow. Theranostics 2018, 8, 4491–4508. [Google Scholar] [CrossRef]
  87. Feng, W.; Wang, R.; Zhou, Y.; Ding, L.; Gao, X.; Zhou, B.; Hu, P.; Chen, Y. Ultrathin Molybdenum Carbide MXene with Fast Biodegradability for Highly Efficient Theory-Oriented Photonic Tumor Hyperthermia. Adv. Funct. Mater. 2019, 29, 1901942. [Google Scholar] [CrossRef]
Materials 13 03083 g001 550
Figure 1. Feature engineering for dataset I obtained for two methods: (A) Random Forest (RF) and (B) Principle Component Analysis (PCA). Both methods are methodological tools that allow simplification of the model.
Figure 1. Feature engineering for dataset I obtained for two methods: (A) Random Forest (RF) and (B) Principle Component Analysis (PCA). Both methods are methodological tools that allow simplification of the model.
Materials 13 03083 g001
Materials 13 03083 g002 550
Figure 2. Ranking of feature importance obtained from RF. The most important are two descriptors: (i) the presence of MxOy, and (ii) Li on the surface, whereas the next four: surface modification, delaminating agent, lateral size, and Cl on the surface are equally important with smaller weights than the previous two. All feature labels are described in Table A2.
Figure 2. Ranking of feature importance obtained from RF. The most important are two descriptors: (i) the presence of MxOy, and (ii) Li on the surface, whereas the next four: surface modification, delaminating agent, lateral size, and Cl on the surface are equally important with smaller weights than the previous two. All feature labels are described in Table A2.
Materials 13 03083 g002
Materials 13 03083 g003 550
Figure 3. Feature engineering for dataset III obtained for two methods: (A) Random Forest (RF) and (B) Principle Component Analysis (PCA). Both methods are methodological tools that allow simplification of the model.
Figure 3. Feature engineering for dataset III obtained for two methods: (A) Random Forest (RF) and (B) Principle Component Analysis (PCA). Both methods are methodological tools that allow simplification of the model.
Materials 13 03083 g003
Materials 13 03083 g004 550
Figure 4. Ranking of feature importance for dataset III. The most important are the first three descriptors, namely MxOy, Li on the surface and surface modifications, respectively. All feature labels are described in Table A2 and Table A3.
Figure 4. Ranking of feature importance for dataset III. The most important are the first three descriptors, namely MxOy, Li on the surface and surface modifications, respectively. All feature labels are described in Table A2 and Table A3.
Materials 13 03083 g004
Table
Table 1. The accuracy score of correct predictions are obtained for unbalanced and balanced data, for the various algorithms employed in this paper. The most accurate results are obtained for balanced data by use of the Support Vector Machine with rbf kernel (regSVM-rbf).
Table 1. The accuracy score of correct predictions are obtained for unbalanced and balanced data, for the various algorithms employed in this paper. The most accurate results are obtained for balanced data by use of the Support Vector Machine with rbf kernel (regSVM-rbf).
Balanced Data
ML AlgorithmsUnbalanced DataWeightSMOTE
regRF0.7470.8260.833
regLOG-l10.7000.7760.783
regLOG-l20.7200.7980.783
regSVM-lin0.8670.7250.808
regSVM-rbf0.6620.9170.933
regSVM-sigmoid0.5550.5430.4042
ERT0.7220.7760.750
Table
Table 2. The metric of accuracy score of correct predictions is obtained, for various of algorithms employed in this paper.
Table 2. The metric of accuracy score of correct predictions is obtained, for various of algorithms employed in this paper.
ML AlgorithmsModel Selection: Weight
regRF0.845
regLOG-l10.845
regLOG-l20.727
regSVM-lin0.781
regSVM-rbf0.876
regSVM-sigmoid0.545
ERT0.793
Table
Table 3. Predicted probability of cytotoxicity score of various MXenes compounds. For each of the compounds different chemical treatment, as well as chemical composition on the surface have been reported. Thus, different probability scores have been obtained. All data taken from high-throughput screening experiments (see last column). The abbreviation Us. denotes Ultrasounds, whereas the the third and fourth columns denote the presence of the M x O y and Li atoms on the surface. Note, that two descriptors etching agent and delaminating agent are denoted here as the synthesis procedure.
Table 3. Predicted probability of cytotoxicity score of various MXenes compounds. For each of the compounds different chemical treatment, as well as chemical composition on the surface have been reported. Thus, different probability scores have been obtained. All data taken from high-throughput screening experiments (see last column). The abbreviation Us. denotes Ultrasounds, whereas the the third and fourth columns denote the presence of the M x O y and Li atoms on the surface. Note, that two descriptors etching agent and delaminating agent are denoted here as the synthesis procedure.
MXene [Ref.]:ProbabilityMxOyLi AtomsSynthesis ProcedureSurface Mod.
Ti3C2 [41]0.18NoYesLiF/HCl; Us.No
Ti3C2 [42]0.05NoNoHFAu
Ti3C2 [43]0.18NoYesLiF/HCl; Us.No
Ti3C2 [44]0.10NoNoHF; Us.No
Ti3C2 [45]0.10NoNoHF; Us.No
Ti3C2 [46]0.06NoYesLiF/HClAPTES + CEA
Ti3C2 [47]0.07NoYesLiF/HCl; Us.DNA, Pt, Pd
Ti3C2 [48]0.05NoNoHF; Us.Ag
Ti3C2 [49]0.07NoYesLiF/HClL-ACC
Ti3C2 [50]0.07NoYesLiF/HCl; Us.PAS
Ti3C2 [51]0.87YesNoLiF/HClNo
Ti3C2 [52]0.88YesNoLiF/HClNo
V2C [50]0.05NoYesLiF/HCl; Us.PAS
V2C [53]0.05NoNoNaF/HClNo
Nb2C [54]0.04NoNoNaF/HClNo
Nb2C [55]0.04NoNoHF; 60∘No
Ti2N [10]0.05NoNoKF/HCl; Us.No
Mo1.33C [56]0.04NoNoHF; TBAOHNo
Ti4N3 [57]0.03NoYesKF/LiF/NaF, TBAOHNo

Share and Cite

MDPI and ACS Style

Marchwiany, M.E.; Birowska, M.; Popielski, M.; Majewski, J.A.; Jastrzębska, A.M. Surface-Related Features Responsible for Cytotoxic Behavior of MXenes Layered Materials Predicted with Machine Learning Approach. Materials 2020, 13, 3083. https://doi.org/10.3390/ma13143083

AMA Style

Marchwiany ME, Birowska M, Popielski M, Majewski JA, Jastrzębska AM. Surface-Related Features Responsible for Cytotoxic Behavior of MXenes Layered Materials Predicted with Machine Learning Approach. Materials. 2020; 13(14):3083. https://doi.org/10.3390/ma13143083

Chicago/Turabian Style

Marchwiany, Maciej E., Magdalena Birowska, Mariusz Popielski, Jacek A. Majewski, and Agnieszka M. Jastrzębska. 2020. "Surface-Related Features Responsible for Cytotoxic Behavior of MXenes Layered Materials Predicted with Machine Learning Approach" Materials 13, no. 14: 3083. https://doi.org/10.3390/ma13143083

APA Style

Marchwiany, M. E., Birowska, M., Popielski, M., Majewski, J. A., & Jastrzębska, A. M. (2020). Surface-Related Features Responsible for Cytotoxic Behavior of MXenes Layered Materials Predicted with Machine Learning Approach. Materials, 13(14), 3083. https://doi.org/10.3390/ma13143083

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop