Citrate capping is one of the most common strategies to achieve the colloidal stability of Au nanoparticles (NPs) with diameters ranging from a few to hundreds of nanometers. Citrate-capped Au nanoparticles (CNPs) represent a step of the synthesis of Au NPs with specific functionalities, as CNPs can be further functionalized via ligand-exchange reactions, leading to the replacement of citrate with other organic ligands. In vitro, CNPs are also used to address the fundamental aspects of NP-membrane interactions, as they can directly interact with cells or model cell membranes. Their affinity for the bilayer is again mediated by the exchange of citrate with lipid molecules. Here, we propose a new computational model of CNPs compatible with the coarse grained Martini force field. The model, which we develop and validate through an extensive comparison with new all-atom molecular dynamics (MD) simulations and UV-vis and Fourier transform infrared spectroscopy data, is aimed at the MD simulation of the interaction between citrate-capped NPs and model phosphatidylcholine lipid membranes. As a test application we show that, during the interaction between a single CNP and a flat planar 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine bilayer, the citrate coating is spontaneously replaced by lipids on the surface of Au NPs, while the NP size and shape determine the final structural configuration of the NP-bilayer complex.