We present a method to calculate (31)P solid-state NMR spectra based on the dynamic input from extended molecular dynamics (MD) simulations. The dynamic information confered by MD simulations is much more comprehensive than the information provided by traditional NMR dynamics models based on, for example, order parameters. Therefore, valuable insight into the dynamics of biomolecules may be achieved by the present method. We have applied this method to study the dynamics of lipid bilayers containing the antimicrobial peptide alamethicin, and we show that the calculated (31)P spectra obtained with input from MD simulations are in agreement with experiments under a large range of different sample conditions, including vesicles and oriented samples with and without peptides. We find that the changes in the (31)P spectra upon addition of peptide stem from lipids with reduced diffusion due to peptide-lipid interactions.