Owing to their atomic thickness and low bending rigidity, suspended two-dimensional (2D) materials are prone to wrinkle formation. Here, we use molecular dynamics (MD) simulations to probe the effect of these wrinkles on the nonlinear elasticity of atomically thin graphene membranes. We observe a stress-strain response that consists of two linear regions that are separated by a transition. It is found that this transition is sharp in membranes where wrinkles are formed by uneven stresses at the boundaries. However, when wrinkles are formed by crystal defects, this nonlinear transition is seen to be more gradual. To capture these effects, we use a phenomenological model based on experimentally measurable quantities. We demonstrate the model's fidelity by fitting it to the MD simulated nonlinear response of many graphene membranes providing evidence that the sharpness of the transition between the linear regions in the stress-strain response is a measure of the type of wrinkles and can be quantified by our model.