In this paper, stochastic predictive computing networks are exploited to investigate the dynamics of the SIS with vaccination impact based epidemic model (SISV-EM) represented by nonlinear systems of stochastic differential equations (SDEs) by exploitation of artificial neural networks (ANNs) with the backpropagated Levenberg-Marquardt technique (BLMT) i.e., (ANNs-BLMT) to approximate the solution behavior. The stochastic nonlinear SISV-EM is governed with three classes: susceptible, infectious, and vaccinated populations. The referenced or target datasets for ANNs-BLMT are constructed by employing Euler-Maruyama (EM) scheme for solving stochastic differential systems in case of sufficiently various nonlinear SISV-EM scenarios by varying the percentage of vaccination for newly born, the coefficient of transmission, the natural mortality rates, the infectious rates of recovery, the rate at which vaccinated people lose their immunity, the rate of death caused by disease, the proportion of vaccinated against susceptible and the white noise in the environment. Based on arbitrary training, testing, and validation samples from the referenced dataset, the ANNs-BLMT provides an approximate solution for the stochastic nonlinear SISV-EM, with significant correlations to the referenced results. Exhaustive simulation-based results using error histograms, mean square errors, and regression analyses further demonstrate that the proposed ANNs-BLMT is efficient, consistent, and accurate for solving SISV-EM. • A two-layer framework of ANNs-BLMT is proposed as an innovative technique based on a stochastic computing paradigm to investigate the dynamics of stochastic nonlinear SISV-EM. • Innovative technique based on a stochastic computing ANNs-BLMT to investigate the dynamics of stochastic nonlinear SISV-EM Reference dataset for ANNs -BLMT is developed with Euler-Maruyama (EM) scheme nonlinear SISV-EM with varying parametersMean square error criteria is used for training of ANNs-BLMT to find approximate solutions to a variety of SISV-EM scenarios Computation of error histogram illustration, regression metrics, and MSE learning curves prove the performance of ANNs-BLMT. • The mean square error criteria are effectively utilized in approximation theory to develop an objective function for training the ANNs-BLMT to determine approximate solutions to a variety of stochastic nonlinear SISV-EM scenarios. • The computation of error histogram illustration, regression metrics, and MSE learning curves significantly improve the performance, precision, and consistency of the ANNs-BLMT for solving the stochastic nonlinear SISV-EM. [ABSTRACT FROM AUTHOR]