In this study, we developed a fractional-order model for polycystic ovarian syndrome issue transmission in women through four compartments with the Mittag-Leffler kernel. For biological feasibility, we discussed the positivity and boundedness of solutions with fixed point mappings in Banach spaces. Sensitivity analysis is performed using different variables as response functions each time, employing Latin Hyper-cube Sampling (LHS) and Partial Rank Correlation Coefficient (PRCC) with 200 iterations. Further, the reproductive number is determined for the model’s biological feasibility with an analysis of parameter impact through Matlab software. The solutions are examined for global stability using the Lyapunov function with a first derivative. Chaos Control will use the regulate for linear responses approach to bring the system to stabilize according to its points of equilibrium. Numerical solutions are derived by using a Newton polynomial in the generalized form of the Mittag-Leffler kernel, which illustrates the effect of the fractional operator and parameters involved in the system. Results are also analyzed with different assumptions of parameter values and initial conditions to observe the different impacts of treatment to overcome the risk of infertility in women. Such type of investigation will be useful to investigate the spread of disease as well as helpful in developing control strategies to overcome the risk of problems and best management care pre-detection of the issue and for treatment.