Molecular Vision 2014; 20:183-199
Received 27 June 2013 | Accepted 05 February 2014 | Published 07 February 2014
1Department of Molecular Ophthalmology, Lions Eye Institute, Perth, Western Australia; 2Centre for Experimental Immunology, Lions Eye Institute, Perth, Western Australia; 3Centre for Ophthalmology and Visual Science, University of Western Australia, Perth, Western Australia; 4EMBL/CRG Systems Biology Research Unit, Centre for Genomic Regulation (CRG), Dr. Aiguader 88, 08003 Barcelona, Spain; 5Universitat Pompeu Fabra (UPF), 08003 Barcelona, Spain
Correspondence to: Elizabeth Rakoczy, Lions Eye Institute, 6 Verdun St, Nedlands, WA 6009, Australia: Phone: +618 9381 0726; FAX:+618 9381 0739, email: firstname.lastname@example.org
Purpose: Following a previous study that demonstrated a correlation between rhodopsin stability and the severity of retinitis pigmentosa (RP), we investigated whether predictions of severity can be improved with a regional analysis of this correlation. The association between changes to the stability of the protein and the relative amount of rhodopsin reaching the plasma membrane was assessed.
Methods: Crystallography-based estimations of mutant rhodopsin stability were compared with descriptions in the scientific literature of the visual function of mutation carriers to determine the extent of associations between rhodopsin stability and clinical phenotype. To test the findings of this analysis, three residues of a green fluorescent protein (GFP) tagged rhodopsin plasmid were targeted with site-directed random mutagenesis to generate mutant variants with a range of stability changes. These plasmids were transfected into HEK-293 cells, and then flow cytometry was used to measure rhodopsin on the cells’ plasma membrane. The GFP signal was used to measure the ratio between this membrane-bound rhodopsin and total cellular rhodopsin. FoldX stability predictions were then compared with the surface staining data and clinical data from the database to characterize the relationship between rhodopsin stability, the severity of RP, and the expression of rhodopsin at the cell surface.
Results: There was a strong linear correlation between the scale of the destabilization of mutant variants and the severity of retinal disease. A correlation was also seen in vitro between stability and the amount of rhodopsin at the plasma membrane. Rhodopsin is drastically reduced on the surface of cells transfected with variants that differ in their inherent stability from the wild-type by more than 2 kcal/mol. Below this threshold, surface levels are closer to those of the wild-type.
Conclusions: There is a correlation between the stability of rhodopsin mutations and disease severity and levels of membrane-bound rhodopsin. Measuring membrane-bound rhodopsin with flow cytometry could improve prognoses for poorly characterized mutations and could provide a platform for measuring the effectiveness of treatments.
The majority of the more than 130 rhodopsin point mutations that have been discovered to date in humans  cause autosomal dominant retinitis pigmentosa (adRP). This degenerative disease of the retina usually begins with night blindness and progresses to peripheral vision loss and in the most severe cases to loss of central vision and total blindness. There is a large variability in the severity of rhodopsin-associated retinitis pigmentosa (RP). Some individuals experience night blindness, poor visual acuity, and reduced visual field soon after birth and can be completely blind by the time they reach their fourth decade of life  whereas others retain useful vision into at least their ninth decade . Severity is heavily influenced by the specific identity of the mutant allele, although environmental, dietary, and genetic cofactors are also likely to contribute .
After rhodopsin is synthesized in the rough endoplasmic reticulum, it is transported to the plasma membrane of the inner segment before traversing the cilia and forming endocytosed vesicles at the base of the outer segments that develop into the nascent discs [5,6]. The extracellular portion of the membrane-bound protein therefore becomes the intradiscal portion in rod discs. Cell culture lines such as HEK-293 and COS-7 have been used intensively to model the mutations that affect the preciliary transport and interactions with the chromophore 11-cis retinal [7-10]. In these models, the rhodopsin that is transported to the plasma membrane does not form discs, and thus, the intradiscal portion of the molecule can also be referred to as the extracellular region in this context.
Shortly after the first rhodopsin mutation was discovered , work by Sung and Nathans et al. [9,10] and the Khorana laboratory  revealed a clear division between two groups of rhodopsin mutations. Sung classified class I mutations as those that resemble the wild-type in yield, efficiency of regeneration with 11-cis-retinal, and in subcellular localization, but do not reach the outer segment [13,14]. Class II mutations in contrast, accumulate at low levels, regenerate poorly with 11-cis-retinal, and are inefficiently transported to the plasma membrane [8,12,15,16]. Mendes later collated evidence from several other laboratories to further extend the biochemical classification  to differentiate mutations that affect endocytosis (class III) , that affect rod-opsin stability and post-translational modification  (class IV), that show an increased rate of transducin activation  (class V), and that induce constitutive opsin activation (class VI) .
The most common reason why class II mutants fail to integrate with the plasma membrane is because they are retained at the endoplasmic reticulum . This organelle is essential for the correct folding of proteins and plays a critical role in numerous other aspects of post-translational processing. The physical association between proteinaceous aggregate and molecular chaperones in cells transfected with class II mutants  is convincing evidence that these mutant rhodopsins are retained by misfolding. A characterization of the way that alterations to rhodopsin stability can lead to aberrant folding is thus critical for a deeper understanding of the mechanisms behind this process. Several recent technologies have the potential to facilitate this.
One such tool is the algorithm FoldX , which was used for a systematic meta-analysis of rhodopsin mutant stability that highlighted a correlation between stability change and disease symptom severity . FoldX has been experimentally validated with 1,088 mutations in 20 proteins within most of the structural environments that occur in proteins  and is currently available as FoldX v3.0 beta3. Previous stability predictions have shown a good correlation with experimental data , and correlations have been identified between stability estimations and the severity of several conditions [24,27,28], including other visual diseases [24,29-31].
The presence of insoluble aggregates in some class II mutations has led to comparisons with other dominant neurodegenerative conditions such as Parkinson’s and Alzheimer’s disease . These misfolding diseases are thought to cause disease by toxic gain of function effects, and there are many similarities between these diseases and rhodopsin-associated RP, including the formation of ubiquitinated aggresomes, the recruitment of cellular chaperones, and the accrual of insoluble aggregates . However, evidence from mice [15,33] and cell cultures  suggests that dominant-negative effects can also result from the destabilizing class II rhodopsin mutations . One possible explanation is that in these cases, the misfolded protein retains an active dimerization motif so the wild-type/mutant heterodimers lose their functionality. Because dimerization in heterozygous individuals can incorporate such heterodimers as well as homodimers, the level of functional wild-type dimer may be reduced below the level that would be seen in the phenotypically normal heterozygote nulls and below the required threshold for normal gene function . The known links between destabilization and misfolding  and between misfolding and cellular localization  led us to investigate whether there is a direct correlation between surface rhodopsin and stability that could contribute to the severity correlation observed for the misfolding mutations that are categorized for rhodopsin as class II mutations.
A more detailed examination of the correlations between stability, structure, and clinical phenotypes in a subgroup of class II mutations is expected to provide unique insight into the traits associated with these mutations. This is an exciting prospect given the potential for improvements in disease prognosis and the possibility of building a basic framework on which to assess amenability to rhodopsin treatments. Our aim in this investigation was to investigate the correlation between the stability of class II rhodopsin mutations and four parameters that can be used to determine ocular health in patients with RP: visual acuity, visual field, night blindness onset, and daytime vision loss onset.
Clinical data was obtained by data-mining the available literature as described previously . Stability calculations were based on the FoldX estimations of the free energy of unfolding (ΔG) of the rhodopsin protein. ΔΔG is the difference between the stability of a mutant rhodopsin variant and the wild-type. Rhodopsin calculations of ΔG were made with the bovine form of the molecule . Bovine rhodopsin differs from human rhodopsin by only three positions within the amino acid sequence: Met216 in human rhodopsin is a leucine in the bovine sequence, human Ser270 is replaced by a glycine, and human Ser297 is replaced by a threonine. Each of these inter-species variations is predicted to have negligible effects on the molecule by FoldX and Sorting Intolerant From Tolerant (SIFT) , and the bovine rhodopsin structure has been used by many authors as a model of mammalian rhodopsin with immediate relevance to human disease owing to its high level of functional and sequence homology [38-40]. Thus, we believed the stability predictions derived from modeling mutations of the bovine structures were likely to closely represent mutations of human rhodopsin. FoldX data for the human disease mutation meta-analysis has been published previously , but the regression analysis in that case was performed on a functional basis. Here, we first grouped the mutations into their respective regions (intradiscal, cytoplasmic, transmembrane helix facing, transmembrane lipid facing, and retinal binding zone) before we analyzed the correlations with phenotype severity.
The other parameters used for our meta-analysis have been outlined previously , but for ease of interpretation, a summary of the inclusion criteria and additional details to clarify the methods are presented here. The age at which visual acuity or visual field was first reported to deteriorate was combined into the single parameter of “onset of daytime vision loss.” In cases where data for the age of visual field or visual acuity loss were not available, daytime vision loss has also been inferred from what were described in the source papers as the onset age of symptoms “other than night blindness,” or when electroretinography readings showed that cones as well as rods were severely reduced. Patients who could perceive only counting fingers, hand motion, or light were all given a visual acuity score of 0. The majority of sources clarified that visual acuity was assessed with best-corrected vision; whereas no sources specified that visual acuity testing was without correction. We therefore included all visual acuity measurements, assuming that they were taken with best corrected vision. Goldmann perimetry tests with isopters of V/4 were used to assess visual fields. Visual field data were either quoted in the original sources as the angle of constriction, the extent of the visual field loss in degrees squared ((o)2) , percentage of the normal visual field [42,43], or by perimetry charts. Visual field references that used “% of normal” were converted to visual field measurements using baselines derived from a previous study of normal peripheral isopter positions . Perimetry chart data were converted to visual field measurements with the software Image J .
Regression analysis (using Excel Analysis ToolPak, Excel 2010; Microsoft, Redmond, WA) of RP clinical severity and the relative stabilities of different rhodopsin mutations was performed on a regional basis, with the molecule separated into five structural regions: intradiscal, transmembrane helix facing, transmembrane lipid facing, retinal binding zone, and cytoplasmic. This was repeated for visual acuity and visual field, but these parameters are known to be influenced by age; therefore, the data were controlled for age by the inclusion of an “age” parameter in a multiple linear regression analysis.
To obtain an objective measure of surface rhodopsin, a plasmid that expresses human rhodopsin fused to a monomeric green fluorescent protein (GFP) was used. GFP tags at the C-terminus do not appear to have any effect on the intracellular transport of rhodopsin within non-ciliated, cultured cells [7,8]; but previous reports have shown that fusion proteins containing rhodopsin and enhanced green fluorescent protein are not trafficked correctly in rod cells [46-49]. We speculated that such aberrant trafficking may be exacerbated by the tendency for GFP to dimerize, and since there are several palmitoylation sites near the C-terminus, we wanted to avoid interference with such critical processes so replaced the dimerizing GFP with monomeric GFP. The Rhodopsin ORF was sub-cloned by Origene from plasmid RC211328 into the vector pCMV6-AC-mGFP to produce the custom-built vector RHO-AC-mGFP (Figure 1).
Site-directed random mutagenesis was performed with random bases at the codons P23, G106, and D190, using a method based on the mutagenesis approach described by Sanchis et al. . The thermal conditions for site-directed random mutagenesis PCR were as follows: 98 °C for 2 min; and then 5 cycles of 98 °C for 10 s, 56 °C−67 °C for 30 s, 72 °C for 1.75 min. This was followed by 20 cycles of 98 °C for 10 s and 72 °C for 30 s or 72 °C for 2.5 min and then a 5 min extension at 72 °C. The 50 µl reaction volume contained 0.2 mM dNTPs, 1 μl F and R primers (5 pmol), 50 ng template, 0.5 µl Phusion DNA Polymerase, Thermo Scientific (Waltham, MA), 1X of the reaction buffer that was supplied with the polymerase. The forward primer sequences were 5ʹ-GTG GAA TCN NKT ACT ACA CGC TCA-3ʹ for D190 mutagenesis; 5ʹ-GCA GCN NKT TCG AGT ACC CAC-3ʹ for P23 mutagenesis; and 5ʹ-GAT ACT TCG TCT TCN NKC CCA C-3ʹ for G106 mutagenesis. The reverse primer for all three reactions was 5ʹ-ACG TCG TGA CTG GGA AAA CC-3ʹ. The three forward primers contained the NNK degenerate base sequence, which randomized the resulting product at the loci of interest. Potential mutant rhodopsin plasmids were sequenced to determine the identity of the mutations.
FoldX data for human mutation analysis had been published previously , but further analysis of the three target residues that underwent site-directed mutagenesis was required. One of the main barriers to accurately predicting protein stability are errors within the structural template. Therefore, one should always model mutations of interest on several template structures. In our earlier analysis, we used five template structures, which mostly resulted in similar energy predictions. We predicted the effects of all possible mutations at the three positions targeted for site-directed mutagenesis, using these original five template structures as well as the respective second chains from four of the crystallographic units. We thus obtained the average and standard deviations (after removing two outliers) from nine PDB structures.
HEK-293 cells at a passage number of between 38 and 50 were seeded at 300 cells/mm2 in DMEM (Gibco® - Life Technologies, Carlsbad, CA) media into wells of a 24-well plate. Several earlier studies used HEK-293 cells for rhodopsin expression studies, and evidence suggests that the majority of rhodopsin has a native conformation in these cells [10,51]. After the cells were incubated for 48 h at 37 °C in 5% CO2, the RHO-AC-mGFP plasmid was transfected into the cells using X-treme gene HP transfection reagent (Roche, Basel, Switzerland) at a ratio of 1 µl transfection reagent to 1 µg of RHO-AC-mGFP DNA in 100 µl of serum-free DMEM media. 50 µl of this transfection mixture was transferred to each well. Untagged rhodopsin controls were transfected using the same conditions as the tagged constructs, but with the rhodopsin plasmid RHO-SC300090 (Origene, Rockville, MD).
Forty hours after transfection, cells were washed once in PBS (136.9 mM NaCl, 2.7 mM KCl, 8.1 mM Na2HPO4, 1.5 mM KH2PO4). They were then disassociated from the vessel surface by covering with 50 µl trypsin and incubated for 4 min at 37 °C. Cells were then resuspended with 150 µl DMEM media and centrifuged at 400 ×g for 4 min, washed by resuspension and centrifugation in PBS, and then resuspended and fixed in 4% paraformaldehyde in PBS for 20 min. After two more PBS washes, a blocking solution containing 3% bovine serum albumin in PBS was added to the cells before anti-rhodopsin RET-P1 monoclonal antibody was added (Sigma-Aldrich, St Louis, MO) at a dilution of 1:2,000 and incubated on ice for 1 h. The RET-P1 antibody recognizes amino acids 4–10 at rhodopsin’s N-terminus and is therefore detected under non-permeabilizing conditions in which no detergents have been added at any stage of sample preparation. The cells were washed twice before APC-Cy7 conjugated goat anti-mouse antibody was added at a dilution of 1:2,000. Further incubation at 30 min on ice was followed by two washes and resuspension in 200 µl PBS for analysis. The Ret-P1 gate was set according to the wavelength at which there were negligible numbers of cells in the secondary antibody control.
The culture conditions and time frames for immunohistochemistry were identical to those used for flow cytometry except that the cells were grown on 0.7 cm2 eight-chamber slides (Becton Dickinson, Franklin Lakes, NJ) and were transfected with 30 µl transfection mix. Washes were performed by adding 500 µl PBS to each chamber, followed by 5 min at room temperature. The secondary antibody for microscopy was donkey anti-mouse conjugated with Alexa Fluor 568 (Molecular Probes, Eigene, OR). Image analysis of fluorescence microscopy photographs was performed using ImageJ  using the following procedure: Fluorescence intensity thresholds for the GFP signal were selected using wild-type images, to limit the regions of interest to cells that were actively producing the RHO-GFP fusion product. The selected areas for each GFP image were then applied to the AF568 antibody image from the equivalent view field. Background grayscale measurements were subtracted from the grayscale measurements of the GFP+ cells, before the red/green ratio was determined.
Confocal microscopy was performed using a Nikon Ti-E inverted microscope with a Nikon A1Si spectral detector. Cells analyzed with confocal microscopy were labeled with biotinylated antimouse CD44 (Biolegend, San Diego, CA) and streptavidin conjugated to Alexa Fluor 594.
Regression analysis of the human clinical data was performed using the Microsoft Excel data analysis toolkit. ANOVA and the Tukey-Kramer post-hoc test were performed using GB-Stat (v9.0).
Our original analysis of the clinical phenotype data identified a correlation between the means of the daytime vision loss onset and ΔΔG when all non-functional residues were analyzed . The new analysis here, in contrast, takes all available data into account individually rather than only the mutational means, and highlights several additional correlations in addition to identifying the region with the strongest correlation. Carriers of mutations within the intradiscal region of rhodopsin (Table 1) retain a highly significant correlation under the new analytical parameters (see Figure 2, Table 2). For every ΔΔG increase of 1 kcal/mol, there was a decrease in the age of onset of night blindness of 0.9 years (p<0.0001, R2=0.3), a decrease in the onset of daytime vision loss of 0.5 years (p=0.008, R2=0.1); a decrease in visual acuity of 0.009 (p=0.007, R2=0.3; Figure 2, Table 3), and a decrease in visual field of 180(o)2 (p=0.003, R2=0.2; Figure 2, Table 4). In contrast to the intradiscal residues, the transmembrane, cytoplasmic, and retinal-zone residues showed no correlation between stability change and visual function.
Although highly significant, the intradiscal residue correlations were associated with low R square values by the standards of other correlative studies due to variability between individuals with the same mutation. These data suggest that there is a genuine correlation between energy change and disease severity, but the increased consistency and strength of correlation in the intradiscal region compared to the other regions should be emphasized.
As expected, the multiple linear regression analysis showed that there was also a strong association between increasing age and decreasing visual acuity (p=5.5×10−7), but the association between age and visual field deterioration was marginal (p=0.05 for the V/4 isopter data).
To determine whether the correlation between rhodopsin stability and RP severity was due to reduced efficiency of its transport to the plasma membrane, we aimed to quantify the levels of surface rhodopsin. Three residues (G106, D190, and P23) were selected for mutagenesis because of their position on the intradiscal region of the protein (Figure 3) and because they are associated with multiple pathogenic mutations in the human population. A series of mutations were generated and the relevant ΔG values calculated (Figure 4).
Cell dissociation methods were optimized for maximum antibody signal intensity. Citric saline was compared with EDTA and trypsin. Previous reports have shown that any one of these substances can be used to dissociate cell monolayers for flow cytometry . Our tests showed that the strongest APC-Cy7 signal was observed in the trypsinized cells.
Population density gradients showed that within cells transfected with wild-type rhodopsin, the intensity of GFP expression correlated closely with the intensity of APC-Cy7 detection (Figure 5A). In mutant variants though, as well as an overall reduction in APC-Cy7 signal intensity, there is a wider range of APC-Cy7 intensities in GFP+ cells (Figure 5B).
Mutations that resulted in a high ΔΔG were associated with a noticeable drop in levels of detectable surface rhodopsin, whereas variants with low ΔΔG had a milder effect (Figure 6A and Table 5). ANOVA of the RET-1/APC-cy7 signal from the GFP+ transfected cells, followed by a Tukey-Kramer post-hoc test showed a significant elevation (p<0.01) of the wild-type (1320.8RFI ±339.9) compared to any of the mutant variants (mutant means ranged from 337.4 RFI ±90.6 SD for G106R to 70.9 RFI ±10.2 SD for P23N). Reanalysis following the normalization of these means to the mutant/wild-type ratio revealed that G106R was significantly elevated (p<0.01) above all of the other mutant variants (G106R=0.259±0.05 SD) and D190N was significantly higher than all but one (D190P; p<0.05) of the other mutations (D190n=0.153±0.02 SD; D190p=0.106 ±0.04 SD; others=0.09±0.02 SD). The difference in the APC-Cy7 signal between D190P and the other more destabilized variants was not significant (Figure 6A,C).
There were some concerns that the drop in surface rhodopsin at high levels of destabilization may have been linked to the trypsinization process. However, a two-way ANOVA comparison between microscopy-based image analysis and flow cytometry showed that although the two methods resulted in significantly different ratios of surface antibody to GFP signal (p<0.01), the interaction effect between the data sets was not significant (p=0.9; see also Figure 6). This suggests that although the different antibodies and detection technologies caused an output shift, the size of the shift remained consistent for each mutant residue and so trypsinization did not affect these data.
Flow cytometry showed that there was a significant weakening of the overall GFP signal in all of the mutant variants compared to the wild type, but that there was little difference in the GFP intensity between the mutants (Figure 6). This would suggest that overall rhodopsin levels were reduced by even the mildest of changes to stability, but that further increases in ΔΔG between 1.94 and 27.44 kcal/mol did not have a compounding effect.
We followed the same strategy for statistical analysis of GFP intensity as we used to analyze APC-Cy7 data. ANOVA with repeated measures and Tukey-Kramer post-hoc analysis showed that the wild-type mean GFP intensity (2444.8 RFI±1025.4 SD, n=5) was significantly higher (p<0.01) than any of the mutants (mutant means ranged from about 886 RFI±453.1 SD for P23K to 1,549 RFI±427.4 SD for P23Y). The mean GFP intensity for all mutants was 1316.5 RFI±658.5 SD). A subsequent comparison of the normalized mutant/wild-type ratios (Figure 6D) revealed a marginally significant difference (p<0.05) between P23K (0.37 ±0.10) and both D190N (0.64 ±0.18 SD) and P23Y (0.67 ±0.13 SD), but there were no significant differences in the GFP signal intensity between any of the other mutants. The differences between the GFP intensities for the different variants were therefore marginal and not as pronounced as they were between the equivalent APC-Cy7 intensities.
Cells transfected with the G106R plasmid, the mutation with the lowest ΔΔG at 1.94 kcal/mol, therefore retained the highest level of membrane-bound rhodopsin of any of the mutant variants. The D190N variant has the second lowest ΔΔG of any mutant variant (2.26 kcal/mol) and flow cytometry demonstrated a drastic drop in surface expression compared to G106R, but higher levels of surface rhodopsin than all of the other mutants.
Microscopy reflected the findings of flow cytometry to show that the least destabilizing mutations exhibited the highest levels of membranous rhodopsin. Surface rhodopsin was not evenly distributed throughout the plasma membrane of HEK-293 cells when transfected with the majority of mutant variants, but was expressed as small intense foci. In contrast, the wild-type, G106R, and D190N transfections demonstrated a more even distribution of surface expression, although D190N surface expression was present at greatly reduced levels. All mutants exhibited some expression of rhodopsin at the surface of the plasma membrane, but every mutant with the exception of G106R was associated with these small foci, even D190N, which also contained the more diffuse membranous protein in some cells.
A stark contrast between the wild-type, G106R, and other mutations was observed with confocal microscopy (Figure 7 and Figure 8). Cells transfected with the wild-type, the two human mutations G106R, D190N, and the most destabilizing mutation (P23Y), were labeled with CD44 to determine the extent of colocalization between the GFP and this plasma membrane marker. As expected, the GFP and CD44 colocalized at high levels in the membranes of cells transfected with wild-type and G106R plasmids. Confocal microscopy also provided a clearer understanding of RHO-GFP localization in cells transfected with the D190N variant. All transfected cells showed high levels of cytoplasmic GFP, but some cells also showed even but low-level colocalization between CD44 and GFP at the plasma membrane.
The observed relationship between destabilization and the clinical data demonstrates the delicate balance of stability that governs the correct function of rhodopsin and its influence on the severity of retinal degeneration. This relationship highlights that in intradiscal residues, the effect of mutation-induced changes to rhodopsin stability has such a powerful impact that other structural/functional severity modifiers (such as interference with glycosylation or interaction motif activity) are masked by the stability change effects.
The high stability-severity correlation in the intradiscal region can be explained by referring to the work of Briscoe et al. , who found that sites containing disease-associated rhodopsin mutations with a lower expected chemical severity (ECS) are on average more critical to the protein’s correct function than those with a high ECS. The intradiscal region has the lowest functionality of all of the regional sub-categories. If a residue has a critical functional role, then factors other than destabilization will become the most important causes of molecular dysfunction. A minor change in structural stability could thus have a drastic effect on a functionally critical residue, and the relationship between ECS or stability and disease severity would be of less relevance to the disease status than the affected function. The respective functions of retinal-binding in the transmembranous region and the cilia trafficking, outer segment transport, and transmolecular interactions in the cytoplasmic region are thus all inhibited by less severe residue changes than are the intradiscal functions. It follows that there is less likely to be a correlation between stability and function in regions containing numerous residues with critical functions.
A hierarchy of functional importance may not provide the only explanation for the more robust correlation in the intradiscal residues, however. The correlation between disease severity and the destabilization level of intradiscal residues in this analysis was much more statistically robust than the correlation observed in the earlier function-excluded data set . One possible explanation for this could be that the small size of the glycosylation sites in the intradiscal region relative to the protein-interaction sites of the other regions means that mutations close to the specific active residues would have a reduced conformational impact on these functions. It is also possible that the intradiscal domain plays a particularly important structural role for the molecule compared to the cytoplasmic domains: The four intradiscal domains are significantly associated with each other , and the central E2 intradiscal domain (see Figure 2) is particularly sensitive to the effects of mutation with 12 of the 26 residues mutated in RP . Whatever the specific cause though, the importance of molecular destabilization relative to residue-specific interaction activity is increased in the intradiscal residues.
Another study has also generated multiple mutations in three residues of a single domain . The transmembrane residues G121, A117, and G114 were each mutagenized seven times to generate a cohort of 21 mutations. These residues are located on helix 3, the sub-domain that contains the counterion to the protonated Schiff’s base at Lys296 that is formed by the binding of 11-cis retinal  and are thus within the retinal binding zone. This study showed a correlation between mutant residue size and the degree of blue-shift in the λmax (wavelength of maximum absorbance) of the pigment. These findings are in agreement with our previous research that showed that the proximity of residues to retinal binding increases the severity of mutations .
Although there have been many previous reports of mutations that cause the internalization of rhodopsin within various cellular models [8-10,32], this is the first time that a quantitative method has shown a direct correlation between rhodopsin stability and its presence at the surface. The correlation is not surprising given that misfolded proteins are retained in the endoplasmic reticulum. However, the discovery of a quantifiable relationship between these two factors suggests that this approach might enable improvements to predictions of mutational effects. Using a transient transfection method rather than a stable transfection means that the rhodopsin sequence will be primarily expressed directly from the plasmid and there will be some expression variability caused by copy number and low-level genome incorporation. However, we considered that the variability incurred by transient transfections would provide mean expression levels that are more consistent for comparisons of different variants than would the genomic incorporation at single sites resulting from stable transfections.
We expected that we might see a linear relationship between FoldX predictions of stability and the surface expression of rhodopsin. Instead, there was a threshold beyond which only minimal rhodopsin was transported toward the plasma membrane, giving rise to a sharp drop in total rhodopsin levels at low ΔΔG. This means that the linear relationship observed between phenotype and stability must be influenced by additional stability-induced effects. Rhodopsin continuously undergoes rapid changes in morphology during the visual cycle, and thus, the nature of rhodopsin stability is more complex than simply maintaining enough stability to avoid misfolding and proteolysis. Some mutant variants are synthesized and transported correctly, but demonstrate aberrant behavior following multiple cycles of retinal isomerization and photobleaching . Pathogenic mutations with moderate misfolding from the intradiscal region can be partially rescued by chaperones in vitro , but these rescued molecules tended to result in an increased level of metarhodopsin-I in relation to metarhodopsin-II, and an accumulation of aberrant photoproduct when subjected to prolonged illumination . The relative stabilities of the different photointermediates are beyond the scope of this paper, but importantly, our method describes aspects of the core biosynthesis stability that allows the molecule to form at the ribosome and subsequently to maintain its conformation, thus allowing its release and transportation from the endoplasmic reticulum to the plasma membrane.
Our approach of integrating a literature-derived database with crystallography-based predictions of molecular stability and with flow cytometry has yielded data that support the view that the degree of rhodopsin destabilization has a powerful influence on disease severity. It has also highlighted a potential mechanism for this phenomenon and has shown that mutations in the intradiscal region are more likely to be amenable to severity predictions based on stability than are mutations in the other regions. There is a large variability in the severity of rhodopsin mutations, even between individuals who harbor the same mutation, and the intradiscal region is no exception. Indeed, the first mutation to be identified in humans (P23H) lies within the intradiscal region and is notoriously variable. Prognostic devices for rhodopsin-associated RP therefore cannot depend exclusively on genotype-based predictions and further characterization of the other contributing elements in this multifactorial disease are critical if reliable and accurate prognoses are to ever be obtained. However, it is important to refine the interpretation of the genetic contribution as far as possible within these confines. Statistical means of severity based on genotype will certainly improve the current methods for estimating the severity of rare mutations. The estimated structural impact of any mutation has been shown to have potential as a valuable tool for practitioners and could be used for risk/benefit analysis of future treatment options.
RM was funded by the Lions Eye Institute department of Molecular Ophthalmology. The authors thank Lori Mull from Origene for valuable advice relating to the design of the RHO-AC-mGFP plasmid. We also thank Sarah Bolt and Julie Marsh from the University of Western Australia Centre for Genetic Epidemiology and Biostatistics for assisting with the analysis of associations between RP severity and FoldX stability estimations. The sequencing team at the Australian Genome Research Facility (Perth node) sequenced the mutant plasmids. Alysia Buckley of the Centre for microscopy, characterization and analysis assisted with confocal microscopy of the transfected cells. We are very grateful to Luis Serrano who provided invaluable training and support for the FoldX analysis.