Large-eddy simulations of atmospheric boundary layers under various stability and surface roughness conditions are performed to investigate the turbulence impact on wind turbines. In particular, the aeroelastic responses of the turbines are studied to characterize the fatigue loading of the turbulence present in the boundary layer and in the wake of the turbines. Two utility-scale 5 MW turbines that are separated by seven rotor diameters are placed in a 3 km by 3 km by 1 km domain. They are subjected to atmospheric turbulent boundary layer flow and data are collected on the structural response of the turbine components. The surface roughness was found to increase the fatigue loads while the atmospheric instability had a small influence. Furthermore, the downstream turbines yielded higher fatigue loads indicating that the turbulent wakes generated from the upstream turbines have significant impact.