International Journal of Biomaterials
Volume 2016 (2016), Article ID 6273414, 10 pages
http://dx.doi.org/10.1155/2016/6273414

Research Article

Developing a Suitable Model for Water Uptake for Biodegradable Polymers Using Small Training Sets

Loreto M. Valenzuela,¹ Doyle D. Knight,² and Joachim Kohn³

¹Department of Chemical and Bioprocess Engineering, Research Center for Nanotechnology and Advanced Materials “CIEN-UC”, Pontificia Universidad Católica de Chile, Vicuña Mackenna 2860, Macul, 7820436 Santiago, Chile
²Department of Mechanical and Aerospace Engineering, Rutgers, The State University of New Jersey, New Brunswick, NJ 08854-8087, USA
³New Jersey Center for Biomaterials, Rutgers, The State University of New Jersey, 145 Bevier Road, Piscataway, NJ 08854, USA

Received 21 December 2015; Accepted 3 April 2016

Academic Editor: Rosalind Labow

Copyright © 2016 Loreto M. Valenzuela et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Prediction of the dynamic properties of water uptake across polymer libraries can accelerate polymer selection for a specific application. We first built semiempirical models using Artificial Neural Networks and all water uptake data, as individual input. These models give very good correlations ( for test set) but very low accuracy on cross-validation sets (less than 19% of experimental points within experimental error). Instead, using consolidated parameters like equilibrium water uptake a good model is obtained ( for test set), with accurate predictions for 50% of tested polymers. The semiempirical model was applied to the 56-polymer library of L-tyrosine-derived polyarylates, identifying groups of polymers that are likely to satisfy design criteria for water uptake. This research demonstrates that a surrogate modeling effort can reduce the number of polymers that must be synthesized and characterized to identify an appropriate polymer that meets certain performance criteria.

1. Introduction

Degradable materials are very important in fabricating biomedical devices. After implantation, they do not need to be removed; rather, under ideal conditions, the implant site repairs itself while the device is resorbed [1]. In comparison, nondegradable materials often need to be surgically removed after their purpose has been achieved, thus subjecting the patient to a second surgery that potentially exposes them to more complications [2]. Degradable devices can be used in a broad range of applications such as vascular stents, vascular bypass grafts, bone fixation devices, and soft tissue replacement scaffolds [3].

Degradable biomaterials have a wide range of requirements depending on the particular clinical application. Parameters such as chemical structure, composition, porosity, and device geometry determine surface and bulk properties of an implant, and thus, they are critical to the selection of the material [4].

One important characteristic of degradable biomaterials is their water uptake versus time, as it is crucial for the determination of how long a polymeric device will reside in the body before erosion leads to the ultimate removal of the device from the implant site [5]. Water uptake affects degradation, swelling, mechanical [6], and adhesive properties [7]; also it determines drug stability [8], drug release profile [9], and biological response [10].

Current methods used to measure water uptake versus time are labor intensive and time consuming. Depending on the polymer, water uptake can take days to weeks to equilibrate [11]. There are potentially very large libraries of polymeric biomaterials, which make it impractical to measure these parameters experimentally for each polymer. For example, a virtual library of about 40,000 polymethacrylates has been described by Kholodovych et al. [12]. This library would clearly be too large for each polymer to be characterized individually by experimental methods.

Computational modeling is a useful tool to minimize the number of experiments needed to characterize a polymer library [13]. Costache et al. [14], Gubskaya [15], and Le et al. [16] published reviews that include the most relevant models currently available for important parameters in biomaterials such as glass transition temperature (), Young’s modulus, air-water contact angle, water uptake, and degradation. Serna et al. (2008) built a model of equilibrium water uptake for 12 aromatic polyamides with very similar levels of water uptake (13.9%–19.1%). They found correlations between the amidic hydrogen charge and the water uptake [17].

Although empirical mathematical modeling has been successfully used to model water uptake for different polymers, all models require parameters that can only be obtained through experimentation. Fick’s diffusion [18], anomalous Fickian diffusion [19], dual-stage Fick’s diffusion [20], power law [18], Weibull equation [21], Langmuir theory [22], and concentration-dependent diffusion coefficient model [23] have been used. Modeling of hydration at the molecular level has been demonstrated using parameters such as free volume redistribution frequency [24], Radial Distribution Functions (RDFs) [25], 3D atomic density maps known as spatial distribution functions [26], and angular distribution functions [27]. Furthermore, from MD simulations, water absorption has been predicted for a single polymer system [28–30].

Prior works by Kholodovych et al. [12, 31], Smith et al. [32–35], Gubskaya et al. [36], and Ghosh et al. [37] showed that it is possible to build computational models of polymer properties for an entire library based upon experimental data for a small subset. In these studies, a polymer library is explored using a combined experimental and computational approach, looking for polymers that fulfill a series of design criteria to be suitable for specific applications. Smith et al. [33, 34] developed semiempirical models using molecular descriptors obtained from two-dimensional polymer structures (i.e., the descriptors were independent of the polymer conformation). These models were able to predict fibrinogen adsorption within experimental error in 38 out of the 45 polymers and rat lung fibroblast proliferation in 41 out of 48 polymers. Pearson correlation coefficient values for these predictions were and , respectively. Gubskaya et al. [36] calculated descriptors from relaxed three-dimensional polymeric structures obtained from Molecular Dynamics (MD) simulations of tetramers in vacuum and implicit water. In this work, Decision Tree Analysis and ANNs were used to predict fibrinogen adsorption with a Pearson coefficient of . The incorporation of three-dimensional descriptors led to important improvements in comparison with previous semiempirical models, increasing the average Pearson correlation coefficient from to .

