Nonlinear mixed effects models structural

Various methods are available to estimate population parameters, but today the nonlinear mixed effects modeling approach is the most common one employed. Population analyses have been performed for mAbs such as basiliximab, daclizu-mab and trastuzumab, as well as several others in development, including clenolixi-mab and sibrotuzumab. Population pharmacokinetic models comprise three submodels the structural the statistical and covariate submodels (Fig. 3.13). Their development and impact for mAbs will be discussed in the following section. 

The number of samples per subject used for this approach is typically small, ranging from one to six. As does the pooled analysis technique, nonlinear mixed-effects modeling approaches analyze the data of all individuals at once, but take the interindividual random effects structure into account. This ensures that confounding correlations and imbalance that may occur in observational data are properly accounted for. 

Data structure analysis is the examination of the raw data for hidden structure, outliers, or leverage observations. This is repeated during the exploratory modeling (and nonlinear mixed effects modeling) steps using case deletion diagnostics (20). This type of analysis is important since outliers or leverage observations may occur in a population PM data set. It is equally important for the reduction of the covariate vector. 

The knowledge discovery basis of PM modeling permits the generation of hypotheses from the relationship discovered during data structure analysis. These relationships can be tested in the nonlinear mixed effects modeling step. It can... 

Nonlinear mixed effects models consist of two components the structural model (which may or may not contain covariates) and the statistical or variance model. The structural model describes the mean response for the population. Similar to a linear mixed effects model, nonlinear mixed effects models can be developed using a hierarchical approach. Data consist of an independent sample of n-subjects with the ith subject having -observations measured at time points t i, t 2, . t n . Let Y be the vector of observations, Y = Y1 1, Yi,2,. ..Ynjl,Yn,2,. ..Yn,ni)T and let s... 

Returning to the previous structural model, a complete nonlinear mixed effects model may be written as... 

Residual variance models were described in detail in the chapter on variance models and weighting. What was detailed in that chapter readily extends to nonlinear mixed effects models. As might be expected, residual variance models model the random, unexplained variability in the regression function f. Hence, the structural model is extended to... 

To get at the question of overall influence, the matrix of structural model parameters and variance components was subjected to principal component analysis. Principal component analysis (PCA) was introduced in the chapter on Nonlinear Mixed Effects Model Theory and transforms a matrix of values to another matrix such that the columns of the transformed matrix are uncorrelated and the first column contains the largest amount of variability, the second column contains the second largest, etc. Hopefully, just the first few principal components contain the majority of the variance in the original matrix. The outcome of PC A is to take X, a matrix of p-variables, and reduce it to a matrix of q-variables (q < p) that contain most of the information within X. In this PC A of the standardized parameters (fixed effects and all variance components), the first three principal components contained 74% of the total variability in the original matrix, so PCA was largely successfully. PCA works best when a high correlation exists between the variables in the original data set. Usually more than 80% variability in the first few components is considered a success. 

An analytical structure-(hyper)polarizability relationship based on a two-state description has also been derived [49]. In this model a parameter MIX is introduced that describes the mixture between the neutral and charge-separated resonance forms of donor-acceptor substituted conjugated molecules. This parameter can be directly related to BLA and can explain solvent effects on the molecular hyperpolarizabilities. NMR studies in solution (e.g. in CDCl3) can give an estimate of the BLA and therefore allow a direct correlation with the nonlinear optical experiments. A similar model introducing a resonance parameter c that can be related to the MIX parameter was also introduced to classify nonlinear optical molecular systems [50,51]. 

The forcefields discussed in Section 2.1 use energy functions which do not take into account quantum effects. Many important processes are intrinsically quantum mechanical, and thus cannot be modelled classically. SA has been used in con-juction with density functional theory, the Schrodinger equation, chemical reaction dynamics, electronic structure studies, and to optimize linear and nonlinear parameters in trial wave functions. This is important because quantum effects are often embedded in an essentially classical system. This has motivated mixing the classical fields with the quantum potentials in simulations known as quantum mechanic/molecular mechanic hybrids. Including quantum effects is important in the study of enzyme reactions, and proton and electron transport studies. 

---