We present an asymptotically optimal importance sampling for Monte Carlo simulation of the Laplace transform of exponential Brownian functionals which plays a prominent role in many disciplines. To this end we utilize the theory of large deviations to reduce finding an asymptotically optimal importance sampling measure to solving a calculus of variations problem. Closed-form solutions are obtained. In addition we also present a path to the test of regularity of optimal drift which is an issue in implementing the proposed method. The performance analysis of the method is provided through the Dothan bond pricing model.