Rhodopsin is the only G protein-coupled receptor (GPCR) whose 3D structure is known; therefore, it serves as a prototype for studies of the GPCR family of proteins. Rhodopsin dysfunction has been linked to misfolding, caused by chemical modifications that affect the naturally occurring disulfide bond between C110 and C187. Here, we identify the structural elements that stabilize rhodopsin by computational analysis of the rhodopsin structure and comparison with data from previous in vitro mutational studies. We simulate the thermal unfolding of rhodopsin by breaking the native-state hydrogen bonds sequentially in the order of their relative strength, using the recently developed Residues most stable under thermal denaturation are part of a core, which is assumed to be important for the formation and stability of folded rhodopsin. This core includes the C110OC187 disulfide bond at the center of residues forming the interface between the transmembrane and the extracellular domains near the retinal binding pocket. Fast mode analysis of rhodopsin using the Gaussian network model also identifies the disulfide bond and the retinal ligand binding pocket to be the most rigid region in rhodopsin. Experiments confirm that 90% of the amino acids predicted by the FIRST method to be part of the core cause misfolding upon mutation. The observed high degree of conservation (78.9%) of this disulfide bond across all GPCR classes suggests that it is critical for the stability and function of GPCRs.network models ͉ membrane protein ͉ folding ͉ G protein-coupled receptor ͉ simulation R hodopsin is the only member of the G protein-coupled receptors (GPCRs), the largest family of cell-surface receptors, whose 3D structure is known (1). The signature motif of the GPCR family is a bundle of seven-transmembrane (TM) helices connected by polypeptide loops that form the cytoplasmic (CP) and the extracellular (EC) domains on opposite sides of the TM domain (Fig. 1). GPCRs perform extremely diverse and vital functions that include responses to light, odor, taste, neurotransmitters, hormones, and a variety of other signals (2). Whereas, in rhodopsin and related visual pigments, the ligand 11-cisretinal (RET) is covalently bound to the apoproteins (opsins), all other GPCRs occur in the ligand-free form, and subsequent binding of appropriate ligand(s) results in their activation. There is a wide variation in the nature of the ligands and their binding modes such as direct binding to the TM domain, the EC domain, or both.Based on pharmacological specificity and sequence conservation, GPCRs are divided into eight classes (3). Although there is no sequence homology between GPCRs in different classes, the seven-TM helix motif is conserved throughout, and all GPCRs share a common topology. They can be grouped into three main classes: receptors related to rhodopsins (class A), secretin receptors (class B), and the metabotropic neurotransmitter receptors (class C). Of these, class A, the largest class, contains Ͼ1,200 distinct members and Ͼ7,000 putat...