Plant metabolic engineering is commonly used in the production of functional foods and quality trait improvement. However, to date, computational model-based approaches have only been scarcely used in this important endeavor, in marked contrast to their prominent success in microbial metabolic engineering. In this study we present a computational pipeline for the reconstruction of fully compartmentalized tissue-specific models of Arabidopsis thaliana on a genome scale. This reconstruction involves automatic extraction of known biochemical reactions in Arabidopsis for both primary and secondary metabolism, automatic gap-filling, and the implementation of methods for determining subcellular localization and tissue assignment of enzymes. The reconstructed tissue models are amenable for constraint-based modeling analysis, and significantly extend upon previous model reconstructions. A set of computational validations (i.e., cross-validation tests, simulations of known metabolic functionalities) and experimental validations (comparison with experimental metabolomics datasets under various compartments and tissues) strongly testify to the predictive ability of the models. The utility of the derived models was demonstrated in the prediction of measured fluxes in metabolically engineered seed strains and the design of genetic manipulations that are expected to increase vitamin E content, a significant nutrient for human health. Overall, the reconstructed tissue models are expected to lay down the foundations for computational-based rational design of plant metabolic engineering. The reconstructed compartmentalized Arabidopsis tissue models are MIRIAM-compliant and are available upon request.C urrent challenges in using plants as factories for bio-energy and nutraceuticals require predesigned and efficient strategies for metabolic engineering (1). Currently, plant metabolic engineering mostly involves trial-and-error approaches, without the utilization of computational modeling procedures to rationally design genetic modifications. The marginal contribution played by metabolic modeling in plants until now stands in marked contrast to its prominent success in microbial metabolic engineering (2, 3). Metabolic network reconstructions were manually reconstructed for dozens of bacterial species (3), and automated approaches recently generated draft models for a total of 130 bacteria (4). A modeling approach, called constraintbased modeling (CBM), serves to analyze the function of such large-scale metabolic networks by solely relying on simple physical-chemical constraints, overcoming the problem of missing enzyme kinetic data (3). Applications of CBM for large-scale microbial networks has proven to be highly successful in predicting metabolic phenotypes in metabolic engineering and many other applications (3).The reconstruction of metabolic network models for multicellular eukaryotes is significantly more challenging than that for bacteria, because of the larger size of the networks, the subcellular compartmentalization of me...