In this study, the interaction potential of water molecule with a graphene oxide (GO) plate containing OH and O groups has been calculated using the M06-2X/6-31g(d,p) level of theory at different orientations and intermolecular distances and fitted to the Born-Huggins-Meyer (BHM) model. There are good agreements between the calculated and the OPLS-AA and Dreiding models, especially for the GO(O)-H2O interactions. To examine the new computed models, we have used the closer potentials to the OPLS-AA and Dreiding models in the molecular dynamics (MD) simulations. We have calculated several properties using the different obtained interaction potentials including average number of hydrogen bonds per water molecule (〈HB〉) between confined water molecules and between water and GO-surfaces, radial distribution function (RDF), self-diffusion coefficient, and angle distribution function of the confined water molecules between the GO plates. Our results showed good agreements between the OPLS-AA and Dreiding models and some calculated potentials. However, some calculated models showed completely different behavior which discussed in details. According to the results, we concluded that the OH-2 and O1-OH2 models show totally better agreement with the famous force fields than the other calculated potentials. This work provides a simple method for the development of new force fields specifically for these types of systems which are in good agreement with the well-known force fields.