The functional dynamics of G protein-coupled receptors (GPCRs) encompasses multiple spatiotemporal scales, ranging from femtoseconds to seconds and Angstroms to micrometers. Computational approaches, often in close collaboration with experimental methods, have been invaluable in unraveling GPCR structure and dynamics at these various hierarchical levels. The binding of natural and synthetic ligands to the wild-type and naturally occurring variant receptors have been analyzed by several computational methods. The activation of receptors from the inactive to the active state has been investigated by atomistic simulations and ongoing work on several receptors will help uncover general and receptor-specific mechanisms. The interaction of GPCRs with complex membranes that contain phospholipids and cholesterol have been probed by coarse-grain methods and shown to directly influence receptor association. In this chapter, we discuss computational approaches that have been successful in analyzing each scale of GPCR dynamics. An overview of these approaches will allow a more judicious choice of the appropriate method. We hope that an appreciation of the power of current computational approaches will encourage more critical collaborations. A comprehensive integration of the different approaches over the entire spatiotemporal scales promises to unravel new facets of GPCR function.