This paper introduces a new 3D chaotic attractor in a gene regulatory network. The proposed model has eighteen parameters. Formulas for characteristic numbers of critical points for three-dimensional systems were considered. We show that the three equilibrium points of the new chaotic 3D system are unstable and deduce that the three-dimensional system exhibits chaotic behavior. The possible outcomes of this 3D model were compared with the results of the Chua circuit. The bifurcation structures of the proposed 3D system are investigated numerically, showing periodic solutions and chaotic solutions. Lyapunov exponents and Kaplan-Yorke dimension are calculated. For calculations, the Wolfram Mathematica is used.