The localized deformation of molecular monolayers constrained between the spherical surfaces of Au nanoparticles is studied by means of molecular dynamics simulations. Alkyl or polyethylene glycol long-chain molecules were homogeneously distributed over the curved Au surface, pushed against each other by repeated cycles of force relaxation and constant-volume equilibration at temperatures increasing from 50 to 300 K before being slowly quenched to near-zero temperature. Plots of minimum configurational energy can be obtained as a function of the nanoparticle distance, according to different directions of approach; therefore, such simulations describe a range of deformations, from perfectly uniaxial compression to a combination of compression and shear. Despite the relative rigidity of molecular backbones, the deformation is always found to be localized at the interface between the opposing molecular monolayers. We find that shorter ligands can be more densely packed on the surface but do no interdigitate upon compression; they respond to the applied force by bending and twisting, thus changing their conformation while remaining disjointed. On the contrary, longer ligands attain lower surface densities and can interprenetrate when the nanoparticles are compressed against each other; such molecules remain rather straight and benefit from the increased overlap to maximize the adhesion by dispersion forces. The apparent Young's and shear moduli of a dense nanostructure, composed of a triangular arrangement of identical MUDA-decorated Au nanoparticles, are found to be smaller than estimates indirectly deduced by atomic-force experiments but quite close to previous computer simulations of molecular monolayers on flat surfaces and of bulk nanoparticle assemblies.