Genetic mutations in X-linked retinoschisis and their role in protein structure/function. ahttp://www.retina-international.org/sci-news/xlrsmut.htm.
We describe a promising approach using
2. Biological factors and electroretinography of XLRS
2.1. Biology of retinoschisin protein: secretion and expression, hypothesized function, underlying mechanisms of functional loss of RS1 protein
Retinoschisin (RS1) is a retinal secreted photoreceptor disulphide-linked oligomeric protein. A number of mutations in RS1 co-segregate with XLRS providing strong evidence that the disease is caused by mutations in the
2.2. Retinal morphology in retinoschisis: photoreceptors disruption and altered photoreceptor/bipolar cells synapse in the RS1-KO mouse
Currently, there are three mouse models of human XLRS. One has been generated by homologus DNA insertion across the endogenous
2.3. Clinical features of XLRS disease: visual acuity, visual field, OCT, fundus appearance, ERG
XLRS classical phenotype is present in two-thirds or more of affected males (George et al., 1996) and consists in an early-onset central vision loss from bilateral foveo-macular cystic cavities involving the inner retina and additional retinal layers (Gregori et al., 2009) that extend from the fovea in a spoke-wheel pattern for approximately 1-1.5 disc diameter (Figure 2).
XLRS patients often have hyperopic refractive error with astigmatism (Roesch et al., 1998; Kato et al., 2001; Vijayasarathy et al., 2009). Visual acuity ranging between 20/20 to 20/200 in young patients, who come to ophthalmological attention because of reading difficulty, due to the macular involvement of the disease. Patients typically exhibit modest severity at a young age (cystic form), with a slow progression until the second decade, when most retain around 20/70 vision. It worsens through the teenage years and stabilizes in adulthood when the atrophic stage can lead to legal blindness (Lesch et al., 2008). However, some XLRS males exhibit unpredictable, severe disease even at very early age, and surprisingly, the ERG can show an unexpected preservation of the scotopic b-wave (Eksandh et al., 2005; Sieving et al., 1999a).
XLRS is a highly penetrant disease condition (Forsius et al., 1973). In some cases subtle retinal changes are difficult to identify by ophthalmoscopy and can be detected only by optical coherence tomography (OCT) scans (Eriksson et al., 2004) (Figure 2). In the late stages of the disease when macular atrophy starts, a hyperfluorescent window defect of the retinal pigment epithelium (RPE) on fluoresceine angiography (FA) can be observed. This feature, however, can mimic retinal abnormalities seen in age-related macular degeneration or Stargardt maculopathy. In such cases, the combination of ERG, FA, and OCT is needed to determine the correct diagnosis. The electronegative scotopic ERG is pathognomonic of XLRS; the FA will not exhibit the dark choroid sign that is typically present in Stargardt’s maculopathy; and the OCT can detect cysts otherwise undetectable by standard ophthalmoscopy. Peripheral retinoschisis can be found in half of affected males. Retinal splitting can involve all the retinal layers through the nerve fiber layer. The consequent abnormal visual signaling causes an dense scotoma in the involved areas (Li et al., 2007). Sometimes peripheral retinal scars, possibly to be considered pre-existing areas of self-limited schisis, can be observed.
2.4. Photoreceptor and bipolar cell function as seen in the full field ERG
The ERG is a useful tool for evaluate outer retinal function and may assist in identifying the primary sight of XLRS pathology. Through the interpretation of the ERG, one can pinpoint disease at the cellular level in outer retinal visual loss and monitor disease progression. The ERG is a good phenotypical measurement in wide spread disease due to its objective nature and ability to measure function of the entire retina. To better understand the pathology behind the characteristic electronegative ERG in XLRS, one must first understand normal retinal physiology. In a dark-adapted retina, free flowing ions move to and from the rod photoreceptor. When rhodopsin in the rod photoreceptor captures a photon of light, it undergoes a conformational change causing a cascade of intracellular events within the photoreceptor outer segment resulting in hyperpolarization of the cell. Rod hyperpolarization causes an increased ionic potential towards the distal retina that results in the a-wave (Breton et al., 1994; Hagins et al., 1970). The rod then signals the on-bipolar cell which depolarizes and causes an influx of ions into the proximal retina. This forms the ERG b-wave. In XLRS, the b-wave amplitude is smaller than the a-wave amplitude (i.e. does not return to pre-flash baseline) and is deemed ‘electronegative’ (Figure 2H). An electronegative ERG suggests dysfunction proximal to the rod outer segment.
2.5. Bipolar cell dysfunction and “electronegative” b-wave in both human and murine model of XLRS
In the retina of the RS1-deficient mouse and the XLRS human the reduced scotopic ERG b-wave amplitude (Figure 2H) suggests an effect on bipolar-cell function, perhaps initiated from an altered photoreceptor/bipolar-cell synapse activity. This implies primary involvement of the inner retina where schisis cavities occur in XLRS. The identification of L-type voltage-gated calcium channels (Shi et al., 2009), Na/K ATPase (alpha3, beta2 isoforms) (Friedrich et al.), the sterile alpha and TIR motif-containing protein SARM1 (Molday, 2007), and phospholipids/Ca2+ (Vijayasarathy et al., 2007) as potential RS1 interacting ligands corroborate a role for RS1 in cell signaling events as well. However, ERG show considerable variation in this disease and ophthalmologists have to be aware about the possibility of relative ERG b-wave preservation especially in presence of atypical or mild manifestation of the disease (Eksandh et al., 2005; Sieving et al., 1999a; Vijayasarathy et al., 2009). These observations may imply that the degree of ERG dysfunctions may not correlate with the extent of fundus involvement and that synaptic integrity might be impaired secondarily over a greater extent of the fundus. In addition, because disruption of the inner segment architecture due to the loss of RS1 might be the basic mechanism that underlies the displacement and disorganization of photoreceptors in RS1-deficient mouse retina, reduced ERG a-wave responses can be found (Kjellstrom et al., 2007b). Clinically, abnormal photoreceptor function has been described in some XLRS subjects (Bradshaw et al., 1999; Kellner et al., 1990) mostly related to the advanced stages of the disease (Peachey et al., 1987) when the outer retinal atrophy may occur.
2.6. The significance of phenotypic b/a-wave ratio
In retina electrophysiology signal transmission between rod photoreceptors and the post-synaptic bipolar cell neurons is measured by the amplitudes of the a-wave and b–wave, respectively (Newman and Odette, 1984; Robson and Frishman, 1995; Stockton and Slaughter, 1989). To account for any compromise of the photoreceptor layer in XLRS, a ratio of photoreceptor to bipolar-cell synaptic transmission compared to photo-transduction quality can be used: the b-wave amplitude to a-wave amplitude ratio (b/a-wave ratio). One can use the b/a-wave ratio to adjust for any photoreceptor dysfunction affect on the b-wave amplitude, and can imply the quantity of dysfunction in bipolar cells. Normally the b-wave amplitude is greater than the a-wave in the dark-adapted ERG for bright stimuli, but in XLRS the b-wave naturally is reduced (Peachey et al., 1987). Thus, a reduced b/a-wave ratio is characteristic of inner retinal disease such as XLRS.
2.7. Natural history and disease progression: role of age
Several reports based on visual acuity loss (Apushkin et al., 2005), fundus changes with schisis cavities modification (Lesch et al., 2008), and retinal function deterioration (Vijayasarathy et al., 2010; Vijayasarathy et al., 2009) have focused on the “progressive” nature of XLRS disease. In detail, visual acuity is mostly impaired starting in grade school (Apushkin et al., 2005) due to macular schisis. However, even if 20% of XLRS children (Roesch et al., 1998) present significant vision loss in early childhood, half of these patients may have poor vision due to complications of the disorder, such as retinal detachment or vitreous hemorrhage. A bimodal age is a characteristic of presentation of XLRS in childhood (George et al., 1996). One group presented strabismus and nystagmus at an age less than 2 years, and the other group presented only poor vision at approximately 7 years of age. Exotropia and esotropia were in the same proportion, and in the group with exotropia the occurrence of poor vision from severe schisis was highest. Rhegmatogenous (tractional retinal detachment) may also be the initial presentation of XLRS. Other XLRS males may present during infancy peripheral bullous elevation of the inner retinal layer accompanied with vitreous hemorrhage (George et al., 1996). In such cases it is difficult to assess a “progressive” trait of the disease.
Retinal appearance at funduscopy, OCT scans and ERG b-wave amplitude reduction (Figure 2, E-H) provide an estimate of the retinal changes in XLRS patients. A substantial change from the “cystic” (< 25 years old) to the “atrophic” (> 25 years old) pattern of the retina with age was recently described (Lesch et al., 2008). We also showed (Vijayasarathy et al., 2009) worsening of retinal function by age in a family harboring a null mutation in exon 5 of the
In addition, disease progression has also been observed in natural history studies of the
2.8. Clinical severity (metanalysis of severity grade systems) and genotype-phenotype correlation in human XLRS
Wide phenotypic heterogeneity is one of the remarkable characteristic of this disease, and it is possibly related to the increasing number of mutations identified in the
Investigators have explored whether the type of mutations on the
Moreover three Japanese XLRS patients were described with retinal detachment or severe macular lesion and low visual acuity at a very young age; two of them harbored deletions of exon 1 or of the boundary region between exon 3 and intron 3 (Shinoda et al., 2000). The authors also reported (Shinoda et al., 2004) a progressive XLRS phenotype in two siblings harboring a deletion of RS1 exon 4 and summarized that severe cases of XLRS are usually associated with upstream deletions (exons 1–3) in the
|Missense||"/110||Change a single nucleotide||Code for a different amino acid||Affects protein folding, sorting, loss of protein|
|Nonsense||9||Codes for a stop codon||Could truncate the protein||Full or partial loss-of-function|
|Splice site||16||A that or a number of nucleotides in the specific site at which of an takes place during the processing of into||The abolishment of the splicing site results in one or more introns remaining in||Production of aberrant|
|Deletions||26||Delete one or more nucleotides in DNA||In-frame deletions or frameshifts||Full or partial loss-of-function|
|Insertions||4||Insert one or more nucleotides into DNA||In-frame deletions or frameshifts||Full or partial loss-of-function|
|Silent||6||Change a single nucleotide||Code for the same amino acid||Native-like protein|
3. Molecular modeling, functional role and severity of pathogenic mutations in human XLRS
Genetic mutations often cause a change in protein expression. As described above several primary mechanisms might be responsible for the loss of RS1 function due to genetic mutations. The first mechanism is the misfolding of the discoidin domain, which changes the putative adhesive properties of the protein. Caused by missense mutations, RS1 misfolding might initiate formation of aberrant proteins in non-native conformation. The second mechanism is the defective disulfide-linked subunit assembly of RS1 into dimers and octamers which results in abnormal protein secretion caused by missense changes affecting reactive thiol groups. The third mechanism is the inability of RS1 to insert into endoplasmic reticulum (ER) membrane as part of the protein secretion process. This is caused by severe changes associated with the insertion or the removal of large fragments in protein molecule, protein truncations, and null mutations.
According to analysis in Paragraph 2.6 and Table 2 the majority of XLRS mutations could be divided in two large groups of (a) missense mutations and (b) non-missense changes with a severe XLRS phenotype. Although a severe role of non-missense changes is confirmed in numerous cases, the severity of XLRS phenotype for missense changes is less obvious. Here we imply molecular modeling of protein structure as described by ERGs clinical severities of XLRS via the computational evaluation of impact of pathogenic missense mutations.
3.1. Protein folding/misfolding in ER and protein folding landscapes
The effect of missense mutations on a protein fold is difficult to understand without some basics of the native protein folding. Over one-third of newly synthesized proteins are translocated into the ER, which ensures their delivery into the secretory pathway. The ER lumen contains molecular chaperones that maintain polypeptide solubility, enzymes that are necessary for post-translational modification of proteins, and factors that directly assist in the folding of newly synthesized polypeptides such as foldases and isomerases. Proteins that not attain their native conformations, because of genetic mutations, cellular stress or other random events, may harm the cell. Therefore, protein systems within the ER lumen mediate quality control, that resulting in the resolution of terminally misfolded proteins from correctly folded proteins and folding intermediates. To compensate for the potentially disastrous consequences of misfolded protein accumulation, ER-retained species are most commonly destroyed.
There are numerous observations of retinoschisins modified by missense mutations and expressed in cell culture which demonstrate the absence of modified RS1 bands from SDS-PAGE or native gels and associated with a loss of protein due to misfolding and rapid degradation (Wang et al., 2002; Wang et al., 2006; Wu and Molday, 2003; Wu et al., 2005). Several major mechanisms lead to degradation of misfolded protein. About 30% of cytosolic proteins carry specific KFERQ-targeting motifs in their sequences allowing for selective targeting of cytosolic proteins to the lysosome for degradation by chaperone-mediated autophagy (Kaushik and Cuervo, 2006, 2008). However, the absence of the specific targeting motif in the RS1 sequence suggests another option. As a result, misfolded or unassembled proteins are retained in the ER bound to chaperones or lectins until they are delivered to the cytosol for degradation in the ubiquitin-proteasome pathway (Goldberg, 2003). Mutant variants with significant changes in protein fold caused by missense mutations or deletions could generate misfolded proteins susceptible to hydrolysis similar to that of human haemoglobin in which about one-fifth of missense proteins undergo rapid degradation (Goldberg, 2003). Therefore, the severe mutant variants are expected to have a markedly perturbed non-native protein fold causing misfolding in the ER and later digestion in the proteasome (Vijayasarathy et al., 2010).
The mechanism by which a polypeptide chain folds to a specific three-dimensional protein structure was recently discussed (Dobson, 2004). In the ER, proteins fold into their native conformations and undergo different post-translational modifications and the formation of disulphide bonds. Native states of proteins almost always correspond to the structures that are most thermodynamically stable under physiological conditions. The primary sequence of protein is the major determinant for proper folding of protein (Anfinsen et al., 1961). The protein folds along several competing pathways into intermediate, non-native structures with decreasing free energies until it achieves a conformation with the lowest energy to form a protein with native interactions (Dobson, 2004). The total number of possible conformations of any polypeptide chain is so large that a systematic search for this particular structure would be impractical. The folding process involves a random search of many protein folds accessible to a polypeptide chain.
The natural fluctuations in the conformation of an unfolded or incompletely folded polypeptide chain cause highly separated residues in the amino-acid sequence to come into contact with one other. In addition, if the energy surface of protein or ‘landscape’ has the correct shape, only a small number of all possible conformers needs to be tested by any given protein molecule during its transition from a random coil to a native structure. The landscape is encoded by the amino-acid sequence of protein. Therefore, natural selection has enabled proteins to change so that polypeptides are able to fold rapidly and efficiently. Such a description is often referred to as the ‘new view’ of protein folding (Dobson, 1999, 2003, 2004). On average, native-like interactions between residues are more determined and the polypeptide chain is finding lowest-energy structure by a process of trial and error. Misfolded or intermediate non-native proteins never achieve this lowest energy minimum. Pathogenic mutations produce a severe or mild effect in a protein structure and cause the protein to stay in an intermediate or even misfolded conformation with a relatively higher free energy.
3.2. Role of cysteine residues and disulphide bonds in protein stability
Cysteine residues normally form disulphide bonds which significantly constrain the number of possible conformers from finding a lowest-energy structure. The disulphide bond or disulphide bridge is formed by the pairing of two thiol groups of cysteines as a result of protein oxidation. A disulfide bond is a strong covalent bond of about 2.05 Å in length connecting two sulfur
The disulfide bond holds two fragments of the protein chain together. This stabilizes the native fold of a protein by lowering entropy of the polypeptide chain and by condensing hydrophobic residues around the bond into local hydrophobic core using hydrophobic interactions. The disulfide bond increases the effective local concentration of hydrophobic residues causing a decrease of the local concentration of water molecules in the hydrophobic core. Since water molecules break up hydrogen bonds, their low concentration stabilizes a protein secondary structure.
The native disulphide bonds are critical for proper protein folding. The native structure of a protein typically has a single disulfide bind. In proteins with more than two cysteines, non-native disulfide species may be created. Non-disulphide species are almost always unfolded. As the number of cysteines increases, the number of non-native species increases significantly. Indeed, RS1 monomer has 10 cysteines which can form up to five disulphide bonds. In general, 10 cysteine residues can form
In the cell in order to avoid non-native disulphides formation, protein disulphide isomerases (PDI) are necessary to catalyze the inter-conversion of disulfide species and accelerate the formation of the native disulfide species. In the ER, PDI has a function of a cellular chaperone with the essential foldase activity to maintain a native protein fold (Noiva, 1999; Tsai et al., 2001). PDI catalyzes oxidation and reshuffling (isomerization) of disulphides in substrate proteins by using the catalytic CxxC motif (Gruber et al., 2006). Current data suggest that PDI docks close to two
Hence, we suggest that genetic mutations which are associated with the cysteine insertion/removal might interfere with suitable PDI binding in the vicinity of missense change. This might cause the appearance of new non-native disulphide bonds previously corrected by the foldase action. The insertion or deletion of cysteine residues should be taken into account when the severity of damage to protein molecule caused by the pathogenic missense change is evaluated.
3.3. Discoidin domains and molecular modeling of RS1 protein structure
Currently more than 72,000 protein structures have been determined by protein crystallography and deposited in the RCSB protein data bank. This list could be significantly expanded if structural genomic approaches such as a sequence homology and molecular modeling are used to generate protein structures. These computational methods were used successfully to predict the effect of thermodynamics experiments, to simulate the dynamics of protein folding and unfolding, and to explain a functional and structural role of elements of protein structure and other areas.
Similar methods could be used to evaluate an impact of missense mutations on protein structure stability and function. We used molecular modeling to predict a structure of mature retinoschisin which is formed by Rs1 and discoidin domains (Sergeev et al., 2010). Discoidin domains have a similar fold with elements of spatial structure in a shape of a β-sandwich, shown schematically in Figure 3 (Baumgartner et al., 1998) as the superposition of three structures demonstrating a binding to different molecules. The galactose, phospholipids, and collagen binding sites are located in shown in the lower part of this figure (area of spikes) suggesting a high conformational flexibility in this area (Carafoli et al., 2009; Macedo-Ribeiro et al., 1999). Moreover, the discoidin domains are able to form a dimeric structure shown by the cyan ribbon (Macedo-Ribeiro et al., 1999). Therefore binding to galactose, membrane, collagen, and a dimer formation could be potentially expected for the discoidin domain of retinoschisin. In general, the discoidin fold is conserved in more than 9 different proteins. The discoidin domain is present in single or multiple copies in extracellular or trans-membrane proteins implicated in cell-cell adhesion, and cell-cell interaction such as the coagulation factors V (FaV) and VIII (FaVIII), milk fat globule, neuropilin and neurexins, and implicated in phospholipid binding (Adams et al., 2004). In different species they have 47 highly conserved residues (~30% sequence identity) including 2 cysteines forming a disulphide bridge and several tryptophans maintaining a stable hydrophobic core of a discoidin domain β-sandwich (Kiedzierska et al., 2007). The structure of mature retinoschisin is shown in Figure 4A. In this structure the Rs1 domain is located in residues 25-64 and the folded, β-sandwich discoidin domain occupies residues 65-224 (Wu and Molday, 2003). The structure of the mature retinoschisin was equilibrated using molecular dynamics (MD). Usually MD is used to obtain an ‘idealized’ or equilibrated protein structure with stereochemistry close to perfect. The main justification of the MD method is that statistical ensemble averages are equal to time averages of the system, so-called the ergodic hypothesis. MD by computer simulation of physical movements by atoms and molecules take into account physical interactions between atoms of the protein structure and molecules of surrounding waters. Whereas it is possible to take "still snapshots" of crystal structures using X-Ray crystallography and NMR, no current experimental technique allows access to all the time scales of motion at atomic resolution.
MD can be used to refine the location of cysteine residues in the RS1 structure which is necessary to understand retinoschisin function (Wang et al., 2002; Wang et al., 2006; Wu and Molday, 2003; Wu et al., 2005). The structure of mature RS1 contains 10 cysteines in positions 38, 40, 42, 59, 63, 83, 110, 142, 219, and 223, which are shared in an even proportion between Rs1 and discoidin domains (Figure 4A). The structure of the discoidin domain and the positions of 5 cysteines were reported in several cases (Fraternali et al., 2003; Wang et al., 2006; Wu and Molday, 2003). The other 5 cysteines are located in the Rs1 domain of mature retinoschisin. In a protein structure, 7 cysteines are predicted to be concentrated in an area close to amino- and carboxy-termini, and 3 cysteines are located close to a spike area shown at the top (terminal extensions) and the bottom (discoidin domain) parts of Figure 4A, respectively. Nearly 70% of the RS1 surface is formed by non-polar groups (see Figures 4B and C, rotation of ~180o) suggesting that non-polar groups (green) could be important for the RS1 normal function.
Although numerous experiments have demonstrated that retinoschisin is secreted as a monomer, dimer, or octamer in model cell systems (Wang et al., 2006; Wu et al., 2005) this protein is secreted
3.4. Computational assessment of severity of protein damage caused by the pathogenic missense change
in the last decade, new approaches using semi-empirical force-fields were developed in molecular modeling. Stability of native protein structure and protein complexes is defined by non-polar interactions, disulphide covalent bonds, salt bridges, and hydrogen bonds. Semi-empirical methods imply simplified forces and potentials are taken into account for these interactions within protein structure. Semi-empirical methods generate values of protein stabilization energies that correlate well with an experimental data for a majority of tested proteins. This approach might help to grade severities of molecular perturbations due to missense changes and explain how a genetic mutation could impact protein structure in quantitative manner.
The free energy of protein stabilization (in kcal mol-1) is a function of several parameters. The first term includes a solvation energy contribution from hydrophobic and polar groups of protein (Eisenberg and McLachlan, 1986) based on a concept of accessible surface area (Lee and Richards, 1971) which could be evaluated using a fast geometry algorithm (Tsodikov et al., 2002). The second term is the energy of hydrogen bonds with regard to simple geometric considerations (Petukhov et al., 1999). Additional terms include Van der Waals energy which takes into account the experimental transfer energies from water to vapor; the electrostatic coulomb free energy term, which describes electrostatic interactions in protein; the crude entropy of protein chain to obtain a measure of free energy and free energy of the steric overlaps between atoms in protein structure (Abagyan and Totrov, 1994; Munoz and Serrano, 1994).
Currently, all these terms are combined in the FoldX force field (Schymkowitz et al., 2005; Schymkowitz et al., 2002). The rapid evaluation of protein stabilization free energy change due to missense mutations could be accessed with this field based on an empirical energy function derived from experimental work on proteins. The predictive power of this approach has been successfully tested with experimental thermodynamics data of 1000 mutant proteins (Guerois et al., 2002) from the ProTherm database (Gromiha et al., 2000; Gromiha et al., 1999; Gromiha et al., 2002). Recently a combination of molecular dynamics folding simulations and FoldX potential has been suggested as an approach of
In addition, we used information about differences between amino acids replacing each other in order to achieve a correct estimate for mutations involving cysteine residues. The amino acids differ by combination of several chemical factors like composition, polarity and molecular volume which correlate best with protein residue substitution frequencies as observed (Grantham, 1974). The Grantham difference for each amino acid pair could be used to characterize a physiochemical differences between wild-type and missense variant amino acid residues (Karchin et al., 2007; Tavtigian et al., 2008).
Therefore, we suggested a molecular grading scale, so-called a computed impact score,
3.5. Molecular grading scale: protein structure-based classification of the pathogenic missense mutations impact in human XLRS
We applied molecular modeling of the RS protein and consider perturbations caused by mutations found in human XLRS affected subjects. Computational analysis has predicted that pathogenic missense mutations could be separated in two groups by using the molecular grading scale.
The ERG b/a-wave ratios were sorted by patient age and by the impact of mutation on protein atomic structure. This analysis was applied independently to ERG b/a-wave ratios as a function of average patient age versus severity of XLRS phenotype (Figure 5A). To minimize the effect of statistical errors and to improve a signal-to-noise ratio in experimental ERG data, we used the method described previously (Sergeev et al., 2010) to a new XLRS cohort of patients. All b/a-wave ratios were divided into groups by age (younger, older) and by mutation severity (low-impact, high-impact) (Figure 5A). The b/a-wave ratio in older (37-43 y/o) XLRS patients was greater than younger (13-17 y/o) for severe missense variants with the b-wave unchanged as the a-wave declined similarly to that observed in the
Differences in phenotypes were verified for low- and high-impact missense changes and severe non-missense mutations (Figure 5B). Low-impact missense changes show average b/a-wave ratio close to 1.0 for a maximum with 15 XLRS patients of different ages (labeled by star in bottom panel). This ratio corresponds to a maximum of 10 patients with an average value of approximately 0.6, which is consistent with more severe XLRS phenotype as expected for high-impact missense variants. In addition, the severity of high impact missense changes is confirmed by identical correspondence to 15 patients with the severe non-missense variants (top panel). Thus, low- and high-impact missense changes in protein structure might be associated with mild and severe changes in the XLRS phenotype.
Earlier we demonstrated that maximum structure perturbations associated with severe XLRS phenotype arise from either the removal or insertion of cysteine residues or from changes in the hydrophobic core (Sergeev et al., 2010). Less severe or mild changes were associated with missense mutations targeting a significant part of protein surface. Our present study, which was uses different datasets, has demonstrated that mild phenotypic changes could be associated with mutations in the water inaccessible residues when polar-to-polar or the hydrophobic-to-hydrophobic changes has appeared.
Our computational analysis has predicted that pathogenic missense mutations could be separated in mild and severe groups of XLRS phenotype by using the molecular grading scale. The results for the expected phenotypes of 109 missense changes which are currently available for XLRS are shown in the Table 3. In order to understand how the phenotype of XLRS disease described by b/a-wave ratios is related to missense changes, we analyzed mutant residue accessibilities (ASA) in the RS1 structure similar to that of shown in the Table 2. All mutations were ranked according to their computed impact. These mutations were separated into two groups, low-impact and high-impact, by their computed severity threshold (~0.4) as previously demonstrated (Sergeev et al., 2010). From this analysis, half of the mild mutations could be considered as mutations which were partially or fully exposed at the protein surface; another half are the mutations predicted to be buried in the hydrophobic core. Mutations from this latter group were variants related to insertion, exchange, or removal of charged residues or replacement with a homologous hydrophobic or polar residue. In only two cases, C223S and C83S, cysteines were affected by a mild missense change. In contrast, the majority of missense variants were predicted to cause severe changes, related to insertion or removal of cysteine residues (23 mutations) and/or buried in the hydrophobic core.
|XLRS phenotype||Missense variants*)|
|Mild||R209H, D158N, R213Q, S73P, E146D, G70S, T185K, H207Q, I133F, A101P, R141H, E72K, G178D, E72D, Q154R, D168H, L113F, H207D, E72Q, E146K, G70A, A211T, G109A, R200H, R209L, E215Q, G109E, G140E, K167N, H222N, N104K, R141G, R102Q, R213W, R102W, P203L, T138A, P192A, A98E, G109R, C223S, W96R, P192T, I199T, P193L, R191P, C83S, A211F, P192S, R102E, E215K, D145H, P192R|
|Severe||C59S, R197P, G70R, I195V, P192L, N179D, L69P, R197H, V76G, G74V, G135V, G140R, L137P, R209P, L216P, R209G, E72A, P193S, R156G, E215V, C219S, I136T, D143V, P68R, C142R, G109W, C110R, L103R, Y89C, F176S, C223R, C110Y, C219G, R209C, N99C, R200C, L127P, R197C, G70D, C219W, R182C, C219R, R141C, W163C, W122C, Y65C, W112C, Y155C, I81N, I81D, W92C, W96C, F108C, C142W, C223Y|
Finally, this study confirms results of our previous work that severe phenotype maximum structure perturbations are related to dramatic changes in a protein hydrophobic core or to the deletion or insertion of cysteine residues affecting in general the stability of protein fold. The kinetics of the disease progression might depend on the degree of the mutation impact on protein stability caused by the pathogenic missense change. Thus, we have to incorporate age, genotype, and molecular modeling into ERG analysis to understand a functional role of pathogenic mutations.
4. Incorporating age and genotype into ERG analysis
In our study, the effect of high-impact missense mutations could be recognized as severe at an earlier age when compared to low-impact missense mutations. In order to recognize the predicted affect of mutations we look to kinetics of the disease progression and incorporate patient age into the analysis of the relationship between genotype and phenotype. An exact understanding of how XLRS disease progresses with age is not clear. A 15-year, longitudinal XLRS study suggests that XLRS is a more stationary disease up to age 35 years, as seen in baseline and follow up ERG testing of 10 subjects. The first examination was performed on subjects between 6 to 15 years of age. A follow up examination was preformed 8 to 14 years later. In both exams, all 10 subjects showed reduced b-wave amplitudes, and no difference was found between the b-wave amplitude of the first or second visit (Kjellstrom et al., 2010).
With these studies in mind, we raised the question of what role XLRS mutation severity plays in the dark-adapted ERG response. We preformed a retrospective cross-sectional study on 73 XLRS subjects who were seen at the National Eye Institute between 2004 and 2010. All subjects underwent ocular examination including a full-field ERG with ISCEV standard, dark-adapted response. All subjects were genotyped for
Our cohort also displayed an earlier aging effect due to mutation than what would be expected from aging alone. The rise and fall of the b-wave amplitude with age in normal populations has been thoroughly studied (Birch and Anderson, 1992). At birth, the maximum response is minimal, but increases approximately 5 fold by age 4 months. By 5 years of age, children have normal adult b-wave amplitudes. The b-wave was not statistically affected from 5 years to 55 years of age, showing less than 25 µV decline per decade. Between ages 45-55 and 55-65, a significant decline of approximately 100 µV b-wave amplitude was seen. Birch’s study showed that subjects ages 68 years have b-wave amplitudes approximately half those of subjects between the ages of 15 to 24 years.
For some time, investigators have questioned the possibility of photoreceptor disease in XLRS. One study found decreased a-wave amplitudes in 1/3 of their XLRS cohort, suggesting that rod photoreceptor dysfunction or loss could be a more primary manifestation of XLRS that first thought (Bradshaw et al., 1999). However, the majority of studies suggest complete or relative preservation of the a-wave in XLRS (Eksandh et al., 2000b; Kellner et al., 1990; Peachey et al., 1987). To look further into the question, Khan et al. modeled the leading edge of the a-wave to assess for pathology in the photo-transduction pathway (Khan et al., 2001). By modeling the first 15 ms of the waveform Khan found no difference between XLRS and normal subjects in sensitivity (S) and photoreceptor maximum amplitude (Kmax) in the majority of XLRS waveforms measured.
4.1. The b/a-wave ratio can be a useful tool to quantify rod bipolar cell dysfunction relationship to rod function
To account for any compromise of the photoreceptor layer in XLRS, a ratio of photoreceptor-to-bipolar cell synaptic transmission compared to photo-transduction quality can be used (b/a-wave ratio). One can use the b/a-wave ratio to adjust for any photoreceptor dysfunction affect on the b-wave amplitude and can amplify the quantity of dysfunction in bipolar cells.
A reduced b/a-ratio is a characteristic of XLRS. The clinical lower limit of normal for the b/a ratio is 1.2 (Sergeev et al., 2010). Given normal b-wave amplitudes previous discussed, one would expect preservation of b/a-ratios in normal subjects until approximately 55 years of age (Birch and Anderson, 1992). In our study, no change in b/a ratio was due to aging alone (n=143 eyes, difference = 0.0039/year, p=0.17; excluding those >55 y/o); however, subjects with low-impact mutations maintained an overall higher b/a-ratio up to 55 years of age than did subjects with high-impact mutations (difference 0.3, p = 0.005, less severe = 54 eyes, more severe = 79 eyes; excluding those >55 y/o; Figure 6).
At first glance, it might appear that our study contradicts the XLRS longitudinal study by Kjellstrom (2010). However, had the longitudinal study identified the genetic mutations of all subjects (only three subjects were genotyped), a decrease in b/a-ratio could have been expected in subjects with severe mutations. Of the three subjects who underwent genotyping, two had severe mutations and had decreased b/a-wave ratio at the follow up visit (ratios decreased from 1.2 to 0.8, and 1.1 to 0.77). According to our data, subjects with low-impact mutations would have only a modest change within the 10 to 15 year followup period.
Our computational analysis has suggested that the disease risk assessment at the atomic level could be implied in clinical studies. The combination of clinical electrophysiology, biology of retinal disease and molecular modeling provide a useful approach for evaluation of possible associations between the predicted structural alteration and/or damage to retinoschisin and the severity of XLRS disease.
We thank Sean Finnegan from the National Eye Institute for his valuable editorial help.