, they do not need to be removed ; rather , under ideal conditions , the implant site repairs itself while the device is resorbed . 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 .
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 .
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 . 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 .
water uptake affects degradation , swelling , mechanical , and adhesive properties ; also it determines drug stability , drug release profile , and biological response .
current methods used to measure water uptake versus time are labor intensive and time consuming .
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 . .
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 .
published reviews that include the most relevant models currently available for important parameters in biomaterials such as glass transition temperature ( tg ) , young 's modulus , air - water contact angle , water uptake , and degradation .
( 2008 ) built a model of equilibrium water uptake for 12 aromatic polyamides with very similar levels of water uptake ( 13.9%19.1% ) .
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 , anomalous fickian diffusion , dual - stage fick 's diffusion , power law , weibull equation , langmuir theory , and concentration - dependent diffusion coefficient model   have been used .
modeling of hydration at the molecular level has been demonstrated using parameters such as free volume redistribution frequency , radial distribution functions ( rdfs ) , 3d atomic density maps known as spatial distribution functions , and angular distribution functions .
furthermore , from md simulations , water absorption has been predicted for a single polymer system [ 2830 ] .
.   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 .
[ 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 0.54  0.12 and 0.54  0.09 , respectively .
gubskaya et al .   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 0.67  0.13 .
the incorporation of three - dimensional descriptors led to important improvements in comparison with previous semiempirical models , increasing the average pearson correlation coefficient from 0.54  0.12 to 0.67  0.13 .
one of the challenges of biomaterials is the change of their interactions and properties over time .
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 .
previously , we built ann models to accurately predict drug release over time on a family of terpolymers   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 , such as bone pins , hernia repair devices , and an antibacterial sleeve that protects recipients of implanted cardiac assist devices from potentially life - threatening infections .
this library , consisting of a - b - type copolymers having an alternating sequence of a diphenol and a diacid , 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 tg and hydrophobicity .
in this study we investigate the effect of polymer backbone and pendent chain on the water uptake profiles of polymer films .
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 .
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 ) .
polymer films were compression molded and annealed at 510c above tg for 20  h before incubation , as described previously .
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 .
briefly , films 1  cm in diameter were incubated in h - radiolabeled water ( 0.2  ci / ml ) at 37c . 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 ( m3h2o ) was calculated using a calibration curve .
water uptake ( wu ) was calculated as the water content relative to the original dry weight ( msample):(1)wu%=100m3h2omsample.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 .
the data - mining package weka ( waikato environment for knowledge analysis )   was used in this study .
the methodology can be summarized in the following steps ( figure 1):polymers were characterized using two - dimensional ( 2d ) descriptors   and three - dimensional ( 3d ) descriptors .descriptors to build the model were selected using correlation based feature selection ( cfs ) , expectation - maximization ( em ) cluster analysis , decision tree analysis , and linear regression.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 .
polymers were characterized using two - dimensional ( 2d ) descriptors   and three - dimensional ( 3d ) descriptors .
descriptors to build the model were selected using correlation based feature selection ( cfs ) , expectation - maximization ( em ) cluster analysis , decision tree analysis , and linear regression .
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 .
descriptors based on the chemical structure of the polymers in implicit water or vacuum incorporating polymer conformation .
two - dimensional descriptors for the entire library of 112 polymers were obtained by smith et al . , using the basic molecular structure derived from the chemical formulae and both the molecular operating environment ( moe , chemical computing group inc . )   and the dragon ( milano chemometrics and qsar research group )   commercial software packages .
three - dimensional descriptors were obtained by gubskaya et al .   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 ( schrdinger )   commercial package with the generalized born / surface area implicit solvent model   and the opls - all atom force field .
although 3d descriptors obtained from tetramers do not capture the realistic structure of large mw polymers , they include very important information about their structure , which allows building more accurate models , as shown previously by gubskaya et al . .
similarly , other authors had previously used monomers   or less than 5 repeating monomeric units   to obtain molecular descriptors . starting with 2,272 descriptors taken from gubskaya et al .   and smith et al . ,
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 .
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 )  
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
the most significant descriptors were selected using a j48 decision tree , selecting descriptors that correctly partition the water uptake values and equilibrium water uptake according to the em cluster analysis . because decision tree analysis can not 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 .
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 .   and calculated by smith et al . .
several authors have shown that an ann model provides more accurate predictions than a linear model [ 5762 ] .
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 .
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
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 .
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 wueq .
2d descriptors include nct , hydrophilic factor , smr_vsa6 , ggi3 , mats3 m , and c-003 .
the hydrophilic factor is calculated from the number of hydrophilic groups ( -oh , -sh , and -nh ) of the molecule   and it was previously used to predict biological response on this polymer library .
smr_vsa6 is a descriptor of subdivided surface area , based on accessible van der waals surface area of each atom , and type of descriptor used before to predict fibrinogen adsorption of this polymer library .
ggi3 is a topological charge descriptor ; similar topological descriptors have been used to predict biological response on polymethacrylate surfaces .
mats3 m is a moran autocorrelation descriptor , which describes the level of correlation between molecules , and it has been used to study protein interactions .
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 .
g2 m is a whim descriptor , which captures relevant 3d information about molecular size , shape , symmetry , and atom distribution with respect to invariant reference frames .
r8p+ is r - getaway descriptor , which accounts for the local aspects of the molecule such as branching , cyclicity , and conformational changes .
mor25 m is a 3d - morse descriptor , which provides structural information of the molecules in the space , 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 , as defined by todeschini and consonni .
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 . .
results in table 4 show that correlation coefficient is not the best indicator of model accuracy . both models
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 wueq is able to predict 67% for training and 50% for test , within experimental variability .
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 wueq 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 wueq 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 ) .
this result is less accurate than predictions of simple physical behaviors such as tg , but it is much more accurate than predictions of more complex processes such as fibrinogen adsorption [ 34 , 36 ] , cell growth , and gene delivery efficiency , where the pearson coefficient for these models was below 0.77 . for each training and test set selection , predictions of equilibrium water uptake
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% ) .
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 .
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 can not 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 .
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
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 ,
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
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 .
 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 ( r2 > 0.78 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 ( r2 = 0.78 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
2. Materials and Methods
3. Results and Discussion
4. Conclusions