Targeted drug delivery using functionalized nanocarriers (NCs) is a strategy in therapeutic and diagnostic applications. In this paper we review the recent development of models at multiple length and time scales and their applications to targeting of antibody functionalized nanocarriers to antigens (receptors) on the endothelial cell (EC) surface. Our mesoscale (100 nm-1 μm) model is based on phenomenological interaction potentials for receptor-ligand interactions, receptor-flexure and resistance offered by glycocalyx. All free parameters are either directly determined from independent biophysical and cell biology experiments or estimated using molecular dynamics simulations. We employ a Metropolis Monte Carlo (MC) strategy in conjunction with the weighted histogram analysis method (WHAM) to compute the free energy landscape (potential of mean force or PMF) associated with the multivalent antigen-antibody interactions mediating the NC binding to EC. The binding affinities (association constants) are then derived from the PMF by computing absolute binding free energy of binding of NC to EC, taking into account the relevant translational and rotational entropy losses of NC and the receptors. We validate our model predictions by comparing the computed binding affinities and PMF to a wide range of experimental measurements, including in vitro cell culture, in vivo endothelial targeting, atomic force microscopy (AFM), and flow chamber experiments. The model predictions agree closely and quantitatively with all types experimental measurements. On this basis, we conclude that our computational protocol represents a quantitative and predictive approach for model driven design and optimization of functionalized NCs in targeted vascular drug delivery.