Analysis of chemical graphs is becoming a major research topic in computational molecular biology due to its potential applications to drug design. One of the major approaches in such a study is inverse quantitative structure activity/property relationships (inverse QSAR/QSPR) analysis, which is to infer chemical structures from given chemical activities/properties. Recently, a novel framework has been proposed for inverse QSAR/QSPR using both artificial neural networks (ANN) and mixed integer linear programming (MILP). This method consists of a prediction phase and an inverse prediction phase. In the first phase, a feature vector f (G) of a chemical graph G is introduced and a prediction function ψ N on a chemical property π is constructed with an ANN N . In the second phase, given a target value y * of the chemical property π, a feature vector x * is inferred by solving an MILP formulated from the trained ANN N so that ψ N (x * ) is close to y * and then a set of chemical structures G * such that f (G * ) = x * is enumerated by a graph search algorithm. The framework has been applied to the case of chemical compounds with cycle index up to 2 so far. The computational results conducted on instances with n non-hydrogen atoms show that a feature vector x * can be inferred for up to around n = 40 whereas graphs G * can be enumerated for up to around n = 15. When applied to the case of chemical acyclic graphs, the maximum computable diameter of G * was around up to around 8. In this paper, we introduce a new characterization of graph structure, called "branch-height" based on which a new MILP formulation and a new graph search algorithm are designed for chemical acyclic graphs. The results of computational experiments using such chemical properties as octanol/water partition coefficient, boiling point and heat of combustion suggest that the proposed method can infer chemical acyclic graphs G * with n = 50 and diameter 30.