An integral part of any systems biology approach is the modelling and simulation of the respective system under investigation. However, the values of many parameters of the system have often not been determined or are not identifiable due to technical experimental difficulties or other constraints. Sensitivity analysis is often employed to quantify the importance of each of the model's parameters in the behaviour of the system. This approach can also be useful in identifying those parts of the system that are most sensitive with the potential of becoming drug targets. A problem of the commonly used methods of sensitivity analysis is that they constitute local methods meaning that they depend directly on the exact parameter space, which in turn is not known exactly. One way to circumvent this problem is to carry out sensitivity analysis over a wide range of values for all parameters, but this is handicapped by expensive computations when the systems are high dimensional. Another approach is to employ global sensitivity analysis, which in this context is mostly based on random sampling methods. In this paper we present an efficient approach that involves using numerical optimizing methods that search a wide region of parameter space for a given model to determine the maximum and minimum values of its metabolic control coefficients. A relevant example for drug development is presented to demonstrate the strategy using the software COPASI.