One of the challenges of biomaterials is the change of their interactions and properties over time [38]. However, all aforementioned models study and predict individual values for each polymer. They do not consider dynamic properties that may change over time. Even Le et al. (2013), who built predictions of phase behavior over time, developed the model using each experimental value as a single input, without considering how the phase behavior changes over time [39]. Previously, we built ANN models to accurately predict drug release over time on a family of terpolymers [40] using molecular descriptors. In this study, we develop and compare models for water uptake over time, first using all individual data separately and then using a global parameter for this property.

Our research has two objectives: (i) the development of computational models for water uptake versus time based upon experimental data from a small subset of polymers in a library and (ii) the application of these models to predict water uptake for an entire library of polymers. The main challenge of this research is to model and predict properties that change over time with particular kinetics using a small set of experimental data. As a model system, a library of L-tyrosine-derived polyarylates was used. Kohn and collaborators used this library to discover promising lead polymers for several medical applications [41], such as bone pins [42], hernia repair devices, and an antibacterial sleeve that protects recipients of implanted cardiac assist devices from potentially life-threatening infections [43].

This library, consisting of A-B-type copolymers having an alternating sequence of a diphenol and a diacid [41], was obtained by copolymerizing 14 tyrosine-derived diphenols with 8 aliphatic diacids in all possible combinations resulting in 112 distinct polymers. Changes in polymer backbone or pendent chain length affect polymer properties such as and hydrophobicity. In this study we investigate the effect of polymer backbone and pendent chain on the water uptake profiles of polymer films.

2. Materials and Methods

2.1. Materials

A subset of the L-tyrosine-derived polyarylates was synthesized as described previously by carbodiimide-mediated solution polycondensation of a diphenol and a diacid at room temperature [44].

2.1.1. Nomenclature

DTR = desaminotyrosyl-tyrosine alkyl ester: R = methyl (M), ethyl (E), iso-propyl (iP), butyl (B), iso-butyl (iB), sec-butyl (sB), hexyl (H), octyl (O), dodecyl (D), benzyl (Bn), 2-(2-ethoxyethoxy)ethyl (G).

HTR = hydroxyacetic acid-tyrosine alkyl ester: R = ethyl (E), hexyl (H), octyl (O).

2.2. Experimental Methods

2.2.1. Film Processing

Polymer films were compression molded and annealed at 5–10°C above for 20 h before incubation, as described previously [11].

2.2.2. Water Uptake

Water uptake was obtained for the selected polymers from the L-tyrosine-derived polyarylates combinatorial library (Table 1) using ³H-labeled water, as described previously [45]. Briefly, films 1 cm in diameter were incubated in ³H-radiolabeled water (0.2 μCi/mL) at 37°C. After 6 h and 12 h and 1, 2, 3, 4, 7, 14, 21, 28, 35, and 42 days, samples were removed from the vial, rinsed with distilled water, blotted dry, and dissolved with 3 mL of tetrahydrofuran (THF) (VWR) and 12 mL of liquid scintillation cocktail (LSC) (Ecolite). Radioactive counts were measured using a scintillation counter (Beckmann 6500), and water content () was calculated using a calibration curve. Water uptake (WU) was calculated as the water content relative to the original dry weight ():

Table 2 lists the estimated values for equilibrium water uptake from the experimental measurements; both this parameter and individual water uptake experimental points were used to build surrogate models for water uptake.

Table 1: Subset of the library of L-tyrosine-derived polyarylates used in this study.

Polymer
(kDa)
(°C)
Polymer set for model
Predictions

Poly(DTO sebacate)

16
●

Poly(DTB adipate)

42
●

Poly(DTO succinate)

43
●

Poly(DTE adipate)

59
●

Poly(DTE glutarate)

64
●

Poly(DTB succinate)

67
●

Poly(HTH sebacate)

23
●

Poly(HTH adipate)

40
●

Poly(DTM sebacate)

45
●

Poly(DTiP adipate)

55
●

Poly(DTM adipate)

67
●

Poly(HTE succinate)
★
78
●

Poly(DTO adipate)

26
●

Poly(DTsBR(+)methyladipate)

45
●

Poly(DTsB R(+) glutarate)

46
●

Poly(DTM R(+) methyladipate)

53
●

Poly(DTBn adipate)

61
●

Poly(HTE adipate)

61
●

Poly(DTO suberate)

21

●

Poly(DTH suberate)

24

●

Poly(HTH suberate)

27

●

Poly(DTO glutarate)

32

●

Poly(DTiB sebacate)

33

●

Poly(DTH R(+) methyladipate)

33

●

Poly(DTH L(−) methyladipate)

33

●

Poly(DTH adipate)

34

●

Poly(DTB R(+) methyladipate)

35

●

Poly(DTB L(−) methyladipate)

35

●

Poly(DTB suberate)

37

●

Poly(DTO diglycolate)

40

●

Poly(DTBn sebacate)

42

●

Poly(DTH glutarate)

43

●

Poly(DTH diglycolate)

45

●

Poly(DTsB L(−) methyladipate)

45

●

Poly(DTsB L(−) glutarate)

46

●

Poly(DTsB R(+) suberate)

46

●

Poly(DTsB L(−) suberate)

46

●

Poly(DTsB R(+) adipate)

50

●

Poly(DTsB L(−) adipate)

50

●

Poly(DTB glutarate)

50

●

