This research article introduces a novel operator termed q-convolution, strategically integrated with foundational principles of q-calculus. Leveraging this innovative operator alongside q-Bernoulli polynomials, a distinctive class of functions emerges, characterized by both analyticity and bi-univalence. The determination of initial coefficients within the Taylor-Maclaurin series for this function class is accomplished, showcasing precise bounds. Additionally, explicit computation of the second Hankel determinant is provided. These pivotal findings, accompanied by their corollaries and implications, not only enrich but also extend previously published results.