This study presents a mathematical model for reversible enzymatic reactions using a hybrid proportional fractional derivative. The model is based on chemical kinetics and enzyme dynamics, representing a two-step substrate-enzyme reversible reaction. Non-linear differential equations are transformed into fractional-order systems using the constant proportional Caputo–Fabrizio (CPCF) and constant proportional Atangana–Baleanu–Caputo (CPABC) operators. The system is simulated using the Laplace–Adomian decomposition method at different fractional orders. Qualitative and quantitative analyses, including boundedness, positivity, unique solution, and feasible concentration, are provided for the proposed model with different hybrid operators. The stability analysis of the proposed scheme is verified using Picard's stable condition through the fixed point theorem. The study highlights the importance of enzyme kinetics in understanding the rates of enzyme-catalyzed reactions and their applications in biochemistry and medicine. It also discusses the use of fractional calculus in modeling complex systems, including enzyme reactions, and the advantages of fractional-order models over traditional integer-order models. The study reviews recent developments in fractional calculus, including the use of different fractional derivatives and operators in modeling biological and chemical processes. The results show that the proposed model can accurately predict the dynamics of reversible enzymatic reactions and provide insights into the behavior of enzymes in metabolic processes. The study also discusses the application of fractional calculus in various fields, including physics, engineering, and biology, and the potential of fractional-order models in improving the accuracy of predictions in complex systems.This study presents a mathematical model for reversible enzymatic reactions using a hybrid proportional fractional derivative. The model is based on chemical kinetics and enzyme dynamics, representing a two-step substrate-enzyme reversible reaction. Non-linear differential equations are transformed into fractional-order systems using the constant proportional Caputo–Fabrizio (CPCF) and constant proportional Atangana–Baleanu–Caputo (CPABC) operators. The system is simulated using the Laplace–Adomian decomposition method at different fractional orders. Qualitative and quantitative analyses, including boundedness, positivity, unique solution, and feasible concentration, are provided for the proposed model with different hybrid operators. The stability analysis of the proposed scheme is verified using Picard's stable condition through the fixed point theorem. The study highlights the importance of enzyme kinetics in understanding the rates of enzyme-catalyzed reactions and their applications in biochemistry and medicine. It also discusses the use of fractional calculus in modeling complex systems, including enzyme reactions, and the advantages of fractional-order models over traditional integer-order models. The study reviews recent developments in fractional calculus, including the use of different fractional derivatives and operators in modeling biological and chemical processes. The results show that the proposed model can accurately predict the dynamics of reversible enzymatic reactions and provide insights into the behavior of enzymes in metabolic processes. The study also discusses the application of fractional calculus in various fields, including physics, engineering, and biology, and the potential of fractional-order models in improving the accuracy of predictions in complex systems.