In vitro rate constants (ka, kdes, ki, and kdeg) were optimized using the GA, originally developed by Holland in the early 1970s. This was typically used as an artificial intelligence algorithm [e.g., it is used for training artificial neural networks (65)] for a more robust optimization of parameters based on evolutionary ideas of natural selection where it does not rely heavily on initial input and often leads to a global minimum (66). For optimization, rate constants were supplied by the GA to produce simulation outputs at particular time points. When fed to the model equations (described above), membrane ([Mem](t)) and cell space ([Cell](t)) simulated outputs were summed at each time point and compared to in vitro measured values using the residual sum of squares described below. To supply model equations (Eqs. 11 and 12) with rate constants, initial vectors (“chromosomes”) composed of rate constants (“genes”) were randomly populated by the GA, fed to the cell kinetic model, calculated for fitness, and underwent selection, crossover, and mutations to maximize diversity and produce better fitness at each iteration (“generation”). Specifically, Eqs. 11 and 12 were first optimized to the calibrated datasets where no degradation was present, holding kdeg at 0. Once optimized, we reconsidered QD raw concentration values containing degradation effects to determine the rate of degradation (kdeg), holding the previously optimized adsorption, desorption, and internalization rates (kads, kdes, and kint) as constant. All simulations were performed in MATLAB v2015b. Parameter optimization was implemented with the GA optimization function from the Optimization Toolbox. Parameters for estimation included the following: initial population, 300; population size, 50; generations, 100; mutation rate, mutation Gaussian; crossover rate, 0.80; selection function, stochastic uniform.

The GA was evaluated using the residual sum squares as the fitness function (Eq. 19)RSS=in(yimi)2(19)where RSS represented the residual sum of squares from model output (mi) at time (i) to measured data (yi) for n time points. Standard error was computed asS=RSSn(20)where S is the standard error, RSS is the residual sum of squares, and n is the total time points.

Model output upper and lower bounds were evaluated at the 95% confidence interval throughCL(95%)=mi±2*S(21)where CL(95%) represented 95% confidence limit. The GA was run for 100 generations, enough to allow convergence at a fitness value representative of measured data.

Note: The content above has been extracted from a research article, so it may not display correctly.



Q&A
Please log in to submit your questions online.
Your question will be posted on the Bio-101 website. We will send your questions to the authors of this protocol and Bio-protocol community members who are experienced with this method. you will be informed using the email address associated with your Bio-protocol account.



We use cookies on this site to enhance your user experience. By using our website, you are agreeing to allow the storage of cookies on your computer.