Poly(DTH succinate)

53

●

Poly(DTM L(−) methyladipate)

53

●

Poly(HTE suberate)

54

●

Poly(DTiP R(+) methyladipate)

54

●

Poly(DTiP L(−) methyladipate)

54

●

Poly(DTM suberate)

55

●

Poly(DTBn R(+) methyladipate)

55

●

Poly(DTBn L(−) methyladipate)

55

●

Poly(DTiB adipate)

56

●

Poly(DTE R(+) methyladipate)

63

●

Poly(DTE L(−) methyladipate)

63

●

Poly(HTE R(+) methyladipate)

63

●

Poly(HTE L(−) methyladipate)

63

●

Poly(DTB diglycolate)

64

●

Poly(DTiB succinate)

75

●

The “” symbol indicates the presence of more than one chiral center in the polymer repeat unit.

Molecular weight ( was measured by THF-GPC (mean value of three different films ± standard deviation (SD)).

The “” symbol indicates the polymers that did not dissolve in THF and, thus, could not be measured, and degradation could not be measured.

Glass transition temperature ( was measured by DSC for the dry polymer before pressing.

Table 2: Equilibrium water uptake for 18 polymers of the L-tyrosine-derived polyarylate library.

Polymer
Equilibrium water uptake (%)

Poly(DTB adipate)

Poly(DTB succinate)

Poly(DTBn adipate)

Poly(DTE adipate)

Poly(DTE glutarate)

Poly(DTiP adipate)

Poly(DTM adipate)

Poly(DTM sebacate)

Poly(DTO adipate)

Poly(DTO sebacate)

Poly(DTO succinate)

Poly(HTE adipate)

Poly(HTE succinate)

Poly(HTH adipate)

Poly(HTH sebacate)

Poly(DTM R(+) methyladipate)

Poly(DTsB R(+) glutarate)

Poly(DTsB R(+) methyladipate)

Polymers are ordered by name used in the descriptor set.

2.3. Computational Methods

The data-mining package WEKA (Waikato Environment for Knowledge Analysis) [46] was used in this study. The methodology can be summarized in the following steps (Figure 1):

(i)Polymers were characterized using two-dimensional (2D) descriptors [32] and three-dimensional (3D) descriptors [36].(ii)Descriptors to build the model were selected using correlation based feature selection (CFS), expectation-maximization (EM) cluster analysis, Decision Tree Analysis, and linear regression.(iii)Either all water uptake experimental data points over time or equilibrium water uptake was used to build the model using ANNs, using 10% for testing and the rest for training.

Figure 1: Scheme of experimental method for surrogate models of water uptake.

2.3.1. Descriptors

The descriptors in this study include “2D” descriptors based on the chemical structure of the polymers [32] and “3D” descriptors based on the chemical structure of the polymers in implicit water or vacuum incorporating polymer conformation [36]. Two-dimensional descriptors for the entire library of 112 polymers were obtained by Smith et al. [34], using the basic molecular structure derived from the chemical formulae and both the Molecular Operating Environment (MOE, Chemical Computing Group Inc.) [47] and the Dragon (Milano Chemometrics and QSAR Research Group) [48] commercial software packages. Three-dimensional descriptors were obtained by Gubskaya et al. [36] for 56 polymers from the polyarylate library. Descriptors were obtained by the Dragon commercial software package using the 3D structures of the tetramers after structure minimization and 1 ns of MD simulations using MacroModel v.8.5 (Schrödinger) [49] commercial package with the generalized Born/surface area implicit solvent model [50] and the OPLS-all atom force field [51]. Although 3D descriptors obtained from tetramers do not capture the realistic structure of large polymers, they include very important information about their structure, which allows building more accurate models, as shown previously by Gubskaya et al. [36]. Similarly, other authors had previously used monomers [52] or less than 5 repeating monomeric units [53] to obtain molecular descriptors.

2.3.2. Descriptor Selection

Starting with 2,272 descriptors taken from Gubskaya et al. [36] and Smith et al. [32], a correlation based feature selection (CFS) was used to reduce the dimensionality of the descriptors for each parameter in study. CFS is a function available in WEKA that evaluates the worth of a subset of attributes (descriptors) by considering the individual predictive ability of each feature along with the degree of redundancy between them. As a result, it selects a subset of attributes that are highly correlated with the parameter while removing irrelevant, redundant, and noisy attributes [54]. A genetic search algorithm was used in conjunction with the CFS, allowing a parallel search of the attribute space and avoiding local optima.

For each model, expectation-maximization (EM) [46] cluster analysis was employed to categorize the polymer property of study (i.e., water uptake and equilibrium water uptake) into three classes (i.e., low, medium, and high). When analyzing all data points for water uptake, both time and water uptake values were included in the cluster analysis.

The most significant descriptors were selected using a J48 Decision Tree [55], selecting descriptors that correctly partition the water uptake values and equilibrium water uptake according to the EM cluster analysis. Because Decision Tree Analysis cannot represent relationships between continuous variables, an additional descriptor was selected by linear regression, that is, the highest weight on the linear regression, for the full training set and the experimental values of water uptake. Time was also included as a descriptor for water uptake with all data points.

2.3.3. Artificial Neural Networks

Linear models are insufficient to capture the complexity of the structure-property-relationships between polymer structure and water uptake profiles. Specifically, we observed that water uptake does not yield a simple correlation with the hydrophilic factor, as defined by Todeschini et al. [56] and calculated by Smith et al. [32].

Several authors have shown that an ANN model provides more accurate predictions than a linear model [57–62]. A multilayer perceptron (MLP) was used to build ANN models for each parameter with the three descriptors selected as explained in Section 2.3.2. Two hidden layers (nodes) were used. Output nodes were unthresholded linear units [46]. Backpropagation by gradient descent was used as MLP learning method. All input variables were scaled to the unit interval while the learning rate and the momentum applied for updating the weights were 0.3 and 0.2, respectively. Training time was set on 1,000 epochs, which showed to be enough for model convergence. To perform cross-validation, 10% of data was separated as test set in each model, in all possible combinations. Randomization of the initial weights and shuffling of the training data were performed by varying the seed for the random number generator. The model obtained with each seed represents a local optimum, based on the initial weights. Thus, running enough seeds and selecting the best model among them would allow finding the global optimum. For the present models, a hundred ANN models were obtained with different seeds, from which the best model in terms of root mean squared error for the training set was selected.

3. Results and Discussion

3.1. Descriptors Selection

Table 3 summarizes the descriptors selected for both models. One 3D descriptor and five 2D descriptors were selected for the model for all time points; two 3D descriptors and one 2D were selected for the model of . 2D descriptors include nCt, hydrophilic factor, SMR_VSA6, GGI3, MATS3m, and C-003. nCt is the number of tertiary carbon atoms (sp3). The hydrophilic factor is calculated from the number of hydrophilic groups (-OH, -SH, and -NH) of the molecule [63] and it was previously used to predict biological response on this polymer library [34]. SMR_VSA6 is a descriptor of subdivided surface area, based on accessible van der Waals surface area of each atom [64], and type of descriptor used before to predict fibrinogen adsorption of this polymer library [35]. GGI3 is a topological charge descriptor; similar topological descriptors have been used to predict biological response on polymethacrylate surfaces [37]. MATS3m is a Moran autocorrelation descriptor, which describes the level of correlation between molecules, and it has been used to study protein interactions [65]. C-003 is the number of CHR3 molecular subfragments, an atom center fragment; it gives information about structural motifs important for the molecular shape and it was used before to predict fibrinogen adsorption on polymethacrylate surfaces [37].

Table 3: Best descriptors and their variability within the training set and within the complete set of 56 polymers.

Model
Descriptor
SD for polymers of the model
SD for the complete library

All data points
Hydrophilic factor
0.246
0.212

SMR_VSA6
0.291
0.242

GGI3
0.227
0.264

MATS3m
0.256
0.273

C-003
0.394
0.478

G2m vacuum
0.231
0.255

nCt
0.287
0.316

Mor25m water
0.212
0.238

R8p+ vacuum
0.243
0.242

3D descriptors include G2m and R8p+ in vacuum and Mor25m in water. G2m is a WHIM descriptor, which captures relevant 3D information about molecular size, shape, symmetry, and atom distribution with respect to invariant reference frames [66]. WHIM descriptors were used to predict fibrinogen adsorption on polymethacrylate surfaces [37]. R8p+ is R-GETAWAY descriptor, which accounts for the local aspects of the molecule such as branching, cyclicity, and conformational changes [67].

Mor25m is a 3D-MoRSE descriptor, which provides structural information of the molecules in the space [68], and it has been suggested that this information is related to the free volume of molecules [69, 70] and, thus, responsible for the ability of the polymer to uptake water. 3D-MoRSE and GETAWAY descriptors have been also correlated with the tendency of a molecule to be solvated by water, measured by the hydrophilic index Hy [71], as defined by Todeschini and Consonni [63]. These types of descriptors encode relevant information of this polymer library that gives information of several physical and chemical processes such as water uptake and even in fibrinogen adsorption as discussed by Gubskaya et al. [36].

3.2. Model for Water Uptake

Results in Table 4 show that correlation coefficient is not the best indicator of model accuracy. Both models present high of training set (>0.92). However, the model using all data presents only 17% or less of predictions within experimental variability, for training and test sets, while the model for is able to predict 67% for training and 50% for test, within experimental variability.

Table 4: Summary of models for water uptake.

Model
n training set
Number of descriptors
training
Within experimental variability (training)
cross-validation (10%)
Within experimental variability (test)

All data points
189
6+ time
0.92
30/189 (16%)
0.83
3/18 (17%)

18
3
0.97
12/18 (67%)
0.78
9/18 (50%)

Results of cross-validation have to be analyzed very carefully when using all data points, because they are not independent of each other. In that case, it is likely to select for cross-validation data that belong to polymers for which there is a large data set in the training set. Thus, depending on how the cross-validation set is selected, different results will be obtained.

On the other hand, the model for obtains its values from several experimental measurements of each polymer after its water content is equilibrated. This gives more representative and reliable experimental data, and it captures more information than single points at the same time of incubation. With this, cross-validation that in this case includes only independent values, considering all possible combinations of leave-two-out (10%) of experimental values, gives accurate predictions in 50% of the cases from test sets, and was correctly classified as high, medium, or low according to the EM cluster analysis previously done, in 83% of the cases. With this, predictions accurately represent the relative order in water uptake of the polymers studied (Figure 2 and Table 4).

Figure 2: Prediction versus experimental values for for polymers as part of training () and test () sets. Black line represents . Values are presented as mean value ± SD of predictions (-error) ± SD of experimental values (-error).

This result is less accurate than predictions of simple physical behaviors such as [14], but it is much more accurate than predictions of more complex processes such as fibrinogen adsorption [34, 36], cell growth [34], and gene delivery efficiency [72], where the Pearson coefficient for these models was below 0.77.

3.3. Predictions of Water Uptake over Rest of the Library

For each training and test set selection, predictions of equilibrium water uptake were made for the rest of the 56-polymer library. As Figure 3 shows, the model predicts low levels of water uptake for polymers containing DTM (0%–14%), DTO (0%–25%), and HTH (5%–14%) (with the exception of methyladipates); low to intermediate levels of water uptake for polymers containing DTBn (18%–34%) and DTE (35%–37%) (with the exception of methyladipates), glutarate (0%–37%), suberate (0%–26%) (with the exception of DTsB), and sebacate (5%–32%) (with the exception of DTiB); intermediate levels of water uptake for succinate-containing polymers (13%–61%); medium to high levels of water uptake for DTiB-containing polymers (82%–120%); high levels of water uptake for DTiP methyladipates (111%–139%); and widely ranging levels of water uptake for DTH (0%–87%), adipate (5%–96%), and methyladipate (31%–139%) polymers. It also predicts that all DTB polymers have low values of water uptake (less than 36%) and only high values of water uptake for DTsB polymers (92%–135%); it predicts low values of water uptake (10%–26%) for HTE polymers (with the exception of methyladipate) and predicts low levels of water uptake for diglycolate polymers (0%–18%).

Figure 3: Predictions of equilibrium water uptake over the remaining 38 polymers of the polymer library. Values are presented as mean value ± SD of the predicted value for each training/test set combination. Polymers are ordered from highest to the lowest water uptake predicted values. Solid lines separate areas of very high, high, medium, and low water uptake polymers.

Some of these predictions would be expected directly from the chemical structure of the polymers, but others would not be easily expected. For example, all DTO polymers would have low water uptake, which is expected from the long pendant chain (8 carbons), while the DTH polymers, with only one carbon less than the DTO, would have water uptake levels from low to medium.

3.4. Limitations of Surrogate Modeling

Limitations of this type of model include the following: (i) it needs experimental data to train the model; (ii) the descriptors give a reference of relevant parameters to the target property, but they cannot explain the mechanism; (iii) experimental measurements must be performed to validate the predictions; (iv) for new polymers outside of the sublibrary, new descriptors must be generated, which is time consuming due to the need for MD simulations. However, this last limitation is only encountered for the first property that you wish to model, for once the descriptors are generated, they can be used to build predictive models for several properties of the polymer library. The obtained models for water uptake can be improved by increasing the size of training set, by generating more meaningful descriptors, such as 3D descriptors in explicit water, by improving the descriptor selection algorithm, and by identifying other surrogate methods.

4. Conclusions

This study describes a new approach to modeling dynamic properties and demonstrates the potential value of this approach. In particular, we developed models for water uptake for a library of polymers using only a small training set and molecular descriptors for all the polymers in the library. We also demonstrate that using a consolidated parameter of water uptake, a dynamic property, gives a more accurate model than using a more conventional approach of all experimental measurement as independent values. By separating time points from one experiment, information about the slope, rate, and progression of the dynamic property is not considered in the model. And, since data points are not independent, accuracy of predictions is compromised.

A surrogate model was built to accurately predict equilibrium water uptake of a polymer sublibrary of 56 L-tyrosine-derived polyarylates using a small training set and only three descriptors selected from a large set of descriptors, calculated from either 2D or 3D structures. Those descriptors included atom counts; 3D information about electron diffraction (3D-MoRSE); and chemical properties of molecular atoms, branching, cyclicity, and conformational changes (GETAWAY). Although these descriptors can be used only for this model in this polymer library, the methodology for selecting descriptors can be applied to any polymer library and/or polymer property.

The model was able to accurately predict low, intermediate, and high levels of water uptake for up to 12 of the 18 polymers. Using this model, predictions were obtained for the rest of the sublibrary. Those predictions can be used primarily as a reference of order of magnitude and ranking of polymers in terms of water uptake.

Finally, having several semiempirical models for different polymer properties such as glass transition temperature, contact angle, fibrinogen adsorption, cell response, water uptake, and degradation for the same polymer library may be used to select a polymer for a specific application. With a known set of design criteria, a group of polymers can be selected from the mentioned models. After this selection, the actual parameters must be measured experimentally, the models must be validated, and the best polymer can be selected to begin the device development process. With this, surrogate modeling of polymer properties may accelerate the discovery and selection of rationally designed materials for a target application.

Disclosure

The content of this paper is solely the responsibility of the authors and does not necessarily represent the official views of the NIH, NIBIB, NCMHD, or CONICYT.

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

This work was supported by RESBIO (Integrated Technology Resource for Polymeric Biomaterials) funded by National Institutes of Health (NIBIB and NCMHD) under Grant P41 EB001046 and by the New Jersey Center for Biomaterials and CONICYT (FONDECYT 11121392).

References

J. Kohn, W. J. Welsh, and D. Knight, “A new approach to the rationale discovery of polymeric biomaterials,” Biomaterials, vol. 28, no. 29, pp. 4171–4177, 2007.

J. M. Lantry, C. S. Roberts, and P. V. Giannoudis, “Operative treatment of scapular fractures: a systematic review,” Injury, vol. 39, no. 3, pp. 271–283, 2008.

L. G. Griffith, “Polymeric biomaterials,” Acta Biomaterialia, vol. 48, no. 1, pp. 263–277, 2000.

N. Angelova and D. Hunkeler, “Rationalizing the design of polymeric biomaterials,” Trends in Biotechnology, vol. 17, no. 10, pp. 409–421, 1999.

M. S. Shoichet, “Polymer scaffolds for biomaterials applications,” Macromolecules, vol. 43, no. 2, pp. 581–591, 2010.

H. Kranz, N. Ubrich, P. Maincent, and R. Bodmeier, “Physicomechanical properties of biodegradable poly(D,L-lactide) and poly(D,L-lactide-co-glycolide) films in the dry and wet states,” Journal of Pharmaceutical Sciences, vol. 89, no. 12, pp. 1558–1566, 2000.

C. Hopkins, P. E. McHugh, N. P. O'Dowd, Y. Rochev, and J. P. McGarry, “A combined computational and experimental methodology to determine the adhesion properties of stent polymer coatings,” Computational Materials Science, vol. 80, pp. 104–112, 2013.

C. Ahlneck and G. Zografi, “The molecular basis of moisture effects on the physical and chemical stability of drugs in the solid state,” International Journal of Pharmaceutics, vol. 62, no. 2-3, pp. 87–95, 1990.

D. Caccavo, S. Cascone, G. Lamberti, and A. A. Barba, “Modeling the drug release from hydrogel-based matrices,” Molecular Pharmaceutics, vol. 12, no. 2, pp. 474–483, 2015.

M. Tanaka, T. Hayashi, and S. Morita, “The roles of water molecules at the biointerface of medical polymers,” Polymer Journal, vol. 45, no. 7, pp. 701–710, 2013.

L. M. Valenzuela, B. Michniak, and J. Kohn, “Variability of water uptake studies of biomedical polymers,” Journal of Applied Polymer Science, vol. 121, no. 3, pp. 1311–1320, 2011.

V. Kholodovych, A. V. Gubskaya, M. Bohrer et al., “Prediction of biological response for large combinatorial libraries of biodegradable polymers: polymethacrylates as a test case,” Polymer, vol. 49, no. 10, pp. 2435–2439, 2008.

D. C. Webster and M. A. R. Meier, “Polymer libraries: preparation and applications,” in Polymer Libraries, vol. 225 of Advances in Polymer Science, pp. 1–15, Springer, Berlin, Germany, 2010.

A. D. Costache, J. Ghosh, D. D. Knight, and J. Kohn, “Computational methods for the development of polymeric biomaterials,” Advanced Engineering Materials, vol. 12, no. 1-2, pp. B3–B17, 2010.

A. V. Gubskaya, “Quantum-chemical descriptors in QSAR/QSPR modeling: achievements, perspectives and trends,” in Quantum Biochem, C. F. Matta, Ed., pp. 693–721, Wiley-VCH, Weinheim, Germany, 2010.

T. Le, V. C. Epa, F. R. Burden, and D. A. Winkler, “Quantitative structure-property relationship modeling of diverse materials properties,” Chemical Reviews, vol. 112, no. 5, pp. 2889–2919, 2012.

F. Serna, F. García, J. L. De La Peña, and J. M. García, “Aromatic polyisophthalamides with mononitro, dinitro and trinitroiminobenzoyl pendant groups,” High Performance Polymers, vol. 20, no. 1, pp. 19–37, 2008.

J. Crank, The Mathematics of Diffusion, Clarendon Press, Oxford, UK, 2nd edition, 1975.

S. Roy, W. X. Xu, S. J. Park, and K. M. Liechti, “Anomalous moisture diffusion in viscoelastic polymers: modeling and testing,” Journal of Applied Mechanics, vol. 67, no. 2, pp. 391–396, 1999.

W. K. Loh, A. D. Crocombe, M. M. A. Wahab, and I. A. Ashcroft, “Modelling anomalous moisture uptake, swelling and thermal characteristics of a rubber toughened epoxy adhesive,” International Journal of Adhesion and Adhesives, vol. 25, no. 1, pp. 1–12, 2005.

W. Weibull, “A statistical distribution of wide applicability,” Journal of Applied Mechanics, vol. 18, pp. 293–297, 1951.

S. Popineau, C. Rondeau-Mouro, C. Sulpice-Gaillet, and M. E. R. Shanahan, “Free/bound water absorption in an epoxy adhesive,” Polymer, vol. 46, no. 24, pp. 10733–10740, 2005.

S. Joannès, L. Mazé, and A. R. Bunsell, “A simple method for modeling the concentration-dependent water sorption in reinforced polymeric materials,” Composites Part B: Engineering, vol. 57, pp. 219–227, 2014.

A. Noorjahan and P. Choi, “Effect of free volume redistribution on the diffusivity of water and benzene in poly(vinyl alcohol),” Chemical Engineering Science, vol. 121, pp. 258–267, 2015.

Y. Tamai, H. Tanaka, and K. Nakanishi, “Molecular dynamics study of polymer-water interaction in hydrogels. 2. Hydrogen-bond dynamics,” Macromolecules, vol. 29, no. 21, pp. 6761–6769, 1996.

A. V. Gubskaya and P. G. Kusalik, “Molecular dynamics simulation study of ethylene glycol, ethylenediamine, and 2-aminoethanol. 1. The local Structure in pure liquids,” Journal of Physical Chemistry A, vol. 108, no. 35, pp. 7151–7164, 2004.

J. Behler, D. W. Price, and M. G. B. Drew, “Water structuring properties of carbohydrates, molecular dynamics studies on 1,5-anhydro-D-fructose,” Physical Chemistry Chemical Physics, vol. 3, no. 4, pp. 588–601, 2001.

M. Canales, D. Aradilla, and C. Alemán, “Water absorption in polyaniline emeraldine base,” Journal of Polymer Science Part B: Polymer Physics, vol. 49, no. 18, pp. 1322–1331, 2011.

T.-X. Xiang and B. D. Anderson, “Distribution and effect of water content on molecular mobility in poly(vinylpyrrolidone) glasses: a molecular dynamics simulation,” Pharmaceutical Research, vol. 22, no. 8, pp. 1205–1214, 2005.

T. X. Xiang and B. D. Anderson, “Water uptake, distribution, and mobility in amorphous poly(D,L-lactide) by molecular dynamics simulation,” Journal of Pharmaceutical Sciences, vol. 103, no. 9, pp. 2759–2771, 2014.

V. Kholodovych, J. R. Smith, D. Knight, S. Abramson, J. Kohn, and W. J. Welsh, “Accurate predictions of cellular response using QSPR: a feasibility test of rational design of polymeric biomaterials,” Polymer, vol. 45, no. 22, pp. 7367–7379, 2004.

J. R. Smith, V. Kholodovych, D. Knight, W. J. Welsh, and J. Kohn, “QSAR models for the analysis of bioresponse data from combinatorial libraries of biomaterials,” QSAR & Combinatorial Science, vol. 24, no. 1, pp. 99–113, 2005.

J. R. Smith, V. Kholodovych, D. Knight, J. Kohn, and W. J. Welsh, “Predicting fibrinogen adsorption to polymeric surfaces in silico: a combined method approach,” Polymer, vol. 46, no. 12, pp. 4296–4306, 2005.

J. R. Smith, A. Seyda, N. Weber, D. Knight, S. Abramson, and J. Kohn, “Integration of combinatorial synthesis, rapid screening, and computational modeling in biomaterials development,” Macromolecular Rapid Communications, vol. 25, no. 1, pp. 127–140, 2004.

J. R. Smith, D. Knight, J. Kohn et al., “Using surrogate modeling in the prediction of fibrinogen adsorption onto polymer surfaces,” Journal of Chemical Information and Computer Sciences, vol. 44, no. 3, pp. 1088–1097, 2004.

A. V. Gubskaya, V. Kholodovych, D. Knight, J. Kohn, and W. J. Welsh, “Prediction of fibrinogen adsorption for biodegradable polymers: integration of molecular dynamics and surrogate modeling,” Polymer, vol. 48, no. 19, pp. 5788–5801, 2007.

J. Ghosh, D. Y. Lewitus, P. Chandra et al., “Computational modeling of in vitro biological responses on polymethacrylate surfaces,” Polymer, vol. 52, no. 12, pp. 2650–2660, 2011.

S. W. Cranford, J. de Boer, C. van Blitterswijk, and M. J. Buehler, “Materiomics: an -omics approach to biomaterials research,” Advanced Materials, vol. 25, no. 6, pp. 802–824, 2013.

T. C. Le, C. E. Conn, F. R. Burden, and D. A. Winkler, “Computational modeling and prediction of the complex time-dependent phase behavior of lyotropic liquid crystals under in meso crystallization conditions,” Crystal Growth and Design, vol. 13, no. 3, pp. 1267–1276, 2013.

A. V. Gubskaya, I. J. Khan, L. M. Valenzuela, Y. V. Lisnyak, and J. Kohn, “Investigating the release of a hydrophobic peptide from matrices of biodegradable polymers: an integrated method approach,” Polymer, vol. 54, no. 15, pp. 3806–3820, 2013.

S. Brocchini, K. James, V. Tangpasuthadol, and J. Kohn, “A combinatorial approach for polymer design,” Journal of the American Chemical Society, vol. 119, no. 19, pp. 4553–4554, 1997.

K. A. Hooper, N. D. Macon, and J. Kohn, “Comparative histological evaluation of new tyrosine-derived polymers and poly (L-lactic acid) as a function of polymer degradation,” Journal of Biomedical Materials Research, vol. 41, no. 3, pp. 443–454, 1998.

H. L. Bloom, L. Constantin, D. Dan et al., “Implantation success and infection in cardiovascular implantable electronic device procedures utilizing an antibacterial envelope,” Pacing and Clinical Electrophysiology, vol. 34, no. 2, pp. 133–142, 2011.

J. Fiordeliso, S. Bron, and J. Kohn, “Design, synthesis, and preliminary characterization of tyrosine-containing polyarylates: new biomaterials for medical applications,” Journal of Biomaterials Science, vol. 5, no. 6, pp. 496–510, 1994.

L. M. Valenzuela, G. Zhang, C. R. Flach et al., “Multiscale analysis of water uptake and erosion in biodegradable polyarylates,” Polymer Degradation and Stability, vol. 97, no. 3, pp. 410–420, 2012.

I. H. Witten and E. Frank, Data Mining: Practical Machine Learning Tools and Techniques with JAVA Implementations, Academic Press, San Diego, Calif, USA, 1st edition, 2000.

Chemical Computing Group, MOE (The Molecular Operating Environment), Chemical Computing Group, Montreal, Canada, 2nd edition, 2003.

R. Todeschini, V. Consonni, A. Mauri, and M. Pavan, Dragon Web Version, 3.0, Milano, Italy, 2003.

Schrödinger, Schrödinger Release 2005: MacroModel, V. 8.5, Schrödinger, New York, NY, USA, 2005.

W. C. Still, A. Tempczyk, R. C. Hawley, and T. Hendrickson, “Semianalytical treatment of solvation for molecular mechanics and dynamics,” Journal of the American Chemical Society, vol. 112, no. 16, pp. 6127–6129, 1990.

W. L. Jorgensen, D. S. Maxwell, and J. Tirado-Rives, “Development and testing of the OPLS all-atom force field on conformational energetics and properties of organic liquids,” Journal of the American Chemical Society, vol. 118, no. 45, pp. 11225–11236, 1996.

A. P. Toropova, A. A. Toropov, V. O. Kudyshkin, D. Leszczynska, and J. Leszczynski, “Optimal descriptors as a tool to predict the thermal decomposition of polymers,” Journal of Mathematical Chemistry, vol. 52, no. 5, pp. 1171–1181, 2014.

P. R. Duchowicz, S. E. Fioressi, D. E. Bacelo, L. M. Saavedra, A. P. Toropova, and A. A. Toropov, “QSPR studies on refractive indices of structurally heterogeneous polymers,” Chemometrics and Intelligent Laboratory Systems, vol. 140, pp. 86–91, 2015.

M. A. Hall, Correlation-Based Feature Selection for Machine Learning, The University of Waikato, 1999.

J. R. Quinlan, Programs for Machine Learning, Morgan Kaufmann Publishers, San Francisco, Calif, USA, 1st edition, 1993.

R. Todeschini, C. Bettiol, G. Giurin, P. Gramatica, P. Miana, and E. Argese, “Modeling and prediction by using WHIM descriptors in QSAR studies: Submitochondrial particles (SMP) as toxicity biosensors of chlorophenols,” Chemosphere, vol. 33, no. 1, pp. 71–79, 1996.

A. Afantitis, G. Melagraki, K. Makridima, A. Alexandridis, H. Sarimveis, and O. Iglessi-Markopoulou, “Prediction of high weight polymers glass transition temperature using RBF neural networks,” Journal of Molecular Structure, vol. 716, no. 1–3, pp. 193–198, 2005.

A. T. Seyhan, G. Tayfur, M. Karakurt, and M. Tanoglu, “Artificial neural network (ANN) prediction of compressive strength of VARTM processed polymer composites,” Computational Materials Science, vol. 34, no. 1, pp. 99–105, 2005.

J. W. Gao, X. Y. Wang, X. B. Li, X. Yu, and H. Wang, “Prediction of polyamide properties using quantum-chemical methods and BP artificial neural networks,” Journal of Molecular Modeling, vol. 12, no. 4, pp. 513–520, 2006.

W. Q. Liu, P. G. Yi, and Z. L. Tang, “QSPR models for various properties of polymethacrylates based on quantum chemical descriptors,” QSAR & Combinatorial Science, vol. 25, no. 10, pp. 936–943, 2006.

F. Gharagheizi, “QSPR analysis for intrinsic viscosity of polymer solutions by means of GA-MLR and RBFNN,” Computational Materials Science, vol. 40, no. 1, pp. 159–167, 2007.

J. Xu, H. Liang, B. Chen, W. Xu, X. Shen, and H. Liu, “Linear and nonlinear QSPR models to predict refractive indices of polymers from cyclic dimer structures,” Chemometrics and Intelligent Laboratory Systems, vol. 92, no. 2, pp. 152–156, 2008.

R. Todeschini and V. Consonni, Handbook of Molecular Descriptors, Wiley-VCH, 1st edition, 2000.

S. A. Wildman and G. M. Crippen, “Prediction of physicochemical parameters by atomic contributions,” Journal of Chemical Information and Computer Sciences, vol. 39, no. 5, pp. 868–873, 1999.

J.-F. Xia, K. Han, and D.-S. Huang, “Sequence-based prediction of protein-protein interactions by means of rotation forest and autocorrelation descriptor,” Protein and Peptide Letters, vol. 17, no. 1, pp. 137–145, 2010.

R. Todeschini and P. Gramatica, “New 3D molecular descriptors: the WHIM theory and QSAR applications,” Perspectives in Drug Discovery and Design, vol. 9, pp. 355–380, 1998.

V. Consonni, R. Todeschini, and M. Pavan, “Structure/response correlations and similarity/diversity analysis by GETAWAY descriptors. 1. Theory of the novel 3D molecular descriptors,” Journal of Chemical Information and Computer Sciences, vol. 42, no. 3, pp. 682–692, 2002.

J. Gasteiger, J. Schuur, P. Selzer, L. Steinhauer, and V. Steinhauer, “Finding the 3D structure of a molecule in its IR spectrum,” Fresenius' Journal of Analytical Chemistry, vol. 359, no. 1, pp. 50–55, 1997.

W. Liu, “Prediction of glass transition temperatures of aromatic heterocyclic polyimides using an ANN model,” Polymer Engineering and Science, vol. 50, no. 8, pp. 1547–1557, 2010.

B. E. Mattioni and P. C. Jurs, “Prediction of glass transition temperatures from monomer and repeat unit structure using computational neural networks,” Journal of Chemical Information and Computer Sciences, vol. 42, no. 2, pp. 232–240, 2002.

Z. Jelcic, “Solvent molecular descriptors on poly(D, L-lactide-co-glycolide) particle size in emulsification—diffusion process,” Colloids and Surfaces A: Physicochemical and Engineering Aspects, vol. 242, no. 1–3, pp. 159–166, 2004.

A. V. Gubskaya, T. O. Bonates, V. Kholodovych et al., “Logical analysis of data in structure-activity investigation of polymeric gene delivery,” Macromolecular Theory and Simulations, vol. 20, no. 4, pp. 275–285, 2011.