This paper presents a simple and novel least squares fitting-based modeling approach for the extraction simple gene regulatory sub- networks from biclusters in microarray time series gene expression data. Preprocessing helps in retaining the strongly interacting gene regulatory pairs. The methodology was applied to public-domain data sets of Yeast and the experimental results were biologically validated based on standard databases and information from literature. © 2009 Springer-Verlag Berlin Heidelberg.