This paper introduces a novel subclass of bi-univalent analytic functions by utilizing a symmetric q-derivative operator in conjunction with Gegenbauer polynomials. Within this newly defined subclass, we derive bounds for the first two Maclaurin coefficients and address the Fekete-Szeg o problem. By varying the parameters in our results between 0 and 1, we obtain a range of new insights and rediscover some previously established results. This approach not only broadens the scope of bi-univalent function theory but also deepens the understanding of coefficient bounds and extremal problems within this context.