This study describes the dynamics of COVID-19 deaths MESHD and infections MESHD via a Monte Carlo approach. The analyses include death MESHD's data from USA, Brazil, Mexico, UK, India and Russia, which comprise the four countries with the highest number of deaths MESHD/ confirmed cases TRANS, as of Aug 07, 2020, according to the WHO. The Gompertz functions were fitted to the data of weekly averaged confirmed deaths MESHD per day by mapping the $\chi^2$ values. The uncertainties, variances and covariances of the model parameters were calculated by propagation. The fitted functions for the average deaths MESHD per day for USA and India have an upward trend, with the former having a higher growth rate and quite huge uncertainties. For Mexico, UK and Russia, the fits are consistent with a slope down pattern. For Brazil we found a subtle trend down, but with significant uncertainties. The USA, UK and India data shown a first peak with a higher growth rate when compared to the second one, demonstrating the benefits of non-pharmaceutical interventions of sanitary measures and social distance flattening the curve. For USA, a third peak seems quite plausible, most likely related with the recent relaxation policies. Brazil's data are satisfactorily described by two highly overlapped Gompertz functions with similar growth rates, suggesting a two-steps process for the pandemic spreading. The 95% CI for the total number of deaths MESHD ($\times 10^3$) predicted by the model for Aug 31, 2020 are 160 to 220, 110 to 130, 59 to 62, 46.6 to 47.3, 54 to 63 and 16.0 to 16.7 for USA, Brazil, Mexico, UK, India and Russia, respectively. Our estimates for the prevalences SERO of infections MESHD are in reasonable agreement with some preliminary reports from serological studies carried out in USA and Brazil. The method represents an effective framework to estimate the line-shape of the infection MESHD curves and the uncertainties of the relevant parameters based on the actual data.