Hybrid nanofluid is an emerging field due to the rapid enhancement of heat transfer and stable nanoparticles in base fluid properties. A three dimensional hybrid nanofluid flow model is constructed over biaxial porous stretching/shrinking sheet with heat transfer, Radiative heat and mass flux (3D-HNF-RHF). Bayesian Regularization technique based on Backpropagated neural networks (BRT-BNNs) is employed to estimate the solution of proposed model. It is the mathematical process which converts nonlinear regression fitting into a well-posed statistical process in the form of ridge regression. This method diminishes length cross validation need and more robust then Backpropagation technique. The proposed flow system 3D-HNF-RHF is transferred to ordinary differential equations (ODEs) possessing physical variations through self-similar transformations. The effect of derived variations such as thermal relaxation parameter, mass flux parameter, Stretching/shrinking parameter, Prandtl number, Skin friction and Nusselt number observed over the velocity and temperature fields. Numerical results of these physical parameters have been presented in tabulated form obtained from a dataset constructed through Homotopy Analysis Method imposed on 3D-HNF-RHF model. Statistical tests through mean square error, histogram curve and regression fitting curves are employed to check the accuracy and convergences of the solution obtained through BRT-BNNs. It is observed that Stretching/shrinking quantity and mass flux parameter slow down the flow rate whereas, increase radiative heat flux upsurges the temperature gradient. The impacts of surface drag force and heat transfer are illustrated through the different graphical illustrations.